Abstract:Accurate free energy estimation is essential for the prediction of RNA secondary structure. The widely Turner model works well for nested structures. For pseudoknotted RNAs, there is no effective free energy rules. We present a new free energy estimation method that goes beyond the Turner model based on the physical consideration that the free energy of loops (including pseudoknotted loops) are mostly determined by the loop length and the end-end distance. Given a specific secondary structure (refer to as topology hereafter), we grow the RNA in three dimensional space and search among all possible arrangement of helices for stable structure. The free energy contribution of helices is estimated from the Turner model, and the free energy of loops is calculated based on our new approach. We assess our free energy model by testing its ability in selecting native structure from decoy topologies. For a testing set of 50-pseudoknotted RNAs, we create thousands of competitive decoy topologies and calculate their free energy. On average, our energy model provides an average of 7-8 percent improvement in term of sensitivity and specificity. In a blind pseudoknot structure prediction without inserting the native topology into the candidate list, our prediction has an average of 0.79 sensitives and 0.75 specificity. This performance is comparable to that of TdFOLD (0.76, 0.79) and better than Hotknots (0.69, 0.72), ILM(0.73, 0.69) and PKNOTS (0.78,0.73). Our free energy estimation method is fundamentally different from the others: it is based on the growing of conformations in three dimensional space and is well suited for predicting pseudoknotted structures. In addition, our method also gives coarse 3D conformation of a given secondary structure, which can lead to accurate 3D ab initio predicted structure upon further refinement.

RNA Secondary Structure Prediction Using a Genetic Algorithm with a Selection Method Based on Free Energy Value and Topological Index

A Hopfield Neural Network Based Algorithm for RNA Secondary Structure Prediction

RNA Secondary Structure Prediction Based on Forest Representation and Genetic Algorithm

RGRNA: prediction of RNA secondary structure based on replacement and growth of stems

RNA folding algorithm based on fuzzy dynamic programming

PSRna: Prediction of small RNA secondary structures based on reverse complementary folding method

Evaluation of RNA secondary structure prediction for both base-pairing and topology

A Fuzzy Dynamic Programming Approach to Predict RNA Secondary Structure

Novel Representation of RNA Secondary Structure Used to Improve Prediction Algorithms.

A dynamic programming algorithm for RNA structure prediction including pseudoknots

GA Combined with Structural Comparison to Improve the Prediction of RNA Secondary Structure from Comparative Sequence Alignment

RNA Secondary Structure Thermodynamics

McGenus: A Monte Carlo algorithm to predict RNA secondary structures with pseudoknots

An Approach For Rna Secondary Structure Prediction Based On Bayesian Network

A Helix-based Minimum Free Energy Algorithm for RNA Secondary Structure Prediction

A New Free Energy Rule Toward Psuedo-knotted RNA Structure Prediction

A Parallel Intelligence Algorithm for RNA Secondary Structure Prediction

Exploring the Efficiency of Deep Graph Neural Networks for RNA Secondary Structure Prediction

Detecting Ribosnitches with RNA Folding Algorithms: a Genome-Wide Benchmark

RNA fold prediction by Monte Carlo in graph space and the statistical mechanics of tertiary interactions

Improved Rna Secondary Structure Prediction by Maximizing Expected Pair Accuracy