The topography of the 5' end of 16-S RNA in the presence and absence of ribosomal proteins S4 and S20

A ribonucleoprotein prepared by strong ribonuclease digestion of a complex of 16-S ribosomal RNA and proteins S4 and S20 from Escherichia coli has been characterized; its nucleotide sequence, the positions of enzyme cuts and the sequence excisions have been placed in the completed sequence of 16-S RNA. The positions and yields of enzyme cuts, and excisions of sequence, are compared with those of various ribonucleoproteins prepared with S4 or S20 alone, and with the ribonuclease-resistant S4 RNA prepared from renatured 16-s RNA in the absence of ribosomal protein. These data yield important information on the topography and organisation of the 5' third of the 16-s RNA which is selectively maintained in its native conformation by the bound proteins; they also provide criteria for testing secondary structural models of this region of 16-S RNA.     

Comparative photocross-linking analysis of the tertiary structures of Escherichia coli and Bacillus subtilis RNase P RNAs.

Bacterial ribonuclease P contains a catalytic RNA subunit that cleaves precursor sequences from the 5' ends of pre-tRNAs. The RNase P RNAs from Bacillus subtilis and Escherichia coli each contain several unique secondary structural elements not present in the other. To understand better how these phylogenetically variable elements affect the global architecture of the ribozyme, photoaffinity cross-linking studies were carried out. Photolysis of photoagents attached at homologous sites in the two RNAs results in nearly identical cross-linking patterns, consistent with the homology of the RNAs and indicating that these RNAs contain a common, core tertiary structure. Distance constraints were used to derive tertiary structure models using a molecular mechanics-based modeling protocol. The resulting superimposition of large sets of equivalent models provides a low resolution (5-10 A) structure for each RNA. Comparison of these structure models shows that the conserved core helices occupy similar positions in space. Variably present helical elements that may play a role in global structural stability are found at the periphery of the core structure. The P5.1 and P15.1 helical elements, unique to the B.subtilis RNase P RNA, and the P6/16/17 helices, unique to the E.coli RNA, occupy similar positions in the structure models and, therefore, may have analogous structural function.     

Hydroxyl radical footprinting of ribosomal proteins on 16S rRNA.

Complexes between 16S rRNA and purified ribosomal proteins, either singly or in combination, were assembled in vitro and probed with hydroxyl radicals generated from free Fe(II)-EDTA. The broad specificity of hydroxyl radicals for attack at the ribose moiety in both single- and double-stranded contexts permitted probing of nearly all of the nucleotides in the 16S rRNA chain. Specific protection of localized regions of the RNA was observed in response to assembly of most of the ribosomal proteins. The locations of the protected regions were in good general agreement with the footprints previously reported for base-specific chemical probes, and with sites of RNA-protein crosslinking. New information was obtained about interaction of ribosomal proteins with 16S rRNA, especially with helical elements of the RNA. In some cases, 5' or 3' stagger in the protection pattern on complementary strands suggests interaction of proteins with the major or minor groove, respectively, of the RNA. These results reinforce and extend previous data on the localization of ribosomal proteins with respect to structural features of 16S rRNA, and offer many new constraints for three-dimensional modeling of the 30S ribosomal subunit.     

Dependence among sites in RNA evolution

Although probabilistic models of genotype (e.g., DNA sequence) evolution have been greatly elaborated, less attention has been paid to the effect of phenotype on the evolution of the genotype. Here we propose an evolutionary model and a Bayesian inference procedure that are aimed at filling this gap. In the model, RNA secondary structure links genotype and phenotype by treating the approximate free energy of a sequence folded into a secondary structure as a surrogate for fitness. The underlying idea is that a nucleotide substitution resulting in a more stable secondary structure should have a higher rate than a substitution that yields a less stable secondary structure. This free energy approach incorporates evolutionary dependencies among sequence positions beyond those that are reflected simply by jointly modeling change at paired positions in an RNA helix. Although there is not a formal requirement with this approach that secondary structure be known and nearly invariant over evolutionary time, computational considerations make these assumptions attractive and they have been adopted in a software program that permits statistical analysis of multiple homologous sequences that are related via a known phylogenetic tree topology. Analyses of 5S ribosomal RNA sequences are presented to illustrate and quantify the strong impact that RNA secondary structure has on substitution rates. Analyses on simulated sequences show that the new inference procedure has reasonable statistical properties. Potential applications of this procedure, including improved ancestral sequence inference and location of functionally interesting sites, are discussed.     

Modeling large RNAs and ribonucleoprotein particles using molecular mechanics techniques.

There is a growing body of low-resolution structural data that can be utilized to devise structural models for large RNAs and ribonucleoproteins. These models are routinely built manually. We introduce an automated refinement protocol to utilize such data for building low-resolution three-dimensional models using the tools of molecular mechanics. In addition to specifying the positions of each nucleotide, the protocol provides quantitative estimates of the uncertainties in those positions, i.e., the resolution of the model. In typical applications, the resolution of the models is about 10-20 A. Our method uses reduced representations and allows us to refine three-dimensional structures of systems as big as the 16S and 23S ribosomal RNAs, which are about one to two orders of magnitude larger than nucleic acids that can be examined by traditional all-atom modeling methods. Nonatomic resolution structural data--secondary structure, chemical cross-links, chemical and enzymatic footprinting patterns, protein positions, solvent accessibility, and so on--are combined with known motifs in RNA structure to predict low-resolution models of large RNAs. These structural constraints are imposed on the RNA chain using molecular mechanics-type potential functions with parameters based on the quality of experimental data. Surface potential functions are used to incorporate shape and positional data from electron microscopy image reconstruction experiments into our models. The structures are optimized using techniques of energy refinement to get RNA folding patterns. In addition to providing a consensus model, the method finds the range of models consistent with the data, which allows quantitative evaluation of the resolution of the model. The method also identifies conflicts in the experimental data. Although our protocol is aimed at much larger RNAs, we illustrate these techniques using the tRNA structure as an example and test-bed.     

Close packing of helices 3 and 12 of 16 S rRNA is required for the normal ribosome function

The along-groove packing motif is a quasi-reciprocal arrangement of two RNA double helices in which a backbone of each helix is closely packed within the minor groove of the other helix. At the center of the inter-helix contact, a GU base pair in one helix packs against a Watson-Crick base pair in the other helix. Here, based on in vivo selection from a combinatorial gene library of 16 S rRNA and on functional characterization of the selected clones, we demonstrate that the normal ribosome performance requires that helices 3 and 12 be closely packed. In some clones the Watson-Crick and GU base pairs exchange in their positions between the two helices, which affects neither the quality of the helix packing, nor the ribosome function. On the other hand, perturbations in the close packing usually lead to a substantial drop in the ribosome activity. The functionality of the clones containing such perturbations may depend on the presence of particular elements in the vicinity of the area of contact between helices 3 and 12. Such cases do not exist in natural 16 S rRNA, and their selection enriches our knowledge of the constraints imposed on the structure of ribosomal RNA in functional ribosomes.     

Xenopus laevis 18S ribosomal RNA: experimental determination of secondary structural elements, and locations of methyl groups in the secondary structure model.

18S ribosomal RNA from X. laevis was subjected to partial digestion with ribonucleases A or T1 under a variety of conditions, and base-paired fragments were isolated. Sequence analysis of the fragments enabled five base-paired secondary structural elements of the 18S RNA to be established. Four of these elements (covering bases 221-256, 713-757, 1494-1555 and 1669-1779) confirm our previous secondary structure predictions, whereas the fifth (comprising bases 1103-1125) represents a phylogenetically conserved "switch" structure, which can also form in prokaryotic 16S RNA. The results are incorporated into a refined model of the 18S RNA secondary structure, which also includes the locations of the many methyl groups in X. laevis 18S RNA. In general the methyl groups occur in non-helical regions, at hairpin loop ends, or at helix boundaries and imperfections. One large cluster of 2'-O-methyl groups occurs in a region of complicated secondary structure in the 5'-one third of the molecule.     

Computer modeling 16 S ribosomal RNA

A three-dimensional structure for 16 S RNA has been produced with a computer protocol that is not dependent on human intervention. This protocol improves upon traditional modeling techniques by using distance geometry to fold the molecule in an objective and reproducible fashion. The method is based on the secondary structure of RNA and treats the molecule as a set of double-stranded helices that are linked by flexible single-strands of variable length. Data derived from chemical cross-linking studies of 16 S RNA and tertiary phylogenetic relationships provide the constraints used to fold the molecule into a compact three-dimensional form. Possibly subjective evaluation of the input data are transformed into verifiable quantitative parameters. Relationships based on general locations within the 30 S subunit or on protein-RNA interactions have been specifically excluded. The resolution of the model exceeds that of electron micrographs and approaches that obtained in preliminary X-ray crystal structures. The model size of 245 x 190 x 140 A is compatible with that of the 30 S subunit as determined by electron microscopy. The volume of the model is 1.87 x 10(6) A which is similar to that of the small subunit in a preliminary X-ray crystal structure. The radius of gyration of the model structure of 76 A is intermediate to that seen for partially denatured and fully folded 16 S RNA. Computer graphics are used to display the results in a manner that maximizes the opportunities for human visual interpretation of the models. A format for displaying the structures has been developed that will make it possible for researchers who have not devoted themselves to ribosomal modeling to comprehend and make use of the information that the models embody. On this basis the computer-generated models are compared with models developed by other researchers and with structural data not included in the folding parameter data set.     

Internal transcribed spacer 1 of the yeast precursor ribosomal RNA. Higher order structure and common structural motifs

The higher order structure of the first internal transcribed spacer between the 18S and the 5.8S rRNA sequences in the Saccharomyces cerevisiae precursor ribosomal RNA has been investigated. Sites of potential base pairing in the RNA region have been determined by using a combination of enzymatic and chemical structure sensitive probes. Data generated have been used to evaluate secondary structure models predicted by minimum free energy calculations. Several alternative suboptimal structures were also evaluated. The derived model contains several stable hairpins. Theoretical secondary structural models for the corresponding RNA region from S. carlsbergensis, S. pombe, N. crassa, X. laevis, and mung bean have also been derived from identical calculations and assumptions. Certain structural motifs appear to be conserved despite extensive divergence in the base sequence. The yeast model should be a useful prototype for investigation of structure and function of precursor ribosomal RNA molecules.     

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.     

The 5S ribosomal RNA of Euglena gracilis cytoplasmic ribosomes is closely homologous to the 5S RNA of the trypanosomatid protozoa.

The complete nucleotide sequence of the major species of cytoplasmic 5S ribosomal RNA of Euglena gracilis has been determined. The sequence is: 5' GGCGUACGGCCAUACUACCGGGAAUACACCUGAACCCGUUCGAUUUCAGAAGUUAAGCCUGGUCAGGCCCAGUUAGUAC UGAGGUGGGCGACCACUUGGGAACACUGGGUGCUGUACGCUUOH3'. This sequence can be fitted to the secondary structural models recently proposed for eukaryotic 5S ribosomal RNAs (1,2). Several properties of the Euglena 5S RNA reveal a close phylogenetic relationship between this organism and the protozoa. Large stretches of nucleotide sequences in predominantly single-stranded regions of the RNA are homologous to that of the trypanosomatid protozoan Crithidia fasticulata. There is less homology when compared to the RNAs of the green alga Chlorella or to the RNAs of the higher plants. The sequence AGAAC near position 40 that is common to plant 5S RNAs is CGAUU in both Euglena and Crithidia. The Euglena 5S RNA has secondary structural features at positions 79-99 similar to that of the protozoa and different from that of the plants. The conclusions drawn from comparative studies of cytochrome c structures which indicate a close phylogenetic relatedness between Euglena and the trypanosomatid protozoa are supported by the comparative data with 5S ribosomal RNAs.     

Investigation of the tertiary folding of Escherichia coli 16S RNA by in situ intra-RNA cross-linking within 30S ribosomal subunits.

Intra-RNA cross-links were introduced into E. coli 30S ribosomal subunits by mild ultraviolet irradiation. The subunits were partially digested with cobra venom nuclease, followed in some cases by a second partial digestion with ribonuclease H in the presence of the hexanucleotide d-(CTTCCC). The cross-linked RNA complexes were separated by two-dimensional gel electrophoresis and the sites of cross-linking analysed by our published procedures. Tertiary structural cross-links in the 16S RNA were identified between positions 31 and 48, between oligonucleotides 1090-1094 and 1161-1164, and between oligonucleotides 1125-1127 and 1280-1281. The first of these imposes a rigid constraint on the relative orientations of helices 3 and 4 of the 16S secondary structure. A further tertiary cross-link (which could not be precisely localised) was found between regions 1-72 and 1020-1095, and secondary structural cross-links were identified between positions 497 and 545-548, and positions 1238-1240 and 1298.     

The 16S rRNA binding site of Thermus thermophilus ribosomal protein S15: comparison with Escherichia coli S15, minimum site and structure.

Binding of Escherichia coli and Thermus thermophilus ribosomal proteins S15 to a 16S ribosomal RNA fragment from T. thermophilus (nt 559-753) has been investigated in detail by extensive deletion analysis, filter-binding assays, gel mobility shift, structure probing, footprinting with chemical, enzymatic, and hydroxyl radical probes. Both S15 proteins recognize two distinct sites. The first one maps in the bottom of helix 638-655/717-734 (H22) and in the three-way junction between helix 560-570/737-747 (H20), helix 571-600/606-634 (H21), and H22. The second is located in a conserved purine-rich region in the center of H22. The first site provides a higher contribution to the free energy of binding than the second one, and both are required for efficient binding. A short RNA fragment of 56 nt containing these elements binds S15 with high affinity. The structure of the rRNA is constrained by the three-way junction and requires both magnesium and S15 to be stabilized. A 3D model, derived by computer modeling with the use of experimental data, suggests that the bound form adopts a Y-shaped conformation, with a quasi-coaxial stacking of H22 on H20, and H21 forming an acute angle with H22. In this model, S15 binds to the shallow groove of the RNA on the exterior side of the Y-shaped structure, making contact with the two sites, which are separated by one helix turn.     

A detailed model of the three-dimensional structure of Escherichia coli 16 S ribosomal RNA in situ in the 30 S subunit

A large body of intra-RNA and RNA-protein crosslinking data, obtained in this laboratory, was used to fold the phylogenetically and experimentally established secondary structure of Escherichia coli 16 S RNA into a three-dimensional model. All the crosslinks were induced in intact 30 S subunits (or in some cases in growing E. coli cells), and the sites of crosslinking were precisely localized on the RNA by oligonucleotide analysis. The RNA-protein crosslinking data (including 28 sites, and involving 13 of the 21 30S ribosomal were used to relate the RNA structure to the distribution of the proteins as determined by neutron scattering. The three-dimensional model of the 16 S RNA has overall dimensions of 220 A x 140 A x 90 A, in good agreement with electron microscopic estimates for the 30 S subunit. The shape of the model is also recognizably the same as that seen in electron micrographs, and the positions in the model of bases localized on the 30 S subunit by immunoelectron microscopy (the 5' and 3' termini, the m7G and m6(2)A residues, and C-1400) correspond closely to their experimentally observed positions. The distances between the RNA-protein crosslink sites in the model correlate well with the distances between protein centres of mass obtained by neutron scattering, only two out of 66 distances falling outside the expected tolerance limits. These two distances both involve protein S13, a protein noted for its anomalous behaviour. A comparison with other experimental information not specifically used in deriving the model shows that it fits well with published data on RNA-protein binding sites, mutation sites on the RNA causing resistance to antibiotics, tertiary interactions in the RNA, and a potential secondary structural "switch". Of the sites on 16 S RNA that have been found to be accessible to chemical modification in the 30 S subunit, 87% are at obviously exposed positions in the model. In contrast, 70% of the sites corresponding to positions that have ribose 2'-O-methylations in the eukaryotic 18 S RNA from Xenopus laevis are at non-exposed (i.e. internal) positions in the model. All nine of the modified bases in the E. coli 16 S RNA itself show a remarkable distribution, in that they form a "necklace" in one plane around the "throat" of the subunit. Insertions in eukaryotic 18 S RNA, and corresponding deletions in chloroplast or mammalian mitochondrial ribosomal RNA relative to E. coli 16 S RNA represent distinct sub-domains in the structure.(ABSTRACT TRUNCATED AT 400 WORDS)     

The lonepair triloop: a new motif in RNA structure

The lonepair triloop (LPTL) is an RNA structural motif that contains a single ("lone") base-pair capped by a hairpin loop containing three nucleotides. The two nucleotides immediately outside of this motif (5' and 3' to the lonepair) are not base-paired to one another, restricting the length of this helix to a single base-pair. Four examples of this motif, along with three tentative examples, were initially identified in the 16S and 23S rRNAs with covariation analysis. An evaluation of the recently determined crystal structures of the Thermus thermophilus 30S and Haloarcula marismortui 50S ribosomal subunits revealed the authenticity for all of these proposed interactions and identified 16 more LPTLs in the 5S, 16S and 23S rRNAs. This motif is found in the T loop in the tRNA crystal structures. The lonepairs are positioned, in nearly all examples, immediately 3' to a regular secondary structure helix and are stabilized by coaxial stacking onto this flanking helix. In all but two cases, the nucleotides in the triloop are involved in a tertiary interaction with another section of the rRNA, establishing an overall three-dimensional function for this motif. Of these 24 examples, 14 occur in multi-stem loops, seven in hairpin loops and three in internal loops. While the most common lonepair, U:A, occurs in ten of the 24 LPTLs, the remaining 14 LPTLs contain seven different base-pair types. Only a few of these lonepairs adopt the standard Watson-Crick base-pair conformations, while the majority of the base-pairs have non-standard conformations. While the general three-dimensional conformation is similar for all examples of this motif, characteristic differences lead to several subtypes present in different structural environments. At least one triloop nucleotide in 22 of the 24 LPTLs in the rRNAs and tRNAs forms a tertiary interaction with another part of the RNA. When a LPTL containing the GNR or UYR triloop sequence forms a tertiary interaction with the first (and second) triloop nucleotide, it recruits a fourth nucleotide to mediate stacking and mimic the tetraloop conformation. Approximately half of the LPTL motifs are in close association with proteins. The majority of these LPTLs are positioned at sites in rRNAs that are conserved in the three phylogenetic domains; a few of these occur in regions of the rRNA associated with ribosomal function, including the presumed site of peptidyl transferase activity in the 23S rRNA.     

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.     

Intra-RNA cross-linking in Escherichia coli 30S ribosomal subunits: selective isolation of cross-linked products by hybridization to specific cDNA fragments.

M13 clones were constructed with cDNA inserts corresponding to specific regions of E. coli ribosomal RNA. The DNA from the clones was immobilized by coupling to diazobenzyloxymethyl cellulose, and was used for the selective isolation by hybridization of cross-linked RNA complexes containing the complementary sequences. Immobilized DNA samples with inserts complementary to four different regions covering bases 735-1384 of the 16S RNA were hybridized with a mixture of 16S RNA fragments generated by partial digestion of 30S subunits that had been cross-linked by ultraviolet irradiation in vivo. After dehybridization, the individual RNA fragments and cross-linked complexes were separated by gel electrophoresis and analysed by our usual procedures. Nine cross-links are described; four of these are hitherto unobserved "secondary structural" cross-links, and one is a new "tertiary structural" cross-link between positions 243-247 and 891-894 of the 16S RNA.     

Prediction of three-dimensional structure of Escherichia coli ribosomal RNA

A model for the tertiary structure of 23S, 16S and 5S ribosomal RNA molecules interacting with three tRNA molecules is presented using the secondary structure models common to E. coli, Z. mays chloroplast, and mammalian mitochondria. This ribosomal RNA model is represented by phosphorus atoms which are separated by 5.9 A in the standard A-form double helix conformation. The accumulated proximity data summarized in Table 1 were used to deduce the most reasonable assembly of helices separated from each other by at least 6.2 A. Straight-line approximation for single strands was adopted to describe the maximum allowed distance between helices. The model of a ribosome binding three tRNA molecules by Nierhaus (1984), the stereochemical model of codon-anticodon interaction by Sundaralingam et al. (1975) and the ribosomal transpeptidation model, forming an alpha-helical nascent polypeptide, by Lim & Spirin (1986), were incorporated in this model. The distribution of chemically modified nucleotides, cross-linked sites, invariant and missing regions in mammalian mitochondrial rRNAs are indicated on the model.     

Comparison of eubacterial and eukaryotic 5S RNA structures: a chemical modification study

The 5S RNAs from Bacillus stearothermophilus and Saccharomyces cerevisiae were probed by nucleotide-specific reagents, with a view to compare and contrast their higher order structures. The progressive unfolding of the RNAs during heating, in the presence and absence of magnesium, was monitored. Evidence was provided for the double-helical segments which occur in the secondary structural models of both RNAs. The results also placed constraints on the possible structuring of the remainder of the RNA and yielded some insight into ways of folding up the molecule. Together with the data from our earlier studies, employing ribonucleases, these results provide a detailed picture of the structuring and topography of the 5S RNAs. The main structural differences between the eubacterial and eukaryotic RNAs occur throughout the loop D/helix IV/loop E/helix V arm; in particular strong evidence is provided for loop D of the eukaryotic RNA being involved in a tertiary interaction.