Molecular Dynamics Simulations of RNA Stem-Loop Folding Using an Atomistic Force Field and a Generalized Born...

Area of Science:

• Computational Biology
• Biophysics
• Molecular Modeling

Background:

• Accurate modeling of ribonucleic acid (RNA) structural dynamics is crucial but challenging.
• RNA stem-loop motifs are fundamental structures with complex folding pathways.

Purpose of the Study:

• To perform de novo folding simulations of diverse RNA stem-loop motifs.
• To assess the accuracy of molecular dynamics simulations with refined force fields and implicit solvent models for RNA folding.

Main Methods:

• Conventional molecular dynamics simulations were used for 26 RNA stem-loops (10–36 residues).
• Simulations employed an atomistic force field refined by the Shaw group and an implicit solvent model by the Simmerling group.
• Structures were analyzed for stem and loop region accuracy using root mean square deviation (RMSD).

Main Results:

• 18 of 18 simple stem-loops folded correctly, retaining native base pairs in stems with RMSD <2 Å for stems.
• 5 of 8 stem-loops with bulges/internal loops were successfully folded, showing stem RMSD of 0.9–4.5 Å.
• Loop region modeling remained challenging, with RMSD around 4 Å for all models.

Conclusions:

• Successful recapitulation of RNA stem folding in fundamental motifs is a significant advancement.
• This study enhances the reliability and accuracy of modeling RNA structural dynamics.
• Further improvements are needed for accurate modeling of RNA loop structures in implicit solvent simulations.