Error and efficiency of replica exchange molecular dynamics simulations.

We derive simple analytical expressions for the error and computational efficiency of replica exchange molecular dynamics (REMD) simulations (and by analogy replica exchange Monte Carlo simulations). The theory applies to the important case of systems whose dynamics at long times is dominated by the slow interconversion between two metastable states. As a specific example, we consider the folding and unfolding of a protein. The efficiency is defined as the rate with which the error in an estimated equilibrium property, as measured by the variance of the estimator over repeated simulations, decreases with simulation time. For two-state systems, this rate is in general independent of the particular property. Our main result is that, with comparable computational resources used, the relative efficiency of REMD and molecular dynamics (MD) simulations is given by the ratio of the number of transitions between the two states averaged over all replicas at the different temperatures, and the number of transitions at the single temperature of the MD run. This formula applies if replica exchange is frequent, as compared to the transition times. High efficiency of REMD is thus achieved by including replica temperatures in which the frequency of transitions is higher than that at the temperature of interest. In tests of the expressions for the error in the estimator, computational efficiency, and the rate of equilibration we find quantitative agreement with the results both from kinetic models of REMD and from actual all-atom simulations of the folding of a peptide in water.

[1]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[2]  R. Levy,et al.  Simple continuous and discrete models for simulating replica exchange simulations of protein folding. , 2008, The journal of physical chemistry. B.

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

[4]  C. Zhang,et al.  Comparison of sampling efficiency between simulated tempering and replica exchange. , 2008, The Journal of chemical physics.

[5]  Gerhard Hummer,et al.  Peptide folding kinetics from replica exchange molecular dynamics. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[6]  G. Hummer,et al.  Thermostat artifacts in replica exchange molecular dynamics simulations. , 2009, Journal of chemical theory and computation.

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

[8]  G. Hummer,et al.  Coarse master equations for peptide folding dynamics. , 2008, The journal of physical chemistry. B.

[9]  Wang,et al.  Replica Monte Carlo simulation of spin glasses. , 1986, Physical review letters.

[10]  A. Szabó,et al.  Model-free approach to the interpretation of nuclear magnetic resonance relaxation in macromolecules. 1. Theory and range of validity , 1982 .

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

[12]  M. Troyer,et al.  Optimized parallel tempering simulations of proteins. , 2006, The Journal of chemical physics.

[13]  Gerhard Hummer,et al.  Diffusion and electrophoretic mobility of single-stranded RNA from molecular dynamics simulations. , 2004, Biophysical journal.

[14]  George H. Weiss,et al.  Theory of single-molecule fluorescence spectroscopy of two-state systems , 1999 .

[15]  Sanford Weisberg,et al.  Computing science and statistics : proceedings of the 30th Symposium on the Interface, Minneapolis, Minnesota, May 13-16, 1998 : dimension reduction, computational complexity and information , 1998 .

[16]  B. Berne,et al.  Replica exchange with solute tempering: a method for sampling biological systems in explicit water. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[17]  K. Sanbonmatsu,et al.  α-Helical stabilization by side chain shielding of backbone hydrogen bonds , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[18]  Daniel J. Sindhikara,et al.  Exchange frequency in replica exchange molecular dynamics. , 2008, The Journal of chemical physics.

[19]  Chris Oostenbrink,et al.  Optimization of replica exchange molecular dynamics by fast mimicking. , 2007, The Journal of chemical physics.

[20]  Paul Tavan,et al.  Optimal temperature ladders in replica exchange simulations , 2009 .

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

[22]  W. Nadler,et al.  Optimized explicit-solvent replica exchange molecular dynamics from scratch. , 2008, The journal of physical chemistry. B.

[23]  Ulrich H E Hansmann,et al.  Dynamics and optimal number of replicas in parallel tempering simulations. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[24]  Paul Tavan,et al.  Efficiency reduction and pseudo-convergence in replica exchange sampling of peptide folding–unfolding equilibria , 2008 .

[25]  A. Berezhkovskii,et al.  Rate Constants for Slow Conformational Transitions and Their Sampling Errors Using Single-Molecule Fluorescence Spectroscopy , 2001 .

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

[27]  A. Szabó,et al.  Electron transfer reaction dynamics in non-Debye solvents , 1998 .

[28]  Martin Gruebele,et al.  Rate−Temperature Relationships in λ-Repressor Fragment λ6-85 Folding† , 2004 .

[29]  Ulrich H E Hansmann,et al.  Generalized ensemble and tempering simulations: a unified view. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[30]  Mark J Abraham,et al.  Ensuring Mixing Efficiency of Replica-Exchange Molecular Dynamics Simulations. , 2008, Journal of chemical theory and computation.