Identification of a tertiary interaction important for cooperative ligand binding by the glycine riboswitch

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.

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.

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.

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

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.

Magnesium controls aptamer-expression platform switching in the SAM-I riboswitch

Investigations of most riboswitches remain confined to the ligand-binding aptamer domain. However, during the riboswitch mediated transcription regulation process, the aptamer domain and the expression platform compete for a shared strand. If the expression platform dominates, an anti-terminator helix is formed, and the transcription process is active (ON state). When the aptamer dominates, transcription is terminated (OFF state). Here, we use an expression platform switching experimental assay and structure-based electrostatic simulations to investigate this ON-OFF transition of the full length SAM-I riboswitch and its magnesium concentration dependence. Interestingly, we find the ratio of the OFF population to the ON population to vary non-monotonically as magnesium concentration increases. Upon addition of magnesium, the aptamer domain pre-organizes, populating the OFF state, but only up to an intermediate magnesium concentration level. Higher magnesium concentration preferentially stabilizes the anti-terminator helix, populating the ON state, relatively destabilizing the OFF state. Magnesium mediated aptamer-expression platform domain closure explains this relative destabilization of the OFF state at higher magnesium concentration. Our study reveals the functional potential of magnesium in controlling transcription of its downstream genes and underscores the importance of a narrow concentration regime near the physiological magnesium concentration ranges, striking a balance between the OFF and ON states in bacterial gene regulation.

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.

Understanding in-line probing experiments by modeling cleavage of nonreactive RNA nucleotides

Ribonucleic acid (RNA) is involved in many regulatory and catalytic processes in the cell. The function of any RNA molecule is intimately related with its structure. In-line probing experiments provide valuable structural data sets for a variety of RNAs and are used to characterize conformational changes in riboswitches. However, the structural determinants that lead to differential reactivities in unpaired nucleotides have not been investigated yet. In this work, we used a combination of theoretical approaches, i.e., classical molecular dynamics simulations, multiscale quantum mechanical/molecular mechanical calculations, and enhanced sampling techniques in order to compute and interpret the differential reactivity of individual residues in several RNA motifs, including members of the most important GNRA and UNCG tetraloop families. Simulations on the multinanosecond timescale are required to converge the related free-energy landscapes. The results for uGAAAg and cUUCGg tetraloops and double helices are compared with available data from in-line probing experiments and show that the introduced technique is able to distinguish between nucleotides of the uGAAAg tetraloop based on their structural predispositions toward phosphodiester backbone cleavage. For the cUUCGg tetraloop, more advanced ab initio calculations would be required. This study is the first attempt to computationally classify chemical probing experiments and paves the way for an identification of tertiary structures based on the measured reactivity of nonreactive nucleotides.

Discovery of a Structurally Unique Small Molecule that Inhibits Protein Synthesis

Identifying and characterizing natural products and synthetic small molecules that inhibit biochemical processes such as ribosomal translation can lead to novel sources of molecular probes and therapeutics. The search for new antibiotics has been invigorated by the increasing burden of drug-resistant bacteria and has identified many clinically essential prokaryote-specific ribosome inhibitors. However, the current cohort of antibiotics is limited with regards to bacterial resistance mechanisms because of structural similarity within classes. From a high-throughput screen for translation inhibitors, we discovered a new compound, T6102, which inhibits bacterial protein synthesis in vitro, inhibits bacterial growth of Bacillus subtilis in vivo, and has a chemical structure that appears to be unique among known classes of translation-inhibiting antibiotics. T6102's unique structure compared to current clinically-utilized antibiotics makes it an exciting new candidate for the development of next-generation antibiotics.

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.

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.

Cooperative binding: a multiple personality

Cooperative binding has been described in many publications and has been related to or defined by several different properties of the binding behavior of the ligand to the target molecule. In addition to the commonly used Hill coefficient, other characteristics such as a sigmoidal shape of the overall titration curve in a linear plot, a change of ligand affinity of the other binding sites when a site of the target molecule becomes occupied, or complex roots of the binding polynomial have been used to define or to quantify cooperative binding. In this work, we analyze how the different properties are related in the most general model for binding curves based on the grand canonical partition function and present several examples which highlight differences between the cooperativity characterizing properties which are discussed. Our results mainly show that among the presented definitions there are not two which fully coincide. Moreover, this work poses the question whether it can make sense to distinguish between positive and negative cooperativity based on the macroscopic binding isotherm only. This article shall emphasize that scientists who investigate cooperative effects in biological systems could help avoiding misunderstandings by stating clearly which kind of cooperativity they discuss.

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.

Consistent global structures of complex RNA states through multidimensional chemical mapping

Accelerating discoveries of non-coding RNA (ncRNA) in myriad biological processes pose major challenges to structural and functional analysis. Despite progress in secondary structure modeling, high-throughput methods have generally failed to determine ncRNA tertiary structures, even at the 1-nm resolution that enables visualization of how helices and functional motifs are positioned in three dimensions. We report that integrating a new method called MOHCA-seq (Multiplexed •OHCleavage Analysis with paired-end sequencing) with mutate-and-map secondary structure inference guides Rosetta 3D modeling to consistent 1-nm accuracy for intricately folded ncRNAs with lengths up to 188 nucleotides, including a blind RNA-puzzle challenge, the lariat-capping ribozyme. This multidimensional chemical mapping (MCM) pipeline resolves unexpected tertiary proximities for cyclic-di-GMP, glycine, and adenosylcobalamin riboswitch aptamers without their ligands and a loose structure for the recently discovered human HoxA9D internal ribosome entry site regulon. MCM offers a sequencing-based route to uncovering ncRNA 3D structure, applicable to functionally important but potentially heterogeneous states.

DOI:http://dx.doi.org/10.7554/eLife.07600.001

Our genetic material, in the form of molecules of DNA, provides instructions for many different processes in our cells. To issue these instructions, particular sections of DNA are copied to make a type of molecule called ribonucleic acid (RNA). Some of these RNA molecules contain instructions to make proteins, but others—known as non-coding RNAs—regulate the activity of genes in cells.

The genetic information within RNA is encoded by the sequence of four different chemical parts called ‘nucleotides’. RNA can exist as a single strand of nucleotides, but the nucleotides can also pair up in specific combinations to form sections of double-stranded RNA. Therefore, a single strand of non-coding RNA can fold into a complex three-dimensional shape that contains loops, twists, and bulges.

The three-dimensional structures of non-coding RNAs are crucial for their roles in cells, but the variety and complexity of shapes that they can form makes it technically difficult to study them. In 2008, researchers developed a new method called MOHCA that can map the positions of nucleotides that are close together in the three-dimensional structure. Highly reactive chemicals are attached to the nucleotides and these can react with, and damage, other nearby nucleotides. By detecting which nucleotides have been damaged, it is possible to map the positions of these nucleotides and decipher the structure of the RNA molecule using computer algorithms.

MOHCA is a promising approach, but the initial methods to find the damaged nucleotides were tedious and required specialized equipment. Now, Cheng, Das et al.—including some of the researchers involved in the 2008 work—have developed an improved version of MOHCA that uses readily available RNA sequencing techniques to find the damaged nucleotides. The RNA sequencing data are then analyzed by a new algorithm in the Rosetta computer modeling software.

Cheng, Das et al. used this newly developed ‘MOHCA-seq’ and Rosetta to reveal the structures of a human non-coding RNA and several other non-coding RNA molecules to a much higher level of detail than before. Together, MOHCA-seq and Rosetta provide a rapid method for researchers to decipher the three-dimensional structure of non-coding RNAs. This method is likely to speed up the analysis of the complex structures of non-coding RNAs. It will be useful in future efforts to work out what roles these RNAs play in cells, including their activity in cancer, neurodegeneration, and other diseases.

DOI:http://dx.doi.org/10.7554/eLife.07600.002

Modeling complex RNA tertiary folds with Rosetta

Reliable modeling of RNA tertiary structures is key to both understanding these structures' roles in complex biological machines and to eventually facilitating their design for molecular computing and robotics. In recent years, a concerted effort to improve computational prediction of RNA structure through the RNA-Puzzles blind prediction trials has accelerated advances in the field. Among other approaches, the versatile and expanding Rosetta molecular modeling software now permits modeling of RNAs in the 100-300 nucleotide size range at consistent subhelical (~1 nm) resolution. Our laboratory's current state-of-the-art methods for RNAs in this size range involve Fragment Assembly of RNA with Full-Atom Refinement (FARFAR), which optimizes RNA conformations in the context of a physically realistic energy function, as well as hybrid techniques that leverage experimental data to inform computational modeling. In this chapter, we give a practical guide to our current workflow for modeling RNA three-dimensional structures using FARFAR, including strategies for using data from multidimensional chemical mapping experiments to focus sampling and select accurate conformations.

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.

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.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

Two glycine riboswitches activate the glycine cleavage system essential for glycine detoxification in Streptomyces griseus

The glycine cleavage (GCV) system catalyzes the oxidative cleavage of glycine into CO2, NH4(+), and a methylene group, which is accepted by tetrahydrofolate (THF) to form N(5),N(10)-methylene-THF. Streptomyces griseus contains gcvP and the gcvT-gcvH operon, which encode three intrinsic components of the GCV system. We identified the transcriptional start sites of gcvTH and gcvP and found putative glycine riboswitches in their 5' untranslated regions (5' UTRs). The ratios of the transcripts of the gcvT and gcvP coding sequences (CDSs) to those of the respective 5' UTRs were significantly higher in the presence of glycine in the wild-type strain. However, the levels of gcvT and gcvP CDS transcripts were not increased by glycine in the respective 5' UTR deletion mutants. A reporter gene assay showed that a transcriptional terminator exists in the 5' UTR of gcvTH. Furthermore, by an in-line probing assay, we confirmed that glycine bound directly to the putative riboswitch RNAs. These results indicate that the S. griseus glycine riboswitches enhance transcriptional read-through to the downstream CDSs, like known glycine riboswitches in other bacteria. We examined the growth of three mutants in which either or both of the gcvTH and gcvP 5' UTRs were deleted. Like the wild-type strain, all mutants grew vigorously in a medium containing 0.9% glucose as a carbon source. However, the mutants showed severely restricted growth in a medium containing 0.9% glucose and 1% glycine, while the wild-type strain grew normally. This indicates that glycine has a growth-inhibitory effect and that the GCV system plays a critical role in glycine detoxification in S. griseus.

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.

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.

Molecular recognition and function of riboswitches

Regulatory mRNAs elements termed riboswitches respond to elevated concentrations of cellular metabolites by modulating expression of associated genes. Riboswitches attain their high metabolite selectivity by capitalizing on the intrinsic tertiary structures of their sensor domains. Over the years, riboswitch structure and folding have been amongst the most researched topics in the RNA field. Most recently, novel structures of single-ligand and cooperative double-ligand sensors have broadened our knowledge of architectural and molecular recognition principles exploited by riboswitches. The structural information has been complemented by extensive folding studies, which have provided several important clues on the formation of ligand-competent conformations and mechanisms of ligand discrimination. These studies have greatly improved our understanding of molecular events in riboswitch-mediated gene expression control and provided the molecular basis for intervention into riboswitch-controlled genetic circuits.

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

Comprised of two aptamers connected by a short nucleotide linker, the glycine riboswitch was the first example of naturally occurring RNA elements reported to bind small organic molecules cooperatively. Earlier works have shown binding of glycine to the second aptamer allows tertiary interactions to be made between the two aptamers, which facilitates binding of a separate glycine molecule to the first aptamer, leading to glycine-binding cooperativity. Prompted by a distinctive protection pattern in the linker region of a minimal glycine riboswitch construct, we have identified a highly conserved (>90%) leader-linker duplex involving leader nucleotides upstream of the previously reported consensus glycine riboswitch sequences. In >50% of the glycine riboswitches, the leader-linker interaction forms a kink-turn motif. Characterization of three glycine ribsowitches showed that the leader-linker interaction improved the glycine-binding affinities by 4.5- to 86-fold. In-line probing and native gel assays with two aptamers in trans suggested synergistic action between glycine-binding and interaptamer interaction during global folding of the glycine riboswitch. Mutational analysis showed that there appeared to be no ligand-binding cooperativity in the glycine riboswitch when the leader-linker interaction is present, and the previously measured cooperativity is simply an artifact of a truncated construct missing the leader sequence.

Automated RNA structure prediction uncovers a kink-turn linker in double glycine riboswitches

The tertiary structures of functional RNA molecules remain difficult to decipher. A new generation of automated RNA structure prediction methods may help address these challenges but have not yet been experimentally validated. Here we apply four prediction tools to a class of double glycine riboswitches that can bind two ligands cooperatively. A novel method (BPPalign), RMdetect, JAR3D, and Rosetta 3D modeling give consistent predictions for a new stem P0 and a kink-turn motif. These elements structure the linker between the RNAs' double aptamers. Chemical mapping on the Fusobacterium nucleatum riboswitch with N-methylisatoic anhydride, dimethyl sulfate and 1-cyclohexyl-3-(2-morpholinoethyl)carbodiimide metho-p-toluenesulfonate probing, mutate-and-map studies, and mutation/rescue experiments all provide strong evidence for the structured linker. Under solution conditions that permit rigorous thermodynamic analysis, disrupting this helix-junction-helix structure gives 120- and 6-30-fold poorer dissociation constants for the RNA's two glycine-binding transitions, corresponding to an overall energetic impact of 4.3 ± 0.5 kcal/mol. Prior biochemical and crystallography studies did not include this critical element due to over-truncation of the RNA. We speculate that several further undiscovered elements are likely to exist in the flanking regions of this and other functional RNAs, and automated prediction tools can play a useful role in their detection and dissection.

Structural basis of cooperative ligand binding by the glycine riboswitch

The glycine riboswitch regulates gene expression through the cooperative recognition of its amino acid ligand by a tandem pair of aptamers. A 3.6 Å crystal structure of the tandem riboswitch from the glycine permease operon of Fusobacterium nucleatum reveals the glycine binding sites and an extensive network of interactions, largely mediated by asymmetric A-minor contacts, that serve to communicate ligand binding status between the aptamers. These interactions provide a structural basis for how the glycine riboswitch cooperatively regulates gene expression.