Replica-Exchange Accelerated Molecular Dynamics (REXAMD) Applied to Thermodynamic Integration

Accelerated molecular dynamics (AMD) is an efficient strategy for accelerating the sampling of molecular dynamics simulations, and observable quantities such as free energies derived on the biased AMD potential can be reweighted to yield results consistent with the original, unmodified potential. In conventional AMD the reweighting procedure has an inherent statistical problem in systems with large acceleration, where the points with the largest biases will dominate the reweighted result and reduce the effective number of data points. We propose a replica exchange of various degrees of acceleration (REXAMD) to retain good statistics while achieving enhanced sampling. The REXAMD method is validated and benchmarked on two simple gas-phase model systems, and two different strategies for computing reweighted averages over a simulation are compared.

[1]  Michael R. Shirts,et al.  Direct calculation of the binding free energies of FKBP ligands. , 2005, The Journal of chemical physics.

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

[3]  J. Mccammon,et al.  Sampling of slow diffusive conformational transitions with accelerated molecular dynamics. , 2007, The Journal of chemical physics.

[4]  G. Bouvignies,et al.  Exploring multiple timescale motions in protein GB3 using accelerated molecular dynamics and NMR spectroscopy. , 2007, Journal of the American Chemical Society.

[5]  Wei Yang,et al.  Simulated scaling method for localized enhanced sampling and simultaneous "alchemical" free energy simulations: a general method for molecular mechanical, quantum mechanical, and quantum mechanical/molecular mechanical simulations. , 2007, The Journal of chemical physics.

[6]  Donald Hamelberg,et al.  Standard free energy of releasing a localized water molecule from the binding pockets of proteins: double-decoupling method. , 2004, Journal of the American Chemical Society.

[7]  J. Mongan,et al.  Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules. , 2004, The Journal of chemical physics.

[8]  Jonathan W. Essex,et al.  The development of replica-exchange-based free-energy methods , 2003 .

[9]  B. Roux,et al.  Absolute binding free energy calculations using molecular dynamics simulations with restraining potentials. , 2006, Biophysical journal.

[10]  Donald Hamelberg,et al.  A statistical analysis of the precision of reweighting-based simulations. , 2008, The Journal of chemical physics.

[11]  Jun Wang,et al.  A Hamiltonian Replica Exchange Approach and Its Application to the Study of Side-Chain Type and Neighbor Effects on Peptide Backbone Conformations. , 2008, Journal of chemical theory and computation.

[12]  Peter A. Kollman,et al.  AMBER, a package of computer programs for applying molecular mechanics, normal mode analysis, molecular dynamics and free energy calculations to simulate the structural and energetic properties of molecules , 1995 .