Abstract:An RNA secondary structure is locally optimal if there is no lower energy structure that can be obtained by the addition or removal of a single base pair, where energy is defined according to the widely accepted Turner nearest neighbor model. Locally optimal structures form kinetic traps, since any evolution away from a locally optimal structure must involve energetically unfavorable folding steps. Here, we present a novel, efficient algorithm to compute the partition function over all locally optimal secondary structures of a given RNA sequence. Our software, RNAlocopt runs in O(n3) time and O(n2) space. Additionally, RNAlocopt samples a user-specified number of structures from the Boltzmann subensemble of all locally optimal structures. We apply RNAlocopt to show that (1) the number of locally optimal structures is far fewer than the total number of structures--indeed, the number of locally optimal structures approximately equal to the square root of the number of all structures, (2) the structural diversity of this subensemble may be either similar to or quite different from the structural diversity of the entire Boltzmann ensemble, a situation that depends on the type of input RNA, (3) the (modified) maximum expected accuracy structure, computed by taking into account base pairing frequencies of locally optimal structures, is a more accurate prediction of the native structure than other current thermodynamics-based methods. The software RNAlocopt constitutes a technical breakthrough in our study of the folding landscape for RNA secondary structures. For the first time, locally optimal structures (kinetic traps in the Turner energy model) can be rapidly generated for long RNA sequences, previously impossible with methods that involved exhaustive enumeration. Use of locally optimal structure leads to state-of-the-art secondary structure prediction, as benchmarked against methods involving the computation of minimum free energy and of maximum expected accuracy. Web server and source code available at http://bioinformatics.bc.edu/clotelab/RNAlocopt/.

K-partite RNA Secondary Structures

A Hopfield Neural Network Based Algorithm for RNA Secondary Structure Prediction

Knotify_V2.0: Deciphering RNA Secondary Structures with H-Type Pseudoknots and Hairpin Loops

Approximation Algorithms for Predicting RNA Secondary Structures with Arbitrary Pseudoknots

Graph-Theoretic Partitioning of RNAs and Classification of Pseudoknots-II

Improved Approximation Algorithms for Predicting RNA Secondary Structures with Arbitrary Pseudoknots

A dynamic programming algorithm for RNA structure prediction including pseudoknots

Deciphering RNA Secondary Structure Prediction: A Probabilistic K-Rook Matching Perspective

[Predicting RNA Secondary Structures Including Pseudoknots by Covariance with Stacking and Minimum Free Energy].

Computing the partition function for kinetically trapped RNA secondary structures

Dual Graph Partitioning Highlights a Small Group of Pseudoknot-Containing RNA Submotifs

Accurate prediction of RNA secondary structure including pseudoknots through solving minimum-cost flow with learned potentials

Concurrent prediction of RNA secondary structures with pseudoknots and local 3D motifs in an integer programming framework

Central and Local Limit Theorems for RNA Structures

CParty: Hierarchically Constrained Partition Function of RNA Pseudoknots

FlexStem: Improving Predictions of RNA Secondary Structures with Pseudoknots by Reducing the Search Space

Cotranscriptional kinetic folding of RNA secondary structures including pseudoknots

Predicting Algorithm of RNA Folding Structure with Pseudoknots.

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

RNA Folding on the 3D Triangular Lattice

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