Riboswitch structure: an internal residue mimicking the purine ligand

The adenine and guanine riboswitches regulate gene expression in response to their purine ligand. X-ray structures of the aptamer moiety of these riboswitches are characterized by a compact fold in which the ligand forms a Watson–Crick base pair with residue 65. Phylogenetic analyses revealed a strict restriction at position 39 of the aptamer that prevents the G39–C65 and A39–U65 combinations, and mutational studies indicate that aptamers with these sequence combinations are impaired for ligand binding. In order to investigate the rationale for sequence conservation at residue 39, structural characterization of the U65C mutant from Bacillus subtilis pbuE adenine riboswitch aptamer was undertaken. NMR spectroscopy and X-ray crystallography studies demonstrate that the U65C mutant adopts a compact ligand-free structure, in which G39 occupies the ligand-binding site of purine riboswitch aptamers. These studies present a remarkable example of a mutant RNA aptamer that adopts a native-like fold by means of ligand mimicking and explain why this mutant is impaired for ligand binding. Furthermore, this work provides a specific insight into how the natural sequence has evolved through selection of nucleotide identities that contribute to formation of the ligand-bound state, but ensures that the ligand-free state remains in an active conformation.     

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

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

Metal-ion binding and metal-ion induced folding of the adenine-sensing riboswitch aptamer domain

Divalent cations are important in the folding and stabilization of complex RNA structures. The adenine-sensing riboswitch controls the expression of mRNAs for proteins involved in purine metabolism by directly sensing intracellular adenine levels. Adenine binds with high affinity and specificity to the ligand binding or aptamer domain of the adenine-sensing riboswitch. The X-ray structure of this domain in complex with adenine revealed an intricate RNA-fold consisting of a three-helix junction stabilized by long-range base-pairing interactions and identified five binding sites for hexahydrated Mg²⁺-ions. Furthermore, a role for Mg²⁺-ions in the ligand-induced folding of this RNA was suggested. Here, we describe the interaction of divalent cations with the RNA–adenine complex in solution as studied by high-resolution NMR spectroscopy. Paramagnetic line broadening, chemical shift mapping and intermolecular nuclear Overhauser effects (NOEs) indicate the presence of at least three binding sites for divalent cations. Two of them are similar to those in the X-ray structure. The third site, which is important for the folding of this RNA, has not been observed previously. The ligand-free state of the RNA is conformationally heterogeneous and contains base-pairing patterns detrimental to ligand binding in the absence of Mg²⁺, but becomes partially pre-organized for ligand binding in the presence of Mg²⁺. Compared to the highly similar guanine-sensing riboswitch, the folding pathway for the adenine-sensing riboswitch aptamer domain is more complex and the influence of Mg²⁺ is more pronounced.     

Thermodynamic and kinetic characterization of ligand binding to the purine riboswitch aptamer domain

Riboswitches are cis-acting genetic regulatory elements found commonly in bacterial mRNAs that consist of a metabolite-responsive aptamer domain coupled to a regulatory switch. Purine riboswitches respond to intracellular concentrations of either adenine or guanine/hypoxanthine to control gene expression. The aptamer domain of the purine riboswitch contains a pyrimidine residue (Y74) that forms a Watson-Crick base-pairing interaction with the bound purine nucleobase ligand that discriminates between adenine and guanine. We sought to understand the structural basis of this specificity and the mechanism of ligand recognition by the purine riboswitch. Here, we present the 2,6-diaminopurine-bound structure of a C74U mutant of the xpt-pbuX guanine riboswitch, along with a detailed thermodynamic and kinetic analysis of nucleobase recognition by both the native and mutant riboswitches. These studies demonstrate clearly that the pyrimidine at position 74 is the sole determinant of purine riboswitch specificity. In addition, the mutant riboswitch binds adenine and adenine derivatives well compared with the guanine-responsive riboswitch. Under our experimental conditions, 2,6-diaminopurine binds the RNA with DeltaH=-40.3 kcal mol(-1), DeltaS=-97.6 cal mol(-1)K(-1), and DeltaG=-10.73 kcal mol(-1). A kinetic determination of the slow rate (0.15 x 10(5)M(-1)s(-1) and 2.1 x 10(5)mM(-1)s(-1) for 2-aminopurine binding the adenine-responsive mutant riboswitch and 7-deazaguanine-binding guanine riboswitch, respectively) of association under varying experimental conditions allowed us to propose a mechanism for ligand recognition by the purine riboswitch. A conformationally dynamic unliganded state for the binding pocket is stabilized first by the Watson-Crick base pairing between the ligand and Y74, and by the subsequent ordering of the J2/3 loop, enclosing the ligand within the three-way junction.     

Guanidine-II aptamer conformations and ligand binding modes through the lens of molecular simulation

Regulation of gene expression via riboswitches is a widespread mechanism in bacteria. Here, we investigate ligand binding of a member of the guanidine sensing riboswitch family, the guanidine-II riboswitch (Gd-II). It consists of two stem–loops forming a dimer upon ligand binding. Using extensive molecular dynamics simulations we have identified conformational states corresponding to ligand-bound and unbound states in a monomeric stem–loop of Gd-II and studied the selectivity of this binding. To characterize these states and ligand-dependent conformational changes we applied a combination of dimensionality reduction, clustering, and feature selection methods. In absence of a ligand, the shape of the binding pocket alternates between the conformation observed in presence of guanidinium and a collapsed conformation, which is associated with a deformation of the dimerization interface. Furthermore, the structural features responsible for the ability to discriminate against closely related analogs of guanidine are resolved. Based on these insights, we propose a mechanism that couples ligand binding to aptamer dimerization in the Gd-II system, demonstrating the value of computational methods in the field of nucleic acids research.     

Binding free energy decomposition and multiple unbinding paths of buried ligands in a PreQ1 riboswitch

Riboswitches are naturally occurring RNA elements that control bacterial gene expression by binding to specific small molecules. They serve as important models for RNA-small molecule recognition and have also become a novel class of targets for developing antibiotics. Here, we carried out conventional and enhanced-sampling molecular dynamics (MD) simulations, totaling 153.5 μs, to characterize the determinants of binding free energies and unbinding paths for the cognate and synthetic ligands of a PreQ₁ riboswitch. Binding free energy analysis showed that two triplets of nucleotides, U6-C15-A29 and G5-G11-C16, contribute the most to the binding of the cognate ligands, by hydrogen bonding and by base stacking, respectively. Mg²⁺ ions are essential in stabilizing the binding pocket. For the synthetic ligands, the hydrogen-bonding contributions of the U6-C15-A29 triplet are significantly compromised, and the bound state resembles the apo state in several respects, including the disengagement of the C15-A14-A13 and A32-G33 base stacks. The bulkier synthetic ligands lead to significantly loosening of the binding pocket, including extrusion of the C15 nucleobase and a widening of the C15-C30 groove. Enhanced-sampling simulations further revealed that the cognate and synthetic ligands unbind in almost opposite directions. Our work offers new insight for designing riboswitch ligands.     

The yjdF riboswitch candidate regulates gene expression by binding diverse azaaromatic compounds

The yjdF motif RNA is an orphan riboswitch candidate that almost exclusively associates with the yjdF protein-coding gene in many bacteria. The function of the YjdF protein is unknown, which has made speculation regarding the natural ligand for this putative riboswitch unusually challenging. By using a structure-probing assay for ligand binding, we found that a surprisingly broad diversity of nitrogen-containing aromatic heterocycles, or “azaaromatics,” trigger near-identical changes in the structures adopted by representative yjdF motif RNAs. Regions of the RNA that undergo ligand-induced structural modulation reside primarily in portions of the putative aptamer region that are highly conserved in nucleotide sequence, as is typical for riboswitches. Some azaaromatic molecules are bound by the RNA with nanomolar dissociation constants, and a subset of these ligands activate riboswitch-mediated gene expression in cells. Furthermore, genetic elements most commonly adjacent to the yjdF motif RNA or to the yjdF protein-coding region are homologous to protein regulators implicated in mitigating the toxic effects of diverse phenolic acids or polycyclic compounds. Although the precise type of natural ligand sensed by yjdF motif RNAs remains unknown, our findings suggest that this riboswitch class might serve as part of a genetic response system to toxic or signaling compounds with chemical structures similar to azaaromatics.     

Structural insights into the interactions of xpt riboswitch with novel guanine analogues: a molecular dynamics simulation study

Ligand recognition in purine riboswitches is a complex process requiring different levels of conformational changes. Recent efforts in the area of purine riboswitch research have focused on ligand analogue binding studies. In the case of the guanine xanthine phosphoribosyl transferase (xpt) riboswitch, synthetic analogues that resemble guanine have the potential to tightly bind and subsequently influence the genetic expression of xpt mRNA in prokaryotes. We have carried out 25 ns Molecular Dynamics (MD) simulation studies of the aptamer domain of the xpt G-riboswitch in four different states: guanine riboswitch in free form, riboswitch bound with its cognate ligand guanine, and with two guanine analogues SJ1 and SJ2. Our work reveals novel interactions of SJ1 and SJ2 ligands with the binding core residues of the riboswitch. The ligands proposed in this work bind to the riboswitch with greater overall stability and lower root mean square deviations and fluctuations compared to guanine ligand. Reporter gene assay data demonstrate that the ligand analogues, upon binding to the RNA, lower the genetic expression of the guanine riboswitch. Our work has important implications for future ligand design and binding studies in the exciting field of riboswitches.     

Core requirements of the adenine riboswitch aptamer for ligand binding

The adenine riboswitch aptamer, the A box, positively regulates gene expression upon adenine binding. To provide insight into structure-function relationships, important for the adenine riboswitch aptamer, we have created alignments for six aptamer sequences that reveal the core requirements. In addition, 2-aminopurine (2AP) binding studies have been used to test the consensus sequence derived from the alignment. Overall, the consensus secondary structure is consistent with 2AP binding studies. However, a position in the core, previously identified as variable, shows restriction in nucleotide sequence. Furthermore, this restriction is found to be related with the ligand specificity of the riboswitch. The implications of this relationship for the riboswitch gene regulation mechanism are discussed.     

Riboswitch Detection Using Profile Hidden Markov Models

Background  

Riboswitches are a type of noncoding RNA that regulate gene expression by switching from one structural conformation to another on ligand binding. The various classes of riboswitches discovered so far are differentiated by the ligand, which on binding induces a conformational switch. Every class of riboswitch is characterized by an aptamer domain, which provides the site for ligand binding, and an expression platform that undergoes conformational change on ligand binding. The sequence and structure of the aptamer domain is highly conserved in riboswitches belonging to the same class. We propose a method for fast and accurate identification of riboswitches using profile Hidden Markov Models (pHMM). Our method exploits the high degree of sequence conservation that characterizes the aptamer domain.     

Conformational capture of the SAM-II riboswitch

Riboswitches are gene regulation elements in mRNA that function by specifically responding to metabolites. Although the metabolite-bound states of riboswitches have proven amenable to structure determination efforts, knowledge of the structural features of riboswitches in their ligand-free forms and their ligand-response mechanisms giving rise to regulatory control is lacking. Here we explore the ligand-induced folding process of the S-adenosylmethionine type II (SAM-II) riboswitch using chemical and biophysical methods, including NMR and fluorescence spectroscopy, and single-molecule fluorescence imaging. The data reveal that the unliganded SAM-II riboswitch is dynamic in nature, in that its stem-loop element becomes engaged in a pseudoknot fold through base-pairing with nucleosides in the 3' overhang containing the Shine-Dalgarno sequence. Although the pseudoknot structure is highly transient in the absence of its ligand, S-adenosylmethionine (SAM), it becomes conformationally restrained upon ligand recognition, through a conformational capture mechanism. These insights provide a molecular understanding of riboswitch dynamics that shed new light on the mechanism of riboswitch-mediated translational regulation.     

Purine sensing by riboswitches.

Structured mRNA elements called riboswitches control gene expression by binding to small metabolites. Over a dozen riboswitch classes have been characterized that target a broad range of molecules and vary widely in size and secondary structure. Four of the known riboswitch classes recognize purines or modified purines. Three of these classes are closely related in conserved sequence and secondary structure, but members of these classes selectively recognize guanine, adenine or 2'-deoxyguanosine. Members of the fourth riboswitch class adopt a distinct structure to form a selective binding pocket for the guanine analogue preQ(1) (7-aminomethyl-7-deazaguanine). All four classes of purine-sensing riboswitches are most likely to recognize their respective metabolites by utilizing a riboswitch residue to make a canonical Watson-Crick base-pair with the ligand. This review will provide a summary of the purine-sensing riboswitches, as well as discuss the complex functions and applications of these RNAs.     

Purine sensing by riboswitches

Structured mRNA elements called riboswitches control gene expression by binding to small metabolites. Over a dozen riboswitch classes have been characterized that target a broad range of molecules and vary widely in size and secondary structure. Four of the known riboswitch classes recognize purines or modified purines. Three of these classes are closely related in conserved sequence and secondary structure, but members of these classes selectively recognize guanine, adenine or 2'-deoxyguanosine. Members of the fourth riboswitch class adopt a distinct structure to form a selective binding pocket for the guanine analogue preQ(1) (7-aminomethyl-7-deazaguanine). All four classes of purine-sensing riboswitches are most likely to recognize their respective metabolites by utilizing a riboswitch residue to make a canonical Watson-Crick base-pair with the ligand. This review will provide a summary of the purine-sensing riboswitches, as well as discuss the complex functions and applications of these RNAs.