410 research outputs found
The Vienna RNA Websuite
The Vienna RNA Websuite is a comprehensive collection of tools for folding, design and analysis of RNA sequences. It provides a web interface to the most commonly used programs of the Vienna RNA package. Among them, we find folding of single and aligned sequences, prediction of RNA–RNA interactions, and design of sequences with a given structure. Additionally, we provide analysis of folding landscapes using the barriers program and structural RNA alignments using LocARNA. The web server together with software packages for download is freely accessible at http://rna.tbi.univie.ac.at/
Analytical description of finite size effects for RNA secondary structures
The ensemble of RNA secondary structures of uniform sequences is studied
analytically. We calculate the partition function for very long sequences and
discuss how the cross-over length, beyond which asymptotic scaling laws apply,
depends on thermodynamic parameters. For realistic choices of parameters this
length can be much longer than natural RNA molecules. This has to be taken into
account when applying asymptotic theory to interpret experiments or numerical
results.Comment: 10 pages, 13 figures, published in Phys. Rev.
Discriminatory power of RNA family models
Motivation: RNA family models group nucleotide sequences that share a common biological function. These models can be used to find new sequences belonging to the same family. To succeed in this task, a model needs to exhibit high sensitivity as well as high specificity. As model construction is guided by a manual process, a number of problems can occur, such as the introduction of more than one model for the same family or poorly constructed models. We explore the Rfam database to discover such problems
R-Coffee: a web server for accurately aligning noncoding RNA sequences
The R-Coffee web server produces highly accurate multiple alignments of noncoding RNA (ncRNA) sequences, taking into account predicted secondary structures. R-Coffee uses a novel algorithm recently incorporated in the T-Coffee package. R-Coffee works along the same lines as T-Coffee: it uses pairwise or multiple sequence alignment (MSA) methods to compute a primary library of input alignments. The program then computes an MSA highly consistent with both the alignments contained in the library and the secondary structures associated with the sequences. The secondary structures are predicted using RNAplfold. The server provides two modes. The slow/accurate mode is restricted to small datasets (less than 5 sequences less than 150 nucleotides) and combines R-Coffee with Consan, a very accurate pairwise RNA alignment method. For larger datasets a fast method can be used (RM-Coffee mode), that uses R-Coffee to combine the output of the three packages which combines the outputs from programs found to perform best on RNA (MUSCLE, MAFFT and ProbConsRNA). Our BRAliBase benchmarks indicate that the R-Coffee/Consan combination is one of the best ncRNA alignment methods for short sequences, while the RM-Coffee gives comparable results on longer sequences. The R-Coffee web server is available at http://www.tcoffee.org
R-Coffee: a web server for accurately aligning noncoding RNA sequences
The R-Coffee web server produces highly accurate multiple alignments of noncoding RNA (ncRNA) sequences, taking into account predicted secondary structures. R-Coffee uses a novel algorithm recently incorporated in the T-Coffee package. R-Coffee works along the same lines as T-Coffee: it uses pairwise or multiple sequence alignment (MSA) methods to compute a primary library of input alignments. The program then computes an MSA highly consistent with both the alignments contained in the library and the secondary structures associated with the sequences. The secondary structures are predicted using RNAplfold. The server provides two modes. The slow/accurate mode is restricted to small datasets (less than 5 sequences less than 150 nucleotides) and combines R-Coffee with Consan, a very accurate pairwise RNA alignment method. For larger datasets a fast method can be used (RM-Coffee mode), that uses R-Coffee to combine the output of the three packages which combines the outputs from programs found to perform best on RNA (MUSCLE, MAFFT and ProbConsRNA). Our BRAliBase benchmarks indicate that the R-Coffee/Consan combination is one of the best ncRNA alignment methods for short sequences, while the RM-Coffee gives comparable results on longer sequences. The R-Coffee web server is available at http://www.tcoffee.org
RNA denaturation: excluded volume, pseudoknots and transition scenarios
A lattice model of RNA denaturation which fully accounts for the excluded
volume effects among nucleotides is proposed. A numerical study shows that
interactions forming pseudoknots must be included in order to get a sharp
continuous transition. Otherwise a smooth crossover occurs from the swollen
linear polymer behavior to highly ramified, almost compact conformations with
secondary structures. In the latter scenario, which is appropriate when these
structures are much more stable than pseudoknot links, probability
distributions for the lengths of both loops and main branches obey scaling with
nonclassical exponents.Comment: 4 pages 3 figure
Conserved RNA secondary structures in viral genomes: a survey
The genomes of RNA viruses often carry conserved RNA structures that perform vital functions during the life cycle of the virus. Such structures can be detected using a combination of structure prediction and co-variation analysis. Here we present results from pilot studies on a variety of viral families performed during bioinformatics computer lab courses in past years
Prediction of locally stable RNA secondary structures for genome-wide surveys
Motivation: Recently novel classes of functional RNAs, most prominently the miRNAs have been discovered, strongly suggesting that further types of functional RNAs are still hidden in the recently completed genomic DNA sequences. Only few techniques are known, however, to survey genomes for such RNA genes. When sufficiently similar sequences are not available for comparative approaches the only known remedy is to search directly for structural features.
Results: We present here efficient algorithms for computing locally stable RNA structures at genome-wide scales. Both the minimum energy structure and the complete matrix of base pairing probabilities can be computed in (N × L2) time and (N + L2) memory in terms of the length N of the genome and the size L of the largest secondary structure motifs of interest. In practice, the 100 Mb of the complete genome of Caenorhabditis elegans can be folded within about half a day on a modern PC with a search depth of L = 100. This is sufficient example for a survey for miRNAs
Zero Temperature Properties of RNA Secondary Structures
We analyze different microscopic RNA models at zero temperature. We discuss
both the most simple model, that suffers a large degeneracy of the ground
state, and models in which the degeneracy has been remove, in a more or less
severe manner. We calculate low-energy density of states using a coupling
perturbing method, where the ground state of a modified Hamiltonian, that
repels the original ground state, is determined. We evaluate scaling exponents
starting from measurements of overlaps and energy differences. In the case of
models without accidental degeneracy of the ground state we are able to clearly
establish the existence of a glassy phase with .Comment: 20 pages including 9 eps figure
- …