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.     

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.     

Dissecting the influence of Mg2+ on 3D architecture and ligand-binding of the guanine-sensing riboswitch aptamer domain

Long-range tertiary interactions determine the three-dimensional structure of a number of metabolite-binding riboswitch RNA elements and were found to be important for their regulatory function. For the guanine-sensing riboswitch of the Bacillus subtilis xpt-pbuX operon, our previous NMR-spectroscopic studies indicated pre-formation of long-range tertiary contacts in the ligand-free state of its aptamer domain. Loss of the structural pre-organization in a mutant of this RNA (G37A/C61U) resulted in the requirement of Mg²⁺ for ligand binding. Here, we investigate structural and stability aspects of the wild-type aptamer domain (Gsw) and the G37A/C61U-mutant (Gsw^loop) of the guanine-sensing riboswitch and their Mg²⁺-induced folding characteristics to dissect the role of long-range tertiary interactions, the link between pre-formation of structural elements and ligand-binding properties and the functional stability. Destabilization of the long-range interactions as a result of the introduced mutations for Gsw^loop or the increase in temperature for both Gsw and Gsw^loop involves pronounced alterations of the conformational ensemble characteristics of the ligand-free state of the riboswitch. The increased flexibility of the conformational ensemble can, however, be compensated by Mg²⁺. We propose that reduction of conformational dynamics in remote regions of the riboswitch aptamer domain is the minimal pre-requisite to pre-organize the core region for specific ligand binding.     

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.     

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.     

Conformational plasticity and dynamic interactions of the N-terminal domain of the chemokine receptor CXCR1

The dynamic interactions between G protein-coupled receptors (GPCRs) and their cognate protein partners are central to several cell signaling pathways. For example, the association of CXC chemokine receptor 1 (CXCR1) with its cognate chemokine, interleukin-8 (IL8 or CXCL8) initiates pathways leading to neutrophil-mediated immune responses. The N-terminal domain of chemokine receptors confers ligand selectivity, but unfortunately the conformational dynamics of this intrinsically disordered region remains unresolved. In this work, we have explored the interaction of CXCR1 with IL8 by microsecond time scale coarse-grain simulations, complemented by atomistic models and NMR chemical shift predictions. We show that the conformational plasticity of the apo-receptor N-terminal domain is restricted upon ligand binding, driving it to an open C-shaped conformation. Importantly, we corroborated the dynamic complex sampled in our simulations against chemical shift perturbations reported by previous NMR studies and show that the trends are similar. Our results indicate that chemical shift perturbation is often not a reporter of residue contacts in such dynamic associations. We believe our results represent a step forward in devising a strategy to understand intrinsically disordered regions in GPCRs and how they acquire functionally important conformational ensembles in dynamic protein-protein interfaces.     

Conformational dynamics of the tetracycline-binding aptamer

The conformational dynamics induced by ligand binding to the tetracycline-binding aptamer is monitored via stopped-flow fluorescence spectroscopy and time-correlated single photon counting experiments. The fluorescence of the ligand is sensitive to changes within the tertiary structure of the aptamer during and after the binding process. In addition to the wild-type aptamer, the mutants A9G, A13U and A50U are examined, where bases important for regulation are changed to inhibit the aptamer's function. Our results suggest a very fast two-step-mechanism for the binding of the ligand to the aptamer that can be interpreted as a binding step followed by a reorganization of the aptamer to accommodate the ligand. Binding to the two direct contact points A13 and A50 was found to occur in the first binding step. The exchange of the structurally important base A9 for guanine induces an enormous deceleration of the overall binding process, which is mainly rooted in an enhancement of the back reaction of the first binding step by several orders of magnitude. This indicates a significant loss of tertiary structure of the aptamer in the absence of the base A9, and underlines the importance of pre-organization on the overall binding process of the tetracycline-binding aptamer.     

Conformational Destabilization of Immunoglobulin G Increases the Low pH Binding Affinity with the Neonatal Fc Receptor

Crystallographic evidence suggests that the pH-dependent affinity of IgG molecules for the neonatal Fc receptor (FcRn) receptor primarily arises from salt bridges involving IgG histidine residues, resulting in moderate affinity at mildly acidic conditions. However, this view does not explain the diversity in affinity found in IgG variants, such as the YTE mutant (M252Y,S254T,T256E), which increases affinity to FcRn by up to 10×. Here we compare hydrogen exchange measurements at pH 7.0 and pH 5.5 with and without FcRn bound with surface plasmon resonance estimates of dissociation constants and FcRn affinity chromatography. The combination of experimental results demonstrates that differences between an IgG and its cognate YTE mutant vary with their pH-sensitive dynamics prior to binding FcRn. The conformational dynamics of these two molecules are nearly indistinguishable upon binding FcRn. We present evidence that pH-induced destabilization in the CH2/3 domain interface of IgG increases binding affinity by breaking intramolecular H-bonds and increases side-chain adaptability in sites that form intermolecular contacts with FcRn. Our results provide new insights into the mechanism of pH-dependent affinity in IgG-FcRn interactions and exemplify the important and often ignored role of intrinsic conformational dynamics in a protein ligand, to dictate affinity for biologically important receptors.     

Folding of the SAM-I riboswitch: A tale with a twist

Riboswitches are ligand-dependent RNA genetic regulators that control gene expression by altering their structures. The elucidation of riboswitch conformational changes before and after ligand recognition is crucial to understand how riboswitches can achieve high ligand binding affinity and discrimination against cellular analogs. The detailed characterization of riboswitch folding pathways suggest that they may use their intrinsic conformational dynamics to sample a large array of structures, some of which being nearly identical to ligand-bound molecules. Some of these structural conformers can be "captured" upon ligand binding, which is crucial for the outcome of gene regulation. Recent studies about the SAM-I riboswitch have revealed unexpected and previously unknown RNA folding mechanisms. For instance, the observed helical twist of the P1 stem upon ligand binding to the SAM-I aptamer adds a new element in the repertoire of RNA strategies for recognition of small metabolites. From an RNA folding perspective, these findings also strongly indicate that the SAM-I riboswitch could achieve ligand recognition by using an optimized combination of conformational capture and induced-fit approaches, a feature that may be shared by other RNA regulatory sequences.     

A Structural Basis for the Recognition of 2′-Deoxyguanosine by the Purine Riboswitch

Riboswitches are noncoding RNA elements that are commonly found in the 5′-untranslated region of bacterial mRNA. Binding of a small-molecule metabolite to the riboswitch aptamer domain guides the folding of the downstream sequence into one of two mutually exclusive secondary structures that directs gene expression. The purine riboswitch family, which regulates aspects of purine biosynthesis and transport, contains three distinct classes that specifically recognize guanine/hypoxanthine, adenine, or 2′-deoxyguanosine (dG). Structural analysis of the guanine and adenine classes revealed a binding pocket that almost completely buries the nucleobase within the core of the folded RNA. Thus, it is somewhat surprising that this family of RNA elements also recognizes dG. We have used a combination of structural and biochemical techniques to understand how the guanine riboswitch could be converted into a dG binder and the structural basis for dG recognition. These studies reveal that a limited number of sequence changes to a guanine-sensing RNA are required to cause a specificity switch from guanine to 2′-deoxyguanosine, and to impart an altered structure for accommodating the additional deoxyribose sugar moiety.     

Role of a heterogeneous free state in the formation of a specific RNA-theophylline complex

The helical regions of RNA are generally very stable, but the single-stranded and loop regions often exist as an ensemble of conformations in solution. The theophylline-binding RNA aptamer forms a very stable structure when bound to the bronchodilator theophylline, but the theophylline binding site is not stably formed in the absence of ligand. The kinetics for theophylline binding were measured here by stopped-flow fluorescence spectroscopy to probe the mechanism for theophylline binding in this RNA aptamer. The kinetic studies showed that formation of the RNA-theophylline complex is over 1000 times slower than a diffusion-controlled rate, and the high affinity of the RNA-theophylline complex arises primarily from a slow dissociation rate for the complex. A theophylline-independent rate was observed for formation of the theophylline-RNA complex at high theophylline concentration, indicating that a conformational change in the RNA is the rate-limiting step in complex formation under these conditions. The RNA-theophylline complex requires divalent metal ions, such as Mg2+, to form a high-affinity complex, and there is a greater than 10000-fold reduction in affinity for theophylline in the absence of Mg2+. This decrease in binding affinity in the absence of Mg2+ results primarily from an increased dissociation rate for the complex. The implications of an ensemble of conformations in the free state of this theophylline-binding RNA are discussed and compared with mechanisms for formation of protein-ligand complexes.     

Structural distinctions between NAD+ riboswitch domains 1 and 2 determine differential folding and ligand binding

Riboswitches are important gene regulatory elements frequently encountered in bacterial mRNAs. The recently discovered nadA riboswitch contains two similar, tandemly arrayed aptamer domains, with the first domain possessing high affinity for nicotinamide adenine dinucleotide (NAD⁺). The second domain which comprises the ribosomal binding site in a putative regulatory helix, however, has withdrawn from detection of ligand-induced structural modulation thus far, and therefore, the identity of the cognate ligand and the regulation mechanism have remained unclear. Here, we report crystal structures of both riboswitch domains, each bound to NAD⁺. Furthermore, we demonstrate that ligand binding to domain 2 requires significantly higher concentrations of NAD⁺ (or ADP retaining analogs) compared to domain 1. Using a fluorescence spectroscopic approach, we further shed light on the structural features which are responsible for the different ligand affinities, and describe the Mg²⁺-dependent, distinct folding and pre-organization of their binding pockets. Finally, we speculate about possible scenarios for nadA RNA gene regulation as a putative two-concentration sensor module for a time-controlled signal that is primed and stalled by the gene regulation machinery at low ligand concentrations (domain 1), and finally triggers repression of translation as soon as high ligand concentrations are reached in the cell (domain 2).     

RNA aptamers that bind L-arginine with sub-micromolar dissociation constants and high enantioselectivity.

A completely randomized RNA pool as well as a degenerate pool comprised of an RNA sequence which binds citrulline with a dissociation constant of 0 muM were used to select for tight binding arginine specific RNA aptamers. A modified in vitro selection scheme, based on affinity chromatography was applied to allow the enrichment of high affinity solution binders. The selection scheme included a negative selection with the non-cognate ligand citrulline, and a heat denaturation step prior to affinity elution with an excess of the cognate ligand arginine. After 20 cycles the majority of the pools bound specifically to the arginine matrix even after denaturation/renaturation in the presence of 20 mM of a non-cognate amino acid. When denatured and eluted in the presence of 20 mM arginine, the selected RNAs quantitatively washed off the column. These RNA aptamers were cloned and sequenced. Equilibrium dialysis performed with the most abundant clone among the selected sequence revealed Kd values of 330 nM for the RNA/arginine affinity, which is nearly a 200-fold improvement over the tightest binding arginine binding RNAs known to date. Arginine recognition by this RNA is highly enantioselectice: L- arginine is bound 12 000-fold better than D-arginine. Chemical modification analysis revealed that the secondary structure of the aptamer might contain a pseudoknot motif. Our tight binding arginine aptamers join a number of natural and in vitro selected RNAs which recognize arginine. The RNAs described here compare in their binding affinity with the tightest binding RNA aptamers for low molecular weight molecules isolated in other in vitro selection experiments.     

Molecular dynamics simulation on the allosteric analysis of the c‐di‐GMP class I riboswitch induced by ligand binding

Riboswitches are RNA molecules that regulate gene expression using conformation change, affected by binding of small molecule ligands. Although a number of ligand‐bound aptamer complex structures have been solved, it is important to know ligand‐free conformations of the aptamers in order to understand the mechanism of specific binding by ligands. In this paper, we use dynamics simulations on a series of models to characterize the ligand‐free and ligand‐bound aptamer domain of the c‐di‐GMP class I (GEMM‐I) riboswitch. The results revealed that the ligand‐free aptamer has a stable state with a folded P2 and P3 helix, an unfolded P1 helix and open binding pocket. The first Mg ions binding to the aptamer is structurally favorable for the successive c‐di‐GMP binding. The P1 helix forms when c‐di‐GMP is successive bound. Three key junctions J1/2, J2/3 and J1/3 in the GEMM‐I riboswitch contributing to the formation of P1 helix have been found. The binding of the c‐di‐GMP ligand to the GEMM‐I riboswitch induces the riboswitch's regulation through the direct allosteric communication network in GEMM‐I riboswitch from the c‐di‐GMP binding sites in the J1/2 and J1/3 junctions to the P1 helix, the indirect ones from those in the J2/3 and P2 communicating to P1 helix via the J1/2 and J1/3 media.     

Mutual Protein-Ligand Conformational Selection Drives cGMP vs. cAMP Selectivity in Protein Kinase G

Protein kinase G (PKG) is a major receptor of cGMP, and controls signaling pathways distinct from those regulated by cAMP. However, the contributions of the two substituents that differentiate cGMP from cAMP (i.e. 6-oxo and 2-NH₂) to the cGMP-versus-cAMP selectivity of PKG remain unclear. Here, using NMR to map how binding affinity and dynamics of the protein and ligand vary along a ligand double-substitution cycle, we show that the contributions of the two substituents to binding affinity are surprisingly non-additive. Such non-additivity stems primarily from mutual protein–ligand conformational selection, whereby not only does the ligand select for a preferred protein conformation upon binding, but also, the protein selects for a preferred ligand conformation. The 6-oxo substituent mainly controls the conformational equilibrium of the bound protein, while the 2-NH₂ substituent primarily controls the conformational equilibrium of the unbound ligand (i.e. syn versus anti). Therefore, understanding the conformational dynamics of both the protein and ligand is essential to explain the cGMP-versus-cAMP selectivity of PKG.     

Solution structure of an informationally complex high-affinity RNA aptamer to GTP

Higher-affinity RNA aptamers to GTP are more informationally complex than lower-affinity aptamers. Analog binding studies have shown that the additional information needed to improve affinity does not specify more interactions with the ligand. In light of those observations, we would like to understand the structural characteristics that enable complex aptamers to bind their ligands with higher affinity. Here we present the solution structure of the 41-nt Class I GTP aptamer (K(d) = 75 nM) as determined by NMR. The backbone of the aptamer forms a reverse-S that shapes the binding pocket. The ligand nucleobase stacks between purine platforms and makes hydrogen bonds with the edge of another base. Interestingly, the local modes of interaction for the Class I aptamer and an RNA aptamer that binds ATP with a K(d) of 6 microM are very much alike. The aptamers exhibit nearly identical levels of binding specificity and fraction of ligand sequestered from the solvent (81%-85%). However, the GTP aptamer is more informationally complex (approximately 45 vs. 35 bits) and has a larger recognition bulge (15 vs. 12 nucleotides) with many more stabilizing base-base interactions. Because the aptamers have similar modes of ligand binding, we conclude that the stabilizing structural elements in the Class I aptamer are responsible for much of the difference in K(d). These results are consistent with the hypothesis that increasing the number of intra-RNA interactions, rather than adding specific contacts to the ligand, is the simplest way to improve binding affinity.     

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.     

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.