Multidimensional Langevin modeling of biomolecular dynamics.

A systematic computational approach to describe the conformational dynamics of biomolecules in reduced dimensionality is presented. The method is based on (i) the decomposition of a high-dimensional molecular dynamics trajectory into a few "system" and (many) "bath" degrees of freedom and (ii) a Langevin simulation of the resulting model. Employing principal component analysis, the dimension of the system is chosen such that it contains all slow large-amplitude motions of the molecule, while the bath coordinates only account for its high-frequency fluctuations. It is shown that a sufficiently large dimension of the model is essential to ensure a clear time scale separation of system and bath variables, which warrants the validity of the memory-free Langevin equation. Applying methods from nonlinear time series analysis, a practical Langevin algorithm is presented which performs a local estimation of the multidimensional Langevin vector fields describing deterministic drift and stochastic driving. Adopting a 800 ns molecular dynamics simulation of the folding of heptaalanine in explicit water, it is shown that a five-dimensional Langevin model correctly reproduces the structure and conformational dynamics of the system. The virtues and limits of the approach are discussed in some detail.

[1]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[2]  Herbert Levine,et al.  Effective stochastic dynamics on a protein folding energy landscape. , 2006, The Journal of chemical physics.

[3]  William Swope,et al.  Describing Protein Folding Kinetics by Molecular Dynamics Simulations. 1. Theory , 2004 .

[4]  Jeremy C. Smith,et al.  Transition Networks for the Comprehensive Characterization of Complex Conformational Change in Proteins. , 2006, Journal of chemical theory and computation.

[5]  B. Jones,et al.  Lecture notes in mathematics: rudiments of Riemann surfaces , 1971 .

[6]  D. Sherrington.  Stochastic Processes in Physics and Chemistry , 1983 .

[7]  Carsten Hartmann,et al.  Data-based parameter estimation of generalized multidimensional Langevin processes. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[8]  García,et al.  Large-amplitude nonlinear motions in proteins. , 1992, Physical review letters.

[9]  I. Jolliffe.  Principal Component Analysis , 2002 .

[10]  Gerhard Stock,et al.  Conformational dynamics of trialanine in water. 2. Comparison of AMBER, CHARMM, GROMOS, and OPLS force fields to NMR and infrared experiments , 2003 .

[11]  D Thirumalai,et al.  Monomer adds to preformed structured oligomers of Aβ-peptides by a two-stage dock–lock mechanism , 2007, Proceedings of the National Academy of Sciences.

[12]  Illia Horenko,et al.  Multiscale Modelling in Molecular Dynamics: Biomolecular Conformations as Metastable States , 2006 .

[13]  H. Berendsen,et al.  Essential dynamics of proteins , 1993, Proteins.

[14]  J. Onuchic,et al.  Theory of protein folding: the energy landscape perspective. , 1997, Annual review of physical chemistry.

[15]  S. Gnanakaran,et al.  Validation of an all-atom protein force field: From dipeptides to larger peptides , 2003 .

[16]  G. Stock,et al.  Conformational Dynamics of Trialanine in Water: A Molecular Dynamics Study , 2002 .

[17]  H Kantz,et al.  Indispensable finite time corrections for Fokker-Planck equations from time series data. , 2001, Physical review letters.

[18]  R. Hegger,et al.  Dihedral angle principal component analysis of molecular dynamics simulations. , 2007, The Journal of chemical physics.

[19]  M. Gruebele.  Protein folding: the free energy surface. , 2002, Current opinion in structural biology.

[20]  H. Schwalbe,et al.  Structure and dynamics of the homologous series of alanine peptides: a joint molecular dynamics/NMR study. , 2007, Journal of the American Chemical Society.

[21]  Graham Richards,et al.  Intermolecular forces , 1978, Nature.

[22]  Gerhard Stock,et al.  How complex is the dynamics of Peptide folding? , 2007, Physical review letters.

[23]  Jeremy C. Smith,et al.  Hierarchical analysis of conformational dynamics in biomolecules: transition networks of metastable states. , 2007, The Journal of chemical physics.

[24]  P. Nguyen,et al.  Energy landscape of a small peptide revealed by dihedral angle principal component analysis , 2004, Proteins.

[25]  Gerhard Stock,et al.  Construction of the free energy landscape of biomolecules via dihedral angle principal component analysis. , 2008, The Journal of chemical physics.

[26]  B. L. de Groot,et al.  Essential dynamics of reversible peptide folding: memory-free conformational dynamics governed by internal hydrogen bonds. , 2001, Journal of molecular biology.

[27]  Oliver F. Lange,et al.  Collective Langevin dynamics of conformational motions in proteins. , 2006, The Journal of chemical physics.

[28]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[29]  Gerhard Stock,et al.  Peptide conformational heterogeneity revealed from nonlinear vibrational spectroscopy and molecular dynamics simulations , 2002 .

[30]  K. Dill,et al.  From Levinthal to pathways to funnels , 1997, Nature Structural Biology.

[31]  N. Go,et al.  Investigating protein dynamics in collective coordinate space. , 1999, Current opinion in structural biology.

[32]  K. Dill,et al.  Automatic discovery of metastable states for the construction of Markov models of macromolecular conformational dynamics. , 2007, The Journal of chemical physics.

[33]  H. Kantz,et al.  Nonlinear time series analysis , 1997 .

[34]  J. A. Hartigan,et al.  A k-means clustering algorithm , 1979 .

[35]  B. Romanowicz,et al.  Anisotropy in the Inner Core: Could It Be Due To Low-Order Convection? , 1996, Science.

[36]  P. Grassberger.  An optimized box-assisted algorithm for fractal dimensions , 1990 .

[37]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[38]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[39]  Paul Tavan,et al.  Extracting Markov Models of Peptide Conformational Dynamics from Simulation Data. , 2005, Journal of chemical theory and computation.

[40]  U. Weiss.  Quantum Dissipative Systems , 1993 .

[41]  Gradisek,et al.  Analysis of time series from stochastic processes , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[42]  John D. Chodera,et al.  Long-Time Protein Folding Dynamics from Short-Time Molecular Dynamics Simulations , 2006, Multiscale Model. Simul..

[43]  H. Berendsen,et al.  Interaction Models for Water in Relation to Protein Hydration , 1981 .

[44]  A Mitsutake,et al.  Generalized-ensemble algorithms for molecular simulations of biopolymers. , 2000, Biopolymers.

[45]  N. Mavromatos,et al.  LECT NOTES PHYS , 2002 .

[46]  M. Karplus,et al.  Collective motions in proteins: A covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations , 1991, Proteins.

[47]  Peter Hamm,et al.  Structure Determination of Trialanine in Water Using Polarization Sensitive Two-Dimensional Vibrational Spectroscopy , 2000 .

[48]  D. Leitner,et al.  Free energy landscape of a biomolecule in dihedral principal component space: Sampling convergence and correspondence between structures and minima , 2007, Proteins.

[49]  R Schweitzer-Stenner,et al.  Dihedral angles of trialanine in D2O determined by combining FTIR and polarized visible Raman spectroscopy. , 2001, Journal of the American Chemical Society.

[50]  P. Deuflhard,et al.  A Direct Approach to Conformational Dynamics Based on Hybrid Monte Carlo , 1999 .

[51]  R. Zwanzig.  Nonequilibrium statistical mechanics , 2001, Physics Subject Headings (PhySH).

[52]  F. Takens.  Detecting strange attractors in turbulence , 1981 .

[53]  Chin-Kun Hu,et al.  Free energy landscape and folding mechanism of a β‐hairpin in explicit water: A replica exchange molecular dynamics study , 2005, Proteins.