Prediction of three-dimensional structure of Escherichia coli ribosomal RNA

The 3D arrangement of the 23 S and 5 S rRNA in the Escherichia coli 50 S ribosomal subunit based on a cryo-electron microscopic reconstruction at 7.5 A resolution

The Escherichia coli 23 S and 5 S rRNA molecules have been fitted helix by helix to a cryo-electron microscopic (EM) reconstruction of the 50 S ribosomal subunit, using an unfiltered version of the recently published 50 S reconstruction at 7.5 A resolution. At this resolution, the EM density shows a well-defined network of fine structural elements, in which the major and minor grooves of the rRNA helices can be discerned at many locations. The 3D folding of the rRNA molecules within this EM density is constrained by their well-established secondary structures, and further constraints are provided by intra and inter-rRNA crosslinking data, as well as by tertiary interactions and pseudoknots. RNA-protein cross-link and foot-print sites on the 23 S and 5 S rRNA were used to position the rRNA elements concerned in relation to the known arrangement of the ribosomal proteins as determined by immuno-electron microscopy. The published X-ray or NMR structures of seven 50 S ribosomal proteins or RNA-protein complexes were incorporated into the EM density. The 3D locations of cross-link and foot-print sites to the 23 S rRNA from tRNA bound to the ribosomal A, P or E sites were correlated with the positions of the tRNA molecules directly observed in earlier reconstructions of the 70 S ribosome at 13 A or 20 A. Similarly, the positions of cross-link sites within the peptidyl transferase ring of the 23 S rRNA from the aminoacyl residue of tRNA were correlated with the locations of the CCA ends of the A and P site tRNA. Sites on the 23 S rRNA that are cross-linked to the N termini of peptides of different lengths were all found to lie within or close to the internal tunnel connecting the peptidyl transferase region with the presumed peptide exit site on the solvent side of the 50 S subunit. The post-transcriptionally modified bases in the 23 S rRNA form a cluster close to the peptidyl transferase area. The minimum conserved core elements of the secondary structure of the 23 S rRNA form a compact block within the 3D structure and, conversely, the points corresponding to the locations of expansion segments in 28 S rRNA all lie on the outside of the structure.

A new model for the three-dimensional folding of Escherichia coli 16 S ribosomal RNA. I. Fitting the RNA to a 3D electron microscopic map at 20 A

Recently published models of the Escherichia coli 70 S ribosome at 20 A resolution, obtained by cryo-electron microscopy (cryo-EM) combined with computerized image processing techniques, exhibit two features that are directly relevant to the in situ three-dimensional folding of the rRNA molecules. First, at this level of resolution many fine structural details are visible, a number of them having dimensions comparable to those of nucleic acid helices. Second, in reconstructions of ribosomes in the pre- and post-translocational states, density can be seen that corresponds directly to the A and P site tRNAs, and to the P and E site tRNAs, respectively, thus enabling the decoding region on the 30 S subunit to be located rather precisely. Accordingly, we have refined our previous model for the 16 S rRNA, based on biochemical evidence, by fitting it to the cryo-EM contour of ribosomes carrying A and P site tRNAs. For this purpose, the most immediately relevant evidence consists of new site-directed cross-linking data in the decoding region, which define sets of contacts between the 16 S rRNA and mRNA, or between 16 S rRNA and tRNA at the A, P and E sites; these contact sites can be correlated directly with the tRNA positions in the EM structure. The model is extended to other parts of the 16 S molecule by fitting individual elements of the well-established secondary structure of the 16 S rRNA into the appropriate fine structural elements of the EM contour, at the same time taking into account other data used in the previous model, such as intra-RNA cross-links within the 16 S rRNA itself. The large body of available RNA-protein cross-linking and foot-printing data is also considered in the model, in order to correlate the rRNA folding with the known distribution of the 30 S ribosomal proteins as determined by neutron scattering and immuno-electron microscopy. The great majority of the biochemical data points involve single-stranded regions of the rRNA, and therefore, in contrast to most previous models, the single-stranded regions are included in our structure, with the help of a specially developed modelling programme, ERNA-3D. This allows the various biochemical data sets to be displayed directly, in this and in the accompanying papers, on diagrams of appropriate parts of the rRNA structure within the cryo-EM contour.

Localization of two domains of a mutant form of Escherichia coli protein L7/L12 that binds the large ribosomal subunit as a single dimer

Escherichia coli ribosomal protein L7/L12 occurs on the large subunit as two dimers: one dimer is extended and comprises the stalk, while the second dimer is folded and occupies a site on the subunit body. A variant protein, in which all 18 amino acids of the flexible hinge region that links separate N-terminal and C-terminal domains of L7/L12 has been deleted, binds the subunit as a single dimer and does not generate stalks that are visible in electron micrographs. Monoclonal antibodies directed against each domain of the protein have been used to localize the variant in electron micrographs of 50S subunits. Both C-terminal domains are seen at a shoulder of the subunit, near its edge as viewed in the most common quasisymmetric projection. N-terminal domains are placed on the subunit body, about 50 A from the C-terminal domains. The antibody to the N-terminal domain also causes dissociation of the variant dimer from the particle and the formation of oligomeric antibody-protein dimer complexes. Similar complexes were seen previously (Olson HM et al (1986) J Biol Chem 261, 6924-6936) when this antibody induced dissociation of one dimer of the native protein. We conclude that the shortened variant most probably occupies the lower-affinity site on the subunit that is normally filled by the stalk dimer.

Incorporation of dinitrophenyl protein L23 into totally reconstituted Escherichia coli 50 S ribosomal subunits and its localization at two sites by immune electron microscopy

Escherichia coli ribosomal protein L23 was derivatized with [3H]2, 4-dinitrofluorobenzene both at the N terminus and at internal lysines. Dinitrophenyl-L23 (DNP-L23) was taken up into 50 S subunits from a reconstitution mixture containing rRNA and total 50 S protein depleted in L23. Unmodified L23 competed with DNP-L23 for uptake, indicating that each protein form bound in an identical or similar position within the subunit. Modified L23, incorporated at a level of 0.7 or 0.4 DNP groups per 50 S, was localized by electron microscopy of subunits complexed with antibodies to dinitrophenol. Antibodies were seen at two major sites with almost equal frequency. One site is beside the central protuberance, in a region previously identified as the peptidyltransferase center. The second location is at the base of the subunit, in the area of the exit site from which the growing peptide leaves the ribosome. Models derived from image reconstruction show hollows or canyons in the subunit and a tunnel that links the transferase and exit sites. Our results indicate that L23 is at the subunit interior, with separate elements of the protein at the subunit surface at or near both ends of this tunnel.

Transfer RNA docking pair model in the ribosomal pre- and post-translocational states

A consensus has been reached that the conformation of the anticodon-codon interactions of two adjacent tRNA molecules on the ribosome is a Sundaralingam-type (S-type). Even if it is kept to the S-type, there are still various possibilities. Various experimental data have been supporting an idea that the conformation of A-site tRNA is different from that of P-site tRNA. Those data as well as the recent result of Brimacombe and co-workers that U20:1 of lupin tRNAmMetbound to the A-site was cross-linked to a region, 875-905, of 23S rRNA in combination with the other recent findings of Nierhaus and co-workers about the spin-contrast method of neutron diffraction of the ribosome and the better accessible nucleotide patterns of phosphorothioated tRNAs on the ribosome have led to a new tRNA docking pair model, in which the highly conserved G18 and G19 of D-loop in A-site tRNA and C56 and C61 of TpsiC-loop in P-site tRNA base pair along with the conventional base pairs of adjacent codon-anticodon interactions. This A-P tRNA pair model can be translocated to the P-E tRNA pair model without changing the conformation except the ACCA termini, keeping the position of the growing nascent polypeptide chain.

Modeling the structure of the ribosome

Considering the size and complexity of the ribosome and the growing body of data from a wide range of experiments on ribosomal structure, it is becoming increasingly important to develop tools that facilitate the development of reliable models for the ribosome. We use a combination of manual and computer-based approaches for building and refining models of the ribosome and other RNA-protein complexes. Our methods are aimed at determining the range of models compatible with the data, making quantitative statements about the positional uncertainties (resolution) of different regions, identifying conflicts in the data, establishing which regions of the ribosome need further experimental exploration, and, where possible, predicting the outcome of future experiments. Our previous low-resolution model for the small subunit of the Escherichia coli ribosome is briefly reviewed, along with progress on atomic resolution modeling of the mRNA-tRNA complex and its interaction with the decoding site of the 16S RNA.

The 70S Escherichia coli ribosome at 23 A resolution: fitting the ribosomal RNA

BACKGROUND

The ribosome--essential for protein synthesis in all organisms--has been an evasive target for structural studies. The best available structures for the 70S Escherichia coli ribosome or its 30S and 50S subunits are based on electron microscopical tilt experiments and are limited in resolution to 28-55 A. The angular reconstitution approach, which exploits the random orientations of particles within a vitreous ice matrix, can be used in conjunction with cryo-electron microscopy to yield a higher-resolution structure.

RESULTS

Our 23 A resolution map of the 70S ribosome elucidates many structural details, such as an extensive system of channels within the 50S subunit and an intersubunit gap ideally shaped to accommodate two transfer RNA molecules. The resolution achieved is sufficient to allow the preliminary fitting of double-helical regions of an earlier three-dimensional ribosomal RNA model.

CONCLUSIONS

Although we are still a long way from attaining an atomic-resolution structure of the ribosome, cryo-electron microscopy, in combination with angular reconstitution, is likely to yield three-dimensional maps with gradually increasing resolution. As exemplified by our current 23 A reconstruction, these maps will lead to progressive refinement of models of the ribosomal RNA.

Structures of small subunit ribosomal RNAs in situ from Escherichia coli and Thermomyces lanuginosus

Small ribosomal subunits from the prokaryote Escherichia coli and the eukaryote Thermomyces lanuginosus were imaged electron spectroscopically, and single particle analysis used to yield three-dimensional reconstructions of the net phosphorus distribution representing the nucleic acid (RNA) backbone. This direct approach showed both ribosomal RNAs to have a three domain structure and other characteristic morphological features. The eukaryotic small ribosomal subunit had a prominent bill present in the head domain, while the prokaryotic subunit had a small vestigial bill. Both ribosomal subunits contained a thick 'collar' central domain which correlates to the site of the evolutionarily conserved ribosomal RNA core, and the location of the majority of ribosomal RNA bases that have been implicated in translation. The reconstruction of the prokaryotic subunit had a prominent protrusion extending from the collar, forming a channel approximately 1.5 nm wide and potentially representing a 'bridge' to the large subunit in the intact monosome. The basal domain of the prokaryotic ribosomal subunit was protein free. In this region of the eukaryotic subunit, there were two basal lobes composed of ribosomal RNA, consistent with previous hypotheses that this is a site for the 'non-conserved core' ribosomal RNA.

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.

Identification of cross-linked amino acids in the protein pair HmaL23-HmaL29 from the 50S ribosomal subunit of the archaebacterium Haloarcula marismortui

50S ribosomal subunits from the extreme halophilic archaebacterium Haloarcula marismortui were treated with the homobifunctional protein-protein cross-linking reagents diepoxybutane (4 A) and dithiobis(succinimidyl propionate) (12 A). The dominant product with both cross-linking reagents was identified on the protein level as HmaL23-HmaL29, which is homologous to the protein pair L23-L29 from Escherichia coli [Walleczek, J., Martin, T., Redl, B., Stöffler-Meilicke, M., & Stöffler, G. (1989) Biochemistry 28, 4099-4105] and from Bacillus stearothermophilus [Brockmöller, J., & Kamp, R. M. (1986) Biol. Chem. Hoppe-Seyler 367, 925-935]. To reveal the exact cross-linking site in HmaL23-HmaL29, the cross-linked complex was purified on a preparative scale by conventional and high-performance liquid chromatography. After endoproteolytic fragmentation of the protein pair, the amino acids engaged in cross-link formation were unambiguously identified by N-terminal sequence analysis and mass spectrometry of the cross-linked peptides. The cross-link is formed between lysine-57 in the C-terminal region of HmaL29 and the alpha-amino group of the N-terminal serine in protein HmaL23, irrespective of the cross-linking reagent. This result demonstrates that the N-terminal region of protein HmaL23 and the C-terminal domain of HmaL29 are highly flexible so that the distance between the two polypeptide chains can vary by at least 8 A. Comparison of our cross-linking results with those obtained with B. stearothermophilus revealed that the fine structure within this ribosomal domain is at least partially conserved.

Electron microscopic visualisation of the 5S rRNA-YL3 complex from Saccharomyces cerevisiae

The complex comprising 5S ribosomal RNA and the ribosomal protein YL3 (5S rRNP) was isolated from yeast (Saccharomyces cerevisiae), and positively contrasted preparations were imaged by transmission electron microscopy. The overall dimensions of the 5S rRNP complex in the micrographs were 10 nm by 6 nm. Three predominant projections were selected from several hundred putative particles for digitisation and computer averaging to yield two-dimensional constructions with reproducible spatial resolutions exceeding 2 nm. The enhanced projection images were compatible with structural models of this complex based on biochemical studies.

The mRNA binding track in the Escherichia coli ribosome for mRNAs of different sequences

Interactions between mRNA and rRNA on the 30S ribosomal subunit or 70S ribosome have been determined by photochemical cross-linking experiments using synthetic mRNA analogs substituted with 4-thiouridine. A set of RNA molecules containing different sequences has been used to determine the extent to which binding contacts are sequence dependent. The 16S rRNA and 23S rRNA nucleotides that form a part of the binding site have been identified by reverse transcription. The nucleotides are U1381, G1338, G1300, G1156, A845, U723, G693, A532, G497, U420, G413/A412, and G436 of 16S rRNA and U887 of 23S rRNA. Several additional nucleotides (U1065 of 23S rRNA and A1227, G818, G524, and G423 of 16S rRNA) are seen for some, but not all, of the mRNAs. Results obtained with two mRNAs containing the Shine-Dalgarno sequence were similar to those obtained with mRNAs lacking the Shine-Dalgarno sequence. Eight of these cross-linking sites were also seen when a mixture of RNA was used in which there are 12 random nucleotides preceding and seven random nucleotides succeeding an AUG codon. These results indicate that to a large extent placement of the mRNA in the ribosome does not depend upon its primary sequence.

Localization of a segment of 16S RNA on the surface of the small ribosomal subunit by immune electron microscopy of complementary oligodeoxynucleotides

Oligonucleotides that complement Escherichia coli 16S ribosomal RNA residues 685-696 and 694-705 have been synthesized so as to incorporate antibody-recognizable markers: a 3'-terminal residue of N6-delta 2-isopentenyladenosine, a 5'-dinitrophenyl group, or both. Each oligonucleotide is able to bind RNA within the small ribosomal subunit, whether free or in 70S ribosomes. Immune electron microscopy places probes at nucleotides 685, 694 and 705 within a single area, at the tip of the subunit platform, very near the position of the 3'-end of the 16S RNA.

The side-by-side model of two tRNA molecules allowing the alpha-helical conformation of the nascent polypeptide during the ribosomal transpeptidation

Lim and Spirin [25] proposed a preferable conformation of the nascent peptide during the ribosomal transpeptidation. Spirin and Lim [26] excluded the possibilities of the side-by-side model proposed by Johnson et al [13] and the three-tRNA binding model (A, P and E sites) of Rheinberger and Nierhaus [3]. However, a slight conformational change at the 3' end regions of both A and P site tRNA molecules can enable the three different tRNA binding models to converge. With a modification of the angles of the ribose rings of both anticodon and mRNA this model can also be related to the model of Sundaralingam et al [19]. In this model of E coli rRNA the 3' end sequence ACCA76 or GCCA76 of P site tRNA is base-paired to UGGU810 of 23S rRNA, while the ACC75 or GCC75 of A site tRNA are base-paired to GGU1621 23S rRNA. The conformation of the A76 of A site tRNA is necessarily different from that of P site tRNA, at least during the course of the transpeptidation. The A76 of A site tRNA overlaps the binding region of puromycin. The C1400 of 16S rRNA in this model is located at a distance of 4 A from the 5' end of the anticodon of P site tRNA [14] and 17 A from the 5' end of the anticodon of A site tRNA [15]. It is also shown that a considerable but reasonable modification in the conformation of the anticodon loops could lead to accommodation of three deacylated tRNA(Phe) molecules at a time on 70S ribosome in the presence of poly(U) as observed experimentally [6]. A sterochemical explanation for the negatively-linked allosteric interactions between the A and E sites is also shown in the present model.

Explicit distance geometry: identification of all the degrees of freedom in a large RNA molecule

An alternative approach to distance geometry ("explicit" distance geometry) is being developed for problems, such as the modeling of RNA folding in the ribosome, where relatively few distances are known. The approach explicitly identifies minimal sets of additional distances that can be added to a distance matrix in order to calculate structures that are consistent with all the known information without distorting the original input data. These additional distances are bounded to the extent possible by the known distances. These explicitly added distances can be treated as degrees of freedom and used to explore the full range of alternative foldings consistent with the original input in an organized way. The present paper establishes that it is practical to explicitly determine such degrees of freedom for even very large RNAs. To demonstrate the feasibility of the approach tRNA was represented as a simple undirected graph containing all relevant information represented in the usual cloverleaf secondary structure and nine base-base tertiary interactions. Using a three atom representation for each residue a total of 206 degrees of freedom are explicitly identified. To accomplish this a graph theoretic approach was used in which a minimal covering cycle basis was determined.

Prediction of the three-dimensional structure of Escherichia coli 30S ribosomal subunit: a molecular mechanics approach

We introduce a computer-assisted procedure for folding large RNA chains into three-dimensional conformations consistent with their secondary structure and other known experimental constraints. The RNA chain is modeled using pseudoatoms at different levels of detail--from a single pseudoatom per helix to a single pseudoatom for each nucleotide. A stepwise procedure is used, starting with a simple representation of the macromolecule that is refined and then extrapolated into higher resolution for further refinement. The procedure is capable of folding different random-walk chains by using energy minimization, allowing generation of a range of conformations consistent with given experimental data. We use this procedure to generate several possible conformations of the 16S RNA in the 30S ribosomal subunit of Escherichia coli by using secondary structure and the neutron-scattering map of the 21 proteins in the small subunit. The RNA chain is modeled using a single pseudoatom per helix. RNA-RNA and RNA-protein crosslinks, reported in current literature, are included in our model. Footprinting data for different ribosomal proteins in the 16S RNA are also used. Several conformations of the 16S RNA are generated and compared to predict gross structural features of the 30S subunit as well as to identify regions of the 16S RNA that cannot be well-defined with current experimental data.

Localization of an oligodeoxynucleotide complementing 16S ribosomal RNA residues 520-531 on the small subunit of Escherichia coli ribosomes: electron microscopy of ribosome-cDNA-antibody complexes

The oligodeoxynucleotide dACCGCGGCTGCT, complementary to Escherichia coli small ribosomal subunit RNA residues 520-531, has been used to probe subunit conformation and to localize the sequence in the subunit. Conditions for binding of the cDNA to 30S subunits were optimized and specificity of the interaction was demonstrated by RNase H cleavage. Three kinds of terminal modification of this cDNA were used to allow its localization by immune electron microscopy. A solid phase support with 5'-dimethoxytrity-N6-delta 2-isopentenyl-adenosine linked to controlled pore glass was synthesized, and used to prepare oligomer with an added 3'-terminal residue of isopentenyl adenosine. cDNA with a 5' primary amine substituent was modified with 1-fluoro-2,4-dinitrobenzene to prepare 5'-dinitrophenyl oligonucleotide, and both modifications together gave doubly-derivatized probes. Immune electron microscopy with antibodies to dinitrophenol, isopentenyl adenosine, or both, was used to place the cDNA on 30S subunits. In each case the probe was placed at a single site at the junction of the head and body of the subunit, near the decoding site and the area in which elongation factor Tu is bound. It is proposed that this segment of ribosomal RNA functions in mRNA binding and orientation.

Distances as degrees of freedom

Tertiary contact distance information of varying resolution for large biological molecules abounds in the literature. The results provided herein develop a framework by which information of this type can be used to reduce the allowable configuration space of a macromolecule. The approach combines graph theory and distance geometry. Large molecules are represented as simple, undirected graphs, with atoms, or groups, as vertices, and distances between them as edges. It is shown that determination of the exact structure of a molecule in three dimensions only requires the specification of all the distances in a single tetrahedron, and four distances to every other atom. This is 4N-10 distances which is a subset of the total N(N-1)/2 unique distances in a molecule consisting of N atoms. This requirement for only 4N-10 distances has serious implications for distance geometry implementations in which all N(N-1)/2 distances are specified by bounded random numbers. Such distance matrices represent overspecified systems which when solved lead to non-obvious distribution of any error caused by inherent contradictions in the input data. It is also shown that numerous valid subsets of 4N-10 distances can be constructed. It is thus possible to tailor a subset of distances using all known distances as degrees of freedom, and thereby reduce the configuration space of the molecule. Simple algebraic relationships are derived that relate sets of distances, and complicated rotations are avoided. These relationships are used to construct minimum, complete sets of distances necessary to specify the exact structure of the entire molecule in three dimensions from incomplete distance information, and to identify sets of inconsistent distances. The method is illustrated for the flexible structural types present in large ribosomal RNAs: 1.) A five-membered ring; 2.) a chemically bonded chain with its ends in contact (i.e., a hairpin loop); 3.) the spatial orientation of two separate molecules, and; 4.) an RNA helix that can have variation in individual base pairs, giving rise to global deviation from standardized helical forms.