Secondary structure of Tetrahymena thermophilia 5S ribosomal RNA as revealed by enzymatic digestion and microdensitometric analysis.

The secondary structure of [32P] end-labeled 5S rRNA from Tetrahymena thermophilia (strain B) has been investigated using the enzymes S1 nuclease, cobra venom ribonuclease and T2 ribonuclease. The results, analyzed by scanning microdensitometry and illustrated by three-dimensional computer graphics, support the secondary structure model of Curtiss and Vournakis for 5S rRNA. Aberrent mobility of certain RNA fragments on sequencing gels was observed as regions of band compression. These regions are postulated to be caused by stable internal base-pairing. The molecule was probed with T2 RNase in neutral (pH 7.5) and acidic (pH 4.5) buffers and only minor structural differences were revealed. One of the helices was found to be susceptible to enzymatic attack by both the single-strand and double-strand specific enzymes. These observations are evidence for the existence of dynamic structural equilibria in 5S rRNA.     

Nucleotide sequence of 5S ribosomal RNA from rainbow trout (Salmo gairdnerii) liver.

Staring from low molecular weight RNA obtained from rainbow trout (Salmo gairdnerii) liver, 5S ribosomal RNA (rRNA) was highly purified by successive chromatography on columns of DEAE-Sephadex A50 and Sephadex G100. Products of complete and partial digestions on this RNA with pancreatic ribonuclease (RNase A) [EC 3.1.4.22] and RNase T [EC 3.1.4.8] were isolated and sequenced by conventional and high-performance liquid chromatography (HPLC) procedures. The nucleotide sequence of this RNA thus established was compared with those of five other vertebrae 5S rRNAs, and the rates of base substitution per site per year were found to be nearly constant in these RNAs. The analyses of the partial digests of the trout 5S rRNA revealed several sites susceptible to RNase attack, which could be accounted for by the secondary structure model for eukaryotic 5S rRNAs proposed by Nishikawa and Takemura (1).     

Structural analysis of prokaryotic and eukaryotic 5S rRNAs by RNase H

The availabilities of single-stranded 5S rRNA regions c, d and d' for base pairing interactions were analyzed by using synthetic DNA oligomers. Hybrid formation was detected by the endonucleolytical mode of the RNA-DNA specific action of RNase H. Provided that the hybrid interaction involved 6 successive base pairs, 5S rRNA loop c nucleotides 42-47 displayed accessibility in Escherichia coli, Bacillus stearothermophilus and Thermus thermophilus 5S rRNAs as well as in eukaryotic 5S rRNAs from Saccharomyces carlsbergensis, Rattus rattus and Equisetum arvense. Investigating eubacterial 5S rRNA regions d and d' (nucleotides 71-76 and 99-105, respectively), susceptibility was observed in E. coli 5S rRNA which, however, decreases in B. stearothermophilus and even more so in T. thermophilus 5S rRNA. For additional evaluation of the data obtained by RNase H cleavage, association constants of the hexanucleotides were determined by equilibrium dialysis at 4 degrees C for B. stearothermophilus 5S rRNA. The results obtained reveal that nucleotides 36-41 of B. stearothermophilus 5S rRNA are inaccessible for Watson-Crick interaction, which suggests that this part of loop c is in a structurally constrained configuration, or buried in the tertiary structure or involved in tertiary interactions.     

Recognition of G-U mismatches by tris(4,7-diphenyl-1,10-phenanthroline)rhodium(III).

The coordination complex tris(4,7-diphenyl-1,10-phenanthroline)rhodium(III) [Rh(DIP)3(3+)], which promotes RNA cleavage upon photoactivation, has been shown to target specifically guanine-uracil (G-U) mismatches in double-helical regions of folded RNAs. Photoactivated cleavage by Rh(DIP)3(3+) has been examined on a series of RNAs that contain G-U mismatches, yeast tRNA(Phe) and yeast tRNA(Asp), as well as on 5S rRNAs from Xenopus oocytes and Escherichia coli. In addition, a "microhelix" was synthesized, which consists of seven base pairs of the acceptor stem of yeast tRNA(Phe) connected by a six-nucleotide loop and contains a mismatch involving residues G4 and U69. A U4.G69 variant of this sequence was also constructed, and cleavage by Rh(DIP)3(3+) was examined. In each of these cases, specific cleavage is observed at the residue which lies to the 3'-side of the wobble-paired U; some cleavage by the rhodium complex is also evident in several structured RNA loops. The remarkable site selectivity for G-U mismatches within double-helical regions is attributed to shape-selective binding by the rhodium complex. This binding furthermore depends upon the orientation of the G-U mismatch, which produces different stacking interactions between the G-U base pair with the Watson-Crick base pair following it on the 5'-side of U compared to the Watson-Crick pair preceding it on the 3'-side of U. Rh(DIP)3(3+) therefore serves as a unique probe of G-U mismatches and may be useful both as a model and in probing RNA-protein interactions as well as in identifying G-U mismatches within double-helical regions of folded RNAs.     

Phylogenetic analysis and evolution of RNase P RNA in proteobacteria.

The secondary structures of the eubacterial RNase P RNAs are being elucidated by a phylogenetic comparative approach. Sequences of genes encoding RNase P RNA from each of the recognized subgroups (alpha, beta, gamma, and delta) of the proteobacteria have now been determined. These sequences allow the refinement, to nearly the base pair level, of the phylogenetic model for RNase P RNA secondary structure. Evolutionary change among the RNase P RNAs was found to occur primarily in four discrete structural domains that are peripheral to a highly conserved core structure. The new sequences were used to examine critically the proposed similarity (C. Guerrier-Takada, N. Lumelsky, and S. Altman, Science 246:1578-1584, 1989) between a portion of RNase P RNA and the "exit site" of the 23S rRNA of Escherichia coli. Phylogenetic comparisons indicate that these sequences are not homologous and that any similarity in the structures is, at best, tenuous.     

Escherichia coli 5S RNA A and B conformers. Characterisation by enzymatic and chemical methods

The structures of the two stable conformers of Escherichia coli 5 S RNA, the and B form, were compared. Information about the structures were obtained using the methods of limited enzymatic hydrolysis and chemical modification of accessible nucleotides. Base-specific modifications were performed for adenosines and cytidines using diethylpyrocarbonate and dimethylsulfate in combination with a strand-scission reaction at the modified site. Base-specific (RNase T1) as well as conformation-specific (nuclease S1, cobra venom nuclease) enzymes were employed for the limited enzymatic hydrolysis. Clear differences in the accessibility of the two 5 S RNA conformers to the enzymes and the chemical reagents were established and the regions with altered reactivities were localized in the 5 S RNA structure. The results are consistent with the disruption of the secondary structural interactions in helix II and partly in helices III and IV during the transition from the A to the B form. (The numbering of the helices is according to the generally accepted Fox and Woese model.) In addition some regions presumably involved in the tertiary structure are distorted. There is evidence, however, for the new formation of structural regions between two distant sites in the 5 S RNA B form. The results enable us to refine the existing 5 S RNA A-form model and provide insight into the structural dynamics that lead to the formation of the 5 S RNA B form.     

Probing the structure of mouse Ehrlich ascites cell 5.8S, 18S and 28S ribosomal RNA in situ.

The secondary structure of mouse Ehrlich ascites 18S, 5.8S and 28S ribosomal RNA in situ was investigated by chemical modification using dimethyl sulphate and 1-cyclohexyl-3-(morpholinoethyl) carbodiimide metho-p-toluene sulphonate. These reagents specifically modify unpaired bases in the RNA. The reactive bases were localized by primer extension followed by gel electrophoresis. The three rRNA species were equally accessible for modification i.e. approximately 10% of the nucleotides were reactive. The experimental data support the theoretical secondary structure models proposed for 18S and 5.8/28S rRNA as almost all modified bases were located in putative single-strand regions of the rRNAs or in helical regions that could be expected to undergo dynamic breathing. However, deviations from the suggested models were found in both 18S and 28S rRNA. In 18S rRNA some putative helices in the 5'-domain were extensively modified by the single-strand specific reagents as was one of the suggested helices in domain III of 28S rRNA. Of the four eukaryote specific expansion segments present in mouse Ehrlich ascites cell 28S rRNA, segments I and III were only partly available for modification while segments II and IV showed average to high modification.     

Improvement of RNA secondary structure prediction using RNase H cleavage and randomized oligonucleotides

RNA secondary structure prediction using free energy minimization is one method to gain an approximation of structure. Constraints generated by enzymatic mapping or chemical modification can improve the accuracy of secondary structure prediction. We report a facile method that identifies single-stranded regions in RNA using short, randomized DNA oligonucleotides and RNase H cleavage. These regions are then used as constraints in secondary structure prediction. This method was used to improve the secondary structure prediction of Escherichia coli 5S rRNA. The lowest free energy structure without constraints has only 27% of the base pairs present in the phylogenetic structure. The addition of constraints from RNase H cleavage improves the prediction to 100% of base pairs. The same method was used to generate secondary structure constraints for yeast tRNA^Phe, which is accurately predicted in the absence of constraints (95%). Although RNase H mapping does not improve secondary structure prediction, it does eliminate all other suboptimal structures predicted within 10% of the lowest free energy structure. The method is advantageous over other single-stranded nucleases since RNase H is functional in physiological conditions. Moreover, it can be used for any RNA to identify accessible binding sites for oligonucleotides or small molecules.     

Resistance of both the 5'- and 3'-domain of isolated Escherichia coli 23S rRNA against digestion with alpha-sarcin

The ribonuclease alpha-sarcin exclusively cleaves the phosphodiester bond after G2661 in the 23S rRNA within 50S subunits, thus inactivating the ribosomes. The resulting alpha-fragment is 243 nucleotides long and contains the 3'-end of the 23S rRNA. The specificity is changed dramatically if isolated 23S rRNA is used as substrate. We have shown previously that 23S rRNA is digested completely except for two fragments, one of which is identical to the alpha-fragment. Here we show that the other fragment comprises the 5'-end of 23S rRNA and contains 385 nucleotides. A similar fragment was obtained when isolated 23S rRNA was digested with RNase A (specific for pyrimidines in single strands). It appears that the 5'-domain (equivalent to 5.8S rRNA of eukaryotic ribosomes) as well as the 3'-domain (equivalent to 4.5S rRNA of chloroplast ribosomes) have a compact and defined tertiary structure in isolated 23S rRNA in contrast to the rRNA region in between. Thus, alpha-sarcin is a convenient tool for detecting compact domains in isolated RNA.     

Conformational determinants of tandem GU mismatches in RNA: insights from molecular dynamics simulations and quantum mechanical calculations

Structure and energetic properties of base pair mismatches in duplex RNA have been the focus of numerous investigations due to their role in many important biological functions. Such efforts have contributed to the development of models for secondary structure prediction of RNA, including the nearest-neighbor model. In RNA duplexes containing GU mismatches, 5'-GU-3' tandem mismatches have a different thermodynamic stability than 5'-UG-3' mismatches. In addition, 5'-GU-3' mismatches in some sequence contexts do not follow the nearest-neighbor model for stability. To characterize the underlying atomic forces that determine the structural and thermodynamic properties of GU tandem mismatches, molecular dynamics (MD) simulations were performed on a series of 5'-GU-3' and 5'-UG-3' duplexes in different sequence contexts. Overall, the MD-derived structural models agree well with experimental data, including local deviations in base step helicoidal parameters in the region of the GU mismatches and the model where duplex stability is associated with the pattern of GU hydrogen bonding. Further analysis of the simulations, validated by data from quantum mechanical calculations, suggests that the experimentally observed differences in thermodynamic stability are dominated by GG interstrand followed by GU intrastrand base stacking interactions that dictate the one versus two hydrogen bonding scenarios for the GU pairs. In addition, the inability of 5'-GU-3' mismatches in different sequence contexts to all fit into the nearest-neighbor model is indicated to be associated with interactions of the central four base pairs with the surrounding base pairs. The results emphasize the role of GG and GU stacking interactions on the structure and thermodynamics of GU mismatches in RNA.     

500-MHz proton homonuclear Overhauser evidence for additional base pair in the common arm of eukaryotic ribosomal 5S RNA: wheat germ

A "common-arm" fragment from wheat germ (Triticum aestivum) 5S RNA has been produced by enzymatic cleavage with RNase T1 and sequenced via autoradiography of electrophoresis gels for the end-labeled fragments obtained by further RNase T1 partial digestion. The existence, base pair composition, and base pair sequence of the common arm are demonstrated for the first time by means of proton 500-MHz nuclear magnetic resonance. From Mg2+ titration, temperature variation, ring current calculations, sequence comparisons, and proton homonuclear Overhauser enhancement experiments, additional base pairs in the common arm of the eukaryotic 5S RNA secondary structure are detected. Two base pairs, G41 X C34 and A42 X U33 in the hairpin loop, could account for the lack of binding between the conserved GAAC segment of 5S RNA and the conserved Watson-Crick-complementary GT psi C segment of tRNAs.     

Identification of a cytidine-specific ribonuclease from chicken liver

Rapid RNA sequencing technology was used to determine if the base specificities of an RNase recently purified from chicken liver would prove useful for RNA sequence analysis. Escherichia coli 5 S [5'-32P]rRNA or yeast 5.8 S [5'-32P]rRNA was digested with the enzyme and this digest, along with digests derived from RNases of known specificity (U2, T1, T2) were subjected to electrophoresis through denaturing polyacrylamide slab gels. Following autoradiography, the banding patterns arising from the activity of each enzyme were compared, and the base specificity of the unknown RNase was established. The chicken liver RNase was found to have a marked preference for phosphodiester bonds containing cytidylic acid residues, a property which should make the enzyme useful for distinguishing between pyrimidines in RNA sequencing.     

A unique secondary folding pattern for 5S RNA corresponds to the lowest energy homologous secondary structure in 17 different prokaryotes.

A general secondary structure is proposed for the 5S RNA of prokaryotic ribosomes, based on helical energy filtering calculations. We have considered all secondary structures that are common to 17 different prokaryotic 5S RNAs and for each 5S sequence calculated the (global) minimum energy secondary structure (300,000 common structures are possible for each sequence). The 17 different minimum energy secondary structures all correspond, with minor differences, to a single, secondary structure model. This is strong evidence that this general 5S folding pattern corresponds to the secondary structure of the functional 5S rRNA. The general 5S secondary structure is forked and in analogy with the cloverleaf of tRNA is named the "wishbone" model. It constant 8 double helical regions; one in the stem, four in the first, or constant arm, and three in the second arm. Four of these double helical regions are present in a model earlier proposed (1) and four additional regions not proposed by them are presented here. In the minimum energy general structure, the four helices in the constant arm are exactly 15 nucleotide pairs long. These helices are stacked in the sequences from gram-positive bacteria and probably stacked in gram-negative sequences as well. In sequences from gram-positive bacteria the length of the constant arm is maintained at 15 stacked pairs by an unusual minimum energy interaction involving a C26-G57 base pair intercalated between two adjacent helical regions.     

Oligonucleotide directed mutagenesis of Escherichia coli 5S ribosomal RNA: construction of mutant and structural analysis.

The ribosomal 5S RNA gene from the rrnB operon of E. coli was mutagenised in vitro using a synthetic oligonucleotide hybridised to M13 ssDNA containing that gene. The oligonucleotide corresponded to the 5S RNA sequence positions 34 to 51 and changed the guanosine at position 41 to a cytidine. The DNA containing the desired mutation was identified by dot blot hybridisation and introduced back into the plasmid pKK 3535 which contains the total rrnB operon in pBR 322. Plasmid coded 5S rRNA was selectively labeled with 32p using a modified maxi-cell system, and the replacement of guanosine G41 by cytidine was confirmed by RNA sequencing. The growth of cells containing mutant 5S rRNA was not altered by the base change, and the 5S rRNA was processed and incorporated into 50S ribosomal subunits and 70S ribosomes. The structure of wildtype and mutant 5S rRNA was compared by chemical modification of accessible guanosines with kethoxal and limited enzymatic digestion using RNase T1 and nuclease S1. These results showed that the wildtype and mutant 5S rRNA do not differ significantly in their structure. Furthermore, the formation, interconversion and stability of the two 5S rRNA A- and B-conformers are unchanged.     

SARNA-Predict: accuracy improvement of RNA secondary structure prediction using permutation-based simulated annealing

Ribonucleic acid (RNA), a single-stranded linear molecule, is essential to all biological systems. Different regions of the same RNA strand will fold together via base pair interactions to make intricate secondary and tertiary structures that guide crucial homeostatic processes in living organisms. Since the structure of RNA molecules is the key to their function, algorithms for the prediction of RNA structure are of great value. In this article, we demonstrate the usefulness of SARNA-Predict, an RNA secondary structure prediction algorithm based on Simulated Annealing (SA). A performance evaluation of SARNA-Predict in terms of prediction accuracy is made via comparison with eight state-of-the-art RNA prediction algorithms: mfold, Pseudoknot (pknotsRE), NUPACK, pknotsRG-mfe, Sfold, HotKnots, ILM, and STAR. These algorithms are from three different classes: heuristic, dynamic programming, and statistical sampling techniques. An evaluation for the performance of SARNA-Predict in terms of prediction accuracy was verified with native structures. Experiments on 33 individual known structures from eleven RNA classes (tRNA, viral RNA, antigenomic HDV, telomerase RNA, tmRNA, rRNA, RNaseP, 5S rRNA, Group I intron 23S rRNA, Group I intron 16S rRNA, and 16S rRNA) were performed. The results presented in this paper demonstrate that SARNA-Predict can out-perform other state-of-the-art algorithms in terms of prediction accuracy. Furthermore, there is substantial improvement of prediction accuracy by incorporating a more sophisticated thermodynamic model (efn2).     

DNA sequence of the 16S rRNA/23S rRNA intercistronic spacer of two rDNA operons of the archaebacterium Methanococcus vannielii

The DNA sequence of the spacer (plus flanking) regions separating the 16S rRNA and 23S rRNA genes of two presumptive rDNA operons of the archaebacterium Methanococcus vannielii was determined. The spacers are 156 and 242 base pairs in size and they share a sequence homology of 49 base pairs following the 3' terminus of the 16S rRNA gene and of about 60 base pairs preceding the 5' end of the 23S rRNA gene. The 242 base pair spacer, in addition contains a sequence which can be transcribed into tRNAAla, whereas no tRNA-like secondary structure can be delineated from the 156 base pair spacer region. Almost complete sequence homology was detected between the end of the 16S rRNA gene and the 3' termini of either Escherichia coli or Halobacterium halobium 16S rRNA, whereas the putative 5' terminal 23S rRNA sequence shared partial homology with E. coli 23S rRNA and eukaryotic 5.8S rRNA.