Introducing titratable water to all-atom molecular dynamics at constant pH.

Recent development of titratable coions has paved the way for realizing all-atom molecular dynamics at constant pH. To further improve physical realism, here we describe a technique in which proton titration of the solute is directly coupled to the interconversion between water and hydroxide or hydronium. We test the new method in replica-exchange continuous constant pH molecular dynamics simulations of three proteins, HP36, BBL, and HEWL. The calculated pKa values based on 10-ns sampling per replica have the average absolute and root-mean-square errors of 0.7 and 0.9 pH units, respectively. Introducing titratable water in molecular dynamics offers a means to model proton exchange between solute and solvent, thus opening a door to gaining new insights into the intricate details of biological phenomena involving proton translocation.

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

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

[3]  Mikael Lund,et al.  Ion pairing as a possible clue for discriminating between sodium and potassium in biological and other complex environments. , 2007, The journal of physical chemistry. B.

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

[5]  G. Voth,et al.  Structure and dynamics of hydronium in the ion channel gramicidin A. , 1996, Biophysical journal.

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

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

[8]  Benoît Roux,et al.  Control of ion selectivity in LeuT: two Na+ binding sites with two different mechanisms. , 2008, Journal of molecular biology.

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

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

[11]  Jessica M J Swanson,et al.  Proton solvation and transport in aqueous and biomolecular systems: insights from computer simulations. , 2007, The journal of physical chemistry. B.

[12]  Emil Alexov,et al.  Role of the protein side‐chain fluctuations on the strength of pair‐wise electrostatic interactions: Comparing experimental with computed pKas , 2002, Proteins.

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

[14]  Yuan Bi.  Studies of the Folding and Stability of the Villin Headpiece Subdomain , 2008 .

[15]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[16]  Arieh Warshel,et al.  Studies of proton translocations in biological systems: simulating proton transport in carbonic anhydrase by EVB-based models. , 2004, Biophysical journal.

[17]  Jennifer L. Knight,et al.  Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent. , 2012, Journal of chemical theory and computation.

[18]  Alexander D. MacKerell,et al.  Polarizable empirical force field for alkanes based on the classical Drude oscillator model. , 2005, The journal of physical chemistry. B.

[19]  H. Grubmüller,et al.  Constant pH Molecular Dynamics in Explicit Solvent with λ-Dynamics , 2011, Journal of chemical theory and computation.

[20]  Carboxyl pK(a) values and acid denaturation of BBL. , 2010, Journal of molecular biology.

[21]  Jana K. Shen,et al.  Charge-leveling and proper treatment of long-range electrostatics in all-atom molecular dynamics at constant pH. , 2012, The Journal of chemical physics.

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

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

[24]  B. Roux,et al.  Simulation of Osmotic Pressure in Concentrated Aqueous Salt Solutions , 2010 .

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

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