An energetically beneficial leader-linker interaction abolishes ligand-binding cooperativity in glycine riboswitches

Cryo-EM reveals remodeling of a tandem riboswitch at 2.9 Å resolution

Riboswitches are non-coding RNA sequences that control cellular processes through ligand binding. Conformational heterogeneity is fundamental to riboswitch functionality, yet this same attribute makes structural characterization of these mRNA elements challenging. Here, we use a combination of molecular dynamics and cryo-electron microscopy to expound the flexible nature of the glycine riboswitch tandem aptamers and characterize diMerent structural populations. We find that Mg2+ partially stabilizes the fully folded state, resulting in one-third of the particles adopting a unique "walking man" conformation, consisting of a rigidified core and two dynamic helices, and two-thirds adopting distinct, partially folded states. Glycine interactions double the relative population of fully folded particles by stabilizing a conserved inter-aptamer Hoogsteen base pair, enabling our capture of a 2.9 Å structure for this RNA-only system. The population data show that glycine and Mg2+ operate synergistically: glycine enhances Mg2+ occupancy, while Mg2+ drives glycine specificity. Our findings indicate that cryo-electron microscopy oMers a promising avenue to characterize RNA folding ensembles.

Linker-Mediated Inactivation of the SAM-II Domain in the Tandem SAM-II/SAM-V Riboswitch

Tandem SAM-II/SAM-V riboswitch belongs to a class of riboswitches found in the marine bacterium 'Candidatus Pelagibacter ubique'. Previous studies have demonstrated that these riboswitches have the potential for digital modulation of gene expression at both the transcriptional and translational levels. In this study, we investigate the conformational changes in the tandem SAM-II/SAM-V riboswitch binding to S-adenosylmethionine (SAM) using selective 2'-hydroxyl acylation analyzed by the primer extension (SHAPE) assay, small-angle X-ray scattering (SAXS), and oligos depressing probing. Our findings reveal that the linker between SAM-II/SAM-V aptamers blocks the SAM response of the SAM-II domain. This result proposes a new mechanism for gene expression regulation, where the ligand-binding functions of tandem riboswitches can be selectively masked or released through a linker.

Bioinformatics and Genomic Analyses of the Suitability of Eight Riboswitches for Antibacterial Drug Targets

Antibiotic resistance (AR) is an acute problem that results in prolonged and debilitating illnesses. AR mortality worldwide is growing and causes a pressing need to research novel mechanisms of action and untested target molecules. This article presents in silico analyses of eight bacterial riboswitches for their suitability for antibacterial drug targets. Most bacterial riboswitches are located in the 5′-untranslated region of messenger RNAs, act as allosteric cis-acting gene control elements, and have not been found in humans before. Sensing metabolites, the riboswitches regulate the synthesis of vital cellular metabolites in various pathogenic bacteria. The analyses performed in this article represent a complete and informative genome-wide bioinformatics analysis of the adequacy of eight riboswitches as antibacterial drug targets in different pathogenic bacteria based on four criteria. Due to the ability of the riboswitch to control biosynthetic pathways and transport proteins of essential metabolites and the presence/absence of alternative biosynthetic pathways, we classified them into four groups based on their suitability for use as antibacterial drug targets guided by our in silico analyses. We concluded that some of them are promising targets for antibacterial drug discovery, such as the PreQ1, MoCo RNA, cyclic-di-GMP I, and cyclic-di-GMP II riboswitches.

Development and characterization of a glycine biosensor system for fine-tuned metabolic regulation in Escherichia coli

Background

In vivo biosensors have a wide range of applications, ranging from the detection of metabolites to the regulation of metabolic networks, providing versatile tools for synthetic biology and metabolic engineering. However, in view of the vast array of metabolite molecules, the existing number and performance of biosensors is far from sufficient, limiting their potential applications in metabolic engineering. Therefore, we developed the synthetic glycine-ON and -OFF riboswitches for metabolic regulation and directed evolution of enzyme in Escherichia coli.

Results

The results showed that a synthetic glycine-OFF riboswitch (glyOFF6) and an increased-detection-range synthetic glycine-ON riboswitch (glyON14) were successfully screened from a library based on the Bacillus subtilis glycine riboswitch using fluorescence-activated cell sorting (FACS) and tetA-based dual genetic selection. The two synthetic glycine riboswitches were successfully used in tunable regulation of lactate synthesis, dynamic regulation of serine synthesis and directed evolution of alanine-glyoxylate aminotransferase in Escherichia coli, respectively. Mutants AGXT22 and AGXT26 of alanine-glyoxylate aminotransferase with an increase of 58% and 73% enzyme activity were obtained by using a high-throughput screening platform based on the synthetic glycine-OFF riboswitch, and successfully used to increase the 5-aminolevulinic acid yield of engineered Escherichia coli.

Conclusions

A synthetic glycine-OFF riboswitch and an increased-detection-range synthetic glycine-ON riboswitch were successfully designed and screened. The developed riboswitches showed broad application in tunable regulation, dynamic regulation and directed evolution of enzyme in E. coli.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12934-022-01779-4.

Architectures and complex functions of tandem riboswitches

Riboswitch architectures that involve the binding of a single ligand to a single RNA aptamer domain result in ordinary dose-response curves that require approximately a 100-fold change in ligand concentration to cover nearly the full dynamic range for gene regulation. However, by using multiple riboswitches or aptamer domains in tandem, these ligand-sensing structures can produce additional, complex gene control outcomes. In the current study, we have computationally searched for tandem riboswitch architectures in bacteria to provide a more complete understanding of the diverse biological and biochemical functions of gene control elements that are made exclusively of RNA. Numerous different arrangements of tandem homologous riboswitch architectures are exploited by bacteria to create more 'digital' gene control devices, which operate over a narrower ligand concentration range. Also, two heterologous riboswitch aptamers are sometimes employed to create two-input Boolean logic gates with various types of genetic outputs. These findings illustrate the sophisticated genetic decisions that can be made by using molecular sensors and switches based only on RNA.

Riboswitch Mechanisms: New Tricks for an Old Dog

Discovered almost twenty years ago, riboswitches turned out to be one of the most common regulatory systems in bacteria, with representatives found in eukaryotes and archaea. Unlike many other regulatory elements, riboswitches are entirely composed of RNA and capable of modulating expression of genes by direct binding of small cellular molecules. While bacterial riboswitches had been initially thought to control production of enzymes and transporters associated with small organic molecules via feedback regulatory circuits, later findings identified riboswitches directing expression of a wide range of genes and responding to various classes of molecules, including ions, signaling molecules, and others. The 5'-untranslated mRNA regions host a vast majority of riboswitches, which modulate transcription or translation of downstream genes through conformational rearrangements in the ligand-sensing domains and adjacent expression-controlling platforms. Over years, the repertoire of regulatory mechanisms employed by riboswitches has greatly expanded; most recent studies have highlighted the importance of alternative mechanisms, such as RNA degradation, for the riboswitch-mediated genetic circuits. This review discusses the plethora of bacterial riboswitch mechanisms and illustrates how riboswitches utilize different features and approaches to elicit various regulatory responses.

Regulation of Glycine Cleavage and Detoxification by a Highly Conserved Glycine Riboswitch in Burkholderia spp

The glycine riboswitch is a known regulatory element that is unique in having two aptamers that are joined by a linker region. In this study, we investigated a glycine riboswitch located in the 5' untranslated region of a glycine cleavage system homolog (gcvTHP) in Burkholderia spp. Structure prediction using the sequence generated a model with a glycine binding pocket composed of base-triple interactions (G62-A64-A86 and G65-U84-C85) that are supported by A/G minor interactions (A17-C60-G88 and G16-C61-G87, respectively) and two ribose-zipper motifs (C11-G12 interacting with A248-A247 and C153-U154 interacting with A79-A78) which had not been previously reported. The capacity of the riboswitch to bind to glycine was experimentally validated by native gel assays and the crucial role of interactions that make up the glycine binding pocket were proven by mutations of A17U and G16C which resulted in conformational differences that may lead to dysfunction. Using glycine supplemented minimal media, we were able to prove that the expression of the gcvTHP genes found downstream of the riboswitch responded to the glycine concentrations introduced thus confirming the role of this highly conserved Burkholderia riboswitch and its associated genes as a putative glycine detoxification system in Burkholderia spp.

Cooperativity and Allostery in RNA Systems

Allostery is among the most basic biological principles employed by biological macromolecules to achieve a biologically active state in response to chemical cues. Although initially used to describe the impact of small molecules on the conformation and activity of protein enzymes, the definition of this term has been significantly broadened to describe long-range conformational change of macromolecules in response to small or large effectors. Such a broad definition could be applied to RNA molecules, which do not typically serve as protein-free cellular enzymes but fold and form macromolecular assemblies with the help of various ligand molecules, including ions and proteins. Ligand-induced allosteric changes in RNA molecules are often accompanied by cooperative interactions between RNA and its ligand, thus streamlining the folding and assembly pathways. This chapter provides an overview of the interplay between cooperativity and allostery in RNA systems and outlines methods to study these two biological principles.

Structural Insights into RNA Dimerization: Motifs, Interfaces and Functions

In comparison with the pervasive use of protein dimers and multimers in all domains of life, functional RNA oligomers have so far rarely been observed in nature. Their diminished occurrence contrasts starkly with the robust intrinsic potential of RNA to multimerize through long-range base-pairing ("kissing") interactions, self-annealing of palindromic or complementary sequences, and stable tertiary contact motifs, such as the GNRA tetraloop-receptors. To explore the general mechanics of RNA dimerization, we performed a meta-analysis of a collection of exemplary RNA homodimer structures consisting of viral genomic elements, ribozymes, riboswitches, etc., encompassing both functional and fortuitous dimers. Globally, we found that domain-swapped dimers and antiparallel, head-to-tail arrangements are predominant architectural themes. Locally, we observed that the same structural motifs, interfaces and forces that enable tertiary RNA folding also drive their higher-order assemblies. These feature prominently long-range kissing loops, pseudoknots, reciprocal base intercalations and A-minor interactions. We postulate that the scarcity of functional RNA multimers and limited diversity in multimerization motifs may reflect evolutionary constraints imposed by host antiviral immune surveillance and stress sensing. A deepening mechanistic understanding of RNA multimerization is expected to facilitate investigations into RNA and RNP assemblies, condensates, and granules and enable their potential therapeutical targeting.

Former orphan riboswitches reveal unexplored areas of bacterial metabolism, signaling, and gene control processes

Comparative sequence analyses have been used to discover numerous classes of structured noncoding RNAs, some of which are riboswitches that specifically recognize small-molecule or elemental ion ligands and influence expression of adjacent downstream genes. Determining the correct identity of the ligand for a riboswitch candidate typically is aided by an understanding of the genes under its regulatory control. Riboswitches whose ligands were straightforward to identify have largely been associated with well-characterized metabolic pathways, such as coenzyme or amino acid biosynthesis. Riboswitch candidates whose ligands resist identification, collectively known as orphan riboswitches, are often associated with genes coding for proteins of unknown function, or genes for various proteins with no established link to one another. The cognate ligands for 16 former orphan riboswitch motifs have been identified to date. The successful pursuit of the ligands for these classes has provided insight into areas of biology that are not yet fully explored, such as ion homeostasis, signaling networks, and other previously underappreciated biochemical or physiological processes. Herein we discuss the strategies and methods used to match ligands with orphan riboswitch classes, and overview the lessons learned to inform and motivate ongoing efforts to identify ligands for the many remaining candidates.

The asymmetry and cooperativity of tandem glycine riboswitch aptamers

Glycine riboswitches utilize both single- and tandem-aptamer architectures. In the tandem system, the relative contribution of each aptamer toward gene regulation is not well understood. To dissect these contributions, the effects of 684 single mutants of a tandem ON switch from Bacillus subtilis were characterized for the wild-type construct and binding site mutations that selectively restrict ligand binding to either the first or second aptamer. Despite the structural symmetry of tandem aptamers, the response to these mutations was frequently asymmetrical. Mutations in the first aptamer often significantly weakened the K_1/2, while several mutations in the second aptamer improved the amplitude. These results demonstrate that this ON switch favors ligand binding to the first aptamer. This is in contrast to the tandem OFF switch variant from Vibrio cholerae, which was previously shown to have preferential binding to its second aptamer. A bioinformatic analysis of tandem glycine riboswitches revealed that the two binding pockets are differentially conserved between ON and OFF switches. Altogether, this indicates that tandem ON switch variants preferentially utilize binding to the first aptamer to promote helical switching, while OFF switch variants favor binding to the second aptamer. The data set also revealed a cooperative glycine response when both binding pockets were maximally stabilized with three GC base pairs. This indicates a cooperative response may sometimes be obfuscated by a difference in the affinities of the two aptamers. This conditional cooperativity provides an additional layer of tunability to tandem glycine riboswitches that adds to their versatility as genetic switches.

Classification of Riboswitch Families Using Block Location-Based Feature Extraction (BLBFE) Method

Purpose: Riboswitches are special non-coding sequences usually located in mRNAs' un-translated regions and regulate gene expression and consequently cellular function. Furthermore, their interaction with antibiotics has been recently implicated. This raises more interest in development of bioinformatics tools for riboswitch studies. Herein, we describe the development and employment of novel block location-based feature extraction (BLBFE) method for classification of riboswitches. Methods: We have already developed and reported a sequential block finding (SBF) algorithm which, without operating alignment methods, identifies family specific sequential blocks for riboswitch families. Herein, we employed this algorithm for 7 riboswitch families including lysine, cobalamin, glycine, SAM-alpha, SAM-IV, cyclic-di-GMP-I and SAH. Then the study was extended toward implementation of BLBFE method for feature extraction. The outcome features were applied in various classifiers including linear discriminant analysis (LDA), probabilistic neural network (PNN), decision tree and k-nearest neighbors (KNN) classifiers for classification of the riboswitch families. The performance of the classifiers was investigated according to performance measures such as correct classification rate (CCR), accuracy, sensitivity, specificity and f-score. Results: As a result, average CCR for classification of riboswitches was 87.87%. Furthermore, application of BLBFE method in 4 classifiers displayed average accuracies of 93.98% to 96.1%, average sensitivities of 76.76% to 83.61%, average specificities of 96.53% to 97.69% and average f-scores of 74.9% to 81.91%. Conclusion: Our results approved that the proposed method of feature extraction; i.e. BLBFE method; can be successfully used for classification and discrimination of the riboswitch families with high CCR, accuracy, sensitivity, specificity and f-score values.

Regulatory context drives conservation of glycine riboswitch aptamers

In comparison to protein coding sequences, the impact of mutation and natural selection on the sequence and function of non-coding (ncRNA) genes is not well understood. Many ncRNA genes are narrowly distributed to only a few organisms, and appear to be rapidly evolving. Compared to protein coding sequences, there are many challenges associated with assessment of ncRNAs that are not well addressed by conventional phylogenetic approaches, including: short sequence length, lack of primary sequence conservation, and the importance of secondary structure for biological function. Riboswitches are structured ncRNAs that directly interact with small molecules to regulate gene expression in bacteria. They typically consist of a ligand-binding domain (aptamer) whose folding changes drive changes in gene expression. The glycine riboswitch is among the most well-studied due to the widespread occurrence of a tandem aptamer arrangement (tandem), wherein two homologous aptamers interact with glycine and each other to regulate gene expression. However, a significant proportion of glycine riboswitches are comprised of single aptamers (singleton). Here we use graph clustering to circumvent the limitations of traditional phylogenetic analysis when studying the relationship between the tandem and singleton glycine aptamers. Graph clustering enables a broader range of pairwise comparison measures to be used to assess aptamer similarity. Using this approach, we show that one aptamer of the tandem glycine riboswitch pair is typically much more highly conserved, and that which aptamer is conserved depends on the regulated gene. Furthermore, our analysis also reveals that singleton aptamers are more similar to either the first or second tandem aptamer, again based on the regulated gene. Taken together, our findings suggest that tandem glycine riboswitches degrade into functional singletons, with the regulated gene(s) dictating which glycine-binding aptamer is conserved.

Characterization and Engineering of a Clostridium Glycine Riboswitch and Its Use To Control a Novel Metabolic Pathway for 5-Aminolevulinic Acid Production in Escherichia coli

A riboswitch, a regulatory RNA that controls gene expression by specifically binding a ligand, is an attractive genetic element for the control of conditional gene expression and metabolic pathways. In this study, we identified a glycine riboswitch located in the 5'-untranslated regions of a glycine:proton symporter gene in Clostridium pasteurianum. The glycine riboswitch is shown to contain two tandem aptamers and to function as an activator of expression of genes fused to its expression platform. Results of singlet aptamer experiments indicated that aptamer-2 has a much higher impact on regulating gene expression than aptamer-1. Further, we successfully obtained synthetic glycine-OFF riboswitches using a dual selection approach, and one of them repressed gene expression up to 10.2-fold with an improved dynamic range. The specific glycine-OFF riboswitch can function as an independent repressor in the presence of glycine, and its repression mechanism is inferred from predicted secondary structure. The selected glycine-OFF riboswitch was used to dynamically control the biosynthesis of 5-aminolevulinic acid (5-ALA) in Escherichia coli with an unnatural 5-ALA synthetic pathway, in which glycine plays a key role. It is demonstrated that the use of a synthetic Clostridium glycine-OFF riboswitch can lead to a significant increase (11%) of 5-ALA in E. coli harboring an unnatural biosynthetic pathway.

Ion-dependent mobility effects of the Fusobacterium nucleatum glycine riboswitch aptamer II via site-directed spin-labeling (SDSL) electron paramagnetic resonance (EPR)

Site-directed spin-labeling (SDSL) with continuous wave electron paramagnetic resonance (cw-EPR) spectroscopy was utilized to probe site-specific changes in backbone dynamics that accompany folding of the isolated 84 nucleotide aptamer II domain of the Fusobacterium nucleatum (FN) glycine riboswitch. Spin-labels were incorporated using splinted ligation strategies. Results show differential dynamics for spin-labels incorporated into the backbone at a base-paired and loop region. Additionally, the addition of a biologically relevant concentration of 5 mM  Mg²⁺, to an RNA solution with 100 mM K⁺, folds and compacts the structure, inferred by a reduction in spin-label mobility. Furthermore, when controlling for ionic strength, Mg²⁺ added to the RNA induces more folding/less flexibility at the two sites than RNA with K⁺ alone. Addition of glycine does not alter the dynamics of this singlet aptamer II, indicating that the full length riboswitch construct may be needed for glycine binding and induced conformational changes. This work adds to our growing understanding of how splinted-ligation SDSL can be utilized to interrogate differential dynamics in large dynamic RNAs, providing insights into how RNA folding and structure is differentially stabilized by monovalent versus divalent cations.

Gene regulation by a glycine riboswitch singlet uses a finely tuned energetic landscape for helical switching

Riboswitches contain structured aptamer domains that, upon ligand binding, facilitate helical switching in their downstream expression platforms to alter gene expression. To fully dissect how riboswitches function requires a better understanding of the energetic landscape for helical switching. Here, we report a sequencing-based high-throughput assay for monitoring in vitro transcription termination and use it to simultaneously characterize the functional effects of all 522 single point mutants of a glycine riboswitch type-1 singlet. Mutations throughout the riboswitch cause ligand-dependent defects, but only mutations within the terminator hairpin alter readthrough efficiencies in the absence of ligand. A comprehensive analysis of the expression platform reveals that ligand binding stabilizes the antiterminator by just 2-3 kcal/mol, indicating that the competing expression platform helices must be extremely close in energy to elicit a significant ligand-dependent response. These results demonstrate that gene regulation by this riboswitch is highly constrained by the energetics of ligand binding and conformational switching. These findings exemplify the energetic parameters of RNA conformational rearrangements driven by binding events.

A Glycine Riboswitch in Streptococcus pyogenes Controls Expression of a Sodium:Alanine Symporter Family Protein Gene

Regulatory RNAs play important roles in the control of bacterial gene expression. In this study, we investigated gene expression regulation by a putative glycine riboswitch located in the 5'-untranslated region of a sodium:alanine symporter family (SAF) protein gene in the group A Streptococcus pyogenes serotype M49 strain 591. Glycine-dependent gene expression mediated by riboswitch activity was studied using a luciferase reporter gene system. Maximal reporter gene expression was observed in the absence of glycine and in the presence of low glycine concentrations. Differences in glycine-dependent gene expression were not based on differential promoter activity. Expression of the SAF protein gene and the downstream putative cation efflux protein gene was investigated in wild-type bacteria by RT-qPCR transcript analyses. During growth in the presence of glycine (≥1 mM), expression of the genes were downregulated. Northern blot analyses revealed premature transcription termination in the presence of high glycine concentrations. Growth in the presence of 0.1 mM glycine led to the production of a full-length transcript. Furthermore, stability of the SAF protein gene transcript was drastically reduced in the presence of glycine. We conclude that the putative glycine riboswitch in S. pyogenes serotype M49 strain 591 represses expression of the SAF protein gene and the downstream putative cation efflux protein gene in the presence of high glycine concentrations. Sequence and secondary structure comparisons indicated that the streptococcal riboswitch belongs to the class of tandem aptamer glycine riboswitches.

In Vivo Behavior of the Tandem Glycine Riboswitch in Bacillus subtilis

In many bacterial species, the glycine riboswitch is composed of two homologous ligand-binding domains (aptamers) that each bind glycine and act together to regulate the expression of glycine metabolic and transport genes. While the structure and molecular dynamics of the tandem glycine riboswitch have been the subject of numerous in vitro studies, the in vivo behavior of the riboswitch remains largely uncharacterized. To examine the proposed models of tandem glycine riboswitch function in a biologically relevant context, we characterized the regulatory activity of mutations to the riboswitch structure in Bacillus subtilis using β-galactosidase assays. To assess the impact disruptions to riboswitch function have on cell fitness, we introduced these mutations into the native locus of the tandem glycine riboswitch within the B. subtilis genome. Our results indicate that glycine does not need to bind both aptamers for regulation in vivo and mutations perturbing riboswitch tertiary structure have the most severe effect on riboswitch function and gene expression. We also find that in B. subtilis, the glycine riboswitch-regulated gcvT operon is important for glycine detoxification.IMPORTANCE The glycine riboswitch is a unique cis-acting mRNA element that contains two tandem homologous glycine-binding domains that act on a single expression platform to regulate gene expression in response to glycine. While many in vitro experiments have characterized the tandem architecture of the glycine riboswitch, little work has investigated the behavior of this riboswitch in vivo In this study, we analyzed the proposed models of tandem glycine riboswitch regulation in the context of its native locus within the Bacillus subtilis genome and examined how disruptions to glycine riboswitch function impact organismal fitness. Our work offers new insights into riboswitch function in vivo and reinforces the potential of riboswitches as novel antimicrobial targets.

Bioinformatic analysis of riboswitch structures uncovers variant classes with altered ligand specificity

Riboswitches are RNAs that form complex, folded structures that selectively bind small molecules or ions. As with certain groups of protein enzymes and receptors, some riboswitch classes have evolved to change their ligand specificity. We developed a procedure to systematically analyze known riboswitch classes to find additional variants that have altered their ligand specificity. This approach uses multiple-sequence alignments, atomic-resolution structural information, and riboswitch gene associations. Among the discoveries are unique variants of the guanine riboswitch class that most tightly bind the nucleoside 2'-deoxyguanosine. In addition, we identified variants of the glycine riboswitch class that no longer recognize this amino acid, additional members of a rare flavin mononucleotide (FMN) variant class, and also variants of c-di-GMP-I and -II riboswitches that might recognize different bacterial signaling molecules. These findings further reveal the diverse molecular sensing capabilities of RNA, which highlights the potential for discovering a large number of additional natural riboswitch classes.

Synthetic Riboswitches: From Plug and Pray toward Plug and Play

In synthetic biology, metabolic engineering, and gene therapy, there is a strong demand for orthogonal or externally controlled regulation of gene expression. Here, RNA-based regulatory devices represent a promising emerging alternative to proteins, allowing a fast and direct control of gene expression, as no synthesis of regulatory proteins is required. Besides programmable ribozyme elements controlling mRNA stability, regulatory RNA structures in untranslated regions are highly interesting for engineering approaches. Riboswitches are especially well suited, as they show a modular composition of sensor and response elements, allowing a free combination of different modules in a plug-and-play-like mode. The sensor or aptamer domain specifically interacts with a trigger molecule as a ligand, modulating the activity of the adjacent response domain that controls the expression of the genes located downstream, in most cases at the level of transcription or translation. In this review, we discuss the recent advances and strategies for designing such synthetic riboswitches based on natural or artificial components and readout systems, from trial-and-error approaches to rational design strategies. As the past several years have shown dramatic development in this fascinating field of research, we can give only a limited overview of the basic riboswitch design principles that is far from complete, and we apologize for not being able to consider every successful and interesting approach described in the literature.

Singlet glycine riboswitches bind ligand as well as tandem riboswitches

The glycine riboswitch often occurs in a tandem architecture, with two ligand-binding domains (aptamers) followed by a single expression platform. Based on previous observations, we hypothesized that "singlet" versions of the glycine riboswitch, which contain only one aptamer domain, are able to bind glycine if appropriate structural contacts are maintained. An initial alignment of 17 putative singlet riboswitches indicated that the single consensus aptamer domain is flanked by a conserved peripheral stem-loop structure. These singlets were sorted into two subtypes based on whether the active aptamer domain precedes or follows the peripheral stem-loop, and an example of each subtype of singlet riboswitch was characterized biochemically. The singlets possess glycine-binding affinities comparable to those of previously published tandem examples, and the conserved peripheral domains form A-minor interactions with the single aptamer domain that are necessary for ligand-binding activity. Analysis of sequenced genomes identified a significant number of singlet glycine riboswitches. Based on these observations, we propose an expanded model for glycine riboswitch gene control that includes singlet and tandem architectures.

Fluorescence-Based Strategies to Investigate the Structure and Dynamics of Aptamer-Ligand Complexes

In addition to the helical nature of double-stranded DNA and RNA, single-stranded oligonucleotides can arrange themselves into tridimensional structures containing loops, bulges, internal hairpins and many other motifs. This ability has been used for more than two decades to generate oligonucleotide sequences, so-called aptamers, that can recognize certain metabolites with high affinity and specificity. More recently, this library of artificially-generated nucleic acid aptamers has been expanded by the discovery that naturally occurring RNA sequences control bacterial gene expression in response to cellular concentration of a given metabolite. The application of fluorescence methods has been pivotal to characterize in detail the structure and dynamics of these aptamer-ligand complexes in solution. This is mostly due to the intrinsic high sensitivity of fluorescence methods and also to significant improvements in solid-phase synthesis, post-synthetic labeling strategies and optical instrumentation that took place during the last decade. In this work, we provide an overview of the most widely employed fluorescence methods to investigate aptamer structure and function by describing the use of aptamers labeled with a single dye in fluorescence quenching and anisotropy assays. The use of 2-aminopurine as a fluorescent analog of adenine to monitor local changes in structure and fluorescence resonance energy transfer (FRET) to follow long-range conformational changes is also covered in detail. The last part of the review is dedicated to the application of fluorescence techniques based on single-molecule microscopy, a technique that has revolutionized our understanding of nucleic acid structure and dynamics. We finally describe the advantages of monitoring ligand-binding and conformational changes, one molecule at a time, to decipher the complexity of regulatory aptamers and summarize the emerging folding and ligand-binding models arising from the application of these single-molecule FRET microscopy techniques.

Engineering and characterization of fluorogenic glycine riboswitches

A set of 12 fluorogenic glycine riboswitches with different thermodynamic and kinetic response properties was engineered. For the design of functional riboswitches, a three-part RNA approach was applied based on the idea of linking a RNA sensor, transmitter and actuator part together. For the RNA sensor and actuator part, we used the tandem glycine aptamer structure from Bacillus subtillis, and fluorogenic aptamer Spinach, respectively. To achieve optimal signal transduction from the sensor to the actuator, a riboswitch library with variable transmitter was screened with a microfluidic large-scale integration chip. This allowed us to establish the complete thermodynamic binding profiles of the riboswitch library. Glycine dissociation constants of the 12 strong fluorescence response riboswitches varied between 99.7 and 570 μM. Furthermore, the kinetic glycine binding (k^on), and dissociation (k^off) rates, and corresponding energy barriers of the 10 strongest fluorescence response riboswitches were determined with the same chip platform. k^on and k^off were in the order of 10⁻³s⁻¹ and 10⁻²s⁻¹, respectively. Conclusively, we demonstrate that systematic screening of synthetic and natural linked RNA parts with microfluidic chip technology is an effective approach to rapidly generate fluorogenic metabolite riboswitches with a broad range of biophysical response properties.

Toward increased concentration sensitivity for continuous wave EPR investigations of spin-labeled biological macromolecules at high fields

High-field, high-frequency electron paramagnetic resonance (EPR) spectroscopy at W-(∼94 GHz) and D-band (∼140 GHz) is important for investigating the conformational dynamics of flexible biological macromolecules because this frequency range has increased spectral sensitivity to nitroxide motion over the 100 ps to 2 ns regime. However, low concentration sensitivity remains a roadblock for studying aqueous samples at high magnetic fields. Here, we examine the sensitivity of a non-resonant thin-layer cylindrical sample holder, coupled to a quasi-optical induction-mode W-band EPR spectrometer (HiPER), for continuous wave (CW) EPR analyses of: (i) the aqueous nitroxide standard, TEMPO; (ii) the unstructured to α-helical transition of a model IDP protein; and (iii) the base-stacking transition in a kink-turn motif of a large 232 nt RNA. For sample volumes of ∼50 μL, concentration sensitivities of 2-20 μM were achieved, representing a ∼10-fold enhancement compared to a cylindrical TE011 resonator on a commercial Bruker W-band spectrometer. These results therefore highlight the sensitivity of the thin-layer sample holders employed in HiPER for spin-labeling studies of biological macromolecules at high fields, where applications can extend to other systems that are facilitated by the modest sample volumes and ease of sample loading and geometry.

Fab Chaperone-Assisted RNA Crystallography (Fab CARC)

Recent discovery of structured RNAs such as ribozymes and riboswitches shows that there is still much to learn about the structure and function of RNAs. Knowledge learned can be employed in both biochemical research and clinical applications. X-ray crystallography gives unparalleled atomic-level structural detail from which functional inferences can be deduced. However, the difficulty in obtaining high-quality crystals and their phasing information make it a very challenging task. RNA crystallography is particularly arduous due to several factors such as RNA's paucity of surface chemical diversity, lability, repetitive anionic backbone, and flexibility, all of which are counterproductive to crystal packing. Here we describe Fab chaperone assisted RNA crystallography (CARC), a systematic technique to increase RNA crystallography success by facilitating crystal packing as well as expediting phase determination through molecular replacement of conserved Fab domains. Major steps described in this chapter include selection of a synthetic Fab library displayed on M13 phage against a structured RNA crystallization target, ELISA for initial choice of binding Fabs, Fab expression followed by protein A affinity then cation exchange chromatography purification, final choice of Fab by binding specificity and affinity as determined by a dot blot assay, and lastly gel filtration purification of a large quantity of chosen Fabs for crystallization.

Biophysical Approaches to Bacterial Gene Regulation by Riboswitches

The last decade has witnessed the discovery of a variety of non-coding RNA sequences that perform a broad range of crucial biological functions. Among these, the ability of certain RNA sequences, so-called riboswitches, has attracted considerable interest. Riboswitches control gene expression in response to the concentration of particular metabolites to which they bind without the need for any protein. These RNA switches not only need to adopt a very specific tridimensional structure to perform their function, but also their sequence has been evolutionary optimized to recognize a particular metabolite with high affinity and selectivity. Thus, riboswitches offer a unique opportunity to get fundamental insights into RNA plasticity and how folding dynamics and ligand recognition mechanisms have been efficiently merged to control gene regulation. Because riboswitch sequences have been mostly found in bacterial organisms controlling the expression of genes associated to the synthesis, degradation or transport of crucial metabolites for bacterial survival, they offer exciting new routes for antibiotic development in an era where bacterial resistance is more than ever challenging conventional drug discovery strategies. Here, we give an overview of the architecture, diversity and regulatory mechanisms employed by riboswitches with particular emphasis on the biophysical methods currently available to characterise their structure and functional dynamics.

Optimizing RNA structures by sequence extensions using RNAcop

A key aspect of RNA secondary structure prediction is the identification of novel functional elements. This is a challenging task because these elements typically are embedded in longer transcripts where the borders between the element and flanking regions have to be defined. The flanking sequences impact the folding of the functional elements both at the level of computational analyses and when the element is extracted as a transcript for experimental analysis. Here, we analyze how different flanking region lengths impact folding into a constrained structure by computing probabilities of folding for different sizes of flanking regions. Our method, RNAcop (RNA context optimization by probability), is tested on known and de novo predicted structures. In vitro experiments support the computational analysis and suggest that for a number of structures, choosing proper lengths of flanking regions is critical. RNAcop is available as web server and stand-alone software via http://rth.dk/resources/rnacop.

Cooperativity, allostery and synergism in ligand binding to riboswitches

Recent progress in identification and characterization of novel types of non-coding RNAs has proven that RNAs carry out a variety of cellular functions ranging from scaffolding to gene expression control. In both prokaryotic and eukaryotic cells, several classes of non-coding RNAs control expression of dozens of genes in response to specific cues. One of the most interesting and outstanding questions in the RNA field is whether regulatory RNAs are capable of employing basic biological concepts, such as allostery and cooperativity, previously attributed to the function of proteins. Aside from regulatory RNAs that form complementary base pairing with their nucleic acid targets, several RNA classes modulate gene expression via molecular mechanisms which can be paralleled to protein-mediated regulation. Among these RNAs are riboswitches, metabolite-sensing non-coding regulatory elements that adopt intrinsic three-dimensional structures and specifically bind various small molecule ligands. These characteristics of riboswitches make them well-suited for complex regulatory responses observed in allosteric and cooperative protein systems. Here we present an overview of the biochemical, genetic, and structural studies of riboswitches with a major focus on complex regulatory mechanisms and biological principles utilized by riboswitches for such genetic modulation.

RNA quaternary structure and global symmetry

Many proteins associate into symmetric multisubunit complexes. Structural analyses suggested that, by contrast, virtually all RNAs with complex 3D structures function as asymmetric monomers. Recent crystal structures revealed that several biological RNAs exhibit global symmetry at the level of their tertiary and quaternary structures. Here we survey known examples of global RNA symmetry, including the true quaternary symmetry of the bacteriophage ϕ29 prohead RNA (pRNA) and the internal pseudosymmetry of the single-chain flavin mononucleotide (FMN), glycine, and cyclic di-AMP (c-di-AMP) riboswitches. For these RNAs, global symmetry stabilizes the RNA fold, coordinates ligand-RNA interactions, and facilitates association with symmetric binding partners.

Ligand binding by the tandem glycine riboswitch depends on aptamer dimerization but not double ligand occupancy

The glycine riboswitch predominantly exists as a tandem structure, with two adjacent, homologous ligand-binding domains (aptamers), followed by a single expression platform. The recent identification of a leader helix, the inclusion of which eliminates cooperativity between the aptamers, has reopened the debate over the purpose of the tandem structure of the glycine riboswitch. An equilibrium dialysis-based assay was combined with binding-site mutations to monitor glycine binding in each ligand-binding site independently to understand the role of each aptamer in glycine binding and riboswitch tertiary interactions. A series of mutations disrupting the dimer interface was used to probe how dimerization impacts ligand binding by the tandem glycine riboswitch. While the wild-type tandem riboswitch binds two glycine equivalents, one for each aptamer, both individual aptamers are capable of binding glycine when the other aptamer is unoccupied. Intriguingly, glycine binding by aptamer-1 is more sensitive to dimerization than glycine binding by aptamer-2 in the context of the tandem riboswitch. However, monomeric aptamer-2 shows dramatically weakened glycine-binding affinity. In addition, dimerization of the two aptamers in trans is dependent on glycine binding in at least one aptamer. We propose a revised model for tandem riboswitch function that is consistent with these results, wherein ligand binding in aptamer-1 is linked to aptamer dimerization and stabilizes the P1 stem of aptamer-2, which controls the expression platform.

The purine riboswitch as a model system for exploring RNA biology and chemistry

Over the past decade the purine riboswitch, and in particular its nucleobase-binding aptamer domain, has emerged as an important model system for exploring various aspects of RNA structure and function. Its relatively small size, structural simplicity and readily observable activity enable application of a wide variety of experimental approaches towards the study of this RNA. These analyses have yielded important insights into small molecule recognition, co-transcriptional folding and secondary structural switching, and conformational dynamics that serve as a paradigm for other RNAs. In this article, the current state of understanding of the purine riboswitch family and how this growing knowledge base is starting to be exploited in the creation of novel RNA devices are examined. This article is part of a Special Issue entitled: Riboswitches.

DNA-rescuable allosteric inhibition of aptamer II ligand affinity by aptamer I element in the shortened Vibrio cholerae glycine riboswitch

Glycine riboswitches contain two aptamers and turn on the expression of downstream genes in bacteria. Although full-length glycine riboswitches were shown to exhibit no glycine-binding cooperativity, the truncated glycine riboswitches were confirmed to bind two glycine molecules cooperatively. Thorough understanding of the ligand-binding cooperativity may shed light on the molecular basis of the cooperativity and help design novel intricate biosensing genetic circuits for application in synthetic biology. A previously proposed sequential model does not readily provide explanation for published data showing a deleterious mutation in the first aptamer inhibiting the glycine binding of the second one. Using the glycine riboswitch from Vibrio cholerae as a model system, we have identified a region in the first aptamer that modulates the second aptamer function especially in the shortened glycine riboswitch. Importantly, this modulation can be rescued by the addition of a complementary oligodeoxynucleotide, demonstrating the feasibility of developing this system into novel genetic circuits that sense both glycine and a DNA signal.

Characterizing the dynamics of the leader-linker interaction in the glycine riboswitch with site-directed spin labeling

Site-directed spin labeling with continuous wave electron paramagnetic resonance (EPR) spectroscopy was utilized to characterize dynamic features of the kink–turn motif formed through a leader–linker interaction in the Vibrio cholerae glycine riboswitch. Efficient incorporation of spin-labels into select sites within the phosphate backbone of the leader–linker region proceeded via splinted ligation of chemically synthesized spin-labeled oligonucleotides to in vitro transcribed larger RNA fragments. The resultant nitroxide EPR line shapes have spectral characteristics consistent with a kink–turn motif and reveal differential backbone dynamics that are modulated by the presence of magnesium, potassium, and glycine.

Specific RNA-binding antibodies with a four-amino-acid code

Numerous large non-coding RNAs are rapidly being discovered, and many of them have been shown to play vital roles in gene expression, gene regulation, and human diseases. Given their often structured nature, specific recognition with an antibody fragment becomes feasible and may help define the structure and function of these non-coding RNAs. As demonstrated for protein antigens, specific antibodies may aid in RNA crystal structure elucidation or the development of diagnostic tools and therapeutic drugs targeting disease-causing RNAs. Recent success and limitation of RNA antibody development has made it imperative to generate an effective antibody library specifically targeting RNA molecules. Adopting the reduced chemical diversity design and further restricting the interface diversity to tyrosines, serines, glycines, and arginines only, we have constructed a RNA-targeting Fab library. Phage display selection and downstream characterization showed that this library yielded high-affinity Fabs for all three RNA targets tested. Using a quantitative specificity assay, we found that these Fabs are highly specific, possibly due to the alternate codon design we used to avoid consecutive arginines in the Fab interface. In addition, the effectiveness of the minimal Fab library may challenge our view of the protein-RNA binding interface and provide a unique solution for future design of RNA-binding proteins.

Site-directed spin-labeling strategies and electron paramagnetic resonance spectroscopy for large riboswitches

Genetic regulation effected by RNA riboswitches is governed by ligand-induced structural reorganization with modulation of RNA conformation and dynamics. Characterization of the conformational states of riboswitches in the presence or absence of salts and ligands is important for understanding how interconversion of riboswitch RNA folding states influences function. The methodology of site-directed spin labeling (SDSL) coupled with electron paramagnetic resonance (EPR) spectroscopy is suitable for such studies, wherein site-specific incorporation of a nitroxide radical spin probe allows for local dynamics and conformational changes to be investigated. This chapter reviews a strategy for SDSL-EPR studies of large riboswitches and uses the full length 232 nucleotide (nt) kink-turn motif-containing Vibrio cholerae (VC) glycine riboswitch as an example. Spin-labeling strategies and the challenges of incorporating spin labels into large riboswitches are reviewed and the approach to overcome these challenges is described. Results are subsequently presented illustrating changes in dynamics within the labeled region of the VC glycine riboswitch as observed using SDSL-EPR.

Glycols modulate terminator stem stability and ligand-dependency of a glycine riboswitch

The Bacillus subtilis glycine riboswitch comprises tandem glycine-binding aptamers and a putative terminator stem followed by the gcvT operon. Gene expression is regulated via the sensing of glycine. However, we found that the riboswitch behaves in a "glycine-independent" manner in the presence of polyethylene glycol (PEG) and ethylene glycol. The effect is related to the formation of a terminator stem within the expression platform under such conditions. The results revealed that increasing PEG stabilized the structure of the terminator stem. By contrast, the addition of ethylene glycol destabilized the terminator stem. PEG and ethylene glycol have opposite effects on transcription as well as on stable terminator stem formation. The glycine-independency of the riboswitch and the effects of such glycols might shed light on the evolution of riboswitches.

Refinement of the Listeria monocytogenes σB regulon through quantitative proteomic analysis

σ(B) is an alternative σ factor that regulates stress response and virulence genes in the foodborne pathogen Listeria monocytogenes. To gain further insight into σ(B)-dependent regulatory mechanisms in L. monocytogenes, we (i) performed quantitative proteomic comparisons between the L. monocytogenes parent strain 10403S and an isogenic ΔsigB mutant and (ii) conducted a meta-analysis of published microarray studies on the 10403S σ(B) regulon. A total of 134 genes were found to be significantly positively regulated by σ(B) at the transcriptomic level with >75 % of these genes preceded by putative σ(B)-dependent promoters; 21 of these 134 genes were also found to be positively regulated by σ(B) through proteomics. In addition, 15 proteins were only found to be positively regulated by σ(B) through proteomics analyses, including Lmo1349, a putative glycine cleavage system protein. The lmo1349 gene is preceded by a 5' UTR that functions as a glycine riboswitch, which suggests regulation of glycine metabolism by σ(B) in L. monocytogenes. Herein, we propose a model where σ(B) upregulates pathways that facilitate biosynthesis and uptake of glycine, which may then activate this riboswitch. Our data also (i) identified a number of σ(B)-dependent proteins that appear to be encoded by genes that are co-regulated by multiple transcriptional regulators, in particular PrfA, and (ii) found σ(B)-dependent genes and proteins to be overrepresented in the 'energy metabolism' role category, highlighting contributions of the σ(B) regulon to L. monocytogenes energy metabolism as well as a role of PrfA and σ(B) interaction in regulating aspects of energy metabolism in L. monocytogenes.

A decade of riboswitches

Riboswitches were discovered in 2002 in bacteria as RNA-based intracellular sensors of vitamin derivatives. During the last decade, naturally occurring RNA sensor elements have been found to bind a range of small metabolites and ions and to exert regulatory control of transcription, translation, splicing, and RNA stability. Extensive biochemical, structural, and genetic studies have established the basic principles underpinning riboswitch function in all three kingdoms of life with implications for developing antibiotics, designing new molecular sensors, and integrating riboswitches into synthetic circuits.

Modulation of quaternary structure and enhancement of ligand binding by the K-turn of tandem glycine riboswitches

Most known glycine riboswitches have two homologous aptamer domains arranged in tandem and separated by a short linker. The two aptamers associate through reciprocal "quaternary" interactions that have been proposed to result in cooperative glycine binding. Recently, the interaptamer linker was found to form helix P0 with a previously unrecognized segment 5' to the first aptamer domain. P0 was shown to increase glycine affinity, abolish cooperativity, and conform to the K-turn motif consensus. We examine the global thermodynamic and structural role of P0 using isothermal titration calorimetry (ITC) and small-angle X-ray scattering (SAXS), respectively. To evaluate the generality of P0 function, we prepared glycine riboswitch constructs lacking and including P0 from Bacillus subtilis, Fusobacterium nucleatum, and Vibrio cholerae. We find that P0 indeed folds into a K-turn, supports partial pre-folding of all three glycine-free RNAs, and is required for ITC observation of glycine binding under physiologic Mg(2+) concentrations. Except for the unusually small riboswitch from F. nucleatum, the K-turn is needed for maximally compacting the glycine-bound states of the RNAs. Formation of a ribonucleoprotein complex between the B. subtilis or the F. nucleatum RNA constructs and the bacterial K-turn binding protein YbxF promotes additional folding of the free riboswitch, and enhances glycine binding. Consistent with the previously reported loss of cooperativity, P0-containing B. subtilis and V. cholerae tandem aptamers bound no more than one glycine molecule per riboswitch. Our results indicate that the P0 K-turn helps organize the quaternary structure of tandem glycine riboswitches, thereby facilitating ligand binding under physiologic conditions.