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/.

A Flexible Stem-Based Local Search Algorithm for Predicting RNA Secondary Structures Including Pseudoknots

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

A Local-stem-search Algorithm to Predict the RNA Secondary Structure

StemP: A fast and deterministic Stem-graph approach for RNA and protein folding prediction

A Hopfield Neural Network Based Algorithm for RNA Secondary Structure Prediction

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

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

A dynamic programming algorithm for RNA structure prediction including pseudoknots

Prediction of RNA secondary structure based on stem region replacement using the RSRNA algorithm

Predicting Model And Algorithm In Rna Folding Structure Including Pseudoknots

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

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

Improved Predicting Algorithm of RNA Pseudoknotted Structure

Predicting Algorithm of RNA Folding Structure with Pseudoknots.

A Local Structural Prediction Algorithm for RNA Triple Helix Structure

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

SparseRNAfolD: optimized sparse RNA pseudoknot-free folding with dangle consideration

Computing the partition function for kinetically trapped RNA secondary structures

An efficient dynamic programming algorithm and implementation for RNA secondary structure prediction

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

A Fuzzy Dynamic Programming Approach to Predict RNA Secondary Structure