Combination of molecular dynamics method and 3D‐RISM theory for conformational sampling of large flexible molecules in solution

We have developed an algorithm for sampling the conformational space of large flexible molecules in solution, which combines the molecular dynamics (MD) method and the three‐dimensional reference interaction site model (3D‐RISM) theory. The solvent‐induced force acting on solute atoms was evaluated as the gradient of the solvation free energy with respect to the solute‐atom coordinates. To enhance the computation speed, we have applied a multiple timestep algorithm based on the RESPA (Reversible System Propagator Algorithm) to the combined MD/3D‐RISM method. By virtue of the algorithm, one can choose a longer timestep for renewing the solvent‐induced force compared with that of the conformational update. To illustrate the present MD/3D‐RISM simulation, we applied the method to a model of acetylacetone in aqueous solution. The multiple timestep algorithm succeeded in enhancing the computation speed by 3.4 times for this model case. Acetylacetone possesses an intramolecular hydrogen‐bonding capability between the hydroxyl group and the carbonyl oxygen atom, and the molecule is significantly stabilized due to this hydrogen bond, especially in gas phase. The intramolecular hydrogen bond was kept intact during almost entire course of the MD simulation in gas phase, while in the aqueous solutions the bond is disrupted in a significant number of conformations. This result qualitatively agrees with the behavior on a free energy barrier lying upon the process for rotating a torsional degree of freedom of the hydroxyl group, where it is significantly reduced in aqueous solution by a cancellation between the electrostatic interaction and the solvation free energy. © 2007 Wiley Periodicals, Inc. J Comput Chem 2008

[1]  Fumio Hirata,et al.  Water molecules in a protein cavity detected by a statistical-mechanical theory. , 2005, Journal of the American Chemical Society.

[2]  Fumio Hirata,et al.  Solvation thermodynamics of protein studied by the 3D-RISM theory , 2004 .

[3]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .

[4]  V. Pande,et al.  Multiplexed-replica exchange molecular dynamics method for protein folding simulation. , 2003, Biophysical journal.

[5]  Fumio Hirata,et al.  A new method to determine electrostatic potential around a macromolecule in solution from molecular wave functions , 2006, J. Comput. Chem..

[6]  B. Montgomery Pettitt,et al.  Application of an extended RISM equation to dipolar and quadrupolar fluids , 1982 .

[7]  B. Montgomery Pettitt,et al.  Integral equation predictions of liquid state structure for waterlike intermolecular potentials , 1982 .

[8]  Yuko Okamoto,et al.  Effects of Side-Chain Charges on α-Helix Stability in C-Peptide of Ribonuclease A Studied by Multicanonical Algorithm , 1999 .

[9]  Anthony K. Felts,et al.  Free Energy Surfaces of -Hairpin and -Helical Peptides Generated by Replica Exchange Molecular Dynamics with the AGBNP Implicit Solvent Model , 2004 .

[10]  Alexander D. MacKerell,et al.  An all-atom empirical energy function for the simulation of nucleic acids , 1995 .

[11]  Yuko Okamoto,et al.  Tertiary Structure Prediction of C-Peptide of Ribonuclease A by Multicanonical Algorithm , 1998, physics/9806017.

[12]  David Chandler,et al.  Optimized Cluster Expansions for Classical Fluids. II. Theory of Molecular Liquids , 1972 .

[13]  C. Brooks,et al.  Comparative Study of the Folding Free Energy Landscape of a Three-Stranded β-Sheet Protein with Explicit and Implicit Solvent Models , 2000 .

[14]  P A Kollman,et al.  Molecular dynamics studies of the HIV-1 TAR and its complex with argininamide. , 2000, Nucleic acids research.

[15]  Fumio Hirata,et al.  Self-consistent description of a metal–water interface by the Kohn–Sham density functional theory and the three-dimensional reference interaction site model , 1999 .

[16]  Marcin Król,et al.  Comparison of various implicit solvent models in molecular dynamics simulations of immunoglobulin G light chain dimer , 2003, J. Comput. Chem..

[17]  Themis Lazaridis,et al.  Implicit solvent simulations of DPC micelle formation. , 2005, The journal of physical chemistry. B.

[18]  Drew Myers,et al.  Surfaces, Interfaces, and Colloids: Principles and Applications , 1991 .

[19]  Fumio Hirata,et al.  An extended rism equation for molecular polar fluids , 1981 .

[20]  Yuko Okamoto,et al.  First-Principle Determination of Peptide Conformations in Solvents: Combination of Monte Carlo Simulated Annealing and RISM Theory , 1998 .

[21]  B. Montgomery Pettitt,et al.  The interionic potential of mean force in a molecular polar solvent from an extended RISM equation , 1983 .

[22]  Noriyuki Yoshii,et al.  A molecular dynamics study of free energy of micelle formation for sodium dodecyl sulfate in water and its size distribution. , 2006, The Journal of chemical physics.

[23]  Haruki Nakamura,et al.  Conformational transition states of a β‐hairpin peptide between the ordered and disordered conformations in explicit water , 2002, Protein science : a publication of the Protein Society.

[24]  Emilio Gallicchio,et al.  Conformational equilibrium of cytochrome P450 BM-3 complexed with N-palmitoylglycine: a replica exchange molecular dynamics study. , 2006, Journal of the American Chemical Society.

[25]  B. Berne,et al.  Can a continuum solvent model reproduce the free energy landscape of a β-hairpin folding in water? , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[26]  M. Karplus,et al.  Understanding beta-hairpin formation. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Yuko Okamoto,et al.  Replica-exchange multicanonical and multicanonical replica-exchange Monte Carlo simulations of peptides. II. Application to a more complex system , 2003 .

[28]  David Chandler,et al.  Free energy functions in the extended RISM approximation , 1985 .

[29]  R. Zhou.  Free energy landscape of protein folding in water: Explicit vs. implicit solvent , 2003, Proteins.

[30]  Fumio Hirata,et al.  Structure of tert-Butyl Alcohol−Water Mixtures Studied by the RISM Theory , 2002 .

[31]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[32]  Yuko Okamoto,et al.  Combination of the Replica-Exchange Monte Carlo Method and the Reference Interaction Site Model Theory for Simulating a Peptide Molecule in Aqueous Solution , 2004 .

[33]  M. Karplus,et al.  Effective energy function for proteins in solution , 1999, Proteins.

[34]  Noriyuki Yoshii,et al.  A molecular dynamics study of surface structure of spherical SDS micelles , 2006 .

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

[36]  Atsushi Suenaga,et al.  Replica-exchange molecular dynamics simulations for a small-sized protein folding with implicit solvent , 2003 .

[37]  Norio Yoshida,et al.  Selective ion-binding by protein probed with the 3D-RISM theory. , 2006, Journal of the American Chemical Society.

[38]  Fumio Hirata,et al.  Solution of three‐dimensional reference interaction site model and hypernetted chain equations for simple point charge water by modified method of direct inversion in iterative subspace , 1999 .

[39]  J. Cordeiro,et al.  Dynamics of meso and thermo citrate synthases with implicit solvation , 2006 .

[40]  Noriyuki Yoshii,et al.  A molecular dynamics study of structural stability of spherical SDS micelle as a function of its size , 2006 .

[41]  C. Brooks.  Computer simulation of liquids , 1989 .

[42]  平田 文男.  Molecular theory of solvation , 2003 .

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

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

[45]  Yuko Okamoto,et al.  Analysis on conformational stability of C-peptide of ribonuclease A in water using the reference interaction site model theory and Monte Carlo simulated annealing , 1999 .

[46]  Konrad Hinsen,et al.  A potential function for computer simulation studies of proton transfer in acetylacetone , 1997, J. Comput. Chem..

[47]  Fumio Hirata,et al.  Analytical energy gradient for the reference interaction site model multiconfigurational self‐consistent‐field method: Application to 1,2‐difluoroethylene in aqueous solution , 1996 .

[48]  Yuko Okamoto,et al.  Effects of the fixed end in single-molecule imaging techniques: A replica-exchange molecular dynamics simulation study , 2005 .

[49]  Yuko Okamoto,et al.  Structures of a peptide fragment of β2‐microglobulin studied by replica‐exchange molecular dynamics simulations – Towards the understanding of the mechanism of amyloid formation , 2005, FEBS letters.