Constant pH Molecular Dynamics Simulations of Nucleic Acids in Explicit Solvent

The nucleosides of adenine and cytosine have pK(a) 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 pK(a) 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 lambda-dynamics. We adopt the new functional form lambda(Nexp) for lambda that was recently developed for multisite lambda-dynamics (MS lambda D) and demonstrate good sampling characteristics in which rapid and frequent transitions between the protonated and unprotonated states at pH = pK(a) are achieved. Our calculated pK(a) values of simple nucleotides are in a good agreement with experimentally measured values, with a mean absolute error of 0.24 pK(a) units. This work demonstrates that CPHMD can be used as a powerful tool to investigate pH-dependent biological properties of RNA macromolecules.