A story: unpaired adenosine bases in ribosomal RNAs

In 1985 an analysis of the Escherichia coli 16 S rRNA covariation-based structure model revealed a strong bias for unpaired adenosines. The same analysis revealed that the majority of the G, C, and U bases were paired. These biases are (now) consistent with the high percentage of unpaired adenosine nucleotides in several structure motifs.An analysis of a larger set of bacterial comparative 16 S and 23 S rRNA structure models has substantiated this initial finding and revealed new biases in the distribution of adenosine nucleotides in loop regions. The majority of the adenosine nucleotides are unpaired, while the majority of the G, C, and U bases are paired in the covariation-based structure model. The unpaired adenosine nucleotides predominate in the middle and at the 3' end of loops, and are the second most frequent nucleotide type at the 5' end of loops (G is the most common nucleotide). There are additional biases for unpaired adenosine nucleotides at the 3' end of loops and adjacent to a G at the 5' end of the helix. The most prevalent consecutive nucleotides are GG, GA, AG, and AA. A total of 70 % of the GG sequences are within helices, while more than 70 % of the AA sequences are unpaired. Nearly 50 % of the GA sequences are unpaired, and approximately one-third of the AG sequences are within helices while another third are at the 3' loop.5' helix junction. Unpaired positions with an adenosine nucleotide in more than 50 % of the sequences at the 3' end of 16 S and 23 S rRNA loops were identified and arranged into the A-motif categories XAZ, AAZ, XAG, AAG, and AAG:U, where G or Z is paired, G:U is a base-pair, and X is not an A and Z is not a G in more than 50 % of the sequences. These sequence motifs were associated with several structural motifs, such as adenosine platforms, E and E-like loops, A:A and A:G pairings at the end of helices, G:A tandem base-pairs, GNRA tetraloop hairpins, and U-turns.     

Diversity of base-pair conformations and their occurrence in rRNA structure and RNA structural motifs

In addition to the canonical base-pairs comprising the standard Watson-Crick (C:G and U:A) and wobble U:G conformations, an analysis of the base-pair types and conformations in the rRNAs in the high-resolution crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits has identified a wide variety of non-canonical base-pair types and conformations. However, the existing nomenclatures do not describe all of the observed non-canonical conformations or describe them with some ambiguity. Thus, a standardized system is required to classify all of these non-canonical conformations appropriately. Here, we propose a new, simple and systematic nomenclature that unambiguously classifies base-pair conformations occurring in base-pairs, base-triples and base-quadruples that are associated with secondary and tertiary interactions. This system is based on the topological arrangement of the two bases and glycosidic bonds in a given base-pair. Base-pairs in the internal positions of regular secondary structure helices usually form with canonical base-pair groups (C:G, U:A, and U:G) and canonical conformations (C:G WC, U:A WC, and U:G Wb). In contrast, non-helical base-pairs outside of regular structure helices usually have non-canonical base-pair groups and conformations. In addition, many non-helical base-pairs are involved in RNA motifs that form a defined set of non-canonical conformations. Thus, each rare non-canonical conformation may be functionally and structurally important. Finally, the topology-based isostericity of base-pair conformations can rationalize base-pair exchanges in the evolution of RNA molecules.     

The Genesis of Ribosome Structure: How a Protein Generates RNA Structure in Real Time

Ribosomal subunit assembly is initiated by the binding of several primary binding proteins. Results from chemical modification studies show that 16S ribosomal RNA undergoes striking structural rearrangements when protein S17 is bound. For the first time, we are able to distinguish and order these structural rearrangements by using time-dependent chemical probing. Initially, protein S17 binds to a portion of helix 11, inducing a kink-turn in that helix that bends helix 7 toward the S17-helix 11 complex in a hairpin-like manner, allowing helix 7 to bind to protein S17. This structural change is rapidly stabilized by interactions at the distal and proximal ends of both RNA helices. Identifying the dynamic nature of interactions between RNA and proteins is not only essential in unraveling ribosome assembly, but also has more general application to all protein-RNA interactions.     

Unique tertiary and neighbor interactions determine conservation patterns of Cis Watson-Crick A/G base-pairs

X-ray, phylogenetic and quantum chemical analysis of molecular interactions and conservation patterns of cis Watson-Crick (W.C.) A/G base-pairs in 16S rRNA, 23S rRNA and other molecules was carried out. In these base-pairs, the A and G nucleotides interact with their W.C. edges with glycosidic bonds oriented cis relative to each other. The base-pair is stabilised by two hydrogen bonds, the C1'-C1' distance is enlarged and the G(N2) amino group is left unpaired. Quantum chemical calculations show that, in the absence of other interactions, the unpaired amino group is substantially non-planar due to its partial sp(3) pyramidalization, while the whole base-pair is internally propeller twisted and very flexible. The unique molecular properties of the cis W.C. A/G base-pairs make them distinct from other base-pairs. They occur mostly at the ends of canonical helices, where they serve as interfaces between the helix and other motifs. The cis W.C. A/G base-pairs play crucial roles in natural RNA structures with salient sequence conservation patterns. The key contribution to conservation is provided by the unpaired G(N2) amino group that is involved in a wide range of tertiary and neighbor contacts in the crystal structures. Many of them are oriented out of the plane of the guanine base and utilize the partial sp(3) pyramidalization of the G(N2). There is a lack of A/G to G/A covariation, which, except for the G(N2) position, would be entirely isosteric. On the contrary, there is a rather frequent occurrence of G/A to G/U covariation, as the G/U wobble base-pair has an unpaired amino group in the same position as the cis W.C. G/A base-pair. The cis W.C. A/G base-pairs are not conserved when there is no tertiary or neighbor interaction. Obtaining the proper picture of the interactions and phylogenetic patterns of the cis W.C. A/G base-pairs requires a detailed analysis of the relation between the molecular structures and the energetics of interactions at a level of single H-bonds and contacts.     

Singly and bifurcated hydrogen-bonded base-pairs in tRNA anticodon hairpins and ribozymes.

The tRNA anticodon loops always comprise seven nucleotides and is involved in many recognition processes with proteins and RNA fragments. We have investigated the nature and the possible interactions between the first (32) and last (38) residues of the loop on the basis of the available sequences and crystal structures. The data demonstrate the conservation of a bifurcated hydrogen bond interaction between residues 32 and 38, located at the stem/loop junction. This interaction leads to the formation of a non-canonical base-pair which is preserved in the known crystal structures of tRNA/synthetase complexes. Among the tRNA and tDNA sequences, 93 % of the 32.38 oppositions can be assigned to two families of isosteric base-pairs, one with a large (86 %) and the other with a much smaller (7 %) population. The remainder (7 %) of the oppositions have been assigned to a third family due to the lack of evidence for assigning them into the first two sets. In all families, the Y32.R38 base-pairs are not isosteric upon reversal (like the sheared G.A or wobble G.U pairs), explaining the strong conservation of a pyrimidine at position 32. Thus, the 32.38 interaction extends the sequence signature of the anticodon loop beyond the conserved U-turn at position 33 and the usually modified purine at position 37. A comparison with other loops containing both a singly hydrogen-bonded base-pair and a U-turn suggests that the 32.38 pair could be involved in the formation of a base triple with a residue in a ribosomal RNA component. It is also observed that two crystal structures of ribozymes (hammerhead and leadzyme) present similar base-pairs at the cleavage site.     

Identification and assignment of base pairs in three helical stems of Torulopsis utilis ribosomal 5S RNA and its RNase T1 and RNase T2 cleaved fragments via 500-MHz proton homonuclear overhauser enhancements

Imino proton resonances in the downfield region (10-14 ppm) of the 500-MHz 1H NMR spectrum of Torulopsis utilis 5S RNA are identified (A X U, G X C, or G X U) and assigned to base pairs in helices I, IV, and V via analysis of homonuclear Overhauser enhancements (NOE) from intact T. utilis 5S RNA, its RNase T1 and RNase T2 digested fragments, and a second yeast (Saccharomyces cerevisiae) 5S RNA whose nucleotide sequence differs at only six residues from that of T. utilis 5S RNA. The near-identical chemical shifts and NOE behavior of most of the common peaks from these four RNAs strongly suggest that helices I, IV, and V retain the same conformation after RNase digestion and that both T. utilis and S. cerevisiae 5S RNAs share a common secondary and tertiary structure. Of the four G X U base pairs identified in the intact 5S RNA, two are assigned to the terminal stem (helix I) and the other two to helices IV and V. Seven of the nine base pairs of the terminal stem have been assigned. Our experimental demonstration of a G X U base pair in helix V supports the 5S RNA secondary structural model of Luehrsen and Fox [Luehrsen, K. R., & Fox, G.E. (1981) Proc. Natl. Acad. Sci. U.S.A. 78, 2150-2154]. Finally, the base-pair proton peak assigned to the terminal G X U in helix V of the RNase T2 cleaved fragment is shifted downfield from that in the intact 5S RNA, suggesting that helices I and V may be coaxial in intact T. utilis 5S RNA.     

RNA structural dynamics: pre-melting and melting transitions in E. coli 5S rRNA

The temperature dependent transition from duplex to a single strand in E. coli 5S ribosomal RNA is a multistep process, and it involves intermediate states. We have analyzed these structural dynamics by chemical modification of cytidines and by single strand specific nuclease digestions. This combined approach led to the characterization of premelting and melting transitions within individual structural segments of the native macromolecule, which we feel may find general application to the structure of biological polyribonucleotides: 1) G-C base pairs at the termini of helices are relatively unstable and they readily undergo premelting transition. 2) Internal G-U/A-U rich stretches of helices exhibit dynamic premelting properties. 3) Hairpin loops have a relatively stronger destabilizing effect than internal loops. 4) Bulge loops destabilize the neighbouring base pairs. 5) Melting of helical segments occurs starting from the destabilizing structures listed above, preferentially from the helix termini. E. coli 5S rRNA has been shown to adopt different conformations. The presence of urea leads to induction of enhancement in the sensitivity for nuclease S1 at several nucleotide positions. The possibility of structural rearrangements will be discussed.     

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.     

Identification and assignment of base pairs in four helical segments of Bacillus megaterium ribosomal 5S RNA and its ribonuclease T1 cleavage fragments by means of 500-MHz proton homonuclear Overhauser enhancements.

Three different fragments of Bacillus megaterium ribosomal 5S RNA have been produced by enzymatic cleavage with ribonuclease T1. Fragment A consists of helices II and III, fragment B contains helix IV, and fragment C contains helix I of the universal 5S rRNA secondary structure. All (eight) imino proton resonances in the downfield region (9-15 ppm) of the 500-MHz proton FT NMR spectrum of fragment B have been identified and assigned as G80.C92-G81.C91-G82.C90-A83.++ +U89-C84.G88 and three unpaired U's (U85, U86, and U87) in helix IV by proton homonuclear Overhauser enhancement connectivities. The secondary structure in helix IV of the prokaryotic loop is completely demonstrated spectroscopically for the first time in any native or enzyme-cleaved 5S rRNA. In addition, G21.C58-A20.U59-G19.C60-A18.U61 in helix II, U32.A46-G31.C47-C30.G48-C29.G49 in helix III, and G4.C112-G5.C111-U6.G110 in the terminal stem (helix I) have been assigned by means of NOE experiments on intact 5S rRNA and its fragments A and C. Base pairs in helices I-IV of the universal secondary structure of B. megaterium 5S RNA are described.     

Structure and function of the conserved 690 hairpin in Escherichia coli 16 S ribosomal RNA: analysis of the stem nucleotides

Nucleotides 680 to 710 of Escherichia coli 16 S rRNA form a distinct structural domain required for ribosome function. The goal of this study was to determine the functional significance of pairing interactions in the 690 region. Two different secondary structures were proposed for this hairpin, based on phylogenetic and chemical modification studies. To study the effect of pairing interactions in the 690 hairpin on ribosome function and to determine which of the proposed secondary structures is biologically significant, we performed an instant-evolution experiment in which the nine nucleotides that form the proposed base-pairs and dangling ends of the 690 stem were randomly mutated, and functional mutant combinations were selected. A total of 96 unique functional mutants were isolated, assayed in vivo, and sequenced. Analysis of these data revealed extensive base-pairing and stacking interactions among the mutated nucleotides. Formation of either a Watson-Crick base-pair or G.U pair between positions 688 and 699 is absolutely required for ribosome function. We also performed NMR studies of a 31-nucleotide RNA which indicate the formation of a functionally important base-pair between nucleotides 688 and 699. Formation of a second base-pair between positions 689 and 698, however, is not essential for ribosome function, but the level of ribosome function correlates with the predicted thermodynamic stability of the nucleotide pairs in these positions. The universally conserved positions G690 and U697 are generally portrayed as forming a G.U mismatch. Our data show co-variation between these positions, but do not support the hypothesis that the G690:U697 pair forms a wobble structure. NMR studies of model 14-nt and 31-nt RNAs support these findings and show that G690 and U697 are involved in unusual stacking interactions but do not form a wobble pair. Preliminary NMR structural analysis reveals that the loop portion of the 690 hairpin folds into a highly structured and novel conformation.     

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.     

Solution probing of metal ion binding by helix 27 from Escherichia coli 16S rRNA

Helix (H)27 from Escherichia coli 16S ribosomal (r)RNA is centrally located within the small (30S) ribosomal subunit, immediately adjacent to the decoding center. Bacterial 30S subunit crystal structures depicting Mg(2+) binding sites resolve two magnesium ions within the vicinity of H27: one in the major groove of the G886-U911 wobble pair, and one within the GCAA tetraloop. Binding of such metal cations is generally thought to be crucial for RNA folding and function. To ask how metal ion-RNA interactions in crystals compare with those in solution, we have characterized, using solution NMR spectroscopy, Tb(3+) footprinting and time-resolved fluorescence resonance energy transfer (tr-FRET), location, and modes of metal ion binding in an isolated H27. NMR and Tb(3+) footprinting data indicate that solution secondary structure and Mg(2+) binding are generally consistent with the ribosomal crystal structures. However, our analyses also suggest that H27 is dynamic in solution and that metal ions localize within the narrow major groove formed by the juxtaposition of the loop E motif with the tandem G894-U905 and G895-U904 wobble pairs. In addition, tr-FRET studies provide evidence that Mg(2+) uptake by the H27 construct results in a global lengthening of the helix. We propose that only a subset of H27-metal ion interactions has been captured in the crystal structures of the 30S ribosomal subunit, and that small-scale structural dynamics afforded by solution conditions may contribute to these differences. Our studies thus highlight an example for differences between RNA-metal ion interactions observed in solution and in crystals.     

The structure of helix III in Xenopus oocyte 5 S rRNA: an RNA stem containing a two-nucleotide bulge

The solution structure of an oligonucleotide containing the helix III sequence from Xenopus oocyte 5 S rRNA has been determined by NMR spectroscopy. Helix III includes two unpaired adenosine residues, flanked on either side by G:C base-pairs, that are required for binding of ribosomal protein L5. The consensus conformation of helix III in the context provided by this oligonucleotide has the two adenosine residues located in the minor groove and stacked upon the 3' flanking guanosine residue, consistent with biochemical studies of free 5 S rRNA in solution. A distinct break in stacking that occurs between the first adenosine residue of the bulge and the flanking 5' guanosine residue exposes the base of the adenosine residue in the minor groove and the base of the guanosine residue in the major groove. The major groove of the helix is widened at the site of the unpaired nucleotides and the helix is substantially bent; nonetheless, the G:C base-pairs flanking the bulge are intact. The data indicate that there may be conformational heterogeneity centered in the bulge region. The corresponding adenosine residues in the Haloarcula marismortui 50 S ribosomal subunit form a dinucleotide platform, which is quite different from the motif seen in solution. Thus, the conformation of helix III probably changes when 5 S rRNA is incorporated into the ribosome.     

A hierarchy of RNA subdomains in assembly of the central domain of the 30 S ribosomal subunit

Beginning with the framework that has been developed for the assembly of the 30 S ribosomal subunit, we have identified a series of RNAs that are minimal binding sites for proteins S15, S6, S18, and S11 in the central domain from Thermus thermophilus. The minimal binding RNA for proteins S15, S6, and S18 consists of helix 22 and three-way junctions at both ends composed of portions of helices 20, 21, and 23. Addition of the remaining portion of helix 23 to this construct results in the minimal site for S11. Surprisingly, almost half of the central domain (helices 24, 25, and 26) is dispensable for binding the central domain proteins. Thus, at least two classes of RNA elements can be identified in ribosomal RNA. A protein-binding core element (such as helices 20, 21, 22, and 23) is required for the association of ribosomal proteins, whereas secondary binding elements (such as helices 24, 25, and 26) associate only with the preformed core RNP complex. Apparently, there may be a hierarchy of ribosomal RNA elements similar to the hierarchy of primary, secondary, and tertiary binding ribosomal proteins.     

5S RNA structure and interaction with transcription factor A. 1. Ribonuclease probe of the structure of 5S RNA from Xenopus laevis oocytes

The structure of Xenopus laevis oocyte (Xlo) 5S ribosomal RNA has been probed with single-strand-specific ribonucleases T1, T2, and A with double-strand-specific ribonuclease V1 from cobra venom. The digestion of 5'- or 3'-labeled renatured 5S RNA samples followed by gel purification of the digested samples allowed the determination of primary cleavage sites. Results of these ribonuclease digestions provide support for the generalized 5S RNA secondary structural model derived from comparative sequence analysis. However, three putative single-stranded regions of the molecule exhibited unexpected V1 cuts, found at C36, U73, U76, and U102. These V1 cuts reflect additional secondary structural features of the RNA including A.G base pairs and support the extended base pairing in the stem containing helices IV and V which was proposed by Stahl et al. [Stahl, D. A., Luehrsen, K. R., Woese, C. R., & Pace, N. R. (1981) Nucleic Acids Res. 9, 6129-6137]. A conserved structure for helix V having a common unpaired uracil residue at Xlo position 84 is proposed for all eukaryotic 5S RNAs. Our results are compared with nuclease probes of other 5S RNAs.     

Groove-backbone interaction in B-DNA. Implication for DNA condensation and recombination

DNA self-fitting is revealed by the study of intermolecular contacts found in the crystal packing of a dodecamer where the helices are locked together by a reciprocal groove-backbone interaction and form a crossed structure. It is proposed that it could be a model for DNA-DNA interaction in several biological processes such as the node of supercoiled DNA and synapsis in recombination. The main topological and symmetrical features of this crossed structure are described and the symmetry-homology relationships are analyzed in the more general case of B-DNA interacting helices. Model-building of Holliday junctions with minimal change from the starting crystal coordinates of the crossed structure leads to at least three different solutions. These various models are compared from the point of view of their symmetry and topology, in the light of their branch migration and resolution properties. In addition, a model for a self-favored reciprocal unwinding mechanism based on the experimentally observed structural alterations, such as the packing-induced opening of G.C base-pairs is proposed. In this model, the phosphate groups of the invading backbone trigger the opening of the base-pairs of the other helix, by pulling cytosine or adenine bases out of the major groove after binding to their amino group.     

An estimate of the nearest neighbor base-pair content of 5S RNA using CD and absorption spectroscopy.

We analyzed the CD and uv absorption spectra of 5S RNA from Escherichia coli using the method developed in the preceding paper. The analysis of spectra of 5S RNA at 20 degrees C in 0.1M NaClO4, 2.5 mM Na+ (phosphate), pH 7.0, and 0.5 mM MgSO4 gave 7 +/- 3.6 A.U base pairs, 25 +/- 3.6 G.C base pairs, and 7.5 +/- 3.6 G.U base pairs. Estimates of nearest neighbor base pairs were more consistent with the Pieler-Erdmann and the Gewirth-Moore secondary structure models than with the Fox-Woese or the Burns-Luoma-Marshall models. We also examined the structure of 5S RNA as a function of temperature. The melting profile exhibited two transitions--one at about 35 degrees C and one above 50 degrees C. Our spectral data showed that helices I and II were stable during the first transition, and agreed with other data that helix III was the most likely helix to have melted. The results from this in-depth study of 5S RNA indicate that our method of analysis should be useful for studying the secondary structures of other small, unmodified RNAs.     

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)