SHAPE directed RNA folding with the ViennaRNA Package

The ViennaRNA Package 2.0 brings powerful dynamic programming algorithms to researchers studying nucleic acid folding. In this post, we explore three SHAPE-guided methods—Deigan, Zarringhalam, and Washietl—that have been integrated into our toolkit to improve predictions of base pair interactions and minimum free energy (MFE) structures for RNA molecules. By combining chemical probing data with in silico modeling, these approaches help capture real-world folding behaviors and enhance the accuracy of computational RNA structure predictions

Bridging Experiment with Computation

Modern experiments generate SHAPE reactivity profiles through selective 2’-hydroxyl acylation analyzed by primer extension, a method that reveals nucleotide flexibility in vitro or in vivo. To translate these reactivities into structural guidance, the Deigan, Zarringhalam, and Washietl methods each adjust the energy parameters used in traditional dynamic programming: Deigan’s model applies pseudo-energies directly to loops, Zarringhalam’s method uses reactivity-specific weights, and Washietl’s approach infers perturbation of the energy parameters to optimally fit the SHAPE signal. Together, these innovations allow the prediction of both the minimum free energy (MFE) structure and the ensemble of probable folds with improved sensitivity to experimental data.

Performance and Benchmarking

To validate these implementations, we benchmarked against a curated set of RNAs with known reference structures, analyzing metrics like sensitivity and positive predictive value for predicted base pairs. The comparisons covered both the MFE structure and the full partition function ensemble, demonstrating that SHAPE-guided predictions outperform purely thermodynamic models, especially for challenging motifs. Detailed supplementary data outline our benchmark strategies, energy parameter settings, and the impact on diverse RNA molecules, confirming the value of integrating experimental evidence into computational pipelines.

Getting Started: Command Line and Source Code

Users can install binary packages for several Linux distributions or build from source code available on our repository. The command line interface supports commands for folding, partition function calculation, and SHAPE-guided modes, with options to specify and adjust energy parameters directly. Whether you’re scripting large-scale screens or running individual analyses, our tools and comprehensive documentation provide everything needed to apply dynamic programming algorithms to your nucleic acid research.

For a deep dive, don’t miss the Supplementary Data, which includes extensive background on each method, parameter details, and full benchmark results for our SHAPE-integrated algorithms. Dive into the code, explore the options, and see how experimental data can guide your next discovery in RNA structure prediction.

Your browser does not support PDFs. Download the PDF



Citation

SHAPE directed RNA folding
Ronny Lorenz, Dominik Luntzer, Ivo L. Hofacker, Peter F. Stadler, Michael T. Wolfinger
Bioinformatics 32: 145–47 (2016) | doi:10.1093/bioinformatics/btv523 | PDF

See also

Predicting RNA Structures from Sequence and Probing Data
Ronny Lorenz, Michael T. Wolfinger, Andrea Tanzer, Ivo L. Hofacker
Methods 103:86–98 (2016) | doi:10.1016/j.ymeth.2016.04.004 | PDF