Reparameterization of the Amber RNA Force Field Non-Bonded Terms

Area of Science:

• Computational Chemistry
• Biomolecular Modeling
• Structural Biology

Background:

• All-atom molecular dynamics simulations are crucial for understanding RNA behavior, including folding and drug interactions.
• Existing RNA force fields, like Amber OL3, have limitations in accurately stabilizing native RNA structures.
• Accurate force fields are essential for reliable prediction of RNA conformational preferences and dynamics.

Purpose of the Study:

• To develop a novel approach for fitting non-bonded parameters in RNA force fields.
• To improve the stability of native RNA structures in molecular dynamics simulations.
• To create a new set of force field parameters, Amber RNA.ROC26, for enhanced RNA modeling.

Main Methods:

• Fitting non-bonded parameters (point charges and Lennard-Jones) using quantum mechanics (QM) interaction energies.
• Utilizing symmetry-adapted perturbation theory (SAPT) with embedded point charges to model electrostatic environments.
• Testing the new Amber RNA.ROC26 parameters on RNA tetraloops and tetramers.

Main Results:

• The Amber RNA.ROC26 parameters significantly improved the stability of native RNA structures compared to Amber OL3.
• Microsecond timescale simulations demonstrated stable native tetraloop conformations, which previously unfolded.
• Conformational preferences of tetramers were enhanced, although A-form helices showed increased internal loop flexibility.

Conclusions:

• The new method for fitting RNA force field parameters to SAPT interaction energies is promising.
• Amber RNA.ROC26 offers improved stability and conformational accuracy for RNA simulations.
• This approach provides a flexible foundation for future revisions and extensions to model RNA-interacting molecules.