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

Exploiting parallelization for RNA secondary structure prediction in cluster

Load Balancing Algorithm in Cluster-based RNA secondary structure Prediction

A Hopfield Neural Network Based Algorithm for RNA Secondary Structure Prediction

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

FPGA-Oriented Fine-Grained Algorithm for RNA Secondary Structure Prediction

Fuzzy Kernel Clustering of RNA Secondary Structure Ensemble Using a Novel Similarity Metric.

A Fuzzy Dynamic Programming Approach to Predict RNA Secondary Structure

Fine-grained parallelism accelerating for RNA secondary structure prediction with pseudoknots based on FPGA

A dynamic programming algorithm for RNA structure prediction including pseudoknots

Parallelization and Optimization of Mfold on Shared Memory System

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

Prediction of the RNA Tertiary Structure Based on a Random Sampling Strategy and Parallel Mechanism

RNA folding algorithm based on fuzzy dynamic programming

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

An algorithm for rapid noncoding RNA sequence-structure alignment

Simultaneous prediction of RNA secondary structure and helix coaxial stacking

Time and Space Efficient RNA-RNA Interaction Prediction via Sparse Folding

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

LinearFold: linear-time approximate RNA folding by 5'-to-3' dynamic programming and beam search

Computing the partition function for kinetically trapped RNA secondary structures

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