Adaptive ligand binding by the purine riboswitch in the recognition of guanine and adenine analogs

Bioinformatic analysis of riboswitch structures uncovers variant classes with altered ligand specificity

Riboswitches are RNAs that form complex, folded structures that selectively bind small molecules or ions. As with certain groups of protein enzymes and receptors, some riboswitch classes have evolved to change their ligand specificity. We developed a procedure to systematically analyze known riboswitch classes to find additional variants that have altered their ligand specificity. This approach uses multiple-sequence alignments, atomic-resolution structural information, and riboswitch gene associations. Among the discoveries are unique variants of the guanine riboswitch class that most tightly bind the nucleoside 2'-deoxyguanosine. In addition, we identified variants of the glycine riboswitch class that no longer recognize this amino acid, additional members of a rare flavin mononucleotide (FMN) variant class, and also variants of c-di-GMP-I and -II riboswitches that might recognize different bacterial signaling molecules. These findings further reveal the diverse molecular sensing capabilities of RNA, which highlights the potential for discovering a large number of additional natural riboswitch classes.

The Emerging Role of RNA as a Therapeutic Target for Small Molecules

Recent advances in understanding different RNAs and unique features of their biology have revealed a wealth of information. However, approaches to identify small molecules that target these newly discovered regulatory elements have been lacking. The application of new biochemical screening and design-based technologies, coupled with a resurgence of interest in phenotypic screening, has resulted in several compelling successes in targeting RNA. A number of recent advances suggest that achieving the long-standing goal of developing drug-like, biologically active small molecules that target RNA is possible. This review highlights advances and successes in approaches to targeting RNA with diverse small molecules, and the potential for these technologies to pave the way to new types of RNA-targeted therapeutics.

Small molecules targeting viral RNA

Highly conserved noncoding RNA (ncRNA) elements in viral genomes and transcripts offer new opportunities to expand the repertoire of drug targets for the development of antiinfective therapy. Ligands binding to ncRNA architectures are able to affect interactions, structural stability or conformational changes and thereby block processes essential for viral replication. Proof of concept for targeting functional RNA by small molecule inhibitors has been demonstrated for multiple viruses with RNA genomes. Strategies to identify antiviral compounds as inhibitors of ncRNA are increasingly emphasizing consideration of drug‐like properties of candidate molecules emerging from screening and ligand design. Recent efforts of antiviral lead discovery for RNA targets have provided drug‐like small molecules that inhibit viral replication and include inhibitors of human immunodeficiency virus (HIV), hepatitis C virus (HCV), severe respiratory syndrome coronavirus (SARS CoV), and influenza A virus. While target selectivity remains a challenge for the discovery of useful RNA‐binding compounds, a better understanding is emerging of properties that define RNA targets amenable for inhibition by small molecule ligands. Insight from successful approaches of targeting viral ncRNA in HIV, HCV, SARS CoV, and influenza A will provide a basis for the future exploration of RNA targets for therapeutic intervention in other viral pathogens which create urgent, unmet medical needs. Viruses for which targeting ncRNA components in the genome or transcripts may be promising include insect‐borne flaviviruses (Dengue, Zika, and West Nile) and filoviruses (Ebola and Marburg). WIREs RNA 2016, 7:726–743. doi: 10.1002/wrna.1373

This article is categorized under: 1

RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems

2

RNA Interactions with Proteins and Other Molecules > Small Molecule–RNA Interactions

3

Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs

Emerging applications of riboswitches - from antibacterial targets to molecular tools

The ability to precisely regulate gene expression is one of the most important features of the living cells as it enables the adaptation and survival in different environmental conditions. The majority of regulatory mechanisms involve protein action, however, multiple genes are controlled by nucleic acids. Among RNA-based regulators, the riboswitches present a large group of specific domains within messenger RNAs able to respond to small metabolites, tRNA, secondary messengers, ions, vitamins or amino acids. A simple, accurate, and efficient mechanism of action as well as easiness in handling and engineering make the riboswitches a potent practical tool in industry, medicine, pharmacy or environmental protection. Hereby, we summarize the current achievements and challenges in designing and practical employment of the riboswitch-based tools.

Watson-Crick-like pairs in CCUG repeats: evidence for tautomeric shifts or protonation

RNA transcripts that include expanded CCUG repeats are associated with myotonic dystrophy type 2. Crystal structures of two CCUG-containing oligomers show that the RNA strands associate into slipped duplexes that contain noncanonical C-U pairs that have apparently undergone tautomeric transition or protonation resulting in an unusual Watson-Crick-like pairing. The overhanging ends of the duplexes interact forming U-U pairs, which also show tautomerism. Duplexes consisting of CCUG repeats are thermodynamically less stable than the trinucleotide repeats involved in the TRED genetic disorders, but introducing LNA residues increases their stability and raises the melting temperature of the studied oligomers by ∼10°C, allowing detailed crystallographic studies. Quantum mechanical calculations were performed to test the possibility of the tautomeric transitions or protonation within the noncanonical pairs. The results indicate that tautomeric or ionic shifts of nucleobases can manifest themselves in biological systems, supplementing the canonical "rules of engagement."

RNA approaches the B-form in stacked single strand dinucleotide contexts

Duplex RNA adopts an A-form structure, while duplex DNA interconverts between the A- and B-forms depending on the environment. The C2'-endo sugar pucker seen in B-form DNA can occur infrequently in ribose sugars as well, but RNA is not understood to assume B-form conformations. Through analysis of over 45,000 stacked single strand dinucleotide (SSD) crystal structure conformations, this study demonstrates that RNA is capable of adopting a wide conformational range between the canonical A- and B-forms at the localized SSD level, including many B-form-like conformations. It does so through C2'-endo ribose conformations in one or both nucleotides, and B-form-like neighboring base stacking patterns. As chemical reactions on nucleic acids involve localized changes in chemical bonds, the understanding of how enzymes distinguish between DNA and RNA nucleotides is altered by the energetic accessibility of these rare B-form-like RNA SSD conformations. The existence of these conformations also has direct implications in parametrization of molecular mechanics energy functions used extensively to model nucleic acid behavior., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 65-82, 2016.

Two Divalent Metal Ions and Conformational Changes Play Roles in the Hammerhead Ribozyme Cleavage Reaction

The hammerhead ribozyme is a self-cleaving RNA broadly dispersed across all kingdoms of life. Although it was the first of the small, nucleolytic ribozymes discovered, the mechanism by which it catalyzes its reaction remains elusive. The nucleobase of G12 is well positioned to be a general base, but it is unclear if or how this guanine base becomes activated for proton transfer. Metal ions have been implicated in the chemical mechanism, but no interactions between divalent metal ions and the cleavage site have been observed crystallographically. To better understand how this ribozyme functions, we have solved crystal structures of wild-type and G12A mutant ribozymes. We observe a pH-dependent conformational change centered around G12, consistent with this nucleotide becoming deprotonated. Crystallographic and kinetic analysis of the G12A mutant reveals a Zn(2+) specificity switch suggesting a direct interaction between a divalent metal ion and the purine at position 12. The metal ion specificity switch and the pH-rate profile of the G12A mutant suggest that the minor imino tautomer of A12 serves as the general base in the mutant ribozyme. We propose a model in which the hammerhead ribozyme rearranges prior to the cleavage reaction, positioning two divalent metal ions in the process. The first metal ion, positioned near G12, becomes directly coordinated to the O6 keto oxygen, to lower the pKa of the general base and organize the active site. The second metal ion, positioned near G10.1, bridges the N7 of G10.1 and the scissile phosphate and may participate directly in the cleavage reaction.

Blu-ray based optomagnetic aptasensor for detection of small molecules

This paper describes an aptamer-based optomagnetic biosensor for detection of a small molecule based on target binding-induced inhibition of magnetic nanoparticle (MNP) clustering. For the detection of a target small molecule, two mutually exclusive binding reactions (aptamer-target binding and aptamer-DNA linker hybridization) are designed. An aptamer specific to the target and a DNA linker complementary to a part of the aptamer sequence are immobilized onto separate MNPs. Hybridization of the DNA linker and the aptamer induces formation of MNP clusters. The target-to-aptamer binding on MNPs prior to the addition of linker-functionalized MNPs significantly hinders the hybridization reaction, thus reducing the degree of MNP clustering. The clustering state, which is thus related to the target concentration, is then quantitatively determined by an optomagnetic readout technique that provides the hydrodynamic size distribution of MNPs and their clusters. A commercial Blu-ray optical pickup unit is used for optical signal acquisition, which enables the establishment of a low-cost and miniaturized biosensing platform. Experimental results show that the degree of MNP clustering correlates well with the concentration of a target small molecule, adenosine triphosphate (ATP) in this work, in the range between 10µM and 10mM. This successful proof-of-concept indicates that our optomagnetic aptasensor can be further developed as a low-cost biosensing platform for detection of small molecule biomarkers in an out-of-lab setting.

Novel riboswitch-binding flavin analog that protects mice against Clostridium difficile infection without inhibiting cecal flora

Novel mechanisms of action and new chemical scaffolds are needed to rejuvenate antibacterial drug discovery, and riboswitch regulators of bacterial gene expression are a promising class of targets for the discovery of new leads. Herein, we report the characterization of 5-(3-(4-fluorophenyl)butyl)-7,8-dimethylpyrido[3,4-b]quinoxaline-1,3(2H,5H)-dione (5FDQD)-an analog of riboflavin that was designed to bind riboswitches that naturally recognize the essential coenzyme flavin mononucleotide (FMN) and regulate FMN and riboflavin homeostasis. In vitro, 5FDQD and FMN bind to and trigger the function of an FMN riboswitch with equipotent activity. MIC and time-kill studies demonstrated that 5FDQD has potent and rapidly bactericidal activity against Clostridium difficile. In C57BL/6 mice, 5FDQD completely prevented the onset of lethal antibiotic-induced C. difficile infection (CDI). Against a panel of bacteria representative of healthy bowel flora, the antibacterial selectivity of 5FDQD was superior to currently marketed CDI therapeutics, with very little activity against representative strains from the Bacteroides, Lactobacillus, Bifidobacterium, Actinomyces, and Prevotella genera. Accordingly, a single oral dose of 5FDQD caused less alteration of culturable cecal flora in mice than the comparators. Collectively, these data suggest that 5FDQD or closely related analogs could potentially provide a high rate of CDI cure with a low likelihood of infection recurrence. Future studies will seek to assess the role of FMN riboswitch binding to the mechanism of 5FDQD antibacterial action. In aggregate, our results indicate that riboswitch-binding antibacterial compounds can be discovered and optimized to exhibit activity profiles that merit preclinical and clinical development as potential antibacterial therapeutic agents.

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.

Visualizing transient Watson-Crick-like mispairs in DNA and RNA duplexes

Rare tautomeric and anionic nucleobases are believed to have fundamental biological roles, but their prevalence and functional importance has remained elusive because they exist transiently, in low abundance, and involve subtle movements of protons that are difficult to visualize. Using NMR relaxation dispersion, we show here that wobble dG•dT and rG•rU mispairs in DNA and RNA duplexes exist in dynamic equilibrium with short-lived, low-populated Watson-Crick-like mispairs that are stabilized by rare enolic or anionic bases. These mispairs can evade Watson-Crick fidelity checkpoints and form with probabilities (10(-3) to 10(-5)) that strongly imply a universal role in replication and translation errors. Our results indicate that rare tautomeric and anionic bases are widespread in nucleic acids, expanding their structural and functional complexity beyond that attainable with canonical bases.

Role of tautomerism in RNA biochemistry

Heterocyclic nucleic acid bases and their analogs can adopt multiple tautomeric forms due to the presence of multiple solvent-exchangeable protons. In DNA, spontaneous formation of minor tautomers has been speculated to contribute to mutagenic mispairings during DNA replication, whereas in RNA, minor tautomeric forms have been proposed to enhance the structural and functional diversity of RNA enzymes and aptamers. This review summarizes the role of tautomerism in RNA biochemistry, specifically focusing on the role of tautomerism in catalysis of small self-cleaving ribozymes and recognition of ligand analogs by riboswitches. Considering that the presence of multiple tautomers of nucleic acid bases is a rare occurrence, and that tautomers typically interconvert on a fast time scale, methods for studying rapid tautomerism in the context of nucleic acids under biologically relevant aqueous conditions are also discussed.

Structural stability, acidity, and halide selectivity of the fluoride riboswitch recognition site

Using static and dynamics DFT methods we show that the Mg(2+)/F(-)/phosphate/water cluster at the center of the fluoride riboswitch is stable by its own and, once assembled, does not rely on any additional factor from the overall RNA fold. Further, we predict that the pKa of the water molecule bridging two Mg cations is around 8.4. We also demonstrate that the halide selectivity of the fluoride riboswitch is determined by the stronger Mg-F bond, which is capable of keeping together the cluster. Replacing F(-) with Cl(-) results in a cluster that is unstable under dynamic conditions. Similar conclusions on the structure and energetics of the cluster in the binding pocket of fluoride-inhibited pyrophosphatase suggest that the peculiarity of fluoride is in its ability to establish much stronger metal-halide bonds.

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.

Halogen bonds in crystal TTF derivatives: an ab initio quantum mechanical study

The stabilisation energies of five ionic and neutral organic crystal structures containing various halogen bonds (I···I, Br···Br, I···Br, I···S and Br···S) were calculated using the DFT-D3 method (B97D/def2-QZVP). Besides them, the ionic I3(-)···I2 and neutral I2···I2, complexes (in the crystal geometries) were also studied. The nature of the bonds was deduced from the electrostatic potential evaluated for all subsystems. In almost all the cases, the σ-hole was positive; it was negative only for the ionic I3(-) system (although more positive than the respective belt value). The strongest halogen bonds were those that involved iodine as a halogen-bond donor and acceptor. Among ionic X···I3(-) and neutral X···I2 and X···Y dimers, the neutral X···I2 complexes were, surprisingly enough, the most stable; the highest stabilisation energy of 13.8 kcal mol(-1) was found for the I2···1,3-dithiole-2-thione-4-carboxylic acid complex. The stabilisation energies of the ionic I3(-)···I2 and neutral I2···1,3-dithiole-2-thione-4-carboxylic acid (20.2 and 20.42 kcal mol(-1), respectively) complexes are very high, which is explained by the favourable geometrical arrangement, allowing the formation of a strong halogen bond. An I···I halogen bond also exists in the neutral I2···I2 complex, having only moderate stabilisation energy (3.9 kcal mol(-1)). This stabilisation energy was, however, shown to be close to that in the optimal gas-phase L-shaped I2···I2 complex. In all the cases, the dispersion energy is important and comparable to electrostatic energy. Only in strong halogen bonds (e.g. I3(-)···I2), the electrostatic energy becomes dominant.

Binding site preorganization and ligand discrimination in the purine riboswitch

The progress of RNA research has suggested a wide variety of RNA molecules as possible targets for pharmaceutical drug molecules. Structure-based computational methods for predicting binding modes and affinities are now important tools in drug discovery, but these methods have mainly been focused on protein targets. Here we employ molecular dynamics free-energy perturbation calculations and the linear interaction energy method to analyze the energetics of ligand binding to purine riboswitches. Calculations are carried out for 14 different purine complexes with the guanine and adenine riboswitches in order to examine their ligand recognition principles. The simulations yield binding affinities in good agreement with experimental data and rationalize the selectivity of the riboswitches for different ligands. In particular, it is found that these receptors have an unusually high degree of electrostatic preorganization for their cognate ligands, and this effect is further quantified by explicit free-energy calculations, which show that the standard electrostatic linear interaction energy parametrization is suboptimal in this case. The adenine riboswitch specifically uses the electrostatic preorganization to discriminate against guanine by preventing the formation of a G-U wobble base pair.

Modular riboswitch toolsets for synthetic genetic control in diverse bacterial species

Ligand-dependent control of gene expression is essential for gene functional analysis, target validation, protein production, and metabolic engineering. However, the expression tools currently available are difficult to transfer between species and exhibit limited mechanistic diversity. Here we demonstrate how the modular architecture of purine riboswitches can be exploited to develop orthogonal and chimeric switches that are transferable across diverse bacterial species, modulating either transcription or translation, to provide tunable activation or repression of target gene expression, in response to synthetic non-natural effector molecules. Our novel riboswitch-ligand pairings are shown to regulate physiologically important genes required for bacterial motility in Escherichia coli and cell morphology in Bacillus subtilis. These findings are relevant for future gene function studies and antimicrobial target validation, while providing new modular and orthogonal regulatory components for deployment in synthetic biology regimes.

Recognition of Watson-Crick base pairs: constraints and limits due to geometric selection and tautomerism

The natural bases of nucleic acids have a strong preference for one tautomer form, guaranteeing fidelity in their hydrogen bonding potential. However, base pairs observed in recent crystal structures of polymerases and ribosomes are best explained by an alternative base tautomer, leading to the formation of base pairs with Watson-Crick-like geometries. These observations set limits to geometric selection in molecular recognition of complementary Watson-Crick pairs for fidelity in replication and translation processes.

Direct observation of multiple tautomers of oxythiamine and their recognition by the thiamine pyrophosphate riboswitch

Structural diversification of canonical nucleic acid bases and nucleotide analogues by tautomerism has been proposed to be a powerful on/off switching mechanism allowing regulation of many biological processes mediated by RNA enzymes and aptamers. Despite the suspected biological importance of tautomerism, attempts to observe minor tautomeric forms in nucleic acid or hybrid nucleic acid-ligand complexes have met with challenges due to the lack of sensitive methods. Here, a combination of spectroscopic, biochemical, and computational tools probed tautomerism in the context of an RNA aptamer-ligand complex; studies involved a model ligand, oxythiamine pyrophosphate (OxyTPP), bound to the thiamine pyrophosphate (TPP) riboswitch (an RNA aptamer) as well as its unbound nonphosphorylated form, oxythiamine (OxyT). OxyTPP, similarly to canonical heteroaromatic nucleic acid bases, has a pyrimidine ring that forms hydrogen bonding interactions with the riboswitch. Tautomerism was established using two-dimensional infrared (2D IR) spectroscopy, variable temperature FTIR and NMR spectroscopies, binding isotope effects (BIEs), and computational methods. All three possible tautomers of OxyT, including the minor enol tautomer, were directly identified, and their distributions were quantitated. In the bound form, BIE data suggested that OxyTPP existed as a 4'-keto tautomer that was likely protonated at the N1'-position. These results also provide a mechanistic framework for understanding the activation of riboswitch in response to deamination of the active form of vitamin B1 (or TPP). The combination of methods reported here revealing the fine details of tautomerism can be applied to other systems where the importance of tautomerism is suspected.

Experimental approaches for measuring pKa's in RNA and DNA

RNA and DNA carry out diverse functions in biology including catalysis, splicing, gene regulation, and storage of genetic information. Interest has grown in understanding how nucleic acids perform such sophisticated functions given their limited molecular repertoire. RNA can fold into diverse shapes that often perturb pKa values and allow it to ionize appreciably under biological conditions, thereby extending its molecular diversity. The goal of this chapter is to enable experimental measurement of pKa's in RNA and DNA. A number of experimental methods for measuring pKa values in RNA and DNA have been developed over the last 10 years, including RNA cleavage kinetics; UV-, fluorescence-, and NMR-detected pH titrations; and Raman crystallography. We begin with general considerations for choosing a pKa assay and then describe experimental conditions, advantages, and disadvantages for these assays. Potential pitfalls in measuring a pKa are provided including the presence of apparent pKa's due to a kinetic pKa or coupled acid- and alkali-promoted RNA unfolding, as well as degradation of RNA, precipitation of metal hydroxides and poor baselines. Use of multiple data fitting procedures and the study of appropriate mutants are described as ways to avoid some of these pitfalls. Application of these experimental methods to RNA and DNA will increase the number of available nucleic acid pKa values in the literature, which should deepen insight into biology and provide benchmarks for pKa calculations. Future directions for measuring pKa's in nucleic acids are discussed.

Probing riboswitch binding sites with molecular docking, focused libraries, and in-line probing assays

Molecular docking calculations combined with chemically focused libraries can bring insight in the exploration of the structure-activity relationships for a series of related compounds against an RNA target. Yet, the in silico engine must be fueled by experimental observations to drive the research into a more effective ligand-discovery path. Here we show how molecular docking predictions can be coupled with in-line probing assays to explore the available chemical and configurational space in a riboswitch binding pocket.

Thermodynamics and kinetics of adaptive binding in the malachite green RNA aptamer

Adaptive binding, the ability of molecules to fold themselves around the structure of a ligand and thereby incorporating it into their three-dimensional fold, is a key feature of most RNA aptamers. The malachite green aptamer (MGA) has been shown to bind several closely related triphenyl dyes with planar and nonplanar structures in this manner. Competitive binding studies using isothermal titration calorimetry and stopped flow kinetics have been conducted with the aim of understanding the adaptive nature of RNA-ligand interaction. The results of these studies reveal that binding of one ligand can reduce the ability of the aptamer pocket to adapt to another ligand, even if this second ligand has a significantly higher affinity to the free aptamer. A similar effect is observed in the presence of Mg(2+) ions which stabilize the binding pocket in a more ligand bound-like conformation.

The promise of riboswitches as potential antibacterial drug targets

Riboswitches represent promising novel RNA structures for developing compounds that artificially regulate gene expression and, thus, bacterial growth. The past years have seen increasing efforts to identify metabolite-analogues which act on riboswitches and which reveal antibacterial activity. Here, we summarize the current inventory of riboswitch-targeting compounds, their characteristics and antibacterial potential.

Loop-loop interaction in an adenine-sensing riboswitch: a molecular dynamics study

Riboswitches are mRNA-based molecules capable of controlling the expression of genes. They undergo conformational changes upon ligand binding, and as a result, they inhibit or promote the expression of the associated gene. The close connection between structural rearrangement and function makes a detailed knowledge of the molecular interactions an important step to understand the riboswitch mechanism and efficiency. We have performed all-atom molecular dynamics simulations of the adenine-sensing add A-riboswitch to study the breaking of the kissing loop, one key tertiary element in the aptamer structure. We investigated the aptamer domain of the add A-riboswitch in complex with its cognate ligand and in the absence of the ligand. The opening of the hairpins was simulated using umbrella sampling using the distance between two loops as the reaction coordinate. A two-step process was observed in all the simulated systems. First, a general loss of stacking and hydrogen bond interactions is seen. The last interactions that break are the two base pairs G37-C61 and G38-C60, but the break does not affect the energy profile, indicating their pivotal role in the tertiary structure formation but not in the structure stabilization. The junction area is partially organized before the kissing loop formation and residue A24 anchors together the loop helices. Moreover, when the distance between the loops is increased, one of the hairpins showed more flexibility by changing its orientation in the structure, while the other conserved its coaxial arrangement with the rest of the structure.

Single-molecule chemical denaturation of riboswitches

To date, single-molecule RNA science has been developed almost exclusively around the effect of metal ions as folding promoters and stabilizers of the RNA structure. Here, we introduce a novel strategy that combines single-molecule Förster resonance energy transfer (FRET) and chemical denaturation to observe and manipulate RNA dynamics. We demonstrate that the competing interplay between metal ions and denaturant agents provides a platform to extract information that otherwise will remain hidden with current methods. Using the adenine-sensing riboswitch aptamer as a model, we provide strong evidence for a rate-limiting folding step of the aptamer domain being modulated through ligand binding, a feature that is important for regulation of the controlled gene. In the absence of ligand, the rate-determining step is dominated by the formation of long-range key tertiary contacts between peripheral stem-loop elements. In contrast, when the adenine ligand interacts with partially folded messenger RNAs, the aptamer requires specifically bound Mg²⁺ ions, as those observed in the crystal structure, to progress further towards the native form. Moreover, despite that the ligand-free and ligand-bound states are indistinguishable by FRET, their different stability against urea-induced denaturation allowed us to discriminate them, even when they coexist within a single FRET trajectory; a feature not accessible by existing methods.

Nucleotides adjacent to the ligand-binding pocket are linked to activity tuning in the purine riboswitch

Direct sensing of intracellular metabolite concentrations by riboswitch RNAs provides an economical and rapid means to maintain metabolic homeostasis. Since many organisms employ the same class of riboswitch to control different genes or transcription units, it is likely that functional variation exists in riboswitches such that activity is tuned to meet cellular needs. Using a bioinformatic approach, we have identified a region of the purine riboswitch aptamer domain that displays conservation patterns linked to riboswitch activity. Aptamer domain compositions within this region can be divided into nine classes that display a spectrum of activities. Naturally occurring compositions in this region favor rapid association rate constants and slow dissociation rate constants for ligand binding. Using X-ray crystallography and chemical probing, we demonstrate that both the free and bound states are influenced by the composition of this region and that modest sequence alterations have a dramatic impact on activity. The introduction of non-natural compositions result in the inability to regulate gene expression in vivo, suggesting that aptamer domain activity is highly plastic and thus readily tunable to meet cellular needs.

Metabolic profiling in Caenorhabditis elegans provides an unbiased approach to investigations of dosage dependent lead toxicity

The nematode, Caenorhabditis elegans (CE), serves as a model system in which to explore the impact of particularly low-levels of lead [250, 500, 1000 and 2000 parts per million (ppm) (1.4 × 10(-6) M to 1.1 × 10(-5) M/nematode)] on specific metabolic pathways and processes. Chromatographic profiles of redox active metabolites are captured through application of high performance liquid chromatography coupled to electrochemical detection (Coularray/HPLC). Principal Component Analysis (PCA: unbiased cluster analysis) and the application of a slicing program, located significant areas of difference occurring within the 2.8-4.58 min section of the chromatograms. It is within this region of the data profiles that known components of the purine pathway reside. Two analytes of unknown structure were detected at 3.5 and 4 min respectively. Alterations in levels of the purine, tryptophan and tyrosine pathway intermediates measured in response to differing concentrations of lead acetate indicate that the effect of lead on these pathways is not linear, yet the ratio of the pathway precursors, tryptophan and tyrosine remains relatively constant. The application of the above combined analytical approaches enhances the value of data generated. Exposure of CE to very low levels of lead produced significant alterations in profiles of electrochemically active compounds. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11306-012-0438-0) contains supplementary material, which is available to authorized users.

Sensing organic molecules by charge transfer through aptamer-target complexes: theory and simulation

Aptamers, i.e., short sequences of RNA and single-stranded DNA, are capable of specificilly binding objects ranging from small molecules over proteins to entire cells. Here, we focus on the structure, stability, dynamics, and electronic properties of oligonucleotides that interact with aromatic or heterocyclic targets. Large-scale molecular dynamics simulations indicate that aromatic rings such as dyes, metabolites, or alkaloides form stable adducts with their oligonucleotide host molecules at least on the simulation time scale. From molecular dynamics snapshots, the energy parameters relevant to Marcus' theory of charge transfer are computed using a modified Su-Schrieffer-Heeger Hamiltonian, permitting an estimate of the charge transfer rates. In many cases, aptamer binding seriously influences the charge transfer kinetics and the charge carrier mobility within the complex, with conductivities up to the nanoampere range for a single complex. We discuss the conductivity properties with reference to potential applications as biosensors.

Structure and mechanism of purine-binding riboswitches

A riboswitch is a non-protein coding sequence capable of directly binding a small molecule effector without the assistance of accessory proteins to regulate expression of the mRNA in which it is embedded. Currently, over 20 different classes of riboswitches have been validated in bacteria with the promise of many more to come, making them an important means of regulating the genome in the bacterial kingdom. Strikingly, half of the known riboswitches recognize effector compounds that contain a purine or related moiety. In the last decade, significant progress has been made to determine how riboswitches specifically recognize these compounds against the background of many other similar cellular metabolites and transduce this signal into a regulatory response. Of the known riboswitches, the purine family containing guanine, adenine and 2'-deoxyguanosine-binding classes are the most extensively studied, serving as a simple and useful paradigm for understanding how these regulatory RNAs function. This review provides a comprehensive summary of the current state of knowledge regarding the structure and mechanism of these riboswitches, as well as insights into how they might be exploited as therapeutic targets and novel biosensors.

Insights into the regulatory landscape of the lysine riboswitch

A prevalent means of regulating gene expression in bacteria is by riboswitches found within mRNA leader sequences. Like protein repressors, these RNA elements must bind an effector molecule with high specificity against a background of other cellular metabolites of similar chemical structure to elicit the appropriate regulatory response. Current crystal structures of the lysine riboswitch do not provide a complete understanding of selectivity as recognition is substantially mediated through main-chain atoms of the amino acid. Using a directed set of lysine analogs and other amino acids, we have determined the relative contributions of the polar functional groups to binding affinity and the regulatory response. Our results reveal that the lysine riboswitch has >1000-fold specificity for lysine over other amino acids. The aptamer is highly sensitive to the precise placement of the ε-amino group and relatively tolerant of alterations to the main-chain functional groups in order to achieve this specificity. At low nucleotide triphosphate (NTP) concentrations, we observe good agreement between the half-maximal regulatory activity (T(50)) and the affinity of the receptor for lysine (K(d)), as well as many of its analogs. However, above 400 μM [NTP], the concentration of lysine required to elicit transcription termination rises, moving into the riboswitch into a kinetic control regime. These data demonstrate that, under physiologically relevant conditions, riboswitches can integrate both effector and NTP concentrations to generate a regulatory response appropriate for global metabolic state of the cell.

Mechanisms for differentiation between cognate and near-cognate ligands by purine riboswitches

Riboswitches are elements in the 5'-untranslated region of mRNAs that regulate gene expression by directly interacting with metabolites related to their own gene products. A remarkable feature of this gene regulation mechanism is the high specificity of riboswitches for their cognate ligands. In this study, we used a combination of static and time-resolved NMR-spectroscopic methods to investigate the mechanisms for ligand specificity in purine riboswitches. We investigate the xpt-aptamer domain from a guanine-responsive riboswitch and the mfl-aptamer domain from a 2'-deoxyguanosine-responsive riboswitch. The xpt-aptamer binds the purine nucleobases guanine/hypoxanthine with high affinity, but, unexpectedly, also the nucleoside 2'-deoxyguanosine. On the other hand, the mfl-aptamer is highly specific for its cognate ligand 2'-deoxyguanosine, and does not bind purine ligands. We addressed the question of aptamer`s ligand specificity by real-time NMR spectroscopy. Our studies of ligand binding and subsequently induced aptamer folding revealed that the xpt-aptamer discriminates against non-cognate ligands by enhanced life-times of the cognate complex compared with non-cognate complexes, whereas the mfl-aptamer rejects non-cognate ligands at the level of ligand association, employing a kinetic proofreading mechanism.

Differential analogue binding by two classes of c-di-GMP riboswitches

The ability of bacteria to adapt to a changing environment is essential for their survival. One mechanism bacteria have evolved to sense environmental cues and translate these signals into phenotypic changes uses the second messenger signaling molecule, cyclic diguanosine monophosphate (c-di-GMP). In addition to several classes of protein receptors, two classes of c-di-GMP-binding riboswitches (class I and class II) have been identified as downstream targets of the second messenger in this signaling pathway. The crystal structures of both riboswitch classes bound to c-di-GMP were previously reported. Here, we further investigate the mechanisms that RNA has evolved for recognition and binding of this second messenger. Using a series of c-di-GMP analogues, we probed the interactions made in the RNA-ligand complex for both classes of riboswitches to identify the most critical elements of c-di-GMP for binding. We found that the structural features of c-di-GMP required for binding differ between these two effectors and that the class II riboswitch is much less discriminatory in ligand binding than the class I riboswitch. These data suggest an explanation for the predicted preferential use of the class I motif over the class II motif in the c-di-GMP signaling pathway.

Influence of ground-state structure and Mg2+ binding on folding kinetics of the guanine-sensing riboswitch aptamer domain

Riboswitch RNAs fold into complex tertiary structures upon binding to their cognate ligand. Ligand recognition is accomplished by key residues in the binding pocket. In addition, it often crucially depends on the stability of peripheral structural elements. The ligand-bound complex of the guanine-sensing riboswitch from Bacillus subtilis, for example, is stabilized by extensive interactions between apical loop regions of the aptamer domain. Previously, we have shown that destabilization of this tertiary loop–loop interaction abrogates ligand binding of the G37A/C61U-mutant aptamer domain (Gsw^loop) in the absence of Mg²⁺. However, if Mg²⁺ is available, ligand-binding capability is restored by a population shift of the ground-state RNA ensemble toward RNA conformations with pre-formed loop–loop interactions. Here, we characterize the striking influence of long-range tertiary structure on RNA folding kinetics and on ligand-bound complex structure, both by X-ray crystallography and time-resolved NMR. The X-ray structure of the ligand-bound complex reveals that the global architecture is almost identical to the wild-type aptamer domain. The population of ligand-binding competent conformations in the ground-state ensemble of Gsw^loop is tunable through variation of the Mg²⁺ concentration. We quantitatively describe the influence of distinct Mg²⁺ concentrations on ligand-induced folding trajectories both by equilibrium and time-resolved NMR spectroscopy at single-residue resolution.

Structural principles of nucleoside selectivity in a 2'-deoxyguanosine riboswitch

Purine riboswitches play an essential role in genetic regulation of bacterial metabolism. This family includes the 2′-deoxyguanosine (dG) riboswitch, involved in feedback control of deoxyguanosine biosynthesis. To understand the principles that define dG selectivity, we determined crystal structures of natural Mesoplasma florum riboswitch bound to cognate dG, as well as non-cognate guanosine, deoxyguanosine monophosphate and guanosine monophosphate. Comparison with related purine riboswitch structures reveals that the dG riboswitch achieves its specificity by modifying key interactions involving the nucleobase and through rearrangement of the ligand-binding pocket, so as to accommodate the additional sugar moiety. In addition, we observe novel conformational changes beyond the junctional binding pocket, extending as far as peripheral loop-loop interactions. It appears that re-engineering riboswitch scaffolds will require consideration of selectivity features dispersed throughout the riboswitch tertiary fold, and that structure-guided drug design efforts targeted to junctional RNA scaffolds need to be addressed within such an expanded framework.

Novel ligands for a purine riboswitch discovered by RNA-ligand docking

The increasing number of RNA crystal structures enables a structure-based approach to the discovery of new RNA-binding ligands. To develop the poorly explored area of RNA-ligand docking, we have conducted a virtual screening exercise for a purine riboswitch to probe the strengths and weaknesses of RNA-ligand docking. Using a standard protein-ligand docking program with only minor modifications, four new ligands with binding affinities in the micromolar range were identified, including two compounds based on molecular scaffolds not resembling known ligands. RNA-ligand docking performed comparably to protein-ligand docking indicating that this approach is a promising option to explore the wealth of RNA structures for structure-based ligand design.

Structure and dynamics of the deoxyguanosine-sensing riboswitch studied by NMR-spectroscopy

The mfl-riboswitch regulates expression of ribonucleotide reductase subunit in Mesoplasma florum by binding to 2′-deoxyguanosine and thereby promoting transcription termination. We characterized the structure of the ligand-bound aptamer domain by NMR spectroscopy and compared the mfl-aptamer to the aptamer domain of the closely related purine-sensing riboswitches. We show that the mfl-aptamer accommodates the extra 2′-deoxyribose unit of the ligand by forming a more relaxed binding pocket than these found in the purine-sensing riboswitches. Tertiary structures of the xpt-aptamer bound to guanine and of the mfl-aptamer bound to 2′-deoxyguanosine exhibit very similar features, although the sequence of the mfl-aptamer contains several alterations compared to the purine-aptamer consensus sequence. These alterations include the truncation of a hairpin loop which is crucial for complex formation in all purine-sensing riboswitches characterized to date. We further defined structural features and ligand binding requirements of the free mfl-aptamer and found that the presence of Mg²⁺ is not essential for complex formation, but facilitates ligand binding by promoting pre-organization of key structural motifs in the free aptamer.

Heterogeneity and dynamics of the ligand recognition mode in purine-sensing riboswitches

High-resolution crystal structures and biophysical analyses of purine-sensing riboswitches have revealed that a network of hydrogen bonding interactions appear to be largey responsible for discrimination of cognate ligands against structurally related compounds. Here we report that by using femtosecond time-resolved fluorescence spectroscopy to capture the ultrafast decay dynamics of the 2-aminopurine base as the ligand, we have detected the presence of multiple conformations of the ligand within the binding pockets of one guanine-sensing and two adenine-sensing riboswitches. All three riboswitches have similar conformational distributions of the ligand-bound state. The known crystal structures represent the global minimum that accounts for 50-60% of the population, where there is no significant stacking interaction between the ligand and bases of the binding pocket, but the hydrogen-bonding cage collectively provides an electronic environment that promotes an ultrafast ( approximately 1 ps) charge transfer pathway. The ligand also samples multiple conformations in which it significantly stacks with either the adenine or the uracil bases of the A21-U75 and A52-U22 base pairs that form the ceiling and floor of the binding pocket, respectively, but favors the larger adenine bases. These alternative conformations with well-defined base stacking interactions are approximately 1-1.5 kcal/mol higher in DeltaG degrees than the global minimum and have distinct charge transfer dynamics within the picosecond to nanosecond time regime. Inside the pocket, the purine ligand undergoes dynamic motion on the low nanosecond time scale, sampling the multiple conformations based on time-resolved anisotropy decay dynamics. These results allowed a description of the energy landscape of the bound ligand with intricate details and demonstrated the elastic nature of the ligand recognition mode by the purine-sensing riboswitches, where there is a dynamic balance between hydrogen bonding and base stacking interactions, yielding the high affinity and specificity by the aptamer domain.

Base pairs and pseudo pairs observed in RNA-ligand complexes

Previously, a geometric nomenclature was proposed in which RNA base pairs were classified by their interaction edges (Watson-Crick, Hoogsteen or sugar-edge) and the glycosidic bond orientations relative to the hydrogen bonds formed (cis or trans). Here, base pairs and pseudo pairs observed in RNA-ligand complexes are classified in a similar manner. Twenty-one basic geometric families are geometrically possible (18 for base pairs formed between a nucleic acid base and a ligand containing heterocycle and 3 families for pseudo pairs). Of those, 16 of them have been observed in X-ray and/or NMR structures.

Multivector fluorescence analysis of the xpt guanine riboswitch aptamer domain and the conformational role of guanine

Purine riboswitches are RNA regulatory elements that control purine metabolism in response to intracellular concentrations of the purine ligands. Conformational changes of the guanine riboswitch aptamer domain induced by guanine binding lead to transcriptional regulation of genes involved in guanine biosynthesis. The guanine riboswitch aptamer domain has three RNA helices designated P1, P2, and P3. An overall model for the Mg²⁺- and guanine-dependent relative orientations and dynamics of P1, P2, and P3 has not been reported, and the conformational role of guanine under physiologically relevant conditions has not been fully elucidated. In this study, an ensemble and single-molecule fluorescence resonance energy transfer (FRET) study was performed on three orthogonally labeled variants of the xpt guanine riboswitch aptamer domain. The combined FRET data support a model in which the unfolded state of the aptamer domain has a highly dynamic P2 helix that switches rapidly between two orientations relative to nondynamic P1 and P3. At ≪1 mM Mg²⁺ (in the presence of a saturating level of guanine) or ≥1 mM Mg²⁺ (in the absence of guanine), the riboswitch starts to adopt a folded conformation in which loop−loop interactions lock P2 and P3 into place. At >5 mM Mg²⁺, further compaction occurs in which P1 more closely approaches P3. Our data help to explain the biological role of guanine as stabilizing the globally folded aptamer domain conformation at physiologically relevant Mg²⁺ concentrations (≤1 mM), whereas in the absence of guanine, much higher Mg²⁺ concentrations are required to induce this folding event.

Predicting loop-helix tertiary structural contacts in RNA pseudoknots

Tertiary interactions between loops and helical stems play critical roles in the biological function of many RNA pseudoknots. However, quantitative predictions for RNA tertiary interactions remain elusive. Here we report a statistical mechanical model for the prediction of noncanonical loop-stem base-pairing interactions in RNA pseudoknots. Central to the model is the evaluation of the conformational entropy for the pseudoknotted folds with defined loop-stem tertiary structural contacts. We develop an RNA virtual bond-based conformational model (Vfold model), which permits a rigorous computation of the conformational entropy for a given fold that contains loop-stem tertiary contacts. With the entropy parameters predicted from the Vfold model and the energy parameters for the tertiary contacts as inserted parameters, we can then predict the RNA folding thermodynamics, from which we can extract the tertiary contact thermodynamic parameters from theory-experimental comparisons. These comparisons reveal a contact enthalpy (DeltaH) of -14 kcal/mol and a contact entropy (DeltaS) of -38 cal/mol/K for a protonated C(+)*(G-C) base triple at pH 7.0, and (DeltaH = -7 kcal/mol, DeltaS = -19 cal/mol/K) for an unprotonated base triple. Tests of the model for a series of pseudoknots show good theory-experiment agreement. Based on the extracted energy parameters for the tertiary structural contacts, the model enables predictions for the structure, stability, and folding pathways for RNA pseudoknots with known or postulated loop-stem tertiary contacts from the nucleotide sequence alone.

Reengineering orthogonally selective riboswitches

The ability to independently control the expression of multiple genes by addition of distinct small-molecule modulators has many applications from synthetic biology, functional genomics, pharmaceutical target validation, through to gene therapy. Riboswitches are relatively simple, small-molecule-dependent, protein-free, mRNA genetic switches that are attractive targets for reengineering in this context. Using a combination of chemical genetics and genetic selection, we have developed riboswitches that are selective for synthetic "nonnatural" small molecules and no longer respond to the natural intracellular ligands. The orthogonal selectivity of the riboswitches is also demonstrated in vitro using isothermal titration calorimetry and x-ray crystallography. The riboswitches allow highly responsive, dose-dependent, orthogonally selective, and dynamic control of gene expression in vivo. It is possible that this approach may be further developed to reengineer other natural riboswitches for application as small-molecule responsive genetic switches in both prokaryotes and eukaryotes.

Directional Weak Intermolecular Interactions: σ-Hole Bonding

The prototypical directional weak interactions, hydrogen bonding and σ-hole bonding (including the special case of halogen bonding) are reviewed in a united picture that depends on the anisotropic nature of the molecular electrostatic potential around the donor atom. Qualitative descriptions of the effects that lead to these anisotropic distributions are given and examples of the importance of σ-hole bonding in crystal engineering and biological systems are discussed.

Role of the adenine ligand on the stabilization of the secondary and tertiary interactions in the adenine riboswitch

Riboswitches are RNA-based genetic control elements that function via a conformational transition mechanism when a specific target molecule binds to its binding pocket. To facilitate an atomic detail interpretation of experimental investigations on the role of the adenine ligand on the conformational properties and kinetics of folding of the add adenine riboswitch, we performed molecular dynamics simulations in both the presence and the absence of the ligand. In the absence of ligand, structural deviations were observed in the J23 junction and the P1 stem. Destabilization of the P1 stem in the absence of ligand involves the loss of direct stabilizing interactions with the ligand, with additional contributions from the J23 junction region. The J23 junction of the riboswitch is found to be more flexible, and the tertiary contacts among the junction regions are altered in the absence of the adenine ligand; results suggest that the adenine ligand associates and dissociates from the riboswitch in the vicinity of J23. Good agreement was obtained with the experimental data with the results indicating dynamic behavior of the adenine ligand on the nanosecond time scale to be associated with the dynamic behavior of hydrogen bonding with the riboswitch. Results also predict that direct interactions of the adenine ligand with U74 of the riboswitch are not essential for stable binding although it is crucial for its recognition. The possibility of methodological artifacts and force-field inaccuracies impacting the present observations was checked by additional molecular dynamics simulations in the presence of 2,6-diaminopurine and in the crystal environment.