Simulated tempering based on global balance or detailed balance conditions: Suwa–Todo, heat bath, and Metropolis algorithms

Simulated tempering (ST) is a useful method to enhance sampling of molecular simulations. When ST is used, the Metropolis algorithm, which satisfies the detailed balance condition, is usually applied to calculate the transition probability. Recently, an alternative method that satisfies the global balance condition instead of the detailed balance condition has been proposed by Suwa and Todo. In this study, ST method with the Suwa–Todo algorithm is proposed. Molecular dynamics simulations with ST are performed with three algorithms (the Metropolis, heat bath, and Suwa–Todo algorithms) to calculate the transition probability. Among the three algorithms, the Suwa–Todo algorithm yields the highest acceptance ratio and the shortest autocorrelation time. These suggest that sampling by a ST simulation with the Suwa–Todo algorithm is most efficient. In addition, because the acceptance ratio of the Suwa–Todo algorithm is higher than that of the Metropolis algorithm, the number of temperature states can be reduced by 25% for the Suwa–Todo algorithm when compared with the Metropolis algorithm. © 2015 Wiley Periodicals, Inc.

[1]  Yoshiharu Mori,et al.  Pressure-Induced Helical Structure of a Peptide Studied by Simulated Tempering Molecular Dynamics Simulations. , 2013, The journal of physical chemistry letters.

[2]  Alan M. Ferrenberg,et al.  Optimized Monte Carlo data analysis. , 1989, Physical Review Letters.

[3]  T. Darden,et al.  A smooth particle mesh Ewald method , 1995 .

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

[5]  Jun S. Liu,et al.  The Multiple-Try Method and Local Optimization in Metropolis Sampling , 2000 .

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

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

[8]  Chun-Jung Chen,et al.  Transformation between α‐helix and β‐sheet structures of one and two polyglutamine peptides in explicit water molecules by replica‐exchange molecular dynamics simulations , 2014, J. Comput. Chem..

[9]  P. Kollman,et al.  Settle: An analytical version of the SHAKE and RATTLE algorithm for rigid water models , 1992 .

[10]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[11]  Michael R. Shirts,et al.  Statistically optimal analysis of samples from multiple equilibrium states. , 2008, The Journal of chemical physics.

[12]  Chin-Kun Hu,et al.  Effect of Taiwan mutation (D7H) on structures of amyloid-β peptides: replica exchange molecular dynamics study. , 2014, The journal of physical chemistry. B.

[13]  Hisashi Okumura,et al.  Dimerization process of amyloid-β(29-42) studied by the Hamiltonian replica-permutation molecular dynamics simulations. , 2014, The journal of physical chemistry. B.

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

[15]  Sanghyun Park,et al.  Comparison of the serial and parallel algorithms of generalized ensemble simulations: an analytical approach. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[16]  A. Gorin,et al.  Accelerated simulated tempering , 2004 .

[17]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

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

[19]  A. Lyubartsev,et al.  New approach to Monte Carlo calculation of the free energy: Method of expanded ensembles , 1992 .

[20]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[21]  B. Berg,et al.  Multicanonical algorithms for first order phase transitions , 1991 .

[22]  Hisashi Okumura,et al.  Hamiltonian replica‐permutation method and its applications to an alanine dipeptide and amyloid‐β(29–42) peptides , 2013, J. Comput. Chem..

[23]  Yuko Okamoto,et al.  Communication: Simulated tempering with fast on-the-fly weight determination. , 2013, The Journal of chemical physics.

[24]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[25]  Yuko Okamoto,et al.  Multidimensional generalized-ensemble algorithms for complex systems. , 2009, The Journal of chemical physics.

[26]  W H Wong,et al.  Dynamic weighting in Monte Carlo and optimization. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Chin-Kun Hu,et al.  New force replica exchange method and protein folding pathways probed by force-clamp technique. , 2007, The Journal of chemical physics.

[28]  Y. Okamoto,et al.  Molecular dynamics, Langevin, and hybrid Monte Carlo simulations in multicanonical ensemble , 1996, physics/9710018.

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

[30]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

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

[32]  H. Okumura,et al.  Molecular dynamics of the structural changes of helical peptides induced by pressure , 2014, Proteins.

[33]  A. Kidera,et al.  Multicanonical Ensemble Generated by Molecular Dynamics Simulation for Enhanced Conformational Sampling of Peptides , 1997 .

[34]  Synge Todo,et al.  Markov chain Monte Carlo method without detailed balance. , 2010, Physical review letters.

[35]  G. Parisi,et al.  Simulated tempering: a new Monte Carlo scheme , 1992, hep-lat/9205018.

[36]  Berg,et al.  Multicanonical ensemble: A new approach to simulate first-order phase transitions. , 1992, Physical review letters.

[37]  Hisashi Okumura,et al.  Replica-Permutation Method with the Suwa-Todo Algorithm beyond the Replica-Exchange Method. , 2012, Journal of chemical theory and computation.

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

[39]  Yuko Okamoto,et al.  Generalized-Ensemble Algorithms for the Isobaric–Isothermal Ensemble , 2010, 1004.2076.

[40]  Ken A Dill,et al.  Use of the Weighted Histogram Analysis Method for the Analysis of Simulated and Parallel Tempering Simulations. , 2007, Journal of chemical theory and computation.