Partial multicanonical algorithm for molecular dynamics and Monte Carlo simulations.

Partial multicanonical algorithm is proposed for molecular dynamics and Monte Carlo simulations. The partial multicanonical simulation samples a wide range of a part of the potential-energy terms, which is necessary to sample the conformational space widely, whereas a wide range of total potential energy is sampled in the multicanonical algorithm. Thus, one can concentrate the effort to determine the weight factor only on the important energy terms in the partial multicanonical simulation. The partial multicanonical, multicanonical, and canonical molecular dynamics algorithms were applied to an alanine dipeptide in explicit water solvent. The canonical simulation sampled the states of P(II), C(5), alpha(R), and alpha(P). The multicanonical simulation covered the alpha(L) state as well as these states. The partial multicanonical simulation also sampled the C(7) (ax) state in addition to the states that were sampled by the multicanonical simulation. In the partial multicanonical simulation, furthermore, backbone dihedral angles phi and psi rotated more frequently than those in the multicanonical and canonical simulations. These results mean that the partial multicanonical algorithm has a higher sampling efficiency than the multicanonical and canonical algorithms.

[1]  Peter A. Kollman,et al.  Ab initio SCF and MP2 calculations on four low-energy conformers of N-acetyl-N'-methylalaninamide , 1992 .

[2]  A. Pohorille,et al.  The development/application of a ‘minimalist’ organic/biochemical molecular mechanic force field using a combination of ab initio calculations and experimental data , 1997 .

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

[4]  Yuko Okamoto,et al.  Theoretical studies of transition states by the multioverlap molecular dynamics methods. , 2006, The Journal of chemical physics.

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

[6]  Yuko Okamoto,et al.  A new generalized-ensemble algorithm: multicanonical–multioverlap algorithm , 2007 .

[7]  Yuko Okamoto,et al.  Multibaric–multithermal ensemble molecular dynamics simulations , 2006, J. Comput. Chem..

[8]  Yuko Okamoto,et al.  Explicit symplectic integrators of molecular dynamics algorithms for rigid-body molecules in the canonical, isobaric-isothermal, and related ensembles. , 2007, The Journal of chemical physics.

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

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

[11]  Paul E. Smith,et al.  THE ALANINE DIPEPTIDE FREE ENERGY SURFACE IN SOLUTION , 1999 .

[12]  Yuko Okamoto,et al.  Monte Carlo simulations in generalized isobaric-isothermal ensembles. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  Yuko Okamoto,et al.  Multibaric–Multithermal Molecular Dynamics Simulation of Alanine Dipeptide in Explicit Water , 2007 .

[14]  Bernd A. Berg,et al.  Generalized ensemble simulations for complex systems , 2001, cond-mat/0110521.

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

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

[17]  Yuko Okamoto,et al.  Multicanonical algorithm combined with the RISM theory for simulating peptides in aqueous solution , 2000 .

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

[19]  Yuko Okamoto,et al.  Generalized-ensemble algorithms: enhanced sampling techniques for Monte Carlo and molecular dynamics simulations. , 2003, Journal of molecular graphics & modelling.

[20]  Yuko Okamoto,et al.  Generalized-ensemble algorithms for molecular dynamics simulations , 2007 .

[21]  Amedeo Caflisch,et al.  Calculation of conformational transitions and barriers in solvated systems: Application to the alanine dipeptide in water , 1999 .

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

[23]  Y. Sugita,et al.  Comparisons of force fields for proteins by generalized-ensemble simulations , 2004 .

[24]  W. L. Jorgensen,et al.  Development and Testing of the OPLS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids , 1996 .

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

[26]  Y. Yoshimoto.  Extended multicanonical method combined with thermodynamically optimized potential: application to the liquid-crystal transition of silicon. , 2006, The Journal of chemical physics.

[27]  P. Weiner,et al.  Computer Simulation of Biomolecular Systems , 1997 .

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

[29]  Berg,et al.  New approach to spin-glass simulations. , 1992, Physical review letters.

[30]  Thomas F. Miller,et al.  Symplectic quaternion scheme for biophysical molecular dynamics , 2002 .

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

[32]  S. Nosé.  A unified formulation of the constant temperature molecular dynamics methods , 1984 .

[33]  R. Friesner,et al.  Evaluation and Reparametrization of the OPLS-AA Force Field for Proteins via Comparison with Accurate Quantum Chemical Calculations on Peptides† , 2001 .

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

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

[36]  Effective sampling in the configurational space of a small peptide by the multicanonical-multioverlap algorithm. , 2007, Physical review. E, Statistical, nonlinear, and soft matter physics.

[37]  H. Okumura,et al.  Monte Carlo simulations in multibaric–multithermal ensemble , 2004 .

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

[39]  T. Kiuchi,et al.  Multi-pH Monte Carlo Simulation of Coil-Globule Transition of Weak Polyelectrolyte , 2005 .

[40]  R A Sayle,et al.  RASMOL: biomolecular graphics for all. , 1995, Trends in biochemical sciences.

[41]  Ronald M. Levy,et al.  Long-time conformational transitions of alanine dipeptide in aqueous solution: Continuous and discrete-state kinetic models , 2004 .

[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]  John D. Chodera,et al.  Long-Time Protein Folding Dynamics from Short-Time Molecular Dynamics Simulations , 2006, Multiscale Model. Simul..

[44]  Y. Okamoto,et al.  Thermodynamics of Helix-Coil Transitions Studied by Multicanonical Algorithms , 1995, chem-ph/9505006.

[45]  Alan M. Ferrenberg,et al.  New Monte Carlo technique for studying phase transitions. , 1988, Physical review letters.

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

[47]  Robert H. Swendsen,et al.  Erratum: ``New Monte Carlo technique for studying phase transitions'' [Phys. Rev. Lett. 61, 2635 (1988)] , 1989 .

[48]  Yuko Okamoto,et al.  The Generalized-Ensemble Approach for Protein Folding Simulations , 1999 .

[49]  Yuko Okamoto,et al.  Multi-overlap molecular dynamics methods for biomolecular systems , 2004 .

[50]  Yuko Okamoto,et al.  Liquid-Gas Phase Transitions Studied by Multibaric-Multithermal Monte Carlo Simulations , 2004 .

[51]  Stephen D. Bond,et al.  The Nosé-Poincaré Method for Constant Temperature Molecular Dynamics , 1999 .

[52]  S. Nosé.  An Improved Symplectic Integrator for Nosé-Poincaré Thermostat , 2001 .

[53]  Yuko Okamoto,et al.  Molecular dynamics simulations in the multibaric–multithermal ensemble , 2004 .

[54]  H. Yoshida.  Construction of higher order symplectic integrators , 1990 .

[55]  Multioverlap simulations for transitions between reference configurations. , 2003, Physical review. E, Statistical, nonlinear, and soft matter physics.

[56]  D. Landau,et al.  Efficient, multiple-range random walk algorithm to calculate the density of states. , 2000, Physical review letters.

[57]  Yuko Okamoto,et al.  Temperature and pressure dependence of alanine dipeptide studied by multibaric-multithermal molecular dynamics simulations. , 2008, The journal of physical chemistry. B.

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