Riboswitches as versatile gene control elements

Regulatory helix plays a key role in genetic ON-OFF switching for the 2'-deoxyguanosine sensing mRNA element

Transcriptional riboswitches, noncoding mRNA elements that operate in cis to regulate gene expression, have a promising potential in medicine, synthetic biology and directed evolution. They bind to cellular metabolites or metal ions with high specificity, leading to conformational rearrangements that facilitate the activation or premature termination of transcription for downstream genes. This elegant mechanism for feedback regulation of metabolic pathways has been identified in prokaryotes and a few in eukaryotes. Our chemical probing of the 2'-deoxyguanosine (2'-dG)-sensing riboswitch demonstrates that the overall conformational state of the full-length riboswitch (dGsw-fl) is unresponsive to the 2'-dG. Although binding proceeds as expected, dGsw-fl exclusively populates an OFF state of transcriptional inhibition. We chemically probed the structure of a known dGsw transcriptional intermediate (dGsw-int) to evaluate the possibility of a cotranscriptional regulatory role. Interestingly, apo dGsw-int adopts an alternative conformation in which a stable anti-terminator helix is formed, leading to an ON state where transcription can proceed. In the presence of 2'-dG, this anti-terminator helix is destabilized to produce a conformation reminiscent of the full-length, OFF-state dGsw. Using a fluorescence quenching assay, we demonstrate that binding 2'-dG to early transcriptional intermediates can inhibit the formation of the anti-terminator helix, locking dGsw in an OFF state. These data suggest that metabolite sensing occurs during a brief window of time between the synthesis of two transcriptional intermediates. Our studies indicate that dGsw does not function as a binary ON-OFF switch, but instead fine-tunes the transcription of downstream genes during RNA synthesis using key intermediates.

Engineering oxypurinol-responsive riboswitches based on bacterial xanthine aptamers for gene expression control in mammalian cell culture

Riboswitch-mediated control of gene expression without the interference of potentially immunogenic proteins is a promising approach for the development of tailor-made tools for biological research and the advancement of gene therapies. However, the current selection of applicable ligands for synthetic riboswitches is limited and strategies have mostly relied on de novo selection of aptamers. Here, we show that the bacterial xanthine I riboswitch aptamer recognizes oxypurinol, the active metabolite of the widely prescribed anti-gout drug allopurinol (Zyloprim®). We have characterized the aptamer/oxypurinol interaction and present a crystal structure of the oxypurinol-bound aptamer, revealing a binding mode similar to that of the cognate ligand xanthine. We then constructed artificial oxypurinol-responsive riboswitches that showed functionality in human cells. By optimizing splicing-based oxypurinol riboswitches using three different strategies, transgene expression could be induced by >100-fold. In summary, we have developed recombinant RNA switches enabling on-demand regulation of gene expression in response to an established and safe drug.

Synthetic Dual-Input Hybrid Riboswitches─Optimized Genetic Regulators in Yeast

Synthetic riboswitches, genetic regulatory elements composed entirely of RNA, have been engineered to control a variety of mechanisms at the level of both transcription and translation in all domains of life. The efficiency of riboswitch regulation can be increased by inserting two of them into an mRNA sequence in close proximity, resulting in a tandem riboswitch. The tandem state results in improved regulation beyond that of a single riboswitch by allowing both binding pockets to contribute to a higher dynamic range. The focus of this study was to create a novel tandem riboswitch design by integrating the binding pockets of two different riboswitches into one continuous structure, thereby creating a dual-input hybrid riboswitch. These hybrids remain compact in size with a shorter sequence length compared to a tandem riboswitch, while taking advantage of the binding pockets and scaffold sequences provided by both parental riboswitches. Through rational design, hybrid constructs derived from the combination of tetracycline-, tobramycin-, neomycin-, and paromomycin-binding riboswitches were engineered that significantly increase the dynamic range (e.g., from 14- to 36-fold for tobramycin) while increasing their expression levels in the absence of ligand (e.g., 28% to 68% expression for tetracycline). This study expands the toolbox of synthetic riboswitches and establishes general design guidelines applicable to similar riboswitches. Additionally, the dual-input state makes hybrid riboswitches an interesting target for the design of genetic regulators following Boolean logic.

Linker-Mediated Inactivation of the SAM-II Domain in the Tandem SAM-II/SAM-V Riboswitch

Tandem SAM-II/SAM-V riboswitch belongs to a class of riboswitches found in the marine bacterium 'Candidatus Pelagibacter ubique'. Previous studies have demonstrated that these riboswitches have the potential for digital modulation of gene expression at both the transcriptional and translational levels. In this study, we investigate the conformational changes in the tandem SAM-II/SAM-V riboswitch binding to S-adenosylmethionine (SAM) using selective 2'-hydroxyl acylation analyzed by the primer extension (SHAPE) assay, small-angle X-ray scattering (SAXS), and oligos depressing probing. Our findings reveal that the linker between SAM-II/SAM-V aptamers blocks the SAM response of the SAM-II domain. This result proposes a new mechanism for gene expression regulation, where the ligand-binding functions of tandem riboswitches can be selectively masked or released through a linker.

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

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

Opportunities for Riboswitch Inhibition by Targeting Co-Transcriptional RNA Folding Events

Antibiotic resistance is a critical global health concern, causing millions of prolonged bacterial infections every year and straining our healthcare systems. Novel antibiotic strategies are essential to combating this health crisis and bacterial non-coding RNAs are promising targets for new antibiotics. In particular, a class of bacterial non-coding RNAs called riboswitches has attracted significant interest as antibiotic targets. Riboswitches reside in the 5'-untranslated region of an mRNA transcript and tune gene expression levels in cis by binding to a small-molecule ligand. Riboswitches often control expression of essential genes for bacterial survival, making riboswitch inhibitors an exciting prospect for new antibacterials. Synthetic ligand mimics have predominated the search for new riboswitch inhibitors, which are designed based on static structures of a riboswitch's ligand-sensing aptamer domain or identified by screening a small-molecule library. However, many small-molecule inhibitors that bind an isolated riboswitch aptamer domain with high affinity in vitro lack potency in vivo. Importantly, riboswitches fold and respond to the ligand during active transcription in vivo. This co-transcriptional folding is often not considered during inhibitor design, and may explain the discrepancy between a low Kd in vitro and poor inhibition in vivo. In this review, we cover advances in riboswitch co-transcriptional folding and illustrate how intermediate structures can be targeted by antisense oligonucleotides-an exciting new strategy for riboswitch inhibitor design.

Evaluating Performance of Different RNA Secondary Structure Prediction Programs Using Self-cleaving Ribozymes

Accurate identification of the correct, biologically relevant RNA structures is critical to understanding various aspects of RNA biology since proper folding represents the key to the functionality of all types of RNA molecules and plays pivotal roles in many essential biological processes. Thus, a plethora of approaches have been developed to predict, identify, or solve RNA structures based on various computational, molecular, genetic, chemical, or physicochemical strategies. Purely computational approaches hold distinct advantages over all other strategies in terms of the ease of implementation, time, speed, cost, and throughput, but they strongly underperform in terms of accuracy that significantly limits their broader application. Nonetheless, the advantages of these methods led to a steady development of multiple in silico RNA secondary structure prediction approaches including recent deep learning-based programs. Here, we compared the accuracy of predictions of biologically relevant secondary structures of dozens of self-cleaving ribozyme sequences using seven in silico RNA folding prediction tools with tasks of varying complexity. We found that while many programs performed well in relatively simple tasks, their performance varied significantly in more complex RNA folding problems. However, in general, a modern deep learning method outperformed the other programs in the complex tasks in predicting the RNA secondary structures, at least based on the specific class of sequences tested, suggesting that it may represent the future of RNA structure prediction algorithms.

The effect of the loop on the thermodynamic and kinetic of single base pair in pseudoknot

RNA pseudoknots are RNA molecules with specialized three-dimensional structures that play important roles in various biological processes. To understand the functions and mechanisms of pseudoknots, it is essential to elucidate their structures and folding pathways. The most fundamental step in RNA folding is the opening and closing of a base pair. The effect of flexible loops on the base pair in pseudoknots remains unclear. In this work, we use molecular dynamics simulations and Markov state model to study the configurations, thermodynamic and kinetic of single base pair in pseudoknots. We find that the presence of the loop leads to a trap state. In addition, the rate-limiting step for the formation of base pair is the disruption of the trap state, rather than the open state to the closed state, which is quite different from the previous studies on non-pseudoknot RNA. For the thermodynamic parameters in pseudoknots, we find that the entropy difference upon opening the base pair between this simulation and the nearest-neighbor model results from the different entropy of different lengths of loop in solution. The thermodynamic parameters of the stack in pseudoknot are close to the nearest-neighbor parameters. The bases on the loop have different distribution patterns in different states, and the slow transition states of the loop are determined by the orientation of the bases.

SPRank─A Knowledge-Based Scoring Function for RNA-Ligand Pose Prediction and Virtual Screening

The growing interest in RNA-targeted drugs underscores the need for computational modeling of interactions between RNA molecules and small compounds. Having a reliable scoring function for RNA-ligand interactions is essential for effective computational drug screening. An ideal scoring function should not only predict the native pose for ligand binding but also rank the affinity of the binding for different ligands. However, existing scoring functions are primarily designed to predict the native binding modes for a given RNA-ligand pair and have not been thoroughly assessed for virtual screening purposes. In this paper, we introduce SPRank, a combination of machine-learning and knowledge-based scoring functions developed through a weighted iterative approach, specifically designed to tackle both binding mode prediction and virtual screening challenges. Our approach incorporates third-party docking software, such as rDock and AutoDock Vina, to sample flexible ligands against an ensemble of RNA structures, capturing the conformational flexibility of both the RNA and the ligand. Through rigorous testing, SPRank demonstrates improved performance compared to the tested scoring functions across four test sets comprising 122, 42, 55, and 71 nucleic acid-ligand complexes. Furthermore, SPRank exhibits improved performance in virtual screening tests targeting the HIV-1 TAR ensemble, which highlights its advantage in drug discovery. These results underscore the advantages of SPRank as a potentially promising tool for the RNA-targeted drug design. The source code of SPRank and the data sets are freely accessible at https://github.com/Vfold-RNA/SPRank.

RNA-Based Sensor Systems for Affordable Diagnostics in the Age of Pandemics

In the era of the COVID-19 pandemic, the significance of point-of-care (POC) diagnostic tools has become increasingly vital, driven by the need for quick and precise virus identification. RNA-based sensors, particularly toehold sensors, have emerged as promising candidates for POC detection systems due to their selectivity and sensitivity. Toehold sensors operate by employing an RNA switch that changes the conformation when it binds to a target RNA molecule, resulting in a detectable signal. This review focuses on the development and deployment of RNA-based sensors for POC viral RNA detection with a particular emphasis on toehold sensors. The benefits and limits of toehold sensors are explored, and obstacles and future directions for improving their performance within POC detection systems are presented. The use of RNA-based sensors as a technology for rapid and sensitive detection of viral RNA holds great potential for effectively managing (dealing/coping) with present and future pandemics in resource-constrained settings.

Utilizing 5' UTR Engineering Enables Fine-Tuning of Multiple Genes within Operons to Balance Metabolic Flux in Bacillus subtilis

The application of synthetic biology tools to modulate gene expression to increase yield has been thoroughly demonstrated as an effective and convenient approach in industrial production. In this study, we employed a high-throughput screening strategy to identify a 5' UTR sequence from the genome of B. subtilis 168. This sequence resulted in a 5.8-fold increase in the expression level of EGFP. By utilizing the 5' UTR sequence to overexpress individual genes within the rib operon, it was determined that the genes ribD and ribAB serve as rate-limiting enzymes in the riboflavin synthesis pathway. Constructing a 5' UTR library to regulate EGFP expression resulted in a variation range in gene expression levels exceeding 100-fold. Employing the same 5' UTR library to regulate the expression of EGFP and mCherry within the operon led to a change in the expression ratio of these two genes by over 10,000-fold. So, employing a 5' UTR library to modulate the expression of the rib operon gene and construct a synthetic rib operon resulted in a 2.09-fold increase in riboflavin production. These results indicate that the 5' UTR sequence identified and characterized in this study can serve as a versatile synthetic biology toolkit for achieving complex metabolic network reconstruction. This toolkit can facilitate the fine-tuning of gene expression to produce target products.

T7 phage-assisted evolution of riboswitches using error-prone replication and dual selection

Leveraging riboswitches, non-coding mRNA fragments pivotal to gene regulation, poses a challenge in effectively selecting and enriching these functional genetic sensors, which can toggle between ON and OFF states in response to their cognate inducers. Here, we show our engineered phage T7, enabling the evolution of a theophylline riboswitch. We have replaced T7's DNA polymerase with a transcription factor controlled by a theophylline riboswitch and have created two types of host environments to propagate the engineered phage. Both types host an error-prone T7 DNA polymerase regulated by a T7 promoter along with another critical gene-either cmk or pifA, depending on the host type. The cmk gene is necessary for T7 replication and is used in the first host type for selection in the riboswitch's ON state. Conversely, the second host type incorporates the pifA gene, leading to abortive T7 infections and used for selection in the riboswitch's OFF state. This dual-selection system, termed T7AE, was then applied to a library of 65,536 engineered T7 phages, each carrying randomized riboswitch variants. Through successive passage in both host types with and without theophylline, we observed an enrichment of phages encoding functional riboswitches that conferred a fitness advantage to the phage in both hosts. The T7AE technique thereby opens new pathways for the evolution and advancement of gene switches, including non-coding RNA-based switches, setting the stage for significant strides in synthetic biology.

Compact RNA sensors for increasingly complex functions of multiple inputs

Designing single molecules that compute general functions of input molecular partners represents a major unsolved challenge in molecular design. Here, we demonstrate that high-throughput, iterative experimental testing of diverse RNA designs crowdsourced from Eterna yields sensors of increasingly complex functions of input oligonucleotide concentrations. After designing single-input RNA sensors with activation ratios beyond our detection limits, we created logic gates, including challenging XOR and XNOR gates, and sensors that respond to the ratio of two inputs. Finally, we describe the OpenTB challenge, which elicited 85-nucleotide sensors that compute a score for diagnosing active tuberculosis, based on the ratio of products of three gene segments. Building on OpenTB design strategies, we created an algorithm Nucleologic that produces similarly compact sensors for the three-gene score based on RNA and DNA. These results open new avenues for diverse applications of compact, single molecule sensors previously limited by design complexity.

Fluorescent Ligand Equilibrium Displacement: A High-Throughput Method for Identification of FMN Riboswitch-Binding Small Molecules

Antibiotic resistance remains a pressing global concern, with most antibiotics targeting the bacterial ribosome or a limited range of proteins. One class of underexplored antibiotic targets is bacterial riboswitches, structured RNA elements that regulate key biosynthetic pathways by binding a specific ligand. We developed a methodology termed Fluorescent Ligand Equilibrium Displacement (FLED) to rapidly discover small molecules that bind the flavin mononucleotide (FMN) riboswitch. FLED leverages intrinsically fluorescent FMN and the quenching effect on RNA binding to create a label-free, in vitro method to identify compounds that can bind the apo population of riboswitch in a system at equilibrium. The response difference between known riboswitch ligands and controls demonstrates the robustness of the method for high-throughput screening. An existing drug discovery library that was screened using FLED resulted in a final hit rate of 0.67%. The concentration response of each hit was determined and revealed a variety of approximate effective concentration values. Our preliminary screening data support the use of FLED to identify small molecules for medicinal chemistry development as FMN riboswitch-targeted antibiotic compounds. This robust, label-free, and cell-free method offers a strong alternative to other riboswitch screening methods and can be adapted to a variety of laboratory setups.

Rock, scissors, paper: How RNA structure informs function

RNA can fold back on itself to adopt a wide range of structures. These range from relatively simple hairpins to intricate 3D folds and can be accompanied by regulatory interactions with both metabolites and macromolecules. The last 50 yr have witnessed elucidation of an astonishing array of RNA structures including transfer RNAs, ribozymes, riboswitches, the ribosome, the spliceosome, and most recently entire RNA structuromes. These advances in RNA structural biology have deepened insight into fundamental biological processes including gene editing, transcription, translation, and structure-based detection and response to temperature and other environmental signals. These discoveries reveal that RNA can be relatively static, like a rock; that it can have catalytic functions of cutting bonds, like scissors; and that it can adopt myriad functional shapes, like paper. We relate these extraordinary discoveries in the biology of RNA structure to the plant way of life. We trace plant-specific discovery of ribozymes and riboswitches, alternative splicing, organellar ribosomes, thermometers, whole-transcriptome structuromes and pan-structuromes, and conclude that plants have a special set of RNA structures that confer unique types of gene regulation. We finish with a consideration of future directions for the RNA structure-function field.

Biological Catalysis and Information Storage Have Relied on N-Glycosyl Derivatives of β-D-Ribofuranose since the Origins of Life

Most naturally occurring nucleotides and nucleosides are N-glycosyl derivatives of β-d-ribose. These N-ribosides are involved in most metabolic processes that occur in cells. They are essential components of nucleic acids, forming the basis for genetic information storage and flow. Moreover, these compounds are involved in numerous catalytic processes, including chemical energy production and storage, in which they serve as cofactors or coribozymes. From a chemical point of view, the overall structure of nucleotides and nucleosides is very similar and simple. However, their unique chemical and structural features render these compounds versatile building blocks that are crucial for life processes in all known organisms. Notably, the universal function of these compounds in encoding genetic information and cellular catalysis strongly suggests their essential role in the origins of life. In this review, we summarize major issues related to the role of N-ribosides in biological systems, especially in the context of the origin of life and its further evolution, through the RNA-based World(s), toward the life we observe today. We also discuss possible reasons why life has arisen from derivatives of β-d-ribofuranose instead of compounds based on other sugar moieties.

"Toehold Switches; a foothold for Synthetic Biology"

Toehold switches are de novo designed riboregulators that contain two RNA components interacting through linear-linear RNA interactions, regulating the gene expression. These are highly versatile, exhibit excellent orthogonality, wide dynamic range, and are highly programmable, so can be used for various applications in synthetic biology. In this review, we summarized and discussed the design characteristics and benefits of toehold switch riboregulators over conventional riboregulators. We also discussed applications and recent advancements of toehold switch riboregulators in various fields like gene editing, DNA nanotechnology, translational repression, and diagnostics (detection of microRNAs and some pathogens). Toehold switches, therefore, furnished advancement in synthetic biology applications in various fields with their prominent features.

Over supplementation with vitamin B12 alters microbe-host interactions in the gut leading to accelerated Citrobacter rodentium colonization and pathogenesis in mice

BACKGROUND

Vitamin B12 supplements typically contain doses that far exceed the recommended daily amount, and high exposures are generally considered safe. Competitive and syntrophic interactions for B12 exist between microbes in the gut. Yet, to what extent excessive levels contribute to the activities of the gut microbiota remains unclear. The objective of this study was to evaluate the effect of B12 on microbial ecology using a B12 supplemented mouse model with Citrobacter rodentium, a mouse-specific pathogen. Mice were fed a standard chow diet and received either water or water supplemented with B12 (cyanocobalamin: ~120 μg/day), which equates to approximately 25 mg in humans. Infection severity was determined by body weight, pathogen load, and histopathologic scoring. Host biomarkers of inflammation were assessed in the colon before and after the pathogen challenge.

RESULTS

Cyanocobalamin supplementation enhanced pathogen colonization at day 1 (P < 0.05) and day 3 (P < 0.01) postinfection. The impact of B12 on gut microbial communities, although minor, was distinct and attributed to the changes in the Lachnospiraceae populations and reduced alpha diversity. Cyanocobalamin treatment disrupted the activity of the low-abundance community members of the gut microbiota. It enhanced the amount of interleukin-12 p40 subunit protein (IL12/23p40; P < 0.001) and interleukin-17a (IL-17A; P < 0.05) in the colon of naïve mice. This immune phenotype was microbe dependent, and the response varied based on the baseline microbiota. The cecal metatranscriptome revealed that excessive cyanocobalamin decreased the expression of glucose utilizing genes by C. rodentium, a metabolic attribute previously associated with pathogen virulence.

CONCLUSIONS

Oral vitamin B12 supplementation promoted C. rodentium colonization in mice by altering the activities of the Lachnospiraceae populations in the gut. A lower abundance of select Lachnospiraceae species correlated to higher p40 subunit levels, while the detection of Parasutterella exacerbated inflammatory markers in the colon of naïve mice. The B12-induced change in gut ecology enhanced the ability of C. rodentium colonization by impacting key microbe-host interactions that help with pathogen exclusion. This research provides insight into how B12 impacts the gut microbiota and highlights potential consequences of disrupting microbial B12 competition/sharing through over-supplementation. Video Abstract.

Genetic switches based on nucleic acid strand displacement

Toehold-mediated strand displacement (TMSD) is an isothermal switching process that enables the sequence-programmable and reversible conversion of DNA or RNA strands between single- and double-stranded conformations or other secondary structures. TMSD processes have already found widespread application in DNA nanotechnology, where they are used to drive DNA-based molecular devices or for the realization of synthetic biochemical computing circuits. Recently, researchers have started to employ TMSD also for the control of RNA-based gene regulatory processes in vivo, in particular in the context of synthetic riboregulators and conditional guide RNAs for CRISPR/Cas. Here, we provide a review over recent developments in this emerging field and discuss the opportunities and challenges for such systems in in vivo applications.

Structural and Energetic Features of Base-Base Stacking Contacts in RNA

Nucleobase π-π stacking is one of the crucial organizing interactions within three-dimensional (3D) RNA architectures. Characterizing the structural variability of these contacts in RNA crystal structures will help delineate their subtleties and their role in determining function. This analysis of different stacking geometries found in RNA X-ray crystal structures is the largest such survey to date; coupled with quantum-mechanical calculations on typical representatives of each possible stacking arrangement, we determined the distribution of stacking interaction energies. A total of 1,735,481 stacking contacts, spanning 359 of the 384 theoretically possible distinct stacking geometries, were identified. Our analysis reveals preferential occurrences of specific consecutive stacking arrangements in certain regions of RNA architectures. Quantum chemical calculations suggest that 88 of the 359 contacts possess intrinsically stable stacking geometries, whereas the remaining stacks require the RNA backbone or surrounding macromolecular environment to force their formation and maintain their stability. Our systematic analysis of π-π stacks in RNA highlights trends in the occurrence and localization of these noncovalent interactions and may help better understand the structural intricacies of functional RNA-based molecular architectures.

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.

Riboswitches as therapeutic targets: promise of a new era of antibiotics

INTRODUCTION

The growth of antibiotic resistance among bacterial pathogens is an impending global threat that can only be averted through the development of novel antibacterial drugs. A promising answer could be the targeting of riboswitches, structured RNA elements found almost exclusively in bacteria.

AREAS COVERED

This review examines the potential of riboswitches as novel antibacterial drug targets. The limited mechanisms of action of currently available antibiotics are summarized, followed by a delineation of the functional mechanisms of riboswitches. We then discuss the potential for developing novel approaches that target paradigmatic riboswitches in the context of their bacterial gene expression machinery.

EXPERT OPINION

We highlight potential advantages of targeting riboswitches in their functional form, embedded within gene expression complexes critical for bacterial survival. We emphasize the benefits of this approach, including potentially higher species specificity and lower side effects.

Nucleic acid strand displacement - from DNA nanotechnology to translational regulation

Nucleic acid strand displacement reactions involve the competition of two or more DNA or RNA strands of similar sequence for binding to a complementary strand, and facilitate the isothermal replacement of an incumbent strand by an invader. The process can be biased by augmenting the duplex comprising the incumbent with a single-stranded extension, which can act as a toehold for a complementary invader. The toehold gives the invader a thermodynamic advantage over the incumbent, and can be programmed as a unique label to activate a specific strand displacement process. Toehold-mediated strand displacement processes have been extensively utilized for the operation of DNA-based molecular machines and devices as well as for the design of DNA-based chemical reaction networks. More recently, principles developed initially in the context of DNA nanotechnology have been applied for the de novo design of gene regulatory switches that can operate inside living cells. The article specifically focuses on the design of RNA-based translational regulators termed toehold switches. Toehold switches utilize toehold-mediated strand invasion to either activate or repress translation of an mRNA in response to the binding of a trigger RNA molecule. The basic operation principles of toehold switches will be discussed as well as their applications in sensing and biocomputing. Finally, strategies for their optimization will be described as well as challenges for their operation in vivo.

5'UTR sequences influence protein levels in Escherichia coli by regulating translation initiation and mRNA stability

A set of 41 synthetic 5'UTRs with different theoretical translation initiation rates were generated to explore the role of 5'UTRs in the regulation of protein levels in Escherichia coli. The roles of the synthetic 5'UTRs in regulating the expression of different reporter genes were analyzed in vivo. Protein levels varied substantially between the different constructs but for most of the 5'UTRs, protein levels were not correlated with theoretical translation initiation rates. Large variations in mRNA concentrations were measured with the different 5'UTRs even though the same concentration of transcription inducer was used in each case. 5'UTRs were also found to strongly affect mRNA stability, and these changes in mRNA stability often contributed to observed differences in mRNA concentration. Unexpectedly, the effect of the 5'UTRs on mRNA half-lives was found to vary depending on the downstream reporter gene. These results clearly demonstrate that 5'UTRs contribute to gene expression regulation at the level of translation initiation and of mRNA stability, to an extent that depends on the nature of the downstream gene.

Sarcin/Ricin Domain RNA Retains Its Structure Better Than A-RNA in Adaptively Biased Molecular Dynamics Simulations

Less than one in thirty of the RNA sequences transcribed in humans are translated into protein. The noncoding RNA (ncRNA) functions in catalysis, structure, regulation, and more. However, for the most part, these functions are poorly characterized. RNA is modular and described by motifs that include helical A-RNA with canonical Watson-Crick base-pairing as well as structures with only noncanonical base pairs. Understanding the structure and dynamics of motifs will aid in deciphering functions of specific ncRNAs. We present computational studies on a standard sarcin/ricin domain (SRD), citrus bark cracking viroid SRD, as well as A-RNA. We have applied enhanced molecular dynamics techniques that construct an inverse free-energy surface (iFES) determined by collective variables that monitor base-pairing and backbone conformation. Each SRD RNA is flanked on each side by A-RNA, allowing comparison of the behavior of these motifs in the same molecule. The RNA iFESs have single peaks, indicating that the combined motifs should denature as a single cohesive unit, rather than by regional melting. Local root-mean-square deviation (RMSD) analysis and communication propensity (CProp, variance in distances between residue pairs) reveal distinct motif properties. Our analysis indicates that the standard SRD is more stable than the viroid SRD, which is more stable than A-RNA. Base pairs at SRD to A-RNA transitions have limited flexibility. Application of CProp reveals extraordinary stiffness of the SRD, allowing residues on opposite sides of the motif to sense each other's motions.

A Riboswitch-Driven Era of New Antibacterials

Riboswitches are structured non-coding RNAs found in the 5′ UTR of important genes for bacterial metabolism, virulence and survival. Upon the binding of specific ligands that can vary from simple ions to complex molecules such as nucleotides and tRNAs, riboswitches change their local and global mRNA conformations to affect downstream transcription or translation. Due to their dynamic nature and central regulatory role in bacterial metabolism, riboswitches have been exploited as novel RNA-based targets for the development of new generation antibacterials that can overcome drug-resistance problems. During recent years, several important riboswitch structures from many bacterial representatives, including several prominent human pathogens, have shown that riboswitches are ideal RNA targets for new compounds that can interfere with their structure and function, exhibiting much reduced resistance over time. Most interestingly, mainstream antibiotics that target the ribosome have been shown to effectively modulate the regulatory behavior and capacity of several riboswitches, both in vivo and in vitro, emphasizing the need for more in-depth studies and biological evaluation of new antibiotics. Herein, we summarize the currently known compounds that target several main riboswitches and discuss the role of mainstream antibiotics as modulators of T-box riboswitches, in the dawn of an era of novel inhibitors that target important bacterial regulatory RNAs.

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.

Structural modifications of berberine and their binding effects towards polymorphic deoxyribonucleic acid structures: A review

Berberine (BBR) is a plant derived quaternary benzylisoquinoline alkaloid, which has been widely used in traditional medicines for a long term. It possesses broad pharmacological effects and is widely applied in clinical. In recent years, the anti-tumor effects of BBR have attracted more and more attention of the researchers. The canonical right-handed double-stranded helical deoxyribonucleic acid (B-DNA) and its polymorphs occur under various environmental conditions and are involved in a plethora of genetic instability-related diseases especially tumor. BBR showed differential binding effects towards various polymorphic DNA structures. But its poor lipophilicity and fast metabolism limited its clinical utility. Structural modification of BBR is an effective approach to improve its DNA binding activity and bioavailability in vivo. A large number of studies dedicated to improving the binding affinities of BBR towards different DNA structures have been carried out and achieved tremendous advancements. In this article, the main achievements of BBR derivatives in polymorphic DNA structures binding researches in recent 20 years were reviewed. The structural modification strategy of BBR, the DNA binding effects of its derivatives, and the structure activity relationship (SAR) analysis have also been discussed.

The cotranscriptional folding landscape for two cyclic di-nucleotide-sensing riboswitches with highly homologous aptamer domains acting either as ON- or OFF-switches

Riboswitches are gene regulatory elements located in untranslated mRNA regions. They bind inducer molecules with high affinity and specificity. Cyclic-di-nucleotide-sensing riboswitches are major regulators of genes for the environment, membranes and motility (GEMM) of bacteria. Up to now, structural probing assays or crystal structures have provided insight into the interaction between cyclic-di-nucleotides and their corresponding riboswitches. ITC analysis, NMR analysis and computational modeling allowed us to gain a detailed understanding of the gene regulation mechanisms for the Cd1 (Clostridium difficile) and for the pilM (Geobacter metallireducens) riboswitches and their respective di-nucleotides c-di-GMP and c-GAMP. Binding capability showed a 25 nucleotide (nt) long window for pilM and a 61 nt window for Cd1. Within this window, binding affinities ranged from 35 μM to 0.25 μM spanning two orders of magnitude for Cd1 and pilM showing a strong dependence on competing riboswitch folds. Experimental results were incorporated into a Markov simulation to further our understanding of the transcriptional folding pathways of riboswitches. Our model showed the ability to predict riboswitch gene regulation and its dependence on transcription speed, pausing and ligand concentration.

Interrogating RNA-Small Molecule Interactions with Structure Probing and Artificial Intelligence-Augmented Molecular Simulations

While there is increasing interest in the study of RNA as a therapeutic target, efforts to understand RNA-ligand recognition at the molecular level lag far behind our understanding of protein-ligand recognition. This problem is complicated due to the more than 10 orders of magnitude in time scales involved in RNA dynamics and ligand binding events, making it not straightforward to design experiments or simulations. Here, we make use of artificial intelligence (AI)-augmented molecular dynamics simulations to directly observe ligand dissociation for cognate and synthetic ligands from a riboswitch system. The site-specific flexibility profiles from our simulations are compared with in vitro measurements of flexibility using selective 2' hydroxyl acylation analyzed by primer extension and mutational profiling (SHAPE-MaP). Our simulations reproduce known relative binding affinity profiles for the cognate and synthetic ligands, and pinpoint how both ligands make use of different aspects of riboswitch flexibility. On the basis of our dissociation trajectories, we also make and validate predictions of pairs of mutations for both the ligand systems that would show differing binding affinities. These mutations are distal to the binding site and could not have been predicted solely on the basis of structure. The methodology demonstrated here shows how molecular dynamics simulations with all-atom force-fields have now come of age in making predictions that complement existing experimental techniques and illuminate aspects of systems otherwise not trivial to understand.

Thermodynamic characterization of naturally occurring RNA pentaloops

RNA folding is hierarchical; therefore, predicting RNA secondary structure from sequence is an intermediate step in predicting tertiary structure. Secondary structure prediction is based on a nearest neighbor model using free energy minimization. To improve secondary structure prediction, all types of naturally occurring secondary structure motifs need to be thermodynamically characterized. However, not all secondary structure motifs are well characterized. Pentaloops, the second most abundant hairpin size, is one such uncharacterized motif. In fact, the current thermodynamic model used to predict the stability of pentaloops was derived from a small data set of pentaloops and from data for other hairpins of different sizes. Here, the most commonly occurring pentaloops were identified and optically melted. New experimental data for 22 pentaloop sequences were combined with previously published data for nine pentaloop sequences. Using linear regression, a pentaloop-specific model was derived. This new model is simpler and more accurate than the current model. The new experimental data and improved model can be incorporated into software that is used to predict RNA secondary structure from sequence.

Nearest neighbor rules for RNA helix folding thermodynamics: improved end effects

Nearest neighbor parameters for estimating the folding stability of RNA secondary structures are in widespread use. For helices, current parameters penalize terminal AU base pairs relative to terminal GC base pairs. We curated an expanded database of helix stabilities determined by optical melting experiments. Analysis of the updated database shows that terminal penalties depend on the sequence identity of the adjacent penultimate base pair. New nearest neighbor parameters that include this additional sequence dependence accurately predict the measured values of 271 helices in an updated database with a correlation coefficient of 0.982. This refined understanding of helix ends facilitates fitting terms for base pair stacks with GU pairs. Prior parameter sets treated 5′GGUC3′ paired to 3′CUGG5′ separately from other 5′GU3′/3′UG5′ stacks. The improved understanding of helix end stability, however, makes the separate treatment unnecessary. Introduction of the additional terms was tested with three optical melting experiments. The average absolute difference between measured and predicted free energy changes at 37°C for these three duplexes containing terminal adjacent AU and GU pairs improved from 1.38 to 0.27 kcal/mol. This confirms the need for the additional sequence dependence in the model.

Crowdsourced RNA design discovers diverse, reversible, efficient, self-contained molecular switches

Internet-based scientific communities promise a means to apply distributed, diverse human intelligence toward previously intractable scientific problems. However, current implementations have not allowed communities to propose experiments to test all emerging hypotheses at scale or to modify hypotheses in response to experiments. We report high-throughput methods for molecular characterization of nucleic acids that enable the large-scale video game–based crowdsourcing of RNA sensor design, followed by high-throughput functional characterization. Iterative design testing of thousands of crowdsourced RNA sensor designs produced near–thermodynamically optimal and reversible RNA switches that act as self-contained molecular sensors and couple five distinct small molecule inputs to three distinct protein binding and fluorogenic outputs. This work suggests a paradigm for widely distributed experimental bioscience.

RNA-ligand molecular docking: advances and challenges

With rapid advances in computer algorithms and hardware, fast and accurate virtual screening has led to a drastic acceleration in selecting potent small molecules as drug candidates. Computational modeling of RNA-small molecule interactions has become an indispensable tool for RNA-targeted drug discovery. The current models for RNA-ligand binding have mainly focused on the docking-and-scoring method. Accurate docking and scoring should tackle four crucial problems: (1) conformational flexibility of ligand, (2) conformational flexibility of RNA, (3) efficient sampling of binding sites and binding poses, and (4) accurate scoring of different binding modes. Moreover, compared with the problem of protein-ligand docking, predicting ligand binding to RNA, a negatively charged polymer, is further complicated by additional effects such as metal ion effects. Thermodynamic models based on physics-based and knowledge-based scoring functions have shown highly encouraging success in predicting ligand binding poses and binding affinities. Recently, kinetic models for ligand binding have further suggested that including dissociation kinetics (residence time) in ligand docking would result in improved performance in estimating in vivo drug efficacy. More recently, the rise of deep-learning approaches has led to new tools for predicting RNA-small molecule binding. In this review, we present an overview of the recently developed computational methods for RNA-ligand docking and their advantages and disadvantages.

High-throughput screening of cell-free riboswitches by fluorescence-activated droplet sorting

Cell-free systems that display complex functions without using living cells are emerging as new platforms to test our understanding of biological systems as well as for practical applications such as biosensors and biomanufacturing. Those that use cell-free protein synthesis (CFPS) systems to enable genetically programmed protein synthesis have relied on genetic regulatory components found or engineered in living cells. However, biological constraints such as cell permeability, metabolic stability, and toxicity of signaling molecules prevent development of cell-free devices using living cells even if cell-free systems are not subject to such constraints. Efforts to engineer regulatory components directly in CFPS systems thus far have been based on low-throughput experimental approaches, limiting the availability of basic components to build cell-free systems with diverse functions. Here, we report a high-throughput screening method to engineer cell-free riboswitches that respond to small molecules. Droplet-sorting of riboswitch variants in a CFPS system rapidly identified cell-free riboswitches that respond to compounds that are not amenable to bacterial screening methods. Finally, we used a histamine riboswitch to demonstrate chemical communication between cell-sized droplets.

Bacterial 2'-Deoxyguanosine Riboswitch Classes as Potential Targets for Antibiotics: A Structure and Dynamics Study

The spread of antibiotic-resistant bacteria represents a substantial health threat. Current antibiotics act on a few metabolic pathways, facilitating resistance. Consequently, novel regulatory inhibition mechanisms are necessary. Riboswitches represent promising targets for antibacterial drugs. Purine riboswitches are interesting, since they play essential roles in the genetic regulation of bacterial metabolism. Among these, class I (2′-dG-I) and class II (2′-dG-II) are two different 2′-deoxyguanosine (2′-dG) riboswitches involved in the control of deoxyguanosine metabolism. However, high affinity for nucleosides involves local or distal modifications around the ligand-binding pocket, depending on the class. Therefore, it is crucial to understand these riboswitches’ recognition mechanisms as antibiotic targets. In this work, we used a combination of computational biophysics approaches to investigate the structure, dynamics, and energy landscape of both 2′-dG classes bound to the nucleoside ligands, 2′-deoxyguanosine, and riboguanosine. Our results suggest that the stability and increased interactions in the three-way junction of 2′-dG riboswitches were associated with a higher nucleoside ligand affinity. Also, structural changes in the 2′-dG-II aptamers enable enhanced intramolecular communication. Overall, the 2′-dG-II riboswitch might be a promising drug design target due to its ability to recognize both cognate and noncognate ligands.

Strategies for Improving Small-Molecule Biosensors in Bacteria

In recent years, small-molecule biosensors have become increasingly important in synthetic biology and biochemistry, with numerous new applications continuing to be developed throughout the field. For many biosensors, however, their utility is hindered by poor functionality. Here, we review the known types of mechanisms of biosensors within bacterial cells, and the types of approaches for optimizing different biosensor functional parameters. Discussed approaches for improving biosensor functionality include methods of directly engineering biosensor genes, considerations for choosing genetic reporters, approaches for tuning gene expression, and strategies for incorporating additional genetic modules.

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

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

Riboswitch RSthiT as a molecular tool in Lactococcus lactis

Previous RNA sequencing has allowed the identification of 129 long 5′ untranslated regions (UTRs) in the Lactococcus lactis MG1363 transcriptome. These sequences potentially harbor cis-acting riboswitches. One of the identified extended 5′ UTRs is a putative thiamine pyrophosphate (TPP) riboswitch. It is located immediately upstream of the thiamine transporter gene thiT (llmg_0334). To confirm this assumption, the 5′-UTR sequence was placed upstream of the gene encoding the superfolder green fluorescent protein (sfGFP), sfgfp, allowing the examination of the expression of sfGFP in the presence or absence of thiamine in the medium. The results show that this sequence indeed represents a thiamine-responsive TPP riboswitch. This RNA-based genetic control device was used to successfully restore the mutant phenotype of an L. lactis strain lacking the major autolysin gene, acmA. The L. lactisthiT TPP riboswitch (RS_thiT) is a useful molecular genetic tool enabling the gradual downregulation of the expression of genes under its control by adjusting the thiamine concentration.

IMPORTANCE The capacity of microbes with biotechnological importance to adapt to and survive under quickly changing industrial conditions depends on their ability to adequately control gene expression. Riboswitches are important RNA-based elements involved in rapid and precise gene regulation. Here, we present the identification of a natural thiamine-responsive riboswitch of Lactococcus lactis, a bacterium used worldwide in the production of dairy products. We used it to restore a genetic defect in an L. lactis mutant and show that it is a valuable addition to the ever-expanding L. lactis genetic toolbox.

Directed Evolution in Drops: Molecular Aspects and Applications

The process of optimizing the properties of biological molecules is paramount for many industrial and medical applications. Directed evolution is a powerful technique for modifying and improving biomolecules such as proteins or nucleic acids (DNA or RNA). Mimicking the mechanism of natural evolution, one can enhance a desired property by applying a suitable selection pressure and sorting improved variants. Droplet-based microfluidic systems offer a high-throughput solution to this approach by helping to overcome the limiting screening steps and allowing the analysis of variants within increasingly complex libraries. Here, we review cases where successful evolution of biomolecules was achieved using droplet-based microfluidics, focusing on the molecular processes involved and the incorporation of microfluidics to the workflow. We highlight the advantages and limitations of these microfluidic systems compared to low-throughput methods and show how the integration of these systems into directed evolution workflows can open new avenues to discover or improve biomolecules according to user-defined conditions.

Regulation of the ThiM riboswitch is facilitated by the trapped structure formed during transcription of the wild-type sequence

The ThiM riboswitch from Escherichia coli is a typical mRNA device that modulates downstream gene expression by sensing TPP. The helix-based RNA folding theory is used to investigate its detailed regulatory behaviors in cells. This RNA molecule is transcriptionally trapped in a state with the unstructured SD sequence in the absence of TPP, which induces downstream gene expression. As a key step to turn on gene expression, formation of this trapped state (the genetic ON state) highly depends on the co-transcriptional folding of its wild-type sequence. Instead of stabilities of the genetic ON and OFF states, the transcription rate, pause, and ligand levels are combined to affect the ThiM riboswitch-mediated gene regulation, which is consistent with a kinetic control model.

Cell-free riboswitches

The emerging community of cell-free synthetic biology aspires to build complex biochemical and genetic systems with functions that mimic or even exceed those in living cells. To achieve such functions, cell-free systems must be able to sense and respond to the complex chemical signals within and outside the system. Cell-free riboswitches can detect chemical signals via RNA-ligand interaction and respond by regulating protein synthesis in cell-free protein synthesis systems. In this article, we review synthetic cell-free riboswitches that function in both prokaryotic and eukaryotic cell-free systems reported to date to provide a current perspective on the state of cell-free riboswitch technologies and their limitations.

The Role of RNA Secondary Structure in Regulation of Gene Expression in Bacteria

Due to the high exposition to changing environmental conditions, bacteria have developed many mechanisms enabling immediate adjustments of gene expression. In many cases, the required speed and plasticity of the response are provided by RNA-dependent regulatory mechanisms. This is possible due to the very high dynamics and flexibility of an RNA structure, which provide the necessary sensitivity and specificity for efficient sensing and transduction of environmental signals. In this review, we will discuss the current knowledge about known bacterial regulatory mechanisms which rely on RNA structure. To better understand the structure-driven modulation of gene expression, we describe the basic theory on RNA structure folding and dynamics. Next, we present examples of multiple mechanisms employed by RNA regulators in the control of bacterial transcription and translation.

Computational based design and tracking of synthetic variants of Porcine circovirus reveal relations between silent genomic information and viral fitness

Viral genomes not only code the protein content, but also include silent, overlapping codes which are important to the regulation of the viral life cycle and affect its evolution. Due to the high density of these codes, their non-modular nature and the complex intracellular processes they encode, the ability of current approaches to decipher them is very limited. We describe the first computational-experimental pipeline for studying the effects of viral silent and non-silent information on its fitness. The pipeline was implemented to study the Porcine Circovirus type 2 (PCV2), the shortest known eukaryotic virus, and includes the following steps: (1) Based on the analyses of 2100 variants of PCV, suspected silent codes were inferred. (2) Five hundred variants of the PCV2 were designed to include various ‘smart’ silent mutations. (3) Using state of the art synthetic biology approaches, the genomes of these five hundred variants were generated. (4) Competition experiments between the variants were performed in Porcine kidney-15 (PK15) cell-lines. (5) The variant titers were analyzed based on novel next-generation sequencing (NGS) experiments. (6) The features related to the titer of the variants were inferred and their analyses enabled detection of various novel silent functional sequence and structural motifs. Furthermore, we demonstrate that 50 of the silent variants exhibit higher fitness than the wildtype in the analyzed conditions.

DRAGoM: Classification and Quantification of Noncoding RNA in Metagenomic Data

Noncoding RNAs (ncRNAs) play important regulatory and functional roles in microorganisms, such as regulation of gene expression, signaling, protein synthesis, and RNA processing. Hence, their classification and quantification are central tasks toward the understanding of the function of the microbial community. However, the majority of the current metagenomic sequencing technologies generate short reads, which may contain only a partial secondary structure that complicates ncRNA homology detection. Meanwhile, de novo assembly of the metagenomic sequencing data remains challenging for complex communities. To tackle these challenges, we developed a novel algorithm called DRAGoM (Detection of RNA using Assembly Graph from Metagenomic data). DRAGoM first constructs a hybrid graph by merging an assembly string graph and an assembly de Bruijn graph. Then, it classifies paths in the hybrid graph and their constituent readsinto differentncRNA families based on both sequence and structural homology. Our benchmark experiments show that DRAGoMcan improve the performance and robustness over traditional approaches on the classification and quantification of a wide class of ncRNA families.

Thiamine pyrophosphate riboswitch regulation: a new possible mechanism involved in the action of nalidixic acid

Objectives

The development of novel antibiotic compounds requires riboswitches; in fact, riboswitches are RNA elements present in the 5′ untranslated region of bacterial mRNA and have a metabolite-binding aptamer domain and an expression platform regulating the expression of vital genes. In the present research, one riboswitch, namely thi-box riboswitch with distinct regulatory mechanisms, was studied. It recognizes Thiamine Pyrophosphates (TPP) regulating TPP-biosynthesis genes in Escherichia coli.

Methods

First, the compounds similar to riboswitch ligands were studied, and their binding with the riboswitch and nucleosides was investigated by molecular docking. Then, compounds containing high binding energy were chosen, and their minimum inhibitory concentration in E. coli was determined by the MIC test. Finally, the binding of compounds to nucleotides and RNA was investigated by measuring the absorbance spectrum through NanoDrop and circular dichroism (CD).

Results

In the thi-box riboswitch, nalidixic acid was found to have the best binding energy (−5.31 kJ/mol), and it inhibited E. coli growth at the minimum inhibitory concentration of 125 μg/mL, and it could bind to ribonucleosides and RNA in vitro.

Conclusions

One possible mechanism involved in the action of nalidixic acid in inhibiting the E. coli growth is to influence thi-box riboswitch.

Molecular crowding and RNA catalysis

RNA enzymes or ribozymes catalyze some of the most important reactions in biology and are thought to have played a central role in the origin and evolution of life on earth. Catalytic function in RNA has evolved in crowded cellular environments that are different from dilute solutions in which most in vitro assays are performed. The presence of molecules such as amino acids, polypeptides, alcohols, and sugars in the cell introduces forces that modify the kinetics and thermodynamics of ribozyme-catalyzed reactions. Synthetic molecules are routinely used in in vitro studies to better approximate the properties of biomolecules under in vivo conditions. This review discusses the various forces that operate within simulated crowded solutions in the context of RNA structure, folding, and catalysis. It also explores ideas about how crowding could have been beneficial to the evolution of functional RNAs and the development of primitive cellular systems in a prebiotic milieu.

RNA Drugs and RNA Targets for Small Molecules: Principles, Progress, and Challenges

RNA-based therapies, including RNA molecules as drugs and RNA-targeted small molecules, offer unique opportunities to expand the range of therapeutic targets. Various forms of RNAs may be used to selectively act on proteins, transcripts, and genes that cannot be targeted by conventional small molecules or proteins. Although development of RNA drugs faces unparalleled challenges, many strategies have been developed to improve RNA metabolic stability and intracellular delivery. A number of RNA drugs have been approved for medical use, including aptamers (e.g., pegaptanib) that mechanistically act on protein target and small interfering RNAs (e.g., patisiran and givosiran) and antisense oligonucleotides (e.g., inotersen and golodirsen) that directly interfere with RNA targets. Furthermore, guide RNAs are essential components of novel gene editing modalities, and mRNA therapeutics are under development for protein replacement therapy or vaccination, including those against unprecedented severe acute respiratory syndrome coronavirus pandemic. Moreover, functional RNAs or RNA motifs are highly structured to form binding pockets or clefts that are accessible by small molecules. Many natural, semisynthetic, or synthetic antibiotics (e.g., aminoglycosides, tetracyclines, macrolides, oxazolidinones, and phenicols) can directly bind to ribosomal RNAs to achieve the inhibition of bacterial infections. Therefore, there is growing interest in developing RNA-targeted small-molecule drugs amenable to oral administration, and some (e.g., risdiplam and branaplam) have entered clinical trials. Here, we review the pharmacology of novel RNA drugs and RNA-targeted small-molecule medications, with a focus on recent progresses and strategies. Challenges in the development of novel druggable RNA entities and identification of viable RNA targets and selective small-molecule binders are discussed. SIGNIFICANCE STATEMENT: With the understanding of RNA functions and critical roles in diseases, as well as the development of RNA-related technologies, there is growing interest in developing novel RNA-based therapeutics. This comprehensive review presents pharmacology of both RNA drugs and RNA-targeted small-molecule medications, focusing on novel mechanisms of action, the most recent progress, and existing challenges.

Riboflow: Using Deep Learning to Classify Riboswitches With ∼99% Accuracy

Riboswitches are cis-regulatory genetic elements that use an aptamer to control gene expression. Specificity to cognate ligand and diversity of such ligands have expanded the functional repertoire of riboswitches to mediate mounting apt responses to sudden metabolic demands and signal changes in environmental conditions. Given their critical role in microbial life, riboswitch characterisation remains a challenging computational problem. Here we have addressed the issue with advanced deep learning frameworks, namely convolutional neural networks (CNN), and bidirectional recurrent neural networks (RNN) with Long Short-Term Memory (LSTM). Using a comprehensive dataset of 32 ligand classes and a stratified train-validate-test approach, we demonstrated the accurate performance of both the deep learning models (CNN and RNN) relative to conventional hyperparameter-optimized machine learning classifiers on all key performance metrics, including the ROC curve analysis. In particular, the bidirectional LSTM RNN emerged as the best-performing learning method for identifying the ligand-specificity of riboswitches with an accuracy >0.99 and macro-averaged F-score of 0.96. An additional attraction is that the deep learning models do not require prior feature engineering. A dynamic update functionality is built into the models to factor for the constant discovery of new riboswitches, and extend the predictive modeling to new classes. Our work would enable the design of genetic circuits with custom-tuned riboswitch aptamers that would effect precise translational control in synthetic biology. The associated software is available as an open-source Python package and standalone resource for use in genome annotation, synthetic biology, and biotechnology workflows.