A disconnect between high-affinity binding and efficient regulation by antifolates and purines in the tetrahydrofolate riboswitch

Assessment of Nucleobase Protomeric and Tautomeric States in Nucleic Acid Structures for Interaction Analysis and Structure-Based Ligand Design

With increasing interest in RNA as a therapeutic and a potential target, the role of RNA structures has become more important. Even slight changes in nucleobases, such as modifications or protomeric and tautomeric states, can have a large impact on RNA structure and function, while local environments in turn affect protonation and tautomerization. In this work, the application of empirical tools for pKa and tautomer prediction for RNA modifications was elucidated and compared with ab initio quantum mechanics (QM) methods and expanded toward macromolecular RNA structures, where QM is no longer feasible. In this regard, the Protonate3D functionality within the molecular operating environment (MOE) was expanded for nucleobase protomer and tautomer predictions and applied to reported examples of altered protonation states depending on the local environment. Overall, observations of nonstandard protomers and tautomers were well reproduced, including structural C+G:C(A) and A+GG motifs, several mismatches, and protonation of adenosine or cytidine as the general acid in nucleolytic ribozymes. Special cases, such as cobalt hexamine-soaked complexes or the deprotonation of guanosine as the general base in nucleolytic ribozymes, proved to be challenging. The collected set of examples shall serve as a starting point for the development of further RNA protonation prediction tools, while the presented Protonate3D implementation already delivers reasonable protonation predictions for RNA and DNA macromolecules. For cases where higher accuracy is needed, like following catalytic pathways of ribozymes, incorporation of QM-based methods can build upon the Protonate3D-generated starting structures. Likewise, this protonation prediction can be used for structure-based RNA-ligand design approaches.

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.

Structural insights into translation regulation by the THF-II riboswitch

In bacteria, expression of folate-related genes is controlled by the tetrahydrofolate (THF) riboswitch in response to specific binding of THF and its derivatives. Recently, a second class of THF riboswitches, named THF-II, was identified in Gram-negative bacteria, which exhibit distinct architecture from the previously characterized THF-I riboswitches found in Gram-positive bacteria. Here, we present the crystal structures of the ligand-bound THF-II riboswitch from Mesorhizobium loti. These structures exhibit a long rod-like fold stabilized by continuous base pair and base triplet stacking across two helices of P1 and P2 and their interconnecting ligand-bound binding pocket. The pterin moiety of the ligand docks into the binding pocket by forming hydrogen bonds with two highly conserved pyrimidines in J12 and J21, which resembles the hydrogen-bonding pattern at the ligand-binding site FAPK in the THF-I riboswitch. Using small-angle X-ray scattering and isothermal titration calorimetry, we further characterized the riboswitch in solution and reveal that Mg2+ is essential for pre-organization of the binding pocket for efficient ligand binding. RNase H cleavage assay indicates that ligand binding reduces accessibility of the ribosome binding site in the right arm of P1, thus down-regulating the expression of downstream genes. Together, these results provide mechanistic insights into translation regulation by the THF-II riboswitch.

Isothermal Titration Calorimetry Analysis of a Cooperative Riboswitch Using an Interdependent-Sites Binding Model

Isothermal titration calorimetry (ITC) is a powerful biophysical tool to characterize energetic profiles of biomacromolecular interactions without any alteration of the underlying chemical structures. In this protocol, we describe procedures for performing, analyzing, and interpreting ITC data obtained from a cooperative riboswitch-ligand interaction.

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

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

Cobalamin Riboswitches Are Broadly Sensitive to Corrinoid Cofactors to Enable an Efficient Gene Regulatory Strategy

In bacteria, many essential metabolic processes are controlled by riboswitches, gene regulatory RNAs that directly bind and detect metabolites. Highly specific effector binding enables riboswitches to respond to a single biologically relevant metabolite. Cobalamin riboswitches are a potential exception because over a dozen chemically similar but functionally distinct cobalamin variants (corrinoid cofactors) exist in nature. Here, we measured cobalamin riboswitch activity in vivo using a Bacillus subtilis fluorescent reporter system and found, among 38 tested riboswitches, a subset responded to corrinoids promiscuously, while others were semiselective. Analyses of chimeric riboswitches and structural models indicate, unlike other riboswitch classes, cobalamin riboswitches indirectly differentiate among corrinoids by sensing differences in their structural conformation. This regulatory strategy aligns riboswitch-corrinoid specificity with cellular corrinoid requirements in a B. subtilis model. Thus, bacteria can employ broadly sensitive riboswitches to cope with the chemical diversity of essential metabolites. IMPORTANCE Some bacterial mRNAs contain a region called a riboswitch which controls gene expression by binding to a metabolite in the cell. Typically, riboswitches sense and respond to a limited range of cellular metabolites, often just one type. In this work, we found the cobalamin (vitamin B₁₂) riboswitch class is an exception, capable of sensing and responding to multiple variants of B₁₂-collectively called corrinoids. We found cobalamin riboswitches vary in corrinoid specificity with some riboswitches responding to each of the corrinoids we tested, while others responding only to a subset of corrinoids. Our results suggest the latter class of riboswitches sense intrinsic conformational differences among corrinoids in order to support the corrinoid-specific needs of the cell. These findings provide insight into how bacteria sense and respond to an exceptionally diverse, often essential set of enzyme cofactors.

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.

Protein-Based Virtual Screening Tools Applied for RNA-Ligand Docking Identify New Binders of the preQ1-Riboswitch

Targeting RNA with small molecules is an emerging field. While several ligands for different RNA targets are reported, structure-based virtual screenings (VSs) against RNAs are still rare. Here, we elucidated the general capabilities of protein-based docking programs to reproduce native binding modes of small-molecule RNA ligands and to discriminate known binders from decoys by the scoring function. The programs were found to perform similar compared to the RNA-based docking tool rDOCK, and the challenges faced during docking, namely, protomer and tautomer selection, target dynamics, and explicit solvent, do not largely differ from challenges in conventional protein-ligand docking. A prospective VS with the Bacillus subtilis preQ₁-riboswitch aptamer domain performed with FRED, HYBRID, and FlexX followed by microscale thermophoresis assays identified six active compounds out of 23 tested VS hits with potencies between 29.5 nM and 11.0 μM. The hits were selected not solely based on their docking score but for resembling key interactions of the native ligand. Therefore, this study demonstrates the general feasibility to perform structure-based VSs against RNA targets, while at the same time it highlights pitfalls and their potential solutions when executing RNA-ligand docking.

The Second Class of Tetrahydrofolate (THF-II) Riboswitches Recognizes the Tetrahydrofolic Acid Ligand via Local Conformation Changes

Riboswitches are regulatory noncoding RNAs found in bacteria, fungi and plants, that modulate gene expressions through structural changes in response to ligand binding. Understanding how ligands interact with riboswitches in solution can shed light on the molecular mechanisms of this ancient regulators. Previous studies showed that riboswitches undergo global conformation changes in response to ligand binding to relay information. Here, we report conformation switching models of the recently discovered tetrahydrofolic acid-responsive second class of tetrahydrofolate (THF-II) riboswitches in response to ligand binding. Using a combination of selective 2′-hydroxyl acylation, analyzed by primer extension (SHAPE) assay, 3D modeling and small-angle X-ray scattering (SAXS), we found that the ligand specifically recognizes and reshapes the THF-II riboswitch loop regions, but does not affect the stability of the P3 helix. Our results show that the THF-II riboswitch undergoes only local conformation changes in response to ligand binding, rearranging the Loop1-P3-Loop2 region and rotating Loop1 from a ~120° angle to a ~75° angle. This distinct conformation changes suggest a unique regulatory mechanism of the THF-II riboswitch, previously unseen in other riboswitches. Our findings may contribute to the fields of RNA sensors and drug design.

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

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

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.

Tying the knot in the tetrahydrofolate (THF) riboswitch: A molecular basis for gene regulation

Effective gene regulation by the tetrahydrofolate riboswitch depends not only on ligand affinity but also on the kinetics of ligand association, which involves two cooperative binding sites. We have determined a 1.9-Å resolution crystal structure of the ligand-free THF riboswitch aptamer. The pseudoknot binding site 'unwinds' in the absence of ligand, whereby the adjacent helical domains (P1, P2, and P3) become disjointed, resulting in rotation and misalignment of the gene-regulatory P1 helix with respect to P3. In contrast, the second binding site at the three-way junction, which is the first to fold, is structurally conserved between apo and holo forms. This suggests a kinetic role for this site, in which binding of the first ligand molecule to the stably folded three-way junction promotes formation of the regulatory pseudoknot site and subsequent binding of the second molecule. As such, these findings provide a molecular basis for both conformational switching and kinetic control.

Occurrence and stability of lone pair-π and OH-π interactions between water and nucleobases in functional RNAs

We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp-π) or the OH-π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1-3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp-π interaction as compared to the OH-π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp-π or the OH-π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.

Parallel Discovery Strategies Provide a Basis for Riboswitch Ligand Design

Riboswitches are mRNA domains that make gene-regulatory decisions upon binding their cognate ligands. Bacterial riboswitches that specifically recognize 5-aminoimidazole-4-carboxamide riboside 5'-monophosphate (ZMP) and 5'-triphosphate (ZTP) regulate genes involved in folate and purine metabolism. Now, we have developed synthetic ligands targeting ZTP riboswitches by replacing the sugar-phosphate moiety of ZMP with various functional groups, including simple heterocycles. Despite losing hydrogen bonds from ZMP, these analogs bind ZTP riboswitches with similar affinities as the natural ligand, and activate transcription more strongly than ZMP in vitro. The most active ligand stimulates gene expression ∼3 times more than ZMP in a live Escherichia coli reporter. Co-crystal structures of the Fusobacterium ulcerans ZTP riboswitch bound to synthetic ligands suggest stacking of their pyridine moieties on a conserved RNA nucleobase primarily determines their higher activity. Altogether, these findings guide future design of improved riboswitch activators and yield insights into how RNA-targeted ligand discovery may proceed.

Small molecule targeting of RNA structures in neurological disorders

Aberrant RNA structure and function operate in neurological disease progression and severity. As RNA contributes to disease pathology in a complex fashion, that is, via various mechanisms, it has become an attractive therapeutic target for small molecules and oligonucleotides. In this review, we discuss the identification of RNA structures that cause or contribute to neurological diseases as well as recent progress toward the development of small molecules that target them, including small molecule modulators of pre-mRNA splicing and RNA repeat expansions that cause microsatellite disorders such as Huntington's disease and amyotrophic lateral sclerosis. The use of oligonucleotide-based modalities is also discussed. There are key differences between small molecule and oligonucleotide targeting of RNA. The former targets RNA structure, while the latter prefers unstructured regions. Thus, some targets will be preferentially targeted by oligonucleotides and others by small molecules.

A unified dinucleotide alphabet describing both RNA and DNA structures

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

High Affinity Binding of N2-Modified Guanine Derivatives Significantly Disrupts the Ligand Binding Pocket of the Guanine Riboswitch

Riboswitches are important model systems for the development of approaches to search for RNA-targeting therapeutics. A principal challenge in finding compounds that target riboswitches is that the effector ligand is typically almost completely encapsulated by the RNA, which severely limits the chemical space that can be explored. Efforts to find compounds that bind the guanine/adenine class of riboswitches with a high affinity have in part focused on purines modified at the C6 and C2 positions. These studies have revealed compounds that have low to sub-micromolar affinity and, in a few cases, have antimicrobial activity. To further understand how these compounds interact with the guanine riboswitch, we have performed an integrated structural and functional analysis of representative guanine derivatives with modifications at the C8, C6 and C2 positions. Our data indicate that while modifications of guanine at the C6 position are generally unfavorable, modifications at the C8 and C2 positions yield compounds that rival guanine with respect to binding affinity. Surprisingly, C2-modified guanines such as N2-acetylguanine completely disrupt a key Watson–Crick pairing interaction between the ligand and RNA. These compounds, which also modulate transcriptional termination as efficiently as guanine, open up a significant new chemical space of guanine modifications in the search for antimicrobial agents that target purine riboswitches.

Responsive self-assembly of tectoRNAs with loop-receptor interactions from the tetrahydrofolate (THF) riboswitch

Naturally occurring RNAs are known to exhibit a high degree of modularity, whereby specific structural modules (or motifs) can be mixed and matched to create new molecular architectures. The modular nature of RNA also affords researchers the ability to characterize individual structural elements in controlled synthetic contexts in order to gain new and critical insights into their particular structural features and overall performance. Here, we characterized the binding affinity of a unique loop–receptor interaction found in the tetrahydrofolate (THF) riboswitch using rationally designed self-assembling tectoRNAs. Our work suggests that the THF loop–receptor interaction has been fine-tuned for its particular role as a riboswitch component. We also demonstrate that the thermodynamic stability of this interaction can be modulated by the presence of folinic acid, which induces a local structural change at the level of the loop–receptor. This corroborates the existence of a THF binding site within this tertiary module and paves the way for its potential use as a THF responsive module for RNA nanotechnology and synthetic biology.

RNA 3D structure prediction guided by independent folding of homologous sequences

BACKGROUND

The understanding of the importance of RNA has dramatically changed over recent years. As in the case of proteins, the function of an RNA molecule is encoded in its tertiary structure, which in turn is determined by the molecule's sequence. The prediction of tertiary structures of complex RNAs is still a challenging task.

RESULTS

Using the observation that RNA sequences from the same RNA family fold into conserved structure, we test herein whether parallel modeling of RNA homologs can improve ab initio RNA structure prediction. EvoClustRNA is a multi-step modeling process, in which homologous sequences for the target sequence are selected using the Rfam database. Subsequently, independent folding simulations using Rosetta FARFAR and SimRNA are carried out. The model of the target sequence is selected based on the most common structural arrangement of the common helical fragments. As a test, on two blind RNA-Puzzles challenges, EvoClustRNA predictions ranked as the first of all submissions for the L-glutamine riboswitch and as the second for the ZMP riboswitch. Moreover, through a benchmark of known structures, we discovered several cases in which particular homologs were unusually amenable to structure recovery in folding simulations compared to the single original target sequence.

CONCLUSION

This work, for the first time to our knowledge, demonstrates the importance of the selection of the target sequence from an alignment of an RNA family for the success of RNA 3D structure prediction. These observations prompt investigations into a new direction of research for checking 3D structure "foldability" or "predictability" of related RNA sequences to obtain accurate predictions. To support new research in this area, we provide all relevant scripts in a documented and ready-to-use form. By exploring new ideas and identifying limitations of the current RNA 3D structure prediction methods, this work is bringing us closer to the near-native computational RNA 3D models.

Pairing interactions between nucleobases and ligands in aptamer:ligand complexes of riboswitches: crystal structure analysis, classification, optimal structures, and accurate interaction energies

In the present work, 67 crystal structures of the aptamer domains of RNA riboswitches are chosen for analysis of the structure and strength of hydrogen bonding (pairing) interactions between nucleobases constituting the aptamer binding pockets and the bound ligands. A total of 80 unique base:ligand hydrogen-bonded pairs containing at least two hydrogen bonds were identified through visual inspection. Classification of these contacts in terms of the interacting edge of the aptamer nucleobase revealed that interactions involving the Watson-Crick edge are the most common, followed by the sugar edge of purines and the Hoogsteen edge of uracil. Alternatively, classification in terms of the chemical constitution of the ligand yields five unique classes of base:ligand pairs: base:base, base:amino acid, base:sugar, base:phosphate, and base:other. Further, quantum mechanical (QM) geometry optimizations revealed that 67 out of 80 pairs exhibit stable geometries and optimal deviations from their macromolecular crystal occurrences. This indicates that these contacts are well-defined RNA aptamer:ligand interaction motifs. QM calculated interaction energies of base:ligand pairs reveal a rich hydrogen bonding landscape, ranging from weak interactions (base:other, -3 kcal/mol) to strong (base:phosphate, -48 kcal/mol) contacts. The analysis was further extended to study the biological importance of base:ligand interactions in the binding pocket of the tetrahydrofolate riboswitch and thiamine pyrophosphate riboswitch. Overall, our study helps in understanding the structural and energetic features of base:ligand pairs in riboswitches, which could aid in developing meaningful hypotheses in the context of RNA:ligand recognition. This can, in turn, contribute toward current efforts to develop antimicrobials that target RNAs.

3dRNA v2.0: An Updated Web Server for RNA 3D Structure Prediction

3D structures of RNAs are the basis for understanding their biological functions. However, experimentally solved RNA 3D structures are very limited in comparison with known RNA sequences up to now. Therefore, many computational methods have been proposed to solve this problem, including our 3dRNA. In recent years, 3dRNA has been greatly improved by adding several important features, including structure sampling, structure ranking and structure optimization under residue-residue restraints. Particularly, the optimization procedure with restraints enables 3dRNA to treat pseudoknots in a new way. These new features of 3dRNA can greatly promote its performance and have been integrated into the 3dRNA v2.0 web server. Here we introduce these new features in the 3dRNA v2.0 web server for the users.

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.

Ligand design for riboswitches, an emerging target class for novel antibiotics

Riboswitches are cis-acting gene regulatory elements and constitute potential targets for new antibiotics. Recent studies in this field have started to explore these targets for drug discovery. New ligands found by fragment screening, design of analogs of the natural ligands or serendipitously by phenotypic screening have shown antibacterial effects in cell assays against a range of bacteria strains and in animal models. In this review, we highlight the most advanced drug design work of riboswitch ligands and discuss the challenges in the field with respect to the development of antibiotics with a new mechanism of action.

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.

Singlet glycine riboswitches bind ligand as well as tandem riboswitches

The glycine riboswitch often occurs in a tandem architecture, with two ligand-binding domains (aptamers) followed by a single expression platform. Based on previous observations, we hypothesized that "singlet" versions of the glycine riboswitch, which contain only one aptamer domain, are able to bind glycine if appropriate structural contacts are maintained. An initial alignment of 17 putative singlet riboswitches indicated that the single consensus aptamer domain is flanked by a conserved peripheral stem-loop structure. These singlets were sorted into two subtypes based on whether the active aptamer domain precedes or follows the peripheral stem-loop, and an example of each subtype of singlet riboswitch was characterized biochemically. The singlets possess glycine-binding affinities comparable to those of previously published tandem examples, and the conserved peripheral domains form A-minor interactions with the single aptamer domain that are necessary for ligand-binding activity. Analysis of sequenced genomes identified a significant number of singlet glycine riboswitches. Based on these observations, we propose an expanded model for glycine riboswitch gene control that includes singlet and tandem architectures.

De Novo Guanine Biosynthesis but Not the Riboswitch-Regulated Purine Salvage Pathway Is Required for Staphylococcus aureus Infection In Vivo

De novo guanine biosynthesis is an evolutionarily conserved pathway that creates sufficient nucleotides to support DNA replication, transcription, and translation. Bacteria can also salvage nutrients from the environment to supplement the de novo pathway, but the relative importance of either pathway during Staphylococcus aureus infection is not known. In S. aureus, genes important for both de novo and salvage pathways are regulated by a guanine riboswitch. Bacterial riboswitches have attracted attention as a novel class of antibacterial drug targets because they have high affinity for small molecules, are absent in humans, and regulate the expression of multiple genes, including those essential for cell viability. Genetic and biophysical methods confirm the existence of a bona fide guanine riboswitch upstream of an operon encoding xanthine phosphoribosyltransferase (xpt), xanthine permease (pbuX), inosine-5'-monophosphate dehydrogenase (guaB), and GMP synthetase (guaA) that represses the expression of these genes in response to guanine. We found that S. aureus guaB and guaA are also transcribed independently of riboswitch control by alternative promoter elements. Deletion of xpt-pbuX-guaB-guaA genes resulted in guanine auxotrophy, failure to grow in human serum, profound abnormalities in cell morphology, and avirulence in mouse infection models, whereas deletion of the purine salvage genes xpt-pbuX had none of these effects. Disruption of guaB or guaA recapitulates the xpt-pbuX-guaB-guaA deletion in vivo In total, the data demonstrate that targeting the guanine riboswitch alone is insufficient to treat S. aureus infections but that inhibition of guaA or guaB could have therapeutic utility.

IMPORTANCE

De novo guanine biosynthesis and purine salvage genes were reported to be regulated by a guanine riboswitch in Staphylococcus aureus We demonstrate here that this is not true, because alternative promoter elements that uncouple the de novo pathway from riboswitch regulation were identified. We found that in animal models of infection, the purine salvage pathway is insufficient for S. aureus survival in the absence of de novo guanine biosynthesis. These data suggest targeting the de novo guanine biosynthesis pathway may have therapeutic utility in the treatment of S. aureus infections.

Recent Advances in Developing Small Molecules Targeting Nucleic Acid

Nucleic acids participate in a large number of biological processes. However, current approaches for small molecules targeting protein are incompatible with nucleic acids. On the other hand, the lack of crystallization of nucleic acid is the limiting factor for nucleic acid drug design. Because of the improvements in crystallization in recent years, a great many structures of nucleic acids have been reported, providing basic information for nucleic acid drug discovery. This review focuses on the discovery and development of small molecules targeting nucleic acids.

Isothermal Titration Calorimetry: Assisted Crystallization of RNA-Ligand Complexes

The success rate of nucleic acids/ligands co-crystallization can be significantly improved by performing preliminary biophysical analyses. Among suitable biophysical approaches, isothermal titration calorimetry (ITC) is certainly a method of choice. ITC can be used in a wide range of experimental conditions to monitor in real time the formation of the RNA- or DNA-ligand complex, with the advantage of providing in addition the complete binding profile of the interaction. Following the ITC experiment, the complex is ready to be concentrated for crystallization trials. This chapter describes a detailed experimental protocol for using ITC as a tool for monitoring RNA/small molecule binding, followed by co-crystallization.

In Vitro Screening and in Silico Modeling of RNA-Based Gene Expression Control

Molecular tools for controlling gene expression are essential for manipulating biological systems. One class of tools includes RNA switches that incorporate RNA-based sensors, known as aptamers. However, most switches reported to date are responsive to toxic molecules or to endogenous metabolites. For effective conditional control, switches must incorporate RNA aptamers that exhibit selectivity against such endogenous metabolites. We report a systematic approach which combines a rapid in vitro assay and an in silico model to support an efficient, streamlined application of aptamers into RNA switches. Model predictions were validated in vivo and demonstrate that the RNA switches enable selective and programmable gene regulation. We demonstrate the method using aptamers that bind the FDA-approved small molecule (6R)-folinic acid, providing access to new molecular targets for gene expression control and much-needed clinically relevant tools for advancing RNA-based therapeutics.

RNA aptamers as genetic control devices: the potential of riboswitches as synthetic elements for regulating gene expression

RNA utilizes many different mechanisms to control gene expression. Among the regulatory elements that respond to external stimuli, riboswitches are a prominent and elegant example. They consist solely of RNA and couple binding of a small molecule ligand to the so-called "aptamer domain" with a conformational change in the downstream "expression platform" which then determines system output. The modular organization of riboswitches and the relative ease with which ligand-binding RNA aptamers can be selected in vitro against almost any molecule have led to the rapid and widespread adoption of engineered riboswitches as artificial genetic control devices in biotechnology and synthetic biology over the past decade. This review highlights proof-of-principle applications to demonstrate the versatility and robustness of engineered riboswitches in regulating gene expression in pro- and eukaryotes. It then focuses on strategies and parameters to identify aptamers that can be integrated into synthetic riboswitches that are functional in vivo, before finishing with a reflection on how to improve the regulatory properties of engineered riboswitches, so that we can not only further expand riboswitch applicability, but also finally fully exploit their potential as control elements in regulating gene expression.

InterRNA: a database of base interactions in RNA structures

A major component of RNA structure stabilization are the hydrogen bonded interactions between the base residues. The importance and biological relevance for large clusters of base interactions can be much more easily investigated when their occurrences have been systematically detected, catalogued and compared. In this paper, we describe the database InterRNA (INTERactions in RNA structures database—http://mfrlab.org/interrna/) that contains records of known RNA 3D motifs as well as records for clusters of bases that are interconnected by hydrogen bonds. The contents of the database were compiled from RNA structural annotations carried out by the NASSAM (http://mfrlab.org/grafss/nassam) and COGNAC (http://mfrlab.org/grafss/cognac) computer programs. An analysis of the database content and comparisons with the existing corpus of knowledge regarding RNA 3D motifs clearly show that InterRNA is able to provide an extension of the annotations for known motifs as well as able to provide novel interactions for further investigations.

Linking aptamer-ligand binding and expression platform folding in riboswitches: prospects for mechanistic modeling and design

The power of riboswitches in regulation of bacterial metabolism derives from coupling of two characteristics: recognition and folding. Riboswitches contain aptamers, which function as biosensors. Upon detection of the signaling molecule, the riboswitch transduces the signal into a genetic decision. The genetic decision is coupled to refolding of the expression platform, which is distinct from, although overlapping with, the aptamer. Early biophysical studies of riboswitches focused on recognition of the ligand by the aptamer-an important consideration for drug design. A mechanistic understanding of ligand-induced riboswitch RNA folding can further enhance riboswitch ligand design, and inform efforts to tune and engineer riboswitches with novel properties. X-ray structures of aptamer/ligand complexes point to mechanisms through which the ligand brings together distal strand segments to form a P1 helix. Transcriptional riboswitches must detect the ligand and form this P1 helix within the timescale of transcription. Depending on the cell's metabolic state and cellular environmental conditions, the folding and genetic outcome may therefore be affected by kinetics of ligand binding, RNA folding, and transcriptional pausing, among other factors. Although some studies of isolated riboswitch aptamers found homogeneous, prefolded conformations, experimental, and theoretical studies point to functional and structural heterogeneity for nascent transcripts. Recently it has been shown that some riboswitch segments, containing the aptamer and partial expression platforms, can form binding-competent conformers that incorporate an incomplete aptamer secondary structure. Consideration of the free energy landscape for riboswitch RNA folding suggests models for how these conformers may act as transition states-facilitating rapid, ligand-mediated aptamer folding.

Cooperativity, allostery and synergism in ligand binding to riboswitches

Recent progress in identification and characterization of novel types of non-coding RNAs has proven that RNAs carry out a variety of cellular functions ranging from scaffolding to gene expression control. In both prokaryotic and eukaryotic cells, several classes of non-coding RNAs control expression of dozens of genes in response to specific cues. One of the most interesting and outstanding questions in the RNA field is whether regulatory RNAs are capable of employing basic biological concepts, such as allostery and cooperativity, previously attributed to the function of proteins. Aside from regulatory RNAs that form complementary base pairing with their nucleic acid targets, several RNA classes modulate gene expression via molecular mechanisms which can be paralleled to protein-mediated regulation. Among these RNAs are riboswitches, metabolite-sensing non-coding regulatory elements that adopt intrinsic three-dimensional structures and specifically bind various small molecule ligands. These characteristics of riboswitches make them well-suited for complex regulatory responses observed in allosteric and cooperative protein systems. Here we present an overview of the biochemical, genetic, and structural studies of riboswitches with a major focus on complex regulatory mechanisms and biological principles utilized by riboswitches for such genetic modulation.

Metal Ion-Mediated Nucleobase Recognition by the ZTP Riboswitch

The ZTP riboswitch is a widespread family of regulatory RNAs that upregulate de novo purine synthesis in response to increased intracellular levels of ZTP or ZMP. As an important intermediate in purine biosynthesis, ZMP also serves as a proxy for the concentration of N10-formyl-tetrahydrofolate, a key component of one-carbon metabolism. Here, we report the structure of the ZTP riboswitch bound to ZMP at a resolution of 1.80 Å. The RNA contains two subdomains brought together through a long-range pseudoknot further stabilized through helix-helix packing. ZMP is bound at the subdomain interface of the RNA through a set of interactions with the base, ribose sugar, and phosphate moieties of the ligand. Unique to nucleobase recognition by RNAs, the Z base is inner-sphere coordinated to a magnesium cation bound by two backbone phosphates. This interaction, along with steric hindrance by the backbone, imparts specificity over chemically similar compounds such as ATP/AMP.

Ligand binding by the tandem glycine riboswitch depends on aptamer dimerization but not double ligand occupancy

The glycine riboswitch predominantly exists as a tandem structure, with two adjacent, homologous ligand-binding domains (aptamers), followed by a single expression platform. The recent identification of a leader helix, the inclusion of which eliminates cooperativity between the aptamers, has reopened the debate over the purpose of the tandem structure of the glycine riboswitch. An equilibrium dialysis-based assay was combined with binding-site mutations to monitor glycine binding in each ligand-binding site independently to understand the role of each aptamer in glycine binding and riboswitch tertiary interactions. A series of mutations disrupting the dimer interface was used to probe how dimerization impacts ligand binding by the tandem glycine riboswitch. While the wild-type tandem riboswitch binds two glycine equivalents, one for each aptamer, both individual aptamers are capable of binding glycine when the other aptamer is unoccupied. Intriguingly, glycine binding by aptamer-1 is more sensitive to dimerization than glycine binding by aptamer-2 in the context of the tandem riboswitch. However, monomeric aptamer-2 shows dramatically weakened glycine-binding affinity. In addition, dimerization of the two aptamers in trans is dependent on glycine binding in at least one aptamer. We propose a revised model for tandem riboswitch function that is consistent with these results, wherein ligand binding in aptamer-1 is linked to aptamer dimerization and stabilizes the P1 stem of aptamer-2, which controls the expression platform.

Binding versus triggering riboswitches

In this issue of Chemistry & Biology, Trausch and Batey report a discrepancy between ligand binding affinity and the effect of transcription termination in a THF riboswitch, raising some important questions about our current understanding of ligand-dependent RNA switches.

The purine riboswitch as a model system for exploring RNA biology and chemistry

Over the past decade the purine riboswitch, and in particular its nucleobase-binding aptamer domain, has emerged as an important model system for exploring various aspects of RNA structure and function. Its relatively small size, structural simplicity and readily observable activity enable application of a wide variety of experimental approaches towards the study of this RNA. These analyses have yielded important insights into small molecule recognition, co-transcriptional folding and secondary structural switching, and conformational dynamics that serve as a paradigm for other RNAs. In this article, the current state of understanding of the purine riboswitch family and how this growing knowledge base is starting to be exploited in the creation of novel RNA devices are examined. This article is part of a Special Issue entitled: Riboswitches.

Gene regulation by structured mRNA elements

The precise temporal and spatial coordination of gene activity, based on the integration of internal and external signals, is crucial for the accurate functioning of all biological processes. Although the basic principles of gene expression were established some 60 years ago, recent research has revealed a surprising complexity in the control of gene activity. Many of these gene regulatory mechanisms occur at the level of the mRNA, including sophisticated gene control tasks mediated by structured mRNA elements. We now know that mRNA folds can serve as highly specific receptors for various types of molecules, as exemplified by metabolite-binding riboswitches, and interfere with pro- and eukaryotic gene expression at the level of transcription, translation, and RNA processing. Gene regulation by structured mRNA elements comprises versatile strategies including self-cleaving ribozymes, RNA-folding-mediated occlusion or presentation of cis-regulatory sequences, and sequestration of trans-acting factors including other RNAs and proteins.

Riboswitch structure and dynamics by smFRET microscopy

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