Molecular sensing by the aptamer domain of the FMN riboswitch: a general model for ligand binding by conformational selection

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.

Advances in machine-learning approaches to RNA-targeted drug design

RNA molecules play multifaceted functional and regulatory roles within cells and have garnered significant attention in recent years as promising therapeutic targets. With remarkable successes achieved by artificial intelligence (AI) in different fields such as computer vision and natural language processing, there is a growing imperative to harness AI's potential in computer-aided drug design (CADD) to discover novel drug compounds that target RNA. Although machine-learning (ML) approaches have been widely adopted in the discovery of small molecules targeting proteins, the application of ML approaches to model interactions between RNA and small molecule is still in its infancy. Compared to protein-targeted drug discovery, the major challenges in ML-based RNA-targeted drug discovery stem from the scarcity of available data resources. With the growing interest and the development of curated databases focusing on interactions between RNA and small molecule, the field anticipates a rapid growth and the opening of a new avenue for disease treatment. In this review, we aim to provide an overview of recent advancements in computationally modeling RNA-small molecule interactions within the context of RNA-targeted drug discovery, with a particular emphasis on methodologies employing ML techniques.

The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

RNA-Small-Molecule Interaction: Challenging the "Undruggable" Tag

RNA targeting, specifically with small molecules, is a relatively new and rapidly emerging avenue with the promise to expand the target space in the drug discovery field. From being "disregarded" as an "undruggable" messenger molecule to FDA approval of an RNA-targeting small-molecule drug Risdiplam, a radical change in perspective toward RNA has been observed in the past decade. RNAs serve important regulatory functions beyond canonical protein synthesis, and their dysregulation has been reported in many diseases. A deeper understanding of RNA biology reveals that RNA molecules can adopt a variety of structures, carrying defined binding pockets that can accommodate small-molecule drugs. Due to its functional diversity and structural complexity, RNA can be perceived as a prospective target for therapeutic intervention. This perspective highlights the proof of concept of RNA-small-molecule interactions, exemplified by targeting of various transcripts with functional modulators. The advent of RNA-oriented knowledge would help expedite drug discovery.

Flipping the script: Understanding riboswitches from an alternative perspective

Riboswitches are broadly distributed regulatory elements most frequently found in the 5'-leader sequence of bacterial mRNAs that regulate gene expression in response to the binding of a small molecule effector. The occupancy status of the ligand-binding aptamer domain manipulates downstream information in the message that instructs the expression machinery. Currently, there are over 55 validated riboswitch classes, where each class is defined based on the identity of the ligand it binds and/or sequence and structure conservation patterns within the aptamer domain. This classification reflects an "aptamer-centric" perspective that dominates our understanding of riboswitches. In this review, we propose a conceptual framework that groups riboswitches based on the mechanism by which RNA manipulates information directly instructing the expression machinery. This scheme does not replace the established aptamer domain-based classification of riboswitches but rather serves to facilitate hypothesis-driven investigation of riboswitch regulatory mechanisms. Based on current bioinformatic, structural, and biochemical studies of a broad spectrum of riboswitches, we propose three major mechanistic groups: (1) "direct occlusion", (2) "interdomain docking", and (3) "strand exchange". We discuss the defining features of each group, present representative examples of riboswitches from each group, and illustrate how these RNAs couple small molecule binding to gene regulation. While mechanistic studies of the occlusion and docking groups have yielded compelling models for how these riboswitches function, much less is known about strand exchange processes. To conclude, we outline the limitations of our mechanism-based conceptual framework and discuss how critical information within riboswitch expression platforms can inform gene regulation.

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.

The Effect of Pseudoknot Base Pairing on Cotranscriptional Structural Switching of the Fluoride Riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of E. coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37 °C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65 °C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

Role of riboflavin biosynthesis gene duplication and transporter in Aeromonas salmonicida virulence in marine teleost fish

Active flavins derived from riboflavin (vitamin B₂) are essential for life. Bacteria biosynthesize riboflavin or scavenge it through uptake systems, and both mechanisms may be present. Because of riboflavin's critical importance, the redundancy of riboflavin biosynthetic pathway (RBP) genes might be present. Aeromonas salmonicida, the aetiological agent of furunculosis, is a pathogen of freshwater and marine fish, and its riboflavin pathways have not been studied. This study characterized the A. salmonicida riboflavin provision pathways. Homology search and transcriptional orchestration analysis showed that A. salmonicida has a main riboflavin biosynthetic operon that includes ribD, ribE1, ribBA, and ribH genes. Outside the main operon, putative duplicated genes ribA, ribB and ribE, and a ribN riboflavin importer encoding gene, were found. Monocistronic mRNA ribA, ribB and ribE2 encode for their corresponding functional riboflavin biosynthetic enzyme. While the product of ribBA conserved the RibB function, it lacked the RibA function. Likewise, ribN encodes a functional riboflavin importer. Transcriptomics analysis indicated that external riboflavin affected the expression of a relatively small number of genes, including a few involved in iron metabolism. ribB was downregulated in response to external riboflavin, suggesting negative feedback. Deletion of ribA, ribB and ribE1 showed that these genes are required for A. salmonicida riboflavin biosynthesis and virulence in Atlantic lumpfish (Cyclopterus lumpus). A. salmonicida riboflavin auxotrophic attenuated mutants conferred low protection to lumpfish against virulent A. salmonicida. Overall, A. salmonicida has multiple riboflavin endowment forms, and duplicated riboflavin provision genes are critical for A. salmonicida infection.

Targeting FMN, TPP, SAM-I, and glmS Riboswitches with Chimeric Antisense Oligonucleotides for Completely Rational Antibacterial Drug Development

Antimicrobial drug resistance has emerged as a significant challenge in contemporary medicine due to the proliferation of numerous bacterial strains resistant to all existing antibiotics. Meanwhile, riboswitches have emerged as promising targets for discovering antibacterial drugs. Riboswitches are regulatory elements in certain bacterial mRNAs that can bind to specific molecules and control gene expression via transcriptional termination, prevention of translation, or mRNA destabilization. By targeting riboswitches, we aim to develop innovative strategies to combat antibiotic-resistant bacteria and enhance the efficacy of antibacterial treatments. This convergence of challenges and opportunities underscores the ongoing quest to revolutionize medical approaches against evolving bacterial threats. For the first time, this innovative review describes the rational design and applications of chimeric antisense oligonucleotides as antibacterial agents targeting four riboswitches selected based on genome-wide bioinformatic analyses. The antisense oligonucleotides are coupled with the cell-penetrating oligopeptide pVEC, which penetrates Gram-positive and Gram-negative bacteria and specifically targets glmS, FMN, TPP, and SAM-I riboswitches in Staphylococcus aureus, Listeria monocytogenes, and Escherichia coli. The average antibiotic dosage of antisense oligonucleotides that inhibits 80% of bacterial growth is around 700 nM (4.5 μg/mL). Antisense oligonucleotides do not exhibit toxicity in human cell lines at this concentration. The results demonstrate that these riboswitches are suitable targets for antibacterial drug development using antisense oligonucleotide technology. The approach is fully rational because selecting suitable riboswitch targets and designing ASOs that target them are based on predefined criteria. The approach can be used to develop narrow or broad-spectrum antibiotics against multidrug-resistant bacterial strains for a short time. The approach is easily adaptive to new resistance using targeting NGS technology.

Capturing heterogeneous conformers of cobalamin riboswitch by cryo-EM

RNA conformational heterogeneity often hampers its high-resolution structure determination, especially for large and flexible RNAs devoid of stabilizing proteins or ligands. The adenosylcobalamin riboswitch exhibits heterogeneous conformations under 1 mM Mg2+ concentration and ligand binding reduces conformational flexibility. Among all conformers, we determined one apo (5.3 Å) and four holo cryo-electron microscopy structures (overall 3.0-3.5 Å, binding pocket 2.9-3.2 Å). The holo dimers exhibit global motions of helical twisting and bending around the dimer interface. A backbone comparison of the apo and holo states reveals a large structural difference in the P6 extension position. The central strand of the binding pocket, junction 6/3, changes from an 'S'- to a 'U'-shaped conformation to accommodate ligand. Furthermore, the binding pocket can partially form under 1 mM Mg2+ and fully form under 10 mM Mg2+ within the bound-like structure in the absence of ligand. Our results not only demonstrate the stabilizing ligand-induced conformational changes in and around the binding pocket but may also provide further insight into the role of the P6 extension in ligand binding and selectivity.

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.

Direct evolution of riboflavin kinase significantly enhance flavin mononucleotide synthesis by design and optimization of flavin mononucleotide riboswitch

Flavin mononucleotide (FMN) is the active form of riboflavin. It has a wide range of application scenarios in the pharmaceutical and food additives. However, there are limitations in selecting generic high-throughput screening platforms that improve the properties of enzymes. First, the biosensor in response to FMN concentration was constructed using the FMN riboswitch and confirmed the function of this sensor. Next, the FMN binding site of the sensor was saturated with a mutation that increased its fluorescence range by approximately 127%. Then, the biosensor and the base editing system based on T7RNAP were combined to construct a platform for rapid mutation and screening of riboflavin kinase gene ribC mutants. The mutants screened using this platform increased the yield of FMN by 8-fold. These results indicate that the high-throughput screening platform can rapidly and effectively improve the activity of target enzymes, and provide a new route for screening industrial enzymes.

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.

Targeting RNA structures with small molecules

RNA adopts 3D structures that confer varied functional roles in human biology and dysfunction in disease. Approaches to therapeutically target RNA structures with small molecules are being actively pursued, aided by key advances in the field including the development of computational tools that predict evolutionarily conserved RNA structures, as well as strategies that expand mode of action and facilitate interactions with cellular machinery. Existing RNA-targeted small molecules use a range of mechanisms including directing splicing - by acting as molecular glues with cellular proteins (such as branaplam and the FDA-approved risdiplam), inhibition of translation of undruggable proteins and deactivation of functional structures in noncoding RNAs. Here, we describe strategies to identify, validate and optimize small molecules that target the functional transcriptome, laying out a roadmap to advance these agents into the next decade.

Long-read single-molecule RNA structure sequencing using nanopore

RNA molecules can form secondary and tertiary structures that can regulate their localization and function. Using enzymatic or chemical probing together with high-throughput sequencing, secondary structure can be mapped across the entire transcriptome. However, a limiting factor is that only population averages can be obtained since each read is an independent measurement. Although long-read sequencing has recently been used to determine RNA structure, these methods still used aggregate signals across the strands to detect structure. Averaging across the population also means that only limited information about structural heterogeneity across molecules or dependencies within each molecule can be obtained. Here, we present Single-Molecule Structure sequencing (SMS-seq) that combines structural probing with native RNA sequencing to provide non-amplified, structural profiles of individual molecules with novel analysis methods. Our new approach using mutual information enabled single molecule structural interrogation. Each RNA is probed at numerous bases enabling the discovery of dependencies and heterogeneity of structural features. We also show that SMS-seq can capture tertiary interactions, dynamics of riboswitch ligand binding, and mRNA structural features.

Computational Modeling of RNA Aptamers: Structure Prediction of the Apo State

RNA aptamers are single-stranded oligonucleotides that bind to specific molecular targets with high affinity and specificity. To design aptamers for new applications, it is critical to understand the ligand binding mechanism in terms of the structure and dynamics of the ligand-bound and apo states. The problem is that most of the NMR or X-ray crystal structures available for RNA aptamers are for ligand-bound states. Available apo state structures, mostly characterized by crystallization under nonphysiological conditions or probed by low resolution techniques, might fail to represent the diverse structural variations of the apo state in solution. Here, we develop an approach to obtain a representative ensemble of apo structures that are based on in silico RNA 3D structure prediction and in vitro experiments that characterize base stacking. Using the neomycin-B aptamer as a case study, an ensemble of structures for the aptamer in the apo (unbound) state are validated and then used to investigate the ligand-binding mechanism for the aptamer in complex with neomycin-B.

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.

SILCS-RNA: Toward a Structure-Based Drug Design Approach for Targeting RNAs with Small Molecules

RNA molecules can act as potential drug targets in different diseases, as their dysregulated expression or misfolding can alter various cellular processes. Noncoding RNAs account for ∼70% of the human genome, and these molecules can have complex tertiary structures that present a great opportunity for targeting by small molecules. In the present study, the site identification by ligand competitive saturation (SILCS) computational approach is extended to target RNA, termed SILCS-RNA. Extensions to the method include an enhanced oscillating excess chemical potential protocol for the grand canonical Monte Carlo calculations and individual simulations of the neutral and charged solutes from which the SILCS functional group affinity maps (FragMaps) are calculated for subsequent binding site identification and docking calculations. The method is developed and evaluated against seven RNA targets and their reported small molecule ligands. SILCS-RNA provides a detailed characterization of the functional group affinity pattern in the small molecule binding sites, recapitulating the types of functional groups present in the ligands. The developed method is also shown to be useful for identification of new potential binding sites and identifying ligand moieties that contribute to binding, granular information that can facilitate ligand design. However, limitations in the method are evident including the ability to map the regions of binding sites occupied by ligand phosphate moieties and to fully account for the wide range of conformational heterogeneity in RNA associated with binding of different small molecules, emphasizing inherent challenges associated with applying computer-aided drug design methods to RNA. While limitations are present, the current study indicates how the SILCS-RNA approach may enhance drug discovery efforts targeting RNAs with small molecules.

Affinity-Based Profiling of the Flavin Mononucleotide Riboswitch

Riboswitches are structural RNA elements that control gene expression. These naturally occurring RNA sensors are of continued interest as antibiotic targets, molecular sensors, and functional elements of synthetic circuits. Here, we describe affinity-based profiling of the flavin mononucleotide (FMN) riboswitch to characterize ligand binding and structural folding. We designed and synthesized photoreactive ligands and used them for photoaffinity labeling. We showed selective labeling of the FMN riboswitch and used this covalent interaction to quantitatively measure ligand binding, which we demonstrate with the naturally occurring antibiotic roseoflavin. We measured conditional riboswitch folding as a function of temperature and cation concentration. Furthermore, combining photoaffinity labeling with reverse transcription revealed ligand binding sites within the aptamer domain with single-nucleotide resolution. The photoaffinity probe was applied to cellular extracts of Bacillus subtilis to demonstrate conditional folding of the endogenous low-abundant ribD FMN riboswitch in biologically derived samples using quantitative PCR. Lastly, binding of the riboswitch-targeting antibiotic roseoflavin to the FMN riboswitch was measured in live bacteria using the photoaffinity probe.

Thoughts on how to think (and talk) about RNA structure

Recent events have pushed RNA research into the spotlight. Continued discoveries of RNA with unexpected diverse functions in healthy and diseased cells, such as the role of RNA as both the source and countermeasure to a severe acute respiratory syndrome coronavirus 2 infection, are igniting a new passion for understanding this functionally and structurally versatile molecule. Although RNA structure is key to function, many foundational characteristics of RNA structure are misunderstood, and the default state of RNA is often thought of and depicted as a single floppy strand. The purpose of this perspective is to help adjust mental models, equipping the community to better use the fundamental aspects of RNA structural information in new mechanistic models, enhance experimental design to test these models, and refine data interpretation. We discuss six core observations focused on the inherent nature of RNA structure and how to incorporate these characteristics to better understand RNA structure. We also offer some ideas for future efforts to make validated RNA structural information available and readily used by all researchers.

Subsite Ligand Recognition and Cooperativity in the TPP Riboswitch: Implications for Fragment-Linking in RNA Ligand Discovery

RNA molecules can show high levels of cooperativity in their global folding and interactions with divalent ions. However, cooperativity at individual ligand-RNA interaction sites remains poorly understood. Here, we investigated the binding of thiamine and methylene diphosphonic acid (MDP, a soluble structural analogue of pyrophosphate) to the thiamine pyrophosphate riboswitch. These ligands each bind weakly at proximal subsites, with 10 μM and 1 mM affinities, respectively. The affinity of MDP moderately improves when thiamine or thiamine-like fragments are pre-bound to the RNA. Covalent linking of thiamine and MDP substantially increases riboswitch binding to a notable high affinity of 20 nM. Crystal structures and single-molecule correlated chemical probing revealed favorable induced fit effects upon binding of individual ligands and, unexpectedly, a substantial thermodynamically unfavorable RNA structural rearrangement upon binding of the linked thiamine-MDP ligand. Thus, linking of two ligands of modest affinity, accompanied by an unfavorable structural rearrangement, still yields a potent linked RNA-binding compound. Since complex ligands often bind riboswitches and other RNAs at proximal subsites, principles derived from this work inform and support fragment-linking strategies for identifying small molecules that interact with RNA specifically and with high affinity.

Deciphering the mechanism of action of 5FDQD and the design of new neutral analogues for the FMN riboswitch: a well-tempered metadynamics simulation study

The FMN riboswitch is a novel drug target for the design of new antibiotics, and efforts have been made to design new charged and uncharged ligands. Uncharged ligands have shown advantages of not requiring any transporter for intracellular transport or proteins for their phosphorylation. 5FDQD (5-(3-(4-fluorophenyl)butyl)-7,8-dimethylpyrido(3,4-b)quinoxaline-1,3(2H,5H)-dione) is a recently reported neutral ligand for the FMN riboswitch active against Clostridium difficile infection in mice. However, the crystal structure of the 5FDQD bound FMN riboswitch is not available, and the mechanism of ligand binding and triggering the function of the riboswitch is not well understood. We have examined 5FDQD for its binding affinity with the FMN riboswitch using the well-tempered metadynamics (WT-MtD) simulation technique. The crystal structure of the FMN riboswitch shows that the FMN interacts with the J4/5 region through the phosphate group with G62; however, the uncharged ligands take advantage of π-π stacking interactions with the same residue of the riboswitch observed from the WT-MtD simulation results. The simulation results show that the presence of fluorine on the phenyl ring in 5FDQD is important to enhance the binding affinity of the neutral ligands with the FMN riboswitch. The WT-MtD results showed that the 1,2-difluoro substitution on the phenyl ring in 5FDQD (FMN-difluoro2) and the 1,3 positions in the phenyl ring (FMN-difluoro1) showed weaker binding energy with the FMN riboswitch compared to 5FDQD. The substitution of another fluorine atom at the 5-position of the phenyl ring (FMN-trifluoro) showed a comparable binding affinity (∼-31.4 kcal mol^-1) to 5FDQD. Electron-donating substitution on the phenyl ring such as the amino group also lowered the binding affinity (-28.8 kcal mol^-1) with the riboswitch compared to 5FDQD. The computed results suggest that the position and nature of substitution in the phenyl ring of the uncharged ligands affect the overall binding and such a delicate balance is important to achieve superior binding affinity with the FMN riboswitch.

RLDOCK method for predicting RNA-small molecule binding modes

RNA molecules play critical roles in cellular functions at the level of gene expression and regulation. The intricate 3D structures and the functional roles of RNAs make RNA molecules ideal targets for therapeutic drugs. The rational design of RNA-targeted drug requires accurate modeling of RNA-ligand interactions. Recently a new computational tool, RLDOCK, was developed to predict ligand binding sites and binding poses. Using an iterative multiscale sampling and search algorithm and a energy-based evaluation of ligand poses, the method enables efficient and accurate predictions for RNA-ligand interactions. Here we present a detailed illustration of the computational procedure for the practical implementation of the RLDOCK method. Using Flavin mononucleotide (FMN) docking to F. nucleatum FMN riboswitch as an example, we illustrate the computational protocol for RLDOCK-based prediction of RNA- ligand interactions. The RLDOCK software is freely accessible at http://https://github.com/Vfold-RNA/RLDOCK.

The potential versatility of RNA catalysis

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

Quantum Mechanics Helps Uncover Atypical Recognition Features in the Flavin Mononucleotide Riboswitch

Estimating the binding energies of small molecules to RNA could help uncover their molecular recognition characteristics and be used to rationally design RNA-targeting chemical probes. Here, we leveraged the ability of the fragment molecular orbital (FMO) method to provide detailed pairwise energetic information to examine the interactions between the aptamer domain of the flavin mononucleotide (FMN)-responsive riboswitch and small-molecule ligands. After developing an efficient protocol for executing high-level FMO calculations on RNA-ligand complexes, we applied our protocol to nine FMN-aptamer-ligand complexes. We then used the results to identify "hot-spots" within the aptamer and decomposed pairwise interactions between the hot-spot residues and the ligands. Interestingly, we found that several of these hot-spot residues interact with the ligands via atypical CH···O hydrogen bonds and anion-π contacts, as well as (face-to-edge) T-shaped π-π interactions. We envision that our results should pave the way for the wider and more prominent use of FMO calculations to study structure-energy relationships in diverse RNA-ligand systems, which in turn may provide a basis for dissecting the molecular recognition characteristics of RNAs.

DrugPred_RNA-A Tool for Structure-Based Druggability Predictions for RNA Binding Sites

RNA is an emerging target for drug discovery. However, like for proteins, not all RNA binding sites are equally suited to be addressed with conventional drug-like ligands. To this end, we have developed the structure-based druggability predictor DrugPred_RNA to identify druggable RNA binding sites. Due to the paucity of annotated RNA binding sites, the predictor was trained on protein pockets, albeit using only descriptors that can be calculated for both RNA and protein binding sites. DrugPred_RNA performed well in discriminating druggable from less druggable binding sites for the protein set and delivered predictions for selected RNA binding sites that agreed with manual assignment. In addition, most drug-like ligands contained in an RNA test set were found in pockets predicted to be druggable, further adding confidence to the performance of DrugPred_RNA. The method is robust against conformational and sequence changes in the binding sites and can contribute to direct drug discovery efforts for RNA targets.

The Perturbed Free-Energy Landscape: Linking Ligand Binding to Biomolecular Folding

The effects of ligand binding on biomolecular conformation are crucial in drug design, enzyme mechanisms, the regulation of gene expression, and other biological processes. Descriptive models such as "lock and key", "induced fit", and "conformation selection" are common ways to interpret such interactions. Another historical model, linked equilibria, proposes that the free-energy landscape (FEL) is perturbed by the addition of ligand binding energy for the bound population of biomolecules. This principle leads to a unified, quantitative theory of ligand-induced conformation change, building upon the FEL concept. We call the map of binding free energy over biomolecular conformational space the "binding affinity landscape" (BAL). The perturbed FEL predicts/explains ligand-induced conformational changes conforming to all common descriptive models. We review recent experimental and computational studies that exemplify the perturbed FEL, with emphasis on RNA. This way of understanding ligand-induced conformation dynamics motivates new experimental and theoretical approaches to ligand design, structural biology and systems biology.

AnnapuRNA: A scoring function for predicting RNA-small molecule binding poses

RNA is considered as an attractive target for new small molecule drugs. Designing active compounds can be facilitated by computational modeling. Most of the available tools developed for these prediction purposes, such as molecular docking or scoring functions, are parametrized for protein targets. The performance of these methods, when applied to RNA-ligand systems, is insufficient. To overcome these problems, we developed AnnapuRNA, a new knowledge-based scoring function designed to evaluate RNA-ligand complex structures, generated by any computational docking method. We also evaluated three main factors that may influence the structure prediction, i.e., the starting conformer of a ligand, the docking program, and the scoring function used. We applied the AnnapuRNA method for a post-hoc study of the recently published structures of the FMN riboswitch. Software is available at https://github.com/filipspl/AnnapuRNA.

Drug development is a lengthy and complicated process, which requires costly experiments on a very large number of chemical compounds. The identification of chemical molecules with desired properties can be facilitated by computational methods. Several methods were developed for computer-aided design of drugs that target protein molecules. However, recently the ribonucleic acid (RNA) emerged as an attractive target for the development of new drugs. Unfortunately, the portfolio of the computer methods that can be applied to study RNA and its interactions with small chemical molecules is very limited. This situation motivated us to develop a new computational method, with which to predict RNA-small molecule interactions. To this end, we collected the information on the statistics of interactions in experimentally determined structures of complexes formed by RNA with small molecules. We then used the statistical data to train machine learning methods aiming to distinguish between RNA-ligand interactions observed experimentally and other interactions that can be observed in theoretical analyses, but are not observed in nature. The resulting method called AnnapuRNA is superior to other similar tools and can be used to predict preferred ligands of RNA molecules and how RNA and small molecules interact with each other.

Riboswitches as Drug Targets for Antibiotics

Riboswitches reside in the untranslated region of RNA and regulate genes involved in the biosynthesis of essential metabolites through binding of small molecules. Since their discovery at the beginning of this century, riboswitches have been regarded as potential antibacterial targets. Using fragment screening, high-throughput screening and rational ligand design guided by X-ray crystallography, lead compounds against various riboswitches have been identified. Here, we review the current status and suitability of the thiamine pyrophosphate (TPP), flavin mononucleotide (FMN), glmS, guanine, and other riboswitches as antibacterial targets and discuss them in a biological context. Further, we highlight challenges in riboswitch drug discovery and emphasis the need to develop riboswitch specific high-throughput screening methods.

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

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

FMN riboswitch aptamer symmetry facilitates conformational switching through mutually exclusive coaxial stacking configurations

Knowledge of both apo and holo states of riboswitches aid in elucidating the various mechanisms of ligand-induced conformational “switching” that underpin their gene-regulating capabilities. Previous structural studies on the flavin mononucleotide (FMN)-binding aptamer of the FMN riboswitch, however, have revealed minimal conformational changes associated with ligand binding that do not adequately explain the basis for the switching behavior. We have determined a 2.7-Å resolution crystal structure of the ligand-free FMN riboswitch aptamer that is distinct from previously reported structures, particularly in the conformation and orientation of the P1 and P4 helices. The nearly symmetrical tertiary structure provides a mechanism by which one of two pairs of adjacent helices (P3/P4 or P1/P6) undergo collinear stacking in a mutually exclusive manner, in the absence or presence of ligand, respectively. Comparison of these structures suggests the stem-loop that includes P4 and L4 is important for maintaining a global conformational state that, in the absence of ligand, disfavors formation of the P1 regulatory helix. Together, these results provide further insight to the structural basis for conformational switching of the FMN riboswitch.

Structural insights into the interactions of flavin mononucleotide (FMN) and riboflavin with FMN riboswitch: a molecular dynamics simulation study

Antibiotics resistance is becoming a serious problem associated with fatalities and suffering patients. New antibiotics that can target the broader spectrum of cellular processes are warranted. One of the recent approaches in this regard is to target the special type of RNA riboswitches in bacteria. In this report, we have explored the mechanistic pathways of ligand-dependent conformational changes of flavin mononucleotide (FMN) riboswitch using molecular dynamics (MD) simulation studies. Cognate ligands FMN and riboflavin (RBF) have shown very different behavior with FMN riboswitch in terms of their role in the gene regulation process. These two ligands have similar scaffold, except the terminal phosphate group in FMN ligand. The MD simulations reveal that the binding of FMN ligand with the riboswitch does not lead to global folding of structure, rather lead to local changes in riboswitch structure. The binding free energy calculated with molecular mechanics Poisson-Boltzmann surface area method suggests the stronger binding of FMN than RBF to the riboswitch and electrostatic energy contributes chiefly to stabilize the complex. Further, the hydrogen bonding analysis identified the key binding site residues G11, G32, G62 of the riboswitch with FMN and RBF. The critical role of the phosphate group in the FMN ligand for binding with the active site of a riboswitch is also borne out in this study. These results unravel the importance of functional groups in natural ligands on designing newer ligands for FMN riboswitch as new antibiotics in the future.Communicated by Ramaswamy H. Sarma.

Polyacetate and Polycarbonate RNA: Acylating Reagents and Properties

Acylation of RNA at 2'-OH groups is widely applied in mapping RNA structure and recently for controlling RNA function. Reactions are described that install the smallest 2-carbon acyl groups on RNA-namely, 2'-O-acetyl and 2'-O-carbonate groups. Hybridization and thermal melting experiments are performed to assess the effects of the acyl groups on duplex formation. Both reagents can be employed at lower concentrations to map RNA secondary structure by reverse transcriptase primer extension (SHAPE) methods.

Regulation of Bacterial Gene Expression by Transcription Attenuation

A wide variety of mechanisms that control gene expression in bacteria are based on conditional transcription termination. Generally, in these mechanisms, a transcription terminator is located between a promoter and a downstream gene(s), and the efficiency of the terminator is controlled by a regulatory effector that can be a metabolite, protein, or RNA. The most common type of regulation involving conditional termination is transcription attenuation, in which the primary regulatory target is an essential element of a single terminator. The terminator can be either intrinsic or Rho dependent, with each presenting unique regulatory targets. Transcription attenuation mechanisms can be divided into five classes based primarily on the manner in which transcription termination is rendered conditional. This review summarizes each class of control mechanisms from a historical perspective, describes important examples in a physiological context and the current state of knowledge, highlights major advances, and discusses expectations of future discoveries.

Using Rosetta for RNA homology modeling

The three-dimensional structures of RNA molecules provide rich and often critical information for understanding their functions, including how they recognize small molecule and protein partners. Computational modeling of RNA 3D structure is becoming increasingly accurate, particularly with the availability of growing numbers of template structures already solved experimentally and the development of sequence alignment and 3D modeling tools to take advantage of this database. For several recent "RNA puzzle" blind modeling challenges, we have successfully identified useful template structures and achieved accurate structure predictions through homology modeling tools developed in the Rosetta software suite. We describe our semi-automated methodology here and walk through two illustrative examples: an adenine riboswitch aptamer, modeled from a template guanine riboswitch structure, and a SAM I/IV riboswitch aptamer, modeled from a template SAM I riboswitch structure.

Evidence for ligandable sites in structured RNA throughout the Protein Data Bank

RNA has attracted considerable attention as a target for small molecules. However, methods to identify, study, and characterize suitable RNA targets have lagged behind strategies for protein targets. One approach that has received considerable attention for protein targets has been to utilize computational analysis to investigate ligandable "pockets" on proteins that are amenable to small molecule binding. These studies have shown that selected physical properties of pockets are important parameters that govern the ability of a structure to bind to small molecules. This work describes a similar analysis to study pockets on all RNAs in the Protein Data Bank (PDB). Using parameters such as buriedness, hydrophobicity, volume, and other properties, the set of all RNAs is analyzed and compared to all proteins. Considerable overlap is observed between the properties of pockets on RNAs and proteins. Thus, many RNAs are capable of populating conformations with pockets that are likely suitable for small molecule binding. Further, principal moment of inertia (PMI) calculations reveal that liganded RNAs exist in diverse structural space, much of which overlaps with protein structural space. Taken together, these results suggest that complex folded RNAs adopt unique structures with pockets that may represent viable opportunities for small molecule targeting.

An NMR View of Protein Dynamics in Health and Disease

Biological molecules are often highly dynamic, and this flexibility can be critical for function. The large range of sampled timescales and the fact that many of the conformers that are continually explored are only transiently formed and sparsely populated challenge current biophysical approaches. Solution nuclear magnetic resonance (NMR) spectroscopy has emerged as a powerful method for characterizing biomolecular dynamics in detail, even in cases where excursions involve short-lived states. Here, we briefly review a number of NMR experiments for studies of biomolecular dynamics on the microsecond-to-second timescale and focus on applications to protein and nucleic acid systems that clearly illustrate the functional relevance of motion in both health and disease.

Duplexed aptamers: history, design, theory, and application to biosensing

Nucleic acid aptamers are single stranded DNA or RNA sequences that specifically bind a cognate ligand. In addition to their widespread use as stand-alone affinity binding reagents in analytical chemistry, aptamers have been engineered into a variety of ligand-specific biosensors, termed aptasensors. One of the most common aptasensor formats is the duplexed aptamer (DA). As defined herein, DAs are aptasensors containing two nucleic acid elements coupled via Watson-Crick base pairing: (i) an aptamer sequence, which serves as a ligand-specific receptor, and (ii) an aptamer-complementary element (ACE), such as a short DNA oligonucleotide, which is designed to hybridize to the aptamer. The ACE competes with ligand binding, such that DAs generate a signal upon ligand-dependent ACE-aptamer dehybridization. DAs possess intrinsic advantages over other aptasensor designs. For example, DA biosensing designs generalize across DNA and RNA aptamers, DAs are compatible with many readout methods, and DAs are inherently tunable on the basis of nucleic acid hybridization. However, despite their utility and popularity, DAs have not been well defined in the literature, leading to confusion over the differences between DAs and other aptasensor formats. In this review, we introduce a framework for DAs based on ACEs, and use this framework to distinguish DAs from other aptasensor formats and to categorize cis- and trans-DA designs. We then explore the ligand binding dynamics and chemical properties that underpin DA systems, which fall under conformational selection and induced fit models, and which mirror classical SN1 and SN2 models of nucleophilic substitution reactions. We further review a variety of in vitro and in vivo applications of DAs in the chemical and biological sciences, including riboswitches and riboregulators. Finally, we present future directions of DAs as ligand-responsive nucleic acids. Owing to their tractability, versatility and ease of engineering, DA biosensors bear a great potential for the development of new applications and technologies in fields ranging from analytical chemistry and mechanistic modeling to medicine and synthetic biology.

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.

Recent advances and future trends of riboswitches: attractive regulatory tools

Bacterial genomes contain a huge amount of different genes. These genes are spatiotemporally expressed to accomplish some required functions within the organism. Inside the cell, any step of gene expression may be modulated at four possible places such as transcription initiation, translation regulation, mRNA stability and protein stability. To achieve this, there is a necessity of strong regulators either natural or synthetic which can fine-tune gene expression regarding the required function. In recent years, riboswitches as metabolite responsive control elements residing in the untranslated regions of certain messenger RNAs, have been known to control gene expression at transcription or translation level. Importantly, these control elements do not prescribe the involvement of protein factors for metabolite binding. However, they own their particular properties to sense intramolecular metabolites (ligands). Herein, we highlighted current important bacterial riboswitches, their applications to support genetic control, ligand-binding domain mechanisms and current progress in synthetic riboswitches.

Structure-Activity Relationship of Flavin Analogues That Target the Flavin Mononucleotide Riboswitch

The flavin mononucleotide (FMN) riboswitch is an emerging target for the development of novel RNA-targeting antibiotics. We previously discovered an FMN derivative, 5FDQD, that protects mice against diarrhea-causing Clostridium difficile bacteria. Here, we present the structure-based drug design strategy that led to the discovery of this fluoro-phenyl derivative with antibacterial properties. This approach involved the following stages: (1) structural analysis of all available free and bound FMN riboswitch structures; (2) design, synthesis, and purification of derivatives; (3) in vitro testing for productive binding using two chemical probing methods; (4) in vitro transcription termination assays; and (5) resolution of the crystal structures of the FMN riboswitch in complex with the most mature candidates. In the process, we delineated principles for productive binding to this riboswitch, thereby demonstrating the effectiveness of a coordinated structure-guided approach to designing drugs against RNA.

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.

Crowding Shifts the FMN Recognition Mechanism of Riboswitch Aptamer from Conformational Selection to Induced Fit

In bacteria, the binding between the riboswitch aptamer domain and ligand is regulated by environmental cues, such as low Mg²⁺ in macrophages during pathogenesis to ensure spatiotemporal expression of virulence genes. Binding was investigated between the flavin mononucleotide (FMN) riboswitch aptamer and its anionic ligand in the presence of molecular crowding agent without Mg²⁺ ion, which mimics pathogenic conditions. Structural, kinetic, and thermodynamic analyses under the crowding revealed more dynamic conformational rearrangements of the FMN riboswitch aptamer compared to dilute Mg²⁺ -containing solution. It is hypothesized that under crowding conditions FMN binds through an induced fit mechanism in contrast to the conformational selection mechanism previously demonstrated in dilute Mg²⁺ solution. Since these two mechanisms involve different conformational intermediates and rate constants, these findings have practical significance in areas such as drug design and RNA engineering.

Discovery of Selective RNA-Binding Small Molecules by Affinity-Selection Mass Spectrometry

Recent advances in understanding the relevance of noncoding RNA (ncRNA) to disease have increased interest in drugging ncRNA with small molecules. The recent discovery of ribocil, a structurally distinct synthetic mimic of the natural ligand of the flavin mononucleotide (FMN) riboswitch, has revealed the potential chemical diversity of small molecules that target ncRNA. Affinity-selection mass spectrometry (AS-MS) is theoretically applicable to high-throughput screening (HTS) of small molecules binding to ncRNA. Here, we report the first application of the Automated Ligand Detection System (ALIS), an indirect AS-MS technique, for the selective detection of small molecule-ncRNA interactions, high-throughput screening against large unbiased small-molecule libraries, and identification and characterization of novel compounds (structurally distinct from both FMN and ribocil) that target the FMN riboswitch. Crystal structures reveal that different compounds induce various conformations of the FMN riboswitch, leading to different activity profiles. Our findings validate the ALIS platform for HTS screening for RNA-binding small molecules and further demonstrate that ncRNA can be broadly targeted by chemically diverse yet selective small molecules as therapeutics.

Real-time crystallographic studies of the adenine riboswitch using an X-ray free-electron laser

Structures of the four reaction states of the adenine riboswitch aptamer domain, including a transient intermediate state were solved by serial femtosecond crystallography. The structures not only demonstrate the use of X-ray free-electron lasers for RNA crystallography but have also proven that transient states can be determined in real time by mix-and-inject crystallography. These results illustrate the structural basis for the ligand-induced conformational changes associated with the molecular 'switch'.

An excited state underlies gene regulation of a transcriptional riboswitch

Riboswitches control gene expression through ligand-dependent structural rearrangements of the sensing aptamer domain. However, we found that the Bacillus cereus fluoride riboswitch aptamer adopts identical tertiary structures in solution with and without ligand. Using chemical-exchange saturation transfer (CEST) NMR spectroscopy, we revealed that the structured ligand-free aptamer transiently accesses a low-populated (∼1%) and short-lived (∼3 ms) excited conformational state that unravels a conserved 'linchpin' base pair to signal transcription termination. Upon fluoride binding, this highly localized, fleeting process is allosterically suppressed, which activates transcription. We demonstrated that this mechanism confers effective fluoride-dependent gene activation over a wide range of transcription rates, which is essential for robust toxicity responses across diverse cellular conditions. These results unveil a novel switching mechanism that employs ligand-dependent suppression of an aptamer excited state to coordinate regulatory conformational transitions rather than adopting distinct aptamer ground-state tertiary architectures, exemplifying a new mode of ligand-dependent RNA regulation.

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.