Consensus Secondary Structure Research Articles

The distributions, mechanisms, and structures of metabolite-binding riboswitches

Phylogenetic analyses revealed insights into the distribution of riboswitch classes in different microbial groups, and structural analyses led to updated aptamer structure models and insights into the mechanism of these non-coding RNA structures.

Complete sequence and secondary structure of the large subunit ribosomal RNA from the harmful unarmored dinoflagellateAkashiwo sanguinea

This is the first report of the complete DNA sequence of the gene encoding the ribosomal large subunit (LSU rDNA, 3336 bp) from the naked gymnodinioid dinoflagellate Akashiwo sanguinea. No introns were found in the LSU rDNA coding region and secondary structures were predicted for both the LSU and 5.8S rRNAs. The predicted LSU structure showed most of the features seen in the consensus secondary structure model proposed for the eukaryotic nuclear LSU rRNAs. However, six helices (C1_1, C1_2, C1_3, D10, D20_1 and H1_2) are not present in the A. sanguinea LSU structure. Particularly, the C branch area (or D2 domain), was extremely reduced compared to the eukaryotic consensus sequence due to nucleotide deletion. Phylogenetic resolution against 12 divergent (D) domains and cores in LSU rDNA showed that the D1, D2 and D12 domains were highly variable and could be used as genetic markers within low taxonomic levels, particularly in the gymnodinioid complex.

Robust prediction of consensus secondary structures using averaged base pairing probability matrices

Recent transcriptomic studies have revealed the existence of a considerable number of non-protein-coding RNA transcripts in higher eukaryotic cells. To investigate the functional roles of these transcripts, it is of great interest to find conserved secondary structures from multiple alignments on a genomic scale. Since multiple alignments are often created using alignment programs that neglect the special conservation patterns of RNA secondary structures for computational efficiency, alignment failures can cause potential risks of overlooking conserved stem structures. We investigated the dependence of the accuracy of secondary structure prediction on the quality of alignments. We compared three algorithms that maximize the expected accuracy of secondary structures as well as other frequently used algorithms. We found that one of our algorithms, called McCaskill-MEA, was more robust against alignment failures than others. The McCaskill-MEA method first computes the base pairing probability matrices for all the sequences in the alignment and then obtains the base pairing probability matrix of the alignment by averaging over these matrices. The consensus secondary structure is predicted from this matrix such that the expected accuracy of the prediction is maximized. We show that the McCaskill-MEA method performs better than other methods, particularly when the alignment quality is low and when the alignment consists of many sequences. Our model has a parameter that controls the sensitivity and specificity of predictions. We discussed the uses of that parameter for multi-step screening procedures to search for conserved secondary structures and for assigning confidence values to the predicted base pairs. The C++ source code that implements the McCaskill-MEA algorithm and the test dataset used in this paper are available at http://www.ncrna.org/papers/McCaskillMEA/. Supplementary data are available at Bioinformatics online.

Measuring covariation in RNA alignments: physical realism improves information measures

The importance of non-coding RNAs is becoming increasingly evident, and often the function of these molecules depends on the structure. It is common to use alignments of related RNA sequences to deduce the consensus secondary structure by detecting patterns of co-evolution. A central part of such an analysis is to measure covariation between two positions in an alignment. Here, we rank various measures ranging from simple mutual information to more advanced covariation measures. Mutual information is still used for secondary structure prediction, but the results of this study indicate which measures are useful. Incorporating more structural information by considering e.g. indels and stacking improves accuracy, suggesting that physically realistic measures yield improved predictions. This can be used to improve both current and future programs for secondary structure prediction. The best measure tested is the RNAalifold covariation measure modified to include stacking. Scripts, data and supplementary material can be found at http://www.binf.ku.dk/Stinus_covariation

Hairpins in a Haystack: recognizing microRNA precursors in comparative genomics data

Recently, genome-wide surveys for non-coding RNAs have provided evidence for tens of thousands of previously undescribed evolutionary conserved RNAs with distinctive secondary structures. The annotation of these putative ncRNAs, however, remains a difficult problem. Here we describe an SVM-based approach that, in conjunction with a non-stringent filter for consensus secondary structures, is capable of efficiently recognizing microRNA precursors in multiple sequence alignments. The software was applied to recent genome-wide RNAz surveys of mammals, urochordates, and nematodes. The program RNAmicro is available as source code and can be downloaded from http://www.bioinf.uni-leipzig/Software/RNAmicro.

Structural and functional insights into the B30.2/SPRY domain

The B30.2/SPRY domain is present in approximately 700 eukaryotic (approximately 150 human) proteins, including medically important proteins such as TRIM5alpha and Pyrin. Nonetheless, the functional role of this modular domain remained unclear. Here, we report the crystal structure of an SPRY-SOCS box family protein GUSTAVUS in complex with Elongins B and C, revealing a highly distorted two-layered beta-sandwich core structure of its B30.2/SPRY domain. Ensuing studies identified one end of the beta-sandwich as the surface interacting with an RNA helicase VASA with a 40 nM dissociation constant. The sequence variation in TRIM5alpha responsible for HIV-1 restriction and most of the mutations in Pyrin causing familial Mediterranean fever map on this surface, implicating the corresponding region in many B30.2/SPRY domains as the ligand-binding site. The amino acids lining the binding surface are highly variable among the B30.2/SPRY domains, suggesting that these domains are protein-interacting modules, which recognize a specific individual partner protein rather than a consensus sequence motif.

The vertebrate E1/U17 small nucleolar ribonucleoprotein particle

Each of the many different box H/ACA ribonucleoprotein particles (RNPs) present in eukaryotes and archaea consists of four common core proteins and one specific H/ACA small RNA, which bears the sequence elements H (ANANNA) and ACA. Most of the H/ACA RNPs are small nucleolar RNPs (snoRNPs), which are localized in nucleoli, and are one of the two major classes of snoRNPs. Most H/ACA RNPs direct pseudouridine synthesis in pre-rRNA and other RNAs. One H/ACA small nucleolar RNA (snoRNA), vertebrate E1/U17 (snR30 in yeast), is required for pre-rRNA cleavage processing that generates mature 18S rRNA. E1 snoRNA is encoded in introns of protein-coding genes, and the evidence suggests that human E1 RNA undergoes uridine insertional RNA editing. The vertebrate E1 RNA consensus secondary structure shows several features that are absent in other box H/ACA snoRNAs. The available UV-induced RNA-protein crosslinking results suggest that the E1 snoRNP is asymmetrical in vertebrate cells, in contrast to other H/ACA snoRNPs. The vertebrate E1 snoRNP in cells is surprisingly complex: (i) E1 RNA contacts directly and specifically several proteins which do not appear to be any of the H/ACA RNP four core proteins; and (ii) multiple E1 RNA sites are needed for E1 snoRNP formation, E1 RNA stability, and E1 RNA-protein direct interactions.

Identification and analysis of ribonuclease P and MRP RNA in a broad range of eukaryotes

RNases P and MRP are ribonucleoprotein complexes involved in tRNA and rRNA processing, respectively. The RNA subunits of these two enzymes are structurally related to each other and play an essential role in the enzymatic reaction. Both of the RNAs have a highly conserved helical region, P4, which is important in the catalytic reaction. We have used a bioinformatics approach based on conserved elements to computationally analyze available genomic sequences of eukaryotic organisms and have identified a large number of novel nuclear RNase P and MRP RNA genes. For MRP RNA for instance, this investigation increases the number of known sequences by a factor of three. We present secondary structure models of many of the predicted RNAs. Although all sequences are able to fold into the consensus secondary structure of P and MRP RNAs, a striking variation in size is observed, ranging from a Nosema locustae MRP RNA of 160 nt to much larger RNAs, e.g. a Plasmodium knowlesi P RNA of 696 nt. The P and MRP RNA genes appear in tandem in some protists, further emphasizing the close evolutionary relationship of these RNAs.

Anomalies de la transcription et diagnostic en génétique constitutionnelle

There is a rapidly growing literature on transcription abnormalities, e.g. differential expression of alleles and the role of some single nucleotide polymorphisms in altering splicing patterns. An average 10% of splicing mutations is reported in the Human Gene Mutation Database but this figure could climb to 50% for some genes such as NF1 or ATM. This paper therefore aims at clarifying some important aspects of transcriptional abnormalities in genetic testing. The main types of alterations are presented, i.e. exonic, intronic and promoter modifications that could modify or create consensus motif and/or secondary structures. DNA, RNA based-diagnostic strategies and in silico tools are then presented and their performances and limitations outlined to build up a picture of the current state of the art.

Large Subunit Mitochondrial rRNA Secondary Structures and Site-Specific Rate Variation in Two Lizard Lineages

A phylogenetic-comparative approach was used to assess and refine existing secondary structure models for a frequently studied region of the mitochondrial encoded large subunit (16S) rRNA in two large lizard lineages within the Scincomorpha, namely the Scincidae and the Lacertidae. Potential pairings and mutual information were analyzed to identify site interactions present within each lineage and provide consensus secondary structures. Many of the interactions proposed by previous models were supported, but several refinements were possible. The consensus structures allowed a detailed analysis of rRNA sequence evolution. Phylogenetic trees were inferred from Bayesian analyses of all sites, and the topologies used for maximum likelihood estimation of sequence evolution parameters. Assigning gamma-distributed relative rate categories to all interacting sites that were homologous between lineages revealed substantial differences between helices. In both lineages, sites within helix G2 were mostly conserved, while those within helix E18 evolved rapidly. Clear evidence of substantial site-specific rate variation (covarion-like evolution) was also detected, although this was not strongly associated with specific helices. This study, in conjunction with comparable findings on different, higher-level taxa, supports the ubiquitous nature of site-specific rate variation in this gene and justifies the incorporation of covarion models in phylogenetic inference.

Exploring consensus mRNA secondary (folding) structure units by stochastic sampling and folding simulation

It is extremely difficult to explore mRNA folding structure by biological experiments. In this report, we use stochastic sampling and folding simulation to test the existence of the stable secondary structural units of mRNA, look for the folding units, and explore the probabilistic stabilization of the units. Using this method, we made simulations for all possible local optimum secondary structures of a single strand mRNA within a certain range, and searched for the common parts of the secondary structures. The consensus secondary structure units (CSSUs) extracted from the above method are mainly hairpins, with a few single strands. These CSSUs suggest that the mRNA folding units could be relatively stable and could perform specific biological function. The significance of these observations for the mRNA folding problem in general is also discussed.

The influence of gapped positions in multiple sequence alignments on secondary structure prediction methods

All currently leading protein secondary structure prediction methods use a multiple protein sequence alignment to predict the secondary structure of the top sequence. In most of these methods, prior to prediction, alignment positions showing a gap in the top sequence are deleted, consequently leading to shrinking of the alignment and loss of position-specific information. In this paper we investigate the effect of this removal of information on secondary structure prediction accuracy. To this end, we have designed SymSSP, an algorithm that post-processes the predicted secondary structure of all sequences in a multiple sequence alignment by (i) making use of the alignment's evolutionary information and (ii) re-introducing most of the information that would otherwise be lost. The post-processed information is then given to a new dynamic programming routine that produces an optimally segmented consensus secondary structure for each of the multiple alignment sequences. We have tested our method on the state-of-the-art secondary structure prediction methods PHD, PROFsec, SSPro2 and JNET using the HOMSTRAD database of reference alignments. Our consensus-deriving dynamic programming strategy is consistently better at improving the segmentation quality of the predictions compared to the commonly used majority voting technique. In addition, we have applied several weighting schemes from the literature to our novel consensus-deriving dynamic programming routine. Finally, we have investigated the level of noise introduced by prediction errors into the consensus and show that predictions of edges of helices and strands are half the time wrong for all the four tested prediction methods.

Coenzyme B12 riboswitches are widespread genetic control elements in prokaryotes.

Recent studies have begun to reveal that numerous fundamental metabolic pathways in bacteria are regulated by riboswitches residing within certain messenger RNAs. These riboswitches selectively bind metabolites and modulate gene expression in response to changing ligand concentrations. Previously, we provided evidence that the btuB mRNAs of Escherichia coli and Salmonella typhimurium each carry a coenzyme B12-dependent riboswitch that causes repressed translation of the encoded cobalamin-transport protein at elevated coenzyme concentrations. Herein, we use a phylogenetic analysis to define a consensus sequence and secondary structure model for the ligand- binding domain of this riboswitch class. RNA structures that conform to this model are widespread in both Gram-positive and Gram-negative organisms. In addition, we find that the 5'-untranslated region (5'-UTR) of the cobalamin biosynthesis (cob) operon of S.typhimurium carries an RNA motif that matches this consensus sequence. Biochemical and genetic characterization of this motif confirms that the RNA directly binds coenzyme B12, and that it likely serves as a genetic control element for regulating expression of the 25-gene operon for cobalamin production in this pathogen.

Less isn't always more.

A favorite experiment of RNA scientists is the “cutdown” or “boundary” experiment where an RNA of interest is truncated from its 5′ or 3′ end and the function of the smaller molecules evaluated. Such experiments often reveal that only a small part of a large RNA molecule is needed for the function of interest. In addition, by reducing the number of alternative conformational states that can plague RNA preparations, truncated RNAs that emerge from boundary experiments are often more active than the original molecule. A common complication in the interpretation of boundary experiments arises from the general way in which RNAs fold into their active structures. As we all know, RNA secondary structure elements often form between two or more parts of the molecule that are quite distant on the primary sequence. A boundary experiment performed on an RNA containing a critical long-range secondary structure element will identify as essential all of the RNA between the 3′-most and 5′-most parts of the element, even if much of the internal sequence is not important for function. This complication can be resolved by comparing several active sequences, deducing a consensus secondary structure, and preparing an appropriate internal deletion. This classic combination of cutdowns and internal deletions has been used to define the minimal RNA binding sequences of many different proteins and small molecules, establish RNA regulatory regions and define catalytic RNA motifs.

Vienna RNA secondary structure server.

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/.

Structural and sequence evolution of U17 small nucleolar RNA (snoRNA) and its phylogenetic congruence in chelonians.

Vertebrate U17 RNA is an intron-encoded H/CA box containing snoRNA, which has been intensively studied in the last decade, though its precise role in ribosome biogenesis is not yet clear. A consensus secondary structure for the U17 RNA molecule has been derived from the comparative sequence and structural evolution analysis of U17 snoRNA among vertebrates. Its phylogenetic congruence above class level has been tested and preliminary data on chelonians suggest that also in this order, U17 snoRNA evolved congruently with phylogeny. We herein extend our analysis to other components of this reptile group. According to the sequence data that have also emerged from chelonians, the U17 RNA molecule can be divided into two main domains: the 5'-variable domain, which presents the sequence motifs capable of base-pairing with the 18S rRNA target and spanning STEM1, -2, and -3, and the 3'-conserved domain, consisting of STEM4. In vertebrates, the latter RNA region shows a high conservation both in secondary structure and in the presence of the three sequence motifs 5'-AUUCCUA-3', 5'-U(G/U)ACU-3', and 5'-AACCC-3'. We tested the phylogenetic congruence of U17 evolution with chelonian relationships: Our results are significantly similar to those emerging from mtDNA and morphological systematics. Some discrepancies (e.g., the position of Platysternon) need to be addressed in greater depth.

Characterization of three novel snoRNAs from Schizosaccharomyces pombe implies the structural and functional diversity of box H/ACA snoRNA family

Site-specific pseudouridylations of rRNAs and snRNAs are directed by small guide RNAs that usually possess conserved H and ACA motifs and the consensus secondary structure consisting of hairpin-hinge-hairpin-tail. In this study, through screening cDNA library, three novel box H/ACA snoRNAs, termed SP12, SP16 and SP20, have been identified from fisson yeast Schizosacchromyces pombe. Each of the small RNAs possesses an ACA-like motif and exhibits secondary structures containing two or three large hairpinlike domains. Although the box H is noncanonical, SP20 can be potentially folded into a typical box H/ACA snoRNA-like structure and is predicted to guide two pseudouridylations, U1208 and U1053, on 18S rRNA. In contrast, SP12 and SP16 display different structures from SP20 and have no antisense elements to rRNAs or snRNAs. Obviously, they are new members of orphan guide snoRNAs whose functions remain elusive. The genes encoded for SP12 and SP16 are both located in long intergenic regions between protein-coding genes. A post-transcriptional processing for SP16 precursor was demonstrated. Interestingly, SP20 gene is nested inside a very hypothetical ORF, but transcribed in a reverse sense. Northern and RT-PCR analyses could not detect the transcript from the hypothetical ORF, and therefore exclude the expression of the spurious gene. This is the first example that a genetic locus predicted for protein-coding gene actually encodes a small non-messenger RNA in S. pombe.

RNAML: a standard syntax for exchanging RNA information.

Analyzing a single data set using multiple RNA informatics programs often requires a file format conversion between each pair of programs, significantly hampering productivity. To facilitate the interoperation of these programs, we propose a syntax to exchange basic RNA molecular information. This RNAML syntax allows for the storage and the exchange of information about RNA sequence and secondary and tertiary structures. The syntax permits the description of higher level information about the data including, but not restricted to, base pairs, base triples, and pseudoknots. A class-oriented approach allows us to represent data common to a given set of RNA molecules, such as a sequence alignment and a consensus secondary structure. Documentation about experiments and computations, as well as references to journals and external databases, are included in the syntax. The chief challenge in creating such a syntax was to determine the appropriate scope of usage and to ensure extensibility as new needs will arise. The syntax complies with the eXtensible Markup Language (XML) recommendations, a widely accepted standard for syntax specifications. In addition to the various generic packages that exist to read and interpret XML formats, an XML processor was developed and put in the open-source MC-Core library for nucleic acid and protein structure computer manipulation.

Homology between O-linked GlcNAc transferases and proteins of the glycogen phosphorylase superfamily

The O-linked GlcNAc transferases (OGTs) are a recently characterized group of largely eukaryotic enzymes that add a single β-N-acetylglucosamine moiety to specific serine or threonine hydroxyls. In humans, this process may be part of a sugar regulation mechanism or cellular signaling pathway that is involved in many important diseases, such as diabetes, cancer, and neurodegeneration. However, no structural information about the human OGT exists, except for the identification of tetratricopeptide repeats (TPR) at the N terminus. The locations of substrate binding sites are unknown and the structural basis for this enzyme’s function is not clear. Here, remote homology is reported between the OGTs and a large group of diverse sugar processing enzymes, including proteins with known structure such as glycogen phosphorylase, UDP-GlcNAc 2-epimerase, and the glycosyl transferase MurG. This relationship, in conjunction with amino acid similarity spanning the entire length of the sequence, implies that the fold of the human OGT consists of two Rossmann-like domains C-terminal to the TPR region. A conserved motif in the second Rossmann domain points to the UDP-GlcNAc donor binding site. This conclusion is supported by a combination of statistically significant PSI-BLAST hits, consensus secondary structure predictions, and a fold recognition hit to MurG. Additionally, iterative PSI-BLAST database searches reveal that proteins homologous to the OGTs form a large and diverse superfamily that is termed GPGTF (glycogen phosphorylase/glycosyl transferase). Up to one-third of the 51 functional families in the CAZY database, a glycosyl transferase classification scheme based on catalytic residue and sequence homology considerations, can be unified through this common predicted fold. GPGTF homologs constitute a substantial fraction of known proteins: 0.4 % of all non-redundant sequences and about 1 % of proteins in the Escherichia coli genome are found to belong to the GPGTF superfamily.

The kink-turn: a new RNA secondary structure motif.

Analysis of the Haloarcula marismortui large ribosomal subunit has revealed a common RNA structure that we call the kink-turn, or K-turn. The six K-turns in H.marismortui 23S rRNA superimpose with an r.m.s.d. of 1.7 A. There are two K-turns in the structure of Thermus thermophilus 16S rRNA, and the structures of U4 snRNA and L30e mRNA fragments form K-turns. The structure has a kink in the phosphodiester backbone that causes a sharp turn in the RNA helix. Its asymmetric internal loop is flanked by C-G base pairs on one side and sheared G-A base pairs on the other, with an A-minor interaction between these two helical stems. A derived consensus secondary structure for the K-turn includes 10 consensus nucleotides out of 15, and predicts its presence in the 5'-UTR of L10 mRNA, helix 78 in Escherichia coli 23S rRNA and human RNase MRP. Five K-turns in 23S rRNA interact with nine proteins. While the observed K-turns interact with proteins of unrelated structures in different ways, they interact with L7Ae and two homologous proteins in the same way.