RNA stability under different combinations of amber force fields and solvation models

Area of Science:

• Computational Chemistry
• Biophysics
• Molecular Modeling

Background:

• Accurate molecular dynamics (MD) simulations of biomolecules require appropriate force field and solvent model selection.
• While extensively studied for proteins, the optimal combination for RNA simulations remains less explored.
• This gap hinders reliable computational analysis of RNA structure and dynamics.

Purpose of the Study:

• To evaluate the performance of different Amber force field parameters (ff98, ff99, ff99bsc0) combined with General Born implicit solvent models (igb1, igb2, igb5).
• To identify the most suitable force field-solvent combinations for maintaining RNA stability and native conformations in MD simulations.
• To provide guidance for future RNA molecular dynamics studies.

Main Methods:

• Utilized standard molecular dynamics (MD) and replica exchange molecular dynamics (REMD) simulations.
• Tested various combinations of Amber force fields (ff98, ff99, ff99bsc0) and General Born implicit solvent models (igb1, igb2, igb5).
• Assessed the stability of RNA hairpin and duplex structures under different simulation conditions.

Main Results:

• Only the combination of the ff98 force field with the igb1 implicit solvent model, and the ff99bsc0 force field with igb1, successfully preserved the native conformations of RNA hairpin and duplex.
• Other tested force field-solvent combinations failed to maintain the structural integrity of the RNA molecules.
• These findings highlight the sensitivity of RNA simulations to parameter choices.

Conclusions:

• The combination of the ff98 Amber force field and the igb1 General Born implicit solvent model is a recommended choice for molecular dynamics simulations of RNA.
• This specific parameter set demonstrated the ability to accurately reproduce and maintain native RNA conformations.
• Further research may explore additional force field-solvent combinations or explicit solvent models for enhanced RNA simulation accuracy.