Nucleic Acids Research, 2003, Vol. 31, No. 13 3429-3431
© 2003 Oxford University Press
Vienna RNA secondary structure server
Institut für Theoretische Chemie und Molekulare Strukturbiologie, Universität Wien, Währingerstraße 17, A-1090 Wien, Austria
*Tel: +43 1 4277 52738; Fax: +43 1 4277 52793; Email: ivo{at}tbi.univie.ac.at
Received February 15, 2003; Revised and Accepted April 5, 2003
| ABSTRACT |
|---|
|
|
|---|
The Vienna RNA secondary structure server provides a web interface to the most frequently used functions of the Vienna RNA software package for the analysis of RNA secondary structures. It currently offers prediction of secondary structure from a single sequence, prediction of the consensus secondary structure for a set of aligned sequences and the design of sequences that will fold into a predefined structure. All three services can be accessed via the Vienna RNA web server at http://rna.tbi.univie.ac.at/.
| INTRODUCTION |
|---|
|
|
|---|
Biomolecules exhibit a close interplay between structure and function. Therefore the growing number of RNA molecules with complex functions, beyond that of encoding proteins, has brought increased demand for RNA structure prediction methods. While prediction of tertiary structure is usually infeasible, the area of RNA secondary structures is an example where computational methods have been highly successful.
The first practical dynamic programming algorithms to predict the optimal secondary structure of an RNA sequence date back over 20 years (1). Since then they have been extended to allow prediction of suboptimal structures (2,3) and thermodynamic ensembles (4), which allow to assign a confidence level or well definedness to the predictions (5).
Recently, several methods have addressed the problem of predicting a consensus structure for a group of related RNA sequences (611). Such conserved structures are of particular interest, since conservation of structure in spite of sequence variation implies that the structure must be functionally important. By enhancing energy rules with sequence covariation these methods also obtain much better prediction accuracies.
The Vienna RNA package (12) is a free software package that implements a variety of algorithms for the prediction and analysis of RNA secondary structures. The package is, however, strongly geared toward Unix command-line users and programmers. For the less computer savvy, or occasional user, it provides neither a point-and-click graphical user interface nor even pre-compiled binaries.
The Vienna RNA web site tries to address these shortcomings by offering access to the most popular features via an easy to use web interface. It consists of three CGI scripts equivalent to the RNAfold, RNAalifold and RNAinverse command line programs, respectively. While the servers have to limit request sizes for performance reasons, they return for each request an equivalent command line invocation. This makes it easier for users to make the transition to locally installed software, should their requirements exceed the limits of the web service.
| THE RNAfold SERVER |
|---|
|
|
|---|
Of the three services, the RNAfold server provides both the most basic and most widely used function. Input consists of a single sequence that has to be typed or pasted into a text field of the input form.
In the simplest case, the server predicts only the minimum free energy (mfe) structure of a single sequence using the classic algorithm of Zuker and Stiegler (1). In addition to mfe folding the server can calculate equilibrium base pairing probabilities via John McCaskill's partition function algorithm (4).
By default the RNA energy parameters of the Turner group (13) are used, but single stranded DNA sequences can be handled as well, by selecting the DNA parameter set provided by John SantaLucia (14).
The fold server output consists of a static html page presenting the predicted mfe structure as a string in bracket notation and links to the plots generated for visualization. Three types of plots can be produced. Firstly, the predicted mfe structure is plotted as a conventional secondary structure graph using the naview layout method (15). The pair probabilities can be visualized in a so-called dot plot: on a square grid of nxn we draw for each possible pair (i, j) a box with area proportional to its probability. Finally, we produce a mountain plot depicting both the predicted mfe and pair probabilities. A mountain plot is an xy-graph that plots the number of base pairs enclosing a sequence position (for pair probabilities the average number of enclosing pairs). See Figure 1 for examples of all three representations.
|
Secondary structure drawing and dot plots are always produced in Postscript format. Postscript is used not only because it gives the highest print quality, but also because it allows the actual data to be embedded in the file, e.g. all pair probabilities are contained in the dot plot in an easy to parse format. On the other hand, Postscript files cannot be used for inline images on web pages and require additional software for viewing (e.g. gsview, http://www.ghostscript.com/).
A suitable alternative is the new standard for Scalable Vector Graphics, SVG (http://www.w3.org/Graphics/SVG). Users with SVG enabled browsers (typically through the use of Adobe's SVG plugin, http://www.adobe.com/svg/) can request structure drawings in SVG, which allows some interactivity such as toggling annotation. Currently the server accepts sequences up to a maximum length of 4000 nt, sequences up to 300 nt will be processed immediately while longer jobs are submitted to a batch queue, in which case the user is notified by email after completion.
| THE Alifold SERVER |
|---|
|
|
|---|
The Alifold service predicts the consensus secondary structure for a set of aligned RNA or DNA sequences by using modified dynamic programming algorithms that add a covariance term to the standard energy model (11), again it supports prediction of mfe structures and pair probabilities. Usage is almost identical to that of the RNAfold service. Instead of typing an input sequence, a precomputed sequence alignment is uploaded via the input form. Currently, only alignments in Clustal format are accepted. The server restricts both the size of the upload and the length of the alignment, current limits being 10 Kb and 2000 nt, respectively.
Results are again visualized in Postscript plots that are enhanced by information on sequence variation. In the structure drawings mutations supporting the predicted structure are marked by circles, in the dot plots and mountain plots, color is used to indicate the number of different pair types. Examples and detailed explanation of these representations can be found on the online help page (http://www.tbi.univie.ac.at/~ivo/RNA/alifoldcgi.html).
| THE INVERSE FOLD SERVER |
|---|
|
|
|---|
Finding sequences that fold into a predefined structure is the inverse of structure prediction problem. Often it is useful to design such sequences, e.g. in order to experimentally test an hypothesis about functional structures. While this is often done manually for very short sequences, it quickly becomes tedious and error prone.
Our inverse folding service treats sequence design as an optimization problem in sequence space that is solved heuristically (12). There are again two variants based on mfe and partition function folding. In the first case we minimize the dissimilarity between the predicted mfe structure and the desired target structure. In the second case we optimize the frequency of the target structure in the thermodynamic ensemble. While the mfe optimization typically yields sequences that are marginally stable, i.e. have many alternative foldings, optimization via the partition function produces sequences with a very strong preference for the target structure.
Input consists simply of the desired structure in bracket notation. The maximum structure length is currently 100 nt. The time needed for the search varies widely depending on the ubiquity of the target structure. Most valid secondary structure strings never occur as mfe structure of some sequence (i.e. many sequence design problems have no solution), while some others are extremely common (for example see 16). Conversely, the number of search steps performed by the algorithm is a good indicator for the frequency of a structure in sequence space.
| FUTURE PLANS |
|---|
|
|
|---|
The Vienna RNA secondary structure server presented here provides only basic access to a subset of the functions in the Vienna RNA software package. Nevertheless they provide a convenient interface for users that need RNA structure prediction only occasionally and a shallow learning curve for those new to the field.
Work is underway to further improve the visualization of the results, e.g. by producing structure drawings annotated with various measures of well-definedness. As SVG enabled browsers become more widespread, a combination of SVG graphics and client side javascript should allow users to explore the predicted structures interactively.
The output web page produced by the server is designed for the interactive user and thus is not ideal for automatic parsing and further processing of the results. To facilitate such interoperation with other programs and web services we plan to offer input and output in a standardized data exchange format. A promising candidate for this is the recently proposed RNAML format (17), an XML based language for the storage of information on RNA sequence and structure.
While the server currently runs on a somewhat dated dual Pentium II 450 MHz machine, the use of a batch queuing system allows jobs to be distributed to other machines should that become necessary.
| ACKNOWLEDGEMENTS |
|---|
This work is supported by the Austrian Fonds zur Förderung der Wissenschaftlichen Forschung, Projects FWF 15893 and P-13545-MAT.
| REFERENCES |
|---|
|
|
|---|
- Zuker,M. and Stiegler,P. (1981) Optimal computer folding of larger RNA sequences using thermodynamics and auxiliary information. Nucleic Acids Res., 9, 133148.
[Abstract/Free Full Text] - Zuker,M. (1989) The use of dynamic programming algorithms in RNA secondary structure prediction. In Waterman,M.S. (ed.), Mathematical Methods for DNA Sequences, CRC Press, Boca Raton, FL, pp. 159184.
- Wuchty,S., Fontana,W., Hofacker,I.L. and Schuster,P. (1999) Complete suboptimal folding of RNA and the stability of secondary structures. Biopolymers, 49, 145165.[CrossRef][ISI][Medline]
- McCaskill,J.S. (1990) The equilibrium partition function and base pair binding probabilities for RNA secondary structure. Biopolymers, 29, 11051119.[CrossRef][ISI][Medline]
- Zuker,M. and Jacobson,A.B. (1995) Well-determined regions in RNA secondary structure prediction: analysis of small subunit ribosomal RNA. Nucleic Acids Res., 23, 27912798.
[Abstract/Free Full Text] - Gorodkin,J., Heyer,L.J. and Stormo,G.D. (1997) Finding the most significant common sequence and structure motifs in a set of RNA sequences. Nucleic Acids Res., 25, 37243732.
[Abstract/Free Full Text] - Hofacker,I.L., Fekete,M., Flamm,C., Huynen,M.A., Rauscher,S., Stolorz,P.E. and Stadler,P.F. (1998) Automatic detection of conserved RNA structure elements in complete RNA virus genomes. Nucleic Acids Res., 26, 38253836.
[Abstract/Free Full Text] - Lück,R., Graf,S. and Steger,G. (1999) Construct: a tool for thermodynamic controlled prediction of conserved secondary structure. Nucleic Acids Res., 27, 42084217.
[Abstract/Free Full Text] - Juan,V. and Wilson,C. (1999) RNA secondary structure prediction based on free energy and phylogenetic analysis. J. Mol. Biol., 289, 935947.[CrossRef][ISI][Medline]
- Knudsen,B. and Hein,J. (1999) RNA secondary structure prediction using stochastic context-free grammars and evolutionary history. Bioinformatics, 15, 446454.
[Abstract/Free Full Text] - Hofacker,I.L., Fekete,M. and Stadler,P.F. (2002) Secondary structure prediction for aligned RNA sequences. J. Mol. Biol., 319, 10591066.[CrossRef][ISI][Medline]
- Hofacker,I.L., Fontana,W., Stadler,P.F., Bonhoeffer,S., Tacker,M. and Schuster,P. (1994) Fast folding and comparison of RNA secondary structures. Monatsh. Chem., 125, 167188.[CrossRef]
- Mathews,D., Sabina,J., Zucker,M. and Turner,H. (1999) Expanded sequence dependence of thermodynamic parameters provides robust prediction of RNA secondary structure. J. Mol. Biol., 288, 911940.[CrossRef][ISI][Medline]
- SantaLucia,J. Jr. (1998) A unified view of polymer, dumbbell, and oligonucleotide DNA nearest-neighbor thermodynamics. Proc. Natl Acad. Sci. USA, 95, 14601465.
[Abstract/Free Full Text] - Bruccoleri,R.E. and Heinrich,G. (1988) An improved algorithm for nucleic acid secondary structure display. CABIOS, 4, 167173.
- Schuster,P., Fontana,W., Stadler,P.F. and Hofacker,I.L. (1994) From sequences to shapes and back: a case study in RNA secondary structures. Proc. Royal Society London B, 255, 279284.
- Waugh,A., Gendron,P., Altman,R., Brown,J., Case,D., Gautheret,D., Harvey,S., Leontis,N., Westbrook,J., Westhof,E. et al. (2002) RNAML: a standard syntax for exchanging RNA information. RNA, 8, 707717.[Abstract]
This article has been cited by other articles:
![]() |
A. Loya, L. Pnueli, Y. Yosefzon, Y. Wexler, M. Ziv-Ukelson, and Y. Arava The 3'-UTR mediates the cellular localization of an mRNA encoding a short plasma membrane protein RNA, July 1, 2008; 14(7): 1352 - 1365. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Xu, X. Zhou, and W. Zhang MicroRNA prediction with a novel ranking algorithm based on random walks Bioinformatics, July 1, 2008; 24(13): i50 - i58. [Abstract] [PDF] |
||||
![]() |
A. R. Gruber, R. Lorenz, S. H. Bernhart, R. Neubock, and I. L. Hofacker The Vienna RNA Websuite Nucleic Acids Res., July 1, 2008; 36(suppl_2): W70 - W74. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Asai, H. Kiryu, M. Hamada, Y. Tabei, K. Sato, H. Matsui, Y. Sakakibara, G. Terai, and T. Mituyama Software.ncrna.org: web servers for analyses of RNA sequences Nucleic Acids Res., July 1, 2008; 36(suppl_2): W75 - W78. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Knudsen and G. Caetano-Anolles NOBAI: a web server for character coding of geometrical and statistical features in RNA structure Nucleic Acids Res., July 1, 2008; 36(suppl_2): W85 - W90. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Moretti, A. Wilm, D. G. Higgins, I. Xenarios, and C. Notredame R-Coffee: a web server for accurately aligning noncoding RNA sequences Nucleic Acids Res., July 1, 2008; 36(suppl_2): W10 - W13. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Torarinsson and S. Lindgreen WAR: Webserver for aligning structural RNAs Nucleic Acids Res., July 1, 2008; 36(suppl_2): W79 - W84. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. M. Chubiz and C. V. Rao Computational design of orthogonal ribosomes Nucleic Acids Res., June 3, 2008; (2008) gkn354v1. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang miRDB: A microRNA target prediction and functional annotation database with a wiki interface RNA, June 1, 2008; 14(6): 1012 - 1017. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Ding, S. Sharma, P. Chalasani, V. V. Demidov, N. E. Broude, and N. V. Dokholyan Ab initio RNA folding by discrete molecular dynamics: From structure prediction to folding mechanisms RNA, June 1, 2008; 14(6): 1164 - 1173. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. A. Glazov, P. A. Cottee, W. C. Barris, R. J. Moore, B. P. Dalrymple, and M. L. Tizard A microRNA catalog of the developing chicken embryo identified by a deep sequencing approach Genome Res., June 1, 2008; 18(6): 957 - 964. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Koch and G. Fuellen A review of bioinformatics education in Germany Brief Bioinform, May 1, 2008; 9(3): 232 - 242. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Wilm, D. G. Higgins, and C. Notredame R-Coffee: a method for multiple alignment of non-coding RNA Nucleic Acids Res., May 1, 2008; 36(9): e52 - e52. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Simmonds, I. Karakasiliotis, D. Bailey, Y. Chaudhry, D. J. Evans, and I. G. Goodfellow Bioinformatic and functional analysis of RNA secondary structure elements among different genera of human and animal caliciviruses Nucleic Acids Res., May 1, 2008; 36(8): 2530 - 2546. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. O. Harmanci, G. Sharma, and D. H. Mathews PARTS: Probabilistic Alignment for RNA joinT Secondary structure prediction Nucleic Acids Res., April 1, 2008; 36(7): 2406 - 2417. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. D. Morin, M. D. O'Connor, M. Griffith, F. Kuchenbauer, A. Delaney, A.-L. Prabhu, Y. Zhao, H. McDonald, T. Zeng, M. Hirst, et al. Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells Genome Res., April 1, 2008; 18(4): 610 - 621. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. D. Morin, G. Aksay, E. Dolgosheina, H. A. Ebhardt, V. Magrini, E. R. Mardis, S. C. Sahinalp, and P. J. Unrau Comparative analysis of the small RNA transcriptomes of Pinus contorta and Oryza sativa Genome Res., April 1, 2008; 18(4): 571 - 584. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Huppi, N. Volfovsky, T. Runfola, T. L. Jones, M. Mackiewicz, S. E. Martin, J. F. Mushinski, R. Stephens, and N. J. Caplen The Identification of MicroRNAs in a Genomically Unstable Region of Human Chromosome 8q24 Mol. Cancer Res., February 1, 2008; 6(2): 212 - 221. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Wang and I. M. El Naqa Prediction of both conserved and nonconserved microRNA targets in animals Bioinformatics, February 1, 2008; 24(3): 325 - 332. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Kiryu, T. Kin, and K. Asai Rfold: an exact algorithm for computing local base pairing probabilities Bioinformatics, February 1, 2008; 24(3): 367 - 373. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Care, C. Bignon, M. C. Pelissier, E. Blanc, B. Canard, and B. Coutard The translation of recombinant proteins in E. coli can be improved by in silico generating and screening random libraries of a -70/+96 mRNA region with respect to the translation initiation codon Nucleic Acids Res., January 17, 2008; 36(1): e6 - e6. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Hertel, I. L. Hofacker, and P. F. Stadler SnoReport: computational identification of snoRNAs with unknown targets Bioinformatics, January 15, 2008; 24(2): 158 - 164. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. H. Buck, J. Santoyo-Lopez, K. A. Robertson, D. S. Kumar, M. Reczko, and P. Ghazal Discrete Clusters of Virus-Encoded MicroRNAs Are Associated with Complementary Strands of the Genome and the 7.2-Kilobase Stable Intron in Murine Cytomegalovirus J. Virol., December 15, 2007; 81(24): 13761 - 13770. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Lindgreen, P. P. Gardner, and A. Krogh MASTR: multiple alignment and structure prediction of non-coding RNAs using simulated annealing Bioinformatics, December 15, 2007; 23(24): 3304 - 3311. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Yousef, S. Jung, A. V. Kossenkov, L. C. Showe, and M. K. Showe Naive Bayes for microRNA target predictions machine learning for microRNA targets Bioinformatics, November 15, 2007; 23(22): 2987 - 2992. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. M. Meyer A practical guide to the art of RNA gene prediction Brief Bioinform, November 1, 2007; 8(6): 396 - 414. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. S. Andersen, A. Lind-Thomsen, B. Knudsen, S. E. Kristensen, J. H. Havgaard, E. Torarinsson, N. Larsen, C. Zwieb, P. Sestoft, J. Kjems, et al. Semiautomated improvement of RNA alignments RNA, November 1, 2007; 13(11): 1850 - 1859. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Matousek, L. Orctova, J. Ptacek, J. Patzak, P. Dedic, G. Steger, and D. Riesner Experimental Transmission of Pospiviroid Populations to Weed Species Characteristic of Potato and Hop Fields J. Virol., November 1, 2007; 81(21): 11891 - 11899. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Shao, C. Y. Chan, A. Maliyekkel, C. E. Lawrence, I. B. Roninson, and Y. Ding Effect of target secondary structure on RNAi efficiency RNA, October 1, 2007; 13(10): 1631 - 1640. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Freyhult, V. Moulton, and P. Clote Boltzmann probability of RNA structural neighbors and riboswitch detection Bioinformatics, August 15, 2007; 23(16): 2054 - 2062. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. R. Gruber, R. Neubock, I. L. Hofacker, and S. Washietl The RNAz web server: prediction of thermodynamically stable and evolutionarily conserved RNA structures Nucleic Acids Res., July 13, 2007; 35(suppl_2): W335 - W338. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Khaladkar, V. Bellofatto, J. T. L. Wang, B. Tian, and B. A. Shapiro RADAR: a web server for RNA data analysis and research Nucleic Acids Res., July 13, 2007; 35(suppl_2): W300 - W304. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Jiang, H. Wu, W. Wang, W. Ma, X. Sun, and Z. Lu MiPred: classification of real and pseudo microRNA precursors using random forest prediction model with combined features Nucleic Acids Res., July 13, 2007; 35(suppl_2): W339 - W344. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Shu, X. Bo, Z. Zheng, and S. Wang RSRE: RNA structural robustness evaluator Nucleic Acids Res., July 13, 2007; 35(suppl_2): W314 - W319. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. T. Souza, M. R. M. Santos, F. M. Lima, N. M. El-Sayed, P. J. Myler, J. C. Ruiz, and J. F. da Silveira New Trypanosoma cruzi Repeated Element That Shows Site Specificity for Insertion Eukaryot. Cell, July 1, 2007; 6(7): 1228 - 1238. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Kiryu, Y. Tabei, T. Kin, and K. Asai Murlet: a practical multiple alignment tool for structural RNA sequences Bioinformatics, July 1, 2007; 23(13): 1588 - 1598. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. U. T. Hellen and S. de Breyne A Distinct Group of Hepacivirus/Pestivirus-Like Internal Ribosomal Entry Sites in Members of Diverse Picornavirus Genera: Evidence for Modular Exchange of Functional Noncoding RNA Elements by Recombination J. Virol., June 1, 2007; 81(11): 5850 - 5863. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. S. Ng and S. K. Mishra De novo SVM classification of precursor microRNAs from genomic pseudo hairpins using global and intrinsic folding measures Bioinformatics, June 1, 2007; 23(11): 1321 - 1330. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Liu, L. Li, J. J. Millership, H. Kang, J. L. Leibowitz, and D. P. Giedroc A U-turn motif-containing stem-loop in the coronavirus 5' untranslated region plays a functional role in replication RNA, May 1, 2007; 13(5): 763 - 780. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Kim, H. H. Gan, and T. Schlick A computational proposal for designing structured RNA pools for in vitro selection of RNAs RNA, April 1, 2007; 13(4): 478 - 492. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Kiryu, T. Kin, and K. Asai Robust prediction of consensus secondary structures using averaged base pairing probability matrices Bioinformatics, February 15, 2007; 23(4): 434 - 441. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. NG Kwang Loong and S. K. Mishra Unique folding of precursor microRNAs: Quantitative evidence and implications for de novo identification RNA, February 1, 2007; 13(2): 170 - 187. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Tian, Z. Pan, and J. Y. Lee Widespread mRNA polyadenylation events in introns indicate dynamic interplay between polyadenylation and splicing Genome Res., February 1, 2007; 17(2): 156 - 165. [Abstract] [Full Text] [PDF] |
||||








