Efficiency reduction and pseudo-convergence in replica exchange sampling of peptide folding–unfolding equilibria

Abstract Replica exchange (RE) molecular dynamics (MD) simulations are frequently applied to sample the folding–unfolding equilibria of β-hairpin peptides in solution, because efficiency gains are expected from this technique. Using a three-state Markov model featuring key aspects of β-hairpin folding we show that RE simulations can be less efficient than conventional techniques. Furthermore we demonstrate that one is easily seduced to erroneously assign convergence to the RE sampling, because RE ensembles can rapidly reach long-lived stationary states. We conclude that typical REMD simulations covering a few tens of nanoseconds are by far too short for sufficient sampling of β-hairpin folding–unfolding equilibria.

[1]  G. Stock,et al.  Molecular dynamics simulation of the structure, dynamics, and thermostability of the RNA hairpins uCACGg and cUUCGg. , 2008, The journal of physical chemistry. B.

[2]  K. Kuczera,et al.  Equilibrium structure and folding of a helix-forming peptide: circular dichroism measurements and replica-exchange molecular dynamics simulations. , 2004, Biophysical journal.

[3]  Martin Karplus,et al.  Aspects of Protein Reaction Dynamics: Deviations from Simple Behavior , 2000 .

[4]  Y. Sugita,et al.  Replica-exchange molecular dynamics method for protein folding , 1999 .

[5]  Hugh Nymeyer,et al.  How Efficient Is Replica Exchange Molecular Dynamics? An Analytic Approach. , 2008, Journal of chemical theory and computation.

[6]  William Swope,et al.  Understanding folding and design: Replica-exchange simulations of ``Trp-cage'' miniproteins , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[7]  F. Rao,et al.  Replica exchange molecular dynamics simulations of reversible folding , 2003 .

[8]  G. Stock,et al.  What NMR Relaxation Can Tell Us about the Internal Motion of an RNA Hairpin:  A Molecular Dynamics Simulation Study. , 2006, Journal of chemical theory and computation.

[9]  F. Gai,et al.  Understanding the key factors that control the rate of β-hairpin folding , 2004 .

[10]  Edward Lyman,et al.  A Second Look at Canonical Sampling of Biomolecules using Replica Exchange Simulation. , 2006, Journal of chemical theory and computation.

[11]  K. Hukushima,et al.  Exchange Monte Carlo Method and Application to Spin Glass Simulations , 1995, cond-mat/9512035.

[12]  K. Sanbonmatsu,et al.  Structure of Met‐enkephalin in explicit aqueous solution using replica exchange molecular dynamics , 2002, Proteins.

[13]  B. Berne,et al.  The free energy landscape for β hairpin folding in explicit water , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[14]  R. Levy,et al.  Simulating replica exchange simulations of protein folding with a kinetic network model , 2007, Proceedings of the National Academy of Sciences.

[15]  F. Rao,et al.  Replica exchange molecular dynamics simulations of amyloid peptide aggregation. , 2004, The Journal of chemical physics.

[16]  V. Pande,et al.  Multiplexed-replica exchange molecular dynamics method for protein folding simulation. , 2003, Biophysical journal.

[17]  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.

[18]  P. Tavan,et al.  Light-triggered β-hairpin folding and unfolding , 2007, Proceedings of the National Academy of Sciences.

[19]  M. Gruebele,et al.  Heterogeneous folding of the trpzip hairpin: full atom simulation and experiment. , 2004, Journal of molecular biology.

[20]  U. Hansmann.  Parallel tempering algorithm for conformational studies of biological molecules , 1997, physics/9710041.

[21]  Eric J. Sorin,et al.  How well can simulation predict protein folding kinetics and thermodynamics? , 2005, Annual review of biophysics and biomolecular structure.