Guanidine-sensing riboswitches: How do they work and what do they regulate?

Atomistic Simulations Reveal Crucial Role of Metal Ions for Ligand Binding in Guanidine-I Riboswitch

Riboswitches are structured ribonucleic acid (RNA) segments that act as specific sensors for small molecules in bacterial metabolism. Due to the flexible nature of these highly charged macromolecules, molecular dynamics simulations are instrumental to investigating the mechanistic details of their regulatory function. In the present study, the guanidine-I riboswitch serves as an example of how atomistic simulations can shed light on the effect of ions on the structure and dynamics of RNA and on ligand binding. Relying on two orthologous crystal structures from different bacterial species, it is demonstrated how the ion setup crucially determines whether the simulation yields meaningful insights into the conformational stability of the RNA, functionally relevant residues and RNA-ligand interactions. The ion setup in this context includes diffuse ions in solution and bound ions associated directly with the RNA, in particular a triad of 2 Mg2+ ions and a K+ ion in close proximity to the guanidinium binding site. A detailed investigation of the binding pocket reveals that the K+ from the ion triad plays a decisive role in stabilizing the ligand binding by stabilizing important localized interactions, which in turn contribute to the overall shape of the folded state of the RNA.

Excited States of apo-Guanidine-III Riboswitch Contribute to Guanidinium Binding through Both Conformational and Induced-Fit Mechanisms

Riboswitches are mRNA segments that regulate gene expression through conformational changes driven by their cognate ligand binding. The ykkC motif forms a riboswitch class that selectively senses a guanidinium ion (Gdm+) and regulates the downstream expression of proteins which aid in the efflux of excess Gdm+ from the cells. The aptamer domain (AD) of the guanidine-III riboswitch forms an H-type pseudoknot with a triple helical domain that binds a Gdm+. We studied the binding of Gdm+ to the AD of the guanidine (ykkC)-III riboswitch using computer simulations to probe the specificity of the riboswitch to Gdm+ binding. We show that Gdm+ binding is a fast process occurring on the nanosecond time scale, with minimal conformational changes to the AD. Using machine learning and Markov-state models, we identified the excited conformational states of the AD, which have a high Gdm+ binding propensity, making the Gdm+ binding landscape complex exhibiting both conformational selection and induced-fit mechanisms. The proposed apo-AD excited states and their role in the ligand-sensing mechanism are amenable to experimental verification. Further, targeting these excited-state conformations in discovering new antibiotics can be explored.

Characterization of Structure and Dynamics of the Guanidine-II Riboswitch from Escherichia coli by NMR Spectroscopy and Small-Angle X-ray Scattering (SAXS)

Riboswitches are regulatory RNA elements that undergo functionally important allosteric conformational switching upon binding of specific ligands. The here investigated guanidine-II riboswitch binds the small cation, guanidinium, and forms a kissing loop-loop interaction between its P1 and P2 hairpins. We investigated the structural changes to support previous studies regarding the binding mechanism. Using NMR spectroscopy, we confirmed the structure as observed in crystal structures and we characterized the kissing loop interaction upon addition of Mg2+ and ligand for the riboswitch aptamer from Escherichia coli. We further investigated closely related mutant constructs providing further insight into functional differences between the two (different) hairpins P1 and P2. Formation of intermolecular interactions were probed by small-angle X-ray scattering (SAXS) and NMR DOSY data. All data are consistent and show the formation of oligomeric states of the riboswitch induced by Mg2+ and ligand binding.

The Biochemical Landscape of Riboswitch Ligands

More than 55 distinct classes of riboswitches that respond to small metabolites or elemental ions have been experimentally validated to date. The ligands sensed by these riboswitches are biased in favor of fundamental compounds or ions that are likely to have been relevant to ancient forms of life, including those that might have populated the "RNA World", which is a proposed biochemical era that predates the evolutionary emergence of DNA and proteins. In the following text, I discuss the various types of ligands sensed by some of the most common riboswitches present in modern bacterial cells and consider implications for ancient biological processes centered on the proven capabilities of these RNA-based sensors. Although most major biochemical aspects of metabolism are represented by known riboswitch classes, there are striking sensory gaps in some key areas. These gaps could reveal weaknesses in the performance capabilities of RNA that might have hampered RNA World evolution, or these could highlight opportunities to discover additional riboswitch classes that sense essential metabolites.

Molecular Dynamics Simulations of the Aptamer Domain of Guanidinium Ion Binding Riboswitch ykkC-III: Structural Insights into the Discrimination of Cognate and Alternate Ligands

Guanidinium ion is a toxic cellular metabolite. The ykkC-III riboswitch, an mRNA stretch, regulates the gene expression by undergoing a conformational change in response to the binding of a free guanidinium ion and thereby plays a potentially important role in alleviating guanidinium toxicity in cells. An experimental crystal structure of the guanidinium-bound aptamer domain of the riboswitch from Thermobifida Fusca revealed the overall RNA architecture and mapped the specific noncovalent interactions that stabilize the ligand within the binding pocket aptamer. However, details of how the aptamer domain discriminates the cognate ligand from its closest structurally analogous physiological metabolites (arginine and urea), and how the binding of cognate ligand arrays information from the aptamer domain to the expression platform for regulating the gene expression, are not well understood. To fill this void, we perform a cumulative of 2 μs all-atom explicit-solvent molecular dynamics (MD) simulations on the full aptamer domain, augmented with quantum-chemical calculations on the ligand-binding pocket, to compare the structural and dynamical details of the guanidinium-bound state with the arginine or urea bound states, as well as the unbound (open) state. Analysis of the ligand-binding pocket reveals that due to unfavorable interactions with the binding-pocket residues, urea cannot bind the aptamer domain and thereby cannot alter the gene expression. Although interaction of the guanidyl moiety of arginine within the binding pocket is either comparable or stronger than the guanidinium ion, additional non-native hydrogen-bonding networks, as well as differences in the dynamical details of the arginine-bound state, explain why arginine cannot transmit the information from the aptamer domain to the expression platform. Based on our simulations, we propose a mechanism of how the aptamer domain communicates with the expression platform. Overall, our work provides interesting insights into the ligand recognition by a specific class of riboswitches and may hopefully inspire future studies to further understand the gene regulation by riboswitches.

Riboregulation in bacteria: From general principles to novel mechanisms of the trp attenuator and its sRNA and peptide products

Gene expression strategies ensuring bacterial survival and competitiveness rely on cis- and trans-acting RNA-regulators (riboregulators). Among the cis-acting riboregulators are transcriptional and translational attenuators, and antisense RNAs (asRNAs). The trans-acting riboregulators are small RNAs (sRNAs) that bind proteins or base pairs with other RNAs. This classification is artificial since some regulatory RNAs act both in cis and in trans, or function in addition as small mRNAs. A prominent example is the archetypical, ribosome-dependent attenuator of tryptophan (Trp) biosynthesis genes. It responds by transcription attenuation to two signals, Trp availability and inhibition of translation, and gives rise to two trans-acting products, the attenuator sRNA rnTrpL and the leader peptide peTrpL. In Escherichia coli, rnTrpL links Trp availability to initiation of chromosome replication and in Sinorhizobium meliloti, it coordinates regulation of split tryptophan biosynthesis operons. Furthermore, in S. meliloti, peTrpL is involved in mRNA destabilization in response to antibiotic exposure. It forms two types of asRNA-containing, antibiotic-dependent ribonucleoprotein complexes (ARNPs), one of them changing the target specificity of rnTrpL. The posttranscriptional role of peTrpL indicates two emerging paradigms: (1) sRNA reprograming by small molecules and (2) direct involvement of antibiotics in regulatory RNPs. They broaden our view on RNA-based mechanisms and may inspire new approaches for studying, detecting, and using antibacterial compounds. This article is categorized under: RNA Interactions with Proteins and Other Molecules > Small Molecule-RNA Interactions RNA Interactions with Proteins and Other Molecules > RNA-Protein Complexes Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Genomic epidemiology of rifampicin ADP-ribosyltransferase (Arr) in the Bacteria domain

Arr is an ADP-ribosyltransferase enzyme primarily reported in association with rifamycin resistance, which has been used to treat tuberculosis in addition to Gram-positive infections and, recently, pan-resistant Gram-negative bacteria. The arr gene was initially identified on the Mycolicibacterium smegmatis chromosome and later on Proteobacteria plasmids. This scenario raised concerns on the distribution and spread of arr, considering the Bacteria domain. Based on 198,082 bacterial genomes/metagenomes, we performed in silico analysis, including phylogenetic reconstruction of Arr in different genomic contexts. Besides, new arr alleles were evaluated by in vitro analysis to assess their association with rifampin resistance phenotype. The arr gene was prevalent in thousands of chromosomes and in hundreds of plasmids from environmental and clinical bacteria, mainly from the phyla Actinobacteria, Proteobacteria, Firmicutes, and Bacteroidetes. Furthermore, this gene was identified in other and new genomic contexts. Interestingly, Arr sequences associated with rifampin resistance were distributed across all phylogeny, indicating that, despite the diversity, their association with rifampin resistance phenotype were maintained. In fact, we found that the key residues were highly conserved. In addition, other analyzes have raised evidence of another Arr function, which is related to guanidine metabolism. Finally, this scenario as a whole also suggested the Actinobacteria phylum as a potential ancestral source of arr within the Bacteria domain.

Guanidine-II aptamer conformations and ligand binding modes through the lens of molecular simulation

Regulation of gene expression via riboswitches is a widespread mechanism in bacteria. Here, we investigate ligand binding of a member of the guanidine sensing riboswitch family, the guanidine-II riboswitch (Gd-II). It consists of two stem–loops forming a dimer upon ligand binding. Using extensive molecular dynamics simulations we have identified conformational states corresponding to ligand-bound and unbound states in a monomeric stem–loop of Gd-II and studied the selectivity of this binding. To characterize these states and ligand-dependent conformational changes we applied a combination of dimensionality reduction, clustering, and feature selection methods. In absence of a ligand, the shape of the binding pocket alternates between the conformation observed in presence of guanidinium and a collapsed conformation, which is associated with a deformation of the dimerization interface. Furthermore, the structural features responsible for the ability to discriminate against closely related analogs of guanidine are resolved. Based on these insights, we propose a mechanism that couples ligand binding to aptamer dimerization in the Gd-II system, demonstrating the value of computational methods in the field of nucleic acids research.

Large Stokes shift fluorescence activation in an RNA aptamer by intermolecular proton transfer to guanine

Fluorogenic RNA aptamers are synthetic functional RNAs that specifically bind and activate conditional fluorophores. The Chili RNA aptamer mimics large Stokes shift fluorescent proteins and exhibits high affinity for 3,5-dimethoxy-4-hydroxybenzylidene imidazolone (DMHBI) derivatives to elicit green or red fluorescence emission. Here, we elucidate the structural and mechanistic basis of fluorescence activation by crystallography and time-resolved optical spectroscopy. Two co-crystal structures of the Chili RNA with positively charged DMHBO⁺ and DMHBI⁺ ligands revealed a G-quadruplex and a trans-sugar-sugar edge G:G base pair that immobilize the ligand by π-π stacking. A Watson-Crick G:C base pair in the fluorophore binding site establishes a short hydrogen bond between the N7 of guanine and the phenolic OH of the ligand. Ultrafast excited state proton transfer (ESPT) from the neutral chromophore to the RNA was found with a time constant of 130 fs and revealed the mode of action of the large Stokes shift fluorogenic RNA aptamer.

Discovery and characterization of a fourth class of guanidine riboswitches

Riboswitches are RNAs that specifically sense a small molecule and regulate genes accordingly. The recent discovery of guanidine-binding riboswitches revealed the biological significance of this compound, and uncovered genes related to its biology. For example, certain sugE genes encode guanidine exporters and are activated by the riboswitches to reduce toxic levels of guanidine in the cell. In order to study guanidine biology and riboswitches, we applied a bioinformatics strategy for discovering additional guanidine riboswitches by searching for new candidate motifs associated with sugE genes. Based on in vitro and in vivo experiments, we determined that one of our six best candidates is a new structural class of guanidine riboswitches. The expression of a genetic reporter was induced 80-fold in response to addition of 5 mM guanidine in Staphylococcus aureus. This new class, called the guanidine-IV riboswitch, reveals additional guanidine-associated protein domains that are extremely rarely or never associated with previously established guanidine riboswitches. Among these protein domains are two transporter families that are structurally distinct from SugE, and could represent novel types of guanidine exporters. These results establish a new metabolite-binding RNA, further validate a bioinformatics method for finding riboswitches and suggest substrate specificities for as-yet uncharacterized transporter proteins.

The structural basis of promiscuity in small multidrug resistance transporters

By providing broad resistance to environmental biocides, transporters from the small multidrug resistance (SMR) family drive the spread of multidrug resistance cassettes among bacterial populations. A fundamental understanding of substrate selectivity by SMR transporters is needed to identify the types of selective pressures that contribute to this process. Using solid-supported membrane electrophysiology, we find that promiscuous transport of hydrophobic substituted cations is a general feature of SMR transporters. To understand the molecular basis for promiscuity, we solved X-ray crystal structures of a SMR transporter Gdx-Clo in complex with substrates to a maximum resolution of 2.3 Å. These structures confirm the family’s extremely rare dual topology architecture and reveal a cleft between two helices that provides accommodation in the membrane for the hydrophobic substituents of transported drug-like cations.

Gdx-Clo is a bacterial transporter from the small multidrug resistance (SMR) family. Here, the authors use solid supported membrane electrophysiology to characterize Gdx-Clo functionally and report crystal structures of Gdx-Clo which confirm the dual topology architecture and offer insight into substrate binding and transport mechanism.

Biochemical Validation of a Fourth Guanidine Riboswitch Class in Bacteria

An intriguing consequence of ongoing riboswitch discovery efforts is the occasional identification of metabolic or toxicity response pathways for unusual ligands. Recently, we reported the experimental validation of three distinct bacterial riboswitch classes that regulate gene expression in response to the selective binding of a guanidinium ion. These riboswitch classes, called guanidine-I, -II, and -III, regulate numerous genes whose protein products include previously misannotated guanidine exporters and enzymes that degrade guanidine via an initial carboxylation reaction. Guanidine is now recognized as the primal substrate of many multidrug efflux pumps that are important for bacterial resistance to certain antibiotics. Guanidine carboxylase enzymes had long been annotated as urea carboxylase enzymes but are now understood to participate in guanidine degradation. Herein, we report the existence of a fourth riboswitch class for this ligand, called guanidine-IV. Members of this class use a novel aptamer to selectively bind guanidine and use an unusual expression platform arrangement that is predicted to activate gene expression when ligand is present. The wide distribution of this abundant riboswitch class, coupled with the striking diversity of other guanidine-sensing RNAs, demonstrates that many bacterial species maintain sophisticated sensory and genetic mechanisms to avoid guanidine toxicity. This finding further highlights the mystery regarding the natural source of this nitrogen-rich chemical moiety.

Riboswitch-Associated Guanidinium-Selective Efflux Pumps Frequently Transmitted on Proteobacterial Plasmids Increase Escherichia coli Biofilm Tolerance to Disinfectants

Members of the small multidrug resistance (SMR) efflux pump family known as SugE (recently renamed Gdx) are known for their narrow substrate selectivity to small guanidinium (Gdm+) compounds and disinfectant quaternary ammonium compounds (QACs). Gdx members have been identified on multidrug resistance plasmids in Gram-negative bacilli, but their functional role remains unclear, as few have been characterized. Here, we conducted a survey of sequenced proteobacterial plasmids that encoded one or more SugE/Gdx sequences in an effort to (i) identify the most frequently represented Gdx member(s) on these plasmids and their sequence diversity, (ii) verify if Gdx sequences possess a Gdm+ riboswitch that regulates their translation similarly to chromosomally encoded Gdx members, and (iii) determine the antimicrobial susceptibility profile of the most predominate Gdx member to various QACs and antibiotics in Escherichia coli strains BW25113 and KAM32. The results of this study determined 14 unique SugE sequences, but only one Gdx sequence, annotated as "SugE(p)," predominated among the >140 plasmids we surveyed. Enterobacterales plasmids carrying sugE(p) possessed a guanidine II riboswitch similar to the upstream region of E. coligdx Cloning and expression of sugE(p), gdx, and emrE sequences into a low-copy-number expression vector (pMS119EH) revealed significant increases in QAC resistance to a limited range of detergent-like QACs only when gdx and sugE(p) transformants were grown as biofilms. These findings suggest that sugE(p) presence on proteobacterial plasmids may be driven by species that frequently encounter Gdm+ and QAC exposure.IMPORTANCE This study characterized the function of antimicrobial-resistant phenotypes attributed to plasmid-encoded guanidinium-selective small multidrug resistance (Gdm/SugE) efflux pumps. These sequences are frequently monitored as biocide resistance markers in antimicrobial resistance surveillance studies. Our findings reveal that enterobacterial gdm sequences transmitted on plasmids possess a guanidine II riboswitch, which restricts transcript translation in the presence of guanidinium. Cloning and overexpression of this gdm sequence revealed that it confers higher resistance to quaternary ammonium compound (QAC) disinfectants (which possess guanidium moieties) when grown as biofilms. Since biofilms are commonly eradicated with QAC-containing compounds, the presence of this gene on plasmids and its biofilm-specific resistance are a growing concern for clinical and food safety prevention measures.

Do the P1 and P2 hairpins of the Guanidine-II riboswitch interact?

Riboswitches regulate genes by adopting different structures in responds to metabolite binding. The guanidine-II riboswitch is the smallest representative of the ykkC class with the mechanism of its function being centred on the idea that its two stem loops P1 and P2 form a kissing hairpin interaction upon binding of guanidinium (Gdm⁺). This mechanism is based on in-line probing experiments with the full-length riboswitch and crystal structures of the truncated stem loops P1 and P2. However, the crystal structures reveal only the formation of the homodimers P1 | P1 and P2 | P2 but not of the proposed heterodimer P1 | P2. Here, site-directed spin labeling (SDSL) in combination with Pulsed Electron–Electron Double Resonance (PELDOR or DEER) is used to study their structures in solution and how they change upon binding of Gdm⁺. It is found that both hairpins adopt different structures in solution and that binding of Gdm⁺ does indeed lead to the formation of the heterodimer but alongside the homodimers in a statistical 1:2:1 fashion. These results do thus support the proposed switching mechanism.

Fundamental studies of functional nucleic acids: aptamers, riboswitches, ribozymes and DNAzymes

This review aims at juxtaposing common versus distinct structural and functional strategies that are applied by aptamers, riboswitches, and ribozymes/DNAzymes. Focusing on recently discovered systems, we begin our analysis with small-molecule binding aptamers, with emphasis on in vitro-selected fluorogenic RNA aptamers and their different modes of ligand binding and fluorescence activation. Fundamental insights are much needed to advance RNA imaging probes for detection of exo- and endogenous RNA and for RNA process tracking. Secondly, we discuss the latest gene expression-regulating mRNA riboswitches that respond to the alarmone ppGpp, to PRPP, to NAD⁺, to adenosine and cytidine diphosphates, and to precursors of thiamine biosynthesis (HMP-PP), and we outline new subclasses of SAM and tetrahydrofolate-binding RNA regulators. Many riboswitches bind protein enzyme cofactors that, in principle, can catalyse a chemical reaction. For RNA, however, only one system (glmS ribozyme) has been identified in Nature thus far that utilizes a small molecule - glucosamine-6-phosphate - to participate directly in reaction catalysis (phosphodiester cleavage). We wonder why that is the case and what is to be done to reveal such likely existing cellular activities that could be more diverse than currently imagined. Thirdly, this brings us to the four latest small nucleolytic ribozymes termed twister, twister-sister, pistol, and hatchet as well as to in vitro selected DNA and RNA enzymes that promote new chemistry, mainly by exploiting their ability for RNA labelling and nucleoside modification recognition. Enormous progress in understanding the strategies of nucleic acids catalysts has been made by providing thorough structural fundaments (e.g. first structure of a DNAzyme, structures of ribozyme transition state mimics) in combination with functional assays and atomic mutagenesis.

Imaginary Ribozymes

Biocatalysis is dominated by protein enzymes, and only a few classes of ribozymes are known to contribute to the task of promoting biochemical transformations. The RNA World theory encompasses the notion that earlier forms of life made use of a much greater diversity of ribozymes and other functional RNAs to guide complex metabolic states long before proteins had emerged in evolution. In recent years, the discoveries of various classes of ribozymes, riboswitches, and other noncoding RNAs in bacteria have provided additional support for the hypothesis that RNA molecules indeed have the catalytic competence to promote diverse chemical reactions without the aid of protein enzymes. Herein, some of the most striking observations made from examinations of natural riboswitches that bind small ligands are highlighted and used as a basis to imagine the characteristics and functions of long-extinct ribozymes from the RNA World.

Structure of an RNA aptamer in complex with the fluorophore tetramethylrhodamine

RNA aptamers-artificially created RNAs with high affinity and selectivity for their target ligand generated from random sequence pools-are versatile tools in the fields of biotechnology and medicine. On a more fundamental level, they also further our general understanding of RNA-ligand interactions e. g. in regard to the relationship between structural complexity and ligand affinity and specificity, RNA structure and RNA folding. Detailed structural knowledge on a wide range of aptamer-ligand complexes is required to further our understanding of RNA-ligand interactions. Here, we present the atomic resolution structure of an RNA-aptamer binding to the fluorescent xanthene dye tetramethylrhodamine. The high resolution structure, solved by NMR-spectroscopy in solution, reveals binding features both common and different from the binding mode of other aptamers with affinity for ligands carrying planar aromatic ring systems such as the malachite green aptamer which binds to the tetramethylrhodamine related dye malachite green or the flavin mononucleotide aptamer.