Abstract:The Kinfold and KFOLD programs for RNA folding kinetics implement the Gillespie algorithm to generate stochastic folding trajectories from an initial structure s to a target structure t, in which each intermediate secondary structure is obtained from its predecessor by the addition, removal or shift of a single base pair. Define MS2 distance between secondary structures s and t to be the minimum path length to refold s to t, where a move from MS2 is applied in each step. We describe algorithms to compute the shortest MS2 folding trajectory between any two given RNA secondary structures. These algorithms include an optimal integer programming (IP) algorithm, an accurate and efficient near-optimal algorithm, a greedy algorithm, a branch-and-bound algorithm, and an optimal algorithm if one allows intermediate structures to contain pseudoknots. Our optimal IP [resp. near-optimal IP] algorithm maximizes [resp. approximately maximizes] the number of shifts and minimizes [resp. approximately minimizes] the number of base pair additions and removals by applying integer programming to (essentially) solve the minimum feedback vertex set (FVS) problem for the RNA conflict digraph, then applies topological sort to tether subtrajectories into the final optimal folding trajectory. We prove NP-hardness of the problem to determine the minimum barrier energy over all possible MS2 folding pathways, and conjecture that computing the MS2 distance between arbitrary secondary structures is NP-hard. Since our optimal IP algorithm relies on the FVS, known to be NP-complete for arbitrary digraphs, we compare the family of RNA conflict digraphs with the following classes of digraphs (planar, reducible flow graph, Eulerian, and tournament) for which FVS is known to be either polynomial time computable or NP-hard. Source code available at <a class="link-external link-http" href="http://bioinformatics.bc.edu/clotelab/MS2distance/" rel="external noopener nofollow">this http URL</a>.

Approximating the Maximum Multiple RNA Interaction Problem

Approximation Algorithms for the Maximum Multiple RNA Interaction Problem.

Approximation Algorithms for Predicting RNA Secondary Structures with Arbitrary Pseudoknots

Improved Approximation Algorithms for Predicting RNA Secondary Structures with Arbitrary Pseudoknots

Improved Predicting Algorithm of RNA Pseudoknotted Structure

BPPart and BPMax: RNA-RNA Interaction Partition Function and Structure Prediction for the Base Pair Counting Model

Predicting Model And Algorithm In Rna Folding Structure Including Pseudoknots

Predicting Algorithm of RNA Folding Structure with Pseudoknots.

Prediction of RNA-RNA interaction structure by centroids in the Boltzmann ensemble

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

Modeling Noncanonical RNA Base Pairs by a Coarse-Grained IsRNA2 Model

Improved Rna Secondary Structure Prediction by Maximizing Expected Pair Accuracy

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

Combinatorial RNA Design Designability and Structure-Approximating Algorithm in Watson-Crick and Nussinov-Jacobson Energy Models

Protein-RNA Interaction Interface Prediction and Design

Non-linear interaction in random matrix models of RNA

A Method for Predicting Rna-Protein Interaction and Interaction Sites

Minimum length RNA folding trajectories

Taking RNA-RNA Interaction to Machine Peak

Combinatorial RNA Design: Designability and Structure-Approximating Algorithm

Homology modeling of Complex Structural RNAs