A Novel Hamiltonian Replica Exchange MD Protocol to Enhance Protein Conformational Space Sampling.

Limited searching in the conformational space is one of the major obstacles for investigating protein dynamics by numerical approaches. For this reason, classical all-atom molecular dynamics (MD) simulations of proteins tend to be confined to local energy minima, particularly when the bulk solvent is treated explicitly. To overcome this problem, we have developed a novel replica exchange protocol that uses modified force-field parameters to treat interparticle nonbonded potentials within the protein and between protein and solvent atoms, leaving unperturbed those relative to solvent-solvent interactions. We have tested the new protocol on the 18-residue-long tip of the P domain of calreticulin in an explicit solvent. With only eight replicas, we have been able to considerably enhance the conformational space sampled during a 100 ns simulation, compared to as many parallel classical molecular dynamics simulations of the same length or to a single one lasting 450 ns. A direct comparison between the various simulations has been possible thanks to the implementation of the weighted histogram analysis method, by which conformations simulated with modified force-field parameters can be assigned different weights. Interatom, inter-residue distances in the structural ensembles obtained with our novel replica exchange approach and by classical MD simulations compare equally well with those derived from NMR data. Rare events, such as unfolding and refolding, occur with reasonable statistical frequency. Visiting of conformations characterized by very small Boltzmann weights is also possible. Despite their low probability, such regions of the conformational space may play an important role in the search for local potential-energy minima and in dynamically controlled functions.

[1]  Y. Sugita,et al.  MOLECULAR DYNAMICS SIMULATIONS OF DNA DIMERS BASED ON REPLICA-EXCHANGE UMBRELLA SAMPLING II: FREE ENERGY ANALYSIS , 2005 .

[2]  Eric J. Sorin,et al.  How well can simulation predict protein folding kinetics and thermodynamics? , 2005, Annual review of biophysics and biomolecular structure.

[3]  A. Schug,et al.  Reproducible protein folding with the stochastic tunneling method. , 2003, Physical review letters.

[4]  Ivano Tavernelli,et al.  Protein dynamics, thermal stability, and free-energy landscapes: a molecular dynamics investigation. , 2003, Biophysical journal.

[5]  Seokmin Shin,et al.  Replica-exchange method using the generalized effective potential. , 2003, Physical review letters.

[6]  K. Wüthrich,et al.  NMR structures of 36 and 73-residue fragments of the calreticulin P-domain. , 2002, Journal of molecular biology.

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

[8]  I. Tavernelli,et al.  The interplay between protein dynamics and frustration of non-bonded interactions as revealed by molecular dynamics simulations , 2001 .

[9]  Berk Hess,et al.  GROMACS 3.0: a package for molecular simulation and trajectory analysis , 2001 .

[10]  X. Daura,et al.  The beta-peptide hairpin in solution: conformational study of a beta-hexapeptide in methanol by NMR spectroscopy and MD simulation. , 2001, Journal of the American Chemical Society.

[11]  R. Riek,et al.  NMR structure of the calreticulin P-domain , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[12]  Yuko Okamoto,et al.  Replica-exchange Monte Carlo method for the isobaric isothermal ensemble , 2001 .

[13]  R. Riek,et al.  Three‐dimensional structure topology of the calreticulin P‐domain based on NMR assignment , 2001, FEBS letters.

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

[15]  Ilya A. Balabin,et al.  Dynamically controlled protein tunneling paths in photosynthetic reaction centers. , 2000, Science.

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

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

[18]  J. Skolnick,et al.  Comparison of three Monte Carlo conformational search strategies for a proteinlike homopolymer model: Folding thermodynamics and identification of low-energy structures , 2000 .

[19]  Yamamoto,et al.  Replica-exchange molecular dynamics simulation for supercooled liquids , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

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

[21]  P. Wolynes,et al.  The physics of protein folding , 1999 .

[22]  K. Hukushima,et al.  Domain-wall free energy of spin-glass models: numerical method and boundary conditions. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[23]  X. Daura,et al.  Peptide Folding: When Simulation Meets Experiment , 1999 .

[24]  Berk Hess,et al.  LINCS: A linear constraint solver for molecular simulations , 1997, J. Comput. Chem..

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

[26]  Shankar Kumar,et al.  Multidimensional free‐energy calculations using the weighted histogram analysis method , 1995, J. Comput. Chem..

[27]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[28]  Wilfred F. van Gunsteren,et al.  A generalized reaction field method for molecular dynamics simulations , 1995 .

[29]  J. Onuchic,et al.  Funnels, pathways, and the energy landscape of protein folding: A synthesis , 1994, Proteins.

[30]  Wilfred F. van Gunsteren,et al.  Consistent dielectric properties of the simple point charge and extended simple point charge water models at 277 and 300 K , 1994 .

[31]  C. Brooks,et al.  Constant-temperature free energy surfaces for physical and chemical processes , 1993 .

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

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

[34]  H. Scheraga,et al.  Application of the diffusion equation method for global optimization to oligopeptides , 1992 .

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

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

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

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

[39]  H. Scheraga,et al.  On the multiple-minima problem in the conformational analysis of molecules: deformation of the potential energy hypersurface by the diffusion equation method , 1989 .

[40]  T. Straatsma,et al.  THE MISSING TERM IN EFFECTIVE PAIR POTENTIALS , 1987 .

[41]  H. Berendsen,et al.  Molecular dynamics with coupling to an external bath , 1984 .

[42]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[43]  L. Ellgaard,et al.  Calnexin, calreticulin, and ERp57 , 2007, Cell Biochemistry and Biophysics.

[44]  Garland R. Marshall,et al.  A potential smoothing algorithm accurately predicts transmembrane helix packing , 1999, Nature Structural Biology.