Predicting the three-dimensional folding of transfer RNA with a computer modeling protocol

We have developed a computer modeling protocol that can be used to predict the three-dimensional folding of a ribonucleic acid on the basis of limited amounts of secondary and tertiary data. This protocol extends the use of distance geometry beyond the domain of NMR data in which it is usually applied. The use of this algorithm to fold the molecule eliminates operator subjectivity and reproducibly predicts the overall dimensions and shape of the transfer RNA molecule. By use of a replacement pseudoatom set based on helical substructures, a series of transfer RNA foldings have been completed that utilize only the primary structure, the phylogenetically deduced secondary structure, and five long-range interactions that were determined without reference to the crystal structure. In a control set of foldings, all the interactions suspected to exist in 1969 have been included. In all cases, the modeling process consistently predicts the global arrangement of the helical domains and to a lesser extent the general path of the backbone of transfer RNA.     

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 three-dimensional structure and function of Escherichia coli ribosomal RNA, as studied by cross-linking techniques

A large number of intra-RNA and RNA-protein cross-link sites have been localized within the 23S RNA from E. coli 50 S ribosomal subunits. These sites, together with other data, are sufficient to constrain the secondary structure of the 23 S molecule into a compact three-dimensional shape. Some of the features of this structure are discussed, in particular, those relating to the orientation of tRNA on the 50 S subunit as studied by site-directed cross-linking techniques. A corresponding model for the 16S RNA within the 30 S subunit has already been described, and here a site-directed cross-linking approach is being used to determine the path followed through the subunit by messenger 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.     

Molecular dynamics investigations of hammerhead ribozyme RNA

The hammerhead ribozyme, a small catalytic RNA molecule, cleaves, in the presence of magnesium ions, a specific phosphodiester bond within its own backbone, leading to 23-cyclic phosphate and 5-OH extremities. In order to study the dynamical flexibility of the hammerhead RNA, we performed molecular dynamics simulations of the solvated crystal structure of an active hammerhead ribozyme, obtained after flash-freezing crystals soaked with magnesium. Because of a careful equilibration protocol and the use of the Ewald summation in calculating the electrostatic interactions, the RNA structure remained close to the crystal structure, as attested by a root-mean-square deviation below 2.5 A after 750 ps of simulation. All Watson-Crick base pairs were intact at the end of the simulations. The tertiary interactions, such as the sheared G.A pairs and the U-turn, important for the stabilisation of the three-dimensional RNA fold, were also retained. The results demonstrate that molecular dynamics simulations can be successfully used to investigate the dynamical behaviour of a ribozyme, thus, opening a road to study the role of transient structural changes involved in ribozyme catalysis.     

Multiple e-Pharmacophore modeling to identify a single molecule that could target both streptomycin and paromomycin binding sites for 30S ribosomal subunit inhibition

The bacterial ribosome is an established target for anti-bacterial therapy since decades. Several inhibitors have already been developed targeting both defined subunits (50S and 30S) of the ribosome. Aminoglycosides and tetracyclines are two classes of antibiotics that bind to the 30S ribosomal subunit. These inhibitors can target multiple active sites on ribosome that have a complex structure. To screen putative inhibitors against 30S subunit of the ribosome, the crystal structures in complex with various known inhibitors were analyzed using pharmacophore modeling approach. Multiple active sites were considered for building energy-based three-dimensional (3D) pharmacophore models. The generated models were validated using enrichment factor on decoy data-set. Virtual screening was performed using the developed 3D pharmacophore models and molecular interaction towards the 30S ribosomal unit was analyzed using the hits obtained for each pharmacophore model. The hits that were common to both streptomycin and paromomycin binding sites were identified. Further, to predict the activity of these hits a robust 2D-QSAR model with good predictive ability was developed using 16 streptomycin analogs. Hence, the developed models were able to identify novel inhibitors that are capable of binding to multiple active sites present on 30S ribosomal subunit.     

Comparison of X-ray crystal structure of the 30S subunit-antibiotic complex with NMR structure of decoding site oligonucleotide-paromomycin complex

Aminoglycoside antibiotics that bind to 16S ribosomal RNA in the aminoacyl-tRNA site (A site) cause misreading of the genetic code and inhibit translocation. Structures of an A site RNA oligonucleotide free in solution and bound to the aminoglycosides paromomycin or gentamicin C1a have been determined by NMR. Recently, the X-ray crystal structure of the entire 30S subunit has been determined, free and bound to paromomycin. Distinct differences were observed in the crystal structure, particularly at A1493. Here, the NMR structure of the oligonucleotide-paromomycin complex was determined with higher precision and is compared with the X-ray crystal structure of the 30S subunit complex. The comparison shows the validity of both structures in identifying critical interactions that affect ribosome function.     

Structural aspects of RbfA action during small ribosomal subunit assembly

Ribosome binding factor A (RbfA) is a bacterial cold shock response protein, required for an efficient processing of the 5' end of the 16S ribosomal RNA (rRNA) during assembly of the small (30S) ribosomal subunit. Here we present a crystal structure of Thermus thermophilus (Tth) RbfA and a three-dimensional cryo-electron microscopic (EM) map of the Tth 30S*RbfA complex. RbfA binds to the 30S subunit in a position overlapping the binding sites of the A and P site tRNAs, and RbfA's functionally important C terminus extends toward the 5' end of the 16S rRNA. In the presence of RbfA, a portion of the 16S rRNA encompassing helix 44, which is known to be directly involved in mRNA decoding and tRNA binding, is displaced. These results shed light on the role played by RbfA during maturation of the 30S subunit, and also indicate how RbfA provides cells with a translational advantage under conditions of cold shock.     

A Geometrical Approach in Folding a Pseudoknot Motif within the E. Coli 16S-RNA

Abstract

Except for tRNA, the tertiary structure of RNA molecules are very little known. The many possibilities in the arrangement of different helices in space and the flexibility in the single- stranded loops that connect the helical regions make the modeling of the tertiary structure of RNA molecule a very complex task. Here, we introduce an approach to fold RNA tertiary structure based only on the information of the secondary structure and the stereochemistry of the molecule. This approach was used to construct an atomic structure of a pseudoknot (bases 500–545) in the E. coli 16S RNA. The resulting structure is a closely packed molecule that is consistent with the predicted secondary structure and stereochemically feasible. This new approach is very general and easily adaptable. Experimental data (e.g., NMR, fluorescence energy transfer, etc.), as they become available, can be incorporated directly into the approach to improve the accuracy of the modeled structure.     

Solution structure of ribosomal protein L16 from Thermus thermophilus HB8

Ribosomal protein L16 is an essential component of the bacterial ribosome. It organizes the architecture of aminoacyl tRNA binding site in the ribosome 50S subunit. The three-dimensional structure of L16 from Thermus thermophilus HB8 was determined by NMR. In solution, L16 forms an alpha+beta sandwich structure combined with two additional beta sheets located at the loop regions connecting the two layers. The terminal regions and a central loop region did not show any specific secondary structure. The structured part of L16 could be superimposed well on the C(alpha) model of L16 determined in the crystal structure of the ribosome 50S subunit. By overlaying the L16 solution structure onto the coordinates of the ribosome crystal structure, we constructed the combined model that represents the ribosome-bound state of L16 in the detailed structure. The model showed that L16 possesses residues in contact with helices 38, 39, 42, 43 and 89 of 23S rRNA and helix 4 of 5S rRNA. This suggests its broad effect on the ribosome architecture. Comparison of L16 with the L10e protein, which is the archaeal counterpart, showed that they share a common fold, but differ in some regions of functional importance, especially in the N-terminal region. All known mutation sites in L16 that confer resistance to avilamycin and evernimicin were positioned so that their side-chains were exposed to solvent in the internal cavity of the ribosome. This suggests the direct participation of L16 as a part of the binding site for antibiotics.     

Topography of 16 S RNA in 30 S subunits and 70 S ribosomes accessibility to cobra venom ribonuclease

A ribonuclease extracted from the venom of the cobra Naja oxiana, which shows an unusual specificity for double-stranded RNA regions, was used to obtain new insight on the topography of Escherichia coli ribosomal 16 S RNA in the 30 S subunit and in the 70 S couple. ³²P-labeled 30 S subunits or reconstituted 70 S tight couples containing ³²P-labeled 16 S RNA have been digested under progressively stronger conditions. The cleavage sites have been precisely localized and the chronology of the hydrolysis process studied.The enzyme cleaves the 16 S RNA within 30 S subunits at 21 different sites, which are not uniformly distributed along the molecule. These results provide valuable information on the 16 S RNA topography and evidence for secondary structure features.The binding of the 50 S subunit markedly reduces the rate of the 16 S RNA hydrolysis and provides protection for several cleavage sites. Four of them are clustered in the 3′-terminal 200 nucleotides of the molecule, one in the middle (at position 772) and one in the 5′ domain (at position 336). Our results provide further evidence that the 3′-terminal and central regions of the RNA chain are close to each other in the ribosome structure and lie at the interface of the two subunits. They also suggest that the 5′ domain is probably not involved exclusively in structure and assembly.     

Higher order structure of the ribosomal 5 S RNA.

The structure of the ribosomal 5 S RNA was examined using Fe(II)-EDTA, a solvent-based reagent that cleaves the phosphodiester backbone of both double- and single-stranded RNA but is restricted by the three-dimensional structure. In the yeast 5 S RNA, cleavages were significantly restricted in six specific regions of the molecule; restrictions in only two of these regions were clearly dependent on a high salt/magnesium ion environment. A comparison of four RNAs of diverse origin revealed strong similarities in the cleavage profiles supporting a highly conserved higher order structure. Taken together with previous studies these data provide a more detailed modeling of the three-dimensional structure.     

Location and characteristics of ribosomal protein binding sites in the 16S RNA of Escherichia coli.

Specific binding sites for five proteins of the Escherichia coli 30S ribosomal subunit have been located within the 16S RNA. The sites are structurally diverse and range in size from 40 to 500 nucleotides; their functional integrity appears to depend upon both the secondary structure and conformation of the RNA molecule. Evidence is presented which indicates that additional proteins interact with the RNA at later stages of subunit assembly.