Abstract:The role of RNA in biology continues to grow but insight into important aspects of RNA behavior is lacking, such as dynamic structural ensembles in different environments, how flexibility is coupled to function, and how function might be modulated by small molecule binding. In the case of proteins, much progress in these areas has been made by complementing experiments with atomistic simulations, but RNA simulation methods and force fields are less mature. It remains challenging to generate stable RNA simulations, even for small systems where well-defined, thermostable structures have been established by experiments. Many different aspects of RNA energetics have been adjusted in force fields, seeking improvements that are transferable across a variety of RNA structural motifs. In this work, we explore the role of weak CH…O interactions, which are ubiquitous in RNA structure but have received less attention in RNA force field development. By comparing data extracted from high-resolution RNA crystal structures to energy profiles from quantum mechanics and force field calculations, we demonstrate that CH…O interactions are overly repulsive in the widely-used Amber RNA force fields. We developed a simple, targeted adjustment of CH…O repulsion that leaves the remainder of the force field unchanged. We then tested the standard and modified force fields using MD simulations with explicit water and salt, amassing over 300 μsec of data for multiple RNA systems containing important features such as presence of loops, base stacking interactions as well as canonical and non-canonical base pairing. Our results demonstrate that the standard force fields lead to reproducible unfolding of the NMR-based structures, as has been reported by others. Including our CH…O adjustment in an otherwise identical protocol dramatically improves the outcome, leading to stable simulations for all RNA systems tested.

Dynamics of Base Pairs with Low Stability in RNA by Solid-State Nuclear Magnetic Resonance Exchange Spectroscopy

Identification of weak non-canonical base pairs around riboswitch-ligand recognition sites by solid-state NMR exchange spectroscopy

Novel Rna Base Pair With Higher Specificity Using Single Selenium Atom

Rational design of hairpin RNA excited states reveals multi-step transitions

Visualizing transient Watson–Crick-like mispairs in DNA and RNA duplexes

Exploring the Molecular Basis of RNA Recognition by the Dimeric RNA-binding Protein Via Molecular Simulation Methods

Adjusting the energy profile for CH–O interactions leads to improved stability of RNA stem-loop structures in MD simulations

Thermodynamic and kinetic analysis of an RNA kissing interaction and its resolution into an extended duplex.

Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch

Interrogating RNA-Small Molecule Interactions with Structure Probing and Artificial Intelligence-Augmented Molecular Simulations

A mechanism for S-adenosyl methionine assisted formation of a riboswitch conformation: a small molecule with a strong arm

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

Local destabilisation of a DNA double helix by a T–T wobble pair

Long-Range Ordered Water Correlations between A–T/C–G Nucleotides

Rupture of DNA Aptamer: new insights from simulations

Phosphorothioate Substitutions in RNA Structure Studied by Molecular Dynamics Simulations, QM/MM Calculations, and NMR Experiments

Bistable secondary structures of small RNAs and their structural probing by comparative imino proton NMR spectroscopy

Extended Structures in RNA Folding Intermediates Are Due to Nonnative Interactions Rather Than Electrostatic Repulsion.

Characterizing RNA oligomers using Stochastic Titration Constant-pH Metadynamics simulations

Structural Insights of Non-Canonical U•U Pair and Hoogsteen Interaction Probed with Se Atom

The Thermodynamic and Kinetic Properties of the dA-rU DNA-RNA Hybrid Base Pair Investigated via Molecular Dynamics Simulations