TIGER2: an improved algorithm for temperature intervals with global exchange of replicas.

An empirical sampling method for molecular simulation based on "temperature intervals with global exchange of replicas" (TIGER2) has been developed to reduce the high demand for computational resources and the low computational efficiency of the conventional replica-exchange molecular dynamics (REMD) method. This new method overcomes the limitation of its previous version, called TIGER, which requires the assumption of constant heat capacity during quenching of replicas from elevated temperatures to the baseline temperature. The robustness of the TIGER2 method is examined by comparing it against a Metropolis Monte Carlo simulation for sampling the conformational distribution of a single butane molecule in vacuum, a REMD simulation for sampling the behavior of alanine dipeptide in explicit solvent, and REMD simulations for sampling the folding behavior of two peptides, (AAQAA)(3) and chignolin, in implicit solvent. The agreement between the results from these conventional sampling methods and the TIGER2 simulations indicates that the TIGER2 algorithm is able to closely approximate a Boltzmann-weighted ensemble of states for these systems but without the limiting assumptions that were required for the original TIGER algorithm. TIGER2 is an efficient replica-exchange sampling method that enables the number of replicas that are used for a replica-exchange simulation to be substantially reduced compared to the conventional REMD method.

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

[2]  E. Mehler.  Self-Consistent, Free Energy Based Approximation To Calculate pH Dependent Electrostatic Effects in Proteins , 1996 .

[3]  A. Nitzan.  Chemical physics: Ultrafast relaxation in water , 1999, Nature.

[4]  Christoph Dellago,et al.  Sampling ensembles of deterministic transition pathways , 1998 .

[5]  Ruhong Zhou,et al.  Hydrophobic aided replica exchange: an efficient algorithm for protein folding in explicit solvent. , 2006, The journal of physical chemistry. B.

[6]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[7]  J. Skinner,et al.  Vibrational spectroscopy of HOD in liquid D2O. III. Spectral diffusion, and hydrogen-bonding and rotational dynamics , 2003 .

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

[9]  G. N. Ramachandran,et al.  Stereochemistry of polypeptide chain configurations. , 1963, Journal of molecular biology.

[10]  Shinya Honda,et al.  10 residue folded peptide designed by segment statistics. , 2004, Structure.

[11]  D. Beveridge,et al.  Exploratory studies of ab initio protein structure prediction: Multiple copy simulated annealing, AMBER energy functions, and a generalized born/solvent accessibility solvation model , 2002, Proteins.

[12]  J. Apostolakis,et al.  Thermodynamics and Kinetics of Folding of Two Model Peptides Investigated by Molecular Dynamics Simulations , 2000 .

[13]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

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

[15]  H. Berendsen,et al.  ALGORITHMS FOR MACROMOLECULAR DYNAMICS AND CONSTRAINT DYNAMICS , 1977 .

[16]  C. Brooks,et al.  Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field. , 2006, Journal of the American Chemical Society.

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

[18]  F Guarnieri,et al.  A self-consistent, microenvironment modulated screened coulomb potential approximation to calculate pH-dependent electrostatic effects in proteins. , 1999, Biophysical journal.

[19]  G. Thorleifsson,et al.  Improved algorithms for simulating crystalline membranes , 1998 .

[20]  Weber,et al.  Computer simulation of local order in condensed phases of silicon. , 1985, Physical review. B, Condensed matter.

[21]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

[22]  E. Stellwagen,et al.  Distribution of Helicity within the Model Peptide Acetyl(AAQAA)3amide , 1994 .

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

[24]  W A Shirley,et al.  Curious structure in “canonical” alanine‐based peptides , 1997, Proteins.

[25]  F. Stillinger,et al.  Molecular Dynamics Study of Liquid Water , 1971 .

[26]  Shinji Saito,et al.  Molecular dynamics simulation of the ice nucleation and growth process leading to water freezing , 2002, Nature.

[27]  Martin Zacharias,et al.  Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing‐potential , 2006, Proteins.

[28]  Xiongwu Wu,et al.  Molecular dynamics simulations of synthetic peptide folding , 1996, Proteins.

[29]  Michael W. Deem,et al.  Strict detailed balance is unnecessary in Monte Carlo simulation , 1999 .

[30]  Joan-Emma Shea,et al.  Folding Landscapes of the Alzheimer Amyloid-β(12-28) Peptide , 2006 .

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

[32]  R. L. Baldwin,et al.  A neutral, water-soluble, .alpha.-helical peptide: the effect of ionic strength on the helix-coil equilibrium , 1991 .

[33]  Car,et al.  Unified approach for molecular dynamics and density-functional theory. , 1985, Physical review letters.

[34]  Erhan Çinlar,et al.  Introduction to stochastic processes , 1974 .

[35]  A. Panagiotopoulos.  Direct determination of phase coexistence properties of fluids by Monte Carlo simulation in a new ensemble , 1987 .

[36]  Charles L. Brooks,et al.  Generalized born model with a simple smoothing function , 2003, J. Comput. Chem..

[37]  R. Levy,et al.  Protein folding pathways from replica exchange simulations and a kinetic network model. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[38]  Ernest L. Mehler,et al.  A general screened Coulomb potential based implicit solvent model: Calculation of secondary structure of small peptides , 2001 .

[39]  S. Nosé.  A molecular dynamics method for simulations in the canonical ensemble , 1984 .

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

[41]  B. Dominy,et al.  Development of a generalized Born model parameterization for proteins and nucleic acids , 1999 .

[42]  Alexander D. MacKerell,et al.  Extending the treatment of backbone energetics in protein force fields: Limitations of gas‐phase quantum mechanics in reproducing protein conformational distributions in molecular dynamics simulations , 2004, J. Comput. Chem..

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

[44]  William L. Jorgensen,et al.  Theoretical studies of medium effects on conformational equilibria , 1983 .

[45]  K. Tominaga,et al.  Dynamical interactions between solute and solvent studied by three-pulse photon echo method , 2005 .

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

[47]  Kentaro Shimizu,et al.  Folding free‐energy landscape of a 10‐residue mini‐protein, chignolin , 2006, FEBS letters.

[48]  A. Young,et al.  Spin glasses and random fields , 1997 .

[49]  James T. Hynes,et al.  A Molecular Jump Mechanism of Water Reorientation , 2006, Science.

[50]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[51]  H. C. Andersen.  Molecular dynamics simulations at constant pressure and/or temperature , 1980 .

[52]  Sergio A Hassan,et al.  Long dynamics simulations of proteins using atomistic force fields and a continuum representation of solvent effects: Calculation of structural and dynamic properties , 2005, Proteins.

[53]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .