key: cord-0867909-5igorndt
authors: Janssen, Stefan; Giegerich, Robert
title: The RNA shapes studio
date: 2015-02-01
journal: Bioinformatics
DOI: 10.1093/bioinformatics/btu649
sha: ed5375ac7d887e311e407121b4c80e8d0ca0afed
doc_id: 867909
cord_uid: 5igorndt

Motivation: Abstract shape analysis, first proposed in 2004, allows one to extract several relevant structures from the folding space of an RNA sequence, preferable to focusing in a single structure of minimal free energy. We report recent extensions to this approach. Results: We have rebuilt the original RNAshapes as a repository of components that allows us to integrate several established tools for RNA structure analysis: RNAshapes, RNAalishapes and pknotsRG, including its recent extension pKiss. As a spin-off, we obtain heretofore unavailable functionality: e. g. with pKiss, we can now perform abstract shape analysis for structures holding pseudoknots up to the complexity of kissing hairpin motifs. The new tool pAliKiss can predict kissing hairpin motifs from aligned sequences. Along with the integration, the functionality of the tools was also extended in manifold ways. Availability and implementation: As before, the tool is available on the Bielefeld Bioinformatics server at http://bibiserv.cebitec.uni-bielefeld.de/rnashapesstudio. Contact: bibi-help@cebitec.uni-bielefeld.de

The framework of algebraic dynamic programming (ADP) allows us to express dynamic programming algorithms for sequence analysis on a high level of abstraction. They are composed from signatures, tree grammars and evaluation algebras (Giegerich et al., 2004a) . Powerful product operations on algebras allow one to derive new types of analysis by the combination of available components, essentially with a single keystroke (Steffen and Giegerich, 2005) . Relying on the recent BELLMAN'S GAP system (Sauthoff et al., 2013) , which implements the ADP framework, we have built a repository of components that allows us to integrate several established tools for RNA structure analysis: RNASHAPES, RNAALISHAPES and PKNOTSRG, including its recent extension PKISS. As a spin-off, we obtain heretofore unavailable functionality: e. g. with PKISS, we can now perform abstract shape analysis for structures holding pseudoknots up to the complexity of kissing hairpin motifs. The new tool PALIKISS can predict kissing hairpin motifs from aligned sequences. Along with the integration, the functionality of the tools was also extended in manifold ways. Figure 1 gives an overview.

1.2.1 Extensions to RNASHAPES It is generally agreed that predicting a single structure of minimal free energy does not adequately capture the subtlety and versatility of RNA structure formation. The RNASHAPES tool introduced the notion of abstract shapes (Giegerich et al., 2004b; 

The work of combines the ideas of RNAALIFOLD and RNASHAPES and performs shape analysis based on pre-aligned RNA sequences. We added the computation of a structure conservation index, different dangling base models, MEA folding and a window mode. RIBOSUM scoring (Bernhart et al., 2008) was added for evaluating sequence similarity. Theis et al. (2010) the ideas of PKNOTSRG (Reeder and Giegerich, 2004) are extended to predict (aside from unknotted structures and H-type pseudoknots) RNA structures that exhibit kissing hairpin motifs in an arbitrarily nested fashion, requiring Oðn 4 Þ time and Oðn 2 Þ space. We added shape analysis, probabilities, different folding strategies and different dangling base models. The -cast option provides comparative prediction of pseudoknotted structures as in the RNACAST approach (Reeder and Giegerich, 2005) . A window mode was also included.

The program PALIKISS allows to predict pseudoknots, including kissing hairpins from aligned sequences. Being composed from the grammars and algebras of the other tools, it inherits all the features and options that make sense for it. *To whom correspondence should be addressed. 

All tools were augmented with utilities to compute folding energy or abstract shape for sequences that are provided with a structure from an outside source, in a way consistent with the tools' energy model. The graphical motif description tool LOCOMOTIF (Reeder et al., 2007) now uses modules from the RNA shapes studio. The KNOTINFRAME (Theis et al., 2008) tool that predicts -1 ribosomal frameshifts has been updated as well.

The Corona virus family contains a frameshift stimulation element (Baranov et al., 2005) , where the frameshift is facilitated by a slippery site together with either an H-type or a K-type pseudoknot. RFAM (Burge et al., 2013) holds the corresponding family model RF00507, although the tools of RFAM cannot explicitly model pseudoknots. PKISS (Theis et al., 2010) in -enforce mode reveals that for 11 family members, minimal free energy structures are H-types, another 11 are K-type pseudoknots and for only one member a purely nested structure has the best energy. The recently sequenced MERS genome (KF958702.1) is annotated with a homologous frameshift site, whereas the structure of the triggering element remains unclear. Structure prediction with PKISS for a 100 bp stretch downstream the slippery site attests a most stable K-type pseudoknot (see Fig. 2) . A second run of PKISS, this time in probability mode, shows that the shape class of this particular K-type pseudoknot has an overwhelming Boltzmann probability of 99%; leaving not much probability mass for any other shape classes.

Programmed ribosomal frameshifting in decoding the sars-cov genome

RNAalifold: improved consensus structure prediction for RNA alignments

Rfam 11.0: 10 years of RNA families

VARNA: Interactive drawing and editing of the RNA secondary structure

A discipline of dynamic programming over sequence data

Abstract shapes of RNA

Lost in folding space? Comparing four variants of the thermodynamic model for RNA secondary structure prediction

Improved RNA secondary structure prediction by maximizing expected pair accuracy

Design, implementation and evaluation of a practical pseudoknot folding algorithm based on thermodynamics

Consensus shapes: an alternative to the Sankoff algorithm for RNA consensus structure prediction

Locomotif: from graphical motif description to RNA motif search

Conservation and occurrence of trans-encoded sRNAs in the rhizobiales

Bellman's GAP -a language and compiler for dynamic programming in sequence analysis

Versatile and declarative dynamic programming using pair algebras

Knotinframe: prediction of -1 ribosomal frameshift events

Prediction of RNA secondary structure including kissing hairpin motifs

Structural analysis of aligned RNAs

Complete probabilistic analysis of RNA shapes

The RNA shapes studio