Base-specific RNA force field improving the dynamics conformation of nucleotide

Area of Science:

• Computational biology
• Biophysics
• Molecular modeling

Background:

• RNA's crucial roles in biological processes necessitate accurate structural and dynamic analysis.
• Traditional experimental methods face limitations in capturing RNA's dynamic conformations.
• Molecular dynamics (MD) simulations are vital for complementing experimental RNA studies.

Purpose of the Study:

• To address limitations in existing RNA force fields, specifically overestimation of base stacking and high ratios of intercalated conformations.
• To develop an optimized RNA force field, BSFF1, for improved accuracy in molecular dynamics simulations.
• To enhance the ability of computational methods to predict RNA structure and dynamics.

Main Methods:

• A two-step optimization strategy was applied to the ff99bsc0χOL3 force field, resulting in BSFF1.
• Adjustments included unbonded parameters of nucleobase heavy atoms and the addition of a ζ/α grid-based energy correction map with reweighting.
• MD simulations were performed on tetranucleotides, single-strand RNA, kink-turn structures, duplexes, riboswitches, and tetraloops.

Main Results:

• BSFF1 significantly reduced the ratio of intercalated conformations in tetranucleotide simulations.
• The BSFF1 force field demonstrated improved accuracy in reproducing conformers for single-strand RNA and kink-turn structures compared to ff99bsc0χOL3.
• BSFF1 successfully stabilized conformers of duplex and riboswitch RNA structures.
• Ab initio folding of tetraloops was achieved, validating the performance of BSFF1.

Conclusions:

• The newly developed BSFF1 force field effectively addresses limitations of previous RNA force fields.
• BSFF1 enhances the accuracy of molecular dynamics simulations for RNA conformational sampling.
• This improved force field will aid in a deeper understanding of RNA structure-function relationships.