Optimization of Umbrella Sampling Replica Exchange Molecular Dynamics by Replica Positioning.

The positioning of sampling windows in an umbrella sampling simulation has an effect on the rate of convergence and computational efficiency. When such simulation is coupled with a Hamiltonian replica exchange setup, we show that such positioning can be optimized for maximal convergence of the results. We present a method for estimating the exchange acceptance ratio (EAR) between two arbitrary positions on a reaction coordinate in umbrella sampling replica exchange (USRE) molecular dynamics (MD). We designed a scoring function to optimize the position of the set of replicas (windows). By maximizing the scoring function, we make EAR the same for all neighbor replica pairs, increasing the efficiency of the method. We tested our algorithm by sampling a torsion for butane in implicit solvent and by studying a salt bridge in explicit solvent. We found that the optimized set of replicas recovers the correct free energy profile much faster than for equally spaced umbrellas.

[1]  R. A. Leibler,et al.  On Information and Sufficiency , 1951 .

[2]  David A. Kofke,et al.  ARTICLES On the acceptance probability of replica-exchange Monte Carlo trials , 2002 .

[3]  B. Berne,et al.  Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[4]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[5]  Danial S. Dashti,et al.  Computing Alchemical Free Energy Differences with Hamiltonian Replica Exchange Molecular Dynamics (H-REMD) Simulations. , 2011, Journal of chemical theory and computation.

[6]  Wei Yang,et al.  Random walk in orthogonal space to achieve efficient free-energy simulation of complex systems , 2008, Proceedings of the National Academy of Sciences.

[7]  A. Iusem,et al.  Full convergence of the steepest descent method with inexact line searches , 1995 .

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

[9]  J. Pablo,et al.  Multicanonical parallel tempering , 2002, cond-mat/0201179.

[10]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[11]  C. Bartels.  Analyzing biased Monte Carlo and molecular dynamics simulations , 2000 .

[12]  C. Predescu,et al.  The incomplete beta function law for parallel tempering sampling of classical canonical systems. , 2003, The Journal of chemical physics.

[13]  R. Zhou.  Free energy landscape of protein folding in water: Explicit vs. implicit solvent , 2003, Proteins.

[14]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

[15]  Daniel J. Sindhikara,et al.  Bad Seeds Sprout Perilous Dynamics: Stochastic Thermostat Induced Trajectory Synchronization in Biomolecules. , 2009, Journal of chemical theory and computation.

[16]  Andrew W. Moore,et al.  Locally Weighted Learning , 1997, Artificial Intelligence Review.

[17]  Edward Lyman,et al.  Resolution exchange simulation. , 2005, Physical review letters.

[18]  A. Roitberg,et al.  pH-replica exchange molecular dynamics in proteins using a discrete protonation method. , 2012, The journal of physical chemistry. B.

[19]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[20]  Hiqmet Kamberaj,et al.  An optimized replica exchange molecular dynamics method. , 2009, The Journal of chemical physics.

[21]  Nadeem A. Vellore,et al.  An improved replica-exchange sampling method: temperature intervals with global energy reassignment. , 2007, The Journal of chemical physics.

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

[23]  Michael P Eastwood,et al.  Minimizing thermodynamic length to select intermediate states for free-energy calculations and replica-exchange simulations. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[25]  Asim Okur,et al.  Improved Efficiency of Replica Exchange Simulations through Use of a Hybrid Explicit/Implicit Solvation Model. , 2006, Journal of chemical theory and computation.

[26]  Walter Thiel,et al.  Bridging the gap between thermodynamic integration and umbrella sampling provides a novel analysis method: "Umbrella integration". , 2005, The Journal of chemical physics.

[27]  J. D. de Pablo,et al.  Optimal allocation of replicas in parallel tempering simulations. , 2005, The Journal of chemical physics.

[28]  S. Takada,et al.  On the Hamiltonian replica exchange method for efficient sampling of biomolecular systems: Application to protein structure prediction , 2002 .

[29]  Gavin E Crooks,et al.  Measuring thermodynamic length. , 2007, Physical review letters.

[30]  Johannes Kästner.  Umbrella integration with higher-order correction terms. , 2012, The Journal of chemical physics.

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

[32]  Michael W Deem,et al.  Parallel tempering: theory, applications, and new perspectives. , 2005, Physical chemistry chemical physics : PCCP.

[33]  J. Neumann.  Proof of the Quasi-Ergodic Hypothesis. , 1932, Proceedings of the National Academy of Sciences of the United States of America.

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

[35]  R. Swendsen,et al.  THE weighted histogram analysis method for free‐energy calculations on biomolecules. I. The method , 1992 .

[36]  Bernard R Brooks,et al.  pH replica‐exchange method based on discrete protonation states , 2011, Proteins.

[37]  Viktor Hornak,et al.  Improving Convergence of Replica-Exchange Simulations through Coupling to a High-Temperature Structure Reservoir. , 2007, Journal of chemical theory and computation.

[38]  D. Kofke,et al.  Selection of temperature intervals for parallel-tempering simulations. , 2005, The Journal of chemical physics.

[39]  V. Hornak,et al.  Modified replica exchange simulation methods for local structure refinement. , 2005, The journal of physical chemistry. B.

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

[41]  Matthias Troyer,et al.  Feedback-optimized parallel tempering Monte Carlo , 2006, cond-mat/0602085.

[42]  Y. Sugita,et al.  Multidimensional replica-exchange method for free-energy calculations , 2000, cond-mat/0009120.

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

[44]  Matthew P Jacobson,et al.  Comparing Conformational Ensembles Using the Kullback-Leibler Divergence Expansion. , 2012, Journal of chemical theory and computation.

[45]  S. Rick,et al.  Replica exchange with dynamical scaling. , 2007, The Journal of chemical physics.

[46]  D. Kofke.  Comment on "The incomplete beta function law for parallel tempering sampling of classical canonical systems" [J. Chem. Phys. 120, 4119 (2004)]. , 2004, The Journal of chemical physics.

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

[48]  F. Escobedo,et al.  Expanded ensemble and replica exchange methods for simulation of protein-like systems , 2003 .

[49]  Daniel J. Sindhikara,et al.  Exchange Often and Properly in Replica Exchange Molecular Dynamics. , 2010, Journal of chemical theory and computation.

[50]  Michael R. Shirts,et al.  Replica exchange and expanded ensemble simulations as Gibbs sampling: simple improvements for enhanced mixing. , 2011, The Journal of chemical physics.

[51]  G. Torrie,et al.  Nonphysical sampling distributions in Monte Carlo free-energy estimation: Umbrella sampling , 1977 .

[52]  D. Case,et al.  Modification of the Generalized Born Model Suitable for Macromolecules , 2000 .

[53]  Andrew J Ballard,et al.  Replica exchange with nonequilibrium switches , 2009, Proceedings of the National Academy of Sciences.

[54]  Y. Sugita,et al.  Free energy calculations for DNA base stacking by replica-exchange umbrella sampling , 2004 .

[55]  J. Kästner.  Umbrella sampling , 2011 .

[56]  A. Roitberg,et al.  Coupling of replica exchange simulations to a non-Boltzmann structure reservoir. , 2007, The journal of physical chemistry. B.