Abstract:The information-encoding molecules RNA and DNA form a combinatorially large set of secondary structures through nucleic acid base pairing. Thermodynamic prediction algorithms predict favoured, or minimum free energy (MFE), secondary structures, and can assign an equilibrium probability to any structure via the partition function: a Boltzman-weighted sum over the set of secondary structures. MFE is NP-hard in the presence pseudoknots, base pairings that violate a restricted planarity condition. However, unpseudoknotted structures are amenable to dynamic programming: for a single DNA/RNA strand there are polynomial time algorithms for MFE and partition function. For multiple strands, the problem is more complicated due to entropic penalties. Dirks et al [SICOMP Review; 2007] showed that for O(1) strands, with N bases, there is a polynomial time in N partition function algorithm, however their technique did not generalise to MFE which they left open. We give the first polynomial time (O(N^4)) algorithm for unpseudoknotted multiple (O(1)) strand MFE, answering the open problem from Dirks et al. The challenge lies in considering rotational symmetry of secondary structures, a feature not immediately amenable to dynamic programming algorithms. Our proof has two main technical contributions: First, a polynomial upper bound on the number of symmetric secondary structures to be considered when computing rotational symmetry penalties. Second, that bound is leveraged by a backtracking algorithm to find the MFE in an exponential space of contenders. Our MFE algorithm has the same asymptotic run time as Dirks et al's partition function algorithm, suggesting efficient handling of rotational symmetry, although higher space complexity. It also seems reasonably tight in the number of strands since Codon, Hajiaghayi & Thachuk [DNA27, 2021] have shown that unpseudoknotted MFE is NP-hard for O(N) strands.

Predicting scheme of RNA folding structure including pseudoknots.

Predicting Algorithm of RNA Folding Structure with Pseudoknots.

Improved Predicting Algorithm of RNA Pseudoknotted Structure

Predicting Model And Algorithm In Rna Folding Structure Including Pseudoknots

A dynamic programming algorithm for RNA structure prediction including pseudoknots

Prediction of Geometrically Feasible Three-Dimensional Structures of Pseudoknotted RNA Through Free Energy Estimation

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

A Hopfield Neural Network Based Algorithm for RNA Secondary Structure Prediction

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

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

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

An Algorithm for RNA Pseudoknots

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

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

A Local Structural Prediction Algorithm for RNA Triple Helix Structure

Predicting 3D structure and stability of RNA pseudoknots in monovalent and divalent ion solutions

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

An efficient algorithm to compute the minimum free energy of interacting nucleic acid strands

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

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

RNA-Puzzles Round II: assessment of RNA structure prediction programs applied to three large RNA structures