The dynamic nature of RNA as key to understanding riboswitch mechanisms

Mechanism of Dual-Site Recognition in a Classic DNA Aptamer

Nucleic acid aptamers possess unique advantages in specific recognition. However, the lack of in-depth investigation into their dynamic recognition mechanisms has restricted their rational design and potential applications in fields such as biosensing and targeted therapy. We herein utilized enhanced sampling molecular dynamics to address affinities of adenosine monophosphate (AMP) to the dual binding sites in the DNA aptamer, focusing on the dynamic recognition mechanism and pathways. The present results indicate that in addition to the widely known intermolecular interactions, inequivalence of chemical environments of the two binding sites leads to slightly higher stability of AMP binding to the site proximal to the aptamer terminus. In the presence of two AMPs captured by the two sites, each binding free energy is enhanced. In particular, an additional hydrogen bond of AMP to A10 is introduced in the dual-site binding complex, which increases the binding energy from -4.25 ± 0.47 to -9.48 ± 0.33 kcal mol-1 in the site close to the loop. For the dual-site recognition process, the free energy landscape and minimum free energy pathway calculations elucidate the crucial role of electrostatic interactions between the AMP phosphate groups and Na+ ions in positively cooperative binding mechanisms.

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

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

Probing and perturbing riboswitch folding using a fluorescent base analogue

Riboswitches are mRNA segments that regulate gene expression in response to ligand binding. The Class I preQ1 riboswitch consists of a stem-loop and an adenine-rich single-stranded tail ("L3"), which adopt a pseudoknot structure upon binding of the ligand preQ1 . We inserted 2-aminopurine (2-AP), a fluorescent analogue of adenine (A), into the riboswitch at six different positions within L3. Here, 2-AP functions both as a spectroscopic probe and as a "mutation" that reveals how alteration of specific A residues impacts the riboswitch. Using fluorescence and circular dichroism spectroscopy, we found that 2-AP decreases the affinity of the riboswitch for preQ1 at all labeling positions tested, although modified and unmodified variants undergo the same global conformational changes at sufficiently high preQ1 concentration. 2-AP substitution is most detrimental to ligand binding at sites proximal to the ligand-binding pocket, while distal labeling sites exhibit the largest impacts on the stability of the L3 domain in the absence of ligand. Insertion of multiple 2-AP residues does not induce significant additional disruptions. Our results show that interactions involving the A residues in L3 play a critical role in ligand recognition by the preQ1 riboswitch and that 2-AP substitution exerts complex and varied impacts on this riboswitch.

Advances in RNA Labeling with Trifluoromethyl Groups

Fluorine labeling of ribonucleic acids (RNA) in conjunction with 19 F NMR spectroscopy has emerged as a powerful strategy for spectroscopic analysis of RNA structure and dynamics, and RNA-ligand interactions. This study presents the first syntheses of 2'-OCF3 guanosine and uridine phosphoramidites, their incorporation into oligoribonucleotides by solid-phase synthesis and a comprehensive study of their properties. NMR spectroscopic analysis showed that the 2'-OCF3 modification is associated with preferential C2'-endo conformation of the U and G ribose in single-stranded RNA. When paired to the complementary strand, slight destabilization of the duplex caused by the modification was revealed by UV melting curve analysis. Moreover, the power of the 2'-OCF3 label for NMR spectroscopy is demonstrated by dissecting RNA pseudoknot folding and its binding to a small molecule. Furthermore, the 2'-OCF3 modification has potential for applications in therapeutic oligonucleotides. To this end, three 2'-OCF3 modified siRNAs were tested in silencing of the BASP1 gene which indicated enhanced performance for one of them. Importantly, together with earlier work, the present study completes the set of 2'-OCF3 nucleoside phosphoramidites to all four standard nucleobases (A, U, C, G) and hence enables applications that utilize the favorable properties of the 2'-OCF3 group without any restrictions in placing the modification into the RNA target sequence.

Phosphonate and Thiasugar Analogues of Glucosamine-6-phosphate: Activation of the glmS Riboswitch and Antibiotic Activity

The glmS riboswitch is a motif found in 5'-untranslated regions of bacterial mRNA that controls the synthesis of glucosamine-6-phosphate (GlcN6P), an essential building block for the bacterial cell wall, by a feedback mechanism. Activation of the glmS riboswitch by GlcN6P mimics interferes with the ability of bacteria to synthesize its cell wall. Accordingly, GlcN6P mimics acting as glmS activators are promising candidates for future antibiotic drugs that may overcome emerging bacterial resistance against established antibiotics. We describe the synthesis of a series of phosphonate mimics of GlcN6P as well as the thiasugar analogue of GlcN6P. The phosphonate mimics differ in their pKa value to answer the question of whether derivatives with a pKa matching that of GlcN6P would be efficient glmS activators. We found that all derivatives activate the riboswitch, however, less efficiently than GlcN6P. This observation can be explained by the missing hydrogen bonds in the case of phosphonates and is valuable information for the design of future GlcN6P mimics. The thiasugar analogue of GlcN6P on the other hand turned out to be a glmS riboswitch activator with the same activity as the natural metabolite GlcN6P. The nonphosphorylated thiasugar displayed antimicrobial activity against certain bacilli. Therefore, the compound is a promising lead structure for the development of future antibiotics with a potentially novel mode of action.

Magnesium Ion-Driven Folding and Conformational Switching Kinetics of Tetracycline Binding Aptamer: Implications for in vivo Riboswitch Engineering

Engineering in vitro selected RNA aptamers into in vivo functional riboswitches represents a long-standing challenge in molecular biology. The highly specific aptamer domain of the riboswitch undergoes a conformational adjustment in response to ligand sensing, which in turn exerts the regulatory function. Besides essential factors like structural complexity and ligand binding kinetics, the active role of magnesium ions in stabilizing RNA tertiary structures and assisting in ligand binding can be a vital criterion. We present spectroscopic studies on the magnesium ion-driven folding of the Tetracycline binding aptamer. Using fluorescent labels, the aptamer pre-folding and subsequent ligand binding is monitored by magnesium titration experiments and time-resolved stopped-flow measurements. A minimum concentration of 0.5 mM magnesium is required to fold into a magnesium ion-stabilized binding-competent state with a preformed binding pocket. Tetracycline binding causes a pronounced conformational change that results in the establishment of the triple helix core motif, and that further propagates towards the closing stem. By a dynamic acquisition of magnesium ions, a kink motif is formed at the intersection of the triple helix and closing stem regions. This ultimately entails a stabilization of the closing stem which is discussed as a key element in the regulatory function of the Tetracycline aptamer.

Photoswitchable molecular glue for RNA: reversible photocontrol of structure and function of the ribozyme

Single-stranded RNA folds into a variety of secondary and higher-order structures. Distributions and dynamics of multiple RNA conformations are responsible for the biological function of RNA. We here developed a photoswitchable molecular glue for RNA, which could reversibly control the association of two unpaired RNA regions in response to light stimuli. The photoswitchable molecular glue, NCTA, is an RNA-binding ligand possessing a photoisomerizable azobenzene moiety. Z-NCTA is an active ligand for the target RNA containing 5'-WGG-3'/5'-WGG-3' (W = U or A) site and stabilizes its hybridized state, while its isomer E-NCTA is not. Photoreversible isomerization of NCTA enabled control of the secondary and tertiary structure of the target RNA. The RNA-cleaving activity of hammerhead ribozyme, where appropriate RNA folding is necessary, could be reversibly regulated by photoirradiation in cells treated with NCTA, demonstrating precise photocontrol of RNA structure and function by the photoswitchable molecular glue.

Pressure pushes tRNALys3 into excited conformational states

Conformational dynamics play essential roles in RNA function. However, detailed structural characterization of excited states of RNA remains challenging. Here, we apply high hydrostatic pressure (HP) to populate excited conformational states of tRNALys3, and structurally characterize them using a combination of HP 2D-NMR, HP-SAXS (HP-small-angle X-ray scattering), and computational modeling. HP-NMR revealed that pressure disrupts the interactions of the imino protons of the uridine and guanosine U-A and G-C base pairs of tRNALys3. HP-SAXS profiles showed a change in shape, but no change in overall extension of the transfer RNA (tRNA) at HP. Configurations extracted from computational ensemble modeling of HP-SAXS profiles were consistent with the NMR results, exhibiting significant disruptions to the acceptor stem, the anticodon stem, and the D-stem regions at HP. We propose that initiation of reverse transcription of HIV RNA could make use of one or more of these excited states.

RPflex: A Coarse-Grained Network Model for RNA Pocket Flexibility Study

RNA regulates various biological processes, such as gene regulation, RNA splicing, and intracellular signal transduction. RNA’s conformational dynamics play crucial roles in performing its diverse functions. Thus, it is essential to explore the flexibility characteristics of RNA, especially pocket flexibility. Here, we propose a computational approach, RPflex, to analyze pocket flexibility using the coarse-grained network model. We first clustered 3154 pockets into 297 groups by similarity calculation based on the coarse-grained lattice model. Then, we introduced the flexibility score to quantify the flexibility by global pocket features. The results show strong correlations between the flexibility scores and root-mean-square fluctuation (RMSF) values, with Pearson correlation coefficients of 0.60, 0.76, and 0.53 in Testing Sets I–III. Considering both flexibility score and network calculations, the Pearson correlation coefficient was increased to 0.71 in flexible pockets on Testing Set IV. The network calculations reveal that the long-range interaction changes contributed most to flexibility. In addition, the hydrogen bonds in the base–base interactions greatly stabilize the RNA structure, while backbone interactions determine RNA folding. The computational analysis of pocket flexibility could facilitate RNA engineering for biological or medical applications.

Thermodynamics and kinetics of an A-U RNA base pair under force studied by molecular dynamics simulations

Mechanical force has been widely used to study RNA folding and unfolding. Understanding how the force affects the opening and closing of a single base pair, which is a basic step for RNA folding and unfolding and a fundamental behavior in some important biological activities, is crucial to understanding the mechanism of RNA folding and unfolding under mechanical force. In this work, we investigated the opening and closing process of an RNA base pair under mechanical force with constant-force stretching molecular dynamics simulations. It was found that high mechanical force results in overstretching, and the open state is a high-energy state. The enthalpy and entropy change of the base-pair opening-closing transition were obtained and the results at low forces were in good agreement with the nearest-neighbor model. The temperature and force dependence of the opening and closing rates were also obtained. The position of the transition state for the base-pair opening-closing transition under mechanical force was determined. The free energy barrier of opening a base pair without force is the enthalpy increase, and the work done by the force from the closed state to the transition state decreases the barrier and increases the opening rate. The free energy barrier of closing the base pair without force results from the entropy loss, and the work done by the force from the open state to the transition state increases the barrier and decreases the closing rate. The transition rates are strongly dependent on the temperature and force, while the transition path times are weakly dependent on force and temperature.

Structure-based investigations of the NAD+-II riboswitch

Riboswitches are conserved non-coding domains in bacterial mRNA with gene regulation function that are essential for maintaining enzyme co-factor metabolism. Recently, the pnuC RNA motif was reported to selectively bind nicotinamide adenine dinucleotide (NAD+), defining a novel class of NAD+ riboswitches (NAD+-II) according to phylogenetic analysis. To reveal the three-dimensional architecture and the ligand-binding mode of this riboswitch, we solved the crystal structure of NAD+-II riboswitch in complex with NAD+. Strikingly and in contrast to class-I riboswitches that form a tight recognition pocket for the adenosine diphosphate (ADP) moiety of NAD+, the class-II riboswitches form a binding pocket for the nicotinamide mononucleotide (NMN) portion of NAD+ and display only unspecific interactions with the adenosine. We support this finding by an additional structure of the class-II RNA in complex with NMN alone. The structures define a novel RNA tertiary fold that was further confirmed by mutational analysis in combination with isothermal titration calorimetry (ITC), and 2-aminopurine-based fluorescence spectroscopic folding studies. Furthermore, we truncated the pnuC RNA motif to a short RNA helical scaffold with binding affinity comparable to the wild-type motif to allude to the potential of engineering the NAD+-II motif for biotechnological applications.

Structure-based insights into recognition and regulation of SAM-sensing riboswitches

Riboswitches are highly conserved RNA elements that located in the 5'-UTR of mRNAs, which undergo real-time structure conformational change to achieve the regulation of downstream gene expression by sensing their cognate ligands. S-adenosylmethionine (SAM) is a ubiquitous methyl donor for transmethylation reactions in all living organisms. SAM riboswitch is one of the most abundant riboswitches that bind to SAM with high affinity and selectivity, serving as regulatory modules in multiple metabolic pathways. To date, seven SAM-specific riboswitch classes that belong to four families, one SAM/SAH riboswitch and one SAH riboswitch have been identified. Each SAM riboswitch family has a well-organized tertiary core scaffold to support their unique ligand-specific binding pocket. In this review, we summarize the current research progress on the distribution, structure, ligand recognition and gene regulation mechanism of these SAM-related riboswitch families, and further discuss their evolutionary prospects and potential applications.

Probing of Fluorogenic RNA Aptamers via Supramolecular Förster Resonance Energy Transfer with a Universal Fluorescent Nucleobase Analog

Fluorogenic RNA aptamers are synthetic RNAs that have been evolved by in vitro selection methods to bind and light up conditionally fluorescent organic ligands. Compared with other probes for RNA detection, they are less invasive than hybridization-based methods (FISH, molecular beacons) and are considerably smaller than fluorescent protein-recruiting systems (MS2, Pumilio variants). Fluorogenic aptamers have therefore found widespread use as genetically encodable tags for RNA detection in live cells and have also been used in combination with riboswitches to construct versatile metabolite sensors for in vitro use. Their success builds on a fundamental understanding of their three-dimensional structure to explain the mechanisms of ligand interaction and to rationally design functional aptamer devices. In this protocol, we describe a supramolecular FRET-based structure probing method for fluorogenic aptamers that exploits distance- and orientation-dependent energy transfer efficiencies between site-specifically incorporated fluorescent nucleoside analogs and non-covalently bound ligands, exemplified by 4-cyanoindol riboside (4CI) and the DMHBI⁺-binding RNA aptamer Chili. This method yields structural restraints that bridge the gap between traditional low-resolution secondary structure probing methods and more elaborate high-resolution methods such as X-ray crystallography and NMR spectroscopy.

AptaMat: a matrix-based algorithm to compare single-stranded oligonucleotides secondary structures

MOTIVATION

Comparing single-stranded nucleic acids (ssNAs) secondary structures is fundamental when investigating their function and evolution and predicting the effect of mutations on their structures. Many comparison metrics exist, although they are either too elaborate or not sensitive enough to distinguish close ssNAs structures.

RESULTS

In this context, we developed AptaMat, a simple and sensitive algorithm for ssNAs secondary structures comparison based on matrices representing the ssNAs secondary structures and a metric built upon the Manhattan distance in the plane. We applied AptaMat to several examples and compared the results to those obtained by the most frequently used metrics, namely the Hamming distance and the RNAdistance, and by a recently developed image-based approach. We showed that AptaMat is able to discriminate between similar sequences, outperforming all the other here considered metrics. In addition, we showed that AptaMat was able to correctly classify 14 RFAM families within a clustering procedure.

AVAILABILITY AND IMPLEMENTATION

The python code for AptaMat is available at https://github.com/GEC-git/AptaMat.git.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Entropy Driving the Mg2+-Induced Folding of TPP Riboswitch RNA

Mg²⁺ is well known to facilitate the structural folding of RNA. However, the thermodynamic and dynamic roles of Mg²⁺ in RNA folding remain elusive. Here, we exploit single-molecule fluorescence resonance energy transfer (smFRET) and isothermal titration calorimetry (ITC) to study the mechanism of Mg²⁺ in facilitating the folding of thiamine pyrophosphate (TPP) riboswitch RNA. The results of smFRET identify that the presence of Mg²⁺ compacts the RNA and enlarges the conformational dispersity among individual RNA molecules, resulting in a large gain of entropy. The compact yet flexible conformations triggered by Mg²⁺ may help the riboswitch recognize its specific ligand and further fold. This is supported by the ITC experiments, in which the Mg²⁺-induced RNA folding is driven by entropy (ΔS) instead of enthalpy (ΔH). Our results complement the understanding of the Mg²⁺-induced RNA folding. The strategy developed in this work can be used to model other RNAs' folding under different conditions.

The promise of cryo-EM to explore RNA structural dynamics

Conformational dynamics are essential to macromolecular function. This is certainly true of RNA, whose ability to undergo programmed conformational dynamics is essential to create and regulate complex biological processes. However, methods to easily and simultaneously interrogate both the structure and conformational dynamics of fully functional RNAs in isolation and in complex with proteins have not historically been available. Due to its ability to image and classify single particles, cryogenic electron microscopy (cryo-EM) has the potential to address this gap and may be particularly amenable to exploring structural dynamics within the three-dimensional folds of biologically active RNAs. We discuss the possibilities and current limitations of applying cryo-EM to simultaneously study RNA structure and conformational dynamics, and present one example that illustrates this (as of yet) not fully realized potential.

RNA-ligand molecular docking: advances and challenges

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

Analysis of the Passage Times for Unfolding/Folding of the Adenine Riboswitch Aptamer

The conformational transitions of the adenosine deaminase A-riboswitch aptamer both with and without ligand binding are investigated within the tenets of the generalized Langevin equation in a complex viscoelastic cellular environment. Steered molecular dynamics (SMD) simulations are performed to evaluate and compare the results of the first passage times (FPTs) with those obtained from the theory for the unfold and fold transitions of the aptamer. The results of the distribution of Kramers's FPT reveal that the unfold-fold transitions are faster and hence more probable as compared to the fold-unfold transitions of the riboswitch aptamer for both ligand-bound and -unbound states. The transition path time is lower than Kramers's FPT for the riboswitch aptamer as the transition path times for the unfold-fold transition of both without and with ligand binding are insensitive to the details of the exact mechanism of the transition events. However, Kramers's FPTs show varied distributions which correspond to different transition pathways, unlike the transition path times. The mean FPT increases with an increase in the complexity of the cellular environment. The results of Kramers's FPT, transition path time distribution, and mean FPT obtained from our calculations qualitatively match with those obtained from the SMD simulations. Analytically derived values of the mean transition path time show good quantitative agreement with those estimated from the single-molecule force spectroscopy experiments for higher barrier heights.

TPP Riboswitch Populates Holo-Form-like Structure Even in the Absence of Cognate Ligand at High Mg2+ Concentration

Riboswitches are noncoding RNA that regulate gene expression by folding into specific three-dimensional structures (holo-form) upon binding by their cognate ligand in the presence of Mg²⁺. Riboswitch functioning is also hypothesized to be under kinetic control requiring large cognate ligand concentrations. We ask the question under thermodynamic conditions, can the riboswitches populate structures similar to the holo-form only in the presence of Mg²⁺ and absence of cognate ligand binding. We addressed this question using thiamine pyrophosphate (TPP) riboswitch as a model system and computer simulations using a coarse-grained model for RNA. The folding free energy surface (FES) shows that with the initial increase in Mg²⁺ concentration ([Mg²⁺]), the aptamer domain (AD) of TPP riboswitch undergoes a barrierless collapse in its dimensions. On further increase in [Mg²⁺], intermediates separated by barriers appear on the FES, and one of the intermediates has a TPP ligand-binding competent structure. We show that site-specific binding of the Mg²⁺ aids in the formation of tertiary contacts. For [Mg²⁺] greater than physiological concentration, AD folds into a structure similar to the crystal structure of the TPP holo-form even in the absence of the TPP ligand. The folding kinetics shows that TPP AD populates an intermediate due to the misalignment of two arms present in the structure, which acts as a kinetic trap, leading to larger folding timescales. The predictions of the intermediate structures from the simulations are amenable for experimental verification.

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

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

An anionic ligand snap-locks a long-range interaction in a magnesium-folded riboswitch

The archetypical transcriptional crcB fluoride riboswitch from Bacillus cereus is an intricately structured non-coding RNA element enhancing gene expression in response to toxic levels of fluoride. Here, we used single molecule FRET to uncover three dynamically interconverting conformations appearing along the transcription process: two distinct undocked states and one pseudoknotted docked state. We find that the fluoride anion specifically snap-locks the magnesium-induced, dynamically docked state. The long-range, nesting, single base pair A40-U48 acts as the main linchpin, rather than the multiple base pairs comprising the pseudoknot. We observe that the proximally paused RNA polymerase further fine-tunes the free energy to promote riboswitch docking. Finally, we show that fluoride binding at short transcript lengths is an early step toward partitioning folding into the docked conformation. These results reveal how the anionic fluoride ion cooperates with the magnesium-associated RNA to govern regulation of downstream genes needed for fluoride detoxification of the cell.

An integrative NMR-SAXS approach for structural determination of large RNAs defines the substrate-free state of a trans-cleaving Neurospora Varkud Satellite ribozyme

The divide-and-conquer strategy is commonly used for protein structure determination, but its applications to high-resolution structure determination of RNAs have been limited. Here, we introduce an integrative approach based on the divide-and-conquer strategy that was undertaken to determine the solution structure of an RNA model system, the Neurospora VS ribozyme. NMR and SAXS studies were conducted on a minimal trans VS ribozyme as well as several isolated subdomains. A multi-step procedure was used for structure determination that first involved pairing refined NMR structures with SAXS data to obtain structural subensembles of the various subdomains. These subdomain structures were then assembled to build a large set of structural models of the ribozyme, which was subsequently filtered using SAXS data. The resulting NMR-SAXS structural ensemble shares several similarities with the reported crystal structures of the VS ribozyme. However, a local structural difference is observed that affects the global fold by shifting the relative orientation of the two three-way junctions. Thus, this finding highlights a global conformational change associated with substrate binding in the VS ribozyme that is likely critical for its enzymatic activity. Structural studies of other large RNAs should benefit from similar integrative approaches that allow conformational sampling of assembled fragments.

Insights into xanthine riboswitch structure and metal ion-mediated ligand recognition

Riboswitches are conserved functional domains in mRNA that mostly exist in bacteria. They regulate gene expression in response to varying concentrations of metabolites or metal ions. Recently, the NMT1 RNA motif has been identified to selectively bind xanthine and uric acid, respectively, both are involved in the metabolic pathway of purine degradation. Here, we report a crystal structure of this RNA bound to xanthine. Overall, the riboswitch exhibits a rod-like, continuously stacked fold composed of three stems and two internal junctions. The binding-pocket is determined by the highly conserved junctional sequence J1 between stem P1 and P2a, and engages a long-distance Watson–Crick base pair to junction J2. Xanthine inserts between a G–U pair from the major groove side and is sandwiched between base triples. Strikingly, a Mg²⁺ ion is inner-sphere coordinated to O6 of xanthine and a non-bridging oxygen of a backbone phosphate. Two further hydrated Mg²⁺ ions participate in extensive interactions between xanthine and the pocket. Our structure model is verified by ligand binding analysis to selected riboswitch mutants using isothermal titration calorimetry, and by fluorescence spectroscopic analysis of RNA folding using 2-aminopurine-modified variants. Together, our study highlights the principles of metal ion-mediated ligand recognition by the xanthine riboswitch.

Microsecond Folding of preQ1 Riboswitch and Its Biological Significance Revealed by Two-Dimensional Fluorescence Lifetime Correlation Spectroscopy

Riboswitches are regulatory elements of bacterial mRNA which function with conformational switching upon binding of specific cellular metabolites. In particular, transcriptional riboswitches regulate gene expression kinetically through the conformational change of the aptamer domain. In this study, we investigate the conformational dynamics and ligand binding mechanisms of the aptamer domain of a transcriptional prequeuosine (preQ₁) riboswitch from Bacillus subtilis using two-dimensional fluorescence lifetime correlation spectroscopy (2D FLCS) with microsecond time resolution. The obtained time-resolved single-molecule data indicate that the aptamer domain undergoes folding/unfolding including three forms, which are attributed to hairpin (O), pseudoknot-like (pF), and H-type pseudoknot (fF) structures. It is found that a cofactor, Mg²⁺, binds only to the fF form with the conformational selection mechanism. In contrast, it is indicated that the ligand, preQ₁, binds to the O form with the induced-fit mechanism and significantly accelerates the microsecond O → pF folding process. It is also shown that the binding with preQ₁ substantially stabilizes the fF form that is generated from the pF form with a long time constant (>10 ms). Combining these results with the results of a former smFRET study on the slower time scale, we obtain an overall picture of the folding/unfolding dynamics of the aptamer domain as well as its energy landscape. On the basis of the picture obtained, we discuss the significance of the microsecond folding/unfolding of the aptamer domain for biological function of the riboswitch and propose the molecular mechanism of the gene expression controlled by the structural dynamics of the aptamer domain.

What defines a synthetic riboswitch? - Conformational dynamics of ciprofloxacin aptamers with similar binding affinities but varying regulatory potentials

Among the many in vitro-selected aptamers derived from SELEX protocols, only a small fraction has the potential to be applied for synthetic riboswitch engineering. Here, we present a comparative study of the binding properties of three different aptamers that bind to ciprofloxacin with similar K_D values, yet only two of them can be applied as riboswitches. We used the inherent ligand fluorescence that is quenched upon binding as the reporter signal in fluorescence titration and in time-resolved stopped-flow experiments. Thus, we were able to demonstrate differences in the binding kinetics of regulating and non-regulating aptamers. All aptamers studied underwent a two-step binding mechanism that suggests an initial association step followed by a reorganization of the aptamer to accommodate the ligand. We show that increasing regulatory potential is correlated with a decreasing back-reaction rate of the second binding step, thus resulting in a virtually irreversible last binding step of regulating aptamers. We suggest that a highly favoured structural adaption of the RNA to the ligand during the final binding step is essential for turning an aptamer into a riboswitch. In addition, our results provide an explanation for the fact that so few aptamers with regulating capacity have been found to date. Based on our data, we propose an adjustment of the selection protocol for efficient riboswitch detection.

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.

In-Cell NMR Spectroscopy of Functional Riboswitch Aptamers in Eukaryotic Cells

We report here the in-cell NMR-spectroscopic observation of the binding of the cognate ligand 2'-deoxyguanosine to the aptamer domain of the bacterial 2'-deoxyguanosine-sensing riboswitch in eukaryotic cells, namely Xenopus laevis oocytes and in human HeLa cells. The riboswitch is sufficiently stable in both cell types to allow for detection of binding of the ligand to the riboswitch. Most importantly, we show that the binding mode established by in vitro characterization of this prokaryotic riboswitch is maintained in eukaryotic cellular environment. Our data also bring important methodological insights: Thus far, in-cell NMR studies on RNA in mammalian cells have been limited to investigations of short (<15 nt) RNA fragments that were extensively modified by protecting groups to limit their degradation in the intracellular space. Here, we show that the in-cell NMR setup can be adjusted for characterization of much larger (≈70 nt) functional and chemically non-modified RNA.

A second riboswitch class for the enzyme cofactor NAD. RNA (NEW YORK, N.Y.) 2021;27:99-105. [PMID: 33087526 PMCID: PMC7749635 DOI: 10.1261/rna.077891.120]

A bacterial noncoding RNA motif almost exclusively associated with pnuC genes was uncovered using comparative sequence analysis. Some PnuC proteins are known to transport nicotinamide riboside (NR), which is a component of the ubiquitous and abundant enzyme cofactor nicotinamide adenine dinucleotide (NAD⁺). Thus, we speculated that the newly found "pnuC motif" RNAs might function as aptamers for a novel class of NAD⁺-sensing riboswitches. RNA constructs that encompass the conserved nucleotides and secondary structure features that define the motif indeed selectively bind NAD⁺, nicotinamide mononucleotide (NMN), and NR. Mutations that disrupt strictly conserved nucleotides of the aptamer also disrupt ligand binding. These bioinformatic and biochemical findings indicate that pnuC motif RNAs are likely members of a second riboswitch class that regulates gene expression in response to NAD⁺ binding.

Structural distinctions between NAD+ riboswitch domains 1 and 2 determine differential folding and ligand binding

Riboswitches are important gene regulatory elements frequently encountered in bacterial mRNAs. The recently discovered nadA riboswitch contains two similar, tandemly arrayed aptamer domains, with the first domain possessing high affinity for nicotinamide adenine dinucleotide (NAD⁺). The second domain which comprises the ribosomal binding site in a putative regulatory helix, however, has withdrawn from detection of ligand-induced structural modulation thus far, and therefore, the identity of the cognate ligand and the regulation mechanism have remained unclear. Here, we report crystal structures of both riboswitch domains, each bound to NAD⁺. Furthermore, we demonstrate that ligand binding to domain 2 requires significantly higher concentrations of NAD⁺ (or ADP retaining analogs) compared to domain 1. Using a fluorescence spectroscopic approach, we further shed light on the structural features which are responsible for the different ligand affinities, and describe the Mg²⁺-dependent, distinct folding and pre-organization of their binding pockets. Finally, we speculate about possible scenarios for nadA RNA gene regulation as a putative two-concentration sensor module for a time-controlled signal that is primed and stalled by the gene regulation machinery at low ligand concentrations (domain 1), and finally triggers repression of translation as soon as high ligand concentrations are reached in the cell (domain 2).

Cooperativity and Allostery in RNA Systems

Allostery is among the most basic biological principles employed by biological macromolecules to achieve a biologically active state in response to chemical cues. Although initially used to describe the impact of small molecules on the conformation and activity of protein enzymes, the definition of this term has been significantly broadened to describe long-range conformational change of macromolecules in response to small or large effectors. Such a broad definition could be applied to RNA molecules, which do not typically serve as protein-free cellular enzymes but fold and form macromolecular assemblies with the help of various ligand molecules, including ions and proteins. Ligand-induced allosteric changes in RNA molecules are often accompanied by cooperative interactions between RNA and its ligand, thus streamlining the folding and assembly pathways. This chapter provides an overview of the interplay between cooperativity and allostery in RNA systems and outlines methods to study these two biological principles.

Aromatic 19 F-13 C TROSY-[19 F, 13 C]-Pyrimidine Labeling for NMR Spectroscopy of RNA

We present the access to [5-19 F, 5-13 C]-uridine and -cytidine phosphoramidites for the production of site-specifically modified RNAs up to 65 nucleotides (nts). The amidites were used to introduce [5-19 F, 5-13 C]-pyrimidine labels into five RNAs-the 30 nt human immunodeficiency virus trans activation response (HIV TAR) 2 RNA, the 61 nt human hepatitis B virus ϵ (hHBV ϵ) RNA, the 49 nt SAM VI riboswitch aptamer domain from B. angulatum, the 29 nt apical stem loop of the pre-microRNA (miRNA) 21 and the 59 nt full length pre-miRNA 21. The main stimulus to introduce the aromatic 19 F-13 C-spin topology into RNA comes from a work of Boeszoermenyi et al., in which the dipole-dipole interaction and the chemical shift anisotropy relaxation mechanisms cancel each other leading to advantageous TROSY properties shown for aromatic protein sidechains. This aromatic 13 C-19 F labeling scheme is now transferred to RNA. We provide a protocol for the resonance assignment by solid phase synthesis based on diluted [5-19 F, 5-13 C]/[5-19 F] pyrimidine labeling. For the 61 nt hHBV ϵ we find a beneficial 19 F-13 C TROSY enhancement, which should be even more pronounced in larger RNAs and will facilitate the NMR studies of larger RNAs. The [19 F, 13 C]-labeling of the SAM VI aptamer domain and the pre-miRNA 21 further opens the possibility to use the biorthogonal stable isotope reporter nuclei in in vivo NMR to observe ligand binding and microRNA processing in a biological relevant setting.

Synthesis of Nucleobase-Modified RNA Oligonucleotides by Post-Synthetic Approach

The chemical synthesis of modified oligoribonucleotides represents a powerful approach to study the structure, stability, and biological activity of RNAs. Selected RNA modifications have been proven to enhance the drug-like properties of RNA oligomers providing the oligonucleotide-based therapeutic agents in the antisense and siRNA technologies. The important sites of RNA modification/functionalization are the nucleobase residues. Standard phosphoramidite RNA chemistry allows the site-specific incorporation of a large number of functional groups to the nucleobase structure if the building blocks are synthetically obtainable and stable under the conditions of oligonucleotide chemistry and work-up. Otherwise, the chemically modified RNAs are produced by post-synthetic oligoribonucleotide functionalization. This review highlights the post-synthetic RNA modification approach as a convenient and valuable method to introduce a wide variety of nucleobase modifications, including recently discovered native hypermodified functional groups, fluorescent dyes, photoreactive groups, disulfide crosslinks, and nitroxide spin labels.

Supramolecular Fluorescence Resonance Energy Transfer in Nucleobase-Modified Fluorogenic RNA Aptamers

RNA aptamers form compact tertiary structures and bind their ligands in specific binding sites. Fluorescence-based strategies reveal information on structure and dynamics of RNA aptamers. Herein, we report the incorporation of the universal emissive nucleobase analog 4-cyanoindole into the fluorogenic RNA aptamer Chili, and its application as a donor for supramolecular FRET to the bound ligands DMHBI+ or DMHBO+ . The photophysical properties of the new nucleobase-ligand-FRET pair revealed structural restraints for the overall RNA aptamer organization and identified nucleotide positions suitable for FRET-based readout of ligand binding. This strategy is generally suitable for binding-site mapping and may also be applied for responsive aptamer devices.

Homogeneous and Universal Detection of Various Targets with a Dual-Step Transduced Toehold Switch Sensor

Toehold switch sensors represent a class of new advances that allow specific targets to trigger in situ expression of a protein reporter. Although they offer unique advantages of a label-free nature and high portability, they generally require repeated sequence design, high expenditure, and laborious optimization of toehold switch sequences according to different targets. To simplify the sensing process further and to improve its practicability, we innovatively introduce a dual-step pre-transduction upon traditional toehold switch sensor. Through two successive toehold-mediated strand-displacement reactions that are initiated, respectively, by a linear and an associative trigger, different DNAs, RNAs, or ligands of functional nucleic acids can be generally transduced into the input of one high-performance toehold switch sensor. This advance greatly increases the target range. Furthermore, the whole process is signal-on, homogeneous, and free of any requirements for complicated operations such as probe labeling, separation, and reconstruction of the toehold switch, being promising and practical even in portable or point-of-care assays.

Alternate RNA Structures

RNA molecules fold into complex three-dimensional structures that sample alternate conformations ranging from minor differences in tertiary structure dynamics to major differences in secondary structure. This allows them to form entirely different substructures with each population potentially giving rise to a distinct biological outcome. The substructures can be partitioned along an existing energy landscape given a particular static cellular cue or can be shifted in response to dynamic cues such as ligand binding. We review a few key examples of RNA molecules that sample alternate conformations and how these are capitalized on for control of critical regulatory functions.

SAM-VI riboswitch structure and signature for ligand discrimination

Riboswitches are metabolite-sensing, conserved domains located in non-coding regions of mRNA that are central to regulation of gene expression. Here we report the first three-dimensional structure of the recently discovered S-adenosyl-L-methionine responsive SAM-VI riboswitch. SAM-VI adopts a unique fold and ligand pocket that are distinct from all other known SAM riboswitch classes. The ligand binds to the junctional region with its adenine tightly intercalated and Hoogsteen base-paired. Furthermore, we reveal the ligand discrimination mode of SAM-VI by additional X-ray structures of this riboswitch bound to S-adenosyl-L-homocysteine and a synthetic ligand mimic, in combination with isothermal titration calorimetry and fluorescence spectroscopy to explore binding thermodynamics and kinetics. The structure is further evaluated by analysis of ligand binding to SAM-VI mutants. It thus provides a thorough basis for developing synthetic SAM cofactors for applications in chemical and synthetic RNA biology.

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

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

Programmable Artificial Cells Using Histamine-Responsive Synthetic Riboswitch

Artificial cells that encapsulate DNA-programmable protein expression machinery are emerging as an attractive platform for studying fundamental cellular properties and applications in synthetic biology. However, interfacing these artificial cells with the complex and dynamic chemical environment remains a major and urgent challenge. We demonstrate that the repertoire of molecules that artificial cells respond to can be expanded by synthetic RNA-based gene switches, or riboswitches. We isolated an RNA aptamer that binds histamine with high affinity and specificity and used it to design robust riboswitches that activate protein expression in the presence of histamine. Finally, the riboswitches were incorporated in artificial cells to achieve controlled release of an encapsulated small molecule and to implement a self-destructive kill-switch. Synthetic riboswitches should serve as modular and versatile interfaces to link artificial cell phenotypes with the complex chemical environment.

Understanding the Contributions of Conformational Changes, Thermodynamics, and Kinetics of RNA-Small Molecule Interactions

The implication of RNA in multiple cellular processes beyond protein coding has revitalized interest in the development of small molecules for therapeutically targeting RNA and for further probing its cellular biology. However, the process of rationally designing such small molecule probes is hampered by the paucity of information about fundamental molecular recognition principles of RNA. In this Review, we summarize two important and often underappreciated aspects of RNA-small molecule recognition: RNA conformational dynamics and the biophysical properties of interactions of small molecules with RNA, specifically thermodynamics and kinetics. While conformational flexibility is often said to impede RNA ligand development, the ability of small molecules to influence the RNA conformational landscape can have a significant effect on the cellular functions of RNA. An analysis of the conformational landscape of RNA and the interactions of individual conformations with ligands can thus guide the development of new small molecule probes, which needs to be investigated further. Additionally, while it is common practice to quantify the binding affinities ( K_a or K_d) of small molecules for biomacromolecules as a measure of their activity, further biophysical characterization of their interaction can provide a deeper understanding. Studies that focus on the thermodynamic and kinetic parameters for interaction between RNA and ligands are next discussed. Finally, this Review provides the reader with a perspective on how such in-depth analysis of biophysical characteristics of the interaction of RNA and small molecules can impact our understanding of these interactions and how they will benefit the future design of small molecule probes.

Enzymatic synthesis of base-modified RNA by T7 RNA polymerase. A systematic study and comparison of 5-substituted pyrimidine and 7-substituted 7-deazapurine nucleoside triphosphates as substrates

We synthesized a small library of eighteen 5-substituted pyrimidine or 7-substituted 7-deazapurine nucleoside triphosphates bearing methyl, ethynyl, phenyl, benzofuryl or dibenzofuryl groups through cross-coupling reactions of nucleosides followed by triphosphorylation or through direct cross-coupling reactions of halogenated nucleoside triphosphates. We systematically studied the influence of the modification on the efficiency of T7 RNA polymerase catalyzed synthesis of modified RNA and found that modified ATP, UTP and CTP analogues bearing smaller modifications were good substrates and building blocks for the RNA synthesis even in difficult sequences incorporating multiple modified nucleotides. Bulky dibenzofuryl derivatives of ATP and GTP were not substrates for the RNA polymerase. In the case of modified GTP analogues, a modified procedure using a special promoter and GMP as initiator needed to be used to obtain efficient RNA synthesis. The T7 RNA polymerase synthesis of modified RNA can be very efficiently used for synthesis of modified RNA but the method has constraints in the sequence of the first three nucleotides of the transcript, which must contain a non-modified G in the +1 position.

Single-Molecule FRET Kinetics of the Mn2+ Riboswitch: Evidence for Allosteric Mg2+ Control of "Induced-Fit" vs "Conformational Selection" Folding Pathways

Gene expression in bacteria is often regulated dynamically by conformational changes in a riboswitch upon ligand binding, a detailed understanding of which is very much in its infancy. For example, the manganese riboswitch is a widespread RNA motif that conformationally responds in regulating bacterial gene expression to micromolar levels of its eponymous ligand, Mn²⁺, but the mechanistic pathways are poorly understood. In this work, we quantitatively explore the dynamic folding behavior of the manganese riboswitch by single-molecule fluorescence resonance energy transfer spectroscopy as a function of cation/ligand conditions. From the detailed analysis of the kinetics, the Mn²⁺ is shown to fold the riboswitch by a "bind-then-fold" (i.e., "induced-fit", IF) mechanism, whereby the ligand binds first and then promotes folding. On the other hand, the data also clearly reveal the presence of a folded yet ligand-free structure predominating due to the addition of physiological Mg²⁺ to a nonselective metal ion binding site. Of particular kinetic interest, such a Mg²⁺ "prefolded" conformation of the riboswitch is shown to exhibit a significantly increased affinity for Mn²⁺ and further stabilization by subsequent binding of the ligand, thereby promoting efficient riboswitch folding by a "fold-then-bind" (i.e., "conformational selection", CS) mechanism. Our results not only demonstrate Mg²⁺-controlled switching between IF and CS riboswitch folding pathways but also suggest a novel heterotropic allosteric control in the manganese riboswitch activity co-regulated by Mg²⁺ binding.

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.

Computational Investigation of RNA A-Bulges Related to the Microtubule-Associated Protein Tau Causing Frontotemporal Dementia and Parkinsonism

Mutations in the human tau gene result in alternative splicing of the tau protein, which causes frontotemporal dementia and Parkinsonism. One disease mechanism is linked to the stability of a hairpin within the microtubule-associated protein tau (MAPT) mRNA, which contains an A-bulge. Here we employ computational methods to investigate the structural and thermodynamic properties of several A-bulge RNAs with different closing base-pairs. We find that the current amber RNA force field has a preference to overstabilize base-triple over stacked states, even though some of the A-bulges are known to prefer stacked states according to NMR studies. We further determined that if the neighboring base-pairs of A-bulges are AU, this situation can lead to base slippage. However, when the 3'-side of the A-bulge has an UA base-pair, the stacked state is stabilized by an extra interaction that is not observed in the other sequences. We suggest that these A-bulge RNA systems could be used as benchmarks to improve the current RNA force fields.