Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent.

The nucleosides of adenine and cytosine have pKa values of 3.50 and 4.08, respectively, and are assumed to be unprotonated under physiological conditions. However, evidence from recent NMR and X-Ray crystallography studies has revealed the prevalence of protonated adenine and cytosine in RNA macromolecules. Such nucleotides with elevated pKa values may play a role in stabilizing RNA structure and participate in the mechanism of ribozyme catalysis. With the work presented here, we establish the framework and demonstrate the first constant pH MD simulations (CPHMD) for nucleic acids in explicit solvent in which the protonation state is coupled to the dynamical evolution of the RNA system via λ-dynamics. We adopt the new functional form λ(Nexp) for λ that was recently developed for Multi-Site λ-Dynamics (MSλD) and demonstrate good sampling characteristics in which rapid and frequent transitions between the protonated and unprotonated states at pH = pKa are achieved. Our calculated pKa values of simple nucleotides are in a good agreement with experimentally measured values, with a mean absolute error of 0.24 pKa units. This work demonstrates that CPHMD can be used as a powerful tool to investigate pH-dependent biological properties of RNA macromolecules.

[1]  Jennifer L. Knight,et al.  Multi-Site λ-dynamics for simulated Structure-Activity Relationship studies. , 2011, Journal of chemical theory and computation.

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

[3]  Miguel Machuqueiro,et al.  Constant-pH molecular dynamics simulations reveal a β-rich form of the human prion protein. , 2010, The journal of physical chemistry. B.

[4]  Charles L. Brooks,et al.  Applying efficient implicit nongeometric constraints in alchemical free energy simulations , 2011, J. Comput. Chem..

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

[6]  Matthias Buck,et al.  Refinement of the primary hydration shell model for molecular dynamics simulations of large proteins , 2009, J. Comput. Chem..

[7]  A. Warshel,et al.  Free energy of charges in solvated proteins: microscopic calculations using a reversible charging process. , 1986, Biochemistry.

[8]  S. Strobel,et al.  Investigation of adenosine base ionization in the hairpin ribozyme by nucleotide analog interference mapping. , 2001, RNA.

[9]  Benoît Roux,et al.  Dominant solvation effects from the primary shell of hydration: Approximation for molecular dynamics simulations , 1995 .

[10]  Bernard R Brooks,et al.  Toward canonical ensemble distribution from self-guided Langevin dynamics simulation. , 2011, The Journal of chemical physics.

[11]  Jennifer A. Doudna,et al.  A conformational switch controls hepatitis delta virus ribozyme catalysis , 2004, Nature.

[12]  Nils G Walter,et al.  Molecular dynamics suggest multifunctionality of an adenine imino group in acid-base catalysis of the hairpin ribozyme. , 2009, RNA.

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

[14]  A. Baptista,et al.  Molecular dynamics at constant pH and reduction potential: application to cytochrome c(3). , 2009, Journal of the American Chemical Society.

[15]  G. Ciccotti,et al.  Numerical Integration of the Cartesian Equations of Motion of a System with Constraints: Molecular Dynamics of n-Alkanes , 1977 .

[16]  D. York,et al.  Quantum mechanical/molecular mechanical simulation study of the mechanism of hairpin ribozyme catalysis. , 2008, Journal of the American Chemical Society.

[17]  K. Réblová,et al.  Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme. , 2010, The journal of physical chemistry. B.

[18]  I. Shih,et al.  A pH-sensitive RNA tertiary interaction affects self-cleavage activity of the HDV ribozymes in the absence of added divalent metal ion. , 2001, Journal of molecular biology.

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

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

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

[22]  Alexander D. MacKerell,et al.  Impact of 2′‐hydroxyl sampling on the conformational properties of RNA: Update of the CHARMM all‐atom additive force field for RNA , 2011, J. Comput. Chem..

[23]  Emil Alexov,et al.  Calculation of pKas in RNA: on the structural origins and functional roles of protonated nucleotides. , 2007, Journal of molecular biology.

[24]  A. Warshel,et al.  Calculations of electrostatic energies in proteins. The energetics of ionized groups in bovine pancreatic trypsin inhibitor. , 1985, Journal of molecular biology.

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

[26]  Arieh Warshel,et al.  Simulation of enzyme reactions using valence bond force fields and other hybrid quantum/classical approaches , 1993 .

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

[28]  L. Scott,et al.  The pH Dependence of Hairpin Ribozyme Catalysis Reflects Ionization of an Active Site Adenine* , 2011, The Journal of Biological Chemistry.

[29]  Jeffrey H. Miller,et al.  Competition between MutY and Mismatch Repair at A · C Mispairs In Vivo , 2003, Journal of bacteriology.

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

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

[32]  Jana K. Shen,et al.  Toward accurate prediction of pKa values for internal protein residues: The importance of conformational relaxation and desolvation energy , 2011, Proteins.

[33]  Maciej Gutowski,et al.  Stabilization of very rare tautomers of uracil by an excess electron. , 2005, Physical chemistry chemical physics : PCCP.

[34]  M. Fedor,et al.  Role of an active site adenine in hairpin ribozyme catalysis. , 2005, Journal of molecular biology.

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

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

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

[38]  Maciej Gutowski,et al.  Stabilization of very rare tautomers of 1-methylcytosine by an excess electron. , 2005, The journal of physical chemistry. A.

[39]  B. Golden,et al.  Direct measurement of a pK(a) near neutrality for the catalytic cytosine in the genomic HDV ribozyme using Raman crystallography. , 2007, Journal of the American Chemical Society.

[40]  Arieh Warshel,et al.  Consistent Calculations of pKa's of Ionizable Residues in Proteins: Semi-microscopic and Microscopic Approaches , 1997 .

[41]  I. Shih,et al.  Involvement of a cytosine side chain in proton transfer in the rate-determining step of ribozyme self-cleavage. , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[42]  Arieh Warshel,et al.  Microscopic and semimicroscopic calculations of electrostatic energies in proteins by the POLARIS and ENZYMIX programs , 1993, J. Comput. Chem..

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

[44]  B. Brooks,et al.  Self-guided Langevin dynamics simulation method , 2003 .

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

[46]  Wei Yang,et al.  Simultaneous escaping of explicit and hidden free energy barriers: application of the orthogonal space random walk strategy in generalized ensemble based conformational sampling. , 2009, The Journal of chemical physics.

[47]  G. Ullmann,et al.  Theoretical investigation of the behavior of titratable groups in proteins , 2006, Photochemical & photobiological sciences : Official journal of the European Photochemistry Association and the European Society for Photobiology.

[48]  P. Bevilacqua,et al.  Driving forces for nucleic acid pK(a) shifting in an A(+).C wobble: effects of helix position, temperature, and ionic strength. , 2010, Biochemistry.

[49]  A. Baptista,et al.  The pH-dependent conformational states of kyotorphin: a constant-pH molecular dynamics study. , 2007, Biophysical journal.

[50]  Wei Yang,et al.  Random walk in orthogonal space to achieve efficient free-energy simulation of complex systems , 2008, Proceedings of the National Academy of Sciences.

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

[52]  D. Patel,et al.  NMR studies of A.C mismatches in DNA dodecanucleotides at acidic pH. Wobble A(anti).C(anti) pair formation. , 1987, The Journal of biological chemistry.

[53]  D. Davis,et al.  Stabilization of the anticodon stem-loop of tRNALys,3 by an A+-C base-pair and by pseudouridine. , 1999, Journal of molecular biology.

[54]  Arieh Warshel,et al.  Using simplified protein representation as a reference potential for all-atom calculations of folding free energy , 1999 .

[55]  J. Piccirilli,et al.  General acid catalysis by the hepatitis delta virus ribozyme , 2005, Nature chemical biology.

[56]  A. Lane,et al.  The contribution of cytosine protonation to the stability of parallel DNA triple helices. , 1998, Journal of molecular biology.

[57]  O. Krechkivska,et al.  Transitions to high tautomeric states can be induced in adenine by interactions with carboxylate and sodium ions: DFT calculation data , 2004 .

[58]  Arieh Warshel,et al.  Microscopic models for quantum mechanical calculations of chemical processes in solutions: LD/AMPAC and SCAAS/AMPAC calculations of solvation energies , 1992 .

[59]  Arieh Warshel,et al.  Multiscale simulations of protein landscapes: Using coarse‐grained models as reference potentials to full explicit models , 2010, Proteins.

[60]  Mark A Olson,et al.  Evaluation of Poisson solvation models using a hybrid explicit/implicit solvent method. , 2005, The journal of physical chemistry. B.

[61]  S. Nakano,et al.  General acid-base catalysis in the mechanism of a hepatitis delta virus ribozyme. , 2000, Science.

[62]  J. Feigon,et al.  Adenine protonation in domain B of the hairpin ribozyme. , 2000, Biochemistry.

[63]  A. Warshel,et al.  Polarization constraints in molecular dynamics simulation of aqueous solutions: The surface constraint all atom solvent (SCAAS) model , 1985 .

[64]  R. Alberty,et al.  The apparent ionization constants of the adenosinephosphates and related compounds. , 1951, The Journal of biological chemistry.

[65]  D. Lilley,et al.  Nucleobase catalysis in the hairpin ribozyme. , 2006, RNA.

[66]  J. Wedekind,et al.  Direct Raman measurement of an elevated base pKa in the active site of a small ribozyme in a precatalytic conformation. , 2009, Journal of the American Chemical Society.

[67]  J. Antosiewicz,et al.  Effects of solute-solvent proton exchange on polypeptide chain dynamics: a constant-pH molecular dynamics study. , 2005, The journal of physical chemistry. B.

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

[69]  L. Cavalieri.  Studies on the Structure of Nucleic Acids. VII. On the Identification of the Isomeric Cytidylic and Adenylic Acids1 , 1953 .

[70]  Charles L. Brooks,et al.  λ‐Dynamics free energy simulation methods , 2009, J. Comput. Chem..

[71]  Arieh Warshel,et al.  A surface constrained all‐atom solvent model for effective simulations of polar solutions , 1989 .

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

[73]  R. M. Izatt,et al.  Sites and thermodynamic quantities associated with proton and metal ion interaction with ribonucleic acid, deoxyribonucleic acid, and their constituent bases, nucleosides, and nucleotides. , 1971, Chemical reviews.

[74]  Charles L. Brooks,et al.  Efficient and Flexible Algorithm for Free Energy Calculations Using the λ-Dynamics Approach , 1998 .

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

[76]  Philip C Bevilacqua,et al.  Mechanistic characterization of the HDV genomic ribozyme: solvent isotope effects and proton inventories in the absence of divalent metal ions support C75 as the general acid. , 2008, Journal of the American Chemical Society.

[77]  Arieh Warshel,et al.  Monte Carlo simulations of proton pumps: on the working principles of the biological valve that controls proton pumping in cytochrome c oxidase. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[78]  Miguel Machuqueiro,et al.  Constant-pH molecular dynamics with ionic strength effects: protonation-conformation coupling in decalysine. , 2006, The journal of physical chemistry. B.

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

[80]  C. Pleij,et al.  Protonation of non-Watson–Crick base pairs and encapsidation of turnip yellow mosaic virus RNA , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[81]  D. Turner,et al.  A CA+ Pair Adjacent to a Sheared GA or AA Pair Stabilizes Size-Symmetric RNA Internal Loops† , 2009, Biochemistry.

[82]  Vijay S Pande,et al.  A smoothly decoupled particle interface: new methods for coupling explicit and implicit solvent. , 2011, The Journal of chemical physics.

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

[84]  E. L. Holbrook,et al.  Structure of an RNA internal loop consisting of tandem C-A+ base pairs. , 1998, Biochemistry.

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

[86]  Eric Westhof,et al.  The non-Watson-Crick base pairs and their associated isostericity matrices. , 2002, Nucleic acids research.

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

[88]  R. Gutell,et al.  Diversity of base-pair conformations and their occurrence in rRNA structure and RNA structural motifs. , 2004, Journal of molecular biology.

[89]  L. Scott,et al.  Direct Measurement of the Ionization State of an Essential Guanine in the Hairpin Ribozyme , 2009, Nature chemical biology.

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

[91]  D. Draper,et al.  Purine-purine mismatches in RNA helices: evidence for protonated G.A pairs and next-nearest neighbor effects. , 1995, Nucleic acids research.

[92]  D. Case,et al.  A novel view of pH titration in biomolecules. , 2001, Biochemistry.

[93]  B. Golden,et al.  A 1.9 A crystal structure of the HDV ribozyme precleavage suggests both Lewis acid and general acid mechanisms contribute to phosphodiester cleavage. , 2010, Biochemistry.

[94]  T. S. Brown,et al.  Catalytic roles for proton transfer and protonation in ribozymes. , 2004, Biopolymers.

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

[96]  S. Strobel,et al.  pH-dependent conformational flexibility within the ribosomal peptidyl transferase center. , 2001, RNA.

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

[98]  J. Chocholousová,et al.  Implicit solvent simulations of DNA and DNA-protein complexes: agreement with explicit solvent vs experiment. , 2006, The journal of physical chemistry. B.

[99]  D. J. Price,et al.  Conformational change of the methionine 20 loop of Escherichia coli dihydrofolate reductase modulates pKa of the bound dihydrofolate , 2007, Protein science : a publication of the Protein Society.

[100]  T. Halgren.  Merck molecular force field. I. Basis, form, scope, parameterization, and performance of MMFF94 , 1996, J. Comput. Chem..

[101]  J. Antosiewicz,et al.  Constant-pH molecular dynamics simulations: a test case of succinic acid , 2004 .

[102]  Mark A Olson,et al.  An efficient hybrid explicit/implicit solvent method for biomolecular simulations , 2004, J. Comput. Chem..