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.     

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.     

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.     

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.     

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

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

Structure and mechanism of purine-binding riboswitches

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

Modeling the noncovalent interactions at the metabolite binding site in purine riboswitches

We present gas phase quantum chemical studies on the metabolite binding interactions in two important purine riboswitches, the adenine and guanine riboswitches, at the B3LYP/6-31G(d,p) level of theory. In order to gain insights into the strucutral basis of their discriminative abilities of regulating gene expression, the structural properties and binding energies for the gas phase optimized geometries of the metabolite bound binding pocket are analyzed and compared with their respective crystal geometries. Kitaura-Morokuma analysis has been carried out to calculate and decompose the interaction energy into various components. NBO and AIM analysis has been carried out to understand the strength and nature of binding of the individual aptamer bases with their respective purine metabolites. The Y74 base, U in case of adenine riboswitch and C in case of guanine riboswitch constitutes the only differentiating element between the two binding pockets. As expected, with W:W cis G:C74 interaction contributing more than 50% of the total binding energy, the interaction energy for metabolite binding as calculated for guanine (-46.43 Kcal/mol) is nearly double compared to the corresponding value for that of adenine (-24.73 Kcal/mol) in the crystal context. Variations in the optimized geometries for different models and comparison of relative contribution to metabolite binding involving four conserved bases reveal the possible role of U47:U51 W:H trans pair in the conformational transition of the riboswitch from the metabolite free to metabolite bound state. Our results are also indicative of significant contributions from stacking and magnesium ion interactions toward cooperativity effects in metabolite recognition.     

Rational design of SAM analogues targeting SAM-II riboswitch aptamer

Riboswitches are noncoding RNA elements embedded in 5′-untranslated region of many bacterial mRNAs regulating gene expression in response to essential metabolites. They are unique from other RNA targets because they have evolved to form specific structural receptors for the purpose of binding small molecular metabolites suggesting that structure-based rational drug design approach may be used in designing metabolite mimics targeting riboswitches. We have developed a fluorescence binding assay for SAM-II riboswitch aptamer and identified an S-adenosylmethionine (SAM) analogue that selectively binds to SAM-II riboswitch aptamer with comparable binding affinity to its native metabolite using structure-based design approach.     

Structural basis for recognition of S-adenosylhomocysteine by riboswitches

S-adenosyl-(L)-homocysteine (SAH) riboswitches are regulatory elements found in bacterial mRNAs that up-regulate genes involved in the S-adenosyl-(L)-methionine (SAM) regeneration cycle. To understand the structural basis of SAH-dependent regulation by RNA, we have solved the structure of its metabolite-binding domain in complex with SAH. This structure reveals an unusual pseudoknot topology that creates a shallow groove on the surface of the RNA that binds SAH primarily through interactions with the adenine ring and methionine main chain atoms and discriminates against SAM through a steric mechanism. Chemical probing and calorimetric analysis indicate that the unliganded RNA can access bound-like conformations that are significantly stabilized by SAH to direct folding of the downstream regulatory switch. Strikingly, we find that metabolites bearing an adenine ring, including ATP, bind this aptamer with sufficiently high affinity such that normal intracellular concentrations of these compounds may influence regulation of the riboswitch.     

Molecular dynamics simulation study of the binding of purine bases to the aptamer domain of the guanine sensing riboswitch

Riboswitches are a novel class of genetic control elements that function through the direct interaction of small metabolite molecules with structured RNA elements. The ligand is bound with high specificity and affinity to its RNA target and induces conformational changes of the RNA''s secondary and tertiary structure upon binding. To elucidate the molecular basis of the remarkable ligand selectivity and affinity of one of these riboswitches, extensive all-atom molecular dynamics simulations in explicit solvent (≈1 μs total simulation length) of the aptamer domain of the guanine sensing riboswitch are performed. The conformational dynamics is studied when the system is bound to its cognate ligand guanine as well as bound to the non-cognate ligand adenine and in its free form. The simulations indicate that residue U51 in the aptamer domain functions as a general docking platform for purine bases, whereas the interactions between C74 and the ligand are crucial for ligand selectivity. These findings either suggest a two-step ligand recognition process, including a general purine binding step and a subsequent selection of the cognate ligand, or hint at different initial interactions of cognate and noncognate ligands with residues of the ligand binding pocket. To explore possible pathways of complex dissociation, various nonequilibrium simulations are performed which account for the first steps of ligand unbinding. The results delineate the minimal set of conformational changes needed for ligand release, suggest two possible pathways for the dissociation reaction, and underline the importance of long-range tertiary contacts for locking the ligand in the complex.     

BODIPY-modified 2'-deoxyguanosine as a novel tool to detect DNA damages

BODIPY-modified 2′-deoxyguanosine was synthesized for use as a detection reagent for genotoxic compounds. BODIPY-FL is a well known fluorescence reagent whose fluorescent light emission diminishes near a guanine base by a photo-induced electron transfer process. We attached BODIPY-Fl to the 5′ position of the deoxyribose moiety of 2′-deoxyguanosine. Although this compound has low fluorescence activity, when depurination by the action of alkylating reagents and dG oxidation by singlet oxygen occurred, the emission of strong fluorescence was observed. BODIPY-dG was found, therefore, to be a very useful tool for selectively detecting DNA damaging activity particularly in natural environmental extracts.     

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.     

Octameric structure of the human bifunctional enzyme PAICS in purine biosynthesis

Phosphoribosylaminoimidazole carboxylase/phosphoribosylaminoimidazole succinocarboxamide synthetase (PAICS) is an important bifunctional enzyme in de novo purine biosynthesis in vertebrate with both 5-aminoimidazole ribonucleotide carboxylase (AIRc) and 4-(N-succinylcarboxamide)-5-aminoimidazole ribonucleotide synthetase (SAICARs) activities. It becomes an attractive target for rational anticancer drug design, since rapidly dividing cancer cells rely heavily on the purine de novo pathway for synthesis of adenine and guanine, whereas normal cells favor the salvage pathway. Here, we report the crystal structure of human PAICS, the first in the entire PAICS family, at 2.8 Å resolution. It revealed that eight PAICS subunits, each composed of distinct AIRc and SAICARs domains, assemble a compact homo-octamer with an octameric-carboxylase core and four symmetric periphery dimers formed by synthetase domains. Based on structural comparison and functional complementation analyses, the active sites of SAICARs and AIRc were identified, including a putative substrate CO₂-binding site. Furthermore, four symmetry-related, separate tunnel systems in the PAICS octamer were found that connect the active sites of AIRc and SAICARs. This study illustrated the octameric nature of the bifunctional enzyme. Each carboxylase active site is formed by structural elements from three AIRc domains, demonstrating that the octamer structure is essential for the carboxylation activity. Furthermore, the existence of the tunnel system implies a mechanism of intermediate channeling and suggests that the quaternary structure arrangement is crucial for effectively executing the sequential reactions. In addition, this study provides essential structural information for designing PAICS-specific inhibitors for use in cancer chemotherapy.     

The plasticity of a structural motif in RNA: Structural polymorphism of a kink turn as a function of its environment

The k-turn is a widespread structural motif that introduces a tight kink into the helical axis of double-stranded RNA. The adenine bases of consecutive G•A pairs are directed toward the minor groove of the opposing helix, hydrogen bonding in a typical A-minor interaction. We show here that the available structures of k-turns divide into two classes, depending on whether N3 or N1 of the adenine at the 2b position accepts a hydrogen bond from the O2′ at the −1n position. There is a coordinated structural change involving a number of hydrogen bonds between the two classes. We show here that Kt-7 can adopt either the N3 or N1 structures depending on environment. While it has the N1 structure in the ribosome, on engineering it into the SAM-I riboswitch, it changes to the N3 structure, resulting in a significant alteration in the trajectory of the helical arms.     

SAM recognition and conformational switching mechanism in the Bacillus subtilis yitJ S box/SAM-I riboswitch

S-box (SAM-I) riboswitches are a widespread class of riboswitches involved in the regulation of sulfur metabolism in Gram-positive bacteria. We report here the 3.0-Å crystal structure of the aptamer domain of the Bacillus subtilis yitJ S-box (SAM-I) riboswitch bound to S-adenosyl-l-methionine (SAM). The RNA folds into two sets of helical stacks spatially arranged by tertiary interactions including a K-turn and a pseudoknot at a four-way junction. The tertiary structure is further stabilized by metal coordination, extensive ribose zipper interactions, and SAM-mediated tertiary interactions. Despite structural differences in the peripheral regions, the SAM-binding core of the B. subtilis yitJ riboswitch is virtually superimposable with the previously determined Thermoanaerobacter tengcongensis yitJ riboswitch structure, suggesting that a highly conserved ligand-recognition mechanism is utilized by all S-box riboswitches. SHAPE (selective 2′-hydroxyl acylation analyzed by primer extension) chemical probing analysis further revealed that the alternative base-pairing element in the expression platform controls the conformational switching process. In the absence of SAM, the apo yitJ aptamer domain folds predominantly into a pre-binding conformation that resembles, but is not identical with, the SAM-bound state. We propose that SAM enters the ligand-binding site through the “J1/2-J3/4” gate and “locks” down the SAM-bound conformation through an induced-fit mechanism. Temperature-dependent SHAPE revealed that the tertiary interaction-stabilized SAM-binding core is extremely stable, likely due to the cooperative RNA folding behavior. Mutational studies revealed that certain modifications in the SAM-binding region result in loss of SAM binding and constitutive termination, which suggests that these mutations lock the RNA into a form that resembles the SAM-bound form in the absence of SAM.     

Role of the Adenine Ligand on the Stabilization of the Secondary and Tertiary Interactions in the Adenine Riboswitch

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