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.

van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA–Protein Complexes

Empirical corrections to the Amber RNA force field with Target Metadynamics

Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics

Simulations of the molecular dynamics of nucleic acids

Structures and Dynamics of DNA Mini-Dumbbells Are Force Field Dependent

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Improved Parameterization of Amine-Carboxyate and Amine-Phosphate Interactions for Molecular Dynamics Simulations Using the CHARMM and AMBER Force Fields

Molecular simulations matching denaturation experiments for N6-Methyladenosine

modXNA: A Modular Approach to Parametrization of Modified Nucleic Acids for Use with Amber Force Fields

Molecular Dynamics Simulations of RNA Motifs to Guide the Architectural Parameters and Design Principles of RNA Nanostructures

RNA-Specific Force Field Optimization with CMAP and Reweighting

RNA force field with accuracy comparable to state-of-the-art protein force fields

Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields

The development of nucleic acids force fields: From an unchallenged past to a competitive future

Z-DNA as a Touchstone for Additive Empirical Force Fields and a Refinement of the Alpha/Gamma DNA torsions for AMBER

Combining simulations and solution experiments as a paradigm for RNA force field refinement

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

Can We Ever Develop an Ideal RNA Force Field? Lessons Learned from Simulations of UUCG RNA Tetraloop and Other Systems.

Understanding nucleic acid structural changes by comparing wide-angle x-ray scattering (WAXS) experiments to molecular dynamics simulations

Ab Initio Nucleic Acid Folding Simulations Using a Physics-Based Atomistic Free Energy Model