Thermodynamic analysis of ligand binding and ligand binding-induced tertiary structure formation by the thiamine pyrophosphate riboswitch

The thi-box riboswitch regulates gene expression in response to the intracellular concentration of thiamine pyrophosphate (TPP) in archaea, bacteria, and eukarya. To complement previous biochemical, genetic, and structural studies of this phylogenetically widespread RNA domain, we have characterized its interaction with TPP by isothermal titration calorimetry. This shows that TPP binding is highly dependent on Mg²⁺ concentration. The dissociation constant decreases from ∼200 nM at 0.5 mM Mg²⁺ concentration to ∼9 nM at 2.5 mM Mg²⁺ concentration. Binding is enthalpically driven, but the unfavorable entropy of binding decreases as Mg²⁺ concentration rises, suggesting that divalent cations serve to pre-organize the RNA. Mutagenesis, biochemical analysis, and a new crystal structure of the riboswitch suggest that a critical element that participates in organizing the riboswitch structure is the tertiary interaction formed between the P3 and L5 regions. This tertiary contact is distant from the TPP binding site, but calorimetric analysis reveals that even subtle mutations in L5 can have readily detectable effects on TPP binding. The thermodynamic signatures of these mutations, namely decreased favorable enthalpy of binding and small effects on entropy of binding, are consistent with the P3–L5 association contributing allosterically to TPP-induced compaction of the RNA.     

Salt bridge exchange binding mechanism between streptavidin and its DNA aptamer – thermodynamics and spectroscopic evidences

Protein‐nucleic acids binding driven by electrostatic interactions typically are characterized by the release of counter ions, and the salt‐inhibited binding association constant (K_a) and the magnitude of exothermic binding enthalpy (ΔH). Here, we report a non‐classical thermodynamics of streptavidin (SA)–aptamer binding in NaCl (140–350 mM) solutions near room temperatures (23–27 °C). By using isothermal titration calorimetry (ITC) and circular dichroism (CD)/fluorescence spectroscopy, we found that the binding was enthalpy driven with a large entropy cost (ΔH ?20.58 kcal mol^?1, TΔS ?10.99 kcal mol^?1, and K_a 1.08 × 10⁷ M^?1 at 140 mM NaCl 25 °C). With the raise of salt concentrations, the ΔH became more exothermic, yet the K_a was almost unchanged (ΔH ?26.29 kcal mol^?1 and K_a 1.50 × 10⁷ M^?1 at 350 mM NaCl 25 °C). The data suggest that no counter Na⁺ was released in the binding. Spectroscopy data suggest that the binding, with a stoichiometry of 2, was accompanied with substantial conformational changes on SA, and the changes were insensitive to the variation of salt concentrations. To account for the non‐classical results, we propose a salt bridge exchange model. The intramolecular binding‐site salt bridge(s) of the free SA and the charged phosphate group of aptamers re‐organize to form the binding complex by forming a new intermolecular salt bridge(s). The salt bridge exchange binding process requires minimum amount of counter ions releasing but dehydration of the contacting surface of SA and the aptamer. The energy required for dehydration is reduced in the case of binding solution with higher salt concentration and account for the higher binding exothermic mainly. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

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.     

Predicting absolute ligand binding free energies to a simple model site

A central challenge in structure-based ligand design is the accurate prediction of binding free energies. Here we apply alchemical free energy calculations in explicit solvent to predict ligand binding in a model cavity in T4 lysozyme. Even in this simple site, there are challenges. We made systematic improvements, beginning with single poses from docking, then including multiple poses, additional protein conformational changes, and using an improved charge model. Computed absolute binding free energies had an RMS error of 1.9 kcal/mol relative to previously determined experimental values. In blind prospective tests, the methods correctly discriminated between several true ligands and decoys in a set of putative binders identified by docking. In these prospective tests, the RMS error in predicted binding free energies relative to those subsequently determined experimentally was only 0.6 kcal/mol. X-ray crystal structures of the new ligands bound in the cavity corresponded closely to predictions from the free energy calculations, but sometimes differed from those predicted by docking. Finally, we examined the impact of holding the protein rigid, as in docking, with a view to learning how approximations made in docking affect accuracy and how they may be improved.     

Evaluation of Staphylococcus aureus DNA aptamer by enzyme‐linked aptamer assay and isothermal titration calorimetry

To monitor the specificity of Staphylococcus aureus aptamer (SA‐31) against its target cell, we used enzyme‐linked aptamer assay. In the presence of target cell, horseradish peroxidase–conjugated streptavidin bound to biotin‐labeled SA‐31 showed specific binding to S   aureus among 3 different bacteria with limit of detection of 10³ colony‐forming unit per milliliter. The apparent K _a was 1.39 μM⁻¹ ± 0.3 μM⁻¹. The binding of SA‐31 to membrane proteins extracted from cell surface was characterized using isothermal titration calorimetry, and the effect of changes in binding temperature and salt concentrations of binding buffer was evaluated based on thermodynamic parameters (K _a, ΔH , and ΔG ). Since binding of aptamer to its targets solely depends on its 3‐dimensional structure under experimental conditions used in selection process, the change in temperature and ion concentration changed the affinity of SA‐31 to its target on surface of bacteria. At 4°C, SA‐31 did not show an affinity to its target with poor heat change upon injection of membrane fraction to aptamer solution. However, the apparent association constants of SA‐31 slightly varied from K _a = 1.56 μM⁻¹ ± 0.69 μM⁻¹ at 25°C to K _a = 1.03 μM⁻¹ ± 0.9 μM⁻¹ at 37°C. At spontaneously occurring exothermic binding reactions, affinities of S  aureus aptamer to its target were also 9.44 μM⁻¹ ± 0.38 μM⁻¹ at 50mM, 1.60 μM⁻¹ ± 0.11 μM⁻¹ at 137mM, and 3.28 μM⁻¹ ± 0.46 μM⁻¹ at 200 mM of salt concentration. In this study, it was demonstrated that enzyme‐linked aptamer assay and isothermal titration calorimetry were useful tools for studying the fundamental binding mechanism between a DNA aptamer and its target on the outer surface of S  aureus .     

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.     

Ligand binding by the tandem glycine riboswitch depends on aptamer dimerization but not double ligand occupancy

The glycine riboswitch predominantly exists as a tandem structure, with two adjacent, homologous ligand-binding domains (aptamers), followed by a single expression platform. The recent identification of a leader helix, the inclusion of which eliminates cooperativity between the aptamers, has reopened the debate over the purpose of the tandem structure of the glycine riboswitch. An equilibrium dialysis-based assay was combined with binding-site mutations to monitor glycine binding in each ligand-binding site independently to understand the role of each aptamer in glycine binding and riboswitch tertiary interactions. A series of mutations disrupting the dimer interface was used to probe how dimerization impacts ligand binding by the tandem glycine riboswitch. While the wild-type tandem riboswitch binds two glycine equivalents, one for each aptamer, both individual aptamers are capable of binding glycine when the other aptamer is unoccupied. Intriguingly, glycine binding by aptamer-1 is more sensitive to dimerization than glycine binding by aptamer-2 in the context of the tandem riboswitch. However, monomeric aptamer-2 shows dramatically weakened glycine-binding affinity. In addition, dimerization of the two aptamers in trans is dependent on glycine binding in at least one aptamer. We propose a revised model for tandem riboswitch function that is consistent with these results, wherein ligand binding in aptamer-1 is linked to aptamer dimerization and stabilizes the P1 stem of aptamer-2, which controls the expression platform.     

Thermodynamic characterization of binding Oxytricha nova single strand telomere DNA with the alpha protein N-terminal domain

The Oxytricha nova telemere binding protein alpha subunit binds single strand DNA and participates in a nucleoprotein complex that protects the very ends of chromosomes. To understand how the N-terminal, DNA binding domain of alpha interacts with DNA we measured the stoichiometry, enthalpy (DeltaH), entropy (DeltaS), and dissociation constant (K(D-DNA)) for binding telomere DNA fragments at different temperatures and salt concentrations using native gel electrophoresis and isothermal titration calorimetry (ITC). About 85% of the total free energy of binding corresponded with non-electrostatic interactions for all DNAs. Telomere DNA fragments d(T(2)G(4)), d(T(4)G(4)), d(G(3)T(4)G(4)), and d(G(4)T(4)G(4)) each formed monovalent protein complexes. In the case of d(T(4)G(4)T(4)G(4)), which has two tandemly repeated d(TTTTTGGGG) telomere motifs, two binding sites were observed. The high-affinity "A site" has a dissociation constant, K(D-DNA(A)) = 13(+/-4) nM, while the low-affinity "B site" is characterized by K(D-DNA(B)) = 5600(+/-600) nM at 25 degrees C. Nucleotide substitution variants verified that the A site corresponds principally with the 3'-terminal portion of d(T(4)G(4)T(4)G(4)). The relative contributions of entropy (DeltaS) and enthalpy (DeltaH) for binding reactions were DNA length-dependent as was heat capacity (DeltaCp). These trends with respect to DNA length likely reflect structural transitions in the DNA molecule that are coupled with DNA-protein association. Results presented here are important for understanding early intermediates and subsequent stages in the assembly of the full telomere nucleoprotein complex and how binding events can prepare the telomere DNA for extension by telomerase, a critical event in telomere biology.     

Calorimetric investigation of phosphorylated and non-phosphorylated peptide ligand binding to the human Grb7-SH2 domain

Grb7 is a member of the Grb7 family of proteins, which also includes Grb10 and Grb14. All three proteins have been found to be overexpressed in certain cancers and cancer cell lines. In particular, Grb7 (along with the receptor tyrosine kinase erbB2) is overexpressed in 20-30% of breast cancers. In general, growth factor receptor bound (Grb) proteins bind to activated membrane-bound receptor tyrosine kinases (RTKs; e.g., the epidermal growth factor receptor, EGFR) through their Src homology 2 (SH2) domains. In particular, Grb7 binds to erbB2 (a.k.a. EGFR2) and may be involved in cell signaling pathways that promote the formation of metastases and inflammatory responses. In previous studies, we reported the solution structure and the backbone relaxation behavior of the Grb7-SH2/erbB2 peptide complex. In this study, isothermal titration calorimetry studies have been completed by measuring the thermodynamic binding parameters of several phosphorylated and non-phosphorylated peptides representative of natural Grb7 receptor ligands as well as ligands developed through combinatorial peptide screening methods. The entirety of these calorimetric studies is interpreted in an effort to describe the specific ligand binding characteristics of the Grb7 protein.     

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.     

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.     

Modulating the folding stability and ligand binding affinity of Pin1 WW domain by proline ring puckering

The pyrrolidine side chain makes proline play a unique role in protein structure and function. The C^γ ring pucker preference and the cis– trans peptidyl bond ratio can be mediated via stereoelectronic effects. Here we used a compact triple‐stranded antiparallel β‐sheet protein, the human Pin1 WW domain, to study the consequences of implanting a preorganized C^γ ring pucker on protein structure and function. The conserved Pro37 is a key residue involved in one hydrophobic core, plays an important role in the WW domain, and adopts a C^γ‐endo ring pucker in the native structure. Pro37 was replaced with C^γ‐exo biased pucker derivatives: (2S,4R)‐4‐hydroxyproline (4R‐Hyp), (2S,4R)‐4‐fluoroproline (4R‐Flp), (2S,4R)‐4‐methoxyproline (4R‐Mop), and C^γ‐endo biased pucker derivatives: (2S,4S)‐4‐hydroxyproline (4S‐hyp), (2S,4S)‐4‐fluoroproline (4S‐flp), (2S,4S)‐4‐methoxyproline (4S‐mop) to examine how a preorganized pucker affects the folding stability and ligand‐binding affinity. Circular dichroism measurements indicate that among the variants, only the one with 4S‐flp substitution (P37flp) is more stable than the wild type, suggesting that the stabilization effects originated from preorganization of the backbone conformation and the hydrophobicity of C? F group. Analysis of ligand‐binding affinity using isothermal titration calorimetry revealed that only P37flp has a stronger ligand affinity than the wild type, showing that 4S‐flp can stabilize the WW domain and increase its ligand affinity. Together we have used 4‐substituted proline derivatives and the WW domain to demonstrate that proline ring puckering can be a key factor in determining the folding stability of a protein but the choice of the derivative groups is also critical. Proteins 2014; 82:67–76. © 2013 Wiley Periodicals, Inc.     

Thermodynamic characterization of OsGID1-gibberellin binding using calorimetry and docking simulations

Gibberellins (GAs) are phytohormones regulating various developmental processes in plants. In rice, the initial GA-signaling events involve the binding of a GA to the soluble GA receptor protein, GID1. Although X-ray structures for certain GID1/GA complexes have recently been determined, an examination of the complexes does not fully clarify how GID1s discriminate among different GAs. Herein, we present a study aimed at defining the types of forces important to binding via a combination of isothermal titration calorimetry (ITC) and computational docking studies that employed rice GID1 (OsGID1), OsGID1 mutants, which were designed to have a decreased possible number of hydrogen bonds with bound GA, and GA variants. We find that, in general, GA binding is enthalpically driven and that a hydrogen bond between the phenolic hydroxyl of OsGID1 Tyr134 and the C-3 hydroxyl of a GA is a defining structural element. A hydrogen-bond network that involves the C-6 carboxyl of a GA that directly hydrogen bonds the hydroxyl of Ser198 and indirectly, via a two-water-molecule network, the phenolic hydroxyl of Tyr329 and the NH of the amide side-chain of Asn255 is also important for GA binding. The binding of OsGID1 by GA(1) is the most enthalpically driven association found for the biologically active GAs evaluated in this study. This observation might be a consequence of a hydrogen bond formed between the hydroxyl at the C-13 position of GA(1) and the main chain carbonyl of OsGID1 Phe245. Our results demonstrate that by combining ITC experiments and computational methods much can be learned about the thermodynamics of ligand/protein binding.     

Identification of a tertiary interaction important for cooperative ligand binding by the glycine riboswitch

The glycine riboswitch has a tandem dual aptamer configuration, where each aptamer is a separate ligand-binding domain, but the aptamers function together to bind glycine cooperatively. We sought to understand the molecular basis of glycine riboswitch cooperativity by comparing sites of tertiary contacts in a series of cooperative and noncooperative glycine riboswitch mutants using hydroxyl radical footprinting, in-line probing, and native gel-shift studies. The results illustrate the importance of a direct or indirect interaction between the P3b hairpin of aptamer 2 and the P1 helix of aptamer 1 in cooperative glycine binding. Furthermore, our data support a model in which glycine binding is sequential; where the binding of glycine to the second aptamer allows tertiary interactions to be made that facilitate binding of a second glycine molecule to the first aptamer. These results provide insight into cooperative ligand binding in RNA macromolecules.     

Kinetics and energetics of ligand binding determined by microcalorimetry: insights into active site mobility in a psychrophilic alpha-amylase

A new microcalorimetric method for recording the kinetic parameters k(cat), K(m) and K(i) of alpha-amylases using polysaccharides and oligosaccharides as substrates is described. This method is based on the heat released by glycosidic bond hydrolysis. The method has been developed to study the active site properties of the cold-active alpha-amylase produced by an Antarctic psychrophilic bacterium in comparison with its closest structural homolog from pig pancreas. It is shown that the psychrophilic alpha-amylase is more active on large macromolecular substrates and that the higher rate constants k(cat) are gained at the expense of a lower affinity for the substrate. The active site is able to accommodate larger inhibitory complexes, resulting in a mixed-type inhibition of starch hydrolysis by maltose. A method for recording the binding enthalpies by isothermal titration calorimetry in a low-affinity system has been developed, allowing analysis of the energetics of weak ligand binding using the allosteric activator chloride. It is shown that the low affinity of the psychrophilic alpha-amylase for chloride is entropically driven. The high enthalpic and entropic contributions of activator binding suggest large structural fluctuations between the free and the bound states of the cold-active enzyme. The kinetic and thermodynamic data for the psychrophilic alpha-amylase indicate that the strictly conserved side-chains involved in substrate binding and catalysis possess an improved mobility, responsible for activity in the cold, and resulting from the disappearance of stabilizing interactions far from the active site.     

Idiosyncratically tuned switching behavior of riboswitch aptamer domains revealed by comparative small-angle X-ray scattering analysis

Riboswitches are structured mRNA elements that regulate gene expression upon binding specific cellular metabolites. It is thought that the highly conserved metabolite-binding domains of riboswitches undergo conformational change upon binding their cognate ligands. To investigate the generality of such a mechanism, we employed small-angle X-ray scattering (SAXS). We probed the nature of the global metabolite-induced response of the metabolite-binding domains of four different riboswitches that bind, respectively, thiamine pyrophosphate (TPP), flavin mononucleotide (FMN), lysine, and S-adenosyl methionine (SAM). We find that each RNA is unique in its global structural response to metabolite. Whereas some RNAs exhibit distinct free and bound conformations, others are globally insensitive to the presence of metabolite. Thus, a global conformational change of the metabolite-binding domain is not a requirement for riboswitch function. It is possible that the range of behaviors observed by SAXS, rather than being a biophysical idiosyncrasy, reflects adaptation of riboswitches to the regulatory requirements of their individual genomic context.     

Crystal structure of Escherichia coli thiamine pyrophosphate-sensing riboswitch in the apo state

1. Download : Download high-res image (195KB)
2. Download : Download full-size image