The nucleotide sequence of ribosomal 5S RNA from Rhizobium meliloti: comparison with 5S rRNA of Agrobacterium

The complete nucleotide sequence of R. meliloti 5S ribosomal RNA has been determined and compared with the already known sequence of A. tumefaciens 5S rRNA (Vandenberghe et al., 1985, Eur. J. Biochem., 149, 537-542) and of other 5S rRNAs from Rodobacteria Alpha-2 (Wolters et al., 1988, Nucleic Acids Res., 16, rl-r70). The differences found at eight positions (23, 73, 83, 72 in helical fragments; 16, 40, 88 in loops; 54 in bulge), which might affect secondary structures of 5S rRNA, are small. Moreover, the sequence analysis specifies both variable and common positions in 5S rRNA secondary structure of Rodobacteria Alpha-2.     

The secondary structure of the protein L1 binding region of ribosomal 23S RNA. Homologies with putative secondary structures of the L11 mRNA and of a region of mitochondrial 16S rRNA.

An heterologous complex was formed between E. coli protein L1 and P. vulgaris 23S RNA. We determined the primary structure of the RNA region which remained associated with protein L1 after RNase digestion of this complex. We also identified the loci of this RNA region which are highly susceptible to T1, S1 and Naja oxiana nuclease digestions respectively. By comparison of these results with those previously obtained with the homologous regions of E. coli and B. stearothermophilus 23S RNAs, we postulate a general structure for the protein L1 binding region of bacterial 23S RNA. Both mouse and human mit 16S rRNAs and Xenopus laevis and Tetrahymena 28S rRNAs contain a sequence similar to the E. coli 23s RNS region preceding the L1 binding site. The region of mit 16S rRNA which follows this sequence has a potential secondary structure bearing common features with the L1-associated region of bacterial 23S rRNA. The 5'-end region of the L11 mRNA also has several sequence potential secondary structures displaying striking homologies with the protein L1 binding region of 23S rRNA and this probably explains how protein L1 functions as a translational repressor. One of the L11 mRNA putative structures bears the features common to both the L1-associated region of bacterial 23S rRNA and the corresponding region of mit 16S rRNA.     

Prediction of secondary structures of 16S and 23S rRNA fromEscherichia coli

Small and large subunits ofEscherichia coli ribosome have three different rRNAs, the sequences of which are known. However, attempts by three groups to predict secondary structures of 16S and 23S rRNAs have certain common limitations namely, these structures are predicted assuming no interactions among various domains of the molecule and only 40% residues are involved in base pairing as against the experimental observation of 60 % residues in base paired state. Recent experimental studies have shown that there is a specific interaction between naked 16S and 23S rRNA molecules. This is significant because we have observed that the regions (oligonucleotides of length 9–10 residues), in 16S rRNA which are complementary to those in 23S rRNA do not have internal complementary sequences. Therefore, we have developed a simple graph theoretical approach to predict secondary structures of 16S and 23S rRNAs. Our method for model building not only uses complete sequence of 16S or 23S rRNA molecule along with other experimental observations but also takes into account the observation that specific recognition is possible through the complementary sequences between 16S and 23S rRNA molecules and, therefore, these parts of the molecules are not used for internal base pairing. The method used to predict secondary structures is discussed. A typical secondary structure of the complex between 16S and 23S rRNA molecules, obtained using our method, is presented and compared Briefly with earlier model Building studies.     

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.     

The 23S ribosomal RNA higher-order structure of Pseudomonas cepacia and other prokaryotes

A 23S ribosomal RNA gene of Pseudomonas cepacia has been cloned and sequenced. A general higher-order structure model based on earlier published models has been derived from comparative analysis of 23S-like rRNAs of eubacteria, archaebacteria, organelles and eukaryotes. Differences between the previous models were carefully analyzed and controversial regions evaluated. Moderately large insertions and deletions have been found at new points in the secondary structure. The analysis of 50 published as well as unpublished 23S rRNA sequences provide additional proof for six of the seven previously suggested tertiary interactions within the 23S rRNA. P. cepacia is the first representative of the beta subgroup of the Proteobacteria phylum whose 23S rRNA has been sequenced. A tree reflecting evolutionary relationships of prokaryotes was constructed. The topology of this tree is in good agreement with the 16S rRNA tree.     

Secondary structures of proteins from the 30S subunit of the Escherichia coli ribosome.

The secondary structures of the proteins S4, S6, S8, S9, S12, S13, S15, S16, S18, S20 and S21 from the subunit of the E. coli ribosome were predicted according to four different methods. From the resultant diagrams indicating regions of helix, turn, extended structure and random coil, average values for the respective secondary structures could be calculated for each protein. Using the known relative distances for residues in the helical, turn and sheet or allowed random conformations, estimates are made of the maximum possible lengths of the proteins in order to correlate these with results obtained from antibody binding studies to the 30S subunit as determined by electron microscopy. The influence of amino acid changes on the predicted secondary structures of proteins from a few selected mutants was studied. The altered residues tend to be structurally conservative or to induce only minimal local changes.     

The 3''-terminal region of bacterial 23S ribosomal RNA: structure and homology with the 3''-terminal region of eukaryotic 28S rRNA and with chloroplast 4.5s rRNA.

The sequence of the 110 nucleotide fragment located at the 3'-end of E.coli, P.vulgaris and A.punctata 23S rRNAs has been determined. The homology between the E.coli and P.vulgaris fragments is 90%, whereas that between the E.coli and A.punctate fragments is only 60%. The three rRNA fragments have sequences compatible with a secondary structure consisting of two hairpins. Using chemical and enzymatic methods recently developed for the study of the secondary structure of RNA, we demonstrated that one of these hairpins and part of the other are actually present in the three 3'-terminal fragments in solution. This supports the existence of these two hairpins in the intact molecule. Indeed, results obtained upon limited digestion of intact 23S RNA with T1 RNase were in good agreement with the existence of these two hairpins. We observed that the primary structures of the 3'-terminal regions of yeast 26S rRNA and X.laevis 28S rRNA are both compatible with a secondary structure similar to that found at the 3'-end of bacterial 23S rRNAs. Furthermore, both tobacco and wheat chloroplast 4.5S rRNAs can also be folded in a similar way as the 3'-terminal region of bacterial 23S rRNA, the 3'-end of chloroplast 4.5S rRNAs being complementary to the 5'-end of chloroplast 23S rRNA. This strongly reinforces the hypothesis that chloroplast 4.5S rRNA originates from the 3'-end of bacterial 23S rRNA and suggests that this rRNA may be base-paired with the 5'-end of chloroplast 23S rRNA. Invariant oligonucleotides are present at identical positions in the homologous secondary structures of E.coli 23S, yeast 26S, X.laevis 28S and wheat and tobacco 4.5S rRNAs. Surprisingly, the sequences of these oligonucleotides are not all conserved in the 3'-terminal regions of A.punctata or even P.vulgaris 23S rRNAs. Results obtained upon mild methylation of E.coli 50S subunits with dimethylsulfate strongly suggest that these invariant oligonucleotides are involved in RNA tertiary structure or in RNA-protein interactions.     

Sequence of the chloroplast 16S rRNA gene and its surrounding regions of Chlamydomonas reinhardii.

The sequence of a 2 kb DNA fragment containing the chloroplast 16S ribosomal RNA gene from Chlamydomonas reinhardii and its flanking regions has been determined. The algal 16S rRNA sequence (1475 nucleotides) and secondary structure are highly related to those found in bacteria and in the chloroplasts of higher plants. In contrast, the flanking regions are very different. In C. reinhardii the 16S rRNA gene is surrounded by AT rich segments of about 180 bases, which are followed by a long stretch of complementary bases separated from each other by 1833 nucleotides. It is likely that these structures play an important role in the folding and processing of the precursor of 16S rRNA. The primary and secondary structures of the binding sites of two ribosomal proteins in the 16SrRNAs of E. coli and C. reinhardii are considerably related.     

Effects of the ribosomal subunit association on the chemical modification of the 16S and 23S RNAs from Escherichia coli

70S ribosomes and 30S and 50S ribosomal subunits from Escherichia coli were modified under non-denaturing conditions with the chemical reagent dimethylsulfate. The ribosomal 23S and 16S RNAs were isolated after the reaction and the last 200 nucleotides from the 3' ends were analyzed for differences in the chemical modification. A number of accessibility changes could be detected for 23S and 16S RNA when 70S ribosomes as opposed to the isolated subunits were modified. In addition to a number of sites which were protected from modification several guanosines showed enhanced reactivities, indicating conformational changes in the ribosomal RNA structures when 30S and 50S subunits associate to a 70S particle. Most of the accessibility changes can be localized in double-helical regions within the secondary structures of the two RNAs. The results confirm the importance of the ribosomal RNAs for ribosomal functions and help to define the RNA domains which constitute the subunit interface of E. coli ribosomes.     

Localization of the binding site for protein S4 on 16 S ribosomal RNA by chemical and enzymatic probing and primer extension

We have examined the effect of binding ribosomal protein S4 to 16 S rRNA on the susceptibility of the RNA to a variety of chemical and enzymatic probes. We have used dimethyl sulfate to probe unpaired adenines (at N-1) and cytosines (at N-3), kethoxal to probe unpaired guanines (at N-1 and N-2) and cobra venom (V1) ribonuclease as a probe of base-paired regions of 16 S rRNA. Sites of attack by the probes were identified by primer extension using synthetic oligodeoxynucleotides. Comparison of probing results for naked and S4-bound rRNA shows: Protein S4 protects a relatively compact region of the 5' domain of 16 S rRNA from chemical and enzymatic attack. This region is bounded by nucleotides 27 to 47 and 394 to 556, and has a secondary structure characterized by the junction of five helical elements. Phylogenetically conserved irregular features (bulged nucleotides, internal loops and flanking unpaired nucleotides) and helical phosphodiester bonds of four of the helices are specifically protected in the S4-RNA complex. We conclude that this is the major, and possibly sole region of contact between 16 S rRNA and S4. Many of the S4-dependent changes mimic those observed on assembly of 16 S rRNA into 30 S ribosomal subunits. Binding of S4 causes enhanced chemical reactivity coupled with protection from V1 nuclease outside the S4 junction region in the 530, 720 and 1140 loops. We interpret these results as indicative of loss of structure, and suggest that S4 binding causes disruption of adventitious pairing in these regions, possibly by stabilizing the geometry of the RNA such that these interactions are prevented from forming.     

The conformation of 16S RNA in the 30S ribosomal subunit from Escherichia coli.

The digestion of E. coli 16S RNA with a single-strand-specific nuclease produced two fractions separable by gel filtration. One fraction was small oligonucleotides, the other, comprising 67.5% of the total RNA, was highly structured double helical fragments of mol. wt. 7,600. There are thus about 44 helical loops of average size corresponding to 12 base pairs in each 16S RNA. 10% of the RNA could be digested from native 30S subunits. Nuclease attack was primarily in the intraloop single-stranded region but two major sites of attack were located in the interloop single-stranded regions. Nuclease digestion of unfolded subunits produced three classes of fragments, two of which, comprising 80% of the total RNA, were identical to fragments from 16S RNA. The third, consisting of 20% RNA, together with an equal weight of peotein, was a resistant core (sedimentation coefficient 7S).     

Protein binding sites on Escherichia coli 16S ribosomal RNA; RNA regions that are protected by proteins S7, S9 and S19, and by proteins S8, S15 and S17.

Selected groups of isolated 14C-labelled proteins from E. coli 30S ribosomal subunits were reconstituted with 32P-labelled 16S RNA, and the reconstituted complexes were partially digested with ribonuclease A. RNA fragments protected by the proteins were separated by gel electrophoresis and subjected to sequence analysis. Complexes containing proteins S7 and S19 protected an RNA region comprising helices 29 to 32, part of helix 41, and helices 42 and 43 of the 16S RNA secondary structure. Addition of protein S9 had no effect. When compared with previous data for proteins S7, S9, S14 and S19, these results suggest that S14 interacts with helix 33, and that S9 and S14 together interact with the loop-end of helix 41. Complexes containing proteins S8, S15 and S17 protected helices 7 to 10 as well as the "S8-S15 binding site" (helices 20, 22 and parts of helices 21 and 23). When protein S15 was omitted, S8 and S18 showed protection of part of helix 44 in addition to the latter regions. The results are discussed in terms of our model for the detailed arrangement of proteins and RNA in the 30S subunit.     

The hydrodynamic shape, conformation, and molecular model of Escherichia coli ribosomal 5 S RNA.

The structure of ribosomal 5 S RNA has been examined using several physical biochemical techniques. Hydrodynamic measurements yield a s020,omega and [eta] of 5.5 x 10(-13) x and 6.9 ml/g, respectively. Other parameters calculated from these values indicate the shape of 5 S RNA is consistent with that of a prolate ellipsoid 160 A in length and 32 A wide. Sedimentation equilibrium results show that 5 S RNA exists as a monomer in the reconstitution buffer with an apparent molecular weight of 44,000. Ultraviolet absorption difference spectra show that approximately 75% of the bases in 5 S RNA are involved in base pairing, and of these base pairs 70% are G-C and 30% are A-U. These results on the overall shape and secondary structure of 5 S RNA have been incorporated with the results of other investigators as to the possible location of single-stranded and double-stranded helical regions, and a molecular model for 5 S RNA is proposed. The molecular model consists of three double helices in the shape of a prolate ellipsoid, with two of the double helical regions at one end of the molecule. The structure is consistent with the available data on the structure and function of 5 S RNA and bears similarity to the molecular model proposed by Osterberg et al. ((1976) Eur. J. Biochem. 68, 481-487) based on small angle x-ray scattering results and the secondary structure proposed by Madison ((1968) Annu. Rev. Biochem. 37, 131-148).     

Directed hydroxyl radical probing of 16S ribosomal RNA in ribosomes containing Fe(II) tethered to ribosomal protein S20.

The 16S ribosomal RNA neighborhood of ribosomal protein S20 has been mapped, in both 30S subunits and 70S ribosomes, using directed hydroxyl radical probing. Cysteine residues were introduced at amino acid positions 14, 23, 49, and 57 of S20, and used for tethering 1-(p-bromoacetamidobenzyl)-Fe(II)-EDTA. In vitro reconstitution using Fe(II)-derivatized S20, together with the remaining small subunit ribosomal proteins and 16S ribosomal RNA (rRNA), yielded functional 30S subunits. Both 30S subunits and 70S ribosomes containing Fe(II)-S20 were purified and hydroxyl radicals were generated from the tethered Fe(II). Hydroxyl radical cleavage of the 16S rRNA backbone was monitored by primer extension. Different cleavage patterns in 16S rRNA were observed from Fe(II) tethered to each of the four positions, and these patterns were not significantly different in 30S and 70S ribosomes. Cleavage sites were mapped to positions 160-200, 320, and 340-350 in the 5' domain, and to positions 1427-1430 and 1439-1458 in the distal end of the penultimate stem of 16S rRNA, placing these regions near each other in three dimensions. These results are consistent with previous footprinting data that localized S20 near these 16S rRNA elements, providing evidence that S20, like S17, is located near the bottom of the 30S subunit.     

Prediction of RNA secondary structure based on helical regions distribution

MOTIVATION: RNAs play an important role in many biological processes and knowing their structure is important in understanding their function. Due to difficulties in the experimental determination of RNA secondary structure, the methods of theoretical prediction for known sequences are often used. Although many different algorithms for such predictions have been developed, this problem has not yet been solved. It is thus necessary to develop new methods for predicting RNA secondary structure. The most-used at present is Zuker's algorithm which can be used to determine the minimum free energy secondary structure. However many RNA secondary structures verified by experiments are not consistent with the minimum free energy secondary structures. In order to solve this problem, a method used to search a group of secondary structures whose free energy is close to the global minimum free energy was developed by Zuker in 1989. When considering a group of secondary structures, if there is no experimental data, we cannot tell which one is better than the others. This case also occurs in combinatorial and heuristic methods. These two kinds of methods have several weaknesses. Here we show how the central limit theorem can be used to solve these problems. RESULTS: An algorithm for predicting RNA secondary structure based on helical regions distribution is presented, which can be used to find the most probable secondary structure for a given RNA sequence. It consists of three steps. First, list all possible helical regions. Second, according to central limit theorem, estimate the occurrence probability of every helical region based on the Monte Carlo simulation. Third, add the helical region with the biggest probability to the current structure and eliminate the helical regions incompatible with the current structure. The above processes can be repeated until no more helical regions can be added. Take the current structure as the final RNA secondary structure. In order to demonstrate the confidence of the program, a test on three RNA sequences: tRNAPhe, Pre-tRNATyr, and Tetrahymena ribosomal RNA intervening sequence, is performed. AVAILABILITY: The program is written in Turbo Pascal 7.0. The source code is available upon request. CONTACT: Wujj@nic.bmi.ac.cn or Liwj@mail.bmi.ac.cn