Sequence elements distal to the ligand binding pocket modulate the efficiency of a synthetic riboswitch

Investigating the Intramolecular Competition of Different RNA Binding Motifs for Neomycin B by Native Top-Down Mass Spectrometry

The ongoing search for small molecule drugs that target ribonucleic acids (RNA) is complicated by a limited understanding of the principles that govern RNA-small molecule interactions. Here we have used stoichiometry-resolved native top-down mass spectrometry (MS) to study the binding of neomycin B to small model hairpin RNAs, an unstructured RNA, and a viral RNA construct. For 15-22 nt model RNAs with hairpin structure, we found that neomycin B binding to hairpin loops relies on interactions with both the nucleobases and the 2'-OH groups, and that a simple 5' or 3' overhang can introduce an additional binding motif. For a 47 nt RNA construct derived from stem IA of the human immunodeficiency virus 1 (HIV-1) rev response element (RRE) RNA, native top-down MS identified four different binding motifs, of which the purine-rich internal loop showed the highest affinity for neomycin B. Stoichiometry-resolved binding site mapping by native top-down MS allows for a new perspective on binding specificity, and has the potential to reveal unexpected principles of small molecule binding to RNA.

Resolving the intricate binding of neomycin B to multiple binding motifs of a neomycin-sensing riboswitch aptamer by native top-down mass spectrometry and NMR spectroscopy

Understanding small molecule binding to RNA can be complicated by an intricate interplay between binding stoichiometry, multiple binding motifs, different occupancies of different binding motifs, and changes in the structure of the RNA under study. Here, we use native top-down mass spectrometry (MS) and nuclear magnetic resonance (NMR) spectroscopy to experimentally resolve these factors and gain a better understanding of the interactions between neomycin B and the 40 nt aptamer domain of a neomycin-sensing riboswitch engineered in yeast. Data from collisionally activated dissociation of the 1:1, 1:2 and 1:3 RNA-neomycin B complexes identified a third binding motif C of the riboswitch in addition to the two motifs A and B found in our previous study, and provided occupancies of the different binding motifs for each complex stoichiometry. Binding of a fourth neomycin B molecule was unspecific according to both MS and NMR data. Intriguingly, all major changes in the aptamer structure can be induced by the binding of the first neomycin B molecule regardless of whether it binds to motif A or B as evidenced by stoichiometry-resolved MS data together with titration data from 1H NMR spectroscopy in the imino proton region. Specific binding of the second and third neomycin B molecules further stabilizes the riboswitch aptamer, thereby allowing for a gradual response to increasing concentrations of neomycin B, which likely leads to a fine-tuning of the cellular regulatory mechanism.

Molecular dynamics in multidimensional space explains how mutations affect the association path of neomycin to a riboswitch

Riboswitches are messenger RNA (mRNA) fragments binding specific small molecules to regulate gene expression. A synthetic N1 riboswitch, inserted into yeast mRNA controls the translation of a reporter gene in response to neomycin. However, its regulatory activity is sensitive to single-point RNA mutations, even those distant from the neomycin binding site. While the association paths of neomycin to N1 and its variants remain unknown, recent fluorescence kinetic experiments indicate a two-step process driven by conformational selection. This raises the question of which step is affected by mutations. To address this, we performed all-atom two-dimensional replica-exchange molecular dynamics simulations for N1 and U14C, U14C[Formula: see text], U15A, and A17G mutants, ensuring extensive conformational sampling of both RNA and neomycin. The obtained neomycin association and binding paths, along with multidimensional free-energy profiles, revealed a two-step binding mechanism, consisting of conformational selection and induced fit. Neomycin binds to a preformed N1 conformation upon identifying a stable upper stem and U-turn motif in the riboswitch hairpin. However, the positioning of neomycin in the binding site occurs at different RNA-neomycin distances for each mutant, which may explain their different regulatory activities. The subsequent induced fit arises from the interactions of the neomycin's N3 amino group with RNA, causing the G9 backbone to rearrange. In the A17G mutant, the critical C6-A17/G17 stacking forms at a closer RNA-neomycin distance compared to N1. These findings together with estimated binding free energies coincide with experiments and elucidate why the A17G mutation decreases and U15A enhances N1 activity in response to neomycin.

The Impact of Second-Shell Nucleotides on Ligand Specificity in Cyclic Dinucleotide Riboswitches

Ligand specificity is an essential requirement for all riboswitches. Some variant riboswitches utilize a common structural motif, yet through subtle sequence differences, they are able to selectively respond to different small molecule ligands and regulate downstream gene expression. These variants discriminate between structurally and chemically similar ligands. Crystal structures provide insight into how specificity is achieved. However, ligand specificity cannot always be explained solely by nucleotides in direct contact with the ligand. The cyclic dinucleotide variant family contains two classes, cyclic-di-GMP and cyclic-AMP-GMP riboswitches, that were distinguished based on the identity of a single nucleotide in contact with the ligand. Here we report a variant riboswitch with a mutation at a second ligand-contacting position that is promiscuous for both cyclic-di-GMP and cyclic-AMP-GMP despite a predicted preference for cyclic-AMP-GMP. A high-throughput mutational analysis, SMARTT, was used to quantitatively assess thousands of sites in the first- and second-shells of ligand contact for impacts on ligand specificity and promiscuity. In addition to nucleotides in direct ligand contact, nucleotides more distal from the binding site, within the J1/2 linker and the terminator helix, were identified that impact ligand specificity. These findings provide an example of how nucleotides outside the ligand binding pocket influence the riboswitch specificity. Moreover, these distal nucleotides could be used to predict promiscuous sequences.

Riboswitches for Controlled Expression of Therapeutic Transgenes Delivered by Adeno-Associated Viral Vectors

Vectors developed from adeno-associated virus (AAV) are powerful tools for in vivo transgene delivery in both humans and animal models, and several AAV-delivered gene therapies are currently approved for clinical use. However, AAV-mediated gene therapy still faces several challenges, including limited vector packaging capacity and the need for a safe, effective method for controlling transgene expression during and after delivery. Riboswitches, RNA elements which control gene expression in response to ligand binding, are attractive candidates for regulating expression of AAV-delivered transgene therapeutics because of their small genomic footprints and non-immunogenicity compared to protein-based expression control systems. In addition, the ligand-sensing aptamer domains of many riboswitches can be exchanged in a modular fashion to allow regulation by a variety of small molecules, proteins, and oligonucleotides. Riboswitches have been used to regulate AAV-delivered transgene therapeutics in animal models, and recently developed screening and selection methods allow rapid isolation of riboswitches with novel ligands and improved performance in mammalian cells. This review discusses the advantages of riboswitches in the context of AAV-delivered gene therapy, the subsets of riboswitch mechanisms which have been shown to function in human cells and animal models, recent progress in riboswitch isolation and optimization, and several examples of AAV-delivered therapeutic systems which might be improved by riboswitch regulation.

Mutations of N1 Riboswitch Affect its Dynamics and Recognition by Neomycin Through Conformational Selection

Short, structured fragments of non-coding mRNA may act as molecular switches upon binding specific ligands, regulating the translation of proteins encoded downstream this mRNA sequence. One switch, called riboswitch N1, is regulated by aminoglycosides such as neomycin. Nucleobase mutations in the apical loop, although distant from the binding pocket, significantly affect neomycin affinity and riboswitch regulatory efficiency. To explain this influence, we conducted molecular dynamics simulations using generalized replica exchange with solute tempering (gREST). Translation assay of a reporter protein in a yeast system shows that mutating A17 to G in the riboswitch apical loop reduces 6-fold the translation regulation efficiency of the mutant. Indeed, simulations of the unbound riboswitch show that G17 frequently stacks with base 7, while base 8 is stabilized towards the binding site in a way that it may interfere with the conformational selection mechanism and decrease riboswitch regulatory activity. In the riboswitch complexes, this single-point A to G mutation disrupts a strong hydrogen bond between nucleotides 5 and 17 and, instead, a new hydrogen bond between residue 17 and neomycin is created. This change forces neomycin to occupy a slightly shifted position in the binding pocket, which increases neomycin flexibility. Our simulations of the U14C mutation suggest that the riboswitch complex with neomycin is more stable if cytosine 14 is protonated. A hydrogen bond between the RNA phosphate and protonated cytosine appears as the stabilizing factor. Also, based on the cell-free translation assay and isothermal titration calorimetry experiments, mutations of nucleotides 14 and 15 affect only slightly the riboswitch ability to bind the ligand and its activity. Indeed, the simulation of the unbound U15A mutant suggests conformations preformed for ligand binding, which may explain slightly higher regulatory activity of this mutant. Overall, our results corroborate the in vivo and in vitro experiments on the N1 riboswitch-neomycin system, detail the relationship between nucleobase mutations and RNA dynamics, and reveal the conformations playing the major role in the conformational selection mechanism.

Structure guided fluorescence labeling reveals a two-step binding mechanism of neomycin to its RNA aptamer

The ability of the cytidine analog Ç_m^f to act as a position specific reporter of RNA-dynamics was spectroscopically evaluated. Ç_m^f-labeled single- and double-stranded RNAs differ in their fluorescence lifetimes, quantum yields and anisotropies. These observables were also influenced by the nucleobases flanking Ç_m^f. This conformation and position specificity allowed to investigate the binding dynamics and mechanism of neomycin to its aptamer N1 by independently incorporating Ç_m^f at four different positions within the aptamer. Remarkably fast binding kinetics of neomycin binding was observed with stopped-flow measurements, which could be satisfactorily explained with a two-step binding. Conformational selection was identified as the dominant mechanism.

Molecular mechanisms for dynamic regulation of N1 riboswitch by aminoglycosides

A synthetic riboswitch N1, inserted into the 5'-untranslated mRNA region of yeast, regulates gene expression upon binding ribostamycin and neomycin. Interestingly, a similar aminoglycoside, paromomycin, differing from neomycin by only one substituent (amino versus hydroxyl), also binds to the N1 riboswitch, but without affecting gene expression, despite NMR evidence that the N1 riboswitch binds all aminoglycosides in a similar way. Here, to explore the details of structural dynamics of the aminoglycoside-N1 riboswitch complexes, we applied all-atom molecular dynamics (MD) and temperature replica-exchange MD simulations in explicit solvent. Indeed, we found that ribostamycin and neomycin affect riboswitch dynamics similarly but paromomycin allows for more flexibility because its complex lacks the contact between the distinctive 6' hydroxyl group and the G9 phosphate. Instead, a transient hydrogen bond of 6'-OH with A17 is formed, which partially diminishes interactions between the bulge and apical loop of the riboswitch, likely contributing to riboswitch inactivity. In many ways, the paromomycin complex mimics the conformations, interactions, and Na+ distribution of the free riboswitch. The MD-derived interaction network helps understand why riboswitch activity depends on aminoglycoside type, whereas for another aminoglycoside-binding site, aminoacyl-tRNA site in 16S rRNA, activity is not discriminatory.

Development of Artificial Riboswitches for Monitoring of Naringenin In Vivo

Microbial strains are considered promising hosts for production of flavonoids because of their rapid growth rate and suitability for large-scale manufacturing. However, productivity and titer of current recombinant strains still do not meet the requirements of industrial processes. Genetically encoded biosensors have been applied for high-throughput screening or dynamic regulation of biosynthetic pathways for enhancing the performance of microbial strains that produce valuable chemicals. Currently, few protein sensor-regulators for flavonoids exist. Unlike the protein-based trans-regulating controllers, riboswitches can respond to their ligands faster and eliminate off-target effects. Here, we developed artificial riboswitches that activate gene expression in response to naringenin, an important flavonoid. RNA aptamers for naringenin were developed using SELEX and cloned upstream of a dual selectable marker gene to construct a riboswitch library. Two in vivo selection routes were applied separately to the library by supplementing naringenin at two different concentrations during enrichments to modulate the operational ranges of the riboswitches. The selected riboswitches were responsive to naringenin and activated gene expression up to 2.91-fold. Operational ranges of the riboswitches were distinguished on the basis of their selection route. Additionally, the specificity of the riboswitches was assessed, and their applicability as dynamic regulators was confirmed. Collectively, the naringenin riboswitches reported in this work will be valuable tools in metabolic engineering of microorganisms for the production of flavonoids.

Synthetic Riboswitches: From Plug and Pray toward Plug and Play

In synthetic biology, metabolic engineering, and gene therapy, there is a strong demand for orthogonal or externally controlled regulation of gene expression. Here, RNA-based regulatory devices represent a promising emerging alternative to proteins, allowing a fast and direct control of gene expression, as no synthesis of regulatory proteins is required. Besides programmable ribozyme elements controlling mRNA stability, regulatory RNA structures in untranslated regions are highly interesting for engineering approaches. Riboswitches are especially well suited, as they show a modular composition of sensor and response elements, allowing a free combination of different modules in a plug-and-play-like mode. The sensor or aptamer domain specifically interacts with a trigger molecule as a ligand, modulating the activity of the adjacent response domain that controls the expression of the genes located downstream, in most cases at the level of transcription or translation. In this review, we discuss the recent advances and strategies for designing such synthetic riboswitches based on natural or artificial components and readout systems, from trial-and-error approaches to rational design strategies. As the past several years have shown dramatic development in this fascinating field of research, we can give only a limited overview of the basic riboswitch design principles that is far from complete, and we apologize for not being able to consider every successful and interesting approach described in the literature.

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.

Riboswitch engineering - making the all-important second and third steps

Synthetic biology uses our understanding of biological systems to develop innovative solutions for challenges in fields as diverse as genetic control and logic devices, bioremediation, materials production or diagnostics and therapy in medicine by designing new biological components. RNA-based elements are key components of these engineered systems. Their structural and functional diversity is ideal for generating regulatory riboswitches that react with many different types of output to molecular and environmental signals. Recent advances have added new sensor and output domains to the existing toolbox, and demonstrated the portability of riboswitches to many different organisms. Improvements in riboswitch design and screens for selecting in vivo active switches provide the means to isolate riboswitches with regulatory properties more like their natural counterparts.

Building a stable RNA U-turn with a protonated cytidine

The U-turn is a classical three-dimensional RNA folding motif first identified in the anticodon and T-loops of tRNAs. It also occurs frequently as a building block in other functional RNA structures in many different sequence and structural contexts. U-turns induce sharp changes in the direction of the RNA backbone and often conform to the 3-nt consensus sequence 5'-UNR-3' (N = any nucleotide, R = purine). The canonical U-turn motif is stabilized by a hydrogen bond between the N3 imino group of the U residue and the 3' phosphate group of the R residue as well as a hydrogen bond between the 2'-hydroxyl group of the uridine and the N7 nitrogen of the R residue. Here, we demonstrate that a protonated cytidine can functionally and structurally replace the uridine at the first position of the canonical U-turn motif in the apical loop of the neomycin riboswitch. Using NMR spectroscopy, we directly show that the N3 imino group of the protonated cytidine forms a hydrogen bond with the backbone phosphate 3' from the third nucleotide of the U-turn analogously to the imino group of the uridine in the canonical motif. In addition, we compare the stability of the hydrogen bonds in the mutant U-turn motif to the wild type and describe the NMR signature of the C+-phosphate interaction. Our results have implications for the prediction of RNA structural motifs and suggest simple approaches for the experimental identification of hydrogen bonds between protonated C-imino groups and the phosphate backbone.