Higher order structure in the 3'-minor domain of small subunit ribosomal RNAs from a gram negative bacterium, a gram positive bacterium and a eukaryote

An experimental approach was used to determine and compare the highest order structure within the 150 to 200 nucleotides at the 3'-ends of the RNAs from the small ribosomal subunits of Escherichia coli, Bacillus stearothermophilus and Saccharomyces cerevisiae. Chemical reagents were employed to establish the degree of stacking and/or accessibility of each adenosine, guanosine and cytidine. The double helices were probed with a cobra venom ribonuclease from Naja naja oxiana, and the relatively unstructured and accessible sequences were localized with the single strand-specific ribonucleases A, T1, T2 and S1. The data enabled the various minimal secondary structural models, proposed for the 3'-regions of the E. coli and S. cerevisiae RNAs, to be critically examined, and to demonstrate that the main common features of these models are correct. The results also reveal the presence and position of additional higher order structure in the renatured free RNA. It can be concluded that a high level of conservation of higher order structure has occurred during the evolution of the gram negative and gram positive eubacteria and the eukaryote in both the double helical regions and the "unstructured" regions. Several unusual structural features were detected. Multiple G X A pairings in two of the putative helices, which are compatible with phylogenetic sequence comparisons, are strongly supported by the occurrence of cobra venom ribonuclease cuts adjacent to, and in one case between, these pairings. Evidence is also provided for the stacking of an A X A pair within a double helix of the yeast RNA. Other special structural features include adenosines bulged out from double helices; such nucleotides, which are hyper-reactive, have been implicated in protein recognition in 5 S ribosomal RNA. The 3'-terminal regions of the RNAs are particularly important for the functioning of the ribosome. They are involved in mRNA, tRNA and ribosomal factor binding. The results reveal that while the functionally important RNA sequences tend to be conserved, they are not always accessible in the free RNA; the pyrimidine-rich "Shine and Dalgarno" sequence, for example, which is involved in mRNA recognition, occurs in a double helix in both eubacterial RNAs.     

RNA secondary structure and compensatory evolution

The classic concept of epistatic fitness interactions between genes has been extended to study interactions within gene regions, especially between nucleotides that are important in maintaining pre-mRNA/mRNA secondary structures. It is shown that the majority of linkage disequilibria found within the Drosophila Adh gene are likely to be caused by epistatic selection operating on RNA secondary structures. A recently proposed method of RNA secondary structure prediction based on DNA sequence comparisons is reviewed and applied to several types of RNAs, including tRNA, rRNA, and mRNA. The patterns of covariation in these RNAs are analyzed based on Kimura's compensatory evolution model. The results suggest that this model describes the substitution process in the pairing regions (helices) of RNA secondary structures well when the helices are evolutionarily conserved and thermodynamically stable, but fails in some other cases. Epistatic selection maintaining pre-mRNA/mRNA secondary structures is compared to weak selective forces that determine features such as base composition and synonymous codon usage. The relationships among these forces and their relative strengths are addressed. Finally, our mutagenesis experiments using the Drosophila Adh locus are reviewed. These experiments analyze long-range compensatory interactions between the 5' and 3' ends of Adh mRNA, the different constraints on secondary structures in introns and exons, and the possible role of secondary structures in RNA splicing.     

Comparative sequence analysis and patterns of covariation in RNA secondary structures

A novel method of RNA secondary structure prediction based on a comparison of nucleotide sequences is described. This method correctly predicts nearly all evolutionarily conserved secondary structures of five different RNAs: tRNA, 5S rRNA, bacterial ribonuclease P (RNase P) RNA, eukaryotic small subunit rRNA, and the 3' untranslated region (UTR) of the Drosophila bicoid (bcd) mRNA. Furthermore, covariations occurring in the helices of these conserved RNA structures are analyzed. Two physical parameters are found to be important determinants of the evolution of compensatory mutations: the length of a helix and the distance between base-pairing nucleotides. For the helices of bcd 3' UTR mRNA and RNase P RNA, a positive correlation between the rate of compensatory evolution and helix length is found. The analysis of Drosophila bcd 3' UTR mRNA further revealed that the rate of compensatory evolution decreases with the physical distance between base-pairing residues. This result is in qualitative agreement with Kimura's model of compensatory fitness interactions, which assumes that mutations occurring in RNA helices are individually deleterious but become neutral in appropriate combinations.     

Comparative photocross-linking analysis of the tertiary structures of Escherichia coli and Bacillus subtilis RNase P RNAs.

Bacterial ribonuclease P contains a catalytic RNA subunit that cleaves precursor sequences from the 5' ends of pre-tRNAs. The RNase P RNAs from Bacillus subtilis and Escherichia coli each contain several unique secondary structural elements not present in the other. To understand better how these phylogenetically variable elements affect the global architecture of the ribozyme, photoaffinity cross-linking studies were carried out. Photolysis of photoagents attached at homologous sites in the two RNAs results in nearly identical cross-linking patterns, consistent with the homology of the RNAs and indicating that these RNAs contain a common, core tertiary structure. Distance constraints were used to derive tertiary structure models using a molecular mechanics-based modeling protocol. The resulting superimposition of large sets of equivalent models provides a low resolution (5-10 A) structure for each RNA. Comparison of these structure models shows that the conserved core helices occupy similar positions in space. Variably present helical elements that may play a role in global structural stability are found at the periphery of the core structure. The P5.1 and P15.1 helical elements, unique to the B.subtilis RNase P RNA, and the P6/16/17 helices, unique to the E.coli RNA, occupy similar positions in the structure models and, therefore, may have analogous structural function.     

Co-conservation of rRNA tetraloop sequences and helix length suggests involvement of the tetraloops in higher-order interactions

Terminal loops containing four nucleotides (tetraloops) are common in structural RNAs, and they frequently conform to one of three sequence motifs, GNRA, UNCG, or CUUG. Here we compare available sequences and secondary structures for rRNAs from bacteria, and we show that helices capped by phylogenetically conserved GNRA loops display a strong tendency to be of conserved length. The simplest interpretation of this correlation is that the conserved GNRA loops are involved in higher-order interactions, intramolecular or intermolecular, resulting in a selective pressure for maintaining the lengths of these helices. A small number of conserved UNCG loops were also found to be associated with conserved length helices, consistent with the possibility that this type of tetraloop also takes part in higher-order interactions.     

Determination of secondary structure in rabbit globin messenger RNA by thermal denaturation.

The secondary structure of highly purified globin messenger RNA has been investigated by alkaline hydrolysis, nuclease digestion, and thermal denaturation. The thermal denaturation properties of globin messenger have been compared to poly(U), poly (A), and a synthetic random sequence RNA copolymer. From these studies it is concluded that globin mRNA contains considerable secondary structure and that the amount of helical structure is greater than that which occurs with a random sequence polyribonucleotide. Globin mRNA contains, by comparison to the secondary structures of native DNA, tRNAs, or 18S rRNA, helices with involve 55-62% of the bases or 58-68% if a correction is made for the 3'-terminal poly(A) segment. The helices of globin mRNA appear to be unique as differences in the NaCl stabilization of this RNA have been noted when compared to other naturally ooccurring and synthetic RNAs. Comparison of the hyperchromicity maxima, obtained at 260 and 280 nm for globin mRNA and 18S rRNA, indicates that the helices of the two RNAs contain similar numbers of G-C base pairs. Differential analysis of NaCl stabilization curves indicate three discrete thermally denaturable helix types in globin mRNA.     

Repeated parallel evolution of minimal rRNAs revealed from detailed comparative analysis

The concept of a minimal ribosomal RNA-containing ribosome, a structure with a minimal set of elements capable of providing protein biosynthesis, is essential for understanding this fundamental cellular process. Nematodes and trypanosomes have minimal mitochondrial rRNAs and detailed reconstructions of their secondary structures indicate that certain conserved helices have been lost in these taxa. In contrast, several recent studies on acariform mites have argued that minimal rRNAs may evolve via shortening of secondary structure elements but not the loss of these elements as shown for trypanosomes and nematodes. Based on extensive structural analysis of chelicerate arthropods, we demonstrate that extremely short rRNAs of acariform mites share certain structural modifications with nematodes and trypanosomes: loss of helices of the GTPase region and divergence in the evolutionarily conserved connecting loop between helices H1648 and H1764 of the large subunit rRNA. These highly concerted parallel modifications indicate that minimal rRNAs were generated under the strong selection that favored or tolerated reductions of helices in particular locations while maintaining the functionality of the rRNA molecules throughout evolution. We also discuss potential evolution of minimal rRNAs and atypical transfer RNAs.     

A computer method for finding common base paired helices in aligned sequences: application to the analysis of random sequences.

We describe a new computer program that identifies conserved secondary structures in aligned nucleotide sequences of related single-stranded RNAs. The program employs a series of hash tables to identify and sort common base paired helices that are located in identical positions in more than one sequence. The program gives information on the total number of base paired helices that are conserved between related sequences and provides detailed information about common helices that have a minimum of one or more compensating base changes. The program is useful in the analysis of large biological sequences. We have used it to examine the number and type of complementary segments (potential base paired helices) that can be found in common among related random sequences similar in base composition to 16S rRNA from Escherichia coli. Two types of random sequences were analyzed. One set consisted of sequences that were independent but they had the same mononucleotide composition as the 16S rRNA. The second set contained sequences that were 80% similar to one another. Different results were obtained in the analysis of these two types of random sequences. When 5 sequences that were 80% similar to one another were analyzed, significant numbers of potential helices with two or more independent base changes were observed. When 5 independent sequences were analyzed, no potential helices were found in common. The results of the analyses with random sequences were compared with the number and type of helices found in the phylogenetic model of the secondary structure of 16S ribosomal RNA. Many more helices are conserved among the ribosomal sequences than are found in common among similar random sequences. In addition, conserved helices in the 16S rRNAs are, on the average, longer than the complementary segments that are found in comparable random sequences. The significance of these results and their application in the analysis of long non-ribosomal nucleotide sequences is discussed.     

Conserved features of Y RNAs revealed by automated phylogenetic secondary structure analysis.

Y RNAs are small 'cytoplasmic' RNAs which are components of the Ro ribonucleoprotein (RNP) complex. The core of this complex, which is found in the cell nuclei of higher eukaryotes as well as the cytoplasm, is composed of a complex between the 60 kDa Ro protein and Y RNAs. Human cells contain four distinct Y RNAs (Y1, Y3, Y4 and Y5), while other eukaryotes contain a variable number of Y RNA homologues. When detected in a particular species, the Ro RNP has been present in every cell type within that particular organism. This characteristic, along with its high conservation among vertebrates, suggests an important function for Ro RNP in cellular metabolism; however, this function has not yet been definitively elucidated. In order to identify conserved features of Y RNA sequences and structures which may be directly involved in Ro RNP function, a phylogenetic comparative analysis of Y RNAs has been performed. Sequences of Y RNA homologues from five vertebrate species have been obtained and, together with previously published Y RNA sequences, used to predict Y RNA secondary structures. A novel RNA secondary structure comparison algorithm, the suboptimal RNA analysis program, has been developed and used in conjunction with available algorithms to find phylogenetically conserved secondary structure models for YI, Y3 and Y4 RNAs. Short, conserved sequences within the Y RNAs have been identified and are invariant among vertebrates, consistent with a direct role for Y RNAs in Ro function. A subset of these are located wholly or partially in looped regions in the Y3 and Y4 RNA predicted model structures, in accord with the possibility that these Y RNAs base pair with other cellular nucleic acids or are sites of interaction between the Ro RNP and other macromolecules.     

Evolutionary relationships amongst archaebacteria. A comparative study of 23 S ribosomal RNAs of a sulphur-dependent extreme thermophile, an extreme halophile and a thermophilic methanogen

The 23 S RNA genes representative of each of the main archaebacterial subkingdoms, Desulfurococcus mobilis an extreme thermophile, Halococcus morrhuae an extreme halophile and Methanobacterium thermoautotrophicum a thermophilic methanogen, were cloned and sequenced. The inferred RNA sequences were aligned with all the available 23 S-like RNAs of other archaebacteria, eubacteria/chloroplasts and the cytoplasm of eukaryotes. Universal secondary structural models containing six major structural domains were refined, and extended, using the sequence comparison approach. Much of the present structure was confirmed but six new helices were added, including one that also exists in the eukaryotic 5.8 S RNA, and extensions were made to several existing helices. The data throw doubt on whether the 5' and 3' ends of the 23 S RNA interact, since no stable helix can form in either the extreme thermophile or the methanogen RNA. A few secondary structural features, specific to the archaebacterial RNAs were identified; two of these were supported by a comparison of the archaebacterial RNA sequences, and experimentally, using chemical and ribonuclease probes. Seven tertiary structural interactions, common to all 23 S-like RNAs, were predicted within unpaired regions of the secondary structural model on the basis of co-variation of nucleotide pairs; two lie in the region of the 23 S RNA corresponding to 5.8 S RNA but they are not conserved in the latter. The flanking sequences of each of the RNAs could base-pair to form long RNA processing stems. They were not conserved in sequence but each exhibited a secondary structural feature that is common to all the archaebacterial stems for both 16 S and 23 S RNAs and constitutes a processing site. Kingdom-specific nucleotides have been identified that are associated with antibiotic binding sites at functional centres in 23 S-like RNAs: in the peptidyl transferase centre (erythromycin-domain V) the archaebacterial RNAs classify with the eukaryotic RNAs; at the elongation factor-dependent GTPase centre (thiostrepton-domain II) they fall with the eubacteria, and at the putative amino acyl tRNA site (alpha-sarcin-domain VI) they resemble eukaryotes. Two of the proposed tertiary interactions offer a structural explanation for how functional coupling of domains II and V occurs at the peptidyl transferase centre. Phylogenetic trees were constructed for the archaebacterial kingdom, and for the other two kingdoms, on the basis of the aligned 23 S-like RNA sequences.(ABSTRACT TRUNCATED AT 400 WORDS)     

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.     

Sequence comparison and secondary structure analysis of the 3' noncoding region of flavivirus genomes reveals multiple pseudoknots.

Sequences of 191 flavivirus RNAs belonging to four sero-groups were used to predict the secondary structure of the 3' noncoding region (3' NCR) directly upstream of the conserved terminal hairpin. In mosquito-borne flavivirus RNAs (n = 164) a characteristic structure element was identified that includes a phylogenetically well-supported pseudoknot. This element is repeated in the dengue and Japanese encephalitis RNAs and centers around the conserved sequences CS2 and RCS2. In yellow fever virus RNAs that contain one CS2 motif, only one copy of this pseudoknotted structure was found. The conserved pseudoknotted element is absent from the 3' NCR of tick-borne virus RNAs, which altogether adopt a secondary structure that is very different from that of mosquito-borne virus RNAs. The strong conservation of the pseudoknot in mosquito-borne flavivirus RNAs implies a stronger relationship between these viruses than concluded from previous secondary structure analyses. The role of the (tandem) pseudoknots in flavivirus replication is discussed.     

Higher-order structure in the 3′-terminal domain VI of the 23 S ribosomal RNAs from Escherichia coli and Bacillus stearothermophilus

An experimental approach was used to determine, and compare, the higher-order structure within domain VI of the 23 S ribosomal RNAs from Escherichia coli and Bacillus stearothermophilus. This domain, which encompasses approximately 300 nucleotides at the 3′ end of the RNAs, consists of two large subdomains. The 5′ subdomain has been conserved during evolution and appears to be functionally important for the binding of the EF-1 · GTP · aminoacyl-tRNA complex in eukaryotes. The 3′ subdomain has diverged widely between eubacteria and eukaryotes and has produced the 4.5 S RNA in the chloroplast ribosomes of flowering plants.The structure of domain VI within the eubacterial RNAs was probed with chemical reagents in order to establish the degree of stacking and/or accessibility of each adenosine, cytidine and guanosine residue; the double-helical segments were localized with the cobra venom ribonuclease from Naja naja oxiana, and the relatively unstructured and accessible sequences were detected with the single-strand-specific ribonucleases A, T₁ and T₂. The data enabled the three secondary structural models, proposed for the E. coli 23 S RNAs, to be examined critically and it was concluded that many of their structural features are correct. Various differences between the models were considered and evidence is provided for additional structuring in the RNA including the stacking of juxtaposed purines into double helices. The 5′ subdomain constitutes a compact and resistant structure whereas the 3′ subdomain is relatively accessible and contains most of the potential protein binding sites. Moreover, comparison of our results with the published results on 4.5 S RNA suggests that the latter forms essentially the same structure as the 3′ subdomain, in contrast to earlier conclusions.A high level of structural conservation has occurred throughout the RNA domain during the evolution of the Gram negative and Gram positive bacteria although the thermophile was generally more stable at base-pairs adjacent to the terminal loops.     

Intracellular localization and unique conserved sequences of three small nucleolar RNAs.

Three human small nucleolar RNAs (snoRNAs), E1, E2 and E3, were reported earlier that have unique sequences, interact directly with unique segments of pre-rRNA in vivo and are encoded in introns of protein genes. In the present report, human and frog E1, E2 and E3 RNAs injected into the cytoplasm of frog oocytes migrated to the nucleus and specifically to the nucleolus. This indicates that the nucleolar and nuclear localization signals of these snoRNAs reside within their evolutionarily conserved segments. Homologs of these snoRNAs from several vertebrates were sequenced and this information was used to develop RNA secondary structure models. These snoRNAs have unique phylogenetically conserved sequences.     

Sequences of three molluscan 5 S ribosomal RNAs confirm the validity of a dynamic secondary structure model.

The collection of known 5 S rRNA primary structures is enriched with the sequences from three mollusca, the snails Helix pomatia and Arion rufus, and the mussel Mytilus edulis. The three sequences can be fitted in a five-helix secondary structure model previously shown (De Wachter et al. (1982) Biochimie 64, 311-329) to apply to all 5 S RNAs regardless of their origin. One of the helices in this model can undergo a bulge-internal loop transition. Within the metazoan kingdom, the dimensions of each helix and loop are rigidly conserved, except for one helix which can comprise either 6 or 7 base pairs.     

Phylogenetic analysis and secondary structure of the Bacillus subtilis bacteriophage RNA required for DNA packaging

An unusual RNA molecule encoded by the Bacillus subtilis bacteriophage phi 29 is a structural component of the viral prohead and is required for the ATP-dependent packaging of DNA. Here we report a model of secondary structure for this prohead RNA developed from a phylogenetic analysis of the primary sequences of prohead RNAs of related phages. Twenty-nine phages related to phi 29 were found to produce prohead RNAs. These RNAs were analyzed by their ability to replace phi 29 RNA in in vitro phage assembly, by Northern blot hybridization with a probe complementary to phi 29 RNA, and by partial and complete sequence analyses. These analyses revealed four quite different sequences ranging in length from 161 to 174 residues. The secondary structure deduced from these sequences, in agreement with earlier observations, indicated that prohead RNA is organized into two domains. The larger 5'-domain (Domain I) is composed of 113-117 residues and contains four helices. Three of these helices appear to be organized into a central stem that is interrupted by two unpaired loops and the fourth helix and loop. The smaller 3'-domain (Domain II) is composed of 40-44 residues and consists of two helices. Domains I and II are separated by 8-13 unpaired residues. Nuclease cleavage occurs readily in this single-stranded joining region, and this cleavage allows the subsequent separation of the two RNA domains. The separated Domain I is fully active in DNA packaging in vitro. The functional significance and biological role of Domain II are unknown. The phylogenetic secondary structure model provides a basis for further analysis of the role of this RNA in bacteriophage morphogenesis.     

Equilibria in 5-S ribosomal RNA secondary structure. Bulges and interior loops in 5-S RNA secondary structure may serve as articulations for a flexible molecule

The basic assumption in this paper is that the secondary structure of a 5-S ribosomal RNA cannot be represented by a single model. We propose that the molecule can adopt, at least within the ribosome, a series of slightly different structures of nearly equal stability. The different structures arise from the existence of ambiguous base-pairing opportunities in bulged helices and the adjacent interior loops. In eubacterial 5-S RNAs there is one such an area, in eukaryotic 5-S RNAs two such areas that can give rise to structural switches. We explain how a change in secondary structure in these areas may influence the relative orientation of the surrounding helices, in other words how bulges and interior loops may serve as articulations and give rise to a flexible tertiary structure.     

Secondary structure of ovalbumin messenger RNA.

The secondary structure of highly purified ovalbumin mRNA was studied by automated thermal denaturation techniques and the data were subjected to computer processing. Comparative studies with 20 natural and synthetic model nucleic acids suggested that the secondary structure of ovalbumin mRNA possesses the following features: the extent of base pairing of ovalbumin mRNA is similar to that found in tRNAs or ribosomal RNAs; the secondary structure of ovalbumin mRNA is more thermolabile than any of the model compounds tested, including the copolymer poly(A-U); ovalbumin mRNA does not have extensive G-C rich stems as found in tRNAs or ribosomal RNAs; the base composition of the double-stranded regions reveals 54% G-C residues which was significantly higher than that noted in the whole molecule (approximately 41.5% G-C). The presence of 46% A-U pairs in short stems of about five base pairs would have a very large destabilizing effect on the secondary structure of ovalbumin mRNA. However, at 0.175 M monovalent cations and 36 degrees C most of the secondary structure of ovalbumin mRNA is preserved. These data suggest that the double-stranded regions in ovalbumin mRNA are of sufficient length to provide the necessary stability for maintaining the open loop regions in an appropriate conformation which may be required for the biological function of ovalbumin mRNA. Furthermore, the lability of the double-stranded regions in ovalbumin mRNA may also be important for the biological function of this mRNA.     

Phylogenetic analysis of the structure of RNase MRP RNA in yeasts

RNase MRP is a ribonucleoprotein enzyme involved in processing precursor rRNA in eukaryotes. To facilitate our structure-function analysis of RNase MRP from Saccharomyces cerevisiae, we have determined the likely secondary structure of the RNA component by a phylogenetic approach in which we sequenced all or part of the RNase MRP RNAs from 17 additional species of the Saccharomycetaceae family. The structure deduced from these sequences contains the helices previously suggested to be common to the RNA subunit of RNase MRP and the related RNA subunit of RNase P, an enzyme cleaving tRNA precursors. However, outside this common region, the structure of RNase MRP RNA determined here differs from a previously proposed universal structure for RNase MRPs. Chemical and enzymatic structure probing analyses were consistent with our revised secondary structure. Comparison of all known RNase MRP RNA sequences revealed three regions with highly conserved nucleotides. Two of these regions are part of a helix implicated in RNA catalysis in RNase P, suggesting that RNase MRP may cleave rRNA using a similar catalytic mechanism.