Structural studies of the purine and SAM binding riboswitches

Structural basis for 2'-deoxyguanosine recognition by the 2'-dG-II class of riboswitches

A recent bioinformatic analysis of well-characterized classes of riboswitches uncovered subgroups unable to bind to the regulatory molecule of the parental class. Within the guanine/adenine class, seven groups of RNAs were identified that deviate from the consensus sequence at one or more of three positions directly involved purine nucleobase recognition, one of which was validated as a second class of 2'-deoxyguanosine riboswitch (called 2'-dG-II). To understand how 2'-dG-II riboswitches recognize their cognate ligand and how they differ from a previously identified class of 2'-deoxyguanosine binding riboswitches, we have solved the crystal structure of a 2'-dG-II aptamer domain bound to 2'-deoxyguanosine. This structure reveals a global architecture similar to other members of the purine riboswitch family, but contains key differences within the ligand binding core. Defining the 2'-dG-II riboswitches is a two-nucleotide insertion in the three-way junction that promotes novel base-base interactions. Unlike 2'-dG-I riboswitches, the 2'-dG-II class only requires local changes to the ligand binding pocket of the guanine/adenine class to achieve a change in ligand preference. Notably, members of the 2'-dG-II family have variable ability to discriminate between 2'-deoxyguanosine and riboguanosine, suggesting that a subset of 2'-dG-II riboswitches may bind either molecule to regulate gene expression.

The structure of the SAM/SAH-binding riboswitch

S-adenosylmethionine (SAM) is a central metabolite since it is used as a methyl group donor in many different biochemical reactions. Many bacteria control intracellular SAM concentrations using riboswitch-based mechanisms. A number of structurally different riboswitch families specifically bind to SAM and mainly regulate the transcription or the translation of SAM-biosynthetic enzymes. In addition, a highly specific riboswitch class recognizes S-adenosylhomocysteine (SAH)—the product of SAM-dependent methyl group transfer reactions—and regulates enzymes responsible for SAH hydrolysis. High-resolution structures are available for many of these riboswitch classes and illustrate how they discriminate between the two structurally similar ligands SAM and SAH. The so-called SAM/SAH riboswitch class binds both ligands with similar affinities and is structurally not yet characterized. Here, we present a high-resolution nuclear magnetic resonance structure of a member of the SAM/SAH-riboswitch class in complex with SAH. Ligand binding induces pseudoknot formation and sequestration of the ribosome binding site. Thus, the SAM/SAH-riboswitches are translational ‘OFF’-switches. Our results establish a structural basis for the unusual bispecificity of this riboswitch class. In conjunction with genomic data our structure suggests that the SAM/SAH-riboswitches might be an evolutionary late invention and not a remnant of a primordial RNA-world as suggested for other riboswitches.

Small-Molecule-Binding Riboswitches

RNA is a versatile biomolecule capable of transferring information, taking on distinct three-dimensional shapes, and reacting to ambient conditions. RNA molecules utilize a wide range of mechanisms to control gene expression. An example of such regulation is riboswitches. Consisting exclusively of RNA, they are able to control important metabolic processes, thus providing an elegant and efficient RNA-only regulation system. Existing across all domains of life, riboswitches appear to represent one of the most highly conserved mechanisms for the regulation of a broad range of biochemical pathways. Through binding of a wide range of small-molecule ligands to their so-called aptamer domain, riboswitches undergo a conformational change in their downstream "expression platform." In consequence, the pattern of gene expression changes, which in turn results in increased or decreased protein production. Riboswitches unite the sensing and transduction of a signal that can directly be coupled to the metabolism of the cell; thus they constitute a very potent regulatory mechanism for many organisms. Highly specific RNA-binding domains not only occur in vivo but can also be evolved by means of the SELEX (systematic evolution of ligands by exponential enrichment) method, which allows in vitro selection of aptamers against almost any ligand. Coupling of these aptamers with an expression platform has led to the development of synthetic riboswitches, a highly active research field of great relevance and immense potential. The aim of this review is to summarize developments in the riboswitch field over the last decade and address key questions of recent research.

Ribo-attenuators: novel elements for reliable and modular riboswitch engineering

Riboswitches are structural genetic regulatory elements that directly couple the sensing of small molecules to gene expression. They have considerable potential for applications throughout synthetic biology and bio-manufacturing as they are able to sense a wide range of small molecules and regulate gene expression in response. Despite over a decade of research they have yet to reach this considerable potential as they cannot yet be treated as modular components. This is due to several limitations including sensitivity to changes in genetic context, low tunability, and variability in performance. To overcome the associated difficulties with riboswitches, we have designed and introduced a novel genetic element called a ribo-attenuator in Bacteria. This genetic element allows for predictable tuning, insulation from contextual changes, and a reduction in expression variation. Ribo-attenuators allow riboswitches to be treated as truly modular and tunable components, thus increasing their reliability for a wide range of applications.

A magnesium-induced triplex pre-organizes the SAM-II riboswitch

Our 13C- and 1H-chemical exchange saturation transfer (CEST) experiments previously revealed a dynamic exchange between partially closed and open conformations of the SAM-II riboswitch in the absence of ligand. Here, all-atom structure-based molecular simulations, with the electrostatic effects of Manning counter-ion condensation and explicit magnesium ions are employed to calculate the folding free energy landscape of the SAM-II riboswitch. We use this analysis to predict that magnesium ions remodel the landscape, shifting the equilibrium away from the extended, partially unfolded state towards a compact, pre-organized conformation that resembles the ligand-bound state. Our CEST and SAXS experiments, at different magnesium ion concentrations, quantitatively confirm our simulation results, demonstrating that magnesium ions induce collapse and pre-organization. Agreement between theory and experiment bolsters microscopic interpretation of our simulations, which shows that triplex formation between helix P2b and loop L1 is highly sensitive to magnesium and plays a key role in pre-organization. Pre-organization of the SAM-II riboswitch allows rapid detection of ligand with high selectivity, which is important for biological function.

Phosphoribosyl Diphosphate (PRPP): Biosynthesis, Enzymology, Utilization, and Metabolic Significance

Phosphoribosyl diphosphate (PRPP) is an important intermediate in cellular metabolism. PRPP is synthesized by PRPP synthase, as follows: ribose 5-phosphate + ATP → PRPP + AMP. PRPP is ubiquitously found in living organisms and is used in substitution reactions with the formation of glycosidic bonds. PRPP is utilized in the biosynthesis of purine and pyrimidine nucleotides, the amino acids histidine and tryptophan, the cofactors NAD and tetrahydromethanopterin, arabinosyl monophosphodecaprenol, and certain aminoglycoside antibiotics. The participation of PRPP in each of these metabolic pathways is reviewed. Central to the metabolism of PRPP is PRPP synthase, which has been studied from all kingdoms of life by classical mechanistic procedures. The results of these analyses are unified with recent progress in molecular enzymology and the elucidation of the three-dimensional structures of PRPP synthases from eubacteria, archaea, and humans. The structures and mechanisms of catalysis of the five diphosphoryltransferases are compared, as are those of selected enzymes of diphosphoryl transfer, phosphoryl transfer, and nucleotidyl transfer reactions. PRPP is used as a substrate by a large number phosphoribosyltransferases. The protein structures and reaction mechanisms of these phosphoribosyltransferases vary and demonstrate the versatility of PRPP as an intermediate in cellular physiology. PRPP synthases appear to have originated from a phosphoribosyltransferase during evolution, as demonstrated by phylogenetic analysis. PRPP, furthermore, is an effector molecule of purine and pyrimidine nucleotide biosynthesis, either by binding to PurR or PyrR regulatory proteins or as an allosteric activator of carbamoylphosphate synthetase. Genetic analyses have disclosed a number of mutants altered in the PRPP synthase-specifying genes in humans as well as bacterial species.

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.

Dynamics of riboswitches: Molecular simulations

Riboswitch RNAs play key roles in bacterial metabolism and represent a promising new class of antibiotic targets for treatment of infectious disease. While many studies of riboswitches have been performed, the exact mechanism of riboswitch operation is still not fully understood at the atomistic level of detail. Molecular dynamics simulations are useful for interpreting existing experimental data and producing predictions for new experiments. Here, a wide range of computational studies on riboswitches is reviewed. By elucidating the key principles of riboswitch operation, computation may aid in the effort to design more specific antibiotics with affinities greater than those of the native ligand. Such a detailed understanding may be required to improve efficacy and reduce side effects. These studies are laying the groundwork for understanding the action mechanism of new compounds that inhibit riboswitch activity. Future directions such as magnesium effects, large-scale conformational changes, expression platforms and co-transcriptional folding are also discussed. This article is part of a Special Issue entitled: Riboswitches.

Modularity of select riboswitch expression platforms enables facile engineering of novel genetic regulatory devices

RNA-based biosensors and regulatory devices have received significant attention for their potential in a broad array of synthetic biology applications. One of the primary difficulties in engineering these molecules is the lack of facile methods to link sensory modules, or aptamers, to readout domains. Such efforts typically require extensive screening or selection of sequences that facilitate interdomain communication. Bacteria have evolved a widespread form of gene regulation known as riboswitches that perform this task with sufficient fidelity to control expression of biosynthetic and transport proteins essential for normal cellular homeostasis. In this work, we demonstrate that select riboswitch readout domains, called expression platforms, are modular in that they can host a variety of natural and synthetic aptamers to create novel chimeric RNAs that regulate transcription both in vitro and in vivo. Importantly, this technique does not require selection of device-specific "communication modules" required to transmit ligand binding to the regulatory domain, enabling rapid engineering of novel functional RNAs.

Evolutionary evidence for alternative structure in RNA sequence co-variation

Sequence conservation and co-variation of base pairs are hallmarks of structured RNAs. For certain RNAs (e.g. riboswitches), a single sequence must adopt at least two alternative secondary structures to effectively regulate the message. If alternative secondary structures are important to the function of an RNA, we expect to observe evolutionary co-variation supporting multiple conformations. We set out to characterize the evolutionary co-variation supporting alternative conformations in riboswitches to determine the extent to which alternative secondary structures are conserved. We found strong co-variation support for the terminator, P1, and anti-terminator stems in the purine riboswitch by extending alignments to include terminator sequences. When we performed Boltzmann suboptimal sampling on purine riboswitch sequences with terminators we found that these sequences appear to have evolved to favor specific alternative conformations. We extended our analysis of co-variation to classic alignments of group I/II introns, tRNA, and other classes of riboswitches. In a majority of these RNAs, we found evolutionary evidence for alternative conformations that are compatible with the Boltzmann suboptimal ensemble. Our analyses suggest that alternative conformations are selected for and thus likely play functional roles in even the most structured of RNAs.

RNA (Ribonucleic Acid) is a messenger of genetic information, master regulator, and catalyst in the cell. To carry out its function, RNA can fold into complex three-dimensional structures. Certain classes of RNAs, called riboswitches, adopt at least two alternative structures to act as a switch. We set out to detect the evolutionary signal for alternative structures in riboswitches as we hypothesize that these RNA sequences must have evolved to allow both conformations. We find that indeed such signals exist when we compare the sequences of riboswitches from multiple species. When we extend this analysis to other RNA regulators in the cell that are not thought of as switches, we detect equivalent evolutionary support for alternative structures. Viewed through the lens of evolutionary structure conservation RNA sequences appear to have adapted to adopt multiple conformations.

Structural and functional studies of S-adenosyl-L-methionine binding proteins: a ligand-centric approach

BACKGROUND

The post-genomic era poses several challenges. The biggest is the identification of biochemical function for protein sequences and structures resulting from genomic initiatives. Most sequences lack a characterized function and are annotated as hypothetical or uncharacterized. While homology-based methods are useful, and work well for sequences with sequence identities above 50%, they fail for sequences in the twilight zone (<30%) of sequence identity. For cases where sequence methods fail, structural approaches are often used, based on the premise that structure preserves function for longer evolutionary time-frames than sequence alone. It is now clear that no single method can be used successfully for functional inference. Given the growing need for functional assignments, we describe here a systematic new approach, designated ligand-centric, which is primarily based on analysis of ligand-bound/unbound structures in the PDB. Results of applying our approach to S-adenosyl-L-methionine (SAM) binding proteins are presented.

RESULTS

Our analysis included 1,224 structures that belong to 172 unique families of the Protein Information Resource Superfamily system. Our ligand-centric approach was divided into four levels: residue, protein/domain, ligand, and family levels. The residue level included the identification of conserved binding site residues based on structure-guided sequence alignments of representative members of a family, and the identification of conserved structural motifs. The protein/domain level included structural classification of proteins, Pfam domains, domain architectures, and protein topologies. The ligand level included ligand conformations, ribose sugar puckering, and the identification of conserved ligand-atom interactions. The family level included phylogenetic analysis.

CONCLUSION

We found that SAM bound to a total of 18 different fold types (I-XVIII). We identified 4 new fold types and 11 additional topological arrangements of strands within the well-studied Rossmann fold Methyltransferases (MTases). This extends the existing structural classification of SAM binding proteins. A striking correlation between fold type and the conformation of the bound SAM (classified as types) was found across the 18 fold types. Several site-specific rules were created for the assignment of functional residues to families and proteins that do not have a bound SAM or a solved structure.

Label-free microfluidic characterization of temperature-dependent biomolecular interactions

We present a microfluidic approach to characterizing temperature-dependent biomolecular interactions. Solvated L-arginine vasopressin (AVP) and its immobilized RNA aptamer (spiegelmer) were allowed to achieve equilibrium binding in a microchip at a series of selected temperatures. Unbound AVP were collected and analyzed with matrix-assisted laser desorption∕ionization mass spectrometry (MALDI-MS), yielding melting curves that reveal highly temperature-dependent zones in which affinity binding (36-45 °C) or dissociation (25-33 °C and 50-65 °C) occurs. Additionally, temperature-dependent binding isotherms were constructed; from these, thermodynamic quantities involved in binding were extracted. The results illustrated a strong change in heat capacity of interaction for this system, suggesting a considerable thermodynamic influence controlling vasopressin-spiegelmer interaction.

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.

Discrimination between closely related cellular metabolites by the SAM-I riboswitch

The SAM-I riboswitch is a cis-acting element of genetic control found in bacterial mRNAs that specifically binds S-adenosylmethionine (SAM). We previously determined the 2.9-A X-ray crystal structure of the effector-binding domain of this RNA element, revealing details of RNA-ligand recognition. To improve this structure, variations were made to the RNA sequence to alter lattice contacts, resulting in a 0.5-A improvement in crystallographic resolution and allowing for a more accurate refinement of the crystallographic model. The basis for SAM specificity was addressed by a structural analysis of the RNA complexed to S-adenosylhomocysteine (SAH) and sinefungin and by measuring the affinity of SAM and SAH for a series of mutants using isothermal titration calorimetry. These data illustrate the importance of two universally conserved base pairs in the RNA that form electrostatic interactions with the positively charged sulfonium group of SAM, thereby providing a basis for discrimination between SAM and SAH.

Monitoring RNA-ligand interactions using isothermal titration calorimetry

Isothermal titration calorimetry (ITC) is a biophysical technique that measures the heat evolved or absorbed during a reaction to report the enthalpy, entropy, stoichiometry of binding, and equilibrium association constant. A significant advantage of ITC over other methods is that it can be readily applied to almost any RNA-ligand complex without having to label either molecule and can be performed under a broad range of pH, temperature, and ionic concentrations. During our application of ITC to investigate the thermodynamic details of the interaction of a variety of compounds with the purine riboswitch, we have explored and optimized experimental parameters that yield the most useful and reproducible results for RNAs. In this chapter, we detail this method using the titration of an adenine-binding RNA with 2,6-diaminopurine (DAP) as a practical example. Our insights should be generally applicable to observing the interactions of a broad range of molecules with structured RNAs.

Determining structures of RNA aptamers and riboswitches by X-ray crystallography

Structural biology plays a central role in gaining a full understanding of the myriad roles of RNA in biology. In recent years, innovative approaches in RNA purification and crystallographic methods have lead to the visualization of an increasing number of unique structures, providing new insights into its function at the atomic level. This article presents general protocols which have streamlined the process of obtaining a homogeneous sample of properly folded and active RNA in high concentrations that crystallizes well in the presence of a suitable heavy-atom for phasing. Of particular importance are approaches toward RNA crystallography that include exploring "construct space" as opposed to "condition space". Moreover, development of a highly flexible method for experimentally phasing RNA crystals may open the door to a relatively simple means of solving these structures.

Riboswitches: emerging themes in RNA structure and function

Riboswitches are RNAs capable of binding cellular metabolites using a diverse array of secondary and tertiary structures to modulate gene expression. The recent determination of the three-dimensional structures of parts of six different riboswitches illuminates common features that allow riboswitches to be grouped into one of two types. Type I riboswitches, as exemplified by the purine riboswitch, are characterized by a single, localized binding pocket supported by a largely pre-established global fold. This arrangement limits ligand-induced conformational changes in the RNA to a small region. In contrast, Type II riboswitches, such as the thiamine pyrophosphate riboswitch, contain binding pockets split into at least two spatially distinct sites. As a result, binding induces both local changes to the binding pocket and global architecture. Similar organizational themes are found in other noncoding RNAs, making it possible to begin to build a hierarchical classification of RNA structure based on the spatial organization of their active sites and associated secondary structural elements.

Ligand-dependent folding of the three-way junction in the purine riboswitch

Riboswitches are highly structured cis-acting elements located in the 5'-untranslated region of messenger RNAs that directly bind small molecule metabolites to regulate gene expression. Structural and biochemical studies have revealed riboswitches experience significant ligand-dependent conformational changes that are coupled to regulation. To monitor the coupling of ligand binding and RNA folding within the aptamer domain of the purine riboswitch, we have chemically probed the RNA with N-methylisatoic anhydride (NMIA) over a broad temperature range. Analysis of the temperature-dependent reactivity of the RNA in the presence and absence of hypoxanthine reveals that a limited set of nucleotides within the binding pocket change their conformation in response to ligand binding. Our data demonstrate that a distal loop-loop interaction serves to restrict the conformational freedom of a significant portion of the three-way junction, thereby promoting ligand binding under physiological conditions.

Non-nearest-neighbor dependence of the stability for RNA bulge loops based on the complete set of group I single-nucleotide bulge loops

Fifty-nine RNA duplexes containing single-nucleotide bulge loops were optically melted in 1 M NaCl, and the thermodynamic parameters DeltaH degrees, DeltaS degrees, DeltaG 37 degrees, and TM for each sequence were determined. Sequences from this study were combined with sequences from previous studies [Longfellow, C. E., et al. (1990) Biochemistry 29, 278-285; Znosko, B. M., et al. (2002) Biochemistry 41, 10406-10417], thus examining all possible group I single-nucleotide bulge loop and nearest-neighbor sequence combinations. The free energy increments at 37 degrees C for the introduction of a group I single-nucleotide bulge loop range between 1.3 and 5.2 kcal/mol. The combined data were used to develop a model for predicting the free energy of a RNA duplex containing a single-nucleotide bulge. For bulge loops with adjacent Watson-Crick base pairs, neither the identity of the bulge nor the nearest-neighbor base pairs had an effect on the influence of the bulge loop on duplex stability. The proposed model for prediction of the stability of a duplex containing a bulged nucleotide was primarily affected by non-nearest-neighbor interactions. The destabilization of the duplex by the bulge was related to the stability of the stems adjacent to the bulge. Specifically, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. The stability of a duplex containing a bulged nucleotide adjacent to a wobble base pair also was primarily affected by non-nearest-neighbor interactions. Again, there was a direct correlation between the destabilization of the duplex and the stability of the less stable duplex stem. However, when one or both of the bulge nearest neighbors was a wobble base pair, the free energy increment for insertion of a bulge loop is dependent upon the position and orientation of the wobble base pair relative the bulged nucleotide. Bulge sequences of the type ((5'UBX)(3'GY)), ((5'GBG)(3'UU)) and ((5'UBU)(3'GG)) are less destabilizing by 0.6 kcal/mol, and bulge sequences of the type ((5'GBX)(3'UY)) and ((5'XBU)(3'YG)) are more destabilizing by 0.4 kcal/mol than bulge loops adjacent to Watson-Crick base pairs.

Mutational analysis of the purine riboswitch aptamer domain

The purine riboswitch is one of a number of mRNA elements commonly found in the 5'-untranslated region capable of controlling expression in a cis-fashion via its ability to directly bind small-molecule metabolites. Extensive biochemical and structural analysis of the nucleobase-binding domain of the riboswitch, referred to as the aptamer domain, has revealed that the mRNA recognizes its cognate ligand using an intricately folded three-way junction motif that completely encapsulates the ligand. High-affinity binding of the purine nucleobase is facilitated by a distal loop-loop interaction that is conserved between both the adenine and guanine riboswitches. To understand the contribution of conserved nucleotides in both the three-way junction and the loop-loop interaction of this RNA, we performed a detailed mutagenic survey of these elements in the context of an adenine-responsive variant of the xpt-pbuX guanine riboswitch from Bacillus subtilis. The varying ability of these mutants to bind ligand as measured by isothermal titration calorimetry uncovered the conserved nucleotides whose identity is required for purine binding. Crystallographic analysis of the bound form of five mutants and chemical probing of their free state demonstrate that the identity of several universally conserved nucleotides is not essential for formation of the RNA-ligand complex but rather for maintaining a binding-competent form of the free RNA. These data show that conservation patterns in riboswitches arise from a combination of formation of the ligand-bound complex, promoting an open form of the free RNA, and participating in the secondary structural switch with the expression platform.

A general strategy to solve the phase problem in RNA crystallography

X-ray crystallography of biologically important RNA molecules has been hampered by technical challenges, including finding heavy-atom derivatives to obtain high-quality experimental phase information. Existing techniques have drawbacks, limiting the rate at which important new structures are solved. To address this, we have developed a reliable means to localize heavy atoms specifically to virtually any RNA. By solving the crystal structures of thirteen variants of the G*U wobble pair cation binding motif, we have identified a version that when inserted into an RNA helix introduces a high-occupancy cation binding site suitable for phasing. This "directed soaking" strategy can be integrated fully into existing RNA crystallography methods, potentially increasing the rate at which important structures are solved and facilitating routine solving of structures using Cu-Kalpha radiation. This method already has been used to solve several crystal structures.

Structural features of metabolite-sensing riboswitches

Riboswitches, metabolite-sensing RNA elements that are present in untranslated regions of the transcripts that they regulate, possess extensive tertiary structure to couple metabolite binding to genetic control. Here we discuss recently published structures from four riboswitch classes and compare these natural RNA structures to those of in-vitro-selected RNA aptamers, which bind ligands similar to those of the riboswitches. In addition, we examine the glmS riboswitch - the first example of a ribozyme-based riboswitch. This RNA provides the latest twist in the riboswitch field and portends exciting advances in the coming years. Our knowledge of the mechanisms underlying genetic regulation by riboswitches has increased mightily in recent years and will continue to grow as new riboswitch classes and ligands are discovered and structurally characterized.