Repeated‐annealing sampling combined with multicanonical algorithm for conformational sampling of bio‐molecules

A novel conformational sampling method (repeated‐annealing sampling method) is proposed to execute an efficient conformational sampling at a reasonable computational cost. In the method, a molecular dynamics simulation is done with repeating an elemental process. An elemental process consists of four subprocesses: high‐temperature run, annealing, room‐temperature run, and fast heating. The sampling is done automatically according to a temperature‐control schedule. The room‐temperature run is treated with the multicanonical algorithm, and the other subprocesses are done with the conventional molecular dynamics algorithm. The method, differing from the generalized ensemble methods recently developed, is not warrantable to give the canonical ensemble because of the nonphysical process in the annealing. However, we observed that the slower the annealing and the longer the high‐temperature run, the closer the sampled conformations to those of the canonical ensemble. A test was performed with tri‐N‐acetyl‐D‐glucosamine in vacuo, and the results were compared with those from the conventional multicanonical simulation. Not only the reweighted canonical distribution function but also the energy landscape were in good agreement with those from the conventional multicanonical simulation. The potential of mean force also showed a fairly good agreement with that from the conventional multicanonical simulation in the room‐temperature region. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1098–1106, 2001

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

[2]  Straub,et al.  Generalized simulated annealing algorithms using Tsallis statistics: Application to conformational optimization of a tetrapeptide. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[3]  A. Kidera,et al.  Enhanced conformational sampling in Monte Carlo simulations of proteins: application to a constrained peptide. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

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

[5]  Robert J. Woods,et al.  Molecular Mechanical and Molecular Dynamic Simulations of Glycoproteins and Oligosaccharides. 1. GLYCAM_93 Parameter Development , 1995 .

[6]  Yuji Sugita,et al.  Replica-exchange multicanonical algorithm and multicanonical replica-exchange method for simulating systems with rough energy landscape , 2000, cond-mat/0009119.

[7]  Christian Bartels,et al.  Determination of equilibrium properties of biomolecular systems using multidimensional adaptive umbrella sampling , 1999 .

[8]  N. Nakajima,et al.  A selectively enhanced multicanonical molecular dynamics method for conformational sampling of peptides in realistic water molecules , 1998 .

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

[10]  Lee,et al.  New Monte Carlo algorithm: Entropic sampling. , 1993, Physical review letters.

[11]  George Chikenji,et al.  Simulation of Lattice Polymers with Multi-Self-Overlap Ensemble , 1998 .

[12]  A. Kidera,et al.  Free energy landscapes of peptides by enhanced conformational sampling. , 2000, Journal of molecular biology.

[13]  Y. Okamoto,et al.  GENERALIZED-ENSEMBLE MONTE CARLO METHOD FOR SYSTEMS WITH ROUGH ENERGY LANDSCAPE , 1997, cond-mat/9710306.

[14]  Christian Bartels,et al.  Multidimensional adaptive umbrella sampling: Applications to main chain and side chain peptide conformations , 1997 .

[15]  Haruki Nakamura,et al.  Presto(protein Engineering Simulator): A Vectorized Molecular Mechanics Program for Biopolymers , 1992, Comput. Chem..

[16]  N. Nakajima,et al.  Enhanced conformational diversity search of CDR‐H3 in antibodies: Role of the first CDR‐H3 residue , 1999, Proteins.

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

[18]  Haruki Nakamura,et al.  The multicanonical weighted histogram analysis method for the free-energy landscape along structural transition paths , 1999 .

[19]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .

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

[21]  M. Karplus,et al.  Solution conformations and thermodynamics of structured peptides: molecular dynamics simulation with an implicit solvation model. , 1998, Journal of molecular biology.

[22]  Yuko Okamoto,et al.  Prediction of peptide conformation by multicanonical algorithm: New approach to the multiple‐minima problem , 1993, J. Comput. Chem..

[23]  Haruki Nakamura,et al.  Flexible docking of a ligand peptide to a receptor protein by multicanonical molecular dynamics simulation , 1997 .

[24]  R. Dwek,et al.  The high degree of internal flexibility observed for an oligomannose oligosaccharide does not alter the overall topology of the molecule. , 1998, European journal of biochemistry.

[25]  Haruki Nakamura,et al.  Two‐component multicanonical Monte Carlo method for effective conformation sampling , 1997 .

[26]  George Chikenji,et al.  MULTI-SELF-OVERLAP ENSEMBLE FOR PROTEIN FOLDING : GROUND STATE SEARCH AND THERMODYNAMICS , 1999, cond-mat/9903003.

[27]  Yuko Okamoto,et al.  An Analysis on Protein Folding Problem by Replica-Exchange Method , 2000 .

[28]  K. Watanabe,et al.  Dissection of protein-carbohydrate interactions in mutant hen egg-white lysozyme complexes and their hydrolytic activity. , 1995, Journal of molecular biology.

[29]  Hesselbo,et al.  Monte Carlo simulation and global optimization without parameters. , 1995, Physical review letters.

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

[31]  Martin Karplus,et al.  Probability Distributions for Complex Systems: Adaptive Umbrella Sampling of the Potential Energy , 1998 .

[32]  T. Yamato,et al.  Molecular Dynamics of a 15-Residue Poly(l-alanine) in Water: Helix Formation and Energetics , 1999 .

[33]  Haruki Nakamura,et al.  β-hairpin folds by molecular dynamics simulations , 2000 .

[34]  William G. Hoover,et al.  High-strain-rate plastic flow studied via nonequilibrium molecular dynamics , 1982 .

[35]  C. Tsallis Possible generalization of Boltzmann-Gibbs statistics , 1988 .

[36]  A. Kidera,et al.  Conformational sampling of CDR-H3 in antibodies by multicanonical molecular dynamics simulation. , 1998, Journal of molecular biology.

[37]  M Karplus,et al.  Characterization of flexible molecules in solution: the RGDW peptide. , 1998, Journal of molecular biology.

[38]  Gary P. Morriss,et al.  The isothermal/isobaric molecular dynamics ensemble , 1983 .

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

[40]  Y. Okamoto,et al.  Numerical comparisons of three recently proposed algorithms in the protein folding problem , 1997 .