Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics

An extension of the constant pH method originally implemented by Mongan et al. (J. Comput. Chem.2004, 25, 2038−2048) is proposed in this study. This adapted version of the method couples the constant pH methodology with the enhanced sampling technique of accelerated molecular dynamics, in an attempt to overcome the sampling issues encountered with current standard constant pH molecular dynamics methods. Although good results were reported by Mongan et al. on application of the standard method to the hen egg-white lysozyme (HEWL) system, residues which possess strong interactions with neighboring groups tend to converge slowly, resulting in the reported inconsistencies for predicted pKa values, as highlighted by the authors. The application of the coupled method described in this study to the HEWL system displays improvements over the standard version of the method, with the improved sampling leading to faster convergence and producing pKa values in closer agreement to those obtained experimentally for the more slowly converging residues.

[1]  D. Case,et al.  Exploring protein native states and large‐scale conformational changes with a modified generalized born model , 2004, Proteins.

[2]  C. Brooks,et al.  Constant‐pH molecular dynamics using continuous titration coordinates , 2004, Proteins.

[3]  R. Rosenfeld.  Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[4]  M. Gilson,et al.  Prediction of pH-dependent properties of proteins. , 1994, Journal of molecular biology.

[5]  P. Kollman,et al.  Simulating proteins at constant pH: An approach combining molecular dynamics and Monte Carlo simulation , 2002, Proteins.

[6]  Tamar Schlick,et al.  New Algorithms for Macromolecular Simulation , 2006 .

[7]  S. Withers,et al.  Catalysis by hen egg-white lysozyme proceeds via a covalent intermediate , 2001, Nature.

[8]  Haruki Nakamura,et al.  Electrostatic forces in two lysozymes: Calculations and measurements of histidine pKa values , 1992, Biopolymers.

[9]  H. Stern.  Molecular simulation with variable protonation states at constant pH. , 2007, The Journal of chemical physics.

[10]  M Karplus,et al.  Improving the accuracy of protein pKa calculations: Conformational averaging versus the average structure , 1998, Proteins.

[11]  Charles L Brooks,et al.  Toward the accurate first-principles prediction of ionization equilibria in proteins. , 2006, Biochemistry.

[12]  John Mongan,et al.  Biomolecular simulations at constant pH. , 2005, Current opinion in structural biology.

[13]  Miguel Machuqueiro,et al.  Acidic range titration of HEWL using a constant‐pH molecular dynamics method , 2008, Proteins.

[14]  Charles L. Brooks,et al.  Generalized born model with a simple smoothing function , 2003, J. Comput. Chem..

[15]  F. Young.  Biochemistry , 1955, The Indian Medical Gazette.

[16]  D. Case,et al.  Constant pH molecular dynamics in generalized Born implicit solvent , 2004, J. Comput. Chem..

[17]  S. Subramaniam,et al.  Explicit solvent models in protein pKa calculations. , 1996, Biophysical journal.

[18]  Donald Hamelberg,et al.  Estimating kinetic rates from accelerated molecular dynamics simulations: alanine dipeptide in explicit solvent as a case study. , 2007, The Journal of chemical physics.

[19]  Charles L Brooks,et al.  Recent advances in implicit solvent-based methods for biomolecular simulations. , 2008, Current opinion in structural biology.

[20]  C Redfield,et al.  Measurement of the individual pKa values of acidic residues of hen and turkey lysozymes by two-dimensional 1H NMR. , 1994, Biophysical journal.

[21]  C. Brooks,et al.  Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field. , 2006, Journal of the American Chemical Society.

[22]  J. Antosiewicz,et al.  Constant-pH molecular dynamics study of protonation-structure relationship in a heptapeptide derived from ovomucoid third domain. , 2004, Physical review. E, Statistical, nonlinear, and soft matter physics.

[23]  C. Tanford,et al.  Interpretation of protein titration curves. Application to lysozyme. , 1972, Biochemistry.

[24]  P. Nguyen,et al.  Energy landscape of a small peptide revealed by dihedral angle principal component analysis , 2004, Proteins.

[25]  Donald Bashford,et al.  Macroscopic electrostatic models for protonation states in proteins. , 2004, Frontiers in bioscience : a journal and virtual library.

[26]  J. Mongan,et al.  Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules. , 2004, The Journal of chemical physics.

[27]  S. Petersen,et al.  Simulation of protein conformational freedom as a function of pH: constant‐pH molecular dynamics using implicit titration , 1997, Proteins.

[28]  A. Yeh,et al.  Rapid and simple protein-stability screens: application to membrane proteins. , 2006, Acta crystallographica. Section D, Biological crystallography.

[29]  C. Brooks,et al.  Constant pH molecular dynamics with proton tautomerism. , 2005, Biophysical journal.

[30]  Gerrit Groenhof,et al.  GROMACS: Fast, flexible, and free , 2005, J. Comput. Chem..

[31]  B Honig,et al.  On the pH dependence of protein stability. , 1993, Journal of molecular biology.

[32]  David A. Case,et al.  Effective Born radii in the generalized Born approximation: The importance of being perfect , 2002, J. Comput. Chem..

[33]  E. Alexov,et al.  Combining conformational flexibility and continuum electrostatics for calculating pK(a)s in proteins. , 2002, Biophysical journal.