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.     

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.     

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.     

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.     

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.     

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.     

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.     

Ligand-induced stabilization of the aptamer terminal helix in the add adenine riboswitch

Riboswitches are structured mRNA elements that modulate gene expression. They undergo conformational changes triggered by highly specific interactions with sensed metabolites. Among the structural rearrangements engaged by riboswitches, the forming and melting of the aptamer terminal helix, the so-called P1 stem, is essential for genetic control. The structural mechanisms by which this conformational change is modulated upon ligand binding mostly remain to be elucidated. Here, we used pulling molecular dynamics simulations to study the thermodynamics of the P1 stem in the add adenine riboswitch. The P1 ligand-dependent stabilization was quantified in terms of free energy and compared with thermodynamic data. This comparison suggests a model for the aptamer folding in which direct P1-ligand interactions play a minor role on the conformational switch when compared with those related to the ligand-induced aptamer preorganization.     

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

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

Ligand-induced conformational capture of a synthetic tetracycline riboswitch revealed by pulse EPR

RNA aptamers are in vitro-selected binding domains that recognize their respective ligand with high affinity and specificity. They are characterized by complex three-dimensional conformations providing preformed binding pockets that undergo conformational changes upon ligand binding. Small molecule-binding aptamers have been exploited as synthetic riboswitches for conditional gene expression in various organisms. In the present study, double electron-electron resonance (DEER) spectroscopy combined with site-directed spin labeling was used to elucidate the conformational transition of a tetracycline aptamer upon ligand binding. Different sites were selected for post-synthetic introduction of either the (1-oxyl-2,2,5,5-tetramethylpyrroline-3-methyl) methanethiosulfonate by reaction with a 4-thiouridine modified RNA or of 4-isocyanato-2,6-tetramethylpiperidyl-N-oxid spin label by reaction with 2'-aminouridine modified RNA. The results of the DEER experiments indicate the presence of a thermodynamic equilibrium between two aptamer conformations in the free state and capture of one conformation upon tetracycline binding.     

Role of ligand binding in structural organization of add A-riboswitch aptamer: a molecular dynamics simulation

The specific binding of ligands is the first step of gene expression or translation regulation by riboswitches. However, understanding the mechanism of the specific binding is still difficult because the tertiary structures of the riboswitch aptamers are available almost only for ligand-bound state at present. In this paper we hope to give some insights into this problem through the studies of the role of ligand-aptamer interaction in the structural organization of add A-riboswitch aptamer, based on the crystal structure of the ligand-bound aptamer. We use all-atom molecular dynamics to simulate the behaviors of the aptamer in ligand-bound, free and mutated states by Amber force field. The results show that the correct paring of the ligand adenine with the nucleotide U74 in the binding pocket is crucial to stabilizing the conformations of the ligand-bound aptamer, especially the helix P1 connecting the expression platform. Our results also suggest that both the nucleotide U74 and U51 may be the key sites of the ligand recognition but the former has much higher probability as the initial docking site. This is in agreement with previous experimental results.     

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.     

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.     

Binding of 30S Ribosome Induces Single-stranded Conformation Within and Downstream of the Expression Platform in a Translational Riboswitch

Translational riboswitches are bacterial gene regulatory elements found in the 5′-untranslated region of mRNAs. They operate through a conformational refolding reaction that is triggered by a concentration change of a modulating small molecular ligand. The translation initiation region (TIR) is either released from or incorporated into base pairing interactions through the conformational switch. Hence, initiation of translation is regulated by the accessibility of the Shine-Dalgarno sequence and start codon. Interaction with the 30S ribosome is indispensable for the structural switch between functional OFF and ON states. However, on a molecular level it is still not fully resolved how the ribosome is accommodated near or at the translation initiation region in the context of translational riboswitches. The standby model of translation initiation postulates a binding site where the mRNA enters the ribosome and where it resides until the initiation site becomes unstructured and accessible. We here investigated the adenine-sensing riboswitch from Vibrio vulnificus. By application of a ¹⁹F labelling strategy for NMR spectroscopy that utilizes ligation techniques to synthesize differentially ¹⁹F labelled riboswitch molecules we show that nucleotides directly downstream of the riboswitch domain are first involved in productive interaction with the 30S ribosomal subunit. Upon the concerted action of ligand and the ribosomal protein rS1 the TIR becomes available and subsequently the 30S ribosome can slide towards the TIR. It will be interesting to see whether this is a general feature in translational riboswitches or if riboswitches exist where this region is structured and represent yet another layer of regulation.