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.

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

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

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

A dynamic programming algorithm for RNA structure prediction including pseudoknots

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

memerna: Sparse RNA Folding Including Coaxial Stacking

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

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

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

Free Energy Landscape and Multiple Folding Pathways of an H-Type RNA Pseudoknot

A Helix-based Minimum Free Energy Algorithm for RNA Secondary Structure Prediction

Predicting Model And Algorithm In Rna Folding Structure Including Pseudoknots

Improved prediction of RNA secondary structure by integrating the free energy model with restraints derived from experimental probing data

RNA Folding on the 3D Triangular Lattice

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

Predicting Algorithm of RNA Folding Structure with Pseudoknots.

A divide-and-conquer approach based on deep learning for long RNA secondary structure prediction: focus on pseudoknots identification

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

Minimum-Free-Energy Distribution of RNA Secondary Structures: Entropic and Thermodynamic Properties of Rare Events

Improved Predicting Algorithm of RNA Pseudoknotted Structure

RNA folding algorithm based on fuzzy dynamic programming