NMR structure determination of the binding site for ribosomal protein S8 from Escherichia coli 16 S rRNA

Structure and Dynamics of the Tetra-A Loop and (A-A)-U Sequence Motif within the Coliphage GA Replicase RNA Operator

The three-dimensional structure of a RNA hairpin containing the RNA operator binding site for bacteriophage GA coat protein is presented. The phage GA operator contains the asymmetric (A-A)-U sequence motif and is capped by a four-adenine (tetra-A) loop. The uridine of the (A-A)-U motif preferentially pairs with the 5'-proximal cross-strand adenine, and the 3'-proximal adenine stacks into the helix. The tetra-A loop is well-ordered with adenine residues 2-4 forming a 3' stack. This loop conformation stands in contrast to the structure of the 5'-AUUA loop of the related phage MS2 operator in which residues 1 and 2 form a 5' stack. The context dependence of the (A-A)-U sequence motif conformation was examined using structures of 76 unique occurrences from the Protein Data Bank. The motif almost always has one adenine bulged and the other adenine adopting an A-U base pair. In the case in which the (A-A)-U motif is flanked by only one Watson-Crick base pair, the adenine adjacent to the flanking base pair tends to bulge; 80% of motifs with a 3' flanking pair have a 3' bulged adenine, and 84% of motifs with a 5' flanking pair have a 5' bulged adenine. The frequencies of 3'- and 5'-proximal adenines bulging are 33 and 67%, respectively, when the (A-A)-U motif is flanked by base pairs on both sides. Although a 3' flanking cytidine correlates (88%) with bulging of the 5'-proximal adenine, no strict dependence on flanking nucleotide identity was identified for the 5' side.

Structure analysis of free and bound states of an RNA aptamer against ribosomal protein S8 from Bacillus anthracis

Several protein-targeted RNA aptamers have been identified for a variety of applications and although the affinities of numerous protein-aptamer complexes have been determined, the structural details of these complexes have not been widely explored. We examined the structural accommodation of an RNA aptamer that binds bacterial r-protein S8. The core of the primary binding site for S8 on helix 21 of 16S rRNA contains a pair of conserved base triples that mold the sugar-phosphate backbone to S8. The aptamer, which does not contain the conserved sequence motif, is specific for the rRNA binding site of S8. The protein-free RNA aptamer adopts a helical structure with multiple non-canonical base pairs. Surprisingly, binding of S8 leads to a dramatic change in the RNA conformation that restores the signature S8 recognition fold through a novel combination of nucleobase interactions. Nucleotides within the non-canonical core rearrange to create a G-(G-C) triple and a U-(A-U)-U quartet. Although native-like S8-RNA interactions are present in the aptamer-S8 complex, the topology of the aptamer RNA differs from that of the helix 21-S8 complex. This is the first example of an RNA aptamer that adopts substantially different secondary structures in the free and protein-bound states and highlights the remarkable plasticity of RNA secondary structure.

The structure of Aquifex aeolicus ribosomal protein S8 reveals a unique subdomain that contributes to an extremely tight association with 16S rRNA

The assembly of ribonucleoprotein complexes occurs under a broad range of conditions, but the principles that promote assembly and allow function at high temperature are poorly understood. The ribosomal protein S8 from Aquifex aeolicus (AS8) is unique in that there is a 41-residue insertion in the consensus S8 sequence. In addition, AS8 exhibits an unusually high affinity for the 16S ribosomal RNA, characterized by a picomolar dissociation constant that is approximately 26,000-fold tighter than the equivalent interaction from Escherichia coli. Deletion analysis demonstrated that binding to the minimal site on helix 21 occurred at the same nanomolar affinity found for other bacterial species. The additional affinity required the presence of a three-helix junction between helices 20, 21, and 22. The crystal structure of AS8 was solved, revealing the helix-loop-helix geometry of the unique AS8 insertion region, while the core of the molecule is conserved with known S8 structures. The AS8 structure was modeled onto the structure of the 30S ribosomal subunit from E. coli, suggesting the possibility that the unique subdomain provides additional backbone and side-chain contacts between the protein and an unpaired base within the three-way junction of helices 20, 21, and 22. Point mutations in the protein insertion subdomain resulted in a significantly reduced RNA binding affinity with respect to wild-type AS8. These results indicate that the AS8-specific subdomain provides additional interactions with the three-way junction that contribute to the extremely tight binding to ribosomal RNA.

Three-dimensional motifs from the SCOR, structural classification of RNA database: extruded strands, base triples, tetraloops and U-turns

Release 2.0.1 of the Structural Classification of RNA (SCOR) database, http://scor.lbl.gov, contains a classification of the internal and hairpin loops in a comprehensive collection of 497 NMR and X-ray RNA structures. This report discusses findings of the classification that have not been reported previously. The SCOR database contains multiple examples of a newly described RNA motif, the extruded helical single strand. Internal loop base triples are classified in SCOR according to their three-dimensional context. These internal loop triples contain several examples of a frequently found motif, the minor groove AGC triple. SCOR also presents the predominant and alternate conformations of hairpin loops, as shown in the most well represented tetraloops, with consensus sequences GNRA, UNCG and ANYA. The ubiquity of the GNRA hairpin turn motif is illustrated by its presence in complex internal loops.

Beyond nucleic acid base pairs: from triads to heptads

Hydrogen-bonded base pairs are an important determinant of nucleic acid structure and function. However, other interactions such as base-base stacking, base-backbone, and backbone-backbone interactions as well as effects exerted by the solvent and by metal or NH(4)(+) ions also have to be taken into account. In addition, hydrogen-bonded base complexes involving more than two bases can occur. With the rapidly increasing number and structural diversity of nucleic acid structures known at atomic detail higher-order hydrogen-bonded base complexes, base polyads, have attracted much interest. This review provides an overview on the occurrence of base polyads in nucleic acid structures and describes computational studies on these nucleic acid building blocks.

Construction and analysis of base-paired regions of the 16S rRNA in the 30S ribosomal subunit determined by constraint satisfaction molecular modelling

Structure models for each of the secondary structure regions from the Escherichia coli 16S rRNA (58 separate elements) were constructed using a constraint satisfaction modelling program to determine which helices deviated from classic A-form geometry. Constraints for each rRNA element included the comparative secondary structure, H-bonding conformations predicted from patterns of base-pair covariation, tertiary interactions predicted from covariation analysis, chemical probing data, rRNA-rRNA crosslinking information, and coordinates from solved structures. Models for each element were built using the MC-SYM modelling algorithm and subsequently were subjected to energy minimization to correct unfavorable geometry. Approximately two-thirds of the structures that result from the input data are very similar to A-form geometry. In the remaining instances, the presence of internal loops and bulges, some sequences (and sequence covariants) and accessory information require deviation from A-form geometry. The structures of regions containing more complex base-pairing arrangements including the central pseudoknot, the 530 region, and the pseudoknot involving base-pairing between G570-U571/A865-C866 and G861-C862/G867-C868 were predicted by this approach. These molecular models provide insight into the connection between patterns of H-bonding, the presence of unpaired nucleotides, and the overall geometry of each element.

Structural features of an influenza virus promoter and their implications for viral RNA synthesis

The influenza A virus, a severe pandemic pathogen, has a segmented RNA genome consisting of eight single-stranded RNA molecules. The 5' and 3' ends of each RNA segment recognized by the influenza A virus RNA-dependent RNA polymerase direct both transcription and replication of the virus's RNA genome. Promoter binding by the viral RNA polymerase and formation of an active open complex are prerequisites for viral replication and proliferation. Here we describe the solution structure of this promoter as solved by multidimensional, heteronuclear magnetic resonance spectroscopy. Our studies show that the viral promoter has a significant dynamic nature and reveal an unusual displacement of an adenosine that forms a novel (A-A) x U motif and a C-A mismatch stacked in a helix. The characterized structural features of the promoter imply that the specificity of polymerase binding results from an internal RNA loop. In addition, an unexpected bending (46 +/- 10 degrees ) near the initiation site suggests the existence of a promoter recognition mechanism similar to that of DNA-dependent RNA polymerase and a possible regulatory function for the terminal structure during open complex formation.

Detailed analysis of RNA-protein interactions within the ribosomal protein S8-rRNA complex from the archaeon Methanococcus jannaschii

The crystal structure of ribosomal protein S8 bound to its target 16 S rRNA from a hyperthermophilic archaeon Methanococcus jannaschii has been determined at 2.6 A resolution. The protein interacts with the minor groove of helix H21 at two sites located one helical turn apart, with S8 forming a bridge over the RNA major groove. The specificity of binding is essentially provided by the C-terminal domain of S8 and the highly conserved nucleotide core, characterized by two dinucleotide platforms, facing each other. The first platform (A595-A596), which is the less phylogenetically and structurally constrained, does not directly contact the protein but has an important shaping role in inducing cross-strand stacking interactions. The second platform (U641-A642) is specifically recognized by the protein. The universally conserved A642 plays a pivotal role by ensuring the cohesion of the complex organization of the core through an array of hydrogen bonds, including the G597-C643-U641 base triple. In addition, A642 provides the unique base-specific interaction with the conserved Ser105, while the Thr106 - Thr107 peptide link is stacked on its purine ring. Noteworthy, the specific recognition of this tripeptide (Thr-Ser-Thr/Ser) is parallel to the recognition of an RNA tetraloop by a dinucleotide platform in the P4-P6 ribozyme domain of group I intron. This suggests a general dual role of dinucleotide platforms in recognition of RNA or peptide motifs. One prominent feature is that conserved side-chain amino acids, as well as conserved bases, are essentially involved in maintaining tertiary folds. The specificity of binding is mainly driven by shape complementarity, which is increased by the hydrophobic part of side-chains. The remarkable similarity of this complex with its homologue in the T. thermophilus 30 S subunit indicates a conserved interaction mode between Archaea and Bacteria.

The energetics of small internal loops in RNA

The energetics of small internal loops are important for prediction of RNA secondary and tertiary structure, selection of drug target sites, and understanding RNA structure-function relationships. Hydrogen bonding, base stacking, electrostatic interactions, backbone distortion, and base-pair size compatibility all contribute to the energetics of small internal loops. Thus, the sequence dependence of these energetics are idiosyncratic. Current approximations for predicting the free energies of internal loops consider size, asymmetry, closing base pairs, and the potential to form GA, GG, or UU pairs. The database of known three-dimensional structures allows for comparison with the models used for predicting stability from sequence.

Factors affecting the thermodynamic stability of small asymmetric internal loops in RNA

Optical melting experiments were used to determine the thermodynamic parameters for oligoribonucleotides containing small asymmetric internal loops. The results show a broad range of thermodynamic stabilities, which depend on loop size, asymmetry, sequence, closing base pairs, and length of helix stems. Imino proton NMR experiments provide evidence for possible hydrogen bonding in GA and UU mismatches in some asymmetric loops. The stabilizing effects of GA, GG, and UU mismatches on the thermodynamic stability of internal loops vary depending on the size and asymmetry of the loop. The dependence of loop stability on Watson-Crick closing base pairs may be explained by an account of hydrogen bonds. Models are presented for approximating the free energy increments of 2 x 3 and 1 x 3 internal loops.

Structure and function of the conserved 690 hairpin in Escherichia coli 16 S ribosomal RNA: analysis of the stem nucleotides

Nucleotides 680 to 710 of Escherichia coli 16 S rRNA form a distinct structural domain required for ribosome function. The goal of this study was to determine the functional significance of pairing interactions in the 690 region. Two different secondary structures were proposed for this hairpin, based on phylogenetic and chemical modification studies. To study the effect of pairing interactions in the 690 hairpin on ribosome function and to determine which of the proposed secondary structures is biologically significant, we performed an instant-evolution experiment in which the nine nucleotides that form the proposed base-pairs and dangling ends of the 690 stem were randomly mutated, and functional mutant combinations were selected. A total of 96 unique functional mutants were isolated, assayed in vivo, and sequenced. Analysis of these data revealed extensive base-pairing and stacking interactions among the mutated nucleotides. Formation of either a Watson-Crick base-pair or G.U pair between positions 688 and 699 is absolutely required for ribosome function. We also performed NMR studies of a 31-nucleotide RNA which indicate the formation of a functionally important base-pair between nucleotides 688 and 699. Formation of a second base-pair between positions 689 and 698, however, is not essential for ribosome function, but the level of ribosome function correlates with the predicted thermodynamic stability of the nucleotide pairs in these positions. The universally conserved positions G690 and U697 are generally portrayed as forming a G.U mismatch. Our data show co-variation between these positions, but do not support the hypothesis that the G690:U697 pair forms a wobble structure. NMR studies of model 14-nt and 31-nt RNAs support these findings and show that G690 and U697 are involved in unusual stacking interactions but do not form a wobble pair. Preliminary NMR structural analysis reveals that the loop portion of the 690 hairpin folds into a highly structured and novel conformation.

Determination of metal ion binding sites within the hairpin ribozyme domains by NMR

Cations play an important role in RNA folding and stabilization. The hairpin ribozyme is a small catalytic RNA consisting of two domains, A and B, which interact in the transition state in an ion-dependent fashion. Here we describe the interaction of mono-, di-, and trivalent cations with the domains of the ribozyme, as studied by homo- and heteronuclear NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping, and intermolecular NOEs indicate that the B domain contains four to five metal binding sites, which bind Mn(2+), Mg(2+), and Co(NH(3))(6)(3+). There is no significant structural change in the B domain upon the addition of Co(NH(3))(6)(3+) or Mg(2+). No specific monovalent ion binding sites exist on the B domain, as determined by (15)NH(4)(+) binding studies. In contrast to the B domain, there are no observable metal ion interactions within the internal loop of the A domain. Model structure calculations of Mn(2+) interactions at two sites within the B domain indicate that the binding sites comprise major groove pockets lined with functional groups oriented so that multiple hydrogen bonds can be formed between the RNA and Mn(H(2)O)(6)(2+) or Co(NH(3))(6)(3+). Site 1 is very similar in geometry to a site within the P4-P6 domain of the Tetrahymena group I intron, while site 2 is unique among known ion binding sites. The site 2 ion interacts with a catalytically essential nucleotide and bridges two phosphates. Due to its location and geometry, this ion may play an important role in the docking of the A and B domains.

RNA-binding properties of the mitochondrial Y-box protein RBP16

We have previously identified a mitochondrial Y-box protein in Trypanosoma brucei that we designated RBP16. The predicted RBP16 amino acid sequence revealed the presence of a cold-shock domain at its N-terminus and a glycine- and arginine-rich C-terminus reminiscent of an RGG RNA-binding motif. Since RBP16 is capable of interacting with different guide RNAs (gRNAs) in vitro and in vivo primarily via the oligo(U) tail, as well as with ribosomal RNAs, possible functions of RBP16 may be in kinetoplastid RNA editing and/or translation. Herein, we report experiments that further define the RNA-binding properties of RBP16. RBP16 forms a single stable complex with the gRNA gA6[14] at low protein concentration, while at higher protein concentration two stable complexes that possibly represent two different conformations are observed. Both complexes are stable at relatively high salt and moderate heparin concentrations indicating that the binding of RBP16 to gA6[14] does not rely primarily on ionic interactions. Phenylglyoxal treatment of the protein indicates that arginine residues are important in RNA binding. The minimal length of RNA sequence necessary for the binding of RBP16 was assessed by gel retardation and UV cross-linking competition assays using oligo(U) ribonucleotides of varying lengths (4-40 nt). Although RBP16 can bind to oligonucleotides as small as U(4), its affinity increases with the length of the oligo(U) ribonucleotide, with a dramatic increase in binding efficiency observed when the length is increased to 10 nt. Gel retardation assays employing T.brucei mRNAs demonstrated that, although it acts as a major binding determinant, a 3' U tail is not an absolute requirement for efficient RBP16-RNA binding. Experiments with oligonucleotides containing U stretches embedded at different positions in oligo(dC) indicated that high-affinity binding requires both a uridine stretch, as well as 5' and 3' non-specific sequences. These results suggest a model for the molecular interactions involved in RBP16-RNA binding.

In vivo selection of functional ribosomes with variations in the rRNA-binding site of Escherichia coli ribosomal protein S8: evolutionary implications

The highly conserved nature of rRNA sequences throughout evolution allows these molecules to be used to build philogenic trees of different species. It is unknown whether the stability of specific interactions and structural features of rRNA reflects an optimal adaptation to a functional task or an evolutionary trap. In the work reported here, we have applied an in vivo selection strategy to demonstrate that unnatural sequences do work as a functional replacement of the highly conserved binding site of ribosomal protein S8. However, growth competition experiments performed between Escherichia coli isolates containing natural and unnatural S8-binding sites showed that the fate of each isolate depended on the growth condition. In exponentially growing cells, one unnatural variant was found to be equivalent to wild type in competition experiments performed in rich media. In culture conditions leading to slow growth, however, cells containing the wild-type sequence were the ultimate winner of the competition, emphasizing that the wild-type sequence is, in fact, the most fit solution for the S8-binding site.

Solution structure and thermodynamics of a divalent metal ion binding site in an RNA pseudoknot

Identification and characterization of a metal ion binding site in an RNA pseudoknot was accomplished using cobalt (III) hexammine, Co(NH3)63+, as a probe for magnesium (II) hexahydrate, Mg(H2O)62+, in nuclear magnetic resonance (NMR) structural studies. The pseudoknot causes efficient -1 ribosomal frameshifting in mouse mammary tumor virus. Divalent metal ions, such as Mg2+, are critical for RNA structure and function; Mg2+preferentially stabilizes the pseudoknot relative to its constituent hairpins. The use of Co(NH3)63+as a substitute for Mg2+was investigated by ultraviolet absorbance melting curves, NMR titrations of the imino protons, and analysis of NMR spectra in the presence of Mg2+or Co (NH3)63+. The structure of the pseudoknot-Co(NH3)63+complex reveals an ion-binding pocket formed by a short, two-nucleotide loop and the major groove of a stem. Co(NH3)63+stabilizes the sharp loop-to-stem turn and reduces the electrostatic repulsion of the phosphates in three proximal strands. Hydrogen bonds are identified between the Co(NH3)63+protons and non-bridging phosphate oxygen atoms, 2' hydroxyl groups, and nitrogen and oxygen acceptors on the bases. The binding site is significantly different from that previously characterized in the major groove surface of tandem G.U base-pairs, but is similar to those observed in crystal structures of a fragment of the 5 S rRNA and the P5c helix of the Tetrahymena thermophila group I intron. Changes in chemical shifts occurred at the same pseudoknot protons on addition of Mg2+as on addition of Co(NH3)63+, indicating that both ions bind at the same site. Ion binding dissociation constants of approximately 0.6 mM and 5 mM (in 200 mM Na+and a temperature of 15 degrees C) were obtained for Co(NH3)63+and Mg2+, respectively, from the change in chemical shift as a function of metal ion concentration. An extensive array of non-sequence-specific hydrogen bond acceptors coupled with conserved structural elements within the binding pocket suggest a general mode of divalent metal ion stabilization of this type of frameshifter pseudoknot. These results provide new thermodynamic and structural insights into the role divalent metal ions play in stabilizing RNA tertiary structural motifs such as pseudoknots.