Evidence for a second class of S-adenosylmethionine riboswitches and other regulatory RNA motifs in alpha-proteobacteria

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

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

Multidimensional Applications and Challenges of Riboswitches in Biosensing and Biotherapy

Riboswitches have received significant attention over the last two decades for their multiple functionalities and great potential for applications in various fields. This article highlights and reviews the recent advances in biosensing and biotherapy. These fields involve a wide range of applications, such as food safety detection, environmental monitoring, metabolic engineering, live cell imaging, wearable biosensors, antibacterial drug targets, and gene therapy. The discovery, origin, and optimization of riboswitches are summarized to help readers better understand their multidimensional applications. Finally, this review discusses the multidimensional challenges and development of riboswitches in order to further expand their potential for novel applications.

A spermidine riboswitch class in bacteria exploits a close variant of an aptamer for the enzyme cofactor S-adenosylmethionine

Natural polyamines such as spermidine and spermine cations have characteristics that make them highly likely to be sensed by riboswitches, such as their general affinity to polyanionic RNA and their broad contributions to cell physiology. Despite previous claims that polyamine riboswitches exist, evidence of their biological functions has remained unconvincing. Here, we report that rare variants of bacterial S-adenosylmethionine-I (SAM-I) riboswitches reject SAM and have adapted to selectively sense spermidine. These spermidine-sensing riboswitch variants are associated with genes whose protein products are directly involved in the production of spermidine and other polyamines. Biochemical and genetic assays demonstrate that representatives of this riboswitch class robustly function as genetic "off" switches, wherein spermidine binding causes premature transcription termination to suppress the expression of polyamine biosynthetic genes. These findings confirm the existence of natural spermidine-sensing riboswitches in bacteria and expand the list of variant riboswitch classes that have adapted to bind different ligands.

Characterization and designing of an SAM riboswitch to establish a high-throughput screening platform for SAM overproduction in Saccharomyces cerevisiae

S-adenosyl- l-methionine (SAM) is a high-value compound widely used in the treatment of various diseases. SAM can be produced through fermentation, but further enhancing the microbial production of SAM requires novel high-throughput screening methods for rapid detection and screening of mutant libraries. In this work, an SAM-OFF riboswitch capable of responding to the SAM concentration was obtained and a high-throughput platform for screening SAM overproducers was established. SAM synthase was engineered by semirational design and directed evolution, which resulted in the SAM2S203F,W164R,T251S,Y285F,S365R mutant with almost twice higher catalytic activity than the parental enzyme. The best mutant was then introduced into Saccharomyces cerevisiae BY4741, and the resulting strain BSM8 produced a sevenfold higher SAM titer in shake-flask fermentation, reaching 1.25 g L-1 . This work provides a reference for designing biosensors to dynamically detect metabolite concentrations for high-throughput screening and the construction of effective microbial cell factories.

Linking folding dynamics and function of SAM/SAH riboswitches at the single molecule level

Riboswitches are regulatory elements found in bacterial mRNAs that control downstream gene expression through ligand-induced conformational changes. Here, we used single-molecule FRET to map the conformational landscape of the translational SAM/SAH riboswitch and probe how co-transcriptional ligand-induced conformational changes affect its translation regulation function. Riboswitch folding is highly heterogeneous, suggesting a rugged conformational landscape that allows for sampling of the ligand-bound conformation even in the absence of ligand. The addition of ligand shifts the landscape, favoring the ligand-bound conformation. Mutation studies identified a key structural element, the pseudoknot helix, that is crucial for determining ligand-free conformations and their ligand responsiveness. We also investigated ribosomal binding site accessibility under two scenarios: pre-folding and co-transcriptional folding. The regulatory function of the SAM/SAH riboswitch involves kinetically favoring ligand binding, but co-transcriptional folding reduces this preference with a less compact initial conformation that exposes the Shine-Dalgarno sequence and takes min to redistribute to more compact conformations of the pre-folded riboswitch. Such slow equilibration decreases the effective ligand affinity. Overall, our study provides a deeper understanding of the complex folding process and how the riboswitch adapts its folding pattern in response to ligand, modulates ribosome accessibility and the role of co-transcriptional folding in these processes.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) predicted over 10 inner-shell Mg2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminishes the selectivity between SAM and SAH. One Mg2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg2+ thus plays vital roles in SAM/SAH riboswitch function.

Magnesium ions mediate ligand binding and conformational transition of the SAM/SAH riboswitch

The SAM/SAH riboswitch binds S-adenosylmethionine (SAM) and S-adenosylhomocysteine (SAH) with similar affinities. Mg 2+ is generally known to stabilize RNA structures by neutralizing phosphates, but how it contributes to ligand binding and conformational transition is understudied. Here, extensive molecular dynamics simulations (totaling 120 μs) identified over 10 inner-shell Mg 2+ ions in the SAM/SAH riboswitch. Six of them line the two sides of a groove to widen it and thereby pre-organize the riboswitch for ligand entry. They also form outer-shell coordination with the ligands and stabilize an RNA-ligand hydrogen bond, which effectively diminish the selectivity between SAM and SAH. One Mg 2+ ion unique to the apo form maintains the Shine-Dalgarno sequence in an autonomous mode and thereby facilitates its release for ribosome binding. Mg 2+ thus plays vital roles in SAM/SAH riboswitch function.

Discovering riboswitches: the past and the future

Riboswitches are structured noncoding RNA domains used by many bacteria to monitor the concentrations of target ligands and regulate gene expression accordingly. In the past 20 years over 55 distinct classes of natural riboswitches have been discovered that selectively sense small molecules or elemental ions, and thousands more are predicted to exist. Evidence suggests that some riboswitches might be direct descendants of the RNA-based sensors and switches that were likely present in ancient organisms before the evolutionary emergence of proteins. We provide an overview of the current state of riboswitch research, focusing primarily on the discovery of riboswitches, and speculate on the major challenges facing researchers in the field.

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.

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.

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.

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.

Variants of the guanine riboswitch class exhibit altered ligand specificities for xanthine, guanine, or 2'-deoxyguanosine

The aptamer portions of previously reported riboswitch classes that sense guanine, adenine, or 2′-deoxyguanosine are formed by a highly similar three-stem junction with distinct nucleotide sequences in the regions joining the stems. The nucleotides in these joining regions form the major features of the selective ligand-binding pocket for each aptamer. Previously, we reported the existence of additional, rare variants of the predominant guanine-sensing riboswitch class that carry nucleotide differences in the ligand-binding pocket, suggesting that these RNAs have further diversified their structures and functions. Herein, we report the discovery and analysis of three naturally occurring variants of guanine riboswitches that are narrowly distributed across Firmicutes. These RNAs were identified using comparative sequence analysis methods, which also revealed that some of the gene associations for these variants are atypical for guanine riboswitches or their previously known natural variants. Binding assays demonstrate that the newfound variant riboswitch representatives recognize xanthine, guanine, or 2′-deoxyguanosine, with the guanine class exhibiting greater discrimination against related purines than the more common guanine riboswitch class reported previously. These three additional variant classes, together with the four previously discovered riboswitch classes that employ the same three-stem junction architecture, reveal how a simple structural framework can be diversified to expand the range of purine-based ligands sensed by RNA.

Incoherent dual regulation by a SAM-II riboswitch controlling translation at a distance

In Sinorhizobium meliloti, the methionine biosynthesis genes metA and metZ are preceded by S-adenosyl-L-methionine (SAM) riboswitches of the SAM-II class. Upon SAM binding, structural changes in the metZ riboswitch were predicted to cause transcriptional termination, generating the sRNA RZ. By contrast, the metA riboswitch was predicted to regulate translation from an AUG1 codon. However, downstream of the metA riboswitch, we found a putative Rho-independent terminator and an in-frame AUG2 codon, which may contribute to metA regulation. We validated the terminator between AUG1 and AUG2, which generates the sRNA RA1 that is processed to RA2. Under high SAM conditions, the activities of the metA and metZ promoters and the steady-state levels of the read-through metA and metZ mRNAs were decreased, while the levels of the RZ and RA2 sRNAs were increased. Under these conditions, the sRNAs and the mRNAs were stabilized. Reporter fusion experiments revealed that the Shine–Dalgarno (SD) sequence in the metA riboswitch is required for translation, which, however, starts 74 nucleotides downstream at AUG2, suggesting a novel translation initiation mechanism. Further, the reporter fusion data supported the following model of RNA-based regulation: Upon SAM binding by the riboswitch, the SD sequence is sequestered to downregulate metA translation, while the mRNA is stabilized. Thus, the SAM-II riboswitches fulfil incoherent, dual regulation, which probably serves to ensure basal metA and metZ mRNA levels under high SAM conditions. This probably helps to adapt to changing conditions and maintain SAM homoeostasis.

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

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

Riboswitch Mechanisms: New Tricks for an Old Dog

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

The potential versatility of RNA catalysis

It is commonly thought that in the early development of life on this planet RNA would have acted both as a store of genetic information and as a catalyst. While a number of RNA enzymes are known in contemporary cells, they are largely confined to phosphoryl transfer reactions, whereas an RNA based metabolism would have required a much greater chemical diversity of catalysis. Here we discuss how RNA might catalyze a wider variety of chemistries, and particularly how information gleaned from riboswitches could suggest how ribozymes might recruit coenzymes to expand their chemical range. We ask how we might seek such activities in modern biology. This article is categorized under: RNA-Based Catalysis > Miscellaneous RNA-Catalyzed Reactions Regulatory RNAs/RNAi/Riboswitches > Riboswitches RNA Structure and Dynamics > RNA Structure, Dynamics and Chemistry.

Methionine biosynthesis pathway genes affect curdlan biosynthesis of Agrobacterium sp. CGMCC 11546 via energy regeneration

Curdlan is a water-insoluble exopolysaccharide produced by Agrobacterium species under nitrogen starvation. The curdlan production in the ΔmdeA, ΔmetA, ΔmetH, and ΔmetZ mutants of methionine biosynthesis pathway of Agrobacterium sp. CGMCC 11546 were significantly impaired. Fermentation profiles of four mutants showed that the consumption of ammonia and sucrose was impaired. Transcriptome analysis of the ΔmetH and ΔmetZ mutants showed that numerous differentially expressed genes involved in the electron transfer chain (ETC) were significantly down-regulated, suggesting that methionine biosynthesis pathway affected the production of energy ATP during the curdlan biosynthesis. Furthermore, metabolomics analysis of the ΔmetH and ΔmetZ mutants showed that ADP and FAD were significantly accumulated, while acetyl-CoA was diminished, suggesting that the impaired curdlan production in the ΔmetH and ΔmetZ mutants might be caused by the insufficient supply of energy ATP. Finally, the addition of both dibasic sodium succinate as a substrate of FAD recycling and methionine significantly restored the curdlan production of four mutants. In conclusion, methionine biosynthesis pathway plays an important role in curdlan biosynthesis in Agrobacterium sp. CGMCC 11546, which affected the sufficient supply of energy ATP from the ETC during the curdlan biosynthesis.

Comprehensive discovery of novel structured noncoding RNAs in 26 bacterial genomes

Comparative sequence analysis methods are highly effective for uncovering novel classes of structured noncoding RNAs (ncRNAs) from bacterial genomic DNA sequence datasets. Previously, we developed a computational pipeline to more comprehensively identify structured ncRNA representatives from individual bacterial genomes. This search process exploits the fact that genomic regions serving as templates for the transcription of structured RNAs tend to be present in longer than average noncoding 'intergenic regions' (IGRs) that are enriched in G and C nucleotides compared to the remainder of the genome. In the present study, we apply this computational pipeline to identify structured ncRNA candidates from 26 diverse bacterial species. Numerous novel structured ncRNA motifs were discovered, including several riboswitch candidates, one whose ligand has been identified and others that have yet to be experimentally validated. Our findings support recent predictions that hundreds of novel ribo-switch classes and other ncRNAs remain undiscovered among the limited number of bacterial species whose genomes have been completely sequenced.

Structural diversity and phylogenetic distribution of valyl tRNA-like structures in viruses

Viruses commonly use specifically folded RNA elements that interact with both host and viral proteins to perform functions important for diverse viral processes. Examples are found at the 3' termini of certain positive-sense ssRNA virus genomes where they partially mimic tRNAs, including being aminoacylated by host cell enzymes. Valine-accepting tRNA-like structures (TLS^Val) are an example that share some clear homology with canonical tRNAs but have several important structural differences. Although many examples of TLS^Val have been identified, we lacked a full understanding of their structural diversity and phylogenetic distribution. To address this, we undertook an in-depth bioinformatic and biochemical investigation of these RNAs, guided by recent high-resolution structures of a TLS^Val We cataloged many new examples in plant-infecting viruses but also in unrelated insect-specific viruses. Using biochemical and structural approaches, we verified the secondary structure of representative TLS^Val substrates and tested their ability to be valylated, confirming previous observations of structural heterogeneity within this class. In a few cases, large stem-loop structures are inserted within variable regions located in an area of the TLS distal to known host cell factor binding sites. In addition, we identified one virus whose TLS has switched its anticodon away from valine, causing a loss of valylation activity; the implications of this remain unclear. These results refine our understanding of the structural and functional mechanistic details of tRNA mimicry and how this may be used in viral infection.

Novel small RNAs expressed by Bartonella bacilliformis under multiple conditions reveal potential mechanisms for persistence in the sand fly vector and human host

Bartonella bacilliformis, the etiological agent of Carrión’s disease, is a Gram-negative, facultative intracellular alphaproteobacterium. Carrión’s disease is an emerging but neglected tropical illness endemic to Peru, Colombia, and Ecuador. B. bacilliformis is spread between humans through the bite of female phlebotomine sand flies. As a result, the pathogen encounters significant and repeated environmental shifts during its life cycle, including changes in pH and temperature. In most bacteria, small non-coding RNAs (sRNAs) serve as effectors that may post-transcriptionally regulate the stress response to such changes. However, sRNAs have not been characterized in B. bacilliformis, to date. We therefore performed total RNA-sequencing analyses on B. bacilliformis grown in vitro then shifted to one of ten distinct conditions that simulate various environments encountered by the pathogen during its life cycle. From this, we identified 160 sRNAs significantly expressed under at least one of the conditions tested. sRNAs included the highly-conserved tmRNA, 6S RNA, RNase P RNA component, SRP RNA component, ffH leader RNA, and the alphaproteobacterial sRNAs αr45 and speF leader RNA. In addition, 153 other potential sRNAs of unknown function were discovered. Northern blot analysis was used to confirm the expression of eight novel sRNAs. We also characterized a Bartonellabacilliformisgroup I intron (BbgpI) that disrupts an un-annotated tRNA_CCU^Arg gene and determined that the intron splices in vivo and self-splices in vitro. Furthermore, we demonstrated the molecular targeting of Bartonellabacilliformissmall RNA 9 (BbsR9) to transcripts of the ftsH, nuoF, and gcvT genes, in vitro.

B. bacilliformis is a bacterial pathogen that is transmitted between humans by phlebotomine sand flies. Bacteria often express sRNAs to fine-tune the production of proteins involved in a wide array of biological processes. We cultured B. bacilliformis in vitro under standard conditions then shifted the pathogen for a period of time to ten distinct environments, including multiple temperatures, pH levels, and infections of human blood and human vascular endothelial cells. After RNA-sequencing, a manual transcriptome search identified 160 putative sRNAs, including seven highly-conserved sRNAs and 153 novel potential sRNAs. We then characterized two of the novel sRNAs, BbgpI and BbsR9. BbgpI is a group I intron (ribozyme) that self-splices and disrupts an unannotated gene coding for a transfer RNA (tRNA_CCU^Arg). BbsR9 is an intergenic sRNA expressed under conditions that simulate the sand fly. We found that BbsR9 targets transcripts of the ftsH, nuoF, and gcvT genes. Furthermore, we determined the specific sRNA-mRNA interactions responsible for BbsR9 binding to its target mRNAs through in vitro mutagenesis and binding assays.

Widespread Protection of RNA Cleavage Sites by a Riboswitch Aptamer that Folds as a Compact Obstacle to Scanning by RNase E

Riboswitches are thought generally to function by modulating transcription elongation or translation initiation. In rare instances, ligand binding to a riboswitch has been found to alter the rate of RNA degradation by directly stimulating or inhibiting nearby cleavage. Here, we show that guanidine-induced pseudoknot formation by the aptamer domain of a guanidine III riboswitch from Legionella pneumophila has a different effect, stabilizing mRNA by protecting distal cleavage sites en masse from ribonuclease attack. It does so by creating a coaxially base-paired obstacle that impedes scanning from a monophosphorylated 5' end to those sites by the regulatory endonuclease RNase E. Ligand binding by other riboswitch aptamers peripheral to the path traveled by RNase E does not inhibit distal cleavage. These findings reveal that a riboswitch aptamer can function independently of any overlapping expression platform to regulate gene expression by acting directly to prolong mRNA longevity in response to ligand binding.

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.

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.

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

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

Interactions between SAM and the 5' UTR mRNA of the sam1 gene regulate translation in S. pombe

The 5' untranslated region (5' UTR) of eukaryotic mRNA plays an important role in translation. Here we report the function of the 5' UTR mRNA of S-adenosylmethionine synthetase (sam1) in translational modulation in the presence of SAM in fission yeast Schizosaccharomyces pombe Reporter assays, binding and chemical probing experiments, and mutational analysis show that the 5' UTR mRNA of sam1 binds to SAM to effect translation. Translational modulation is dependent on a tertiary structure transition in the RNA upon SAM binding. The characterization of such an RNA that is directly associated with an essential metabolic process in eukaryotes provides additional evidence that ligand binding by RNAs plays an important role in eukaryotic gene regulation.

In vivo evolutionary engineering of riboswitch with high-threshold for N-acetylneuraminic acid production

Riboswitches with desired properties, such as sensitivity, threshold, dynamic range, is important for its application. However, the property change of a natural riboswitch is difficult due to the lack of the understanding of aptamer ligand binding properties and a proper screening method for both rational and irrational design. In this study, an effective method to change the threshold of riboswitch was established in vivo based on growth coupled screening by combining both positive and negative selections. The feasibility of the method was verified by the model library. Using this method, an N-acetylneuraminic acid (NeuAc) riboswitch was evolved and modified riboswitches with high threshold and large dynamic range were obtained. Then, using a new NeuAc riboswitch, both ribosome binding sites and key gene in NeuAc biosynthesis pathway were optimized. The highest NeuAc production of 14.32 g/l that has been reported using glucose as sole carbon source was obtained.

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.

The structure of the SAM/SAH-binding riboswitch

S-adenosylmethionine (SAM) is a central metabolite since it is used as a methyl group donor in many different biochemical reactions. Many bacteria control intracellular SAM concentrations using riboswitch-based mechanisms. A number of structurally different riboswitch families specifically bind to SAM and mainly regulate the transcription or the translation of SAM-biosynthetic enzymes. In addition, a highly specific riboswitch class recognizes S-adenosylhomocysteine (SAH)—the product of SAM-dependent methyl group transfer reactions—and regulates enzymes responsible for SAH hydrolysis. High-resolution structures are available for many of these riboswitch classes and illustrate how they discriminate between the two structurally similar ligands SAM and SAH. The so-called SAM/SAH riboswitch class binds both ligands with similar affinities and is structurally not yet characterized. Here, we present a high-resolution nuclear magnetic resonance structure of a member of the SAM/SAH-riboswitch class in complex with SAH. Ligand binding induces pseudoknot formation and sequestration of the ribosome binding site. Thus, the SAM/SAH-riboswitches are translational ‘OFF’-switches. Our results establish a structural basis for the unusual bispecificity of this riboswitch class. In conjunction with genomic data our structure suggests that the SAM/SAH-riboswitches might be an evolutionary late invention and not a remnant of a primordial RNA-world as suggested for other riboswitches.

Disruption of the OLE ribonucleoprotein complex causes magnesium toxicity in Bacillus halodurans

OLE RNAs represent an unusual class of bacterial noncoding RNAs common in Gram-positive anaerobes. The OLE RNA of the alkaliphile Bacillus halodurans is highly expressed and naturally interacts with at least two RNA-binding proteins called OapA and OapB. The phenotypes of the corresponding knockouts include growth inhibition when exposed to ethanol or other short-chain alcohols or when incubated at modestly reduced temperatures (e.g. 20°C). Intriguingly, the OapA 'PM1' mutant, which carries two amino acid changes to a highly conserved region, yields a dominant-negative phenotype that causes more severe growth defects under these same stress conditions. Herein, we report that the PM1 strain also exhibits extreme sensitivity to elevated Mg²⁺ concentrations, beginning as low as 2 mM. Suppressor mutants predominantly map to genes for aconitate hydratase and isocitrate dehydrogenase, which are expected to alter cellular citrate concentrations. Citrate reduces the severity of the Mg²⁺ toxicity phenotype, but neither the genomic mutations nor the addition of citrate to the medium overcomes ethanol toxicity or temperature sensitivity. These findings reveal that OLE RNA and its protein partners are involved in biochemical responses under several stress conditions, wherein the unusual sensitivity to Mg²⁺ can be independently suppressed by specific genomic mutations.

Structural prediction of RNA switches using conditional base-pair probabilities

An RNA switch triggers biological functions by toggling between two conformations. RNA switches include bacterial riboswitches, where ligand binding can stabilize a bound structure. For RNAs with only one stable structure, structural prediction usually just requires a straightforward free energy minimization, but for an RNA switch, the prediction of a less stable alternative structure is often computationally costly and even problematic. The current sampling-clustering method predicts stable and alternative structures by partitioning structures sampled from the energy landscape into two clusters, but it is very time-consuming. Instead, we predict the alternative structure of an RNA switch from conditional probability calculations within the energy landscape. First, our method excludes base pairs related to the most stable structure in the energy landscape. Then, it detects stable stems (“seeds”) in the remaining landscape. Finally, it folds an alternative structure prediction around a seed. While having comparable riboswitch classification performance, the conditional-probability computations had fewer adjustable parameters, offered greater predictive flexibility, and were more than one thousand times faster than the sampling step alone in sampling-clustering predictions, the competing standard. Overall, the described approach helps traverse thermodynamically improbable energy landscapes to find biologically significant substructures and structures rapidly and effectively.

Polyamine biosynthesis and biological roles in rhizobia

Polyamines are ubiquitous molecules containing two or more amino groups that fulfill varied and often essential physiological and regulatory roles in all organisms. In the symbiotic nitrogen-fixing bacteria known as rhizobia, putrescine and homospermidine are invariably produced while spermidine and norspermidine synthesis appears to be restricted to the alfalfa microsymbiont Sinorhizobium meliloti. Studies with rhizobial mutants deficient in the synthesis of one or more polyamines have shown that these compounds are important for growth, stress resistance, motility, exopolysaccharide production and biofilm formation. In this review, we describe these studies and examine how polyamines are synthesized and regulated in rhizobia.

mRNA dynamics and alternative conformations adopted under low and high arginine concentrations control polyamine biosynthesis in Salmonella

Putrescine belongs to the large group of polyamines, an essential class of metabolites that exists throughout all kingdoms of life. The Salmonella speF gene encodes an inducible ornithine decarboxylase that produces putrescine from ornithine. Putrescine can be also synthesized from arginine in a parallel metabolic pathway. Here, we show that speF expression is controlled at multiple levels through regulatory elements contained in a long leader sequence. At the heart of this regulation is a short open reading frame, orf34, which is required for speF production. Translation of orf34 interferes with Rho-dependent transcription termination and helps to unfold an inhibitory RNA structure sequestering speF ribosome-binding site. Two consecutive arginine codons in the conserved domain of orf34 provide a third level of speF regulation. Uninterrupted translation of orf34 under conditions of high arginine allows the formation of a speF mRNA structure that is degraded by RNase G, whereas ribosome pausing at the consecutive arginine codons in the absence of arginine enables the formation of an alternative structure that is resistant to RNase G. Thus, the rate of ribosome progression during translation of the upstream ORF influences the dynamics of speF mRNA folding and putrescine production. The identification of orf34 and its regulatory functions provides evidence for the evolutionary conservation of ornithine decarboxylase regulatory elements and putrescine production.

Polyamines are widely distributed in nature, they bind nucleic acids and proteins and although their exact mechanism of action is not clear, their effect on fundamental cellular functions is well documented. The canonical biosynthesis pathway of polyamines is conserved and begins with speF encoding ornithine decarboxylase, an inducible enzyme that produces putrescine from ornithine. Putrescine can also be produced from arginine in an alternative metabolic pathway. Here, we show that the rate of ribosome progression during translation of a short ORF (ORF34) upstream of speF influences the dynamics of speF mRNA folding and thus putrescine production. Uninterrupted translation of orf34 carrying two consecutive arginine codons, under conditions of high arginine, results in the formation of a speF mRNA structure that is degraded by RNase G, whereas ribosomes slow-down at the consecutive arginine codons in the absence of arginine enables the formation of an alternative structure that is unsusceptible to RNase G and thus results in putrescine production. The study of Salmonella speF regulation provides evidence that, despite variations in the mechanistic details, RNA-based regulation of putrescine biosynthesis and ornithine decarboxylase is conserved from bacteria to mammals.

Large Noncoding RNAs in Bacteria

Bacterial noncoding RNA (ncRNA) classes longer than 200 nucleotides are rare but are responsible for performing some of the most fundamental tasks in living cells. RNAs such as 16S and 23S rRNA, group I and group II introns, RNase P ribozymes, transfer-messenger RNAs, and coenzyme B₁₂ riboswitches are diverse in structure and accomplish biochemical functions that rival the activities of proteins. Over the last decade, a number of new classes of large ncRNAs have been uncovered in bacteria. A total of 21 classes with no established functions have been identified through the use of bioinformatics search strategies. Based on precedents for bacterial large ncRNAs performing sophisticated functions, it seems likely that some of these structured ncRNAs also will prove to carry out complex functions. Thus, determining their roles will provide a better understanding of fundamental biological processes. A few studies have produced data that provide clues to the purposes of some of these recently found classes, but the true functions of most classes remain mysterious.

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.

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 dual regulatory circuit consisting of S-adenosylmethionine decarboxylase protein and its reaction product controls expression of the paralogous activator prozyme in Trypanosoma brucei

Polyamines are essential for cell growth of eukaryotes including the etiologic agent of human African trypanosomiasis (HAT), Trypanosoma brucei. In trypanosomatids, a key enzyme in the polyamine biosynthetic pathway, S-adenosylmethionine decarboxylase (TbAdoMetDC) heterodimerizes with a unique catalytically-dead paralog called prozyme to form the active enzyme complex. In higher eukaryotes, polyamine metabolism is subject to tight feedback regulation by spermidine-dependent mechanisms that are absent in trypanosomatids. Instead, in T. brucei an alternative regulatory strategy based on TbAdoMetDC prozyme has evolved. We previously demonstrated that prozyme protein levels increase in response to loss of TbAdoMetDC activity. Herein, we show that prozyme levels are under translational control by monitoring incorporation of deuterated leucine into nascent prozyme protein. We furthermore identify pathway factors that regulate prozyme mRNA translation. We find evidence for a regulatory feedback mechanism in which TbAdoMetDC protein and decarboxylated AdoMet (dcAdoMet) act as suppressors of prozyme translation. In TbAdoMetDC null cells expressing the human AdoMetDC enzyme, prozyme levels are constitutively upregulated. Wild-type prozyme levels are restored by complementation with either TbAdoMetDC or an active site mutant, suggesting that TbAdoMetDC possesses an enzyme activity-independent function that inhibits prozyme translation. Depletion of dcAdoMet pools by three independent strategies: inhibition/knockdown of TbAdoMetDC, knockdown of AdoMet synthase, or methionine starvation, each cause prozyme upregulation, providing independent evidence that dcAdoMet functions as a metabolic signal for regulation of the polyamine pathway in T. brucei. These findings highlight a potential regulatory paradigm employing enzymes and pseudoenzymes that may have broad implications in biology.

Trypanosoma brucei is a single-celled eukaryotic pathogen and the causative agent of human African trypanosomiasis (HAT). Polyamines are organic polycations that are essential for growth in T. brucei to facilitate protein translation and to maintain redox homeostasis. The pathway is the target of eflornithine, a current frontline therapy for treatment of HAT. Polyamine biosynthetic enzymes are regulated at multiple levels in mammals (e.g. transcription, translation and protein turnover), but in contrast, T. brucei lacks these mechanisms. Instead in T. brucei a central enzyme in polyamine metabolism called AdoMetDC must form a complex with a sister protein (termed a pseudoenzyme) to be active. Herein, we show that cellular levels of this sister protein we call prozyme are in turn feedback regulated by both AdoMetDC and by its reaction product in response to cell treatments that reduce pathway output. This regulatory paradigm highlights how pseudoenzymes can evolve to play an important role in metabolic pathway regulation and in organismal fitness.

SuhB, a small non-coding RNA involved in catabolite repression of tetralin degradation genes in Sphingopyxis granuli strain TFA

Global dRNA-seq analysis of transcription start sites combined with in silico annotation using Infernal software revealed the expression of 91 putative non-coding sRNA in Sphingopyxis granuli TFA cells grown on different carbon sources. Excluding housekeeping sRNAs, only one additional sRNA, which belongs to the Rfam SuhB family (RF00519), was detected by Infernal but with an incorrect size according to the experimental results. SuhB is highly conserved across the Sphingopyxis genus. Expression data revealed that SuhB is present in rapidly growing TFA cells. A suhB deletion mutant exhibited de-repression of tetralin degradation (thn) gene expression and higher amounts of their LysR-type activator, ThnR, under conditions of carbon catabolite repression (CCR). Interaction between SuhB and the 5'UTR of thnR mRNA was demonstrated in vitro. Moreover, co-immunoprecipitation experiments, combined with fluorescence measurements of gfp fusions to the 5'UTR of thnR mRNA and the phenotype of an hfq deletion mutant, suggest the involvement of Hfq in this interaction. Taken together, these data support an Hfq-mediated repressive role for SuhB, on ThnR mRNA translation that prevents thn gene induction. SuhB, which is a highly conserved sRNA in the Sphingopyxis genus, is the first identified element directly involved in CCR of thn gene expression in S. granuli strain TFA.

Classification of riboswitch sequences using k-mer frequencies

Riboswitches are non-coding RNAs that regulate gene expression by altering the structural conformation of mRNA transcripts. Their regulation mechanism might be exploited for interesting biomedical applications such as drug targets and biosensors. A major challenge consists in accurately identifying metabolite-binding RNA switches which are structurally complex and diverse. In this regard, we investigated the classification of 16 riboswitch families using supervised learning algorithms trained solely with sequence-based features. We generated a reduced feature set and proposed a visual representation to explore its components. We induced Support Vector Machine, Random Forest, Naive Bayes, J48, and HyperPipes classifiers with our proposed feature set and tested their performance over independent data. Our best multi-class classifier achieved F-measure values of 0.996 and 0.966 in the training and test phases, respectively, outperforming those of a previous approach. When compared against BLAST, our best classifiers yielded competitive results. This work shows that the classifiers trained with our sequence-based feature set accurately discriminate riboswitches.

Structure and ligand binding of the SAM-V riboswitch

SAM-V is one of the class of riboswitches that bind S-adenosylmethione, regulating gene expression by controlling translation. We have solved the crystal structure of the metY SAM-V riboswitch bound to its SAM ligand at 2.5 Å resolution. The RNA folds as an H-type pseudoknot, with a major-groove triple helix in which resides the SAM ligand binding site. The bound SAM adopts an elongated conformation aligned with the axis of the triple helix, and is held at either end by hydrogen bonding to the adenine and the amino acid moieties. The central sulfonium cation makes electrostatic interactions with an U:A.U base triple, so conferring specificity. We propose a model in which SAM binding leads to association of the triplex third strand that stabilizes a short helix and occludes the ribosome binding site. Thus the new structure explains both ligand specificity and the mechanism of genetic control.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

Challenges of ligand identification for the second wave of orphan riboswitch candidates

Orphan riboswitch candidates are noncoding RNA motifs whose representatives are believed to function as genetic regulatory elements, but whose target ligands have yet to be identified. The study of certain orphans, particularly classes that have resisted experimental validation for many years, has led to the discovery of important biological pathways and processes once their ligands were identified. Previously, we highlighted details for four of the most common and intriguing orphan riboswitch candidates. This facilitated the validation of riboswitches for the signaling molecules c-di-AMP, ZTP, and ppGpp, the metal ion Mn2+, and the metabolites guanidine and PRPP. Such studies also yield useful linkages between the ligands sensed by the riboswitches and numerous biochemical pathways. In the current report, we describe the known characteristics of 30 distinct classes of orphan riboswitch candidates - some of which have remained unsolved for over a decade. We also discuss the prospects for uncovering novel biological insights via focused studies on these RNAs. Lastly, we make recommendations for experimental objectives along the path to finding ligands for these mysterious RNAs.

SAM-VI RNAs selectively bind S-adenosylmethionine and exhibit similarities to SAM-III riboswitches

Five distinct riboswitch classes that regulate gene expression in response to the cofactor S-adenosylmethionine (SAM) or its metabolic breakdown product S-adenosylhomocysteine (SAH) have been reported previously. Collectively, these SAM- or SAH-sensing RNAs constitute the most abundant collection of riboswitches, and are found in nearly every major bacterial lineage. Here, we report a potential sixth member of this pervasive riboswitch family, called SAM-VI, which is predominantly found in Bifidobacterium species. SAM-VI aptamers selectively bind the cofactor SAM and strongly discriminate against SAH. The consensus sequence and structural model for SAM-VI share some features with the consensus model for the SAM-III riboswitch class, whose members are mainly found in lactic acid bacteria. However, there are sufficient differences between the two classes such that current bioinformatics methods separately cluster representatives of the two motifs. These findings highlight the abundance of RNA structures that can form to selectively recognize SAM, and showcase the ability of RNA to utilize diverse strategies to perform similar biological functions.

RNAseq analysis of α-proteobacterium Gluconobacter oxydans 621H

Background

The acetic acid bacterium Gluconobacter oxydans 621H is characterized by its exceptional ability to incompletely oxidize a great variety of carbohydrates in the periplasm. The metabolism of this α-proteobacterium has been characterized to some extent, yet little is known about its transcriptomes and related data. In this study, we applied two different RNAseq approaches. Primary transcriptomes enriched for 5′-ends of transcripts were sequenced to detect transcription start sites, which allow subsequent analysis of promoter motifs, ribosome binding sites, and 5´-UTRs. Whole transcriptomes were sequenced to identify expressed genes and operon structures.

Results

Sequencing of primary transcriptomes of G. oxydans revealed 2449 TSSs, which were classified according to their genomic context followed by identification of promoter and ribosome binding site motifs, analysis of 5´-UTRs including validation of predicted cis-regulatory elements and correction of start codons. 1144 (41%) of all genes were found to be expressed monocistronically, whereas 1634 genes were organized in 571 operons. Together, TSSs and whole transcriptome data were also used to identify novel intergenic (18), intragenic (328), and antisense transcripts (313).

Conclusions

This study provides deep insights into the transcriptional landscapes of G. oxydans. The comprehensive transcriptome data, which we made publicly available, facilitate further analysis of promoters and other regulatory elements. This will support future approaches for rational strain development and targeted gene expression in G. oxydans. The corrections of start codons further improve the high quality genome reference and support future proteome analysis.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4415-x) contains supplementary material, which is available to authorized users.

An account of evolutionary specialization: the AbcR small RNAs in the Rhizobiales

The AbcR small RNAs (sRNAs) are a fascinating example of two highly conserved sRNAs that differ tremendously at the functional level among organisms. From their transcriptional activation to their regulatory capabilities, the AbcR sRNAs exhibit varying characteristics in three well-studied bacteria belonging to the Rhizobiales order: the plant symbiont Sinorhizobium meliloti, the plant pathogen Agrobacterium tumefaciens, and the animal pathogen Brucella abortus. This review outlines the similarities and differences of the AbcR sRNAs between each of these organisms, and discusses reasons as to why this group of sRNAs has diverged in their genetic organization and regulatory functions across species. In the end, this review will shed light on how regulatory systems, although seemingly conserved among bacteria, can vary based on the environmental niche and lifestyle of an organism.

Structure of the Guanidine III Riboswitch

Riboswitches are structural elements found in mRNA molecules that couple small-molecule binding to regulation of gene expression, usually by controlling transcription or translation. We have determined high-resolution crystal structures of the ykkC guanidine III riboswitch from Thermobifida fusca. The riboswitch forms a classic H-type pseudoknot that includes a triple helix that is continuous with a central core of conserved nucleotides. These form a left-handed helical ramp of inter-nucleotide interactions, generating the guanidinium cation binding site. The ligand is hydrogen bonded to the Hoogsteen edges of two guanine bases. The binding pocket has a side opening that can accommodate a small side chain, shown by structures with bound methylguanidine, aminoguanidine, ethylguanidine, and agmatine. Comparison of the new structure with those of the guanidine I and II riboswitches reveals that evolution generated three different structural solutions for guanidine binding and subsequent gene regulation, although with some common elements.

RNA-mediated signal perception in pathogenic bacteria

Bacterial pathogens encounter several different environments during an infection, many of them possibly being detrimental. In order to sense its surroundings and adjust the gene expression accordingly, different regulatory schemes are undertaken. With these, the bacterium appropriately can differentiate between various environmental cues to express the correct virulence factor at the appropriate time and place. An attractive regulator device is RNA, which has an outstanding ability to alter its structure in response to external stimuli, such as metabolite concentration or alterations in temperature, to control its downstream gene expression. This review will describe the function of riboswitches and thermometers, with a particular emphasis on regulatory RNAs being important for bacterial pathogenicity. WIREs RNA 2017, 8:e1429. doi: 10.1002/wrna.1429 For further resources related to this article, please visit the WIREs website.

Riboswitch-Mediated Gene Regulation: Novel RNA Architectures Dictate Gene Expression Responses

Riboswitches are RNA elements that act on the mRNA with which they are cotranscribed to modulate expression of that mRNA. These elements are widely found in bacteria, where they have a broad impact on gene expression. The defining feature of riboswitches is that they directly recognize a physiological signal, and the resulting shift in RNA structure affects gene regulation. The majority of riboswitches respond to cellular metabolites, often in a feedback loop to repress synthesis of the enzymes used to produce the metabolite. Related elements respond to the aminoacylation status of a specific tRNA or to a physical parameter, such as temperature or pH. Recent studies have identified new classes of riboswitches and have revealed new insights into the molecular mechanisms of signal recognition and gene regulation. Application of structural and biophysical approaches has complemented previous genetic and biochemical studies, yielding new information about how different riboswitches operate.

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.