Role of conserved nucleotides in building the 16S rRNA binding site of E. coli ribosomal protein S8

Ancestral Sequence Reconstruction of the Ribosomal Protein uS8 and Reduction of Amino Acid Usage to a Smaller Alphabet

Understanding the origin and early evolution of proteins is important for unveiling how the RNA world developed into an RNA-protein world. Because the composition of organic molecules in the Earth's primitive environment was plausibly not as diverse as today, the number of different amino acids used in early protein synthesis is likely to be substantially less than the current 20 proteinogenic residues. In this study, we have explored the thermal stability and RNA binding of ancestral variants of the ribosomal protein uS8 constructed from a reduced-alphabet of amino acids. First, we built a phylogenetic tree based on the amino acid sequences of uS8 from multiple extant organisms and used the tree to infer two plausible amino acid sequences corresponding to the last bacterial common ancestor of uS8. Both ancestral proteins were thermally stable and bound to an RNA fragment. By eliminating individual amino acid letters and monitoring thermal stability and RNA binding in the resulting proteins, we reduced the size of the amino acid set constituting one of the ancestral proteins, eventually finding that convergent sequences consisting of 15- or 14-amino acid alphabets still folded into stable structures that bound to the RNA fragment. Furthermore, a simplified variant reconstructed from a 13-amino-acid alphabet retained affinity for the RNA fragment, although it lost conformational stability. Collectively, RNA-binding activity may be achieved with a subset of the current 20 amino acids, raising the possibility of a simpler composition of RNA-binding proteins in the earliest stage of protein evolution.

Deducing putative ancestral forms of GNRA/receptor interactions from the ribosome

Stable RNAs rely on a vast repertoire of long-range interactions to assist in the folding of complex cellular machineries such as the ribosome. The universally conserved L39/H89 interaction is a long-range GNRA-like/receptor interaction localized in proximity to the peptidyl transferase center of the large subunit of the ribosome. Because of its central location, L39/H89 likely originated at an early evolutionary stage of the ribosome and played a significant role in its early function. However, L39/H89 self-assembly is impaired outside the ribosomal context. Herein, we demonstrate that structural modularity principles can be used to re-engineer L39/H89 to self-assemble in vitro. The new versions of L39/H89 improve affinity and loop selectivity by several orders of magnitude and retain the structural and functional features of their natural counterparts. These versions of L39/H89 are proposed to be ancestral forms of L39/H89 that were capable of assembling and folding independently from proteins and post-transcriptional modifications. This work demonstrates that novel RNA modules can be rationally designed by taking advantage of the modular syntax of RNA. It offers the prospect of creating new biochemical models of the ancestral ribosome and increases the tool kit for RNA nanotechnology and synthetic biology.

Modular Organization of Residue-Level Contacts Shapes the Selection Pressure on Individual Amino Acid Sites of Ribosomal Proteins

Understanding the molecular evolution of macromolecular complexes in the light of their structure, assembly, and stability is of central importance. Here, we address how the modular organization of native molecular contacts shapes the selection pressure on individual residue sites of ribosomal complexes. The bacterial ribosomal complex is represented as a residue contact network where nodes represent amino acid/nucleotide residues and edges represent their van der Waals interactions. We find statistically overrepresented native amino acid-nucleotide contacts (OaantC, one amino acid contacts one or multiple nucleotides, internucleotide contacts are disregarded). Contact number is defined as the number of nucleotides contacted. Involvement of individual amino acids in OaantCs with smaller contact numbers is more random, whereas only a few amino acids significantly contribute to OaantCs with higher contact numbers. An investigation of structure, stability, and assembly of bacterial ribosome depicts the involvement of these OaantCs in diverse biophysical interactions stabilizing the complex, including high-affinity protein-RNA contacts, interprotein cooperativity, intersubunit bridge, packing of multiple ribosomal RNA domains, etc. Amino acid-nucleotide constituents of OaantCs with higher contact numbers are generally associated with significantly slower substitution rates compared with that of OaantCs with smaller contact numbers. This evolutionary rate heterogeneity emerges from the strong purifying selection pressure that conserves the respective amino acid physicochemical properties relevant to the stabilizing interaction with OaantC nucleotides. An analysis of relative molecular orientations of OaantC residues and their interaction energetics provides the biophysical ground of purifying selection conserving OaantC amino acid physicochemical properties.

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.

Evolutionary relationships among Chlamydophila abortus variant strains inferred by rRNA secondary structure-based phylogeny

The evolutionary relationships among known Chlamydophila abortus variant strains including the LLG and POS, previously identified as being highly distinct, were investigated based on rRNA secondary structure information. PCR-amplified overlapping fragments of the 16S, 16S-23S intergenic spacer (IS), and 23S domain I rRNAs were subjected to cloning and sequencing. Secondary structure analysis revealed the presence of transitional single nucleotide variations (SNVs), two of which occurred in loops, while seven in stem regions that did not result in compensatory substitutions. Notably, only two SNVs, in 16S and 23S, occurred within evolutionary variable regions. Maximum likelihood and Bayesian phylogeny reconstructions revealed that C. abortus strains could be regarded as representing two distinct lineages, one including the “classical” C. abortus strains and the other the “LLG/POS variant”, with the type strain B577^T possibly representing an intermediate of the two lineages. The two C. abortus lineages shared three unique (apomorphic) characters in the 23S domain I and 16S-23S IS, but interestingly lacked synapomorphies in the 16S rRNA. The two lineages could be distinguished on the basis of eight positions; four of these comprised residues that appeared to be signature or unique for the “classical” lineage, while three were unique for the “LLG/POS variant”. The U277 (E. coli numbering) signature character, corresponding to a highly conserved residue of the 16S molecule, and the unique G681 residue, conserved in a functionally strategic region also of 16S, are the most pronounced attributes (autapomorphies) of the “classical” and the “LLG/POS variant” lineages, respectively. Both lineages were found to be descendants of a common ancestor with the Prk/Daruma C. psittaci variant. Compared with the “classical”, the “LLG/POS variant” lineage has retained more ancestral features. The current rRNA secondary structure-based analysis and phylogenetic inference reveal new insights into how these two C. abortus lineages have differentiated during their evolution.

Affinity of ribosomal protein S8 from mesophilic and (hyper)thermophilic archaea and bacteria for 16S rRNA correlates with the growth temperatures of the organisms

The ribosomal protein S8 plays a pivotal role in the assembly of the 30S ribosomal subunit. Using filter binding assays, S8 proteins from mesophilic, and (hyper)thermophilic species of the archaeal genus Methanococcus and from the bacteria Escherichia coli and Thermus thermophilus were tested for their affinity to their specific 16S rRNA target site. S8 proteins from hyperthermophiles exhibit a 100-fold and S8 from thermophiles exhibit a 10-fold higher affinity than their mesophilic counterparts. Thus, there is a striking correlation of affinity of S8 proteins for their specific RNA binding site and the optimal growth temperatures of the respective organisms. The stability of individual rRNA-protein complexes might modulate the stability of the ribosome, providing a maximum of thermostability and flexibility at the growth temperature of the organism.

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 location of protein S8 and surrounding elements of 16S rRNA in the 70S ribosome from combined use of directed hydroxyl radical probing and X-ray crystallography

Ribosomal protein S8, which is essential for the assembly of the central domain of 16S rRNA, is one of the most thoroughly studied RNA-binding proteins. To map its surrounding RNA in the ribosome, we carried out directed hydroxyl radical probing of 16S rRNA using Fe(II) tethered to nine different positions on the surface of protein S8 in 70S ribosomes. Hydroxyl radical-induced cleavage was observed near the classical S8-binding site in the 620 stem, and flanking the other S8-footprinted regions of the central domain at the three-helix junction near position 650 and the 825 and 860 stems. In addition, cleavage near the 5' terminus of 16S rRNA, in the 300 region of its 5' domain, and in the 1070 region of its 3'-major domain provide information about the proximity to S8 of RNA elements not directly involved in its binding. These data, along with previous footprinting and crosslinking results, allowed positioning of protein S8 and its surrounding RNA elements in a 7.8-A map of the Thermus thermophilus 70S ribosome. The resulting model is in close agreement with the extensive body of data from previous studies using protein-protein and protein-RNA crosslinking, chemical and enzymatic footprinting, and genetics.

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.

Ribosomal RNA is the target for oxazolidinones, a novel class of translational inhibitors

Oxazolidinones are antibacterial agents that act primarily against gram-positive bacteria by inhibiting protein synthesis. The binding of oxazolidinones to 70S ribosomes from Escherichia coli was studied by both UV-induced cross-linking using an azido derivative of oxazolidinone and chemical footprinting using dimethyl sulphate. Oxazolidinone binding sites were found on both 30S and 50S subunits, rRNA being the only target. On 16S rRNA, an oxazolidinone footprint was found at A864 in the central domain. 23S rRNA residues involved in oxazolidinone binding were U2113, A2114, U2118, A2119, and C2153, all in domain V. This region is close to the binding site of protein L1 and of the 3' end of tRNA in the E site. The mechanism of action of oxazolidinones in vitro was examined in a purified translation system from E. coli using natural mRNA. The rate of elongation reaction of translation was decreased, most probably because of an inhibition of tRNA translocation, and the length of nascent peptide chains was strongly reduced. Both binding sites and mode of action of oxazolidinones are unique among the antibiotics known to act on the ribosome.

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

Many cellular processes involve the preferential interaction of an RNA molecule with a specific protein. A detailed analysis of the individual protein and RNA components of these interactions can provide unique insights into the structural features important for protein-RNA recognition. Ribosomal protein S8 of Escherichia coli plays a key role in 30 S ribosomal subunit assembly through its interaction with 16 S rRNA. The binding site for protein S8 comprises a portion of helix 21, nucleotides G588 to G604 and C634 to C651. This region forms a base-paired helix that is interrupted by a non-Watson-Crick segment composed of nine phylogenetically conserved nucleotides. We have investigated the detailed structure of the conserved segment and the interaction of this region with metal ions using NMR spectroscopy. Twenty-four of the 40 calculated structures converged to similar conformations and were grouped into two families. The main difference between the families is the orientation of the base of U641. The rms deviation between the heavy-atoms of the ten lowest-energy structures is 1.24 A. The orientations of the G597.C643 base-pair and A595.(A596.U644) base-triple within the conserved core have been defined and appear to extend the proximal segment of helix 21 into the phylogenetically conserved core. The base of A642 terminates this helix by stacking against C643 and the base of U641 forms hydrogen bonds with core nucleotides. The conserved core also contains a Mg2+-binding site that promotes stabilization of the secondary and tertiary structure elements of the core. A model for the interaction of S8 with its RNA-binding site is proposed.

Crystal structure of ribosomal protein S8 from Thermus thermophilus reveals a high degree of structural conservation of a specific RNA binding site

S8 is one of the core ribosomal proteins. It binds to 16 S RNA with high affinity and independently of other ribosomal proteins. It also acts as a translational repressor in Escherichia coli by binding to its own mRNA. The structure of Thermus thermophilus S8 has been determined by the method of multiple isomorphous replacement at 2.9 A resolution and refined to a crystallographic R-factor of 16.2% (Rfree 27.5%). The two domains of the structure have an alpha/beta fold and are connected by a long protruding loop. The two molecules in the asymmetric unit of the crystal interact through an extensive hydrophobic core and form a tightly associated dimer, while symmetry-related molecules form a joint beta-sheet of mixed type. This type of protein-protein interaction could be realized within the ribosomal assembly. A comparison of the structures of T. thermophilus and Bacillus stearothermophilus S8 shows that the interdomain loop is eight residues longer in the former and reveals high structural conservation of an extensive region, located in the C-terminal domain. From mutational studies this region was proposed earlier to be involved in specific interaction with RNA. On the basis of these data and on the comparison of the two structures of S8, it is proposed that the three-dimensional structure of specific RNA binding sites in ribosomal proteins is highly conserved among different species.

A new model for the three-dimensional folding of Escherichia coli 16 S ribosomal RNA. II. The RNA-protein interaction data

The map of the mass centres of the 21 proteins from the Escherichia coli 30 S ribosomal subunit, as determined by neutron scattering, was fitted to a cryoelectron microscopic (cryo-EM) model at a resolution of 20 A of 70 S ribosomes in the pre-translocational state, carrying tRNA molecules at the A and P sites. The fit to the 30 S moiety of the 70 S particles was accomplished with the help of the well-known distribution of the ribosomal proteins in the head, body and side lobe regions of the 30 S subunit, as determined by immuno electron microscopy (IEM). Most of the protein mass centres were found to lie close to the surface (or even outside) of the cryo-EM contour of the 30 S subunit, supporting the idea that the ribosomal proteins are arranged peripherally around the rRNA. The ribosomal protein distribution was then compared with the corresponding model for the 16 S rRNA, fitted to the same EM contour (described in an accompanying paper), in order to analyse the mutual compatibility of the arrangement of proteins and rRNA in terms of the available RNA-protein interaction data. The information taken into account included the hydroxyl radical and base foot-printing data from Noller's laboratory, and our own in situ cross-linking results. Proteins S1 and S14 were not considered, due to the lack of RNA-protein data. Among the 19 proteins analysed, 12 (namely S2, S4, S5, S7, S8, S9, S10, S11, S12, S15, S17 and S21) showed a fit to the rRNA model that varied from being excellent to at least acceptable. Of the remaining 7, S3 and S13 showed a rather poor fit, as did S18 (which is considered in combination with S6 in the foot-printing experiments). S16 was difficult to evaluate, as the foot-print data for this protein cover a large area of the rRNA. S19 and S20 showed a bad fit in terms of the neutron map, but their foot-print and cross-link sites were clustered into compact groups in the rRNA model in those regions of the 30 S subunit where these proteins have respectively been located by IEM studies.

The structures and relative stabilities of d(G x G) reverse Hoogsteen, d(G x T) reverse wobble, and d(G x C) reverse Watson-Crick base-pairs in DNA crystals

We have solved the structures of the homoduplex d(Gm5CGCGCG)2, and the heteroduplexes d(GCGCGCG)/d(TCGCGCG) and d(GCGCGCG)/d(CCGCGCG). The structures form six base-pairs of identical Z-DNA duplexes with single nucleotides overhanging at the 5'-ends. The overhanging nucleotide from one strand remains stacked and sandwiched between the blunt-ends of two adjacent Z-DNA duplexes, while the overhanging base of the opposing strand is extra-helical. The stacked and the extra-helical bases from adjacent duplexes pair to form a distorted d(G x G) reverse Hoogsteen base-pair in the d(Gm5CGCGCG)2 homoduplex, and d(G x T) reverse wobble and d(G x C) reverse Watson-Crick base-pairs in the d(GCGCGCG)/d(TCGCGCG) and d(GCGCGCG)/d(CCGCGCG) heteroduplexes, respectively. Interestingly, only the d(G,T) and d(G x C) base-pairs were observed in the heteroduplexes, suggesting that both the d(G x T) reverse wobble and d(G x C) reverse Watson-Crick base-pairs are more stable in this crystal environment than the d(G x G) reverse Hoogsteen base-pair. To estimate the relative stability of the three types of reverse base-pairs, crystals were grown using various mixtures of sequences and their strand compositions analyzed by mass spectrometry. The d(G x C) reverse Watson-Crick base-pair was estimated to be more stable by approximately 1.5 kcal/mol and the d(G x T) reverse wobble base-pair more stable by approximately 0.5 kcal/mol than the d(G x G) reverse Hoogsteen base-pair. The step during crystallization responsible for discriminating between the strands in the crystal is highly cooperative, suggesting that it occurs during the initial nucleating event of crystal growth.

Identification and sequence analysis of contact sites between ribosomal proteins and rRNA in Escherichia coli 30 S subunits by a new approach using matrix-assisted laser desorption/ionization-mass spectrometry combined with N-terminal microsequencing

Cross-linked peptide-oligoribonucleotide complexes derived from distinct regions of the rRNA and individual ribosomal proteins of the 30 S ribosomal subunits from Escherichia coli were isolated and purified. Cross-linking sites at the amino acid and nucleotide level were determined by N-terminal amino acid sequence analysis in combination with matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS). MALDI-MS analysis performed subsequent to a partial alkaline hydrolysis of cross-linked peptide-oligoribonucleotide complexes allowed for the first time the cross-linked rRNA moiety to be sequenced by this technique. In this manner Lys-44 in S4 was determined to be cross-linked to the oligoribonucleotide at positions 1531-1542 on the 16 S RNA (whereby either U-1541 or A-1542 is the actual cross-link site), Lys-75 in S7 to positions 1374-1379 (C-1378 cross-linked), Met-114 in S7 to 1234-1241 (U-1240 cross-linked), Lys-55 in S8 to 651-654 (U-653 cross-linked), and Lys-29 in S17 to 629-633 (U-632 cross-linked). The novel approach applied here promises to be useful for similar studies on other known protein.RNA complexes.

Structural features of the binding site for ribosomal protein S8 in Escherichia coli 16S rRNA defined using NMR spectroscopy

Ribosomal protein S8 of Escherichia coli plays a key role in 30S ribosomal subunit assembly through its interaction with 16S rRNA. S8 also participates in the translational regulation of ribosomal protein expression through its interaction with spc operon mRNA. The binding site for protein S8 within the 16S rRNA encompasses nucleotides G588 to G604 and C634 to C651 and is composed of two base paired helical regions that flank a phylogenetically conserved core element containing nine residues. We have investigated the structure of the rRNA binding site for S8 both in the free state and in the presence of protein using NMR spectroscopy. The integrity of the two helical segments has been verified, and the presence of G597 x C643 and A596 x U644 base pairs within the conserved core, predicted from comparative analysis, have been confirmed. In addition, we have identified a base triple within the core that is composed of residues A595 x (A596 x U644). The NMR data suggest that S8-RNA interaction is accomplished without significant changes in the RNA. Nonetheless, S8 binding promotes formation of the U598 x A640 base pair and appears to stabilize the G597 x C643 and A596 x U644 base pairs.

Structural evidence for specific S8-RNA and S8-protein interactions within the 30S ribosomal subunit: ribosomal protein S8 from Bacillus stearothermophilus at 1.9 A resolution

BACKGROUND

Prokaryotic ribosomal protein S8 is an important RNA-binding protein that occupies a central position within the small ribosomal subunit. It interacts extensively with 16S rRNA and is crucial for the correct folding of the central domain of the rRNA. S8 also controls the synthesis of several ribosomal proteins by binding to mRNA. It binds specifically to very similar sites in the two RNA molecules.

RESULTS

S8 is divided into two tightly associated domains and contains three regions that are proposed to interact with other ribosomal components: two potential RNA-binding sites, and a hydrophobic patch that may interact with a complementary hydrophobic region of S5. The N-terminal domain fold is found in several proteins including two that bind double-stranded DNA.

CONCLUSIONS

These multiple RNA-binding sites are consistent with the role of S8 in organizing the central domain and agree with the latest models of the 16S RNA which show that the S8 location coincides with a region of complicated nucleic-acid structure. The presence in a wide variety of proteins of a region homologous to the N-terminal domain supports the idea that ribosomal proteins must represent some of the earliest protein molecules.

Determining the conformation of RNAs in solution. Application to a retroviral system: structure of the HIV-1 primer binding site region and effect of tRNA(3Lys) binding

RNAs play a crucial and central role in a large variety of biological functions obviously linked to the wide variety of structures that they can adopt. Understanding the function of RNAs thus requires the knowledge of their two- and three-dimensional structures. We describe in detail the way to access the secondary structure of RNAs, by combining sequence comparison, secondary structure prediction by computer and, mainly, experimental data obtained by probing with chemicals and ribonucleases. These approaches were used to investigate secondary structure of the region containing the primer binding site of HIV-1 genomic RNA either free or involved in the binary complex with the replication primer tRNA(3Lys).

The 16S ribosomal RNA mutation database (16SMDB)

The 16S ribosomal RNA mutation database (16SMDB) provides a list of mutated positions in 16S ribosomal RNA from Escherichia coli and the identity of each alteration. Information provided for each mutation includes: (i) a brief description of the phenotype(s) associated with each mutation; (ii) whether a mutant phenotype has been detected by in vivo or in vitro methods; (iii) relevant literature citations. The database is available via ftp and on the World Wide Web.