A widespread riboswitch candidate that controls bacterial genes involved in molybdenum cofactor and tungsten cofactor metabolism

Shared functions of Fe-S cluster assembly and Moco biosynthesis

Molybdenum cofactor (Moco) biosynthesis is a complex process that involves the coordinated function of several proteins. In the recent years it has become evident that the availability of Fe-S clusters play an important role for the biosynthesis of Moco. First, the MoaA protein binds two [4Fe-4S] clusters per monomer. Second, the expression of the moaABCDE and moeAB operons is regulated by FNR, which senses the availability of oxygen via a functional [4Fe-4S] cluster. Finally, the conversion of cyclic pyranopterin monophosphate to molybdopterin requires the availability of the L-cysteine desulfurase IscS, which is an enzyme involved in the transfer of sulfur to various acceptor proteins with a main role in the assembly of Fe-S clusters. In this review, we dissect the dependence of the production of active molybdoenzymes in detail, starting from the regulation of gene expression and further explaining sulfur delivery and Fe-S cluster insertion into target enzymes. Further, Fe-S cluster assembly is also linked to iron availability. While the abundance of selected molybdoenzymes is largely decreased under iron-limiting conditions, we explain that the expression of the genes is dependent on an active FNR protein. FNR is a very important transcription factor that represents the master-switch for the expression of target genes in response to anaerobiosis. Moco biosynthesis is further directly dependent on the presence of ArcA and also on an active Fur protein.

Iron limitation indirectly reduces the Escherichia coli torCAD operon expression by a reduction of molybdenum cofactor availability

The expression of most molybdoenzymes in Escherichia coli has so far been revealed to be regulated by anaerobiosis and requires the presence of iron, based on the necessity of the transcription factor FNR to bind one [4Fe-4S] cluster. One exception is trimethylamine-N-oxide reductase encoded by the torCAD operon, which has been described to be expressed independently from FNR. In contrast to other alternative anaerobic respiratory systems, the expression of the torCAD operon was shown not to be completely repressed by the presence of dioxygen. To date, the basis for the O2-dependent expression of the torCAD operon has been related to the abundance of the transcriptional regulator IscR, which represses the transcription of torS and torT, and is more abundant under aerobic conditions than under anaerobic conditions. In this study, we reinvestigated the regulation of the torCAD operon and its dependence on the presence of iron and identified a novel regulation that depends on the presence of the bis-molybdopterin guanine dinucleotide (bis-MGD) molybdenum cofactor . We confirmed that the torCAD operon is directly regulated by the heme-containing protein TorC and is indirectly regulated by ArcA and by the availability of iron via active FNR and Fur, both regulatory proteins that influence the synthesis of the molybdenum cofactor. Furthermore, we identified a novel regulation mode of torCAD expression that is dependent on cellular levels of bis-MGD and is not used by other bis-MGD-containing enzymes like nitrate reductase.IMPORTANCEIn bacteria, molybdoenzymes are crucial for anaerobic respiration using alternative electron acceptors. FNR is a very important transcription factor that represents the master switch for the expression of target genes in response to anaerobiosis. Only Escherichia coli trimethylamine-N-oxide (TMAO) reductase escapes this regulation by FNR. We identified that the expression of TMAO reductase is regulated by the amount of bis-molybdopterin guanine dinucleotide (bis-MGD) cofactor synthesized by the cell itself, representing a novel regulation pathway for the expression of an operon coding for a molybdoenzyme. Furthermore, TMAO reductase gene expression is indirectly regulated by the presence of iron, which is required for the production of the bis-MGD cofactor in the cell.

Development of Cupriavidus necator H16 as a host for heterologous production of formate dehydrogenase I of Methylorubrum extorquens: Possibilities and limitations

The discovery of formate dehydrogenase (Me-FDH1) from Methylorubrum extorquens has provided an avenue for sustainable CO2 fixation and utilization. However, the mass production of Me-FDH1 is challenging due to the presence of its unique tungsto-bis-metalopterin guanine dinucleotide (W-bis-MGD) cofactor, limiting its practical applications. In this study, C. necator H16 is proposed as a host for the large-scale production of Me-FDH1, utilizing fructose as a carbon source and its inherent machinery for cofactor synthesis. In a minimal salt medium, C. necator H16 could produce active Me-FDH1, which exhibited a specific activity of 80 to 100 U/mg for CO2 conversion to formate. In fed batch bioreactor experiments, approximately 50 g CDW/L (cell dry weight/L) and 10,000 U/L Me-FDH1 were achieved within 50 h. This study highlights C. necator H16 as the recombinant host for Me-FDH1, paving the way for the future development of efficient mass-production methods for this crucial enzyme.

G·U base pairing motifs in long non-coding RNAs

Long non-coding RNAs (lncRNAs) are recently-discovered transcripts involved in gene expression regulation and associated with diseases. Despite the unprecedented molecular complexity of these transcripts, recent studies of the secondary and tertiary structure of lncRNAs are starting to reveal the principles of lncRNA structural organization, with important functional implications. It therefore starts to be possible to analyze lncRNA structures systematically. Here, using a set of prototypical and medically-relevant lncRNAs of known secondary structure, we specifically catalogue the distribution and structural environment of one of the first-identified and most frequently occurring non-canonical Watson-Crick interactions, the G·U base pair. We compare the properties of G·U base pairs in our set of lncRNAs to those of the G·U base pairs in other well-characterized transcripts, like rRNAs, tRNAs, ribozymes, and riboswitches. Furthermore, we discuss how G·U base pairs in these targets participate in establishing interactions with proteins or miRNAs, and how they enable lncRNA tertiary folding by forming intramolecular or metal-ion interactions. Finally, by identifying highly-G·U-enriched regions of yet unknown function in our target lncRNAs, we provide a new rationale for future experimental investigation of these motifs, which will help obtain a more comprehensive understanding of lncRNA functions and molecular mechanisms in the future.

Riboswitches, from cognition to transformation

Riboswitches are functional RNA elements that regulate gene expression by directly detecting metabolites. Twenty years have passed since it was first discovered, researches on riboswitches are becoming increasingly standardized and refined, which could significantly promote people's cognition of RNA function as well. Here, we focus on some representative orphan riboswitches, enumerate the structural and functional transformation and artificial design of riboswitches including the coupling with ribozymes, hoping to attain a comprehensive understanding of riboswitch research.

The structural features of the ligand-free moaA riboswitch and its ion-dependent folding

Riboswitches are structural elements of mRNA involved in the regulation of gene expression by responding to specific cellular metabolites. To fulfil their regulatory function, riboswitches prefold into an active state, the so-called binding competent form, that guarantees metabolite binding and allows a consecutive refolding of the RNA. Here, we describe the folding pathway to the binding competent form as well as the ligand free structure of the moaA riboswitch of E. coli. This RNA proposedly responds to the molybdenum cofactor (Moco), a highly oxygen-sensitive metabolite, essential in the carbon and sulfur cycles of eukaryotes. K⁺- and Mg²⁺-dependent footprinting assays and spectroscopic investigations show a high degree of structure formation of this RNA already at very low ion-concentrations. Mg²⁺ facilitates additionally a general compaction of the riboswitch towards its proposed active structure. We show that this fold agrees with the earlier suggested secondary structure which included also a long-range tetraloop/tetraloop-receptor like interaction. Metal ion cleavage assays revealed specific Mg²⁺-binding pockets within the moaA riboswitch. These Mg²⁺ binding pockets are good indicators for the potential Moco binding site, since in riboswitches, Mg²⁺ was shown to be necessary to bind phosphate-carrying metabolites. The importance of the phosphate and of other functional groups of Moco is highlighted by binding assays with tetrahydrobiopterin, the reduced and oxygen-sensitive core moiety of Moco. We demonstrate that the general molecular shape of pterin by its own is insufficient for the recognition by the riboswitch.

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.

Antisense-acting riboswitches: A poorly characterized yet important model of transcriptional regulation in prokaryotic organisms

Riboswitches are RNA elements involved in regulating genes that participate in the biosynthesis or transport of essential metabolites. They are characterized by their ability to recognize their target molecules with high affinity and specificity. Riboswitches are commonly cotranscribed with their target genes and are located at the 5' end of their transcriptional units. To date, only two exceptional cases of riboswitches being situated at the 3' end and transcribing in the antisense direction of their regulated genes have been described. The first case involves a SAM riboswitch located at the 3' end of the ubiG-mccB-mccA operon in Clostridium acetobutylicum involved in converting methionine to cysteine. The second case concerns a Cobalamin riboswitch in Listeria monocytogenes that regulates the transcription factor PocR related to this organism's pathogenic process. In almost a decade since the first descriptions of antisense-acting riboswitches, no new examples have been described. In this work, we performed a computational analysis to identify new examples of antisense-acting riboswitches. We found 292 cases in which, according to the available information, we infer that the expected regulation of the riboswitch is consistent with the signaling molecule it senses and the metabolic function of the regulated gene. The metabolic implications of this novel type of regulation are thoroughly discussed.

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.

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

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

Molecular dynamics simulation on the Thermosinus carboxydivorans pfl ZTP riboswitch by ligand binding

Riboswitches are RNA molecules that can regulate gene expression which is affected by ligand-binding during cotranscriptional folding process. However, the role of ligand during the folding is still unclear. In this study, the pfl domain of Thermosinus carboxydivorans ZTP riboswitch was discussed. The ligand is molecule ZMP. We mainly analyzed the change of ZMP-free and ZMP-bound aptamer domain by the dynamics simulation method. Structural features by calculating their RMSD, RMSF, etc. are analyzed. The results demonstrate that the binding domain require the presence of ZMP to maintain a stable fold. It also suggested that ZMP specificly binding to ZTP can generate more hydrogen bonds in the binding domain. Through the calculation of binding free energy decomposition of each nucleotide, molecule ZMP was found to promote the recognition and binding process of ligands by controlling some special nucleotides in the process of ligand binding. At last, the dynamical correlation and components of conformational motions were both applied to explore the effect of molecule ZMP to ZTP riboswitch. In general, ZMP can effectively affect the motions of the pfl riboswitch and facilitate the folding process of the ZTP riboswitch.These results may provide some new ideas for structural changes in riboswitches and their cotranscriptional folding process.

The Biochemical Landscape of Riboswitch Ligands

More than 55 distinct classes of riboswitches that respond to small metabolites or elemental ions have been experimentally validated to date. The ligands sensed by these riboswitches are biased in favor of fundamental compounds or ions that are likely to have been relevant to ancient forms of life, including those that might have populated the "RNA World", which is a proposed biochemical era that predates the evolutionary emergence of DNA and proteins. In the following text, I discuss the various types of ligands sensed by some of the most common riboswitches present in modern bacterial cells and consider implications for ancient biological processes centered on the proven capabilities of these RNA-based sensors. Although most major biochemical aspects of metabolism are represented by known riboswitch classes, there are striking sensory gaps in some key areas. These gaps could reveal weaknesses in the performance capabilities of RNA that might have hampered RNA World evolution, or these could highlight opportunities to discover additional riboswitch classes that sense essential metabolites.

Proteomic analysis of Rhodospirillum rubrum after carbon monoxide exposure reveals an important effect on metallic cofactor biosynthesis

Some carboxydotrophs like Rhodospirillum rubrum are able to grow with CO as their sole source of energy using a Carbone monoxide dehydrogenase (CODH) and an Energy conserving hydrogenase (ECH) to perform anaerobically the so called water-gas shift reaction (WGSR) (CO + H₂O → CO₂ + H₂). Several studies have focused at the biochemical and biophysical level on this enzymatic system and a few OMICS studies on CO metabolism. Knowing that CO is toxic in particular due to its binding to heme iron atoms, and is even considered as a potential antibacterial agent, we decided to use a proteomic approach in order to analyze R. rubrum adaptation in term of metabolism and management of the toxic effect. In particular, this study allowed highlighting a set of proteins likely implicated in ECH maturation, and important perturbations in term of cofactor biosynthesis, especially metallic cofactors. This shows that even this CO tolerant microorganism cannot avoid completely CO toxic effects associated with its interaction with metallic ions. SIGNIFICANCE: This proteomic study highlights the fact that even in a microorganism able to handle carbon monoxide and in some way detoxifying it via the intrinsic action of the carbon monoxide dehydrogenase (CODH), CO has important effects on metal homeostasis, metal cofactors and metalloproteins. These effects are direct or indirect via transcription regulation, and amplified by the high interdependency of cofactors biosynthesis.

Harnessing Escherichia coli for Bio-Based Production of Formate under Pressurized H2 and CO2 Gases

Escherichia coli is a Gram-negative bacterium that is a workhorse for biotechnology. The organism naturally performs a mixed-acid fermentation under anaerobic conditions where it synthesizes formate hydrogenlyase (FHL-1). The physiological role of the enzyme is the disproportionation of formate into H₂ and CO₂. However, the enzyme has been observed to catalyze hydrogenation of CO₂ given the correct conditions, and so it has possibilities in bio-based carbon capture and storage if it can be harnessed as a hydrogen-dependent CO₂ reductase (HDCR). In this study, an E. coli host strain was engineered for the continuous production of formic acid from H₂ and CO₂ during bacterial growth in a pressurized batch bioreactor. Incorporation of tungsten, in place of molybdenum, in FHL-1 helped to impose a degree of catalytic bias on the enzyme. This work demonstrates that it is possible to couple cell growth to simultaneous, unidirectional formate production from carbon dioxide and develops a process for growth under pressurized gases. IMPORTANCE Greenhouse gas emissions, including waste carbon dioxide, are contributing to global climate change. A basket of solutions is needed to steadily reduce emissions, and one approach is bio-based carbon capture and storage. Here, we present our latest work on harnessing a novel biological solution for carbon capture. The Escherichia coli formate hydrogenlyase (FHL-1) was engineered to be constitutively expressed. Anaerobic growth under pressurized H₂ and CO₂ gases was established, and aqueous formic acid was produced as a result. Incorporation of tungsten into the enzyme in place of molybdenum proved useful in poising FHL-1 as a hydrogen-dependent CO₂ reductase (HDCR).

Selection of a self-cleaving ribozyme activated in a chemically and thermally denaturing environment

A new self-cleaving ribozyme was obtained from in vitro selection, displaying site-specific cleavage activity under denaturing conditions, such as high temperatures up to 95 °C, and denaturing solvents (20 M formamide). Adding salt such as Mg2+ also inhibited its activity. The conserved ribozyme sequence was found in the genome of several extremophiles, suggesting its potential biological relevance. This study provides an example of a ribozyme working under exotic conditions, which may expand the application of ribozymes in non-biological environments.

Coenzymes and Their Role in the Evolution of Life

The evolution of coenzymes, or their impact on the origin of life, is fundamental for understanding our own existence. Having established reasonable hypotheses about the emergence of prebiotic chemical building blocks, which were probably created under palaeogeochemical conditions, and surmising that these smaller compounds must have become integrated to afford complex macromolecules such as RNA, the question of coenzyme origin and its relation to the evolution of functional biochemistry should gain new impetus. Many coenzymes have a simple chemical structure and are often nucleotide-derived, which suggests that they may have coexisted with the emergence of RNA and may have played a pivotal role in early metabolism. Based on current theories of prebiotic evolution, which attempt to explain the emergence of privileged organic building blocks, this Review discusses plausible hypotheses on the prebiotic formation of key elements within selected extant coenzymes. In combination with prebiotic RNA, coenzymes may have dramatically broadened early protometabolic networks and the catalytic scope of RNA during the evolution of life.

Regulatory roles of Escherichia coli 5' UTR and ORF-internal RNAs detected by 3' end mapping

Many bacterial genes are regulated by RNA elements in their 5´ untranslated regions (UTRs). However, the full complement of these elements is not known even in the model bacterium Escherichia coli. Using complementary RNA-sequencing approaches, we detected large numbers of 3´ ends in 5´ UTRs and open reading frames (ORFs), suggesting extensive regulation by premature transcription termination. We documented regulation for multiple transcripts, including spermidine induction involving Rho and translation of an upstream ORF for an mRNA encoding a spermidine efflux pump. In addition to discovering novel sites of regulation, we detected short, stable RNA fragments derived from 5´ UTRs and sequences internal to ORFs. Characterization of three of these transcripts, including an RNA internal to an essential cell division gene, revealed that they have independent functions as sRNA sponges. Thus, these data uncover an abundance of cis- and trans-acting RNA regulators in bacterial 5´ UTRs and internal to ORFs.

In most organisms, specific segments of a cell’s genetic information are copied to form single-stranded molecules of various sizes and purposes. Each of these RNA molecules, as they are known, is constructed as a chain that starts at the 5´ end and terminates at the 3´ end.

Certain RNAs carry the information present in a gene, which provides the instructions that a cell needs to build proteins. Some, however, are ‘non-coding’ and instead act to fine-tune the activity of other RNAs. These regulatory RNAs can be separate from the RNAs they control, or they can be embedded in the very sequences they regulate; new evidence also shows that certain regulatory RNAs can act in both ways.

Many regulatory RNAs are yet to be catalogued, even in simple, well-studied species such as the bacterium Escherichia coli. Here, Adams et al. aimed to better characterize the regulatory RNAs present in E. coli by mapping out the 3´ ends of every RNA molecule in the bacterium.

This revealed many new regulatory RNAs and offered insights into where these sequences are located. For instance, the results show that several of these RNAs were embedded within RNA produced from larger genes. Some were nested in coding RNAs, and were parts of a longer RNA sequence that is adjacent to the protein coding segment. Others, however, were present within the instructions that code for a protein.

The work by Adams et al. reveals that regulatory RNAs can be located in unexpected places, and provides a method for identifying them. This can be applied to other types of bacteria, in particular in species with few known RNA regulators.

A rare bacterial RNA motif is implicated in the regulation of the purF gene whose encoded enzyme synthesizes phosphoribosylamine

The Fibro-purF motif is a putative structured noncoding RNA domain that was discovered previously in species of Fibrobacter by using comparative sequence analysis methods. An updated bioinformatics search yielded a total of only 30 unique-sequence representatives, exclusively found upstream of the purF gene that codes for the enzyme amidophosphoribosyltransferase. This enzyme synthesizes the compound 5-phospho-D-ribosylamine (PRA), which is the first committed step in purine biosynthesis. The consensus model for Fibro-purF motif RNAs includes a predicted three-stem junction that carries numerous conserved nucleotide positions within the regions joining the stems. This architecture appears to be of sufficient size and complexity for the formation of the ligand-binding aptamer portion of a riboswitch. In this study, we conducted biochemical analyses of a representative Fibro-purF motif RNA to confirm that the RNA generally folds according to the predicted consensus model. However, due to the instability of PRA, binding of this ligand candidate by the RNA could not be directly assessed. Genetic analyses were used to demonstrate that Fibro-purF motif RNAs regulate gene expression in accordance with predicted PRA concentrations. These findings indicate that Fibro-purF motif RNAs are genetic regulation elements that likely suppress PRA biosynthesis when sufficient levels of this purine precursor are present.

The coenzyme/protein pair and the molecular evolution of life

Covering: up to 2020What was first? Coenzymes or proteins? These questions are archetypal examples of causal circularity in living systems. Classically, this "chicken-and-egg" problem was discussed for the macromolecules RNA, DNA and proteins. This report focuses on coenzymes and cofactors and discusses the coenzyme/protein pair as another example of causal circularity in life. Reflections on the origin of life and hypotheses on possible prebiotic worlds led to the current notion that RNA was the first macromolecule, long before functional proteins and hence DNA. So these causal circularities of living systems were solved by a time travel into the past. To tackle the "chicken-and-egg" problem of the protein-coenzyme pair, this report addresses this problem by looking for clues (a) in the first hypothetical biotic life forms such as protoviroids and the last unified common ancestor (LUCA) and (b) in considerations and evidence of the possible prebiotic production of amino acids and coenzymes before life arose. According to these considerations, coenzymes and cofactors can be regarded as very old molecular players in the origin and evolution of life, and at least some of them developed independently of α-amino acids, which here are evolutionarily synonymous with proteins. Discussions on "chicken-and-egg" problems open further doors to the understanding of evolution.

Diverse Mechanisms and Circuitry for Global Regulation by the RNA-Binding Protein CsrA

The carbon storage regulator (Csr) or repressor of stationary phase metabolites (Rsm) system of Gammaproteobacteria is among the most complex and best-studied posttranscriptional regulatory systems. Based on a small RNA-binding protein, CsrA and homologs, it controls metabolism, physiology, and bacterial lifestyle decisions by regulating gene expression on a vast scale. Binding of CsrA to sequences containing conserved GGA motifs in mRNAs can regulate translation, RNA stability, riboswitch function, and transcript elongation. CsrA governs the expression of dozens of transcription factors and other regulators, further expanding its influence on cellular physiology, and these factors can participate in feedback to the Csr system. Expression of csrA itself is subject to autoregulation via translational inhibition and indirect transcriptional activation. CsrA activity is controlled by small noncoding RNAs (sRNAs), CsrB and CsrC in Escherichia coli, which contain multiple high affinity CsrA binding sites that compete with those of mRNA targets. Transcription of CsrB/C is induced by certain nutrient limitations, cellular stresses, and metabolites, while these RNAs are targeted for degradation by the presence of a preferred carbon source. Consistent with these findings, CsrA tends to activate pathways and processes that are associated with robust growth and repress stationary phase metabolism and stress responses. Regulatory loops between Csr components affect the signaling dynamics of the Csr system. Recently, systems-based approaches have greatly expanded our understanding of the roles played by CsrA, while reinforcing the notion that much remains to be learned about the Csr system.

Nucleic acid therapeutics: a focus on the development of aptamers

INTRODUCTION

Aptamers provide exciting opportunities for the development of specific and targeted therapeutic approaches.

AREAS COVERED

In this review, the authors discuss different therapeutic options available with nucleic acids, including aptamers, focussing on similarities and differences between them. The authors concentrate on case studies with specific aptamers, which exemplify their distinct advantages. The reasons for failure, wherever available, are deliberated upon. Attempts to accelerate the in vitro selection process have been discussed. Challenges with aptamers in terms of their specificity and targeted delivery and strategies to overcome these are described. Examples of precise regulation of systemic half-life of aptamers using antidotes are discussed.

EXPERT OPINION

Despite their nontoxic nature, a variety of reasons limit the therapeutic potential of aptamers in the clinic. The analysis of adverse effects observed with the pegnivacogin/anivamersen pair has highlighted the need to screen for preexisting PEG antibodies in any clinical trial involving pegylated molecules. Surprisingly, and promisingly, the ability of nucleic acid therapeutics to breach the blood brain barrier seems achievable. The recognition of specific motifs, e.g. G-quadruplex in thrombin-binding aptamers, or a 'nucleation' zone while designing aptamer-antidote pairs, is likely to accelerate the discovery of therapeutically efficacious molecules.

A novel riboswitch classification based on imbalanced sequences achieved by machine learning

Riboswitch, a part of regulatory mRNA (50-250nt in length), has two main classes: aptamer and expression platform. One of the main challenges raised during the classification of riboswitch is imbalanced data. That is a circumstance in which the records of a sequences of one group are very small compared to the others. Such circumstances lead classifier to ignore minority group and emphasize on majority ones, which results in a skewed classification. We considered sixteen riboswitch families, to be in accord with recent riboswitch classification work, that contain imbalanced sequences. The sequences were split into training and test set using a newly developed pipeline. From 5460 k-mers (k value 1 to 6) produced, 156 features were calculated based on CfsSubsetEval and BestFirst function found in WEKA 3.8. Statistically tested result was significantly difference between balanced and imbalanced sequences (p < 0.05). Besides, each algorithm also showed a significant difference in sensitivity, specificity, accuracy, and macro F-score when used in both groups (p < 0.05). Several k-mers clustered from heat map were discovered to have biological functions and motifs at the different positions like interior loops, terminal loops and helices. They were validated to have a biological function and some are riboswitch motifs. The analysis has discovered the importance of solving the challenges of majority bias analysis and overfitting. Presented results were generalized evaluation of both balanced and imbalanced models, which implies their ability of classifying, to classify novel riboswitches. The Python source code is available at https://github.com/Seasonsling/riboswitch.

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

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

The biosynthesis of the molybdenum cofactors in Escherichia coli

The biosynthesis of the molybdenum cofactor (Moco) is highly conserved among all kingdoms of life. In all molybdoenzymes containing Moco, the molybdenum atom is coordinated to a dithiolene group present in the pterin-based 6-alkyl side chain of molybdopterin (MPT). In general, the biosynthesis of Moco can be divided into four steps in in bacteria: (i) the starting point is the formation of the cyclic pyranopterin monophosphate (cPMP) from 5'-GTP, (ii) in the second step the two sulfur atoms are inserted into cPMP leading to the formation of MPT, (iii) in the third step the molybdenum atom is inserted into MPT to form Moco and (iv) in the fourth step bis-Mo-MPT is formed and an additional modification of Moco is possible with the attachment of a nucleotide (CMP or GMP) to the phosphate group of MPT, forming the dinucleotide variants of Moco. This review presents an update on the well-characterized Moco biosynthesis in the model organism Escherichia coli including novel discoveries from the recent years.

Evidence that the nadA motif is a bacterial riboswitch for the ubiquitous enzyme cofactor NAD. RNA (NEW YORK, N.Y.) 2019;25:1616-1627. [PMID: 31467147 PMCID: PMC6859854 DOI: 10.1261/rna.072538.119]

The nadA motif is a riboswitch candidate present in various Acidobacteria species that was previously identified by bioinformatic analysis of bacterial DNA data sets. More than 100 unique representatives have been identified exclusively upstream of nadA genes, which code for an enzyme in the biosynthetic pathway of the ubiquitous coenzyme NAD⁺ The architecture of nadA motif RNAs suggests they use structurally similar tandem ligand-binding aptamer domains to control translation initiation. Biochemical analyses reveal that the first domain selectively binds ligands carrying an adenosine 5'-diphosphate (5' ADP) moiety, including NAD⁺ and its reduced form, NADH. Genetic analyses indicate that a tandem nadA motif RNA suppresses gene expression when NAD⁺ is abundant, and that both aptamer domains are required for maximal gene regulation. However, we have not observed selective binding of the nicotinamide moiety of NAD⁺ or binding by the second putative aptamer in vitro, despite sequence and structural similarities between the tandem domains.

Iron-Dependent Regulation of Molybdenum Cofactor Biosynthesis Genes in Escherichia coli

Molybdenum cofactor (Moco) biosynthesis is a complex process that involves the coordinated function of several proteins. In recent years it has become obvious that the availability of iron plays an important role in the biosynthesis of Moco. First, the MoaA protein binds two [4Fe-4S] clusters per monomer. Second, the expression of the moaABCDE and moeAB operons is regulated by FNR, which senses the availability of oxygen via a functional [4Fe-4S] cluster. Finally, the conversion of cyclic pyranopterin monophosphate to molybdopterin requires the availability of the l-cysteine desulfurase IscS, which is a shared protein with a main role in the assembly of Fe-S clusters. In this report, we investigated the transcriptional regulation of the moaABCDE operon by focusing on its dependence on cellular iron availability. While the abundance of selected molybdoenzymes is largely decreased under iron-limiting conditions, our data show that the regulation of the moaABCDE operon at the level of transcription is only marginally influenced by the availability of iron. Nevertheless, intracellular levels of Moco were decreased under iron-limiting conditions, likely based on an inactive MoaA protein in addition to lower levels of the l-cysteine desulfurase IscS, which simultaneously reduces the sulfur availability for Moco production.IMPORTANCE FNR is a very important transcriptional factor that represents the master switch for the expression of target genes in response to anaerobiosis. Among the FNR-regulated operons in Escherichia coli is the moaABCDE operon, involved in Moco biosynthesis. Molybdoenzymes have essential roles in eukaryotic and prokaryotic organisms. In bacteria, molybdoenzymes are crucial for anaerobic respiration using alternative electron acceptors. This work investigates the connection of iron availability to the biosynthesis of Moco and the production of active molybdoenzymes.

Systematic Evaluation of Genetic and Environmental Factors Affecting Performance of Translational Riboswitches

Since their discovery, riboswitches have been attractive tools for the user-controlled regulation of gene expression in bacterial systems. Riboswitches facilitate small molecule mediated fine-tuning of protein expression, making these tools of great use to the synthetic biology community. However, the use of riboswitches is often restricted due to context dependent performance and limited dynamic range. Here, we report the drastic improvement of a previously developed orthogonal riboswitch achieved through in vivo functional selection and optimization of flanking coding and noncoding sequences. The behavior of the derived riboswitches was mapped under a wide array of growth and induction conditions, using a structured Design of Experiments approach. This approach successfully improved the maximal protein expression levels 8.2-fold relative to the original riboswitches, and the dynamic range was improved to afford riboswitch dependent control of 80-fold. The optimized orthogonal riboswitch was then integrated downstream of four endogenous stress promoters, responsive to phosphate starvation, hyperosmotic stress, redox stress, and carbon starvation. These responsive stress promoter-riboswitch devices were demonstrated to allow for tuning of protein expression up to ∼650-fold in response to both environmental and cellular stress responses and riboswitch dependent attenuation. We envisage that these riboswitch stress responsive devices will be useful tools for the construction of advanced genetic circuits, bioprocessing, and protein expression.

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.

The regulation of Moco biosynthesis and molybdoenzyme gene expression by molybdenum and iron in bacteria

The regulation of the operons involved in Moco biosynthesis is dependent on the availability of Fe–S clusters in the cell.

Genome-scale analysis of Acetobacterium bakii reveals the cold adaptation of psychrotolerant acetogens by post-transcriptional regulation

Acetogens synthesize acetyl-CoA via CO2 or CO fixation, producing organic compounds. Despite their ecological and industrial importance, their transcriptional and post-transcriptional regulation has not been systematically studied. With completion of the genome sequence of Acetobacterium bakii (4.28-Mb), we measured changes in the transcriptome of this psychrotolerant acetogen in response to temperature variations under autotrophic and heterotrophic growth conditions. Unexpectedly, acetogenesis genes were highly up-regulated at low temperatures under heterotrophic, as well as autotrophic, growth conditions. To mechanistically understand the transcriptional regulation of acetogenesis genes via changes in RNA secondary structures of 5'-untranslated regions (5'-UTR), the primary transcriptome was experimentally determined, and 1379 transcription start sites (TSS) and 1100 5'-UTR were found. Interestingly, acetogenesis genes contained longer 5'-UTR with lower RNA-folding free energy than other genes, revealing that the 5'-UTRs control the RNA abundance of the acetogenesis genes under low temperature conditions. Our findings suggest that post-transcriptional regulation via RNA conformational changes of 5'-UTRs is necessary for cold-adaptive acetogenesis.

Coordinated regulation of nitrogen fixation and molybdate transport by molybdenum

Biological nitrogen fixation, the reduction of chemically inert dinitrogen to bioavailable ammonia, is a central process in the global nitrogen cycle highly relevant for life on earth. N₂ reduction to NH₃ is catalyzed by nitrogenases exclusively synthesized by diazotrophic prokaryotes. All diazotrophs have a molybdenum nitrogenase containing the unique iron-molybdenum cofactor FeMoco. In addition, some diazotrophs encode one or two alternative Mo-free nitrogenases that are less efficient at reducing N₂ than Mo-nitrogenase. To permit biogenesis of Mo-nitrogenase and other molybdoenzymes when Mo is scarce, bacteria synthesize the high-affinity molybdate transporter ModABC. Generally, Mo supports expression of Mo-nitrogenase genes, while it represses production of Mo-free nitrogenases and ModABC. Since all three nitrogenases and ModABC can reach very high levels at suitable Mo concentrations, tight Mo-mediated control saves considerable resources and energy. This review outlines the similarities and differences in Mo-responsive regulation of nitrogen fixation and molybdate transport in diverse diazotrophs.

The Molybdenum Storage Protein: A soluble ATP hydrolysis-dependent molybdate pump

A continuous FeMo cofactor supply for nitrogenase maturation is ensured in Azotobacter vinelandii by developing a cage-like molybdenum storage protein (MoSto) capable to store ca. 120 molybdate molecules ( <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:msubsup><mml:mtext>MoO</mml:mtext> <mml:mrow><mml:mn>4</mml:mn></mml:mrow> <mml:mrow><mml:mn>2</mml:mn> <mml:mo>-</mml:mo></mml:mrow> </mml:msubsup> </mml:math> ) as discrete polyoxometalate (POM) clusters. To gain mechanistic insight into this process, MoSto was characterized by Mo and ATP/ADP content, structural, and kinetic analysis. We defined three functionally relevant states specified by the presence of both ATP/ADP and POM clusters (MoSto_funct ), of only ATP/ADP (MoSto_basal ) and of neither ATP/ADP nor POM clusters (MoSto_zero ), respectively. POM clusters are only produced when ATP is hydrolyzed to ADP and phosphate. V_max was ca. 13 μmol_phosphate ·min^-1 ·mg^-1 and K_m for molybdate and ATP/Mg²⁺ in the low micromolar range. ATP hydrolysis presumably proceeds at subunit α, inferred from a highly occupied α-ATP/Mg²⁺ and a weaker occupied β-ATP/no Mg²⁺ -binding site found in the MoSto_funct structure. Several findings indicate that POM cluster storage is separated into a rapid ATP hydrolysis-dependent molybdate transport across the protein cage wall and a slow molybdate assembly induced by combined auto-catalytic and protein-driven processes. The cage interior, the location of the POM cluster depot, is locked in all three states and thus not rapidly accessible for molybdate from the outside. Based on V_max , the entire Mo storage process should be completed in less than 10 s but requires, according to the molybdate content analysis, ca. 15 min. Long-time incubation of MoSto_basal with nonphysiological high molybdate amounts implicates an equilibrium in and outside the cage and POM cluster self-formation without ATP hydrolysis. DATABASES: The crystal structures MoSto in the MoSto-F6, MoSto-F7, MoSto_basal , MoSto_zero , and MoSto-F1_vitro states were deposited to PDB under the accession numbers PDB 6GU5, 6GUJ, 6GWB, 6GWV, and 6GX4.

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.

ykkC riboswitches employ an add-on helix to adjust specificity for polyanionic ligands

The ykkC family of bacterial riboswitches combines several widespread classes that have similar secondary structures and consensus motifs but control different genes in response to different cellular metabolites. Here we report the crystal structures of two distinct ykkC riboswitches specifically bound to their cognate ligands ppGpp, a second messenger involved in stress response, or PRPP, a precursor in purine biosynthesis. Both RNAs adopt similar structures and contain a conserved core previously observed in the guanidine-specific ykkC riboswitch. However, ppGpp and PRPP riboswitches uniquely employ an additional helical element that joins the ends of the ligand-sensing domains and creates a tunnel for direct and Mg²⁺-mediated binding of ligands. Mutational and footprinting experiments highlight the importance of conserved nucleotides forming the tunnel and long-distance contacts for ligand binding and genetic response. Our work provides new insights into the specificity of riboswitches and gives a unique opportunity for future studies of RNA evolution.

Bacterial riboswitches and RNA thermometers: Nature and contributions to pathogenesis

Bacterial pathogens are always challenged by fluctuations of chemical and physical parameters that pose serious threats to cellular integrity and metabolic status. Sudden deprivation of nutrients or key metabolites, changes in surrounding pH, and temperature shifts are the most important examples of such parameters. To elicit a proper response to such fluctuations, bacterial cells coordinate the expression of parameter-relevant genes. Although protein-mediated control of gene expression is well appreciated since many decades, RNA-based regulation has been discovered in early 2000s as a parallel level of regulation. Small regulatory RNAs have emerged as one of the most widespread and important gene regulatory systems in bacteria with rare representatives found in Archaea and Eukarya. Riboswitches and thermosensors are cis-encoded RNA regulatory elements that employ different mechanisms to regulate the expression of related genes controlling key metabolic pathways and genes of temperature relevant proteins including virulence factors. The extent of RNA contributions to gene regulation is not completely known even in well-studied models such E. coli and B. subtilis. In depth understanding of riboswitches is promising for opportunity to discover a narrow spectrum antibacterial drugs that target riboswitches of essential metabolic pathways.

Engineering riboswitch in L. major: From prediction to conceptualization

Post-transcriptional gene regulation is a vital process to regulate expression of the key genes in the eukaryotic cell. Such processes are essential for pathogens which reside inside the host cell. One such pathogen is Leishmania major, which causes cutaneous leishmaniasis. The parasite lives inside the macrophages of mammalian host (mostly human). Inside the macrophage, Leishmania genes show complex host-pathogen interaction regulating a plethora of gene expression. Till date, most of the studies have shown this kind of regulation with respect to the host macrophages. Here, based on an extensive in silico analysis, we have hypothesized a novel Theophylline binding riboswitch mediated post-transcriptional regulation of a gene i.e. RNA Polymerase III subunit1 (Lmjf_09_1060), an essential gene for the parasite's survival both in its promastigote as well as in its amastigote form. Later, we have conceptualized the working of the identified putative Theophylline binding riboswitch cassette in in vitro using E. coli based reporter assay, wherein, a reporter gene (eGFP) is used instead of RNA Polymerase III subunit1 gene and apparently have shown the downregulation of the reporter gene (eGFP) expression under the influence of in silico identified Theophylline binding riboswitch.

Tandem riboswitches form a natural Boolean logic gate to control purine metabolism in bacteria

Gene control systems sometimes interpret multiple signals to set the expression levels of the genes they regulate. In rare instances, ligand-binding riboswitch aptamers form tandem arrangements to approximate the function of specific two-input Boolean logic gates. Here, we report the discovery of riboswitch aptamers for phosphoribosyl pyrophosphate (PRPP) that naturally exist either in singlet arrangements, or occur in tandem with guanine aptamers. Tandem guanine-PRPP aptamers can bind the target ligands, either independently or in combination, to approximate the function expected for an IMPLY Boolean logic gate to regulate transcription of messenger RNAs for de novo purine biosynthesis in bacteria. The existence of sophisticated all-RNA regulatory systems that sense two ancient ribonucleotide derivatives to control synthesis of RNA molecules supports the hypothesis that RNA World organisms could have managed a complex metabolic state without the assistance of protein regulatory factors.

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.

Riboswitch regulation in cyanobacteria is independent of their habitat adaptations

Cyanobacteria are one of the ancient bacterial species occupying a variety of habitats with diverse metabolic preferences. RNA regulators like riboswitches play significant role in controlling the gene expression in prokaryotes. The taxonomic distribution of riboswitches suggests that they might be one of the oldest mechanisms of gene control system. In this paper, we analyzed the distribution of different riboswitch families in various cyanobacterial genomes. It was observed that only four riboswitch classes were abundant in cyanobacteria, B12-element (Cob)/AdoCbl/AdoCbl-variant riboswitch being the most abundant. The analysis suggests that riboswitch mode of regulation is present in cyanobacterial species irrespective of their habitat types. A large number of unidentified genes regulated by riboswitches listed in this analysis indicate the wide range of targets for these riboswitch families. The analysis revealed a large number of genes regulated by riboswitches which may assist in elaborating the diversity among the cyanobacterial species.

Using RNA Sequence and Structure for the Prediction of Riboswitch Aptamer: A Comprehensive Review of Available Software and Tools

RNA molecules are essential players in many fundamental biological processes. Prokaryotes and eukaryotes have distinct RNA classes with specific structural features and functional roles. Computational prediction of protein structures is a research field in which high confidence three-dimensional protein models can be proposed based on the sequence alignment between target and templates. However, to date, only a few approaches have been developed for the computational prediction of RNA structures. Similar to proteins, RNA structures may be altered due to the interaction with various ligands, including proteins, other RNAs, and metabolites. A riboswitch is a molecular mechanism, found in the three kingdoms of life, in which the RNA structure is modified by the binding of a metabolite. It can regulate multiple gene expression mechanisms, such as transcription, translation initiation, and mRNA splicing and processing. Due to their nature, these entities also act on the regulation of gene expression and detection of small metabolites and have the potential to helping in the discovery of new classes of antimicrobial agents. In this review, we describe software and web servers currently available for riboswitch aptamer identification and secondary and tertiary structure prediction, including applications.

Multi-metal-dependent nucleic acid enzymes

Nucleic acid enzymes require metal ions for activity, and many recently discovered enzymes can use multiple metals, either binding to the scissile phosphate or also playing an allosteric role.

Determination of the Genome and Primary Transcriptome of Syngas Fermenting Eubacterium limosum ATCC 8486

Autotrophic conversion of CO₂ to value-added biochemicals has received considerable attention as a sustainable route to replace fossil fuels. Particularly, anaerobic acetogenic bacteria are naturally capable of reducing CO₂ or CO to various metabolites. To fully utilize their biosynthetic potential, an understanding of acetogenesis-related genes and their regulatory elements is required. Here, we completed the genome sequence of the syngas fermenting Eubacterium limosum ATCC 8486 and determined its transcription start sites (TSS). We constructed a 4.4 Mb long circular genome with a GC content of 47.2% and 4,090 protein encoding genes. To understand the transcriptional and translational regulation, the primary transcriptome was augmented, identifying 1,458 TSSs containing a high pyrimidine (T/C) and purine nucleotide (A/G) content at the −1 and +1 position, respectively, along with 1,253 5′-untranslated regions, and principal promoter elements such as −10 (TATAAT) and −35 (TTGACA), and Shine-Dalgarno motifs (GGAGR). Further analysis revealed 93 non-coding RNAs, including one for potential transcriptional regulation of the hydrogenase complex via interaction with molybdenum or tungsten cofactors, which in turn controls formate dehydrogenase activity of the initial step of Wood-Ljungdahl pathway. Our results provide comprehensive genomic information for strain engineering to enhance the syngas fermenting capacity of acetogenic bacteria.

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.

Allosteric pathways in tetrahydrofolate sensing riboswitch with dynamics correlation network

Riboswitches are cis-acting genetic control elements. Due to their fundamental importance in bacteria gene regulation, they have been proposed as antibacterial drug targets. Tetrahydrofolate (THF) is an essential cofactor of one-carbon transfer reactions and downregulates the expression of downstream genes. However, information on how to transfer from the binding site of THF to the expression platform is still unavailable. Herein, a nucleotide/nucleotide dynamics correlation network based on an all-atom molecular dynamic simulation was used to reveal the regulation mechanism of a THF-sensing riboswitch. Shortest pathway analysis based on the network illustrates that there is an allosteric pathway through the P2 helix to the pseudoknot, then to the P1 helix in the THF-riboswitch. Thus the hypothesis of "THF-binding induced allosteric switching" was proposed and evaluated using THF and pseudoknot weakened experiments. Furthermore, a possible allosteric pathway of C30-C31-G33-A34-G35-G36-G37-A38-G48-G47-U46-A90-U91-C92-G93-C94-G95-C96 was identified and confirmed through the perturbation of the network. The proposed allosteric mechanism and the underlying allosteric pathway provide fundamental insights for the regulation of THF sensing riboswitches.

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.