title: Nonlinear vs. linear biasing in Trp-cage folding simulations.
abstract: Biased simulations have great potential for the study of slow processes, including protein folding. Atomic motions in molecules are nonlinear, which suggests that simulations with enhanced sampling of collective motions traced by nonlinear dimensionality reduction methods may perform better than linear ones. In this study, we compare an unbiased folding simulation of the Trp-cage miniprotein with metadynamics simulations using both linear (principle component analysis) and nonlinear (Isomap) low dimensional embeddings as collective variables. Folding of the mini-protein was successfully simulated in 200 ns simulation with linear biasing and non-linear motion biasing. The folded state was correctly predicted as the free energy minimum in both simulations. We found that the advantage of linear motion biasing is that it can sample a larger conformational space, whereas the advantage of nonlinear motion biasing lies in slightly better resolution of the resulting free energy surface. In terms of sampling efficiency, both methods are comparable.
authors: Department of Biochemistry and Microbiology, University of Chemistry and Technology, Prague, Technická 3, Prague 6 166 28, Czech Republic.
authors: Vojtěch
authors: V
authors: Spiwok
authors: Department of Biochemistry and Microbiology, University of Chemistry and Technology, Prague, Technická 3, Prague 6 166 28, Czech Republic.
authors: Pavel
authors: P
authors: Oborský
authors: Institute of Computer Science, Masaryk University, Botanická 554/68a, 602 00 Brno, Czech Republic.
authors: Jana
authors: J
authors: Pazúriková
authors: Institute of Computer Science, Masaryk University, Botanická 554/68a, 602 00 Brno, Czech Republic.
authors: Aleš
authors: A
authors: Křenek
authors: Department of Biochemistry and Microbiology, University of Chemistry and Technology, Prague, Technická 3, Prague 6 166 28, Czech Republic.
authors: Blanka
authors: B
authors: Králová
mesh_terms: D003198:Computer Simulation; D016014:Linear Models; D008956:Models, Chemical; D009038:Motion; D017711:Nonlinear Dynamics; D017510:Protein Folding; D011506:Proteins / Q000737:chemistry; D012997:Solvents / Q000737:chemistry
chemical_list: D011506:Proteins; D012997:Solvents