Search

Bacterial genomes encode a plethora of small RNAs (sRNAs), which are heterogeneous in size, structure and function. Most sRNAs act as post-transcriptional regulators by means of specific base pairing interactions with the 5'-untranslated region of mRNA transcripts, thereby modifying the stability of the target transcript and/or its ability to be translated. Here, we present RNApredator, a web server for the prediction of sRNA targets. The user can choose from a set of over 2155 genomes and plasmids from 1183 bacterial species. RNApredator then uses a dynamic programming approach, RNAplex, to compute putative targets. Compared to web servers with a similar task, RNApredator takes the accessibility of the target during the target search into account, improving the specificity of the predictions. Furthermore, enrichment in Gene Ontology terms, cellular pathways as well as changes in accessibilities along the target sequence can be done in fully automated post-processing steps. The predictive performance of the underlying dynamic programming approach RNAplex is similar to that of more complex methods, but needs at least three orders of magnitude less time to complete. RNApredator is available at http://rna.tbi.univie.ac.at/RNApredator.

Motivation: Small nucleolar RNAs (snoRNAs) are an abundant class of non-coding RNAs with a wide variety of cellular functions including chemical modification of RNA, telomere maintanance, pre-rRNA processing and regulatory activities in alternative splicing. The main role of box C/D snoRNAs is to determine the targets for 2'-O-ribose methylation, which is important for rRNA maturation and splicing regulation of some mRNAs. The targets are still unknown, however, for many 'orphan' snoRNAs. While a fast and efficient target predictor for box H/ACA snoRNAs is available, no comparable tool exists for box C/D snoRNAs, even though they bind to their targets in a much less complex manner., Results: PLEXY is a dynamic programming algorithm that computes thermodynamically optimal interactions of a box C/D snoRNA with a putative target RNA. Implemented as scanner for large input sequences and equipped with filters on the duplex structure, PLEXY is an efficient and reliable tool for the prediction of box C/D snoRNA target sites., Availability: The perl script PLEXY is freely available at http://www.bioinf.uni-leipzig.de/Software/PLEXY.

A synergistic combination of two next-generation sequencing platforms with a detailed comparative BAC physical contig map provided a cost-effective assembly of the genome sequence of the domestic turkey (Meleagris gallopavo). Heterozygosity of the sequenced source genome allowed discovery of more than 600,000 high quality single nucleotide variants. Despite this heterozygosity, the current genome assembly (∼1.1 Gb) includes 917 Mb of sequence assigned to specific turkey chromosomes. Annotation identified nearly 16,000 genes, with 15,093 recognized as protein coding and 611 as non-coding RNA genes. Comparative analysis of the turkey, chicken, and zebra finch genomes, and comparing avian to mammalian species, supports the characteristic stability of avian genomes and identifies genes unique to the avian lineage. Clear differences are seen in number and variety of genes of the avian immune system where expansions and novel genes are less frequent than examples of gene loss. The turkey genome sequence provides resources to further understand the evolution of vertebrate genomes and genetic variation underlying economically important quantitative traits in poultry. This integrated approach may be a model for providing both gene and chromosome level assemblies of other species with agricultural, ecological, and evolutionary interest., Competing Interests: The authors have declared that no competing interests exist.

Motivation: Small nucleolar RNAs are an abundant class of non-coding RNAs that guide chemical modifications of rRNAs, snRNAs and some mRNAs. In the case of many 'orphan' snoRNAs, the targeted nucleotides remain unknown, however. The box H/ACA subclass determines uridine residues that are to be converted into pseudouridines via specific complementary binding in a well-defined secondary structure configuration that is outside the scope of common RNA (co-)folding algorithms., Results: RNAsnoop implements a dynamic programming algorithm that computes thermodynamically optimal H/ACA-RNA interactions in an efficient scanning variant. Complemented by an support vector machine (SVM)-based machine learning approach to distinguish true binding sites from spurious solutions and a system to evaluate comparative information, it presents an efficient and reliable tool for the prediction of H/ACA snoRNA target sites. We apply RNAsnoop to identify the snoRNAs that are responsible for several of the remaining 'orphan' pseudouridine modifications in human rRNAs, and we assign a target to one of the five orphan H/ACA snoRNAs in Drosophila., Availability: The C source code of RNAsnoop is freely available at http://www.tbi.univie.ac.at/ -htafer/RNAsnoop

RNA interference, mediated by small interfering RNAs (siRNAs), is a powerful tool for investigation of gene functions and it is increasingly being used as a therapeutic agent. However, not all siRNAs are equally potent - although simple rules for the selection of good siRNAs were proposed early on, siRNAs are still plagued with widely fluctuating efficiency. Recently, new design tools that incorporate both the structural features of the targeted RNAs and the sequence features of the siRNAs have substantially improved the efficacy of siRNAs. In this chapter, we present the algorithms behind these accessibility-aided tools and show how to design efficient siRNAs with their help.

In RNAi (RNA interference), siRNAs (small interfering RNAs) are loaded into the RISC (RNA-induced silencing complex), which then mediates endonucleolytic cleavage of complementary target RNAs. Although RNAi has become one of the most powerful tools in molecular biology to assess gene function, there remains a great number of ineffective siRNAs. It is already known that the assembly and activation of RISC is a crucial determinant of RNAi activity, but downstream effects such as target accessibility have not been analysed extensively. Therefore we assessed the effect of target site accessibility and found that it significantly improves the potency of siRNAs. Similarly, miRNAs (microRNAs) act by repressing protein synthesis through imperfect base-pairing to the 3'-UTR (untranslated region) of target mRNAs. We found that predicted target sites reside in regions of high accessibility and tested whether this criterion could be used in the search of functional miRNA targets. In addition, we performed reporter gene assays to test whether accessibility correlates with measured mRNA suppression levels. The results of our initial study suggest that secondary structures might add a so far underrepresented layer of complexity in the recognition of RNA targets by miRNAs.

Motivation: Regulatory RNAs often unfold their action via RNA-RNA interaction. Transcriptional gene silencing by means of siRNAs and miRNA as well as snoRNA directed RNA editing rely on this mechanism. Additionally ncRNA regulation in bacteria is mainly based upon RNA duplex formation. Finding putative target sites for newly discovered ncRNAs is a lengthy task as tools for cofolding RNA molecules like RNAcofold and RNAup are too slow for genome-wide search. Tools like RNAhybrid that neglects intramolecular interactions have runtimes proportional to O(m x n), albeit with a large prefactor. Still in many cases the need for even faster methods exists., Results: We present a new program, RNAplex, especially designed to quickly find possible hybridization sites for a query RNA in large RNA databases. RNAplex uses a slightly different energy model which reduces the computational time by a factor 10-27 compared to RNAhybrid. In addition a length penalty allows to focus the target search on short highly stable interactions., Availability: RNAplex can be downloaded at http://www.tbi.univie.ac.at/~htafer/, Supplementary Information: Supplementary data are available at Bioinformatics online.

Small-interfering RNAs (siRNAs) assemble into RISC, the RNA-induced silencing complex, which cleaves complementary mRNAs. Despite their fluctuating efficacy, siRNAs are widely used to assess gene function. Although this limitation could be ascribed, in part, to variations in the assembly and activation of RISC, downstream events in the RNA interference (RNAi) pathway, such as target site accessibility, have so far not been investigated extensively. In this study we present a comprehensive analysis of target RNA structure effects on RNAi by computing the accessibility of the target site for interaction with the siRNA. Based on our observations, we developed a novel siRNA design tool, RNAxs, by combining known siRNA functionality criteria with target site accessibility. We calibrated our method on two data sets comprising 573 siRNAs for 38 genes, and tested it on an independent set of 360 siRNAs targeting four additional genes. Overall, RNAxs proves to be a robust siRNA selection tool that substantially improves the prediction of highly efficient siRNAs.

RNA interference (RNAi) has become an important tool to study and utilize gene silencing by introducing short interfering RNA (siRNA). In order to predict the most efficient siRNAs, a new software tool, RNA Workbench (RNAWB), has been designed and is freely available (after registration) on http://www.rnaworkbench.com. In addition to the standard selection rules, RNAWB includes the possibility of statistical analyses of the applied selection rules (criteria). The role of RNA secondary structures in the RNA interference process as well as the application of sequence rules are discussed to show the applicability of the software.

7SK RNA is a highly abundant noncoding RNA in mammalian cells whose function in transcriptional regulation has only recently been elucidated. Despite its highly conserved sequence throughout vertebrates, all attempts to discover 7SK RNA homologues in invertebrate species have failed so far. Here we report on a combined experimental and computational survey that succeeded in discovering 7SK RNAs in most of the major deuterostome clades and in two protostome phyla: mollusks and annelids. Despite major efforts, no candidates were found in any of the many available ecdysozoan genomes, however. The additional sequence data confirm the evolutionary conservation and hence functional importance of the previously described 3' and 5' stem-loop motifs, and provide evidence for a third, structurally well-conserved domain.

Background: Reliable prediction of RNA-RNA binding energies is crucial, e.g. for the understanding on RNAi, microRNA-mRNA binding and antisense interactions. The thermodynamics of such RNA-RNA interactions can be understood as the sum of two energy contributions: (1) the energy necessary to 'open' the binding site and (2) the energy gained from hybridization., Methods: We present an extension of the standard partition function approach to RNA secondary structures that computes the probabilities Pu[i, j] that a sequence interval [i, j] is unpaired., Results: Comparison with experimental data shows that Pu[i, j] can be applied as a significant determinant of local target site accessibility for RNA interference (RNAi). Furthermore, these quantities can be used to rigorously determine binding free energies of short oligomers to large mRNA targets. The resource consumption is comparable with a single partition function computation for the large target molecule. We can show that RNAi efficiency correlates well with the binding energies of siRNAs to their respective mRNA target., Availability: RNAup will be distributed as part of the Vienna RNA Package, www.tbi.univie.ac.at/~ivo/RNA/

Background: RNA has been recognized as a key player in cellular regulation in recent years. In many cases, non-coding RNAs exert their function by binding to other nucleic acids, as in the case of microRNAs and snoRNAs. The specificity of these interactions derives from the stability of inter-molecular base pairing. The accurate computational treatment of RNA-RNA binding therefore lies at the heart of target prediction algorithms., Methods: The standard dynamic programming algorithms for computing secondary structures of linear single-stranded RNA molecules are extended to the co-folding of two interacting RNAs., Results: We present a program, RNAcofold, that computes the hybridization energy and base pairing pattern of a pair of interacting RNA molecules. In contrast to earlier approaches, complex internal structures in both RNAs are fully taken into account. RNAcofold supports the calculation of the minimum energy structure and of a complete set of suboptimal structures in an energy band above the ground state. Furthermore, it provides an extension of McCaskill's partition function algorithm to compute base pairing probabilities, realistic interaction energies, and equilibrium concentrations of duplex structures.

On the basis of sequence-specific resonance assignments for the complete polypeptide backbone and most of the amino acid side chains by heteronuclear nuclear magnetic resonance (NMR) spectroscopy, the urea-unfolded form of the outer membrane protein X (OmpX) from Escherichia coli has been structurally characterized. (1)H-(1)H nuclear Overhauser effects (NOEs), dispersion of the chemical shifts, amide proton chemical shift temperature coefficients, amide proton exchange rates, and (15)N[(1)H]-NOEs show that OmpX in 8 M urea at pH 6.5 is globally unfolded, but adopts local nonrandom conformations in the polypeptide segments of residues 73-82 and 137-145. For these two regions, numerous medium-range and longer-range NOEs were observed, which were used as the input for structure calculations of these polypeptide segments with the program DYANA. The segment 73-82 forms a quite regular helical structure, with only loosely constrained amino acid side chains. In the segment 137-145, the tryptophan residue 140 forms the core of a small hydrophobic cluster. Both nonrandom structures are present with an abundance of about 25% of the protein molecules. The sequence-specific NMR assignment and the physicochemical characterization of urea-denatured OmpX presented in this paper are currently used as a platform for investigations of the folding mechanism of this integral membrane protein.