Structural basis for recognition of S-adenosylhomocysteine by riboswitches

Modeling riboswitches: The impact of concentration of Mg2+ ions on the folding of the SAM-II riboswitch

Riboswitches are regulatory elements present in bacterial messenger RNA acting as sensors of small molecules and consequently playing a vital role in bacterial gene regulation. The SAM-II riboswitch is a class of riboswitches that recognizes S-adenosyl methionine. It has been previously shown that the presence of Mg2+ ions stabilizes the preexisting minor state of the riboswitch, which is structurally characterized having a nucleated pseudoknot, leading to the increase of its probability. In this study, an analytical equilibrium model is developed to describe the impact of Mg2+ ion concentration on the folding of the SAM-II riboswitch, thus linking RNA folding and tertiary interactions energetics to ligand binding. The model enables quantitative predictions of equilibrium folding intermediates as a function of concentration of Mg2+ ions.

A microscale thermophoresis-based enzymatic RNA methyltransferase assay enables the discovery of DNMT2 inhibitors

RNA methyltransferases (MTases) have recently become increasingly important in drug discovery. Yet, most frequently utilized RNA MTase assays are limited in their throughput and hamper this rapidly evolving field of medicinal chemistry. This study developed a microscale thermophoresis (MST)-based split aptamer assay for enzymatic MTase investigations, improving current methodologies by offering a non-proprietary, cost-effective, and highly sensitive approach. Our findings demonstrate the assay's effectiveness across different RNA MTases, including inhibitor characterization of METTL3/14, DNMT2, NSUN2, and S. aureus TrmD, enabling future drug discovery efforts. Using this concept, a pilot screening on the cancer drug target DNMT2 discovered several hit compounds with micromolar potency.

A general strategy for engineering GU base pairs to facilitate RNA crystallization

X-ray crystallography is a fundamental technique that provides atomic-level insights into RNA structures. However, obtaining crystals of RNA structures diffracting to high resolution is challenging. We introduce a simple strategy to enhance the resolution limit of RNA crystals by the selective substitution of Watson-Crick pairs by GU pairs within RNA sequences. Our approach has successfully yielded high-resolution structures for eight unique RNA crystals. Notably, six instances showed marked resolution enhancement upon GC/AU to GU base pair substitution, with two cases achieving high-resolution structures from initially poor data. In one case, reverting GU to GC base pairs also improved resolution. Our method facilitated the first structural determinations of the Long Interspersed Nuclear Element-1 and Olfactory Receptor family 4 subfamily K member 15 ribozymes, the 2'-deoxyguanosine-III riboswitch and the Broccoli RNA aptamer. The placement of GU base pairs within the first 5' helical stem of any given RNA species, or in one peripheral stem, is shown to be sufficient. These results offer a simple and effective approach for designing sequences or selecting sequences from homologous sequences, for high-resolution RNA structure determination.

Knotty is nice: Metabolite binding and RNA-mediated gene regulation by the preQ1 riboswitch family

Riboswitches sense specific cellular metabolites, leading to messenger RNA conformational changes that regulate downstream genes. Here, we review the three known prequeosine1 (preQ1) riboswitch classes, which encompass five gene-regulatory motifs derived from distinct consensus models of folded RNA pseudoknots. Structural and functional analyses reveal multiple gene-regulation strategies ranging from partial occlusion of the ribosome-binding Shine-Dalgarno sequence (SDS), SDS sequestration driven by kinetic or thermodynamic folding pathways, direct preQ1 recognition by the SDS, and complete SDS burial with in the riboswitch architecture. Family members can also induce elemental transcriptional pausing, which depends on ligand-mediated pseudoknot formation. Accordingly, preQ1 family members provide insight into a wide range of gene-regulatory tactics as well as a diverse repertoire of chemical approaches used to recognize the preQ1 metabolite. From a broader perspective, future challenges for the field will include the identification of new riboswitches in mRNAs that do not possess an SDS or those that induce ligand-dependent transcriptional pausing. When choosing an antibacterial target, the field must also consider how well a riboswitch accommodates mutations. Investigation of riboswitches in their natural context will also be critical to elucidate how RNA-mediated gene regulation influences organism fitness, thus providing a firm foundation for antibiotic development.

Riboswitch Mechanisms for Regulation of P1 Helix Stability

Riboswitches are highly structured RNA regulators of gene expression. Although found in all three domains of life, they are particularly abundant and widespread in bacteria, including many human pathogens, thus making them an attractive target for antimicrobial development. Moreover, the functional versatility of riboswitches to recognize a myriad of ligands, including ions, amino acids, and diverse small-molecule metabolites, has enabled the generation of synthetic aptamers that have been used as molecular probes, sensors, and regulatory RNA devices. Generally speaking, a riboswitch consists of a ligand-sensing aptamer domain and an expression platform, whose genetic control is achieved through the formation of mutually exclusive secondary structures in a ligand-dependent manner. For most riboswitches, this involves formation of the aptamer's P1 helix and the regulation of its stability, whose competing structure turns gene expression ON/OFF at the level of transcription or translation. Structural knowledge of the conformational changes involving the P1 regulatory helix, therefore, is essential in understanding the structural basis for ligand-induced conformational switching. This review provides a summary of riboswitch cases for which ligand-free and ligand-bound structures have been determined. Comparative analyses of these structures illustrate the uniqueness of these riboswitches, not only in ligand sensing but also in the various structural mechanisms used to achieve the same end of regulating switch helix stability. In all cases, the ligand stabilizes the P1 helix primarily through coaxial stacking interactions that promote helical continuity.

Thermodynamic compensation to temperature extremes in B. subtilis vs T. maritima lysine riboswitches

T. maritima and B. subtilis are bacteria that inhabit significantly different thermal environments, ∼80 vs. ∼40°C, yet employ similar lysine riboswitches to aid in the transcriptional regulation of the genes involved in the synthesis and transport of amino acids. Despite notable differences in G-C basepair frequency and primary sequence, the aptamer moieties of each riboswitch have striking similarities in tertiary structure, with several conserved motifs and long-range interactions. To explore genetic adaptation in extreme thermal environments, we compare the kinetic and thermodynamic behaviors in T. maritima and B. subtilis lysine riboswitches via single-molecule fluorescence resonance energy transfer analysis. Kinetic studies reveal that riboswitch folding rates increase with lysine concentration while the unfolding rates are independent of lysine. This indicates that both riboswitches bind lysine through an induced-fit ("bind-then-fold") mechanism, with lysine binding necessarily preceding conformational changes. Temperature-dependent van't Hoff studies reveal qualitative similarities in the thermodynamic landscapes for both riboswitches in which progression from the open, lysine-unbound state to both transition states (‡) and closed, lysine-bound conformations is enthalpically favored yet entropically penalized, with comparisons of enthalpic and entropic contributions extrapolated to a common [K+] = 100 mM in quantitative agreement. Finally, temperature-dependent Eyring analysis reveals the TMA and BSU riboswitches to have remarkably similar folding/unfolding rate constants when extrapolated to their respective (40 and 80°C) environmental temperatures. Such behavior suggests a shared strategy for ligand binding and aptamer conformational change in the two riboswitches, based on thermodynamic adaptations in number of G-C basepairs and/or modifications in tertiary structure that stabilize the ligand-unbound conformation to achieve biocompetence under both hyperthermophilic and mesothermophilic conditions.

Development and validation of a generic methyltransferase enzymatic assay based on an SAH riboswitch

Methylation of proteins and nucleic acids plays a fundamental role in epigenetic regulation, and discovery of methyltransferase (MT) inhibitors is an area of intense activity. Because of the diversity of MTs and their products, assay methods that detect S-adenosylhomocysteine (SAH) - the invariant product of S-adenosylmethionine (SAM)-dependent methylation reactions - offer some advantages over methods that detect specific methylation events. However, direct, homogenous detection of SAH requires a reagent capable of discriminating between SAH and SAM, which differ by a single methyl group. Moreover, MTs are slow enzymes and many have submicromolar affinities for SAM; these properties translate to a need for detection of SAH at low nanomolar concentrations in the presence of excess SAM. To meet these needs, we leveraged the exquisite molecular recognition properties of a naturally occurring SAH-sensing RNA aptamer, or riboswitch. By splitting the riboswitch into two fragments, such that SAH binding induces assembly of a trimeric complex, we engineered sensors that transduce binding of SAH into positive fluorescence polarization (FP) and time resolved Förster resonance energy transfer (TR-FRET) signals. The split riboswitch configuration, called the AptaFluor™ SAH Methyltransferase Assay, allows robust detection of SAH (Z' > 0.7) at concentrations below 10 nM, with overnight signal stability in the presence of typical MT assay components. The AptaFluor assay tolerates diverse MT substrates, including histones, nucleosomes, DNA and RNA, and we demonstrated its utility as a robust, enzymatic assay method for several methyltransferases with SAM Km values < 1 µM. The assay was validated for HTS by performing a pilot screen of 1,280 compounds against the SARS-CoV-2 RNA capping enzyme, nsp14. By enabling direct, homogenous detection of SAH at low nanomolar concentrations, the AptaFluor assay provides a universal platform for screening and profiling MTs at physiologically relevant SAM concentrations.

Structure-based insights into recognition and regulation of SAM-sensing riboswitches

Riboswitches are highly conserved RNA elements that located in the 5'-UTR of mRNAs, which undergo real-time structure conformational change to achieve the regulation of downstream gene expression by sensing their cognate ligands. S-adenosylmethionine (SAM) is a ubiquitous methyl donor for transmethylation reactions in all living organisms. SAM riboswitch is one of the most abundant riboswitches that bind to SAM with high affinity and selectivity, serving as regulatory modules in multiple metabolic pathways. To date, seven SAM-specific riboswitch classes that belong to four families, one SAM/SAH riboswitch and one SAH riboswitch have been identified. Each SAM riboswitch family has a well-organized tertiary core scaffold to support their unique ligand-specific binding pocket. In this review, we summarize the current research progress on the distribution, structure, ligand recognition and gene regulation mechanism of these SAM-related riboswitch families, and further discuss their evolutionary prospects and potential applications.

Decoding the Identification Mechanism of an SAM-III Riboswitch on Ligands through Multiple Independent Gaussian-Accelerated Molecular Dynamics Simulations

S-Adenosyl-l-methionine (SAM)-responsive riboswitches play a central role in the regulation of bacterial gene expression at the level of transcription attenuation or translation inhibition. In this study, multiple independent Gaussian-accelerated molecular dynamics simulations were performed to decipher the identification mechanisms of SAM-III (S_MK) on ligands SAM, SAH, and EEM. The results reveal that ligand binding highly affects the structural flexibility, internal dynamics, and conformational changes of SAM-III. The dynamic analysis shows that helices P3 and P4 as well as two junctions J23 and J24 of SAM-III are highly susceptible to ligand binding. Analyses of free energy landscapes suggest that ligand binding induces different free energy profiles of SAM-III, which leads to the difference in identification sites of SAM-III on ligands. The information on ligand-nucleotide interactions not only uncovers that the π-π, cation-π, and hydrogen bonding interactions drive identification of SAM-III on the three ligands but also reveals that different electrostatic properties of SAM, SAH, and EEM alter the active sites of SAM-III. Meanwhile, the results also verify that the adenine group of SAM, SAH, and EEM is well recognized by conserved nucleotides G7, A29, U37, A38, and G48. We expect that this study can provide useful information for understanding the applications of SAM-III in chemical, synthetic RNA biology, and biomedical fields.

Lysine-Dependent Entropy Effects in the B. subtilis Lysine Riboswitch: Insights from Single-Molecule Thermodynamic Studies

Riboswitches play an important role in RNA-based sensing/gene regulation control for many bacteria. In particular, the accessibility of multiple conformational states at physiological temperatures allows riboswitches to selectively bind a cognate ligand in the aptamer domain, which triggers secondary structural changes in the expression platform, and thereby "switching" between on or off transcriptional or translational states for the downstream RNA. The present work exploits temperature-controlled, single-molecule total internal reflection fluorescence (TIRF) microscopy to study the thermodynamic landscape of such ligand binding/folding processes, specifically for the Bacillus subtilis lysine riboswitch. The results confirm that the riboswitch folds via an induced-fit (IF) mechanism, in which cognate lysine ligand first binds to the riboswitch before structural rearrangement takes place. The transition state to folding is found to be enthalpically favored (ΔH_fold^‡ < 0), yet with a free-energy barrier that is predominantly entropic (-TΔS_fold^‡ > 0), which results in folding (unfolding) rate constants strongly dependent (independent) of lysine concentration. Analysis of the single-molecule kinetic "trajectories" reveals this rate constant dependence of k_fold on lysine to be predominantly entropic in nature, with the additional lysine conferring preferential advantage to the folding process by the presence of ligands correctly oriented with respect to the riboswitch platform. By way of contrast, van't Hoff analysis reveals enthalpic contributions to the overall folding thermodynamics (ΔH⁰) to be surprisingly constant and robustly independent of lysine concentration. The results demonstrate the crucial role of hydrogen bonding between the ligand and riboswitch platform but with only a relatively modest fraction (45%) of the overall enthalpy change needed to access the transition state and initiate transcriptional switching.

Structures of artificially designed discrete RNA nanoarchitectures at near-atomic resolution

Although advances in nanotechnology have enabled the construction of complex and functional synthetic nucleic acid–based nanoarchitectures, high-resolution discrete structures are lacking because of the difficulty in obtaining good diffracting crystals. Here, we report the design and construction of RNA nanostructures based on homooligomerizable one-stranded tiles for x-ray crystallographic determination. We solved three structures to near-atomic resolution: a 2D parallelogram, a 3D nanobracelet unexpectedly formed from an RNA designed for a nanocage, and, eventually, a bona fide 3D nanocage designed with the guidance of the two previous structures. Structural details of their constituent motifs, such as kissing loops, branched kissing loops, and T-junctions, that resemble natural RNA motifs and resisted x-ray determination are revealed, providing insights into those natural motifs. This work unveils the largely unexplored potential of crystallography in gaining high-resolution feedback for nanoarchitectural design and suggests a route to investigate RNA motif structures by configuring them into nanoarchitectures.

RNA stabilization by a poly(A) tail 3'-end binding pocket and other modes of poly(A)-RNA interaction

Polyadenylate [poly(A)] tail addition to the 3' end of a wide range of RNAs is a highly conserved modification that plays a central role in cellular RNA function. Elements for nuclear expression (ENEs) are cis-acting RNA elements that stabilize poly(A) tails by sequestering them in RNA triplex structures. A crystal structure of a double ENE from a rice hAT transposon messenger RNA complexed with poly(A)₂₈ at a resolution of 2.89 angstroms reveals multiple modes of interaction with poly(A), including major-groove triple helices, extended minor-groove interactions with RNA double helices, a quintuple-base motif that transitions poly(A) from minor-groove associations to major-groove triple helices, and a poly(A) 3'-end binding pocket. Our findings both expand the repertoire of motifs involved in long-range RNA interactions and provide insights into how polyadenylation can protect an RNA's extreme 3' end.

Imaginary Ribozymes

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

Analysis of a preQ1-I riboswitch in effector-free and bound states reveals a metabolite-programmed nucleobase-stacking spine that controls gene regulation

Riboswitches are structured RNA motifs that recognize metabolites to alter the conformations of downstream sequences, leading to gene regulation. To investigate this molecular framework, we determined crystal structures of a preQ1-I riboswitch in effector-free and bound states at 2.00 Å and 2.65 Å-resolution. Both pseudoknots exhibited the elusive L2 loop, which displayed distinct conformations. Conversely, the Shine-Dalgarno sequence (SDS) in the S2 helix of each structure remained unbroken. The expectation that the effector-free state should expose the SDS prompted us to conduct solution experiments to delineate environmental changes to specific nucleobases in response to preQ1. We then used nudged elastic band computational methods to derive conformational-change pathways linking the crystallographically-determined effector-free and bound-state structures. Pathways featured: (i) unstacking and unpairing of L2 and S2 nucleobases without preQ1-exposing the SDS for translation and (ii) stacking and pairing L2 and S2 nucleobases with preQ1-sequestering the SDS. Our results reveal how preQ1 binding reorganizes L2 into a nucleobase-stacking spine that sequesters the SDS, linking effector recognition to biological function. The generality of stacking spines as conduits for effector-dependent, interdomain communication is discussed in light of their existence in adenine riboswitches, as well as the turnip yellow mosaic virus ribosome sensor.

Crystal structure and ligand-induced folding of the SAM/SAH riboswitch

While most SAM riboswitches strongly discriminate between SAM and SAH, the SAM/SAH riboswitch responds to both ligands with similar apparent affinities. We have determined crystal structures of the SAM/SAH riboswitch bound to SAH, SAM and other variant ligands at high resolution. The riboswitch forms an H-type pseudoknot structure with coaxial alignment of the stem–loop helix (P1) and the pseudoknot helix (PK). An additional three base pairs form at the non-open end of P1, and the ligand is bound at the interface between the P1 extension and the PK helix. The adenine nucleobase is stacked into the helix and forms a trans Hoogsteen–Watson–Crick base pair with a uridine, thus becoming an integral part of the helical structure. The majority of the specific interactions are formed with the adenosine. The methionine or homocysteine chain lies in the groove making a single hydrogen bond, and there is no discrimination between the sulfonium of SAM or the thioether of SAH. Single-molecule FRET analysis reveals that the riboswitch exists in two distinct conformations, and that addition of SAM or SAH shifts the population into a stable state that likely corresponds to the form observed in the crystal. A model for translational regulation is presented whereby in the absence of ligand the riboswitch is largely unfolded, lacking the PK helix so that translation can be initiated at the ribosome binding site. But the presence of ligand stabilizes the folded conformation that includes the PK helix, so occluding the ribosome binding site and thus preventing the initiation of translation.

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.

A unified dinucleotide alphabet describing both RNA and DNA structures

By analyzing almost 120 000 dinucleotides in over 2000 nonredundant nucleic acid crystal structures, we define 96+1 diNucleotide Conformers, NtCs, which describe the geometry of RNA and DNA dinucleotides. NtC classes are grouped into 15 codes of the structural alphabet CANA (Conformational Alphabet of Nucleic Acids) to simplify symbolic annotation of the prominent structural features of NAs and their intuitive graphical display. The search for nontrivial patterns of NtCs resulted in the identification of several types of RNA loops, some of them observed for the first time. Over 30% of the nearly six million dinucleotides in the PDB cannot be assigned to any NtC class but we demonstrate that up to a half of them can be re-refined with the help of proper refinement targets. A statistical analysis of the preferences of NtCs and CANA codes for the 16 dinucleotide sequences showed that neither the NtC class AA00, which forms the scaffold of RNA structures, nor BB00, the DNA most populated class, are sequence neutral but their distributions are significantly biased. The reported automated assignment of the NtC classes and CANA codes available at dnatco.org provides a powerful tool for unbiased analysis of nucleic acid structures by structural and molecular biologists.

Exploring the RNA-Recognition Mechanism Using Supervised Molecular Dynamics (SuMD) Simulations: Toward a Rational Design for Ribonucleic-Targeting Molecules?

Although proteins have represented the molecular target of choice in the development of new drug candidates, the pharmaceutical importance of ribonucleic acids has gradually been growing. The increasing availability of structural information has brought to light the existence of peculiar three-dimensional RNA arrangements, which can, contrary to initial expectations, be recognized and selectively modulated through small chemical entities or peptides. The application of classical computational methodologies, such as molecular docking, for the rational development of RNA-binding candidates is, however, complicated by the peculiarities characterizing these macromolecules, such as the marked conformational flexibility, the singular charges distribution, and the relevant role of solvent molecules. In this work, we have thus validated and extended the applicability domain of SuMD, an all-atoms molecular dynamics protocol that allows to accelerate the sampling of molecular recognition events on a nanosecond timescale, to ribonucleotide targets of pharmaceutical interest. In particular, we have proven the methodological ability by reproducing the binding mode of viral or prokaryotic ribonucleic complexes, as well as that of artificially engineered aptamers, with an impressive degree of accuracy.

In silico prediction and expression profile analysis of small non-coding RNAs in Herbaspirillum seropedicae SmR1

BACKGROUND

Herbaspirillum seropedicae is a diazotrophic bacterium from the β-proteobacteria class that colonizes endophytically important gramineous species, promotes their growth through phytohormone-dependent stimulation and can express nif genes and fix nitrogen inside plant tissues. Due to these properties this bacterium has great potential as a commercial inoculant for agriculture. The H. seropedicae SmR1 genome is completely sequenced and annotated but despite the availability of diverse structural and functional analysis of this genome, studies involving small non-coding RNAs (sRNAs) has not yet been done. We have conducted computational prediction and RNA-seq analysis to select and confirm the expression of sRNA genes in the H. seropedicae SmR1 genome, in the presence of two nitrogen independent sources and in presence of naringenin, a flavonoid secreted by some plants.

RESULTS

This approach resulted in a set of 117 sRNAs distributed in riboswitch, cis-encoded and trans-encoded categories and among them 20 have Rfam homologs. The housekeeping sRNAs tmRNA, ssrS and 4.5S were found and we observed that a large number of sRNAs are more expressed in the nitrate condition rather than the control condition and in the presence of naringenin. Some sRNAs expression were confirmed in vitro and this work contributes to better understand the post transcriptional regulation in this bacterium.

CONCLUSIONS

H. seropedicae SmR1 express sRNAs in the presence of two nitrogen sources and/or in the presence of naringenin. The functions of most of these sRNAs remains unknown but their existence in this bacterium confirms the evidence that sRNAs are involved in many different cellular activities to adapt to nutritional and environmental changes.

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

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

SAM-VI riboswitch structure and signature for ligand discrimination

Riboswitches are metabolite-sensing, conserved domains located in non-coding regions of mRNA that are central to regulation of gene expression. Here we report the first three-dimensional structure of the recently discovered S-adenosyl-L-methionine responsive SAM-VI riboswitch. SAM-VI adopts a unique fold and ligand pocket that are distinct from all other known SAM riboswitch classes. The ligand binds to the junctional region with its adenine tightly intercalated and Hoogsteen base-paired. Furthermore, we reveal the ligand discrimination mode of SAM-VI by additional X-ray structures of this riboswitch bound to S-adenosyl-L-homocysteine and a synthetic ligand mimic, in combination with isothermal titration calorimetry and fluorescence spectroscopy to explore binding thermodynamics and kinetics. The structure is further evaluated by analysis of ligand binding to SAM-VI mutants. It thus provides a thorough basis for developing synthetic SAM cofactors for applications in chemical and synthetic RNA biology.

Identification of Structural Motifs Using Networks of Hydrogen-Bonded Base Interactions in RNA Crystallographic Structures

RNA structural motifs can be identified using methods that analyze base–base interactions and the conformation of a structure’s backbone; however, these approaches do not necessarily take into consideration the hydrogen bonds that connect the bases or the networks of inter-connected hydrogen-bonded bases that are found in RNA structures. Large clusters of RNA bases that are tightly inter-connected by a network of hydrogen bonds are expected to be stable and relatively rigid substructures. Such base arrangements could therefore be present as structural motifs in RNA structures, especially when there is a requirement for a highly stable support platform or substructure to ensure the correct folding and spatial maintenance of functional sites that partake in catalysis or binding interactions. In order to test this hypothesis, we conducted a search in available RNA crystallographic structures in the Protein Data Bank database using queries that searched for profiles of bases inter-connected by hydrogen bonds. This method of searching does not require to have prior knowledge of the arrangement being searched. Our search results identified two clusters of six bases that are inter-connected by a network of hydrogen bonds. These arrangements of base sextuples have never been previously reported, thus making this the first report that proposes them as novel RNA tertiary motifs.

RNA-ALIS: Methodology for screening soluble RNAs as small molecule targets using ALIS affinity-selection mass spectrometry

Recent advances resulting from the completion of the human genome have shown that RNA has the promise to be a target for small molecule drugs, and therefore represents a previously unexploited class of targets for novel human therapeutics. We recently reported the adaptation of an affinity selection mass spectrometry screening technique, termed ALIS (Automatic Ligand Identification System), to screen and characterize a variety of RNA species from both prokaryotic and eukaryotic sources. We demonstrated that the ALIS technique, which had previously been used for protein targets, was also compatible for screening, ranking and characterizing small molecule ligands for RNA targets. We present here a detailed description of the use of ALIS for screening and characterizing ligands for RNA and discuss issues of validating and testing RNA for use in the ALIS system. We have also further elaborated on issues of RNA stability and testing in the ALIS system and demonstrate that the affinity-selection screening system has the potential to be a general solution for label-free screening and characterization of small molecule drug candidates for RNA targets.

Small-Molecule-Binding Riboswitches

RNA is a versatile biomolecule capable of transferring information, taking on distinct three-dimensional shapes, and reacting to ambient conditions. RNA molecules utilize a wide range of mechanisms to control gene expression. An example of such regulation is riboswitches. Consisting exclusively of RNA, they are able to control important metabolic processes, thus providing an elegant and efficient RNA-only regulation system. Existing across all domains of life, riboswitches appear to represent one of the most highly conserved mechanisms for the regulation of a broad range of biochemical pathways. Through binding of a wide range of small-molecule ligands to their so-called aptamer domain, riboswitches undergo a conformational change in their downstream "expression platform." In consequence, the pattern of gene expression changes, which in turn results in increased or decreased protein production. Riboswitches unite the sensing and transduction of a signal that can directly be coupled to the metabolism of the cell; thus they constitute a very potent regulatory mechanism for many organisms. Highly specific RNA-binding domains not only occur in vivo but can also be evolved by means of the SELEX (systematic evolution of ligands by exponential enrichment) method, which allows in vitro selection of aptamers against almost any ligand. Coupling of these aptamers with an expression platform has led to the development of synthetic riboswitches, a highly active research field of great relevance and immense potential. The aim of this review is to summarize developments in the riboswitch field over the last decade and address key questions of recent research.

Long-Range Interactions in Riboswitch Control of Gene Expression

Riboswitches are widespread RNA motifs that regulate gene expression in response to fluctuating metabolite concentrations. Known primarily from bacteria, riboswitches couple specific ligand binding and changes in RNA structure to mRNA expression in cis. Crystal structures of the ligand binding domains of most of the phylogenetically widespread classes of riboswitches, each specific to a particular metabolite or ion, are now available. Thus, the bound states-one end point-have been thoroughly characterized, but the unbound states have been more elusive. Consequently, it is less clear how the unbound, sensing riboswitch refolds into the ligand binding-induced output state. The ligand recognition mechanisms of riboswitches are diverse, but we find that they share a common structural strategy in positioning their binding sites at the point of the RNA three-dimensional fold where the residues farthest from one another in sequence meet. We review how riboswitch folds adhere to this fundamental strategy and propose future research directions for understanding and harnessing their ability to specifically control gene expression.

A magnesium-induced triplex pre-organizes the SAM-II riboswitch

Our 13C- and 1H-chemical exchange saturation transfer (CEST) experiments previously revealed a dynamic exchange between partially closed and open conformations of the SAM-II riboswitch in the absence of ligand. Here, all-atom structure-based molecular simulations, with the electrostatic effects of Manning counter-ion condensation and explicit magnesium ions are employed to calculate the folding free energy landscape of the SAM-II riboswitch. We use this analysis to predict that magnesium ions remodel the landscape, shifting the equilibrium away from the extended, partially unfolded state towards a compact, pre-organized conformation that resembles the ligand-bound state. Our CEST and SAXS experiments, at different magnesium ion concentrations, quantitatively confirm our simulation results, demonstrating that magnesium ions induce collapse and pre-organization. Agreement between theory and experiment bolsters microscopic interpretation of our simulations, which shows that triplex formation between helix P2b and loop L1 is highly sensitive to magnesium and plays a key role in pre-organization. Pre-organization of the SAM-II riboswitch allows rapid detection of ligand with high selectivity, which is important for biological function.

Rational design of a synthetic mammalian riboswitch as a ligand-responsive -1 ribosomal frame-shifting stimulator

Metabolite-responsive RNA pseudoknots derived from prokaryotic riboswitches have been shown to stimulate −1 programmed ribosomal frameshifting (PRF), suggesting −1 PRF as a promising gene expression platform to extend riboswitch applications in higher eukaryotes. However, its general application has been hampered by difficulty in identifying a specific ligand-responsive pseudoknot that also functions as a ligand-dependent -1 PRF stimulator. We addressed this problem by using the −1 PRF stimulation pseudoknot of SARS-CoV (SARS-PK) to build a ligand-dependent −1 PRF stimulator. In particular, the extra stem of SARS-PK was replaced by an RNA aptamer of theophylline and designed to couple theophylline binding with the stimulation of −1 PRF. Conformational and functional analyses indicate that the engineered theophylline-responsive RNA functions as a mammalian riboswitch with robust theophylline-dependent −1 PRF stimulation activity in a stable human 293T cell-line. Thus, RNA–ligand interaction repertoire provided by in vitro selection becomes accessible to ligand-specific −1 PRF stimulator engineering using SARS-PK as the scaffold for synthetic biology application.

An Adaptive Defect Weighted Sampling Algorithm to Design Pseudoknotted RNA Secondary Structures

Computational design of RNA sequences that fold into targeted secondary structures has many applications in biomedicine, nanotechnology and synthetic biology. An RNA molecule is made of different types of secondary structure elements and an important RNA element named pseudoknot plays a key role in stabilizing the functional form of the molecule. However, due to the computational complexities associated with characterizing pseudoknotted RNA structures, most of the existing RNA sequence designer algorithms generally ignore this important structural element and therefore limit their applications. In this paper we present a new algorithm to design RNA sequences for pseudoknotted secondary structures. We use NUPACK as the folding algorithm to compute the equilibrium characteristics of the pseudoknotted RNAs, and describe a new adaptive defect weighted sampling algorithm named Enzymer to design low ensemble defect RNA sequences for targeted secondary structures including pseudoknots. We used a biological data set of 201 pseudoknotted structures from the Pseudobase library to benchmark the performance of our algorithm. We compared the quality characteristics of the RNA sequences we designed by Enzymer with the results obtained from the state of the art MODENA and antaRNA. Our results show our method succeeds more frequently than MODENA and antaRNA do, and generates sequences that have lower ensemble defect, lower probability defect and higher thermostability. Finally by using Enzymer and by constraining the design to a naturally occurring and highly conserved Hammerhead motif, we designed 8 sequences for a pseudoknotted cis-acting Hammerhead ribozyme. Enzymer is available for download at https://bitbucket.org/casraz/enzymer.

In Vitro and In Vivo Enzyme Activity Screening via RNA-Based Fluorescent Biosensors for S-Adenosyl-l-homocysteine (SAH)

High-throughput enzyme activity screens are essential for target characterization and drug development, but few assays employ techniques or reagents that are applicable to both in vitro and live cell settings. Here, we present a class of selective and sensitive fluorescent biosensors for S-adenosyl-l-homocysteine (SAH) that provide a direct "mix and go" activity assay for methyltransferases (MTases), an enzyme class that includes several cancer therapeutic targets. Our riboswitch-based biosensors required an alternate inverted fusion design strategy, but retained full selectivity for SAH over its close structural analogue, the highly abundant methylation cofactor S-adenosyl-l-methionine (SAM). The level of ligand selectivity for these fluorescent biosensors exceeded that of commercial antibodies for SAH and proved critical to cellular applications, as we employed them to measure methylthioadenosine nucleosidase (MTAN) activity in live Escherichia coli. In particular, we were able to monitor in vivo increase of SAH levels upon chemical inhibition of MTAN using flow cytometry, which demonstrates high-throughput, single cell measurement of an enzyme activity associated with the biosynthesis of quorum sensing signal AI-2. Thus, this study presents RNA-based fluorescent biosensors as promising molecular reagents for high-throughput enzymatic assays that successfully bridge the gap between in vitro and in vivo applications.

Soaking Hexammine Cations into RNA Crystals to Obtain Derivatives for Phasing Diffraction Data

Solving a novel RNA structure by x-ray crystallography requires a means to obtain initial phase estimates. This is a challenge because many of the tools available for solving protein structures are not available for RNA. We have developed a reliable means to use hexammine cations to address this challenge. The process involves engineering the RNA to introduce a reliable hexammine binding site into the structure, then soaking crystals of these RNAs with an iridium (III) or cobalt (III) compound in a "directed soaking" strategy. Diffraction data obtained from these crystals then can be used in SAD or MAD phasing. In many cases, suitable derivatives can be obtained by soaking the hexammine into RNA crystals that have not been engineered. Considerations for using this method and example protocols are presented.

Recent advances and versatility of MAGE towards industrial applications

The genome engineering toolkit has expanded significantly in recent years, allowing us to study the functions of genes in cellular networks and assist in over-production of proteins, drugs, chemicals and biofuels. Multiplex automated genome engineering (MAGE) has been recently developed and gained more scientific interest towards strain engineering. MAGE is a simple, rapid and efficient tool for manipulating genes simultaneously in multiple loci, assigning genetic codes and integrating non-natural amino acids. MAGE can be further expanded towards the engineering of fast, robust and over-producing strains for chemicals, drugs and biofuels at industrial scales.

Rapid and high-throughput construction of microbial cell-factories with regulatory noncoding RNAs

Due to global crises such as pollution and depletion of fossil fuels, sustainable technologies based on microbial cell-factories have been garnering great interest as an alternative to chemical factories. The development of microbial cell-factories is imperative in cutting down the overall manufacturing cost. Thus, diverse metabolic engineering strategies and engineering tools have been established to obtain a preferred genotype and phenotype displaying superior productivity. However, these tools are limited to only a handful of genes with permanent modification of a genome and significant labor costs, and this is one of the bottlenecks associated with biofactory construction. Therefore, a groundbreaking rapid and high-throughput engineering tool is needed for efficient construction of microbial cell-factories. During the last decade, copious small noncoding RNAs (ncRNAs) have been discovered in bacteria. These are involved in substantial regulatory roles like transcriptional and post-transcriptional gene regulation by modulating mRNA elongation, stability, or translational efficiency. Because of their vulnerability, ncRNAs can be used as another layer of conditional control over gene expression without modifying chromosomal sequences, and hence would be a promising high-throughput tool for metabolic engineering. Here, we review successful design principles and applications of ncRNAs for high-throughput metabolic engineering or physiological studies of diverse industrially important microorganisms.

Secondary structural entropy in RNA switch (Riboswitch) identification

BACKGROUND

RNA regulatory elements play a significant role in gene regulation. Riboswitches, a widespread group of regulatory RNAs, are vital components of many bacterial genomes. These regulatory elements generally function by forming a ligand-induced alternative fold that controls access to ribosome binding sites or other regulatory sites in RNA. Riboswitch-mediated mechanisms are ubiquitous across bacterial genomes. A typical class of riboswitch has its own unique structural and biological complexity, making de novo riboswitch identification a formidable task. Traditionally, riboswitches have been identified through comparative genomics based on sequence and structural homology. The limitations of structural-homology-based approaches, coupled with the assumption that there is a great diversity of undiscovered riboswitches, suggests the need for alternative methods for riboswitch identification, possibly based on features intrinsic to their structure. As of yet, no such reliable method has been proposed.

RESULTS

We used structural entropy of riboswitch sequences as a measure of their secondary structural dynamics. Entropy values of a diverse set of riboswitches were compared to that of their mutants, their dinucleotide shuffles, and their reverse complement sequences under different stochastic context-free grammar folding models. Significance of our results was evaluated by comparison to other approaches, such as the base-pairing entropy and energy landscapes dynamics. Classifiers based on structural entropy optimized via sequence and structural features were devised as riboswitch identifiers and tested on Bacillus subtilis, Escherichia coli, and Synechococcus elongatus as an exploration of structural entropy based approaches. The unusually long untranslated region of the cotH in Bacillus subtilis, as well as upstream regions of certain genes, such as the sucC genes were associated with significant structural entropy values in genome-wide examinations.

CONCLUSIONS

Various tests show that there is in fact a relationship between higher structural entropy and the potential for the RNA sequence to have alternative structures, within the limitations of our methodology. This relationship, though modest, is consistent across various tests. Understanding the behavior of structural entropy as a fairly new feature for RNA conformational dynamics, however, may require extensive exploratory investigation both across RNA sequences and folding models.

Crystal structures of S-adenosylhomocysteine hydrolase from the thermophilic bacterium Thermotoga maritima

S-adenosylhomocysteine (SAH) hydrolase catalyzes the reversible hydrolysis of SAH into adenosine and homocysteine by using NAD(+) as a cofactor. The enzyme from Thermotoga maritima (tmSAHH) has great potentials in industrial applications because of its hyperthermophilic properties. Here, two crystal structures of tmSAHH in complex with NAD(+) show both open and closed conformations despite the absence of bound substrate. Each subunit of the tetrameric enzyme is composed of three domains, namely the catalytic domain, the NAD(+)-binding domain and the C-terminal domain. The NAD(+) binding mode is clearly observed and a substrate analogue can also be modeled into the active site, where two cysteine residues in mesophilic enzymes are replaced by serine and threonine in tmSAHH. Notably, the C-terminal domain of tmSAHH lacks the second loop region of mesophilic SAHH, which is important in NAD(+) binding, and thus exposes the bound cofactor to the solvent. The difference explains the higher NAD(+) requirement of tmSAHH because of the reduced affinity. Furthermore, the feature of missing loop is consistently observed in thermophilic bacterial and archaeal SAHHs, and may be related to their thermostability.

Single-molecule conformational dynamics of a biologically functional hydroxocobalamin riboswitch

Riboswitches represent a family of highly structured regulatory elements found primarily in the leader sequences of bacterial mRNAs. They function as molecular switches capable of altering gene expression; commonly, this occurs via a conformational change in a regulatory element of a riboswitch that results from ligand binding in the aptamer domain. Numerous studies have investigated the ligand binding process, but little is known about the structural changes in the regulatory element. A mechanistic description of both processes is essential for deeply understanding how riboswitches modulate gene expression. This task is greatly facilitated by studying all aspects of riboswitch structure/dynamics/function in the same model system. To this end, single-molecule fluorescence resonance energy transfer (smFRET) techniques have been used to directly observe the conformational dynamics of a hydroxocobalamin (HyCbl) binding riboswitch (env8HyCbl) with a known crystallographic structure.1 The single-molecule RNA construct studied in this work is unique in that it contains all of the structural elements both necessary and sufficient for regulation of gene expression in a biological context. The results of this investigation reveal that the undocking rate constant associated with the disruption of a long-range kissing-loop (KL) interaction is substantially decreased when the ligand is bound to the RNA, resulting in a preferential stabilization of the docked conformation. Notably, the formation of this tertiary KL interaction directly sequesters the Shine-Dalgarno sequence (i.e., the ribosome binding site) via base-pairing, thus preventing translation initiation. These results reveal that the conformational dynamics of this regulatory switch are quantitatively described by a four-state kinetic model, whereby ligand binding promotes formation of the KL interaction. The results of complementary cell-based gene expression experiments conducted in Escherichia coli are highly correlated with the smFRET results, suggesting that KL formation is directly responsible for regulating gene expression.

Folding of SAM-II riboswitch explored by replica-exchange molecular dynamics simulation

Riboswitches are cis-acting RNA fragments that function via a conformational transition mechanism when a specific target molecule binds to its binding pocket, representing an inviting new class of biomolecular target for the development of antibiotics. To understand the folding mechanism of SAM-II riboswitch, occurring predominantly in proteobacteria, a 100ns replica-exchange molecular dynamics simulation in explicit solvent is performed. Our results show that this RNA pseudoknot has multiple folding pathways, and various intermediate structures. The resultant riboswitch conformational transition map is well consistent with the recent fluorescence measurement, which confirms the dynamical properties of this pseudoknot. Moreover, a novel transition pathway is predicted. The global folding dynamics is mainly coupled with the helix rather than the loop region. The potential folding pathways of the riboswitch presented here should lead to a deeper understanding of the folding mechanism of the riboswitch, as well as the conformational change of RNA pseudoknot.

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.

Structure and function of pseudoknots involved in gene expression control

Natural RNA molecules can have a high degree of structural complexity but even the most complexly folded RNAs are assembled from simple structural building blocks. Among the simplest RNA elements are double-stranded helices that participate in the formation of different folding topologies and constitute the major fraction of RNA structures. One common folding motif of RNA is a pseudoknot, defined as a bipartite helical structure formed by base-pairing of the apical loop in the stem-loop structure with an outside sequence. Pseudoknots constitute integral parts of the RNA structures essential for various cellular activities. Among many functions of pseudoknotted RNAs is feedback regulation of gene expression, carried out through specific recognition of various molecules. Pseudoknotted RNAs autoregulate ribosomal and phage protein genes in response to downstream encoded proteins, while many metabolic and transport genes are controlled by cellular metabolites interacting with pseudoknotted RNA elements from the riboswitch family. Modulation of some genes also depends on metabolite-induced messenger RNA (mRNA) cleavage performed by pseudoknotted ribozymes. Several regulatory pseudoknots have been characterized biochemically and structurally in great detail. These studies have demonstrated a plethora of pseudoknot-based folds and have begun uncovering diverse molecular principles of the ligand-dependent gene expression control. The pseudoknot-mediated mechanisms of gene control and many unexpected and interesting features of the regulatory pseudoknots have significantly advanced our understanding of the genetic circuits and laid the foundation for modulation of their outcomes.

Structure, composition and metagenomic profile of soil microbiomes associated to agricultural land use and tillage systems in Argentine Pampas

Agriculture is facing a major challenge nowadays: to increase crop production for food and energy while preserving ecosystem functioning and soil quality. Argentine Pampas is one of the main world producers of crops and one of the main adopters of conservation agriculture. Changes in soil chemical and physical properties of Pampas soils due to different tillage systems have been deeply studied. Still, not much evidence has been reported on the effects of agricultural practices on Pampas soil microbiomes. The aim of our study was to investigate the effects of agricultural land use on community structure, composition and metabolic profiles on soil microbiomes of Argentine Pampas. We also compared the effects associated to conventional practices with the effects of no-tillage systems. Our results confirmed the impact on microbiome structure and composition due to agricultural practices. The phyla Verrucomicrobia, Plactomycetes, Actinobacteria, and Chloroflexi were more abundant in non cultivated soils while Gemmatimonadetes, Nitrospirae and WS3 were more abundant in cultivated soils. Effects on metabolic metagenomic profiles were also observed. The relative abundance of genes assigned to transcription, protein modification, nucleotide transport and metabolism, wall and membrane biogenesis and intracellular trafficking and secretion were higher in cultivated fertilized soils than in non cultivated soils. We also observed significant differences in microbiome structure and taxonomic composition between soils under conventional and no- tillage systems. Overall, our results suggest that agronomical land use and the type of tillage system have induced microbiomes to shift their life-history strategies. Microbiomes of cultivated fertilized soils (i.e. higher nutrient amendment) presented tendencies to copiotrophy while microbiomes of non cultivated homogenous soils appeared to have a more oligotrophic life-style. Additionally, we propose that conventional tillage systems may promote copiotrophy more than no-tillage systems by decreasing soil organic matter stability and therefore increasing nutrient availability.

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.

Insights into ligand binding to PreQ1 Riboswitch Aptamer from molecular dynamics simulations

Riboswitches play roles in transcriptional or translational regulation through specific ligand binding of their aptamer domains. Although a number of ligand-bound aptamer complex structures have been solved, it is important to know ligand-free conformations of the aptamers in order to understand the mechanism of specific binding by ligands. In this paper, preQ₁ riboswitch aptamer domain from Bacillus subtilis is studied by overall 1.5 μs all-atom molecular dynamics simulations We found that the ligand-free aptamer has a stable state with a folded P1-L3 and open binding pocket. The latter forms a cytosine-rich pool in which the nucleotide C19 oscillates between close and open positions, making it a potential conformation for preQ₁ entrance. The dynamic picture further suggests that the specific recognition of preQ₁ by the aptamer domain is not only facilitated by the key nucleotide C19 but also aided and enhanced by other cytosines around the binding pocket. These results should help to understand the details of preQ₁ binding.

Themes and variations in riboswitch structure and function

The complexity of gene expression control by non-coding RNA has been highlighted by the recent progress in the field of riboswitches. Discovered a decade ago, riboswitches represent a diverse group of non-coding mRNA regions that possess a unique ability to directly sense cellular metabolites and modulate gene expression through formation of alternative metabolite-free and metabolite-bound conformations. Such protein-free metabolite sensing domains utilize sophisticated three-dimensional folding of RNA molecules to discriminate between a cognate ligand from related compounds so that only the right ligand would trigger a genetic response. Given the variety of riboswitch ligands ranging from small cations to large coenzymes, riboswitches adopt a great diversity of structures. Although many riboswitches share structural principles to build metabolite-competent folds, form precise ligand-binding pockets, and communicate a ligand-binding event to downstream regulatory regions, virtually all riboswitch classes possess unique features for ligand recognition, even those tuned to recognize the same metabolites. Here we present an overview of the biochemical and structural research on riboswitches with a major focus on common principles and individual characteristics adopted by these regulatory RNA elements during evolution to specifically target small molecules and exert genetic responses. This article is part of a Special Issue entitled: Riboswitches.

Computational analysis of riboswitch-based regulation

Advances in computational analysis of riboswitches in the last decade have contributed greatly to our understanding of riboswitch regulatory roles and mechanisms. Riboswitches were originally discovered as part of the sequence analysis of the 5'-untranslated region of mRNAs in the hope of finding novel gene regulatory sites, and the existence of structural RNAs appeared to be a spurious phenomenon. As more riboswitches were discovered, they illustrated the diversity and adaptability of these RNA regulatory sequences. The fact that a chemically monotonous molecule like RNA can discern a wide range of substrates and exert a variety of regulatory mechanisms was subsequently demonstrated in diverse genomes and has hastened the development of sophisticated algorithms for their analysis and prediction. In this review, we focus on some of the computational tools for riboswitch detection and secondary structure prediction. The study of this simple yet efficient form of gene regulation promises to provide a more complete picture of a world that RNA once dominated and allows rational design of artificial riboswitches. This article is part of a Special Issue entitled: Riboswitches.

Riboswitch structure and dynamics by smFRET microscopy

Riboswitches are structured noncoding RNA elements that control the expression of their embedding messenger RNAs by sensing the intracellular concentration of diverse metabolites. As the name suggests, riboswitches are dynamic in nature so that studying their inherent conformational dynamics and ligand-mediated folding is important for understanding their mechanism of action. Single-molecule fluorescence energy transfer (smFRET) microscopy is a powerful and versatile technique for studying the folding pathways and intra- and intermolecular dynamics of biological macromolecules, especially RNA. The ability of smFRET to monitor intramolecular distances and their temporal evolution make it a particularly insightful tool for probing the structure and dynamics of riboswitches. Here, we detail the general steps for using prism-based total internal reflection fluorescence microscopy for smFRET studies of the structure, dynamics, and ligand-binding mechanisms of riboswitches.

Riboswitches in eubacteria sense the second messenger c-di-AMP

Cyclic di-adenosine monophosphate (c-di-AMP) is a recently discovered bacterial second messenger implicated in the control of cell wall metabolism, osmotic stress responses and sporulation. However, the mechanisms by which c-di-AMP triggers these physiological responses have remained largely unknown. Notably, a candidate riboswitch class called ydaO associates with numerous genes involved in these same processes. Although a representative ydaO motif RNA recently was reported to weakly bind ATP, we report that numerous members of this noncoding RNA class selectively respond to c-di-AMP with subnanomolar affinity. Our findings resolve the mystery regarding the primary ligand for this extremely common riboswitch class and expose a major portion of the super-regulon of genes that are controlled by the widespread bacterial second messenger c-di-AMP.

Engineering modular 'ON' RNA switches using biological components

Riboswitches are cis-acting regulatory elements broadly distributed in bacterial mRNAs that control a wide range of critical metabolic activities. Expression is governed by two distinct domains within the mRNA leader: a sensory 'aptamer domain' and a regulatory 'expression platform'. Riboswitches have also received considerable attention as important tools in synthetic biology because of their conceptually simple structure and the ability to obtain aptamers that bind almost any conceivable small molecule using in vitro selection (referred to as SELEX). In the design of artificial riboswitches, a significant hurdle has been to couple the two domains enabling their efficient communication. We previously demonstrated that biological transcriptional 'OFF' expression platforms are easily coupled to diverse aptamers, both biological and SELEX-derived, using simple design rules. Here, we present two modular transcriptional 'ON' riboswitch expression platforms that are also capable of hosting foreign aptamers. We demonstrate that these biological parts can be used to facilely generate artificial chimeric riboswitches capable of robustly regulating transcription both in vitro and in vivo. We expect that these modular expression platforms will be of great utility for various synthetic biological applications that use RNA-based biosensors.

Tuning a riboswitch response through structural extension of a pseudoknot

Structural and dynamic features of RNA folding landscapes represent critical aspects of RNA function in the cell and are particularly central to riboswitch-mediated control of gene expression. Here, using single-molecule fluorescence energy transfer imaging, we explore the folding dynamics of the preQ1 class II riboswitch, an upstream mRNA element that regulates downstream encoded modification enzymes of queuosine biosynthesis. For reasons that are not presently understood, the classical pseudoknot fold of this system harbors an extra stem-loop structure within its 3'-terminal region immediately upstream of the Shine-Dalgarno sequence that contributes to formation of the ligand-bound state. By imaging ligand-dependent preQ1 riboswitch folding from multiple structural perspectives, we reveal that the extra stem-loop strongly influences pseudoknot dynamics in a manner that decreases its propensity to spontaneously fold and increases its responsiveness to ligand binding. We conclude that the extra stem-loop sensitizes this RNA to broaden the dynamic range of the ON/OFF regulatory switch.

Exploiting preQ(1) riboswitches to regulate ribosomal frameshifting

Knowing the molecular details of the interaction between riboswitch aptamers and their corresponding metabolites is important to understand gene expression. Here we report on a novel in vitro assay to study preQ(1) riboswitch aptamers upon binding of 7-aminomethyl-7-deazaguanine (preQ(1)). The assay is based on the ability of the preQ(1) aptamer to fold, upon ligand binding, into a pseudoknotted structure that is capable of stimulating -1 ribosomal frameshifting (-1 FS). Aptamers from three different species were found to induce between 7% and 20% of -1 FS in response to increasing preQ(1) levels, whereas preQ(1) analogues were 100-1000-fold less efficient. In depth mutational analysis of the Fusobacterium nucleatum aptamer recapitulates most of the structural details previously identified for preQ(1) aptamers from other bacteria by crystallography and/or NMR spectroscopy. In addition to providing insight into the role of individual nucleotides of the preQ(1) riboswitch aptamer in ligand binding, the presented system provides a valuable tool to screen small molecules against bacterial riboswitches in a eukaryotic background.