Continuous Constant pH Molecular Dynamics in Explicit Solvent with pH-Based Replica Exchange.

A computational tool that offers accurate pKa values and atomically detailed knowledge of protonation-coupled conformational dynamics is valuable for elucidating mechanisms of energy transduction processes in biology, such as enzyme catalysis and electron transfer as well as proton and drug transport. Toward this goal we present a new technique of embedding continuous constant pH molecular dynamics within an explicit-solvent representation. In this technique we make use of the efficiency of the generalized-Born (GB) implicit-solvent model for estimating the free energy of protein solvation while propagating conformational dynamics using the more accurate explicit-solvent model. Also, we employ a pH-based replica exchange scheme to significantly enhance both protonation and conformational state sampling. Benchmark data of five proteins including HP36, NTL9, BBL, HEWL, and SNase yield an average absolute deviation of 0.53 and a root mean squared deviation of 0.74 from experimental data. This level of accuracy is obtained with 1 ns simulations per replica. Detailed analysis reveals that explicit-solvent sampling provides increased accuracy relative to the previous GB-based method by preserving the native structure, providing a more realistic description of conformational flexibility of the hydrophobic cluster, and correctly modeling solvent mediated ion-pair interactions. Thus, we anticipate that the new technique will emerge as a practical tool to capture ionization equilibria while enabling an intimate view of ionization coupled conformational dynamics that is difficult to delineate with experimental techniques alone.

[1]  P. Hünenberger,et al.  Explicit-solvent molecular dynamics simulation at constant pH: Methodology and application to small amines , 2001 .

[2]  Asim Okur,et al.  Improved Efficiency of Replica Exchange Simulations through Use of a Hybrid Explicit/Implicit Solvation Model. , 2006, Journal of chemical theory and computation.

[3]  B. García-Moreno E.,et al.  Stabilization of internal charges in a protein: water penetration or conformational change? , 2004, Biophysical journal.

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

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

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

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

[8]  B. Kuhlman,et al.  pKa values and the pH dependent stability of the N-terminal domain of L9 as probes of electrostatic interactions in the denatured state. Differentiation between local and nonlocal interactions. , 1999, Biochemistry.

[9]  A. Warshel.  Calculations of enzymatic reactions: calculations of pKa, proton transfer reactions, and general acid catalysis reactions in enzymes. , 1981, Biochemistry.

[10]  Ulrich H E Hansmann,et al.  Generalized ensemble and tempering simulations: a unified view. , 2006, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  Angel E García,et al.  Molecular dynamics study of water penetration in staphylococcal nuclease , 2005, Proteins.

[12]  Charles L. Brooks,et al.  λ‐dynamics: A new approach to free energy calculations , 1996 .

[13]  P. Hünenberger,et al.  pH-Dependent Stability of a Decalysine α-Helix Studied by Explicit-Solvent Molecular Dynamics Simulations at Constant pH , 2004 .

[14]  Charles L Brooks,et al.  Exploring atomistic details of pH-dependent peptide folding , 2006, Proceedings of the National Academy of Sciences.

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

[16]  Lauren Wickstrom,et al.  Secondary structure bias in generalized Born solvent models: comparison of conformational ensembles and free energy of solvent polarization from explicit and implicit solvation. , 2007, The journal of physical chemistry. B.

[17]  R. Norel,et al.  Electrostatic aspects of protein-protein interactions. , 2000, Current opinion in structural biology.

[18]  Gerhard Hummer,et al.  Free Energy of Ionic Hydration , 1996 .

[19]  C. Soares,et al.  Constant-pH molecular dynamics using stochastic titration , 2002 .

[20]  Jana Khandogin.  MODELING PROTONATION EQUILIBRIA IN BIOLOGICAL MACROMOLECULES , 2009 .

[21]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[22]  Charles L. Brooks,et al.  New analytic approximation to the standard molecular volume definition and its application to generalized Born calculations , 2003, J. Comput. Chem..

[23]  J. Mccammon,et al.  Calculating pKa values in enzyme active sites , 2003, Protein science : a publication of the Protein Society.

[24]  A. Warshel.  Electrostatic basis of structure-function correlation in proteins , 1981 .

[25]  Charles Tanford,et al.  [84] Examination of titration behavior , 1967 .

[26]  C. Pace,et al.  pK values of the ionizable groups of proteins , 2006, Protein science : a publication of the Protein Society.

[27]  Charles L Brooks,et al.  Linking folding with aggregation in Alzheimer's β-amyloid peptides , 2007, Proceedings of the National Academy of Sciences.

[28]  Bernard R. Brooks,et al.  Removal of pressure and free energy artifacts in charged periodic systems via net charge corrections to the Ewald potential , 1998 .

[29]  Charles L Brooks,et al.  Folding intermediate in the villin headpiece domain arises from disruption of a N-terminal hydrogen-bonded network. , 2007, Journal of the American Chemical Society.

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

[31]  Jana K. Shen,et al.  Uncovering specific electrostatic interactions in the denatured states of proteins. , 2010, Biophysical journal.

[32]  Jana K. Shen,et al.  Predicting pKa values with continuous constant pH molecular dynamics. , 2009, Methods in enzymology.

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

[34]  Asim Okur,et al.  Evaluation of Salt Bridge Structure and Energetics in Peptides Using Explicit, Implicit, and Hybrid Solvation Models. , 2008, Journal of chemical theory and computation.

[35]  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..

[36]  P. Beroza,et al.  Application of a pairwise generalized Born model to proteins and nucleic acids: inclusion of salt effects , 1999 .

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

[38]  B. Brooks,et al.  Constant pressure molecular dynamics simulation: The Langevin piston method , 1995 .

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

[40]  Jana K. Shen,et al.  Probing the strand orientation and registry alignment in the propagation of amyloid fibrils. , 2010, Biochemistry.

[41]  C. Castañeda,et al.  Molecular determinants of the pKa values of Asp and Glu residues in staphylococcal nuclease , 2009, Proteins.

[42]  P E Wright,et al.  Electrostatic calculations of side-chain pK(a) values in myoglobin and comparison with NMR data for histidines. , 1993, Biochemistry.

[43]  Damien Farrell,et al.  Remeasuring HEWL pKa values by NMR spectroscopy: Methods, analysis, accuracy, and implications for theoretical pKa calculations , 2011, Proteins.

[44]  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.

[45]  D. Beglov,et al.  Atomic Radii for Continuum Electrostatics Calculations Based on Molecular Dynamics Free Energy Simulations , 1997 .

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

[47]  Vincent A Voelz,et al.  Unfolded-state dynamics and structure of protein L characterized by simulation and experiment. , 2010, Journal of the American Chemical Society.

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

[49]  Charles L Brooks,et al.  Implicit modeling of nonpolar solvation for simulating protein folding and conformational transitions. , 2008, Physical chemistry chemical physics : PCCP.

[50]  A. Fersht,et al.  Downhill versus barrier-limited folding of BBL 1: energetic and structural perturbation effects upon protonation of a histidine of unusually low pKa. , 2009, Journal of molecular biology.

[51]  W. Nadler,et al.  Optimized explicit-solvent replica exchange molecular dynamics from scratch. , 2008, The journal of physical chemistry. B.

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

[53]  V. Hornak,et al.  Investigation of Salt Bridge Stability in a Generalized Born Solvent Model. , 2006, Journal of chemical theory and computation.

[54]  Adrian E Roitberg,et al.  Constant pH replica exchange molecular dynamics in biomolecules using a discrete protonation model. , 2010, Journal of chemical theory and computation.

[55]  Michael Feig,et al.  MMTSB Tool Set: enhanced sampling and multiscale modeling methods for applications in structural biology. , 2004, Journal of molecular graphics & modelling.

[56]  Jana K. Shen.  A method to determine residue-specific unfolded-state pKa values from analysis of stability changes in single mutant cycles. , 2010, Journal of the American Chemical Society.