Multiscale Modelling in Molecular Dynamics: Biomolecular Conformations as Metastable States

We report on a novel approach to the automatic identification of metastable states from long term simulation of complex molecular systems. The new approach is based on a hierarchical concept of metastability: metastable states are understood as subsets of state or configuration space from which the dynamics exits only very rarely; subsets with the smallest exit probabilities are of most interest, their further decomposition then may reveal subsets from which exiting is less but comparably difficult for the system under investigation. The article gives a survey of the theoretical foundation of the approach and its algorithmic realization that generalizes the well-known concept of Hidden Markov Models. The performance of the resulting algorithm are illustrated by application to a 100 ns simulation of penta-alanine with explicit water. We demonstrate the resulting metastable states allow to reveal the conformation dynamics of the moelcule.

[1]  A. Bovier,et al.  Metastability in stochastic dynamics of disordered mean-field models , 1998, cond-mat/9811331.

[2]  Marcus Weber Improved Perron Cluster Analysis , 2003 .

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

[4]  Peter Deuflhard,et al.  Transfer Operator Approach to Conformational Dynamics in Biomolecular Systems , 2001 .

[5]  W. E,et al.  Finite temperature string method for the study of rare events. , 2002, Journal of Physical Chemistry B.

[6]  P. Deuflhard,et al.  Identification of almost invariant aggregates in reversible nearly uncoupled Markov chains , 2000 .

[7]  David Chandler,et al.  Finding transition pathways: throwing ropes over rough mountain passes, in the dark , 1998 .

[8]  E. Vanden-Eijnden,et al.  String method for the study of rare events , 2002, cond-mat/0205527.

[9]  B. Berne,et al.  Novel methods of sampling phase space in the simulation of biological systems. , 1997, Current opinion in structural biology.

[10]  P. Deuflhard,et al.  Hierarchical Uncoupling-Coupling of Metastable Conformations , 2002 .

[11]  Illia Horenko,et al.  Macroscopic Dynamics of Complex Metastable Systems: Theory, Algorithms, and Application to B-DNA , 2008, SIAM J. Appl. Dyn. Syst..

[12]  M. Karplus,et al.  Multiple conformational states of proteins: a molecular dynamics analysis of myoglobin. , 1987, Science.

[13]  Van Nostrand,et al.  Error Bounds for Convolutional Codes and an Asymptotically Optimum Decoding Algorithm , 1967 .

[14]  P. Deuflhard,et al.  Robust Perron cluster analysis in conformation dynamics , 2005 .

[15]  G. N. Ramachandran,et al.  Conformation of polypeptides and proteins. , 1968, Advances in protein chemistry.

[16]  O. Junge,et al.  On the Approximation of Complicated Dynamical Behavior , 1999 .

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

[18]  Stephen D. Bond,et al.  The Nosé-Poincaré Method for Constant Temperature Molecular Dynamics , 1999 .

[19]  Eric Vanden Eijnden,et al.  Metastability, conformation dynamics, and transition pathways in complex systems , 2004 .

[20]  Jeff A. Bilmes,et al.  A gentle tutorial of the em algorithm and its application to parameter estimation for Gaussian mixture and hidden Markov models , 1998 .

[21]  Alexander Fischer,et al.  Identification of biomolecular conformations from incomplete torsion angle observations by hidden markov models , 2007, J. Comput. Chem..

[22]  G. Singleton,et al.  ASYMPTOTICALLY EXACT ESTIMATES FOR METASTABLE MARKOV SEMIGROUPS , 1984 .

[23]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[24]  Illia Horenko,et al.  Reduced Stochastic Models for Complex Molecular Systems , 2006 .

[25]  Wilhelm Huisinga,et al.  On Conformational Dynamics induced by Langevin Processes , 2000 .

[26]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Alexander Fischer,et al.  Automated Model Reduction for Complex Systems Exhibiting Metastability , 2006, Multiscale Model. Simul..

[28]  E. Davies,et al.  Metastable States of Symmetric Markov Semigroups I , 1982 .

[29]  S. Meyn,et al.  Phase transitions and metastability in Markovian and molecular systems , 2004 .

[30]  P. Lezaud Chernoff and Berry-Esséen inequalities for Markov processes , 2001 .

[31]  W. Huisinga,et al.  Metastability and Dominant Eigenvalues of Transfer Operators , 2006 .