RNA Conformational Sampling: II. Arbitrary Length Multinucleotide Loop Closure.

Area of Science:

• Computational Biology
• Structural Biology
• Biophysics

Background:

• Accurate modeling of RNA structures is crucial for understanding their function.
• Simulating large-scale conformational changes in RNA remains computationally challenging.
• Existing methods struggle with efficiently closing RNA segments of arbitrary length.

Purpose of the Study:

• To generalize the inverse kinematic solution for RNA loop closure.
• To develop a novel Monte Carlo algorithm for efficient sampling of large-scale RNA conformations.
• To demonstrate the algorithm's utility in simulating complex RNA rearrangements.

Main Methods:

• Generalization of inverse kinematics using virtual coordinates (RETO) to reduce problem complexity.
• Formulation and implementation of a new Monte Carlo algorithm for atomistic RNA simulations.
• Exclusive use of torsion angle updates for moving arbitrary length loops.
• Integration with conventional Monte Carlo moves for enhanced conformational sampling.

Main Results:

• A generalized inverse kinematic solution reduces the closure problem from 6x6 to 4x4 constraints.
• A novel Monte Carlo algorithm efficiently samples large-scale RNA conformations.
• The algorithm successfully models RNA segments of arbitrary length without disturbing existing structures.
• Demonstrated utility in simulating the unfolding of a full-length hammerhead ribozyme.

Conclusions:

• The generalized closure solution and new Monte Carlo algorithm significantly improve the efficiency of sampling large-scale RNA conformational changes.
• This approach provides a powerful tool for studying RNA dynamics and complex rearrangements.
• The method is applicable to RNA segments of any length and nucleotide composition.