Abstract:Recently, several experimental techniques have emerged for probing RNA structures based on high-throughput sequencing. However, most secondary structure prediction tools that incorporate probing data are designed and optimized for particular types of experiments. For example, RNAstructure-Fold is optimized for SHAPE data, while SeqFold is optimized for PARS data. Here, we report a new RNA secondary structure prediction method, restrained MaxExpect (RME), which can incorporate multiple types of experimental probing data and is based on a free energy model and an MEA (maximizing expected accuracy) algorithm. We first demonstrated that RME substantially improved secondary structure prediction with perfect restraints (base pair information of known structures). Next, we collected structure-probing data from diverse experiments (e.g. SHAPE, PARS and DMS-seq) and transformed them into a unified set of pairing probabilities with a posterior probabilistic model. By using the probability scores as restraints in RME, we compared its secondary structure prediction performance with two other well-known tools, RNAstructure-Fold (based on a free energy minimization algorithm) and SeqFold (based on a sampling algorithm). For SHAPE data, RME and RNAstructure-Fold performed better than SeqFold, because they markedly altered the energy model with the experimental restraints. For high-throughput data (e.g. PARS and DMS-seq) with lower probing efficiency, the secondary structure prediction performances of the tested tools were comparable, with performance improvements for only a portion of the tested RNAs. However, when the effects of tertiary structure and protein interactions were removed, RME showed the highest prediction accuracy in the DMS-accessible regions by incorporating in vivo DMS-seq data.

Quantitative DMS mapping for automated RNA secondary structure inference

A quantitative framework for structural interpretation of DMS reactivity

Toward a Next-Generation Atlas of RNA Secondary Structure

Simultaneous mapping of RNA secondary, tertiary, and quaternary structure in living cells by multi-site DMS probing

Modeling RNA secondary structure folding ensembles using SHAPE mapping data

High-throughput determination of RNA tertiary contact thermodynamics by quantitative DMS chemical mapping

Web-based platform for analysis of RNA folding from high throughput chemical probing data

Isoform-specific RNA structure determination using Nano-DMS-MaP

STRUCT: a statistical approach to identify RNA secondary structures from raw sequencing data, bypassing multiple sequence alignment

Differential Analysis of RNA Structure Probing Experiments at Nucleotide Resolution: Uncovering Regulatory Functions of RNA Structure

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

A new reagent for in vivo structure probing of RNA G and U residues that improves RNA structure prediction alone and combined with DMS

Correcting a SHAPE-directed RNA structure by a mutate-map-rescue approach

RNA-MobiSeq: Deep mutational scanning and mobility-based selection for RNA structure inference

Two-dimensional chemical mapping for non-coding RNAs

SHAPE-Seq 2.0: systematic optimization and extension of high-throughput chemical probing of RNA secondary structure with next generation sequencing

Modern Semiempirical Electronic Structure Methods and Machine Learning Potentials for Drug Discovery: Conformers, Tautomers, and Protonation States

RNA Structure Probing Reveals the Structural Basis of Dicer Binding and Cleavage.

Characteristic chemical probing patterns of loop motifs improve prediction accuracy of RNA secondary structures.

RNA tertiary structure and conformational dynamics revealed by BASH MaP

Conceptualizing Ribosomal mRNA Unwinding Utilizing DMS MaP-Seq