Simulating RNA folding kinetics on approximated energy landscapes

Area of Science:

• Computational biology
• Biophysics
• Molecular biology

Background:

• RNA function is often dictated by folding kinetics, not solely by sequence or lowest free-energy state.
• Understanding RNA folding pathways is crucial for elucidating structure-function relationships.

Purpose of the Study:

• To present a general computational approach for simulating RNA folding kinetics.
• To extract population kinetics, folding rates, and intermediate substructure formation.
• To enable the study of larger RNA molecules and their kinetic behavior.

Main Methods:

• Building an approximate folding energy landscape map.
• Solving the master equation on the energy map to analyze population kinetics.
• Employing map-based Monte Carlo simulations to extract folding pathways.

Main Results:

• The approximate energy map captures major features of the complete landscape.
• The method scales effectively to simulate kinetics of RNAs over 200 nucleotides.
• Accurate computation of kinetics-based functional rates for ColE1 RNAII and MS2 phage RNAs, showing excellent experimental agreement.

Conclusions:

• The developed computational method provides a scalable and accurate approach to simulate RNA folding kinetics.
• This technique offers valuable insights into RNA dynamics and function, complementing traditional structure-based analyses.
• The findings validate the method's utility for studying complex biological systems like wild-type and mutant RNAs.