Abstract:RNA molecules exhibit various biological functions intrinsically dependent on their diverse ecosystem of highly flexible structures. This flexibility arises from complex hydrogen-bonding networks defined by canonical and non-canonical base pairs that require protonation events to stabilize or perturb these interactions. Constant pH molecular dynamics (CpHMD) methods provide a reliable framework to explore the conformational and protonation space of dynamic structures and for robust calculations of pH-dependent properties, such as the pK$_\mathrm{a}$ of titrable sites. Despite growing biological evidence concerning pH regulation of certain motifs and in biotechnological applications, pH-sensitive in silico methods have rarely been applied to nucleic acids. In this work, we extended the stochastic titration CpHMD method to include RNA parameters from the standard $\chi$OL3 AMBER force field and highlighted its capability to depict titration events of nucleotides in single-stranded RNAs. We validated the method using trimers and pentamers with a single central titrable site while integrating a well-tempered metadynamics approach into the st-CpHMD methodology (CpH-MetaD) using PLUMED. This approach enhanced the convergence of the conformational landscape and enabled more efficient sampling of protonation-conformation coupling. Our pK$_\mathrm{a}$ estimates agree with experimental data, validating the method's ability to reproduce electrostatic changes around a titrable nucleobase in single-stranded RNA. These findings provided molecular insight into intramolecular phenomena, such as nucleobase stacking and phosphate interactions, that dictate the experimentally observed pK$_\mathrm{a}$ shifts between different strands. Overall, this work validates both the st-CpHMD and the metadynamics integration as reliable tools for studying biologically relevant RNA systems.

Influence of ion and hydration atmospheres on RNA structure and dynamics: insights from advanced theoretical and computational methods

Unravelling Mg$^{2+}$-RNA binding with atomistic molecular dynamics

Theory and simulations for RNA folding in mixtures of monovalent and divalent cations

Impact of RNA Melting on Hydrating Water Structure Mapped by Femtosecond 2D-IR Spectroscopy

Change of Hydration Patterns upon RNA Melting Probed by Excitations of Phosphate Backbone Vibrations

Counterions in RNA structure: Structural bioinformatics analysis to identify the role of Mg2+ ions in base pair formation

Molecular origin of distinct hydration dynamics in double helical DNA and RNA sequences

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

Elucidating the influence of RNA modifications and Magnesium ions on tRNA conformational dynamics in : Insights from Replica Exchange Molecular Dynamics simulations

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

Ion-mediated RNA structural collapse: effect of spatial confinement

Exploring RNA Destabilization Mechanisms in Biomolecular Condensates through Atomistic Simulations

Modeling Loop Composition and Ion Concentration Effects in RNA Hairpin Folding Stability

Non-specific binding of Na$^+$ and Mg$^{2+}$ to RNA determined by force spectroscopy methods

Principles of ion binding to RNA inferred from the analysis of a 1.55 Å resolution bacterial ribosome structure – Part I: Mg

An Intermediate Resolution Model of RNA Dynamics and Phase Separation with Explicit Mg

Elucidating the influence of RNA modifications and Magnesium ions on tRNA Phe conformational dynamics in S. cerevisiae : Insights from Replica Exchange Molecular Dynamics simulations

Magnesium Contact Ions Stabilize the Tertiary Structure of Transfer RNA: Electrostatics Mapped by Two-Dimensional Infrared Spectra and Theoretical Simulations

Molecular simulations to investigate the impact of N6-methylation in RNA recognition: Improving accuracy and precision of binding free energy prediction

Exploring RNA structure and dynamics through enhanced sampling simulations

Characterizing RNA oligomers using Stochastic Titration Constant-pH Metadynamics simulations