Enhancing Conformation and Protonation State Sampling of Hen Egg White Lysozyme Using pH Replica Exchange Molecular Dynamics.

We evaluate the efficiency of the pH replica exchange molecular dynamics (pH-REMD) method proposed by Itoh et al. (Proteins2011, 79, 3420-3436) by using it to predict the pKa values of the titratable residues in hen egg white lysozyme (HEWL). pKa values predicted using pH-REMD converge significantly faster than those calculated using constant pH molecular dynamics (CpHMD). Furthermore, increasing the frequency between exchange attempts in pH-REMD simulations improves protonation and conformational state sampling. By enabling the simulation to sample both conformational and protonation states more rapidly, pH-REMD simulations provide valuable insight into the pH-dependence of HEWL that the CpHMD simulations failed to capture. We present an efficient and highly scalable implementation of pH-REMD as an attractive enhancement to traditional CpHMD methods.

[1]  Peter J. Artymiuk,et al.  The Structures of the Monoclinic and Orthorhombic Forms of Hen Egg-White Lysozyme at 6 Angstroms Resolution , 1981 .

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

[3]  Kay Hamacher,et al.  Information theoretical measures to analyze trajectories in rational molecular design , 2007, J. Comput. Chem..

[4]  C. Anfinsen,et al.  SOME RELATIONSHIPS OF STRUCTURE TO FUNCTION IN RIBONUCLEASE , 1959, Annals of the New York Academy of Sciences.

[5]  Daniel J. Sindhikara,et al.  Bad Seeds Sprout Perilous Dynamics: Stochastic Thermostat Induced Trajectory Synchronization in Biomolecules. , 2009, Journal of chemical theory and computation.

[6]  César Augusto F. de Oliveira,et al.  Coupling Constant pH Molecular Dynamics with Accelerated Molecular Dynamics , 2010, Journal of chemical theory and computation.

[7]  Gregory A. Grothaus,et al.  A simple clustering algorithm can be accurate enough for use in calculations of pKs in macromolecules , 2006, Proteins.

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

[9]  A. Cornish-Bowden,et al.  The pH-dependence of pepsin-catalysed reactions. , 1969, The Biochemical journal.

[10]  D. Bashford,et al.  Electrostatic calculations of the pKa values of ionizable groups in bacteriorhodopsin , 1992 .

[11]  Rebecca C. Wade,et al.  Improving the Continuum Dielectric Approach to Calculating pKas of Ionizable Groups in Proteins , 1996 .

[12]  Daniel J. Sindhikara,et al.  Exchange Often and Properly in Replica Exchange Molecular Dynamics. , 2010, Journal of chemical theory and computation.

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

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

[15]  Sarah L. Williams,et al.  Progress in the prediction of pKa values in proteins , 2011, Proteins.

[16]  A. Roitberg,et al.  pH-replica exchange molecular dynamics in proteins using a discrete protonation method. , 2012, The journal of physical chemistry. B.

[17]  M. Karplus,et al.  pKa's of ionizable groups in proteins: atomic detail from a continuum electrostatic model. , 1990, Biochemistry.

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

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

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

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

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

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

[24]  Daniel J. Sindhikara,et al.  Exchange frequency in replica exchange molecular dynamics. , 2008, The Journal of chemical physics.

[25]  Matthew P Jacobson,et al.  Comparing Conformational Ensembles Using the Kullback-Leibler Divergence Expansion. , 2012, Journal of chemical theory and computation.

[26]  E. Alexov,et al.  Incorporating protein conformational flexibility into the calculation of pH-dependent protein properties. , 1997, Biophysical journal.

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

[28]  A. Baptista,et al.  Is the prediction of pKa values by constant‐pH molecular dynamics being hindered by inherited problems? , 2011, Proteins.

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

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

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

[32]  Michael R. Shirts,et al.  Replica exchange and expanded ensemble simulations as Gibbs sampling: simple improvements for enhanced mixing. , 2011, The Journal of chemical physics.

[33]  William Swope,et al.  Understanding folding and design: Replica-exchange simulations of ``Trp-cage'' miniproteins , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[34]  Jana K. Shen,et al.  Continuous Constant pH Molecular Dynamics in Explicit Solvent with pH-Based Replica Exchange. , 2011, Journal of chemical theory and computation.

[35]  Jan H. Jensen,et al.  PROPKA3: Consistent Treatment of Internal and Surface Residues in Empirical pKa Predictions. , 2011, Journal of chemical theory and computation.

[36]  Ulrich H E Hansmann,et al.  Folding proteins by first-passage-times-optimized replica exchange. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[37]  Danial S. Dashti,et al.  Computing Alchemical Free Energy Differences with Hamiltonian Replica Exchange Molecular Dynamics (H-REMD) Simulations. , 2011, Journal of chemical theory and computation.

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

[39]  Bernard R Brooks,et al.  pH replica‐exchange method based on discrete protonation states , 2011, Proteins.

[40]  C. Tanford,et al.  Theory of Protein Titration Curves. I. General Equations for Impenetrable Spheres , 1957 .

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

[42]  Karsten Schwan,et al.  Keeneland: Bringing Heterogeneous GPU Computing to the Computational Science Community , 2011, Computing in Science & Engineering.

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