pH-Dependent Dynamics of Complex RNA Macromolecules

The role of pH-dependent protonation equilibrium in modulating RNA dynamics and function is one of the key unanswered questions in RNA biology. Molecular dynamics (MD) simulations can provide insight into the mechanistic roles of protonated nucleotides, but it is only capable of modeling fixed protonation states and requires prior knowledge of the key residue's protonation state. Recently, we developed a framework for constant pH molecular dynamics simulations (CPHMDMS lambda D) of nucleic acids, where the nucleotides' protonation states are modeled as dynamic variables that are coupled to the structural dynamics of the RNA. In the present study, we demonstrate the application of CPHMDMS lambda D to the lead-dependent ribozyme; establishing the validity of this approach for modeling complex RNA structures. We show that CPHMDMS lambda D) accurately predicts the direction of the plc shifts and reproduces experimentally measured microscopic pK(a) values with an average unsigned error of 1.3 pK(a) units. The effects of coupled titration states in RNA structures are modeled, and the importance of conformation sampling is highlighted. The general accuracy of CPHMDMS lambda D simulations in reproducing pH-dependent observables reported in this work demonstrates that constant pH simulations provide a powerful tool to investigate pH-dependent processes in nucleic acids.