The ligand-free state of the TPP riboswitch: a partially folded RNA structure

Supercoiling-dependent DNA binding: quantitative modeling and applications to bulk and single-molecule experiments

DNA stores our genetic information and is ubiquitous in applications, where it interacts with binding partners ranging from small molecules to large macromolecular complexes. Binding is modulated by mechanical strains in the molecule and can change local DNA structure. Frequently, DNA occurs in closed topological forms where topology and supercoiling add a global constraint to the interplay of binding-induced deformations and strain-modulated binding. Here, we present a quantitative model with a straight-forward numerical implementation of how the global constraints introduced by DNA topology modulate binding. We focus on fluorescent intercalators, which unwind DNA and enable direct quantification via fluorescence detection. Our model correctly describes bulk experiments using plasmids with different starting topologies, different intercalators, and over a broad range of intercalator and DNA concentrations. We demonstrate and quantitatively model supercoiling-dependent binding in a single-molecule assay, where we directly observe the different intercalator densities going from supercoiled to nicked DNA. The single-molecule assay provides direct access to binding kinetics and DNA supercoil dynamics. Our model has broad implications for the detection and quantification of DNA, including the use of psoralen for UV-induced DNA crosslinking to quantify torsional tension in vivo, and for the modulation of DNA binding in cellular contexts.

Crystal structure of Escherichia coli thiamine pyrophosphate-sensing riboswitch in the apo state

The thiamine pyrophosphate (TPP)-sensing riboswitch is one of the earliest discovered and most widespread riboswitches. Numerous structural studies have been reported for this riboswitch bound with various ligands. However, the ligand-free (apo) structure remains unknown. Here, we report a 3.1 Å resolution crystal structure of Escherichia coli TPP riboswitch in the apo state, which exhibits an extended, Y-shaped conformation further supported by small-angle X-ray scattering data and driven molecular dynamics simulations. The loss of ligand interactions results in helical uncoiling of P5 and disruption of the key tertiary interaction between the sensory domains. Opening of the aptamer propagates to the gene-regulatory P1 helix and generates the key conformational flexibility needed for the switching behavior. Much of the ligand-binding site at the three-way junction is unaltered, thereby maintaining a partially preformed pocket. Together, these results paint a dynamic picture of the ligand-induced conformational changes in TPP riboswitches that confer conditional gene regulation.

Exploring the structure, function of thiamine pyrophosphate riboswitch, and designing small molecules for antibacterial activity

During the last decade, riboswitches emerged as new small-molecule sensing RNA in bacteria. Thiamine pyrophosphate (TPP) riboswitch is widely distributed and occurs in plants, bacteria, fungi, and archaea. Extensive biochemical, structural, and genetic studies have been carried out to elucidate the recognition mechanism of TPP riboswitches. However, a comprehensive report summarizing all information on recognition principles and newly designed ligands for TPP riboswitch is scarce in the literature. This review gives a comprehensive understanding of the TPP riboswitch's structure, mechanism, and methods applied to design ligands for the TPP riboswitch. The ligand-bound TPP riboswitch was studied with various experimental and theoretical techniques to elucidate the conformational dynamics. The mutation studies shed light on the significance of pyrimidine sensing helix for the binding of ligands. Further, the structure-activity relationship study and fragment-based approach lead to the development of ligands with Kd values at the sub-micromolar level. However, there is a need to design more potent inhibitors for TPP riboswitch for therapeutic applications. The recent advancements in ligand design highlight the TPP riboswitch as a promising target for developing new antibiotics. This article is categorized under: Regulatory RNAs/RNAi/Riboswitches > Riboswitches Regulatory RNAs/RNAi/Riboswitches > Regulatory RNAs.

Tools shaping drug discovery and development

Spectroscopic, scattering, and imaging methods play an important role in advancing the study of pharmaceutical and biopharmaceutical therapies. The tools more familiar to scientists within industry and beyond, such as nuclear magnetic resonance and fluorescence spectroscopy, serve two functions: as simple high-throughput techniques for identification and purity analysis, and as potential tools for measuring dynamics and structures of complex biological systems, from proteins and nucleic acids to membranes and nanoparticle delivery systems. With the expansion of commercial small-angle x-ray scattering instruments into the laboratory setting and the accessibility of industrial researchers to small-angle neutron scattering facilities, scattering methods are now used more frequently in the industrial research setting, and probe-less time-resolved small-angle scattering experiments are now able to be conducted to truly probe the mechanism of reactions and the location of individual components in complex model or biological systems. The availability of atomic force microscopes in the past several decades enables measurements that are, in some ways, complementary to the spectroscopic techniques, and wholly orthogonal in others, such as those related to nanomechanics. As therapies have advanced from small molecules to protein biologics and now messenger RNA vaccines, the depth of biophysical knowledge must continue to serve in drug discovery and development to ensure quality of the drug, and the characterization toolbox must be opened up to adapt traditional spectroscopic methods and adopt new techniques for unraveling the complexities of the new modalities. The overview of the biophysical methods in this review is meant to showcase the uses of multiple techniques for different modalities and present recent applications for tackling particularly challenging situations in drug development that can be solved with the aid of fluorescence spectroscopy, nuclear magnetic resonance spectroscopy, atomic force microscopy, and small-angle scattering.

Molecular dynamics simulation on the Thermosinus carboxydivorans pfl ZTP riboswitch by ligand binding

Riboswitches are RNA molecules that can regulate gene expression which is affected by ligand-binding during cotranscriptional folding process. However, the role of ligand during the folding is still unclear. In this study, the pfl domain of Thermosinus carboxydivorans ZTP riboswitch was discussed. The ligand is molecule ZMP. We mainly analyzed the change of ZMP-free and ZMP-bound aptamer domain by the dynamics simulation method. Structural features by calculating their RMSD, RMSF, etc. are analyzed. The results demonstrate that the binding domain require the presence of ZMP to maintain a stable fold. It also suggested that ZMP specificly binding to ZTP can generate more hydrogen bonds in the binding domain. Through the calculation of binding free energy decomposition of each nucleotide, molecule ZMP was found to promote the recognition and binding process of ligands by controlling some special nucleotides in the process of ligand binding. At last, the dynamical correlation and components of conformational motions were both applied to explore the effect of molecule ZMP to ZTP riboswitch. In general, ZMP can effectively affect the motions of the pfl riboswitch and facilitate the folding process of the ZTP riboswitch.These results may provide some new ideas for structural changes in riboswitches and their cotranscriptional folding process.

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

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

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.

Ensemble Switching Unveils a Kinetic Rheostat Mechanism of the Eukaryotic Thiamine Pyrophosphate Riboswitch

Thiamine pyrophosphate (TPP) riboswitches regulate thiamine metabolism by inhibiting the translation of enzymes essential to thiamine synthesis pathways upon binding to thiamine pyrophosphate in cells across all domains of life. Recent work on the Arabidopsis thaliana TPP riboswitch suggests a multi-step TPP binding process involving multiple riboswitch configurational ensembles and that Mg2+ dependence underlies the mechanism of TPP recognition and subsequent transition to the expression-inhibiting state of the aptamer domain followed by changes in the expression platform. However, details of the relationship between TPP riboswitch conformational changes and interactions with TPP and Mg2+ ¬¬in the aptamer domain constituting this mechanism are unknown. Therefore, we integrated single-molecule multiparameter fluorescence and force spectroscopy with atomistic molecular dynamics simulations and found that conformational transitions within the aptamer domain's sensor helices associated with TPP and Mg2+ ligand binding occurred between at least five different ensembles on timescales ranging from µs to ms. These dynamics are orders of magnitude faster than the 10 second-timescale folding kinetics associated with expression-state switching in the switch sequence. Together, our results show that a TPP and Mg2+ dependent mechanism determines dynamic configurational state ensemble switching of the aptamer domain's sensor helices that regulates the stability of the switch helix, which ultimately may lead to the expression-inhibiting state of the riboswitch. Additionally, we propose that two pathways exist for ligand recognition and that this mechanism underlies a kinetic rheostat-like behavior of the Arabidopsis thaliana TPP riboswitch.

Molecular level insights into the inhibition of gene expression by thiamine pyrophosphate (TPP) analogs for TPP riboswitch: A well-tempered metadynamics simulations study

Riboswitches are metabolite sensing aptamer domains present in non-coding regions in RNA and act as gene-regulating elements. Thiamine pyrophosphate (TPP) riboswitch is evolved as a new target for developing antibiotics against many pathogenic bacteria. The earlier reports suggest that the modification of the pyrophosphate group in the ligand molecule can enhance gene expression. In this work, we have examined the binding affinity and efficacy of TPP and two recently reported ligands, CH₂-TPP, and CF₂-TPP, using Well-tempered metadynamics (WT-MtD) simulations. The experimental in vitro assays show that both TPP and CH₂-TPP repress the gene expression to the same extent. The calculated binding energies correlate well with the experimental study and show the same trend of binding affinity of ligands for the TPP riboswitch. The root mean square fluctuation profiles suggest that the CH₂-TPP and TPP trigger higher fluctuations in P1 and L3 region, and such fluctuations in the P1 region is involved in the gene regulation process. The metal ion mediated contact of TPP ligand with pyrophosphate binding helix is found to be critical in the gene regulation process. The simulation results corroborate the experimental observations that the role of conformational changes occurring in different riboswitch regions upon ligand binding is essential to repress the gene expression process. This work sheds light on the subtle change in the ligand structure that can induce a more considerable impact on binding affinity and efficacy of ligands with riboswitch.

Finely tuned conformational dynamics regulate the protective function of the lncRNA MALAT1 triple helix

Nucleic acid triplexes may regulate many important biological processes. Persistent accumulation of the oncogenic 7-kb long noncoding RNA MALAT1 is dependent on an unusually long intramolecular triple helix. This triplex structure is positioned within a conserved ENE (element for nuclear expression) motif at the lncRNA 3′ terminus and protects the entire transcript from degradation in a polyA-independent manner. A requisite 3′ maturation step leads to triplex formation though the precise mechanism of triplex folding remains unclear. Furthermore, the contributions of several peripheral structural elements to triplex formation and protective function have not been determined. We evaluated the stability, conformational fluctuations, and function of this MALAT1 ENE triple helix (M1^TH) protective element using in vitro mutational analyses coupled with biochemical and biophysical characterizations. Using fluorescence and UV melts, FRET, and an exonucleolytic decay assay we define a concerted mechanism for triplex formation and uncover a metastable, dynamic triplex population under near-physiological conditions. Structural elements surrounding the triplex regulate the dynamic M1^TH conformational variability, but increased triplex dynamics lead to M1^TH degradation. Taken together, we suggest that finely tuned dynamics may be a general mechanism regulating triplex-mediated functions.

Construction of Boolean logic gates based on dual-vector circuits of multiple gene regulatory elements

Gene circuits are constructed to run complex logical operations for the precise regulation of biological metabolic processes. At present, the implementation of most genetic circuits is based on the regulatory mechanism of various circuit components, but we hope to realize complex logic gates through biological metabolic pathways of organisms. In this study, we matched the regulatory elements of different functional mechanisms to build a Boolean logic gate model by means of a dual-vector circuit. In Escherichia coli, we made 12 circuit logic gate modules and validated the functions of four of the logic gates, including "AND", "NAND", "OR" and "NOR" by the expression and analysis of a reporter gene. The inputs were converted into outputs by an intermediate product of the host metabolism. The results indicated that these logic gate circuits had the expected efficacy and regulatory characteristics. Our study provides new ideas for designing genetic circuits and precisely controlling metabolic pathways.

TPP riboswitch aptamer: Role of Mg2+ ions, ligand unbinding, and allostery

Riboswitches are non-coding RNAs that regulate gene expression in response to the binding of metabolites. Their abundance in bacteria makes them ideal drug targets. The prokaryotic thiamine pyrophosphate (TPP) riboswitch regulates gene expression in a wide range of bacteria by undergoing conformational changes in response to the binding of TPP. Although an experimental structure for the aptamer domain of the riboswitch is now available, details of the conformational changes that occur during the binding of the ligand, and the factors that govern these conformational changes, are still not clear. This study employs microsecond-scale molecular dynamics simulations to provide insights into the functioning of the riboswitch aptamer in atomistic detail. A mechanism for the transmission of conformational changes from the ligand-binding site to the P1 switch helix is proposed. Mg²⁺ ions in the binding site play a critical role in anchoring the ligand to the riboswitch. Finally, modeling the egress of TPP from the binding site reveals a two-step mechanism for TPP unbinding. Findings from this study can motivate the design of future studies aimed at modulating the activity of this drug target.

Gold nanocrystal labels provide a sequence-to-3D structure map in SAXS reconstructions

Small-angle x-ray scattering (SAXS) is a powerful technique to probe the structure of biological macromolecules and their complexes under virtually arbitrary solution conditions, without the need for crystallization. While it is possible to reconstruct molecular shapes from SAXS data ab initio, the resulting electron density maps have a resolution of ~1 nm and are often insufficient to reliably assign secondary structure elements or domains. We show that SAXS data of gold-labeled samples significantly enhance the information content of SAXS measurements, allowing the unambiguous assignment of macromolecular sequence motifs to specific locations within a SAXS structure. We first demonstrate our approach for site-specifically internally and end-labeled DNA and an RNA motif. In addition, we present a protocol for highly uniform and site-specific labeling of proteins with small (~1.4 nm diameter) gold particles and apply our method to the signaling protein calmodulin. In all cases, the position of the small gold probes can be reliably identified in low-resolution electron density maps. Enhancing low-resolution measurements by site-selective gold labeling provides an attractive approach to aid modeling of a large range of macromolecular systems.

Single-molecule FRET studies on the cotranscriptional folding of a thiamine pyrophosphate riboswitch

Because RNAs fold as they are being synthesized, their transcription rate can affect their folding. Here, we report the results of single-molecule fluorescence studies that characterize the ligand-dependent cotranscriptional folding of the Escherichia coli thiM riboswitch that regulates translation. We found that the riboswitch aptamer folds into the "off" conformation independent of its ligand, but switches to the "on" conformation during transcriptional pausing near the translational start codon. Ligand binding maintains the riboswitch in the off conformation during transcriptional pauses. We expect our assay will permit the controlled study of the two main physical mechanisms that regulate cotranscriptional folding: transcriptional pausing and transcriptional speed.

Three-dimensional structure of the 3'X-tail of hepatitis C virus RNA in monomeric and dimeric states

The 3'X domain is a 98-nt region located at the 3' end of hepatitis C virus genomic RNA that plays essential functions in the viral life cycle. It contains an absolutely conserved, 16-base palindromic sequence that promotes viral RNA dimerization, overlapped with a 7-nt tract implicated in a distal contact with a nearby functional sequence. Using small angle X-ray scattering measurements combined with model building guided by NMR spectroscopy, we have studied the stoichiometry, structure, and flexibility of domain 3'X and two smaller subdomain sequences as a function of ionic strength, and obtained a three-dimensional view of the full-length domain in its monomeric and dimeric states. In the monomeric form, the 3'X domain adopted an elongated conformation containing two SL1' and SL2' double-helical stems stabilized by coaxial stacking. This structure was significantly less flexible than that of isolated subdomain SL2' monomers. At higher ionic strength, the 3'X scattering envelope nearly doubled its size, reflecting the formation of extended homodimers containing an antiparallel SL2' duplex flanked by coaxially stacked SL1' helices. Formation of these dimers could initialize and/or regulate the packaging of viral RNA genomes into virions.

Cotranscriptional folding of a riboswitch at nucleotide resolution

RNAs can begin to fold immediately as they emerge from RNA polymerase. During cotranscriptional folding, interactions between nascent RNAs and ligands are able to direct the formation of alternative RNA structures, a feature exploited by noncoding RNAs called riboswitches to make gene-regulatory decisions. Despite their importance, cotranscriptional folding pathways have yet to be uncovered with sufficient resolution to reveal how cotranscriptional folding governs RNA structure and function. To access cotranscriptional folding at nucleotide resolution, we extended selective 2'-hydroxyl acylation analyzed by primer-extension sequencing (SHAPE-seq) to measure structural information of nascent RNAs during transcription. Using cotranscriptional SHAPE-seq, we determined how the cotranscriptional folding pathway of the Bacillus cereus crcB fluoride riboswitch undergoes a ligand-dependent bifurcation that delays or promotes terminator formation via a series of coordinated structural transitions. Our results directly link cotranscriptional RNA folding to a genetic decision and establish a framework for cotranscriptional analysis of RNA structure at nucleotide resolution.

Conformational Changes and Flexibility of DNA Devices Observed by Small-Angle X-ray Scattering

Self-assembled DNA origami nanostructures enable the creation of precisely defined shapes at the molecular scale. Dynamic DNA devices that are capable of switching between defined conformations could afford completely novel functionalities for diagnostic, therapeutic, or engineering applications. Developing such objects benefits strongly from experimental feedback about conformational changes and 3D structures, ideally in solution, free of potential biases from surface attachment or labeling. Here, we demonstrate that small-angle X-ray scattering (SAXS) can quantitatively resolve the conformational changes of a DNA origami two-state switch device as a function of the ionic strength of the solution. In addition, we show how SAXS data allow for refinement of the predicted idealized three-dimensional structure of the DNA object using a normal mode approach based on an elastic network model. The results reveal deviations from the idealized design geometries that are otherwise difficult to resolve. Our results establish SAXS as a powerful tool to investigate conformational changes and solution structures of DNA origami and we anticipate our methodology to be broadly applicable to increasingly complex DNA and RNA devices.

Tetracycline determines the conformation of its aptamer at physiological magnesium concentrations

Synthetic riboswitches are versatile tools for the study and manipulation of biological systems. Yet, the underlying mechanisms governing its structural properties and regulation under physiological conditions are poorly studied. We performed spectroscopic and calorimetric experiments to explore the folding kinetics and thermodynamics of the tetracycline-binding aptamer, which can be employed as synthetic riboswitch, in the range of physiological magnesium concentrations. The dissociation constant of the ligand-aptamer complex was found to strongly depend on the magnesium concentration. At physiological magnesium concentrations, tetracycline induces a significant conformational shift from a compact, but heterogeneous intermediate state toward the completely formed set of tertiary interactions defining the regulation-competent structure. Thus, the switching functionality of the tetracycline-binding aptamer appears to include both a conformational rearrangement toward the regulation-competent structure and its thermodynamic stabilization.

Progress and Current Challenges in Modeling Large RNAs

Recent breakthroughs in next-generation sequencing technologies have led to the discovery of several classes of non-coding RNAs (ncRNAs). It is now apparent that RNA molecules are not only just carriers of genetic information but also key players in many cellular processes. While there has been a rapid increase in the number of ncRNA sequences deposited in various databases over the past decade, the biological functions of these ncRNAs are largely not well understood. Similar to proteins, RNA molecules carry out a function by forming specific three-dimensional structures. Understanding the function of a particular RNA therefore requires a detailed knowledge of its structure. However, determining experimental structures of RNA is extremely challenging. In fact, RNA-only structures represent just 1% of the total structures deposited in the PDB. Thus, computational methods that predict three-dimensional RNA structures are in high demand. Computational models can provide valuable insights into structure-function relationships in ncRNAs and can aid in the development of functional hypotheses and experimental designs. In recent years, a set of diverse RNA structure prediction tools have become available, which differ in computational time, input data and accuracy. This review discusses the recent progress and challenges in RNA structure prediction methods.

Consistent global structures of complex RNA states through multidimensional chemical mapping

Accelerating discoveries of non-coding RNA (ncRNA) in myriad biological processes pose major challenges to structural and functional analysis. Despite progress in secondary structure modeling, high-throughput methods have generally failed to determine ncRNA tertiary structures, even at the 1-nm resolution that enables visualization of how helices and functional motifs are positioned in three dimensions. We report that integrating a new method called MOHCA-seq (Multiplexed •OHCleavage Analysis with paired-end sequencing) with mutate-and-map secondary structure inference guides Rosetta 3D modeling to consistent 1-nm accuracy for intricately folded ncRNAs with lengths up to 188 nucleotides, including a blind RNA-puzzle challenge, the lariat-capping ribozyme. This multidimensional chemical mapping (MCM) pipeline resolves unexpected tertiary proximities for cyclic-di-GMP, glycine, and adenosylcobalamin riboswitch aptamers without their ligands and a loose structure for the recently discovered human HoxA9D internal ribosome entry site regulon. MCM offers a sequencing-based route to uncovering ncRNA 3D structure, applicable to functionally important but potentially heterogeneous states.

DOI:http://dx.doi.org/10.7554/eLife.07600.001

Our genetic material, in the form of molecules of DNA, provides instructions for many different processes in our cells. To issue these instructions, particular sections of DNA are copied to make a type of molecule called ribonucleic acid (RNA). Some of these RNA molecules contain instructions to make proteins, but others—known as non-coding RNAs—regulate the activity of genes in cells.

The genetic information within RNA is encoded by the sequence of four different chemical parts called ‘nucleotides’. RNA can exist as a single strand of nucleotides, but the nucleotides can also pair up in specific combinations to form sections of double-stranded RNA. Therefore, a single strand of non-coding RNA can fold into a complex three-dimensional shape that contains loops, twists, and bulges.

The three-dimensional structures of non-coding RNAs are crucial for their roles in cells, but the variety and complexity of shapes that they can form makes it technically difficult to study them. In 2008, researchers developed a new method called MOHCA that can map the positions of nucleotides that are close together in the three-dimensional structure. Highly reactive chemicals are attached to the nucleotides and these can react with, and damage, other nearby nucleotides. By detecting which nucleotides have been damaged, it is possible to map the positions of these nucleotides and decipher the structure of the RNA molecule using computer algorithms.

MOHCA is a promising approach, but the initial methods to find the damaged nucleotides were tedious and required specialized equipment. Now, Cheng, Das et al.—including some of the researchers involved in the 2008 work—have developed an improved version of MOHCA that uses readily available RNA sequencing techniques to find the damaged nucleotides. The RNA sequencing data are then analyzed by a new algorithm in the Rosetta computer modeling software.

Cheng, Das et al. used this newly developed ‘MOHCA-seq’ and Rosetta to reveal the structures of a human non-coding RNA and several other non-coding RNA molecules to a much higher level of detail than before. Together, MOHCA-seq and Rosetta provide a rapid method for researchers to decipher the three-dimensional structure of non-coding RNAs. This method is likely to speed up the analysis of the complex structures of non-coding RNAs. It will be useful in future efforts to work out what roles these RNAs play in cells, including their activity in cancer, neurodegeneration, and other diseases.

DOI:http://dx.doi.org/10.7554/eLife.07600.002

Predicting RNA 3D structure using a coarse-grain helix-centered model

A 3D model of RNA structure can provide information about its function and regulation that is not possible with just the sequence or secondary structure. Current models suffer from low accuracy and long running times and either neglect or presume knowledge of the long-range interactions which stabilize the tertiary structure. Our coarse-grained, helix-based, tertiary structure model operates with only a few degrees of freedom compared with all-atom models while preserving the ability to sample tertiary structures given a secondary structure. It strikes a balance between the precision of an all-atom tertiary structure model and the simplicity and effectiveness of a secondary structure representation. It provides a simplified tool for exploring global arrangements of helices and loops within RNA structures. We provide an example of a novel energy function relying only on the positions of stems and loops. We show that coupling our model to this energy function produces predictions as good as or better than the current state of the art tools. We propose that given the wide range of conformational space that needs to be explored, a coarse-grain approach can explore more conformations in less iterations than an all-atom model coupled to a fine-grain energy function. Finally, we emphasize the overarching theme of providing an ensemble of predicted structures, something which our tool excels at, rather than providing a handful of the lowest energy structures.

Emerging applications of small angle solution scattering in structural biology

Small angle solution X-ray and neutron scattering recently resurfaced as powerful tools to address an array of biological problems including folding, intrinsic disorder, conformational transitions, macromolecular crowding, and self or hetero-assembling of biomacromolecules. In addition, small angle solution scattering complements crystallography, nuclear magnetic resonance spectroscopy, and other structural methods to aid in the structure determinations of multidomain or multicomponent proteins or nucleoprotein assemblies. Neutron scattering with hydrogen/deuterium contrast variation, or X-ray scattering with sucrose contrast variation to a certain extent, is a convenient tool for characterizing the organizations of two-component systems such as a nucleoprotein or a lipid-protein assembly. Time-resolved small and wide-angle solution scattering to study biological processes in real time, and the use of localized heavy-atom labeling and anomalous solution scattering for applications as FRET-like molecular rulers, are amongst promising newer developments. Despite the challenges in data analysis and interpretation, these X-ray/neutron solution scattering based approaches hold great promise for understanding a wide variety of complex processes prevalent in the biological milieu.

Global analysis of riboswitches by small-angle X-ray scattering and calorimetry

Riboswitches are phylogenetically widespread non-coding mRNA domains that directly bind cellular metabolites and regulate transcription, translation, RNA stability or splicing via alternative RNA structures modulated by ligand binding. The details of ligand recognition by many riboswitches have been elucidated using X-ray crystallography and NMR. However, the global dynamics of riboswitch-ligand interactions and their thermodynamic driving forces are less understood. By compiling the work of many laboratories investigating riboswitches using small-angle X-ray scattering (SAXS) and isothermal titration calorimetry (ITC), we uncover general trends and common themes. There is a pressing need for community-wide consensus experimental conditions to allow results of riboswitch studies to be compared rigorously. Nonetheless, our meta-analysis reveals considerable diversity in the extent to which ligand binding reorganizes global riboswitch structures. It also demonstrates a wide spectrum of enthalpy-entropy compensation regimes across riboswitches that bind a diverse set of ligands, giving rise to a relatively narrow range of physiologically relevant free energies and ligand affinities. From the strongly entropy-driven binding of glycine to the predominantly enthalpy-driven binding of c-di-GMP to their respective riboswitches, these distinct thermodynamic signatures reflect the versatile strategies employed by RNA to adapt to the chemical natures of diverse ligands. Riboswitches have evolved to use a combination of long-range tertiary interactions, conformational selection, and induced fit to work with distinct ligand structure, charge, and solvation properties. This article is part of a Special Issue entitled: Riboswitches.

Analysis of riboswitch structure and ligand binding using small-angle X-ray scattering (SAXS)

Small-angle X-ray scattering (SAXS) is a powerful tool for examining the global conformation of riboswitches in solution, and how this is modulated by binding of divalent cations and small molecule ligands. SAXS experiments, which typically require only minutes per sample, directly yield two quantities describing the size and shape of the RNA: the radius of gyration (Rg) and the maximum linear dimension (Dmax). Examination of these quantities can reveal if a riboswitch undergoes cation-induced compaction. Comparison of the Rg and Dmax values between samples containing different concentrations of ligand reveals the overall structural response of the riboswitch to ligand. The Kratky plot (a graphical representation that emphasizes the higher-resolution SAXS data) and the P(r) plot or pair-probability distribution (an indirect Fourier transform, or power spectrum of the data) can provide additional evidence of riboswitch conformational changes. Simulation methods have been developed for generating three-dimensional reconstructions consistent with the one-dimensional SAXS data. These low-resolution molecular envelopes can aid in deciphering the relative helical arrangement within the RNA.

A personal perspective on chemistry-driven RNA research

In this mini review, we discuss how our understanding of ribonucleic acid (RNA) properties becomes significantly deepened when a broad range of modern chemical and biophysical methods is applied. We span our perspective from RNA solid-phase synthesis and site-specific labeling to single-molecule fluorescence-resonance-energy-transfer imaging and NMR spectroscopy approaches to explore the dynamics of RNA over a broad timescale. We then move on to Fourier-transform-ion-cyclotron-resonance mass spectrometry (FT-ICR-MS) as a powerful technique for RNA sequencing and modification analysis. The novel methodological developments are discussed for selected biological systems that include the thiamine-pyrophosphate riboswitch, HIV and ribosomal A-site RNA, and transfer RNA.

Wide-angle X-ray solution scattering for protein-ligand binding: multivariate curve resolution with Bayesian confidence intervals

A new way to use wide-angle x-ray solution scattering to study protein-ligand binding is presented. First, scattering patterns are measured at different protein and ligand concentrations. Multivariate curve resolution based on singular value decomposition and global analysis is applied to estimate the binding affinities and reference patterns (i.e., the scattering patterns of individual components). As validated by simulation, Bayesian confidence intervals provide accurate uncertainty estimates for the binding free energies and reference patterns. Experimental results from several protein-ligand systems demonstrate the feasibility of the approach, which promises to expand the role of wide-angle x-ray scattering as a quantitative biophysical tool.

Nucleic acid structure characterization by small angle X-ray scattering (SAXS)

Small angle X-ray scattering (SAXS) is a powerful method for investigating macromolecular structure in solution. SAXS data provide information about the size and shape of a molecule with a resolution of ∼2 to 3 nm. SAXS is particularly useful for the investigation of nucleic acids, which scatter X-rays strongly due to the electron-rich phosphate backbone. Therefore, SAXS has become an increasingly popular method for modeling nucleic acid structures, an endeavor made tractable by the highly regular helical nature of nucleic acid secondary structures. Recently, SAXS was used in combination with NMR to filter and refine all-atom models of a U2/U6 small nuclear RNA complex. In this unit, general protocols for sample preparation, data acquisition, and data analysis and processing are given. Additionally, examples of correctly and incorrectly processed SAXS data and expected results are provided.

Metabolite recognition principles and molecular mechanisms underlying riboswitch function

Riboswitches are mRNA elements capable of modulating gene expression in response to specific binding by cellular metabolites. Riboswitches exert their function through the interplay of alternative ligand-free and ligand-bound conformations of the metabolite-sensing domain, which in turn modulate the formation of adjacent gene expression controlling elements. X-ray crystallography and NMR spectroscopy have determined three-dimensional structures of virtually all the major riboswitch classes in the ligand-bound state and, for several riboswitches, in the ligand-free state. The resulting spatial topologies have demonstrated the wide diversity of riboswitch folds and revealed structural principles for specific recognition by cognate metabolites. The available three-dimensional information, supplemented by structure-guided biophysical and biochemical experimentation, has led to an improved understanding of how riboswitches fold, what RNA conformations are required for ligand recognition, and how ligand binding can be transduced into gene expression modulation. These studies have greatly facilitated the dissection of molecular mechanisms underlying riboswitch action and should in turn guide the anticipated development of tools for manipulating gene regulatory circuits.

Folding and ligand recognition of the TPP riboswitch aptamer at single-molecule resolution

Thiamine pyrophosphate (TPP)-sensitive mRNA domains are the most prevalent riboswitches known. Despite intensive investigation, the complex ligand recognition and concomitant folding processes in the TPP riboswitch that culminate in the regulation of gene expression remain elusive. Here, we used single-molecule fluorescence resonance energy transfer imaging to probe the folding landscape of the TPP aptamer domain in the absence and presence of magnesium and TPP. To do so, distinct labeling patterns were used to sense the dynamics of the switch helix (P1) and the two sensor arms (P2/P3 and P4/P5) of the aptamer domain. The latter structural elements make interdomain tertiary contacts (L5/P3) that span a region immediately adjacent to the ligand-binding site. In each instance, conformational dynamics of the TPP riboswitch were influenced by ligand binding. The P1 switch helix, formed by the 5' and 3' ends of the aptamer domain, adopts a predominantly folded structure in the presence of Mg(2+) alone. However, even at saturating concentrations of Mg(2+) and TPP, the P1 helix, as well as distal regions surrounding the TPP-binding site, exhibit an unexpected degree of residual dynamics and disperse kinetic behaviors. Such plasticity results in a persistent exchange of the P3/P5 forearms between open and closed configurations that is likely to facilitate entry and exit of the TPP ligand. Correspondingly, we posit that such features of the TPP aptamer domain contribute directly to the mechanism of riboswitch-mediated translational regulation.

The expression platform and the aptamer: cooperativity between Mg2+ and ligand in the SAM-I riboswitch

Riboswitch operation involves the complex interplay between the aptamer domain and the expression platform. During transcription, these two domains compete against each other for shared sequence. In this study, we explore the cooperative effects of ligand binding and Magnesium interactions in the SAM-I riboswitch in the context of aptamer collapse and anti-terminator formation. Overall, our studies show the apo-aptamer acts as (i) a pre-organized aptamer competent to bind ligand and undergo structural collapse and (ii) a conformation that is more accessible to anti-terminator formation. We show that both Mg(2+) ions and SAM are required for a collapse transition to occur. We then use competition between the aptamer and expression platform for shared sequence to characterize the stability of the collapsed aptamer. We find that SAM and Mg(2+) interactions in the aptamer are highly cooperative in maintaining switch polarity (i.e. aptamer 'off-state' versus anti-terminator 'on-state'). We further show that the aptamer off-state is preferentially stabilized by Mg(2+) and similar divalent ions. Furthermore, the functional switching assay was used to select for phosphorothioate interference, and identifies potential magnesium chelation sites while characterizing their coordinated role with SAM in aptamer stabilization. In addition, we find that Mg(2+) interactions with the apo-aptamer are required for the full formation of the anti-terminator structure, and that higher concentrations of Mg(2+) (>4 mM) shift the equilibrium toward the anti-terminator on-state even in the presence of SAM.

Fingerprinting noncanonical and tertiary RNA structures by differential SHAPE reactivity

Many RNA structures are composed of simple secondary structure elements linked by a few critical tertiary interactions. SHAPE chemistry has made interrogation of RNA dynamics at single-nucleotide resolution straightforward. However, de novo identification of nucleotides involved in tertiary interactions remains a challenge. Here we show that nucleotides that form noncanonical or tertiary contacts can be detected by comparing information obtained using two SHAPE reagents, N-methylisatoic anhydride (NMIA) and 1-methyl-6-nitroisatoic anhydride (1M6). Nucleotides that react preferentially with NMIA exhibit slow local nucleotide dynamics and usually adopt the less common C2'-endo ribose conformation. Experiments and first-principles calculations show that 1M6 reacts preferentially with nucleotides in which one face of the nucleobase allows an unhindered stacking interaction with the reagent. Differential SHAPE reactivities were used to detect noncanonical and tertiary interactions in four RNAs with diverse structures and to identify preformed noncanonical interactions in partially folded RNAs. Differential SHAPE reactivity analysis will enable experimentally concise, large-scale identification of tertiary structure elements and ligand binding sites in complex RNAs and in diverse biological environments.

Allosteric tertiary interactions preorganize the c-di-GMP riboswitch and accelerate ligand binding

Cyclic diguanylate (c-di-GMP) is a bacterial second messenger important for physiologic adaptation and virulence. Class-I c-di-GMP riboswitches are phylogenetically widespread and thought to mediate pleiotropic genetic responses to the second messenger. Previous studies suggest that the RNA aptamer domain switches from an extended free state to a compact, c-di-GMP-bound conformation in which two helical stacks dock side-by-side. Single molecule fluorescence resonance energy transfer (smFRET) experiments now reveal that the free RNA exists in four distinct populations that differ in dynamics in the extended and docked conformations. In the presence of c-di-GMP and Mg(2+), a stably docked population (>30 min) becomes predominant. smFRET mutant analysis demonstrates that tertiary interactions distal to the c-di-GMP binding site strongly modulate the RNA population structure, even in the absence of c-di-GMP. These allosteric interactions accelerate ligand recognition by preorganizing the RNA, favoring rapid c-di-GMP binding.

Molecular recognition and function of riboswitches

Regulatory mRNAs elements termed riboswitches respond to elevated concentrations of cellular metabolites by modulating expression of associated genes. Riboswitches attain their high metabolite selectivity by capitalizing on the intrinsic tertiary structures of their sensor domains. Over the years, riboswitch structure and folding have been amongst the most researched topics in the RNA field. Most recently, novel structures of single-ligand and cooperative double-ligand sensors have broadened our knowledge of architectural and molecular recognition principles exploited by riboswitches. The structural information has been complemented by extensive folding studies, which have provided several important clues on the formation of ligand-competent conformations and mechanisms of ligand discrimination. These studies have greatly improved our understanding of molecular events in riboswitch-mediated gene expression control and provided the molecular basis for intervention into riboswitch-controlled genetic circuits.

Probing riboswitch-ligand interactions using thiamine pyrophosphate analogues

The Escherichia coli thiM riboswitch forms specific contacts with its natural ligand, thiamine pyrophosphate (TPP or thiamine diphosphate), allowing it to generate not only nanomolar binding affinity, but also a high degree of discrimination against similar small molecules. A range of synthetic TPP analogues have been used to probe each of the riboswitch-ligand interactions. The results show that the pyrimidine-sensing helix of thiM is exquisitely tuned to select for TPP by recognising the H-bonding donor and acceptors around its aminopyrimidine ring and also by forming π-stacking interactions that may be sensitive to the electronics of the ring. The central thiazolium ring of TPP appears to be more important for ligand recognition than previously thought. It may contribute to binding via long-range electrostatic interactions and/or by exerting an electron withdrawing effect on the pyrimidine ring, allowing its presence to be sensed indirectly and thereby allowing discrimination between thiamine (and its phosphate esters) and other aminopyrimidines found in vivo. The pyrophosphate moiety is essential for submicromolar binding affinity, but unexpectedly, it does not appear to be strictly necessary for modulation of gene expression.

Multiple conformations of SAM-II riboswitch detected with SAXS and NMR spectroscopy

Riboswitches are a newly discovered large family of structured functional RNA elements that specifically bind small molecule targets out of a myriad of cellular metabolites to modulate gene expression. Structural studies of ligand-bound riboswitches by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy have provided insights into detailed RNA-ligand recognition and interactions. However, the structures of ligand-free riboswitches remain poorly characterized. In this study, we have used a variety of biochemical, biophysical and computational techniques including small-angle X-ray scattering and NMR spectroscopy to characterize the ligand-free and ligand-bound forms of SAM-II riboswitch. Our data demonstrate that the RNA adopts multiple conformations along its folding pathway and suggest that the RNA undergoes marked conformational changes upon Mg(2+) compaction and S-adenosylmethionine (SAM) metabolite binding. Further studies indicated that Mg(2+) ion is not essential for the ligand binding but can stabilize the complex by facilitating loop/stem interactions. In the presence of millimolar concentration of Mg(2+) ion, the RNA samples a more compact conformation. This conformation is near to, but distinct from, the native fold and competent to bind the metabolite. We conclude that the formation of various secondary and tertiary structural elements, including a pseudoknot, occur to sequester the putative Shine-Dalgarno sequence of the RNA only after metabolite binding.

Rationalization and prediction of selective decoding of pseudouridine-modified nonsense and sense codons

A stop or nonsense codon is an in-frame triplet within a messenger RNA that signals the termination of translation. One common feature shared among all three nonsense codons (UAA, UAG, and UGA) is a uridine present at the first codon position. It has been recently shown that the conversion of this uridine into pseudouridine (Ψ) suppresses translation termination, both in vitro and in vivo. Furthermore, decoding of the pseudouridylated nonsense codons is accompanied by the incorporation of two specific amino acids in a nonsense codon-dependent fashion. Ψ differs from uridine by a single N¹H group at the C5 position; how Ψ suppresses termination and, more importantly, enables selective decoding is poorly understood. Here, we provide molecular rationales for how pseudouridylated stop codons are selectively decoded. Our analysis applies crystal structures of ribosomes in varying states of translation to consider weakened interaction of Ψ with release factor; thermodynamic and geometric considerations of the codon-anticodon base pairs to rank and to eliminate mRNA-tRNA pairs; the mechanism of fidelity check of the codon-anticodon pairing by the ribosome to evaluate noncanonical codon-anticodon base pairs and the role of water. We also consider certain tRNA modifications that interfere with the Ψ-coordinated water in the major groove of the codon-anticodon mini-helix. Our analysis of nonsense codons enables prediction of potential decoding properties for Ψ-modified sense codons, such as decoding ΨUU potentially as Cys and Tyr. Our results provide molecular rationale for the remarkable dynamics of ribosome decoding and insights on possible reprogramming of the genetic code using mRNA modifications.

Determining RNA three-dimensional structures using low-resolution data

Knowing the 3-D structure of an RNA is fundamental to understand its biological function. Nowadays X-ray crystallography and NMR spectroscopy are systematically applied to newly discovered RNAs. However, the application of these high-resolution techniques is not always possible, and thus scientists must turn to lower resolution alternatives. Here, we introduce a pipeline to systematically generate atomic resolution 3-D structures that are consistent with low-resolution data sets. We compare and evaluate the discriminative power of a number of low-resolution experimental techniques to reproduce the structure of the Escherichia coli tRNA(VAL) and P4-P6 domain of the Tetrahymena thermophila group I intron. We test single and combinations of the most accessible low-resolution techniques, i.e. hydroxyl radical footprinting (OH), methidiumpropyl-EDTA (MPE), multiplexed hydroxyl radical cleavage (MOHCA), and small-angle X-ray scattering (SAXS). We show that OH-derived constraints are accurate to discriminate structures at the atomic level, whereas EDTA-based constraints apply to global shape determination. We provide a guide for choosing which experimental techniques or combination of thereof is best in which context. The pipeline represents an important step towards high-throughput low-resolution RNA structure determination.

Folding energy landscape of the thiamine pyrophosphate riboswitch aptamer

Riboswitches are motifs in the untranslated regions (UTRs) of RNA transcripts that sense metabolite levels and modulate the expression of the corresponding genes for metabolite import, export, synthesis, or degradation. All riboswitches contain an aptamer: an RNA structure that, upon binding ligand, folds to expose or sequester regulatory elements in the adjacent sequence through alternative nucleotide pairing. The coupling between ligand binding and aptamer folding is central to the regulatory mechanisms of thiamine pyrophosphate (TPP) riboswitches and has not been fully characterized. Here, we show that TPP aptamer folding can be decomposed into ligand-independent and -dependent steps that correspond to the formation of secondary and tertiary structures, respectively. We reconstructed the full energy landscape for folding of the wild-type (WT) aptamer and measured perturbations of this landscape arising from mutations or ligand binding. We show that TPP binding proceeds in two steps, from a weakly to a strongly bound state. Our data imply a hierarchical folding sequence, and provide a framework for understanding molecular mechanism throughout the TPP riboswitch family.

Riboswitches: discovery of drugs that target bacterial gene-regulatory RNAs

Riboswitches are messenger RNA (mRNA) domains that regulate gene function in response to the intracellular concentration of a variety of metabolites and second messengers. They control essential genes in many pathogenic bacteria, thus representing an inviting new class of biomolecular target for the development of antibiotics and chemical-biological tools. In this Account, we briefly review the discovery of riboswitches in the first years of the 21st century and their ensuing characterization over the past decade. We then discuss the progress achieved so far in using riboswitches as a focus for drug discovery, considering both the value of past serendipity and the particular challenges that confront current researchers. Five mechanisms of gene regulation have been determined for riboswitches. Most bacterial riboswitches modulate either transcription termination or translation initiation in response to ligand binding. All known examples of eukaryotic riboswitches, and some bacterial riboswitches, control gene expression by alternative splicing. The glmS riboswitch, which is widespread in Gram-positive bacteria, is a catalytic RNA activated by ligand binding: its self-cleavage destabilizes the mRNA of which it is part. Finally, one example of a trans-acting riboswitch is known. Three-dimensional structures have been determined for representatives of 13 structurally distinct riboswitch classes, providing atomic-level insight into their mechanisms of ligand recognition. While cellular and viral RNAs have attracted widespread interest as potential drug targets, riboswitches show special promise due to the diversity of small-molecule recognition strategies that are on display in their ligand-binding pockets. Moreover, riboswitches have evolved to recognize small-molecule ligands, which is unique among known structured RNA domains. Structural and biochemical advances in the study of riboswitches provide an impetus for the development of methods for the discovery of novel riboswitch activators and inhibitors. Recent rational drug design efforts focused on select riboswitch classes have yielded a small number of candidate antibiotic compounds, including one active in a mouse model of Staphylococcus aureus infection. The development of high-throughput methods suitable for riboswitch-specific drug discovery is ongoing. A fragment-based screening approach employing equilibrium dialysis that may be generically useful has demonstrated early success. Riboswitch-mediated gene regulation is widely employed by bacteria; however, only the thiamine pyrophosphate-responsive riboswitch has thus far been found in eukaryotes. Thus, riboswitches are particularly attractive as targets for antibacterials. Indeed, antimicrobials with previously unknown mechanisms have been found to function by binding riboswitches and causing aberrant gene expression.

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

Understanding the nature of the free state of riboswitch aptamers is important for illuminating common themes in gene regulation by riboswitches. Prior evidence indicated the flavin mononucleotide (FMN)-binding riboswitch aptamer adopted a 'bound-like' structure in absence of FMN, suggesting only local conformational changes upon ligand binding. In the scope of pinpointing the general nature of such changes at the nucleotide level, we performed SHAPE mapping experiments using the aptamer domain of two phylogenetic variants, both in absence and in presence of FMN. We also solved the crystal structures of one of these domains both free (3.3 Å resolution) and bound to FMN (2.95 Å resolution). Our comparative study reveals that structural rearrangements occurring upon binding are restricted to a few of the joining regions that form the binding pocket in both RNAs. This type of binding event with minimal structural perturbations is reminiscent of binding events by conformational selection encountered in other riboswitches and various RNAs.

Atomic-scale characterization of conformational changes in the preQ₁ riboswitch aptamer upon ligand binding

Riboswitches are mRNA structural elements that act as intracellular sensors of small-molecule metabolites. By undergoing conformational changes capable of modulating translation or terminating transcription, riboswitches are able to play a role in regulating the concentration of essential metabolites in the cell. Computer-guided fluorescence experiments were carried out to interrogate molecular dynamics and conformational changes in the minimal riboswitch aptamer that binds 7-aminomethyl-7-deazaguanine (preQ₁). Our combined experimental results and computational analysis suggest that the preQ₁ riboswitch apo form is structured but shows no evidence of a ligand-binding pocket. Simulations of the apo and bound forms indicate a large conformational change is triggered by the breaking of the Watson-Crick base pairing of nucleotides G11 and C31 upon preQ₁ removal, followed by collapse of the pocket due to interfering π-stacking. Computational predictions of local aptamer dynamics were validated by fluorescence experiments employing 2-aminopurine substitutions. In-line probing reactions confirmed that fluorophore-labeled riboswitches retain similar higher-order structural features as the unlabeled aptamer upon ligand binding, although their affinity for the ligand is reduced by the introduction of the fluorescent reporter.

Comparison of a preQ1 riboswitch aptamer in metabolite-bound and free states with implications for gene regulation

Riboswitches are RNA regulatory elements that govern gene expression by recognition of small molecule ligands via a high affinity aptamer domain. Molecular recognition can lead to active or attenuated gene expression states by controlling accessibility to mRNA signals necessary for transcription or translation. Key areas of inquiry focus on how an aptamer attains specificity for its effector, the extent to which the aptamer folds prior to encountering its ligand, and how ligand binding alters expression signal accessibility. Here we present crystal structures of the preQ(1) riboswitch from Thermoanaerobacter tengcongensis in the preQ(1)-bound and free states. Although the mode of preQ(1) recognition is similar to that observed for preQ(0), surface plasmon resonance revealed an apparent K(D) of 2.1 ± 0.3 nm for preQ(1) but a value of 35.1 ± 6.1 nm for preQ(0). This difference can be accounted for by interactions between the preQ(1) methylamine and base G5 of the aptamer. To explore conformational states in the absence of metabolite, the free-state aptamer structure was determined. A14 from the ceiling of the ligand pocket shifts into the preQ(1)-binding site, resulting in "closed" access to the metabolite while simultaneously increasing exposure of the ribosome-binding site. Solution scattering data suggest that the free-state aptamer is compact, but the "closed" free-state crystal structure is inadequate to describe the solution scattering data. These observations are distinct from transcriptional preQ(1) riboswitches of the same class that exhibit strictly ligand-dependent folding. Implications for gene regulation are discussed.

Comparison of solution and crystal structures of preQ1 riboswitch reveals calcium-induced changes in conformation and dynamics

Riboswitches regulate gene expression via specific recognition of cognate metabolites by their aptamer domains, which fold into stable conformations upon ligand binding. However, the recently reported solution and crystal structures of the Bacillus subtilis preQ₁ riboswitch aptamer show small but significant differences, suggesting that there may be conformational heterogeneity in the ligand-bound state. We present a structural and dynamic characterization of this aptamer by solution NMR spectroscopy. The aptamer−preQ₁ complex is intrinsically flexible in solution, with two regions that undergo motions on different time scales. Three residues move in concert on the micro-to-millisecond time scale and may serve as the lid of the preQ₁-binding pocket. Several Ca²⁺ ions are present in the crystal structure, one of which binds with an affinity of 47 ± 2 μM in solution to a site that is formed only upon ligand binding. Addition of Ca²⁺ to the aptamer−preQ₁ complex in solution results in conformational changes that account for the differences between the solution and crystal structures. Remarkably, the Ca²⁺ ions present in the crystal structure, which were proposed to be important for folding and ligand recognition, are not required for either in solution.

Insights into metalloregulation by M-box riboswitch RNAs via structural analysis of manganese-bound complexes

The M-box riboswitch couples intracellular magnesium levels to expression of bacterial metal transport genes. Structural analyses on other riboswitch RNA classes, which typically respond to a small organic metabolite, have revealed that ligand recognition occurs through a combination of base-stacking, electrostatic, and hydrogen-bonding interactions. In contrast, the M-box RNA triggers a change in gene expression upon association with an undefined population of metals, rather than responding to only a single ligand. Prior biophysical experimentation suggested that divalent ions associate with the M-box RNA to promote a compacted tertiary conformation, resulting in sequestration of a short sequence tract otherwise required for downstream gene expression. Electrostatic shielding from loosely associated metals is undoubtedly an important influence during this metal-mediated compaction pathway. However, it is also likely that a subset of divalent ions specifically occupies cation binding sites and promotes proper positioning of functional groups for tertiary structure stabilization. To better elucidate the role of these metal binding sites, we resolved a manganese-chelated M-box RNA complex to 1.86 Å by X-ray crystallography. These data support the presence of at least eight well-ordered cation binding pockets, including several sites that had been predicted by biochemical studies but were not observed in prior structural analysis. Overall, these data support the presence of three metal-binding cores within the M-box RNA that facilitate a network of long-range interactions within the metal-bound, compacted conformation.

Clustering to identify RNA conformations constrained by secondary structure

RNA often folds hierarchically, so that its sequence defines its secondary structure (helical base-paired regions connected by single-stranded junctions), which subsequently defines its tertiary fold. To preserve base-pairing and chain connectivity, the three-dimensional conformations that RNA can explore are strongly confined compared to when secondary structure constraints are not enforced. Using three examples, we studied how secondary structure confines and dictates an RNA's preferred conformations. We made use of Macromolecular Conformations by SYMbolic programming (MC-Sym) fragment assembly to generate RNA conformations constrained by secondary structure. Then, to understand the correlations between different helix placements and orientations, we robustly clustered all RNA conformations by employing unique methods to remove outliers and estimate the best number of conformational clusters. We observed that the preferred conformation (as judged by largest cluster size) for each type of RNA junction molecule tested is consistent with its biological function. Further, the improved quality of models in our pruned datasets facilitates subsequent discrimination using scoring functions based either on statistical analysis (knowledge based) or experimental data.

Recognition of S-adenosylmethionine by riboswitches

Riboswitches are regulatory elements commonly found in the 5' leader sequences of bacterial mRNAs that bind cellular metabolites to direct expression at either the transcriptional or translational level. The effectors of these RNAs are chemically diverse, including nucleobases and nucleosides, amino acids, cofactors, and second messenger molecules. Over the last few years, a number of structures have revealed the architectural means by which RNA creates binding pockets of high affinity and specificity for these compounds. For most effectors, there is a single class of associated riboswitches. However, eight individual classes of S-adenosylmethionine (SAM) and/or S-adenosylhomocysteine (SAH) responsive riboswitches that control various aspects of sulfur metabolism have been validated, revealing a diverse set of solutions to the recognition of these ubiquitous metabolites. This review focuses upon the structures of RNAs that bind SAM and SAH and how they discriminate between these compounds.