Comprehensive Comparison and Critical Assessment of RNA-Specific Force Fields.

Abstract

Molecular dynamics simulations play a pivotal role in elucidating the dynamic behaviors of RNA structures, offering a valuable complement to traditional methods such as nuclear magnetic resonance or X-ray. Despite this, the current precision of RNA force fields lags behind that of protein force fields. In this work, we systematically compared the performance of four RNA force fields (ff99bsc0χOL3, AMBERDES, ff99OL3_CMAP1, AMBERMaxEnt) across diverse RNA structures. Our findings highlight significant challenges in maintaining stability, particularly with regard to cross-strand and cross-loop hydrogen bonds. Furthermore, we observed the limitations in accurately describing the conformations of nonhelical structural motif, terminal nucleotides, and also base pairing and base stacking interactions by the tested RNA force fields. The identified deficiencies in existing RNA force fields provide valuable insights for subsequent force field development. Concurrently, these findings offer recommendations for selecting appropriate force fields in RNA simulations.