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.     

A functional pseudoknot in 16S ribosomal RNA.

Several lines of evidence indicate that the universally conserved 530 loop of 16S ribosomal RNA plays a crucial role in translation, related to the binding of tRNA to the ribosomal A site. Based upon limited phylogenetic sequence variation, Woese and Gutell (1989) have proposed that residues 524-526 in the 530 hairpin loop are base paired with residues 505-507 in an adjoining bulge loop, suggesting that this region of 16S rRNA folds into a pseudoknot structure. Here, we demonstrate that Watson-Crick interactions between these nucleotides are essential for ribosomal function. Moreover, we find that certain mild perturbations of the structure, for example, creation of G-U wobble pairs, generate resistance to streptomycin, an antibiotic known to interfere with the decoding process. Chemical probing of mutant ribosomes from streptomycin-resistant cells shows that the mutant ribosomes have a reduced affinity for streptomycin, even though streptomycin is thought to interact with a site on the 30S subunit that is distinct from the 530 region. Data from earlier in vitro assembly studies suggest that the pseudoknot structure is stabilized by ribosomal protein S12, mutations in which have long been known to confer streptomycin resistance and dependence.     

Interaction between 16S ribosomal RNA and ribosomal protein S12: differential effects of paromomycin and streptomycin.

Strains containing a series of restrictive and non-restrictive mutations in ribosomal protein S12 have been transformed with plasmids carrying the rrnB operon with mutations at positions 1409 and 1491 in 16S rRNA. The effects of the double-mutant constructs have been measured by growth rate, paromomycin and streptomycin sensitivity, resistance and dependence. The results demonstrate a functional interaction between the 1409-1491 region of rRNA and ribosomal protein S12.     

Interaction of proteins S16, S17 and S20 with 16 S ribosomal RNA

We have used rapid chemical probing methods to examine the effect of assembly of ribosomal proteins S16, S17 and S20 on the reactivity of individual residues of 16 S rRNA. Protein S17 strongly protects a compact region of the RNA between positions 245 and 281, a site previously assigned to binding of S20. Protein S20 also protects many of these same positions, albeit more weakly than S17. Strong S20-dependent protections are seen elsewhere in the 5' domain, most notably at positions 108, and in the 160-200 and 330 loop regions. Enenpectedly, S20 also causes protection of several bases in the 1430-1450 region, in the 3' minor domain. In the presence of the primary binding proteins S4, S8 and S20, we observe a variety of effects that result from assembly of the secondary binding protein S16. Most strongly protected are nucleotides around positions 50, 120, 300 to 330 and 360 in the 5' domain, and positions 606 to 630 in the central domain. In addition, numerous nucleotides in the 5' and central domains exhibit enhanced reactivity in response to S16. Interestingly, the strength of the S20-dependent effects in the 1430-1450 region is attenuated in the presence of S4 + S8 + S20, and restored in the presence of S4 + S8 + S20 + S16. Finally, the previously observed rearrangement of the 300 region stem-loop that occurs during assembly is shown to be an S16-dependent event. We discuss these findings with respect to assignment of RNA binding sites for these proteins, and in regard to the co-operativity of ribosome assembly.     

Nuclease S1 mapping of 16S ribosomal RNA in ribosomes

Escherichia coli 16S rRNA and 16S-like rRNAs from other species have several universally conserved sequences which are believed to be single-stranded in ribosomes. The quantitative disposition of these sequences within ribosomes is not known. Here we describe experiments designed to explore the availability of universal 16S rRNA sequences for hybridization with DNA probes in 30S particles and 70S ribosomes. Unlike previous investigations, quantitative data on the accessibility of DNA probes to the conserved portions of 16S rRNA within ribosomes was acquired. Uniquely, the experimental design also permitted investigation of cooperative interactions involving portions of conserved 16S rRNA. The basic strategy employed ribosomes, 30S subunits, and 16S rRNAs, which were quantitatively analyzed for hybridization efficiency with synthetic DNA in combination with nuclease S1. In deproteinated E. coli 16S rRNA and 30S subunits, the regions 520-530, 1396-1404, 1493-1504, and 1533-1542 are all single-stranded and unrestricted for hybridization to short synthetic DNAs. However, the quantitative disposition of the sequences in 70S ribosomes varies with each position. In 30S subunits there appear to be no cooperative interactions between the 16S rRNA universal sequences investigated.     

A Comparison Between Elastic Network Interpolation and MD Simulation of 16S Ribosomal RNA

Abstract

In this paper a coarse-grained method called elastic network interpolation (ENI) is used to generate feasible transition pathways between two given conformations of the core central domain of 16S Ribosomal RNA (16S rRNA). The two given conformations are the extremes generated by a molecular dynamics (MD) simulation, which differ from each other by 10Å in root-mean-square deviation (RMSD). It takes only several hours to build an ENI pathway on a 1.5GHz Pentium with 512 MB memory, while the MD takes several weeks on high-performance multi-processor servers such as the SGI ORIGIN 2000/2100. It is shown that multiple ENI pathways capture the essential anharmonic motions of millions of timesteps in a particular MD simulation. A coarse-grained normal mode analysis (NMA) is performed on each intermediate ENI conformation, and the lowest 1% of the normal modes (representing about 40 degrees of freedom (DOF)) are used to parameterize fluctuations. This combined ENI/NMA method captures all intermediate conformations in the MD run with 1.5Å RMSD on average. In addition, if we restrict attention to the time interval of the MD run between the two extreme conformations, the RMSD between the closest ENI/NMA pathway and the MD results is about 1Å. These results may serve as a paradigm for reducedDOF dynamic simulations of large biological macromolecules as well as a method for the reduced-parameter interpretation of massive amounts of MD data.     

Complementary binding of oligonucleotides with 16S RNA and ribosomal ribonucleoproteins

The accessibility of single-stranded sequences of 16S RNA in free state and in ribonucleoprotein particles (RNP) to complementary binding with isoplith fractions of oligonucleotides was studied. RNP had different protein composition and corresponded to intermediate stages of E. coli 30S subunit assembly in vitro. Gel-filtration was used to detect the most strong binding. It was found that S4 essentially inhibited the hexamer binding to RNA. Core proteins bound to 16S RNA strongly increased the shielding of single-stranded regions while split proteins insignificantly changed the hexamer binding. Nevertheless evidence is presented that split proteins might also interact directly with 16S RNA in the 30S subunit.     

Scanning tunnelling microscopy of 16S ribosomal RNA in water

The scanning tunnelling microscope has been used to image 16S ribosomal RNA molecules in water electrophoretically deposited on graphite surface. Two kinds of images have been obtained: images showing aggregates of 16S ribosomal RNA molecules similar to those obtained from DNA solutions and others showing individual 16S ribosomal RNA molecules. An interesting characteristic of these images, recorded in constant current mode, is that the 16S ribosomal RNA molecules appear to be located below the graphite surface. The morphology and several structural parameters of the molecules were consistent with the data obtained from electron microscopy.     

The Escherichia coli 30S ribosomal subunit; an optimized three-dimensional fit between the ribosomal proteins and the 16S RNA.

We have generated a computerized fit between the 3-dimensional map of the E.coli 30S ribosomal proteins, as determined by neutron scattering, and the recently published 3-dimensional model for the 16S RNA. To achieve this, the framework of coordinates for RNA-protein cross-link sites on the phosphate backbone in the RNA model was related to the corresponding framework of coordinates for the mass centres of the proteins by a least squares fitting procedure. The resulting structure, displayed on a computer graphics system, gives the first complete picture of the E.coli 30S ribosomal subunit showing both the proteins and the double-helical regions of the RNA. The root mean square distance between cross-link sites and protein centres is 32 A. The position of the mass centre of the combined double-helical regions was calculated from the model and compared with the position of the mass centre of the complete set of proteins. The two centres are displaced relative to one another by 20 A in the model structure, in good agreement with the experimental value of 25 A found by neutron scattering.     

Mutations in ribosomal proteins S4 and S12 influence the higher order structure of 16 S ribosomal RNA

We have studied the effects of protein mutations on the higher order structure of 16 S rRNA in Escherichia coli ribosomes, using a set of structure-sensitive chemical probes. Ten mutant strains were studied, which contained alterations in ribosomal proteins S4 and S12, including double mutants containing both altered S4 and S12. Two ribosomal ambiguity (ram) S4 mutant strains, four streptomycin resistant (SmR) S12 mutant strains, one streptomycin pseudodependent (SmP) S12 mutant strain, one streptomycin dependent (SmD) S12 mutant strain and two streptomycin independent (Sm1) double mutants (containing both-SmD and ram mutations) were probed and compared to an isogenic wild-type strain. In ribosomes from strains containing S4 ram mutations, nucleotides A8 and A26 become more reactive to dimethyl sulfate (DMS) at their N-1 positions. In ribosomes from strains bearing the SmD allele, A908, A909, A1413 and G1487 are significantly less reactive to chemical probes. These same effects are observed when the S4 and S12 mutations are present simultaneously in the double mutants. An interesting correlation is found between the reactivity of A908 and the miscoding potential of SmR, SmD, SmP and wild-type ribosomes; the reactivity of A908 increases as the translational error frequency of the ribosomes increases. In the case of ram ribosomes, the reactivity of A908 resembles that of wild-type, unless tRNA is bound, in which case it becomes hyper-reactive. Similarly, streptomycin has little effect on A908 in wild-type ribosomes unless tRNA is bound, in which case its reactivity increases to resemble that of ram ribosomes with bound tRNA. Finally, interaction of streptomycin with SmP and SmD ribosomes causes the reactivity of A908 to increase to near-wild-type levels. A simple model is proposed, in which the reactivity of A908 reflects the position of an equilibrium between two conformational states of the 30 S subunit, one of which is DMS-reactive, and the other DMS-unreactive. In this model, the balance between these two states would be influenced by proteins S4 and S12. Mutations in S12 generally cause a shift toward the unreactive conformer, and in the case of SmD and SmP ribosomes, this shift can be suppressed phenotypically by streptomycin, ram mutations in protein S4 cause a shift toward the reactive conformer, but only when tRNA is bound. This suggests that the opposing effects of these two classes of mutations influence the proof-reading process by somewhat different mechanisms.     

Base-pairing between distant regions of the Escherichia coli 16 S ribosomal RNA in solution.

Intramolecular crosslinks have been introduced into Escherichia coli 16 S ribosomal RNA in aqueous solution by irradiation in the presence of hydroxymethyl-trimethylpsoralen. When the crosslinked RNA is denatured and examined in the electron microscope the most striking features are a variety of large open loops. In addition, because the crosslinked molecules are shortened compared to non-crosslinked molecules, there are likely to be small hairpins not resolved by the present technique. The sizes and positions of 11 loop classes have been determined and oriented on the molecule. The frequency of occurrence of the different classes of loops depends on the crosslinking conditions. When the crosslinking is done in solutions containing Mg²⁺, at least four of the loop classes appear with greater frequency than they do in 3.5 mm-NaCl. The loops presumably arise because complementary sequences separated by long intervening regions are being crosslinked. These base-pairing interactions between residues distant in the primary structure appear to be prominent features of the secondary structure of rRNA in solution.     

Mutations in 16S ribosomal RNA disrupt antibiotic--RNA interactions.

Two of six mutations at a base-paired site in Escherichia coli 16S rRNA confer resistance to nine different aminoglycoside antibiotics in vivo. Chemical probing of mutant and wild-type ribosomes in the presence of paromomycin indicates that interactions between the antibiotic and 16S rRNA in mutant ribosomes are disrupted. The altered interactions measured in vitro correlate precisely with resistance seen in vivo and may be attributable to specific structural changes observed in the mutant rRNA.     

Chromosomal loci for 16S ribosomal RNA in Escherichia coli

Summary Genetic loci for 16S ribosomal RNA (rRNA) on the Escherichia coli chromosome were determined using the K-sequence, a characteristic oligonucleotide of strain K12, as a genetic marker. Oligonucleotide analyses of 16S rRNA from various recombinants between strain K12 and strain B(H) showed that the loci for 16S rRNA containing the K-sequence were near the metB locus which was at 77 min. on the chromosome map.