Prediction of the Recognition Sites on 16S and 23S rRNAs from E.coli for the Formation of 16S-23S rRNA Complex

Abstract

Interactions between RNA molecules have been postulated to play an important role in the assembly of ribosomes. Using the sequence analysis and the search of continuous complementary regions on 16S rRNA and 23S rRNA, the recognition sites involved in the formation of ribosome of E.coli are postulated. The number of postulated sites was narrowed down by taking available experimental data. The suggestive evidence for correct postulation is obtained from sequence comparison studies of 16S and 23S rRNAs from various species. The sites 891–899 and 1195–1203 on 16S rRNA along with the corresponding complementary sites 1904–1912 and 760–768 on 23S rRNA are predicted to be the most probable candidates for the sites of recognition between 16S and 23S rRNAs. The possibility of the involvement of the additional site 630–638 on 16S rRNA with its complementary site 2031–2039 on 23S rRNA cannot be ruled out.     

The accuracy of ribosomal RNA comparative structure models

The determination of the 16S and 23S rRNA secondary structure models was initiated shortly after the first complete 16S and 23S rRNA sequences were determined in the late 1970s. The structures that are common to all 16S rRNAs and all 23S rRNAs were determined using comparative methods from the analysis of thousands of rRNA sequences. Twenty-plus years later, the 16S and 23S rRNA comparative structure models have been evaluated against the recently determined high-resolution crystal structures of the 30S and 50S ribosomal subunits. Nearly all of the predicted covariation-based base pairs, including the regular base pairs and helices, and the irregular base pairs and tertiary interactions, were present in the 30S and 50S crystal structures.     

Prediction of the recognition sites on 16S and 23S rRNAs from E. coli for the formation of 16S-23S rRNA complex

Interactions between RNA molecules have been postulated to play an important role in the assembly of ribosomes. Using the sequence analysis and the search of continuous complementary regions on 16S rRNA and 23S rRNA, the recognition sites involved in the formation of ribosome of E. coli are postulated. The number of postulated sites was narrowed down by taking available experimental data. The suggestive evidence for correct postulation is obtained from sequence comparison studies of 16S and 23S rRNAs from various species. The sites 891-899 and 1195-1203 on 16S rRNA along with the corresponding complementary sites 1904-1912 and 760-768 on 23S rRNA are predicted to be the most probable candidates for the sites of recognition between 16S and 23S rRNAs. The possibility of the involvement of the additional site 630-638 on 16S rRNA with its complementary site 2031-2039 on 23S rRNA cannot be ruled out.     

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.     

tRNA genes are found between 16S and 23S rRNA genes in Bacillus subtilis.

There are at least nine, and probably ten, ribosomal RNA gene sets in the genome of Bacillus subtilis. Each gene set contains sequences complementary to 16S, 23S and 5S rRNAs. We have determined the nucleotide sequences of two DNA fragments which each contain 165 base pairs of the 16S rRNA gene, 191 base pairs of the 23S rRNA gene, and the spacer region between them. The smaller space region is 164 base pairs in length and the larger one includes an additional 180 base pairs. The extra nucleotides could be transcribed in tRNAIIe and tRNA Ala sequences. Evidence is also presented for the existence of a second spacer region which also contains tRNAIIe and tRNA Ala sequences. No other tRNAs appear to be encoded in the spacer regions between the 16S and 23S rRNA genes. Whereas the nucleotide sequences corresponding to the 16S rRNA, 23S rRNA and the spacer tRNAs are very similar to those of E. coli, the sequences between these structural genes are very different.     

Computer analysis of potential stem structures of rRNA operons in various procaryote genomes

The processing of 16S rRNA and 23S rRNA by RNase III in E.coli is known to involve stem structures formed by both ends of the rRNA. Indeed, complementary nucleotide sequences are usually found at both ends of 16S rRNA and 23S rRNA. However, whether or not this phenomenon exists in various other bacteria has not yet been adequately studied. We have conducted computer analyses of potential stem structures of rRNA operons in 12 bacterial and 3 archaeal genomes, and compared characteristics of the stem structures among these species. We systematically computed free energy values by exhaustively 'annealing' sequences around the 5' end and sequences around the 3' end of both 16S rRNA and 23S rRNA genes, in order to predict potential stem structures.The results suggest that rRNAs in most species form stem structures at both ends. Some species, such as A.aeolicus, seem to form unusually stable stem structures. On the other hand, some rRNAs, such as rRNAs of D.radiodurans, seem not to form solid stem structures. This suggests that rRNA processing in those species must employ a reliable targeting mechanism other than recognizing stem structures by RNase III.     

The secondary structure of human 28S rRNA: The structure and evolution of a mosaic rRNA gene

Summary We have determined the secondary structure of the human 28S rRNA molecule based on comparative analysis of available eukaryotic cytoplasmic and prokaryotic large-rRNA gene sequences. Examination of large-rRNA sequences of both distantly and closely related species has enabled us to derive a structure that accounts both for highly conserved sequence tracts and for previously unanalyzed variable-sequence tracts that account for the evolutionary differences in size among the large rRNAs.Human 28S rRNA is composed of two different types of sequence tracts: conserved and variable. They differ in composition, degree of conservation, and evolution. The conserved regions demonstrate a striking constancy of size and sequence. We have confirmed that the conserved regions of large-rRNA molecules are capable of forming structures that are superimposable on one another. The variable regions contain the sequences responsible for the 83% increase in size of the human large-rRNA molecule over that ofEscherichia coli. Their locations in the gene are maintained during evolution. They are G+C rich and largely nonhomologous, contain simple repetitive sequences, appear to evolve by frequent recombinational events, and are capable of forming large, stable hairpins.The secondary-structure model presented here is in close agreement with existing prokaryotic 23S rRNA secondary-structure models. The introduction of this model helps resolve differences between previously proposed prokaryotic and eukaryotic large-rRNA secondary-structure models.     

Secondary structure model for bacterial 16S ribosomal RNA: phylogenetic, enzymatic and chemical evidence.

We have derived a secondary structure model for 16S ribosomal RNA on the basis of comparative sequence analysis, chemical modification studies and nuclease susceptibility data. Nucleotide sequences of the E. coli and B. brevis 16S rRNA chains, and of RNAse T1 oligomer catalogs from 16S rRNAs of over 100 species of eubacteria were used for phylogenetic comparison. Chemical modification of G by glyoxal, A by m-chloroperbenzoic acid and C by bisulfite in naked 16S rRNA, and G by kethoxal in active and inactive 30S ribosomal subunits was taken as an indication of single stranded structure. Further support for the structure was obtained from susceptibility to RNases A and T1. These three approaches are in excellent agreement. The structure contains fifty helical elements organized into four major domains, in which 46 percent of the nucleotides of 16S rRNA are involved in base pairing. Phylogenetic comparison shows that highly conserved sequences are found principally in unpaired regions of the molecule. No knots are created by the structure.     

Enhancement of translation by the epsilon element is independent of the sequence of the 460 region of 16S rRNA

The epsilon enhancer element is a pyrimidine-rich sequence that increases expression of T7 gene 10 and a number of Escherichia coli mRNAs during initiation of translation and inhibits expression of the recF mRNA during elongation. Based on its complementarity to the 460 region of 16S rRNA, it has been proposed that epsilon exerts its enhancer activity by base pairing to this complementary rRNA sequence. We have tested this model of enhancer action by constructing mutations in the 460 region of 16S rRNA and examining expression of epsilon-containing CAT reporter genes and recF-lacZ fusions in strains expressing the mutant rRNAs. Replacement of the 460 E.coli stem-loop with that of Salmonella enterica serovar Typhimurium or a stem-loop containing a reversal of all 8 bp in the helical region produced fully functional rRNAs with no apparent effect on cell growth or expression of any epsilon-containing mRNA. Our experiments confirm the reported effects of the epsilon elements on gene expression but show that these effects are independent of the sequence of the 460 region of 16S rRNA, indicating that epsilon-rRNA base pairing does not occur.     

In vivo analysis of plant RNA structure: soybean 18S ribosomal and ribulose-1,5-bisphosphate carboxylase small subunit RNAs

A method to investigate the structure of RNA molecules within intact plant tissues has been developed. The RNA structures are analyzed using dimethyl sulfate (DMS), which modifies substituents of adenine and cytosine residues within single-stranded regions of RNA molecules. Reactive sites are identified by primer extension analysis. Using this procedure, an analysis of the secondary structure of the cytoplasmic 18S ribosomal RNA in soybean seedling leaves has been completed. DMS modification data are in good agreement with the phylogenetic structure predicted for soybean 18S rRNA. However, there are a few notable exceptions where residues thought to be involved in double-stranded regions in all 18S rRNAs are strongly modified in soybean leaf samples. These data taken together with the phylogenetic structure suggest that alternate structures may exist in vivo.The further applicability of this technique is demonstrated by comparing the modification pattern obtained in vivo to that obtained in vitro for a particular mRNA molecule encoding the small subunit of ribulose-1,5-bisphosphate carboxylase. The results obtained are compared to a predicted minimum energy secondary structure. The data indicate that the conformation of RNA molecules within the cell may not be reflected in a structural analysis of purified mRNA molecules.     

Model-Free RNA Sequence and Structure Alignment Informed by SHAPE Probing Reveals a Conserved Alternate Secondary Structure for 16S rRNA

Discovery and characterization of functional RNA structures remains challenging due to deficiencies in de novo secondary structure modeling. Here we describe a dynamic programming approach for model-free sequence comparison that incorporates high-throughput chemical probing data. Based on SHAPE probing data alone, ribosomal RNAs (rRNAs) from three diverse organisms – the eubacteria E. coli and C. difficile and the archeon H. volcanii – could be aligned with accuracies comparable to alignments based on actual sequence identity. When both base sequence identity and chemical probing reactivities were considered together, accuracies improved further. Derived sequence alignments and chemical probing data from protein-free RNAs were then used as pseudo-free energy constraints to model consensus secondary structures for the 16S and 23S rRNAs. There are critical differences between these experimentally-informed models and currently accepted models, including in the functionally important neck and decoding regions of the 16S rRNA. We infer that the 16S rRNA has evolved to undergo large-scale changes in base pairing as part of ribosome function. As high-quality RNA probing data become widely available, structurally-informed sequence alignment will become broadly useful for de novo motif and function discovery.     

Intermolecular base-paired interaction between complementary sequences present near the 3'' ends of 5S rRNA and 18S (16S) rRNA might be involved in the reversible association of ribosomal subunits.

Highly conserved sequences present at an identical position near the 3' ends of eukaryotic and prokaryotic 5S rRNAs are complementary to the 5' strand of the m2(6)A hairpin structure near the 3' ends of 18S rRNA and 16S rRNA, respectively. The extent of base-pairing and the calculated stabilities of the hybrids that can be constructed between 5S rRNAs and the small ribosomal subunit RNAs are greater than most, if not all, RNA-RNA interactions that have been implicated in protein synthesis. The existence of complementary sequences in 5S rRNA and small ribosomal subunit RNA, along with the previous observation that there is very efficient and selective hybridization in vitro between 5S and 18S rRNA, suggests that base-pairing between 5S rRNA in the large ribosomal subunit and 18S (16S) rRNA in the small ribosomal subunit might be involved in the reversible association of ribosomal subunits. Structural and functional evidence supporting this hypothesis is discussed.     

一种预测单链RNA分子二级结构的计算机算法

本文给出了一个利用已知能量数据构成具有最小自由能的单链RNA分子二级结构的计算机算法,并给出了此算法的可行性证明和应用实例。     

The sequence of the 5.8 S ribosomal RNA of the crustacean Artemia salina. With a proposal for a general secondary structure model for 5.8 S ribosomal RNA.

We report the primary structure of 5.8 S rRNA from the crustacean Artemia salina. The preparation shows length heterogeneity at the 5'-terminus, but consists of uninterrupted RNA chains, in contrast to some insect 5.8 S rRNAs, which consist of two chains of unequal length separated in the gene by a short spacer. The sequence was aligned with those of 11 other 5.8 S rRNAs and a general secondary structure model derived. It has four helical regions in common with the model of Nazar et al. (J. Biol. Chem. 250, 8591-8597 (1975)), but for a fifth helix a different base pairing scheme was found preferable, and the terminal sequences are presumed to bind to 28 S rRNA instead of binding to each other. In the case of yeast, where both the 5.8 S and 26 S rRNA sequences are known, the existence of five helices in 5.8 S rRNA is shown to be compatible with a 5.8 S - 26 S rRNA interaction model.