Multicanonical ab inito QM/MM molecular dynamics simulation of a peptide in an aqueous environment

We developed a multicanonical ab initio QM/MM molecular dynamics simulation method to enhance conformational sampling of biomolecules in an aqueous environment. We applied this method to an alanine dipeptide immersed in a sphere of explicit water molecules. The peptide and the water molecules were treated by the QM method at the HF/3‐21G level and by the MM method, respectively. The van der Waals interactions between the peptide and the water molecules were calculated at the MM level, while the electrostatic interaction terms between them were incorporated into the QM Hamiltonian to account for the effect of the solvent on the electronic structure of the peptide. The simulation was performed for 1 ns, and a free‐energy map was calculated with respect to the peptide conformation. All the conformations (C5, PII, C7eq, and αR) that have been experimentally suggested to exist in solution formed basins on the free‐energy surface. Analysis of the water distribution revealed that the αR conformation was stabilized by the interaction between the large electric dipole moments of this peptide conformation and the water electric dipole moments, whereas the PII conformation was stabilized by the formation of characteristic hydrogen bonds with the water molecules. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010

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

[2]  J. Hermans,et al.  Comparison of a QM/MM force field and molecular mechanics force fields in simulations of alanine and glycine “dipeptides” (Ace‐Ala‐Nme and Ace‐Gly‐Nme) in water in relation to the problem of modeling the unfolded peptide backbone in solution , 2003, Proteins.

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

[4]  A. Kidera,et al.  Generalized form of the conserved quantity in constant-temperature molecular dynamics , 2002 .

[5]  V. Barone,et al.  Quantum Calculation of Molecular Energies and Energy Gradients in Solution by a Conductor Solvent Model , 1998 .

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

[7]  C. Breneman,et al.  Determining atom‐centered monopoles from molecular electrostatic potentials. The need for high sampling density in formamide conformational analysis , 1990 .

[8]  Takashi Imai,et al.  Temperature and pressure effects on conformational equilibria of alanine dipeptide in aqueous solution. , 2004, Biopolymers.

[9]  L. Onsager.  Electric Moments of Molecules in Liquids , 1936 .

[10]  Alexander D. MacKerell,et al.  Improved treatment of the protein backbone in empirical force fields. , 2004, Journal of the American Chemical Society.

[11]  J. Thar,et al.  When is a molecule properly solvated by a continuum model or in a cluster ansatz? A first-principles simulation of alanine hydration. , 2008, The journal of physical chemistry. B.

[12]  Fei Zhang.  Operator-splitting integrators for constant-temperature molecular dynamics , 1997 .

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

[14]  L-alanine in a droplet of water: a density-functional molecular dynamics study. , 2007, The journal of physical chemistry. B.

[15]  Tohru Terada,et al.  A method for evaluating multicanonical potential function without iterative refinement: Application to conformational sampling of a globular protein in water , 2003 .

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

[17]  Giovanni Scalmani,et al.  Energies, structures, and electronic properties of molecules in solution with the C‐PCM solvation model , 2003, J. Comput. Chem..

[18]  James C. Weisshaar,et al.  Conformational Analysis of Alanine Dipeptide from Dipolar Couplings in a Water-Based Liquid Crystal , 2003 .

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

[20]  Haruki Nakamura,et al.  Peptide free‐energy profile is strongly dependent on the force field: Comparison of C96 and AMBER95 , 2000 .

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

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

[23]  Christoph Weise,et al.  Do Bridging Water Molecules Dictate the Structure of a Model Dipeptide in Aqueous Solution , 2000 .

[24]  J. Pople,et al.  Self-consistent molecular orbital methods. 21. Small split-valence basis sets for first-row elements , 2002 .

[25]  Matthew T. Eddy,et al.  Structure of the Alanine Dipeptide in Condensed Phases Determined by 13C NMR , 2004 .

[26]  Thom Vreven,et al.  Combining Quantum Mechanics Methods with Molecular Mechanics Methods in ONIOM. , 2006, Journal of chemical theory and computation.

[27]  Yong Duan,et al.  Solvation effects on alanine dipeptide: A MP2/cc‐pVTZ//MP2/6‐31G** study of (Φ, Ψ) energy maps and conformers in the gas phase, ether, and water , 2004, J. Comput. Chem..

[28]  Minhaeng Cho,et al.  Classical and quantum mechanical/molecular mechanical molecular dynamics simulations of alanine dipeptide in water: comparisons with IR and vibrational circular dichroism spectra. , 2008, The Journal of chemical physics.

[29]  T. Terada,et al.  A multicanonical ab initio molecular dynamics method: Application to conformation sampling of alanine tripeptide , 2006 .

[30]  Sándor Suhai,et al.  Theoretical Study of Aqueous N-Acetyl-l-alanine N‘-Methylamide: Structures and Raman, VCD, and ROA Spectra , 1998 .

[31]  K. Kopple,et al.  Solvent-dependent conformational distributions of some dipeptides , 1980 .

[32]  Yuto Komeiji,et al.  Change in a protein's electronic structure induced by an explicit solvent: An ab initio fragment molecular orbital study of ubiquitin , 2007, J. Comput. Chem..

[33]  Simona Golic Grdadolnik,et al.  Determination of conformational preferences of dipeptides using vibrational spectroscopy. , 2008, The journal of physical chemistry. B.

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

[35]  Ross C Walker,et al.  Implementation of the SCC-DFTB method for hybrid QM/MM simulations within the amber molecular dynamics package. , 2007, The journal of physical chemistry. A.

[36]  V. Hornak,et al.  Comparison of multiple Amber force fields and development of improved protein backbone parameters , 2006, Proteins.

[37]  Matteo Dal Peraro,et al.  Solute-solvent charge transfer in aqueous solution. , 2005, Chemphyschem : a European journal of chemical physics and physical chemistry.