The S(MK) box is a new SAM-binding RNA for translational regulation of SAM synthetase

We have identified the S(MK) box as a conserved RNA motif in the 5' untranslated leader region of metK (SAM synthetase) genes in lactic acid bacteria, including Enterococcus, Streptococcus and Lactococcus species. This RNA element bound SAM in vitro, and binding of SAM caused an RNA structural rearrangement that resulted in sequestration of the Shine-Dalgarno (SD) sequence. Mutations that disrupted pairing between the SD region and a sequence complementary to the SD blocked SAM binding, whereas compensatory mutations that restored pairing restored SAM binding. The Enterococcus faecalis S(MK) box conferred translational repression of a lacZ reporter when cells were grown under conditions where SAM pools are elevated, and mutations that blocked SAM binding resulted in loss of repression, demonstrating that the S(MK) box is functional in vivo. The S(MK) box therefore represents a new SAM-binding riboswitch distinct from the previously identified S box RNAs.     

Discrimination between Closely Related Cellular Metabolites by the SAM-I Riboswitch

The SAM-I riboswitch is a cis-acting element of genetic control found in bacterial mRNAs that specifically binds S-adenosylmethionine (SAM). We previously determined the 2.9-Å X-ray crystal structure of the effector-binding domain of this RNA element, revealing details of RNA-ligand recognition. To improve this structure, variations were made to the RNA sequence to alter lattice contacts, resulting in a 0.5-Å improvement in crystallographic resolution and allowing for a more accurate refinement of the crystallographic model. The basis for SAM specificity was addressed by a structural analysis of the RNA complexed to S-adenosylhomocysteine (SAH) and sinefungin and by measuring the affinity of SAM and SAH for a series of mutants using isothermal titration calorimetry. These data illustrate the importance of two universally conserved base pairs in the RNA that form electrostatic interactions with the positively charged sulfonium group of SAM, thereby providing a basis for discrimination between SAM and SAH.     

An mRNA structure that controls gene expression by binding S-adenosylmethionine

Riboswitches are metabolite-binding RNA structures that serve as genetic control elements for certain messenger RNAs. These RNA switches have been identified in all three kingdoms of life and are typically responsible for the control of genes whose protein products are involved in the biosynthesis, transport or utilization of the target metabolite. Herein, we report that a highly conserved RNA domain found in bacteria serves as a riboswitch that responds to the coenzyme S-adenosylmethionine (SAM) with remarkably high affinity and specificity. SAM riboswitches undergo structural reorganization upon introduction of SAM, and these allosteric changes regulate the expression of 26 genes in Bacillus subtilis. This and related findings indicate that direct interaction between small metabolites and allosteric mRNAs is an important and widespread form of genetic regulation in bacteria.     

Common themes and differences in SAM recognition among SAM riboswitches

The recent discovery of short cis-acting RNA elements termed riboswitches has caused a paradigm shift in our understanding of genetic regulatory mechanisms. The three distinct superfamilies of S-adenosyl-l-methionine (SAM) riboswitches are the most commonly found riboswitch classes in nature. These RNAs represent three independent evolutionary solutions to achieve specific SAM recognition. This review summarizes research on 1) modes of gene regulatory mechanisms, 2) common themes and differences in ligand recognition, and 3) ligand-induced conformational dynamics among SAM riboswitch families. The body of work on the SAM riboswitch families constitutes a useful primer to the topic of gene regulatory RNAs as a whole.     

Structural basis for recognition of S-adenosylhomocysteine by riboswitches

S-adenosyl-(L)-homocysteine (SAH) riboswitches are regulatory elements found in bacterial mRNAs that up-regulate genes involved in the S-adenosyl-(L)-methionine (SAM) regeneration cycle. To understand the structural basis of SAH-dependent regulation by RNA, we have solved the structure of its metabolite-binding domain in complex with SAH. This structure reveals an unusual pseudoknot topology that creates a shallow groove on the surface of the RNA that binds SAH primarily through interactions with the adenine ring and methionine main chain atoms and discriminates against SAM through a steric mechanism. Chemical probing and calorimetric analysis indicate that the unliganded RNA can access bound-like conformations that are significantly stabilized by SAH to direct folding of the downstream regulatory switch. Strikingly, we find that metabolites bearing an adenine ring, including ATP, bind this aptamer with sufficiently high affinity such that normal intracellular concentrations of these compounds may influence regulation of the riboswitch.     

The yjdF riboswitch candidate regulates gene expression by binding diverse azaaromatic compounds

The yjdF motif RNA is an orphan riboswitch candidate that almost exclusively associates with the yjdF protein-coding gene in many bacteria. The function of the YjdF protein is unknown, which has made speculation regarding the natural ligand for this putative riboswitch unusually challenging. By using a structure-probing assay for ligand binding, we found that a surprisingly broad diversity of nitrogen-containing aromatic heterocycles, or “azaaromatics,” trigger near-identical changes in the structures adopted by representative yjdF motif RNAs. Regions of the RNA that undergo ligand-induced structural modulation reside primarily in portions of the putative aptamer region that are highly conserved in nucleotide sequence, as is typical for riboswitches. Some azaaromatic molecules are bound by the RNA with nanomolar dissociation constants, and a subset of these ligands activate riboswitch-mediated gene expression in cells. Furthermore, genetic elements most commonly adjacent to the yjdF motif RNA or to the yjdF protein-coding region are homologous to protein regulators implicated in mitigating the toxic effects of diverse phenolic acids or polycyclic compounds. Although the precise type of natural ligand sensed by yjdF motif RNAs remains unknown, our findings suggest that this riboswitch class might serve as part of a genetic response system to toxic or signaling compounds with chemical structures similar to azaaromatics.     

Evidence for a second class of S-adenosylmethionine riboswitches and other regulatory RNA motifs in alpha-proteobacteria

Background

Riboswitches are RNA elements in the 5' untranslated leaders of bacterial mRNAs that directly sense the levels of specific metabolites with a structurally conserved aptamer domain to regulate expression of downstream genes. Riboswitches are most common in the genomes of low GC Gram-positive bacteria (for example, Bacillus subtilis contains examples of all known riboswitches), and some riboswitch classes seem to be restricted to this group.

Results

We used comparative sequence analysis and structural probing to identify five RNA elements (serC, speF, suhB, ybhL, and metA) that reside in the intergenic regions of Agrobacterium tumefaciens and many other α-proteobacteria. One of these, the metA motif, is found upstream of methionine biosynthesis genes and binds S-adenosylmethionine (SAM). This natural aptamer most likely functions as a SAM riboswitch (SAM-II) with a consensus sequence and structure that is distinct from the class of SAM riboswitches (SAM-I) predominantly found in Gram-positive bacteria. The minimal functional SAM-II aptamer consists of fewer than 70 nucleotides, which form a single stem and a pseudoknot. Despite its simple architecture and lower affinity for SAM, the SAM-II aptamer strongly discriminates against related compounds.

Conclusion

SAM-II is the only metabolite-binding riboswitch class identified so far that is not found in Gram-positive bacteria, and its existence demonstrates that biological systems can use multiple RNA structures to sense a single chemical compound. The two SAM riboswitches might be 'RNA World' relics that were selectively retained in certain bacterial lineages or new motifs that have emerged since the divergence of the major bacterial groups.     

Identification of 22 candidate structured RNAs in bacteria using the CMfinder comparative genomics pipeline

We applied a computational pipeline based on comparative genomics to bacteria, and identified 22 novel candidate RNA motifs. We predicted six to be riboswitches, which are mRNA elements that regulate gene expression on binding a specific metabolite. In separate studies, we confirmed that two of these are novel riboswitches. Three other riboswitch candidates are upstream of either a putative transporter gene in the order Lactobacillales, citric acid cycle genes in Burkholderiales or molybdenum cofactor biosynthesis genes in several phyla. The remaining riboswitch candidate, the widespread Genes for the Environment, for Membranes and for Motility (GEMM) motif, is associated with genes important for natural competence in Vibrio cholerae and the use of metal ions as electron acceptors in Geobacter sulfurreducens. Among the other motifs, one has a genetic distribution similar to a previously published candidate riboswitch, ykkC/yxkD, but has a different structure. We identified possible non-coding RNAs in five phyla, and several additional cis-regulatory RNAs, including one in ε-proteobacteria (upstream of purD, involved in purine biosynthesis), and one in Cyanobacteria (within an ATP synthase operon). These candidate RNAs add to the growing list of RNA motifs involved in multiple cellular processes, and suggest that many additional RNAs remain to be discovered.     

Atomic-level insights into metabolite recognition and specificity of the SAM-II riboswitch

Although S-adenosylhomocysteine (SAH), a metabolic by-product of S-adenosylmethionine (SAM), differs from SAM only by a single methyl group and an overall positive charge, SAH binds the SAM-II riboswitch with more than 1000-fold less affinity than SAM. Using atomistic molecular dynamics simulations, we investigated the molecular basis of such high selectivity in ligand recognition by SAM-II riboswitch. The biosynthesis of SAM exclusively generates the (S,S) stereoisomer, and (S,S)-SAM can spontaneously convert to the (R,S) form. We, therefore, also examined the effects of (R,S)-SAM binding to SAM-II and its potential biological function. We find that the unfavorable loss in entropy in SAM-II binding is greater for (S,S)- and (R,S)-SAM than SAH, which is compensated by stabilizing electrostatic interactions with the riboswitch. The positively charged sulfonium moiety on SAM acts as the crucial anchor point responsible for the formation of key ionic interactions as it fits favorably in the negatively charged binding pocket. In contrast, SAH, with its lone pair of electrons on the sulfur, experiences repulsion in the binding pocket of SAM-II and is enthalpically destabilized. In the presence of SAH, similar to the unbound riboswitch, the pseudoknot structure of SAM-II is not completely formed, thus exposing the Shine-Dalgarno sequence. Unlike SAM, this may further facilitate ribosomal assembly and translation initiation. Our analysis of the conformational ensemble sampled by SAM-II in the absence of ligands and when bound to SAM or SAH reveals that ligand binding follows a combination of conformational selection and induced-fit mechanisms.     

Riboswitches that sense S-adenosylmethionine and S-adenosylhomocysteine.

Numerous riboswitches have been discovered that specifically recognize metabolites and modulate gene expression. Each riboswitch class is defined either by the consensus sequence and structural features of its metabolite-binding aptamer domain, or by the distinct metabolite that the aptamer recognizes. Several distinct classes of riboswitches that respond to S-adenosylmethionine (SAM or AdoMet) have been discovered. Representatives of these classes have been shown to strongly discriminate against S-adenosylhomocystenine (SAH or AdoHcy), which is the metabolic byproduct produced when SAM is used as a cofactor for methylation reactions. However, a distinct class of riboswitches that selectively binds SAH, and strongly discriminates against SAM, also has been discovered. Herein we compare the features of SAM and SAH riboswitches, which help showcase the enormous structural diversity that RNA can harness to form precision genetic switches for compounds that are critical for fundamental metabolic processes.     

Prediction of signal recognition particle RNA genes

We describe a method for prediction of genes that encode the RNA component of the signal recognition particle (SRP). A heuristic search for the strongly conserved helix 8 motif of SRP RNA is combined with covariance models that are based on previously known SRP RNA sequences. By screening available genomic sequences we have identified a large number of novel SRP RNA genes and we can account for at least one gene in every genome that has been completely sequenced. Novel bacterial RNAs include that of Thermotoga maritima, which, unlike all other non-gram-positive eubacteria, is predicted to have an Alu domain. We have also found the RNAs of Lactococcus lactis and Staphylococcus to have an unusual UGAC tetraloop in helix 8 instead of the normal GNRA sequence. An investigation of yeast RNAs reveals conserved sequence elements of the Alu domain that aid in the analysis of these RNAs. Analysis of the human genome reveals only two likely genes, both on chromosome 14. Our method for SRP RNA gene prediction is the first convenient tool for this task and should be useful in genome annotation.     

Confirmation of a second natural preQ1 aptamer class in Streptococcaceae bacteria

Bioinformatics searches of eubacterial genomes have yielded many riboswitch candidates where the identity of the ligand is not immediately obvious on examination of associated genes. One of these motifs is found exclusively in the family Streptococcaceae within the 5' untranslated regions (UTRs) of genes encoding the hypothetical membrane protein classified as COG4708 or DUF988. While the function of this protein class is unproven, a riboswitch binding the queuosine biosynthetic intermediate pre-queuosine(1) (preQ(1)) has been identified in the 5' UTR of homologous genes in many Firmicute species of bacteria outside of Streptococcaceae. Here we show that a representative of the COG4708 RNA motif from Streptococcus pneumoniae R6 also binds preQ(1). Furthermore, representatives of this RNA have structural and molecular recognition characteristics that are distinct from those of the previously described preQ(1) riboswitch class. PreQ(1) is the second metabolite for which two or more distinct classes of natural aptamers exist, indicating that natural aptamers utilizing different structures to bind the same metabolite may be more common than is currently known. Additionally, the association of preQ(1) binding RNAs with most genes encoding proteins classified as COG4708 strongly suggests that these proteins function as transporters for preQ(1) or another queuosine biosynthetic intermediate.     

Folding of the SAM aptamer is determined by the formation of a K-turn-dependent pseudoknot

The S-adenosylmethionine (SAM) riboswitch is one of the most recurrent riboswitches found in bacteria and has three known different natural aptamers. The Bacillus subtilis yitJ SAM riboswitch aptamer is organized around a four-way junction which is characterized by the presence of a pseudoknot and a K-turn motif. By replacing the adenine involved in a Watson-Crick base pair at position 138 in the core region of the aptamer with the fluorescent analogue 2-aminopurine (2AP), we show that the ligand-induced reorganization of the aptamer strongly attenuates 2AP fluorescence. The fluorescence quenching process is specific to SAM on the basis of the observation that the structural analogue S-adenosylhomocysteine does not promote a similar effect. We find that the pseudoknot is important for the reorganization of the core domain and that the K-turn motif also has a marked influence on the core domain reorganization, most probably through its important role in pseudoknot formation. Finally, we show that SAM riboswitch ligand binding is facilitated by the L7Ae K-turn binding protein, which suggests that K-turn motifs may be protein anchor sites used by riboswitches to promote RNA folding.     

Coenzyme B12 riboswitches are widespread genetic control elements in prokaryotes

Recent studies have begun to reveal that numerous fundamental metabolic pathways in bacteria are regulated by riboswitches residing within certain messenger RNAs. These riboswitches selectively bind metabolites and modulate gene expression in response to changing ligand concentrations. Previously, we provided evidence that the btuB mRNAs of Escherichia coli and Salmonella typhimurium each carry a coenzyme B12-dependent riboswitch that causes repressed translation of the encoded cobalamin-transport protein at elevated coenzyme concentrations. Herein, we use a phylogenetic analysis to define a consensus sequence and secondary structure model for the ligand- binding domain of this riboswitch class. RNA structures that conform to this model are widespread in both Gram-positive and Gram-negative organisms. In addition, we find that the 5'-untranslated region (5'-UTR) of the cobalamin biosynthesis (cob) operon of S.typhimurium carries an RNA motif that matches this consensus sequence. Biochemical and genetic characterization of this motif confirms that the RNA directly binds coenzyme B12, and that it likely serves as a genetic control element for regulating expression of the 25-gene operon for cobalamin production in this pathogen.