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

Vojtech Mlynsky,Petra Kuehrova,Martin Pykal,Miroslav Krepl,Petr Stadlbauer,Michal Otyepka,Pavel Banas,Jiri Sponer

DOI: https://doi.org/10.1101/2024.10.08.617195

2024-10-11

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.

Biophysics

What problem does this paper attempt to address?

The paper aims to address the following issues:

1. **Evaluating the performance of RNA force fields**: The paper assesses the performance of commonly used RNA force fields in simulating the UUCG tetraloop (TL) system through extensive molecular dynamics (MD) simulations. The study finds that most force fields fail to maintain the native structure and tend to form alternative loop conformations.

2. **Identifying force fields that can stabilize the UUCG tetraloop structure**: Despite the poor performance of most force fields, the paper mentions that three of the latest force fields (OL3 CP–gHBfix21, DESAMBER, and OL3 R2.7) successfully maintain the native structure of the UUCG tetraloop over long timescales (10×20 microseconds).

3. **Evaluating folding free energy**: To further assess these force fields, the researchers conducted enhanced sampling folding simulations and calculated the folding free energy (ΔG°fold) of the 8-mer UUCG tetraloop. The results show significant differences in the ΔG°fold values predicted by different force fields, with the value from OL3 CP–gHBfix21 being the closest to experimental data.

4. **Exploring the limitations of RNA force fields**: Through standard simulations of other RNA and DNA systems, the study reveals performance issues with the DESAMBER and OL3 R2.7 force fields in certain scenarios, highlighting the complexity of accurately simulating RNA dynamics and pointing out that developing a force field that reliably works across various RNA motifs remains highly challenging.

In summary, the main goal of the paper is to provide valuable insights for the future development of force fields by thoroughly evaluating the performance of existing RNA force fields in simulating specific RNA structures, such as the UUCG tetraloop.