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

Crystal structure of Escherichia coli thiamine pyrophosphate-sensing riboswitch in the apo state

The thiamine pyrophosphate (TPP)-sensing riboswitch is one of the earliest discovered and most widespread riboswitches. Numerous structural studies have been reported for this riboswitch bound with various ligands. However, the ligand-free (apo) structure remains unknown. Here, we report a 3.1 Å resolution crystal structure of Escherichia coli TPP riboswitch in the apo state, which exhibits an extended, Y-shaped conformation further supported by small-angle X-ray scattering data and driven molecular dynamics simulations. The loss of ligand interactions results in helical uncoiling of P5 and disruption of the key tertiary interaction between the sensory domains. Opening of the aptamer propagates to the gene-regulatory P1 helix and generates the key conformational flexibility needed for the switching behavior. Much of the ligand-binding site at the three-way junction is unaltered, thereby maintaining a partially preformed pocket. Together, these results paint a dynamic picture of the ligand-induced conformational changes in TPP riboswitches that confer conditional gene regulation.

Ornate, large, extremophilic (OLE) RNA forms a kink turn necessary for OapC protein recognition and RNA function

Ornate, large, extremophilic (OLE) RNAs represent a class of noncoding RNAs prevalent in Gram-positive, extremophilic/anaerobic bacterial species. OLE RNAs (∼600 nt), whose precise biochemical functions remain mysterious, form an intricate secondary structure interspersed with regions of highly conserved nucleotides. In the alkali-halophilic bacterium Bacillus halodurans, OLE RNA is a component of a ribonucleoprotein (RNP) complex involving at least two proteins named OapA and OapB, but additional components may exist that could point to functional roles for the RNA. Disruption of the genes for either OLE RNA, OapA, or OapB result in the inability of cells to overcome cold, alcohol, or Mg²⁺ stresses. In the current study, we used in vivo crosslinking followed by OLE RNA isolation to identify the protein YbxF as a potential additional partner in the OLE RNP complex. Notably, a mutation in the gene for this same protein was also reported to be present in a strain wherein the complex is nonfunctional. The B. halodurans YbxF (herein renamed OapC) is homologous to a bacterial protein earlier demonstrated to bind kink turn (k-turn) RNA structural motifs. In vitro RNA-protein binding assays reveal that OLE RNA forms a previously unrecognized k-turn that serves as the natural binding site for YbxF/OapC. Moreover, B. halodurans cells carrying OLE RNAs with disruptive mutations in the k-turn exhibit phenotypes identical to cells lacking functional OLE RNP complexes. These findings reveal that the YbxF/OapC protein of B. halodurans is important for the formation of a functional OLE RNP complex.

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.

Structure and ligand binding of the ADP-binding domain of the NAD+ riboswitch

The nadA motif is the first known NAD+-dependent riboswitch, comprising two similar tandem bulged stem-loop structures. We have determined the structure of the 5' domain 1 of the riboswitch. It has three coaxial helical segments, separated by an ACANCCCC bulge and by an internal loop, with a tertiary contact between them that includes two C:G base pairs. We have determined the structure with a number of ligands related to NADH, but in each case only the ADP moiety is observed. The adenosine adopts an anti conformation, forms multiple hydrogen bonds across the width of the sugar edge of the penultimate C:G base pair of the helix preceding the bulge, and the observed contacts have been confirmed by mutagenesis and calorimetry. Two divalent metal ions play a key structural role at the narrow neck of the bulge. One makes direct bonding contacts to the diphosphate moiety, locking it into position. Thus the nucleobase, ribose, and phosphate groups of the ADP moiety are all specifically recognized by the RNA. The NAD+ riboswitch is modular. Domain 1 is an ADP binding domain that may be ancient and could potentially be used in combination with other ligand binding motifs such as CoA.

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.

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.

The kink-turn in the structural biology of RNA

The kink-turn (k-turn) is a widespread structural motif found in functional RNA species. It typically comprises a three-nucleotide bulge followed by tandem trans sugar edge-Hoogsteen G:A base pairs. It introduces a sharp kink into the axis of duplex RNA, juxtaposing the minor grooves. Cross-strand H-bonds form at the interface, accepted by the conserved adenine nucleobases of the G:A basepairs. Alternative acceptors for one of these divides the k-turns into two conformational classes N3 and N1. The base pair that follows the G:A pairs (3b:3n) determines which conformation is adopted by a given k-turn. k-turns often mediate tertiary contacts in folded RNA species and frequently bind proteins. Common k-turn binding proteins include members of the L7Ae family, such as the human 15·5k protein. A recognition helix within these proteins binds in the widened major groove on the outside of the k-turn, that makes specific H-bonds with the conserved guanine nucleobases of the G:A pairs. L7Ae binds with extremely high affinity, and single-molecule data are consistent with folding by conformational selection. The standard, simple k-turn can be elaborated in a variety of ways, that include the complex k-turns and the k-junctions. In free solution in the absence of added metal ions or protein k-turns do not adopt the tightly-kinked conformation. They undergo folding by the binding of proteins, by the formation of tertiary contacts, and some (but not all) will fold on the addition of metal ions. Whether or not folding occurs in the presence of metal ions depends on local sequence, including the 3b:3n position, and the -1b:-1n position (5' to the bulge). In most cases -1b:-1n = C:G, so that the 3b:3n position is critical since it determines both folding properties and conformation. In general, the selection of these sequence matches a given k-turn to its biological requirements. The k-turn structure is now very well understood, to the point at which they can be used as a building block for the formation of RNA nano-objects, including triangles and squares.

Finely tuned conformational dynamics regulate the protective function of the lncRNA MALAT1 triple helix

Nucleic acid triplexes may regulate many important biological processes. Persistent accumulation of the oncogenic 7-kb long noncoding RNA MALAT1 is dependent on an unusually long intramolecular triple helix. This triplex structure is positioned within a conserved ENE (element for nuclear expression) motif at the lncRNA 3′ terminus and protects the entire transcript from degradation in a polyA-independent manner. A requisite 3′ maturation step leads to triplex formation though the precise mechanism of triplex folding remains unclear. Furthermore, the contributions of several peripheral structural elements to triplex formation and protective function have not been determined. We evaluated the stability, conformational fluctuations, and function of this MALAT1 ENE triple helix (M1^TH) protective element using in vitro mutational analyses coupled with biochemical and biophysical characterizations. Using fluorescence and UV melts, FRET, and an exonucleolytic decay assay we define a concerted mechanism for triplex formation and uncover a metastable, dynamic triplex population under near-physiological conditions. Structural elements surrounding the triplex regulate the dynamic M1^TH conformational variability, but increased triplex dynamics lead to M1^TH degradation. Taken together, we suggest that finely tuned dynamics may be a general mechanism regulating triplex-mediated functions.

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.

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.

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.

Rapid RNA-ligand interaction analysis through high-information content conformational and stability landscapes

The structure and biological properties of RNAs are a function of changing cellular conditions, but comprehensive, simultaneous investigation of the effect of multiple interacting environmental variables is not easily achieved. We have developed an efficient, high-throughput method to characterize RNA structure and thermodynamic stability as a function of multiplexed solution conditions using Förster resonance energy transfer (FRET). In a single FRET experiment using conventional quantitative PCR instrumentation, 19,400 conditions of MgCl₂, ligand and temperature are analysed to generate detailed empirical conformational and stability landscapes of the cyclic diguanylate (c-di-GMP) riboswitch. The method allows rapid comparison of RNA structure modulation by cognate and non-cognate ligands. Landscape analysis reveals that kanamycin B stabilizes a non-native, idiosyncratic conformation of the riboswitch that inhibits c-di-GMP binding. This demonstrates that allosteric control of folding, rather than direct competition with cognate effectors, is a viable approach for pharmacologically targeting riboswitches and other structured RNA molecules.

The structure and biological properties of RNAs are a function of changing cellular conditions. Here, Baird et al. report a high-throughput Förster resonance energy transfer (FRET) method to rapidly compare RNA structure modulation by cognate and non-cognate ligands across multiplexed solution conditions.

Recognition of the bacterial alarmone ZMP through long-distance association of two RNA subdomains

The bacterial alarmone 5-aminoimidazole-4-carboxamide riboside 5'-triphosphate (ZTP), derived from the monophosphorylated purine precursor ZMP, accumulates during folate starvation. ZTP regulates genes involved in purine and folate metabolism through a cognate riboswitch. The linker connecting this riboswitch’s two sub-domains varies in length by over 100 nucleotides. We report the co-crystal structure of the Fusobacterium ulcerans riboswitch bound to ZMP, which spans the two sub-domains whose interface also comprises a pseudoknot and ribose zipper. The riboswitch recognizes the carboxamide oxygen of ZMP through an unprecedented inner-sphere coordination with a Mg²⁺ ion. We demonstrate that the affinity of the riboswitch for ZMP is modulated by the linker length. Notably, ZMP can bind to the two sub-domains together even when synthesized as separate RNAs. The ZTP riboswitch demonstrates how specific small-molecule binding can drive association of distant non-coding RNA domains to regulate gene expression.

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.

Crystal structure of a c-di-AMP riboswitch reveals an internally pseudo-dimeric RNA

Cyclic diadenosine monophosphate (c-di-AMP) is a second messenger that is essential for growth and homeostasis in bacteria. A recently discovered c-di-AMP-responsive riboswitch controls the expression of genes in a variety of bacteria, including important pathogens. To elucidate the molecular basis for specific binding of c-di-AMP by a gene-regulatory mRNA domain, we have determined the co-crystal structure of this riboswitch. Unexpectedly, the structure reveals an internally pseudo-symmetric RNA in which two similar three-helix-junction elements associate head-to-tail, creating a trough that cradles two c-di-AMP molecules making quasi-equivalent contacts with the riboswitch. The riboswitch selectively binds c-di-AMP and discriminates exquisitely against other cyclic dinucleotides, such as c-di-GMP and cyclic-AMP-GMP, via interactions with both the backbone and bases of its cognate second messenger. Small-angle X-ray scattering experiments indicate that global folding of the riboswitch is induced by the two bound cyclic dinucleotides, which bridge the two symmetric three-helix domains. This structural reorganization likely couples c-di-AMP binding to gene expression.

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.

A critical base pair in k-turns that confers folding characteristics and correlates with biological function

Kink turns (k-turns) are widespread elements in RNA that mediate tertiary contacts by kinking the helical axis. We have found that the ability of k-turns to undergo ion-induced folding is conferred by a single base pair that follows the conserved A·G pairs, that is, the 3b·3n position. A Watson–Crick pair leads to an inability to fold in metal ions alone, while 3n=G or 3b=C (but not both) permits folding. Crystallographic study reveals two hydrated metal ions coordinated to O6 of G3n and G2n of Kt-7. Removal of either atom impairs Mg²⁺-induced folding in solution. While SAM-I riboswitches have 3b·3n sequences that would predispose them to ion-induced folding, U4 snRNA are strongly biased to an inability to such folding. Thus riboswitch sequences allow folding to occur independently of protein binding, while U4 should remain unfolded until bound by protein. The empirical rules deduced for k-turn folding have strong predictive value.

The k-turn is a widespread RNA element that adopts a kinked structure that mediates tertiary contacts and frequently binds specific proteins. Here, McPhee et al. show that the ability of a given k-turn to fold in the presence of metal ions alone—or to otherwise require protein binding—is attributable to a specific base pair.

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.

Primary transcriptome map of the hyperthermophilic archaeon Thermococcus kodakarensis

Background

Prokaryotes have relatively small genomes, densely-packed with protein-encoding sequences. RNA sequencing has, however, revealed surprisingly complex transcriptomes and here we report the transcripts present in the model hyperthermophilic Archaeon, Thermococcus kodakarensis, under different physiological conditions.

Results

Sequencing cDNA libraries, generated from RNA isolated from cells under different growth and metabolic conditions has identified >2,700 sites of transcription initiation, established a genome-wide map of transcripts, and consensus sequences for transcription initiation and post-transcription regulatory elements. The primary transcription start sites (TSS) upstream of 1,254 annotated genes, plus 644 primary TSS and their promoters within genes, are identified. Most mRNAs have a 5'-untranslated region (5'-UTR) 10 to 50 nt long (median = 16 nt), but ~20% have 5'-UTRs from 50 to 300 nt long and ~14% are leaderless. Approximately 50% of mRNAs contain a consensus ribosome binding sequence. The results identify TSS for 1,018 antisense transcripts, most with sequences complementary to either the 5'- or 3'-region of a sense mRNA, and confirm the presence of transcripts from all three CRISPR loci, the RNase P and 7S RNAs, all tRNAs and rRNAs and 69 predicted snoRNAs. Two putative riboswitch RNAs were present in growing but not in stationary phase cells. The procedure used is designed to identify TSS but, assuming that the number of cDNA reads correlates with transcript abundance, the results also provide a semi-quantitative documentation of the differences in T. kodakarensis genome expression under different growth conditions and confirm previous observations of substrate-dependent specific gene expression. Many previously unanticipated small RNAs have been identified, some with relative low GC contents (≤50%) and sequences that do not fold readily into base-paired secondary structures, contrary to the classical expectations for non-coding RNAs in a hyperthermophile.

Conclusion

The results identify >2,700 TSS, including almost all of the primary sites of transcription initiation upstream of annotated genes, plus many secondary sites, sites within genes and sites resulting in antisense transcripts. The T. kodakarensis genome is small (~2.1 Mbp) and tightly packed with protein-encoding genes, but the transcriptomes established also contain many non-coding RNAs and predict extensive RNA-based regulation in this model Archaeon.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-684) contains supplementary material, which is available to authorized users.

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.

The K-turn motif in riboswitches and other RNA species

The kink turn is a widespread structure motif that introduces a tight bend into the axis of duplex RNA. This generally functions to mediate tertiary interactions, and to serve as a specific protein binding site. K-turns or closely related structures are found in at least seven different riboswitch structures, where they function as key architectural elements that help generate the ligand binding pocket. This article is part of a Special Issue entitled: Riboswitches.

Global analysis of riboswitches by small-angle X-ray scattering and calorimetry

Riboswitches are phylogenetically widespread non-coding mRNA domains that directly bind cellular metabolites and regulate transcription, translation, RNA stability or splicing via alternative RNA structures modulated by ligand binding. The details of ligand recognition by many riboswitches have been elucidated using X-ray crystallography and NMR. However, the global dynamics of riboswitch-ligand interactions and their thermodynamic driving forces are less understood. By compiling the work of many laboratories investigating riboswitches using small-angle X-ray scattering (SAXS) and isothermal titration calorimetry (ITC), we uncover general trends and common themes. There is a pressing need for community-wide consensus experimental conditions to allow results of riboswitch studies to be compared rigorously. Nonetheless, our meta-analysis reveals considerable diversity in the extent to which ligand binding reorganizes global riboswitch structures. It also demonstrates a wide spectrum of enthalpy-entropy compensation regimes across riboswitches that bind a diverse set of ligands, giving rise to a relatively narrow range of physiologically relevant free energies and ligand affinities. From the strongly entropy-driven binding of glycine to the predominantly enthalpy-driven binding of c-di-GMP to their respective riboswitches, these distinct thermodynamic signatures reflect the versatile strategies employed by RNA to adapt to the chemical natures of diverse ligands. Riboswitches have evolved to use a combination of long-range tertiary interactions, conformational selection, and induced fit to work with distinct ligand structure, charge, and solvation properties. This article is part of a Special Issue entitled: Riboswitches.

New molecular engineering approaches for crystallographic studies of large RNAs

Crystallization of RNAs with complex three-dimensional architectures remains a formidable experimental challenge. We review a number of successful heuristics involving engineering of the target RNAs to facilitate crystal contact formation, such as those that enabled the crystallization and structure determination of the cognate tRNA complexes of RNase P holoenzyme and the Stem I domain of the T-box riboswitch. Recently, RNA-targeted antibody Fab fragments and Kink-turn binding proteins have joined the ranks of successful chaperones for RNA crystallization. Lastly, we review the use of structured RNAs to facilitate crystallization of RNA-binding proteins and other RNAs.

Structure of a rare non-standard sequence k-turn bound by L7Ae protein

Kt-23 from Thelohania solenopsae is a rare RNA kink turn (k-turn) where an adenine replaces the normal guanine at the 2n position. L7Ae is a member of a strongly conserved family of proteins that bind a range of k-turn structures in the ribosome, box C/D and H/ACA small nucleolar RNAs and U4 small nuclear RNA. We have solved the crystal structure of T. solenopsae Kt-23 RNA bound to Archeoglobus fulgidus L7Ae protein at a resolution of 2.95 Å. The protein binds in the major groove displayed on the outer face of the k-turn, in a manner similar to complexes with standard k-turn structures. The k-turn adopts a standard N3 class conformation, with a single hydrogen bond from A2b N6 to A2n N3. This contrasts with the structure of the same sequence located in the SAM-I riboswitch, where it adopts an N1 structure, showing the inherent plasticity of k-turn structure. This potentially can affect any tertiary interactions in which the RNA participates.

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.

Co-crystal structure of a T-box riboswitch stem I domain in complex with its cognate tRNA

In Gram-positive bacteria, T-box riboswitches regulate the expression of aminoacyl-tRNA synthetases and other proteins in response to fluctuating transfer RNA aminoacylation levels under various nutritional states. T-boxes reside in the 5'-untranslated regions of the messenger RNAs they regulate, and consist of two conserved domains. Stem I contains the specifier trinucleotide that base pairs with the anticodon of cognate tRNA. 3' to stem I is the antiterminator domain, which base pairs with the tRNA acceptor end and evaluates its aminoacylation state. Despite high phylogenetic conservation and widespread occurrence in pathogens, the structural basis of tRNA recognition by this riboswitch remains ill defined. Here we demonstrate that the ~100-nucleotide T-box stem I is necessary and sufficient for specific, high-affinity (dissociation constant (Kd) ~150 nM) tRNA binding, and report the structure of Oceanobacillus iheyensis glyQ stem I in complex with its cognate tRNA at 3.2 Å resolution. Stem I recognizes the overall architecture of tRNA in addition to its anticodon, something accomplished by large ribonucleoproteins such as the ribosome, or proteins such as aminoacyl-tRNA synthetases, but is unprecedented for a compact mRNA domain. The C-shaped stem I cradles the L-shaped tRNA, forming an extended (1,604 Å(2)) intermolecular interface. In addition to the specifier-anticodon interaction, two interdigitated T-loops near the apex of stem I stack on the tRNA elbow in a manner analogous to those of the J11/12-J12/11 motif of RNase P and the L1 stalk of the ribosomal E-site. Because these ribonucleoproteins and T-boxes are unrelated, this strategy to recognize a universal tRNA feature probably evolved convergently. Mutually induced fit of stem I and the tRNA exploiting the intrinsic flexibility of tRNA and its conserved post-transcriptional modifications results in high shape complementarity, which in addition to providing specificity and affinity, globally organizes the T-box to orchestrate tRNA-dependent transcription regulation.

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.