Conserved unpaired adenine residues are important for ordered structures of 5S ribosomal RNA. An infrared study of the secondary and tertiary structure of Thermus thermophilus 5S rRNA

An improved set of infrared calibration spectra for the determination of G X C and A X U base pairs leads to 32 +/- 3 G X C (+ G X U) and 4 +/- 1 A X U base pairs for Thermus thermophilus 5S RNA in the presence and absence of Mg2+. These results give further support for the consensus secondary structure of 5S RNA recently proposed by several groups. T. thermophilus 5S RNA shows, in the presence of Mg2+, a distinct two-step thermal melting of its ordered structure. Based on new data about the stacking dependence of infrared intensities of unpaired ribonucleotides the spectral changes of the low-temperature transition should be explained by melting of stacked arrangements of unpaired bases and/or non-standard base pairs. Striking is the reduction in A stacking, which is not related to the melting of A X U base pairs, indicating the importance of the mostly conserved unpaired adenines for the Mg2+ stabilized higher-order structures especially within internal loops of 5S RNA.     

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.     

Template-Independent Repair of the 3′ End of Cucumber Mosaic Virus Satellite RNA Controlled by RNAs 1 and 2 of Helper Virus

RNA viruses which do not have a poly(A) tail or a tRNA-like structure for the protection of their vulnerable 3′ termini may have developed a different strategy to maintain their genome integrity. We provide evidence that deletions of up to 7 nucleotides from the 3′ terminus of cucumber mosaic cucumovirus (CMV) satellite RNA (satRNA) were repaired in planta in the presence of the helper virus (HV) CMV. Sequence comparison of 3′-end-repaired satRNA progenies, and of satRNA and HV RNA, suggested that the repair was not dependent on a viral template. The 3′ end of CMV satRNA lacking the last three cytosines was not repaired in planta in the presence of tomato aspermy cucumovirus (TAV), although TAV is an efficient helper for the replication of CMV satRNA. With use of pseudorecombinants constructed by the interchange of RNAs 1 and 2 of TAV and CMV, evidence was provided that the 3′-end repair was controlled by RNAs 1 and 2 of CMV, which encode subunits of the viral RNA replicase. These results, and the observation of short repeated sequences close to the 3′ terminus of repaired molecules, suggest that the HV replicase maintains the integrity of the satRNA genome, playing a role analogous to that of cellular telomerases.     

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.     

A method for estimating the nearest neighbor base-pair content of RNAs using CD and absorption spectroscopy.

CD and absorption spectra are sensitive to the secondary structure of RNAs. By fitting the spectra contained in our basis set to the CD and absorption spectra of an RNA of known sequence, we could determine the fractions of base pairs, the fractions of each of the nearest neighbor base pairs, and the fractions of the single-stranded nucleotides in that RNA. The basis set included 58 CD and 58 absorption spectra. The fitting was done with a guided selection routine. The estimated error was about 0.05 for predicting the fractions of the nearest neighbor base pairs, 0.06 for predicting the fractions of A.U, G.C, and G.U base pairs, and 0.04 for predicting the fractions of the single-stranded nucleotides.     

A new pathway expressed during a distinct stage of Drosophila development for the removal of dUMP residues in DNA

In view of removing lesions in DNA produced by the deamination of cytosine to uracil, uracil-DNA glycosylases were anticipated to be ubiquitous. However, an analogous activity in Drosophila melanogaster was not detected. Instead, a nuclease was identified that acts specifically upon DNA containing uracil. The cleavage of uracil-containing DNA by the nuclease generates acid-soluble oligonucleotides in a reaction which can be inhibited by pretreatment of the DNA with Escherichia coli uracil-DNA glycosylase. Uracil-containing DNA with either A:U base pairs or G:U base pairs were susceptible to cleavage by the nuclease, whereas other damaged DNA substrates were not. The nuclease activity is transient and appears only in third instar larvae, with other developmental stages of Drosophila lacking significant levels of the nuclease.     

Isoalloxazine derivatives promote photocleavage of natural RNAs at G.U base pairs embedded within helices.

We have recently shown that isoalloxazine derivatives are able to photocleave RNA specifically at G.U base pairs embedded within a helical stack. The reaction involves the selective molecular recognition of G.U base pairs by the isoalloxazine ring and the removal of one nucleoside downstream of the uracil residue. Divalent metal ions are absolutely required for cleavage. Here we extend our studies to complex natural RNA molecules with known secondary and tertiary structures, such as tRNAs and a group I intron (td). G.U pairs were cleaved in accordance with the phylogenetically and experimentally derived secondary and tertiary structures. Tandem G.U pairs or certain G.U pairs located at a helix extremity were not affected. These new cleavage data, together with the RNA crystal structure, allowed us to perform molecular dynamics simulations to provide a structural basis for the observed specificity. We present a stable structural model for the ternary complex of the G. U-containing helical stack, the isoalloxazine molecule and a metal ion. This model provides significant new insight into several aspects of the cleavage phenomenon, mechanism and specificity for G. U pairs. Our study shows that in large natural RNAs a secondary structure motif made of an unusual base pair can be recognized and cleaved with high specificity by a low molecular weight molecule. This photocleavage reaction thus opens up the possibility of probing the accessibility of G.U base pairs, which are endowed with specific structural and functional roles in numerous structured and catalytic RNAs and interactions of RNA with proteins, in folded RNAs.     

The G x U wobble base pair. A fundamental building block of RNA structure crucial to RNA function in diverse biological systems

The G x U wobble base pair is a fundamental unit of RNA secondary structure that is present in nearly every class of RNA from organisms of all three phylogenetic domains. It has comparable thermodynamic stability to Watson-Crick base pairs and is nearly isomorphic to them. Therefore, it often substitutes for G x C or A x U base pairs. The G x U wobble base pair also has unique chemical, structural, dynamic and ligand-binding properties, which can only be partially mimicked by Watson-Crick base pairs or other mispairs. These features mark sites containing G x U pairs for recognition by proteins and other RNAs and allow the wobble pair to play essential functional roles in a remarkably wide range of biological processes.     

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.     

Higher order structure of chloroplastic 5S ribosomal RNA from spinach

The secondary and tertiary structure of chloroplastic 5S ribosomal RNA from spinach was investigated by the use of several chemical and enzymatic structure probes. The four bases were monitored at one of their Watson-Crick base-pairing positions with dimethyl sulfate [at A(N1) and C(N3)] and with 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate [at G(N1) and U(N3)]. Position N7 of purines was probed with diethyl pyrocarbonate (adenines) and with dimethyl sulfate (guanines). Ethylnitrosourea was used to probe phosphate involved in tertiary interaction or in cation coordination. In order to estimate the degree of stability of helices, the various chemical reagents were employed under "native" conditions (300 mM KCl and 20 mM magnesium at 37 degrees C), under "semidenaturing" conditions [1 mM ethylenediaminetetraacetic acid (EDTA) at 37 degrees C], and under denaturing conditions (1 mM EDTA at 90 degrees C). Unstructured regions were also tested with single-strand-specific nucleases T1, U2, and S1 and double-stranded or stacked regions with RNase V1 from cobra Naja naja oxiana venom. The results confirm the existence of the five helices and the two external loops proposed in the consensus model of 5S rRNA. However, the regions depicted as unpaired internal loops appear to be folded into a more complex conformation. A three-dimensional model derived from the present data and graphic modeling for a region encompassing helix IV, helix V, loop D, and loop E (nucleotides 70-110) is proposed. Nucleotides in the so-called loop E (73-79/100-106) display unusual features: Noncanonical base pairs (A-A and A-G) are formed, and three nucleotides (C75, U78, and U105) are bulging out. This region adopts an unwound and extended conformation that can be well suited for tertiary interactions or for protein binding. Several bases and phosphates candidate for the tertiary folding of the RNA were also identified.     

Identification of three G.U base pairs in Bacillus subtilis ribosomal 5S RNA via 500-MHz proton homonuclear Overhauser enhancements

Three distinct G.U base pairs in Bacillus subtilis 5S RNA have been identified via homonuclear Overhauser enhancements (NOE) of their low-field (9-15 ppm) proton Fourier transform nuclear magnetic resonances at 11.75 T. With these G.U resonances as starting points, short segments of NOE connectivity can be established. One G.U-G.C-G.C segment (most probably G4.C112-G5.C111-U6.G110) can definitely be assigned to the terminal helix. The existence of at least part of the terminal helical stem of the secondary structure of a Gram-positive bacterial 5S RNA has thus been established for the first time by direct experimental observation. Addition of Mg2+ produces almost no conformational changes in the terminal stem but results in major conformational changes elsewhere in the structure, as reflected by changes in the 1H 500-MHz low-field NMR spectrum. Assignment of the two remaining G.U base pairs will require further experiments (e.g., enzymatic-cleavage fragments). Finally, the implications of these results for analysis of RNA secondary structure are discussed.     

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.     

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.     

Satellite RNA of cucumber mosaic virus forms a secondary structure with partial 3''-terminal homology to genomal RNAs.

Sat-RNA is one of several replicating satellite RNAs which have been isolated from RNA encapsidated in cucumber mosaic virus (CMV) and which are totally dependent on CMV for replication. The 336 residue sequence of Sat-RNA obtained using the dideoxynucleotide chain termination and partial enzymic digestion procedures shows only a few short stretches (up to 11 residues) of sequence homology with one of the three CMV genomal RNAs so far sequenced. Sat-RNA has 88% sequence homology with another, previously sequenced, satellite RNA of CMV, CARNA 5. Analysis of partial digests of 5'- or 3' -32P-Sat-RNA with nuclease S1 or RNase T1 under non-denaturing conditions showed that only about 10% of the residues in Sat-RNA were cleaved. Further data on base-paired segments of Sat-RNA were obtained using digestion with RNase T1 followed by electrophoretic fractionation of the resulting fragments under both non-denaturing and denaturing conditions. On the basis of this data, a complete secondary structure model is proposed for Sat-RNA with 52% of its residues involved in base pairs. A prominent hairpin at the 3'-terminus of Sat-RNA shows considerable sequence and structural homology with parts of the 3'-terminal tRNA-like structure of the CMV genomal RNAs.     

Computational analysis of hydrogen bonds in protein-RNA complexes for interaction patterns

Structural analysis of protein-RNA complexes is labor-intensive, yet provides insight into the interaction patterns between a protein and RNA. As the number of protein-RNA complex structures reported has increased substantially in the last few years, a systematic method is required for automatically identifying interaction patterns. This paper presents a computational analysis of the hydrogen bonds in the most representative set of protein-RNA complexes. The analysis revealed several interesting interaction patterns. (1) While residues in the beta-sheets favored unpaired nucleotides, residues in the helices showed no preference and residues in turns favored paired nucleotides. (2) The backbone hydrogen bonds were more dominant than the base hydrogen bonds in the paired nucleotides, but the reverse was observed in the unpaired nucleotides. (3) The protein-RNA complexes contained more paired nucleotides than unpaired nucleotides, but the unpaired nucleotides were observed more frequently interacting with the proteins. And (4) Arg-U, Thr-A, Lys-A, and Asn-U were the most frequently observed pairs. The interaction patterns discovered from the analysis will provide us with useful information in predicting the structure of the RNA binding protein and the structure of the protein binding RNA.     

Sequence dependence for the energetics of dangling ends and terminal base pairs in ribonucleic acid

Stability increments of terminal unpaired nucleotides (dangling ends) and terminal base pairs on the core helixes AUGCAU and UGCGCA are reported. Enthalpy, entropy, and free energy changes of helix formation were measured spectrophotometrically for 18 oligoribonucleotides containing the core sequences. The results indicate 3' dangling purines add more stability than 3' dangling pyrimidines. In most cases, the additional stability from a 3' dangling end on an AU base pair is less than that on a GC base pair [Freier, S.M., Burger, B.J., Alkema, D., Neilson, T., & Turner, D.H. (1985) Biochemistry 22, 6198-6206]. The sequence dependence provides a test for the importance of dangling ends for various RNA interactions. Correlations are suggested with codon context effects and with the three-dimensional structure of yeast phenylalanine transfer RNA. In the latter case, all terminal unpaired nucleotides having stability increments more favorable than -1 kcal/mol are stacked on the adjacent base pair. All terminal unpaired nucleotides having stability increments less favorable than -0.3 kcal/mol are not stacked on the adjacent base pair. In several cases, this lack of stacking is associated with a turn in the sugar-phosphate backbone. This suggests stability increments measured on oligoribonucleotides may be useful for predicting tertiary structure in large RNA molecules. Comparison of the stability increments for terminal dangling ends and base pairs, and of terminal GC and AU base pairs, indicates the free energy increment associated with forming a hydrogen bond can be about -1 kcal/mol of hydrogen bond.     

Divalent metal ion binding to a conserved wobble pair defining the upstream site of cleavage of group I self-splicing introns.

The upstream site of cleavage of all group I self-splicing introns is identified by an absolutely conserved U.G base pair. Although a wobble C.A pair can substitute the U.G pair, all other combinations of nucleotides at this position abolish splicing, suggesting that it is an unusual RNA structure, rather than sequence, that is recognized by the catalytic intron core. RNA enzymes are metalloenzymes, and divalent metal ion binding may be an important requirement for splice site recognition and catalysis. The paramagnetic broadening of NMR resonances upon manganese binding at specific sites was used to probe the interaction between divalent metal ions and an oligonucleotide model of a group I intron ribozyme substrate. Unlike previous studies in which only imino proton resonances were monitored, we have used isotopically labelled RNA and a set of complete spectral assignments to identify the location of the divalent metal binding site with much greater detail than previously possible. Two independent metal binding sites were identified for this oligonucleotide. A first metal binding site is located in the major groove of the three consecutive G.C base pairs at the end of double helical stem. A second site is found in the major groove of the RNA double helix in the vicinity of the U.G base pair. These results suggest that metal ion coordination (or a metal bridge) and tertiary interactions identified biochemically, may be used by group I intron ribozymes for substrate recognition.