Hierarchical analysis of RNA secondary structures with pseudoknots based on sections

Area of Science:

• Computational Biology
• Molecular Biology
• Bioinformatics

Background:

• Predicting RNA structures with pseudoknots is computationally intensive.
• Existing methods have high time complexity, limiting their application.

Purpose of the Study:

• To develop a computationally efficient method for RNA pseudoknot prediction.
• To maintain high accuracy in predicting complex RNA structures.

Main Methods:

• A hierarchical approach dividing RNA into sections of unpaired bases.
• Utilizing a nearest-neighbor energy model with dynamic programming for section interactions.
• Analyzing minimum free energy (MFE) gain for section pairs.

Main Results:

• The algorithm achieves a time complexity of O(N^3), a significant improvement over O(N^6).
• Biologically relevant pseudoknots are concentrated in the top 3% of section pairs ranked by MFE gain.
• High prediction accuracy was achieved for 2-clusters (sensitivity > 0.9, PPV > 0.8).
• Asymmetric prediction accuracy for 3-clusters suggests sequential co-transcriptional folding.

Conclusions:

• The hierarchical method offers substantial computational advantages for RNA pseudoknot prediction.
• Findings suggest complex pseudoknots form via sequential folding, not just global energy minimization.
• This provides new insights into RNA folding dynamics and structure prediction.