An analysis of G-U base pair occurrence in eukaryotic 5S rRNAs

The structure-function relationship in RNA molecules is a key to understanding of the expression of genetic information. Various types of RNA play crucial roles at almost every step of protein biosynthesis. In recent years, it has been shown that one of the most important structural elements in RNA is a wobble pair G-U. In this paper, we present for the first time an analysis of the distribution of G-U pairs in eukaryotic 5S ribosomal RNAs. Interestingly, the G-U pair in 5S rRNA species is predominantly found in two intrahelical regions of the stems I and V and at the junction of helix IV and loop A. The distribution of G-U pairs and the nature of adjacent bases suggests their possible role as a recognition site in interactions with other components of protein biosynthesis machinery.     

Stacking of Crick Wobble pair and Watson-Crick pair: stability rules of G-U pairs at ends of helical stems in tRNAs and the relation to codon-anticodon Wobble interaction.

The occurrence of the noncomplementary G-U base pair at the end of a helix is found to be governed by stacking interactions. As a rule, a G-U pair with G on the 5'-side of a Watson-Crick base pair exhibits strikingly greater stacking overlap with the Watson-Crick base pair than a G-U pair on the 3'-side of a Watson-Crick base pair. The former arrangement is expected to be more stable and indeed is observed 29 times out of 32 in the known transfer RNA molecules. In accordance with this rule, the major wobble base pairs G-U or I-U in codon-anticodon interactions have G or I on the 5'-side of the anticodon. Similarly, in initiator tRNAs, this rule is obeyed where now the G is the first letter of the codon (5'-side). In the situation where U is in the wobble position of the anticodon, it is usually substituted at C(5) andmay also have a 2-thio group and it can read one to four codons depending on its modifications. A G at the wobble position of the anticodon can recognize the two codons ending with U or C and modification of G (unless it is I) does not change its reading properties.     

Is there a special function for U.G basepairs in ribosomal RNA?

U.G basepairs are well-established elements of RNA structure. The geometry of this pair is different, however, from classical Watson-Crick basepairs. This leads to an unusual stacking of the basepair: overlap with the basepair at the 5' side of the U (and the 3' side of the G) is strong (stacked) while it is weak with the basepair on the other side (destacked). The closure of an RNA helix by a U.G pair will be energetically unfavourable when the U residue occupies the 5' end. In transfer RNA there is a strong selection against a 'destacked' U.G pair at helix ends. In the 16S rRNA model of Escherichia coli there are 72 U.G pairs of which 36 or 22 occupy a helix end, depending on how such an end is defined. There is a slight preference for 'stacked' U.G's in these positions. It is remarkable, however, that of 13 very conserved U.G pairs in the 16S (-like) rRNA, 7 occur at helix ends and that 5 of these have the 'destacked' configuration. It is suggested that these pairs, if they exist at all in a hydrogen-bounded form, are stabilized by co-axial stacking with other helices or by interaction with a protein.     

Nucleotides in 16S rRNA that are required in unmodified form for features recognized by ribosomal protein S8.

Nucleotides in 16S rRNA which are required in unmodified form for specific recognition of ribosomal protein S8 from Escherichia coli were identified using a damage-selection experimental approach. Prior to complex formation with S8, 16S rRNA was treated under fully denaturing conditions with either diethyl pyrocarbonate or 25% hydrazine. Following separation of bound from unbound fragments of RNA, those associated with S8 were analyzed for their content of modified bases by treatment with aniline. Nucleotides found to be consistently unmodified in such fragments were located near the base of a stable helix (encompassing bases 581-656) or near the apex of the helix on the 3' proximal side. A minor S8 ribonucleoprotein particle was found to contain fragments which extended in the 3' direction to position 671.     

An NMR structural study of deaminated base pairs in DNA.

The structurally aberrant base pairs TG, UG and TI may occur in DNA as a consequence of deamination of 5-methylcytosine, cytosine and adenine respectively. Results of NMR spectroscopic studies are reported here for these deaminated base pairs in a model seven base pair long oligonucleotide duplex. We find that in all three cases, the DNA helix is a normal B form and both mispaired bases are intrahelical and hydrogen bonded with one another in a wobble geometry. Similarly, in all three cases, all sugars are found to be normal C2' endo in conformation. Symmetric structural perturbations are observed in the helix twist on the 3' side of the mispaired pyrimidine and on the 5' side of the mispaired purine. In all three cases, the amino group of the G residue on the 3' side of the mispaired pyrimidine shows hindered rotation. Although less thermodynamically stable than helices containing only Watson-Crick base pairs, these helices melt normally from the ends and not from the mispair outwards.     

The effects of stem I and loop A on the processing of 5 S rRNA from Drosophila melanogaster.

The 135-nucleotide Drosophila melanogaster 5 S RNA precursor is processed by removal of 15 nucleotides from its 3' end before incorporation into the large ribosomal subunit. Mature 5 S RNA consists of five helical stem-loops; stem IV and part of V are dispensable, whereas stem III and the 1/118 G-C base pair closest to the processing site at nucleotide 120 are required for processing (Preiser, P., and Levinger, L. (1991) J. Biol. Chem. 266, 7509-7516; Preiser, P., and Levinger, L. (1991) J. Biol. Chem. 266, 23602-23605). We have investigated the effects of stem I and loop A transversions, transitions, selected additions and deletions on 5 S RNA processing. Stem I single substitutions generally prevent processing, whereas compensatory double substitutions restore a range of processing rates. Proximal to the processing site, stem I double substitutions inhibit processing. In the distal portion of stem I and loop A, the processing effect of paired sequence changes varies widely in an irregular pattern. The 7/112 GU pair and nucleotide 13A least tolerate sequence changes; several mutations clustered close to the stem I-loop A boundary stimulate processing. We interpret these results in terms of the RNA helix path and possible RNA-protein contacts.     

The influence of the thymine C5 methyl group on spontaneous base pair breathing in DNA

Sequences of four or more AT base pairs without a 5'-TA-3' step, so-called A-tracts, influence the global properties of DNA by causing curvature of the helix axis if phased with the helical repeat and also influence nucleosome packaging. Hence it is interesting to understand this phenomenon on the molecular level, and numerous studies have been devoted to investigations of dynamical and structural features of A-tract DNA. It was early observed that anomalously slow base pair-opening kinetics were a striking physical property unique to DNA A-tracts (Leroy, J. L., Charretier, E., Kochoyan, M., and Gueron, M. (1988) Biochemistry 27, 8894-8898). Furthermore, a strong correlation between DNA curvature and anomalously slow base pair-opening dynamics was found. In the present work it is shown, using imino proton exchange measurements by NMR spectroscopy that the main contribution to the dampening of the base pair-opening fluctuations in A-tracts comes from the C5 methylation of the thymine base. Because the methyl group has been shown to have a very limited effect on the DNA curvature as well as the structure of the DNA helix, the thymine C5 methyl group stabilizes the helix directly. Empirical potential energy calculations show that methylation of the tract improves the stacking energy of a base pair with its neighbors in the tract by 3-4 kcal/mol.     

Crosslinking of transcription factor TFIIIA to ribosomal 5S RNA from X. laevis by trans-diamminedichloroplatinum (II)

Trans-diamminnedichloroplatinum (II) was used to induce reversible crosslinks between 5S rRNA and TFIIIA within the 7S RNP particle from X. laevis immature oocyte. The crosslinked fragments have been unambiguously identified. These fragments exclusively arise from three RNA regions centered around the hinge region at the junction of the three helical domains. Major crosslinking sites are located in region 9-21 (comprising loops A and helix II) and region 54-71 (comprising loop B, helices II and V). A minor site is also found in the 3' part of helix I and helix V (region 100-120). Our results point to the crucial role of the junction region and of the three-dimensional folding of the RNA in the recognition of the 5S rRNA by TFIIIA.     

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.     

Do mRNA and rRNA binding sites of E.coli ribosomal protein S15 share common structural determinants?

Escherichia coli ribosomal protein S15 recognizes two RNA targets: a three-way junction in 16S rRNA and a pseudoknot structure on its own mRNA. Binding to mRNA occurs when S15 is expressed in excess over its rRNA target, resulting in an inhibition of translation start. The sole apparent similarity between the rRNA and mRNA targets is the presence of a G-U/G-C motif that contributes only modestly to rRNA binding but is essential for mRNA. To get more information on the structural determinants used by S15 to bind its mRNA target as compared to its rRNA site, we used site-directed mutagenesis, substitution by nucleotide analogs, footprinting experiments on both RNA and protein, and graphic modeling. The size of the mRNA-binding site could be reduced to 45 nucleotides, without loss of affinity. This short RNA preferentially folds into a pseudoknot, the formation of which depends on magnesium concentration and temperature. The size of the loop L2 that bridges the two stems of the pseudoknot through the minor groove could not be reduced below nine nucleotides. Then we showed that the pseudoknot recognizes the same side of S15 as 16S rRNA, although shielding a smaller surface area. It turned out that the G-U/G-C motif is recognized from the minor groove in both cases, and that the G-C pair is recognized in a very similar manner. However, the wobble G-U pair of the mRNA is not directly contacted by S15, as in rRNA, but is most likely involved in building a precise conformation of the RNA, essential for binding. Otherwise, unique specific features are utilized, such as the three-way junction in the case of 16S rRNA and the looped out A(-46) for the mRNA pseudoknot.     

Different conformational forms of Escherichia coli and rat liver 5S rRNA revealed by Pb(II)-induced hydrolysis.

Different stable forms of Escherichia coli and rat liver 5S rRNA have been probed by Pb(II)-induced hydrolysis. In the native A forms of 5S rRNA, Pb2+ reveal single-stranded RNA stretches and regions of increased conformational flexibility or distorted by the presence of bulged nucleotides. Hydrolysis of urea/EDTA-treated E. coli 5S rRNA (B form) shows the presence of two strong helical domains; helix A retained from the A form and a helix composed of RNA regions G33-C42 and G79-C88. Other RNA regions resistant to hydrolysis may be involved in alternative base pairing, causing conformational heterogeneity of that form. Pb(II)-induced hydrolysis distinguishes two different forms of rat liver 5S rRNA; the native A form and the form obtained by renaturation of 5S rRNA in the presence of EDTA. Pb(II)-hydrolysis data suggest that both forms are highly structured. In the latter form, the orientation of the bulged C66 is changed with respect to helix B. At the same time, a new helical segment is possibly formed, composed of nucleotides from helix C and loop c on one side and from helix E and loop d' on the other.     

Assignment of phosphorus-31 and nonexchangeable proton resonances in a symmetrical 14 base pair lac pseudooperator DNA fragment

The 31P chemical shifts of all 13 phosphates and the chemical shifts of nearly all of the non-exchangeable protons of a symmetrical 14 base pair lac pseudooperator DNA fragment have been assigned by regiospecific labeling with oxygen-17 and two-dimensional NMR techniques. At 22 degrees C, 8 of the 13 phosphorus resonances can distinctly be resolved while the remaining 5 resonances occur in two separate overlapping regions. The 31P chemical shifts of this particular 14 base pair oligonucleotide do not follow the general observation that the more internal the phosphate is located within the oligonucleotide sequence the more upfield the 31P resonance occurs, as shown from other 31P assignment studies. Failure of this general rule is believed to be a result of helical distortions that occur along the oligonucleotide double helix, on the basis of the analysis of Callidine [Callidine, C.R. (1982) J. Mol. Biol. 161, 343-352]. Notable exceptions to the phosphate position relationship are 5'-Py-Pu-3' dinucleotide sequences, which resonate at a lower field strength than expected in agreement with similar results as reported by Ott and Eckstein [Ott, J., & Eckstein, F. (1985) Biochemistry 24, 253]. A reasonable correlation exists between 31P chemical shifts values of the 14-mer and the helical twist sum function of Calladine. The most unusual 31P resonance occurs most upfield in the 31P spectrum, which has been assigned to the second phosphate position (5'-GpT-3') from the 5' end. This unusual chemical shift may be the result of the predicted large helical twist angle that occurs at this position in the 14-mer sequence. Further, it is believed that the large helical twist represents a unique structural feature responsible for optimum binding contact between lac repressor protein and this 14-mer lac pseudooperator segment. Assignments of proton resonances were made from two-dimensional 1H-1H nuclear Overhauser effect (NOESY) connectivities in a sequential manner applicable to right-handed B-DNA, in conjunction with two-dimensional homonuclear and heteronuclear J-correlated spectroscopies (1H-1H COSY and 31P-1H HETCOR). Most nonexchangeable base proton and deoxyribose proton (except for some unresolved H4', H5', and H5" protons) resonances were assigned.     

A second base pair interaction between U3 small nucleolar RNA and the 5'-ETS region is required for early cleavage of the yeast pre-ribosomal RNA

In eukaryotes, U3 snoRNA is essential for pre-rRNA maturation. Its 5'-domain was found to form base pair interactions with the 18S and 5'-ETS parts of the pre-rRNA. In Xenopus laevis, two segments of U3 snoRNA form base-pair interactions with the 5'-ETS region and only one of them is essential to the maturation process. In Saccharomyces cerevisiae, two similar U3 snoRNA-5' ETS interactions are possible; but, the functional importance of only one of them had been tested. Surprisingly, this interaction, which corresponds to the non-essential one in X. laevis, is essential for cell growth and pre-rRNA maturation in yeast. In parallel with [Dutca et al. (2011) The initial U3 snoRNA:pre-rRNA base pairing interaction required for pre-18S rRNA folding revealed by in vivo chemical probing. Nucleic Acids Research, 39, 5164-5180], here we show, that the second possible 11-bp long interaction between the 5' domain of S. cerevisiae U3 snoRNA and the pre-rRNA 5'-ETS region (helix VI) is also essential for pre-rRNA processing and cell growth. Compensatory mutations in one-half of helix VI fully restored cell growth. Only a partial restoration of growth was obtained upon extension of compensatory mutations to the entire helix VI, suggesting sequence requirement for binding of specific proteins. Accordingly, we got strong evidences for a role of segment VI in the association of proteins Mpp10, Imp4 and Imp3.     

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.     

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.     

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.     

Nearest neighbor parameters for inosine x uridine pairs in RNA duplexes

An enzyme family known as adenosine deaminases that act on RNA (ADARs) catalyzes adenosine deamination in RNA. ADARs act on RNA that is largely double-stranded and convert adenosine to inosine, resulting, in many cases, in an I x U pair. Thermodynamic parameters derived from optical melting studies are reported for a series of 14 oligoribonucleotides containing single I x U pairs adjacent to Watson-Crick pairs. In order to determine unique linearly independent nearest neighbor parameters for I x U pairs, four duplexes containing 3'-terminal I x U pairs and four duplexes containing 5'-terminal I x U pairs have also been thermodynamically characterized. This data was combined with previously published data of seven duplexes containing internal, terminal, or tandem I x U pairs from Strobel et al. [Strobel, S. A., Cech, T. R., Usman, N., and Beigelman, L. (1994) Biochemistry 33, 13824-13838] and Serra et al. [Serra, M. J., Smolter, P. E., and Westhof, E. (2004) Nucleic Acids Res. 32, 1824-1828]. On average, a duplex with an internal I x U pair is 2.3 kcal/mol less stable than the same duplex with an A-U pair, however, a duplex with a terminal I x U pair is 0.8 kcal/mol more stable than the same duplex with an A-U pair. Although isosteric with a G-U pair, on average, a duplex with an internal I x U pair is 1.9 kcal/mol less stable than the same duplex with a G-U pair, however, a duplex with a terminal I x U pair is 0.9 kcal/mol more stable than the same duplex with a G-U pair. Duplexes with tandem I x U pairs are on average 5.9 and 3.8 kcal/mol less stable than the same duplex with tandem A-U or tandem G-U pairs, respectively. Using the combined thermodynamic data and a complete linear least-squares fitting routine, nearest neighbor parameters for all nearest neighbor combinations of I x U pairs and an additional parameter for terminal I x U pairs have been derived.     

Solution structure of an A-tract DNA bend

The solution structure of a DNA dodecamer d(GGCAAAAAACGG)/d(CCGTTTTTTGCC) containing an A-tract has been determined by NMR spectroscopy with residual dipolar couplings. The structure shows an overall helix axis bend of 19 degrees in a geometry consistent with solution and gel electrophoresis experiments. Fourteen degrees of the bending occurs in the GC regions flanking the A-tract. The remaining 5 degrees is spread evenly over its six AT base-pairs. The A-tract is characterized by decreasing minor groove width from the 5' to the 3' direction along the A strand. This is a result of propeller twist in the AT pairs and the increasing negative inclination of the adenine bases at the 3' side of the run of adenine bases. The four central thymine bases all have negative inclination throughout the A-tract with an average value of -6.1 degrees. Although this negative inclination makes the geometry of the A-tract different from all X-ray structures, the proton on N6 of adenine and the O4 of thymine one step down the helix are within distance to form bifurcated hydrogen bonds. The 5' bend of 4 degrees occurs at the junction between the GC flank and the A-tract through a combination of tilt and roll. The larger 3' bend, 10 degrees, occurs in two base steps: the first composed of tilt, -4.1 degrees, and the second a combination of tilt, -4.2 degrees, and roll, 6.0 degrees. This second step is a direct consequence of the change in inclination between an adjacent cytosine base, which has an inclination of -12 degrees, and the next base, a guanine, which has 3 degrees inclination. This bend is a combination of tilt and roll. The large change in inclination allows the formation of a hydrogen bond between the protons of N4 of the 3' cytosine and the O6 of the next 3' base, a guanine, stabilizing the roll component in the bend. These structural features differ from existing models for A-tract bends.For comparison, we also determined the structure of the control sequence, d(GGCAAGAAACGG)/d(CCGTTTCTTGCC), with an AT to GC transition in the center of the A-tract. This structure has no negative inclination in most of the bases within the A-tract, resulting in a bend of only 9 degrees. When ligated in phase, the control sequence has nearly normal mobility in gel electrophoresis experiments.     

A small angle x-ray scattering study of a fragment derived from E. coli 5S RNA.

The structure of a 62 base nuclease resistant fragment of E. coli 5S RNA (bases 1-11, 69-87, 89-120) has been examined by small angle x-ray scattering. The results obtained are indistinguishable from those expected if this oligonucleotide complex were a perfect RNA double helix of about 30 base pairs. These results indicate that this portion of 5S RNA is in a configuration which is approximately double helical, even though proper base pairing is possible over only half its length.