Abstract:Molecular dynamics (MD) simulations are an important and well-established tool for investigating RNA structural dynamics, but their accuracy relies heavily on the quality of the employed force field ( ). In this work, we present a comprehensive evaluation of widely used pair-additive and polarizable RNA s using the challenging UUCG tetraloop (TL) benchmark system. Extensive standard MD simulations, initiated from the NMR structure of the 14-mer UUCG TL, revealed that most s did not maintain the native state, instead favoring alternative loop conformations. Notably, three very recent variants of pair-additive s, OL3 —gHBfix21, DESAMBER, and OL3 , successfully preserved the native structure over a 10 × 20 μs timescale. To further assess these s, we performed enhanced sampling folding simulations of the shorter 8-mer UUCG TL, starting from the single-stranded conformation. Estimated folding free energies (ΔG° ) varied significantly among these three s, with values of 0.0 ± 0.6 kcal/mol, 2.4 ± 0.8 kcal/mol, and 7.4 ± 0.2 kcal/mol for OL3 —gHBfix21, DESAMBER, and OL3 , respectively. The ΔG° value from OL3 —gHBfix21 was closest to experimental data, while the higher ΔG° values from DESAMBER and OL3 were unexpected, suggesting an over- or underestimation of key interactions within the folded and mis(un)folded ensembles. These discrepancies led us to further test DESAMBER and OL3 s on additional RNA and DNA systems, where further performance issues were observed. Our results emphasize the complexity of accurately modeling RNA dynamics and suggest that creating an RNA capable of reliably performing across a wide range of RNA motifs remains extremely challenging. In conclusion, our study provides valuable insights into the capabilities of current RNA s and highlights key areas for future development.

Simulations of the molecular dynamics of nucleic acids

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

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

Molecular dynamics analysis of biomolecular systems including nucleic acids

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics

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

Extensive Evaluation of Force Fields for G-Quadruplexes

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

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

Molecular Dynamics Simulations of Protein RNA Complexes by Using an Advanced Electrostatic Model

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

Molecular simulations of proteins: From simplified physical interactions to complex biological phenomena

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

Tumuc1: A New Accurate DNA Force Field Consistent with High-Level Quantum Chemistry

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

Systematic Validation of Protein Force Fields against Experimental Data

Molecular Dynamics Simulation for All

Biomolecular Simulation: A Computational Microscope for Molecular Biology

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

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

Development of models and methods for the molecular simulation of large systems and molecules