Evaluation of linked protonation effects in protein binding reactions using isothermal titration calorimetry.

A theoretical development in the evaluation of proton linkage in protein binding reactions by isothermal titration calorimetry (ITC) is presented. For a system in which binding is linked to protonation of an ionizable group on a protein, we show that by performing experiments as a function of pH in buffers with varying ionization enthalpy, one can determine the pK(a)'s of the group responsible for the proton linkage in the free and the liganded states, the protonation enthalpy for this group in these states, as well as the intrinsic energetics for ligand binding (delta H(o), delta S(o), and delta C(p)). Determination of intrinsic energetics in this fashion allows for comparison with energetics calculated empirically from structural information. It is shown that in addition to variation of the ligand binding constant with pH, the observed binding enthalpy and heat capacity change can undergo extreme deviations from their intrinsic values, depending upon pH and buffer conditions.     

Effect of protein backbone folding on the stability of protein-ligand complexes

The role played by the degree of folding of protein backbones in explaining the binding energetics of protein-ligand interactions has been studied. We analyzed the protein/peptide interactions in the RNase-S system in which amino acids at two positions of the peptide S have been mutated. The global degree of folding of the protein S correlates in a significant way with the free energy and enthalpy of the protein-peptide interactions. A much better correlation is found with the local contribution to the degree of folding of one amino acid residue: Thr36. This residue is shown to have a destabilizing interaction with Lys41, which interacts directly with peptide S. Another system, consisting of the interactions of small organic molecules with HIV-1 protease was also studied. In this case, the global change in the degree of folding of the protease backbone does not explain the binding energetics of protein-ligand interactions. However, a significant correlation is observed between the free energy of binding and the contribution of two amino acid residues in the HVI-1 protease: Gly49 and Ile66. In general, it was observed that the changes in the degree of folding are not restricted to the binding site of the protein chain but are distributed along the whole protein backbone. This study provides a basis for further consideration of the degree of folding as a parameter for empirical structural parametrizations of the binding energetics of protein folding and binding.     

Quantitative study of protein association at picomolar concentrations: the lambda phage cl repressor

A method has been developed for radiolabeling the lambda cl repressor to a specific activity sufficiently high to permit accurate quantitation of the protein in the picomolar range of concentration. Procedures are described whereby the labeled protein can be used for accurate quantitative study of the energetics of repressor assembly by large zone analytical gel chromatography. This methodology is applicable to other systems in which the stoichiometry and energetics of tightly associating DNA binding proteins are currently difficult to measure.     

Streptavidin-biotin binding energetics

The high affinity energetics in the streptavidin-biotin system provide an excellent model system for studying how proteins balance enthalpic and entropic components to generate an impressive overall free energy for ligand binding. We review here concerted site-directed mutagenesis, biophysical, and computational studies of aromatic and hydrogen bonding interaction energetics between streptavidin and biotin. These results also have provided insight into how streptavidin builds a large activation barrier to dissociation by managing the enthalpic and entropic activation components. Finally, we review recent studies of the biotin dissociation pathway that address the fundamental question of how ligands exit protein binding pockets.     

Thermodynamics of zinc binding to hepatitis C virus NS3 protease: A folding by binding event

The hepatitis C virus (HCV) nonstructural protein 3 (NS3) protease is responsible for the processing of the non‐structural region of the viral precursor polyprotein in infected hepatic cells. HCV NS3 is a zinc‐dependent serine protease. The zinc ion, which is bound far away from the active site and considered to have a structural role, is essential for the structural integrity of the protein; furthermore, the ion is required for the hydrolytic activity. Consequently, the NS3 zinc binding site has been considered for a long time as a possible target for drug discovery. As a first step towards this goal, the energetics of the NS3‐zinc interaction and its effect on the NS3 conformation must be established and discussed. The thermodynamic characterization of zinc binding to NS3 protease by isothermal titration calorimetry and spectroscopy is presented here. Spectroscopic and calorimetric results suggest that a considerable conformational change in the protein is coupled to zinc binding. The energetics of the conformational change is comparable to that of the folding of a protein of similar size. Therefore, zinc binding to NS3 protease can be considered as a “folding by binding” event. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

Intrinsic Thermodynamics and Structure Correlation of Benzenesulfonamides with a Pyrimidine Moiety Binding to Carbonic Anhydrases I,II, VII,XII, and XIII

The early stage of drug discovery is often based on selecting the highest affinity lead compound. To this end the structural and energetic characterization of the binding reaction is important. The binding energetics can be resolved into enthalpic and entropic contributions to the binding Gibbs free energy. Most compound binding reactions are coupled to the absorption or release of protons by the protein or the compound. A distinction between the observed and intrinsic parameters of the binding energetics requires the dissection of the protonation/deprotonation processes. Since only the intrinsic parameters can be correlated with molecular structural perturbations associated with complex formation, it is these parameters that are required for rational drug design. Carbonic anhydrase (CA) isoforms are important therapeutic targets to treat a range of disorders including glaucoma, obesity, epilepsy, and cancer. For effective treatment isoform-specific inhibitors are needed. In this work we investigated the binding and protonation energetics of sixteen [(2-pyrimidinylthio)acetyl]benzenesulfonamide CA inhibitors using isothermal titration calorimetry and fluorescent thermal shift assay. The compounds were built by combining four sulfonamide headgroups with four tailgroups yielding 16 compounds. Their intrinsic binding thermodynamics showed the limitations of the functional group energetic additivity approach used in fragment-based drug design, especially at the level of enthalpies and entropies of binding. Combined with high resolution crystal structural data correlations were drawn between the chemical functional groups on selected inhibitors and intrinsic thermodynamic parameters of CA-inhibitor complex formation.     

A designed four-alpha-helix bundle that binds the volatile general anesthetic halothane with high affinity

The structural features of volatile anesthetic binding sites on proteins are being examined with the use of a defined model system consisting of a four-alpha-helix bundle scaffold with a hydrophobic core. Previous work has suggested that introducing a cavity into the hydrophobic core improves anesthetic binding affinity. The more polarizable methionine side chain was substituted for a leucine, in an attempt to enhance the dispersion forces between the ligand and the protein. The resulting bundle variant has an improved affinity (K(d) = 0.20 +/- 0.01 mM) for halothane binding, compared with the leucine-containing bundle (K(d) = 0.69 +/- 0.06 mM). Photoaffinity labeling with (14)C-halothane reveals preferential labeling of the W15 residue in both peptides, supporting the view that fluorescence quenching by bound anesthetic reports both the binding energetics and the location of the ligand in the hydrophobic core. The rates of amide hydrogen exchange were similar for the two bundles, suggesting that differences in binding affinity were not due to changes in protein stability. Binding of halothane to both four-alpha-helix bundle proteins stabilized the native folded conformations. Molecular dynamics simulations of the bundles illustrate the existence of the hydrophobic core, containing both W15 residues. These results suggest that in addition to packing defects, enhanced dispersion forces may be important in providing higher affinity anesthetic binding sites. Alternatively, the effect of the methionine substitution on halothane binding energetics may reflect either improved access to the binding site or allosteric optimization of the dimensions of the binding pocket. Finally, preferential stabilization of folded protein conformations may represent a fundamental mechanism of inhaled anesthetic action.     

Binding of TmHU to single dsDNA as observed by optical tweezers

Optical tweezers are employed to study the action of the histone-like protein from Thermotoga maritima (TmHU) on DNA at a single molecule level. Binding and disruption of TmHU to and from DNA are found to take place in discrete steps of 4-5 nm length and a net binding enthalpy of about 16kBT. This is in reasonable agreement with a microscopic model that estimates the extension of the binding sites of the protein and evaluates the energetics mainly for bending of the DNA in the course of interaction.     

Exploring the impact of polyproline II (PII) conformational bias on the binding of peptides to the SEM-5 SH3 domain

The left-handed polyproline II helical structure (P(II)) is observed to be a dominant conformation in the disordered states of protein and small polypeptide chains, even when no prolines are present in the sequence. Recently, in work by Ferreon and Hilser, the energetics associated with Ala and Gly substitutions at a surface exposed proline site were determined calorimetrically by measuring the binding energetics of Sos peptide variants to the C-terminal Src Homology 3 domain of SEM-5. The results were interpreted as a significant conformational bias toward the bound conformation (i.e., P(II)), even when the ligand is unbound. That study was not able to determine, however, whether the conformational bias of the peptides could be explained in terms other than that of a P(II) preference. Here, we test, using a computer algorithm based on the hard sphere collision (HSC) model, the notion of whether a bias in the unbound states of the peptide ligands is specific for the P(II) conformation, or if a bias to any other region of (phi, psi) space can also result in the same observed binding energetics. The results of these computer simulations indicate that, of the regions of (phi, psi) modeled for bias in the small peptides, only the bias to the P(II) conformation, and at rates of bias similar to the experimentally observed rates, quantitatively reproduced the experimental binding energetics.     

Exploring and designing protein function with restricted diversity

Combinatorial libraries with restricted diversity can be used to rapidly map binding energetics across protein interfaces. Shotgun scanning strategies have been used for alanine scanning and for alternative mutagenesis schemes that provide high-resolution functional views of binding interfaces. In addition, synthetic antibodies have been derived from na?ve libraries restricted to a binary code to explore the minimal requirements for molecular recognition. These studies shed light on the underlying principles governing molecular recognition, and provide rapid yet quantitative alternatives to conventional biophysical methods for exploring protein structure and function.     

Energetics of substrate binding and catalysis by class 1 (glycosylhydrolase family 47) alpha-mannosidases involved in N-glycan processing and endoplasmic reticulum quality control

Nascent glycoproteins are subject to quality control in the lumen of the endoplasmic reticulum (ER) where they can either be effectively folded with the aid of a collection of ER chaperones or they can be targeted for disposal in a process known as ER-associated degradation. Initiation of the ER disposal process involves selective trimming of N-glycans by ER alpha-mannosidase I and subsequent recognition by the ER degradation-enhancing alpha-mannosidase-like protein family of lectins, both members of glycosylhydrolase family 47. The kinetics and energetics of substrate binding and catalysis by members of this family were investigated here by the analysis of wild type and mutant forms of human ER alpha-mannosidase I. The contributions of several amino acid residues and an enzyme-associated Ca(2+) ion to substrate binding and catalysis were demonstrated by a combination of surface plasmon resonance and enzyme kinetic analyses. One mutant, E330Q, shown previously to alter general acid function within the catalytic site, resulted in an enzyme that possessed increased glycan binding affinity but compromised glycan hydrolysis. This mutant protein was used in a series of glycan binding studies with a library of mannose-containing ligands to examine the energetics of Man(9)GlcNAc(2) substrate interactions. These studies provide a framework for understanding the nature of the unusual substrate interactions within the family 47 mannosidases involved in glycan maturation and ER-associated glycoprotein degradation.     

On the link between conformational changes,ligand binding and heat capacity

BackgroundConformational changes coupled to ligand binding constitute the structural and energetics basis underlying cooperativity, allostery and, in general, protein regulation. These conformational rearrangements are associated with heat capacity changes. ITC is a unique technique for studying binding interactions because of the simultaneous determination of the binding affinity and enthalpy, and for providing the best estimates of binding heat capacity changes.Scope of reviewStill controversial issues in ligand binding are the discrimination between the “conformational selection model” and the “induced fit model”, and whether or not conformational changes lead to temperature dependent apparent binding heat capacities. The assessment of conformational changes associated with ligand binding by ITC is discussed. In addition, the “conformational selection” and “induced fit” models are reconciled, and discussed within the context of intrinsically (partially) unstructured proteins.Major conclusionsConformational equilibrium is a major contribution to binding heat capacity changes. A simple model may explain both conformational selection and induced fit scenarios. A temperature-independent binding heat capacity does not necessarily indicate absence of conformational changes upon ligand binding. ITC provides information on the energetics of conformational changes associated with ligand binding (and other possible additional coupled equilibria).General significancePreferential ligand binding to certain protein states leads to an equilibrium shift that is reflected in the coupling between ligand binding and additional equilibria. This represents the structural/energetic basis of the widespread dependence of ligand binding parameters on temperature, as well as pH, ionic strength and the concentration of other chemical species. This article is part of a Special Issue entitled Microcalorimetry in the BioSciences — Principles and Applications, edited by Fadi Bou-Abdallah.     

Energetics for displacing a single chain from the surface of microcrystalline cellulose into the active site of Acidothermus cellulolyticus Cel5A

A series of molecular mechanics calculations were used to analyze the energetics for moving a single polysaccharide chain from the surface of microcrystalline cellulose into the binding cleft of the Cel5A cellulase from Acidothermus cellulolyticus. A build-up procedure was used to model the placement of a 12-residue oligosaccharide chain along the surface of the enzyme, using as a guide the four residues of the tetrasaccharide substrate co-crystallized with the protein in the crystallographic structure determination. The position of this 12-residue oligosaccharide was used to orient the enzyme properly above two different surfaces of cellulose 1beta, the (1,0,0) and the (1,1,0) faces of the crystal. Constrained molecular dynamics simulations were then used to pull a target chain directly below the enzyme up out of the crystal surface and into the binding groove. The energetics for this process were favorable for both faces, with the step face being more favorable than the planar face, implying that this surface could be hydrolyzed more readily.     

Energetics-based protein profiling on a proteomic scale: identification of proteins resistant to proteolysis

Native states of proteins are flexible, populating more than just the unique native conformation. The energetics and dynamics resulting from this conformational ensemble are inherently linked to protein function and regulation. Proteolytic susceptibility is one feature determined by this conformational energy landscape. As an attempt to investigate energetics of proteins on a proteomic scale, we challenged the Escherichia coli proteome with extensive proteolysis and determined which proteins, if any, have optimized their energy landscape for resistance to proteolysis. To our surprise, multiple soluble proteins survived the challenge. Maltose binding protein, a survivor from thermolysin digestion, was characterized by in vitro biophysical studies to identify the physical origin of proteolytic resistance. This experimental characterization shows that kinetic stability is responsible for the unusual resistance in maltose binding protein. The biochemical functions of the identified survivors suggest that many of these proteins may have evolved extreme proteolytic resistance because of their critical roles under stressed conditions. Our results suggest that under functional selection proteins can evolve extreme proteolysis resistance by modulating their conformational energy landscapes without the need to invent new folds, and that proteins can be profiled on a proteomic scale according to their energetic properties by using proteolysis as a structural probe.     

Energetics of RNA binding by the West Nile virus RNA triphosphatase

The West Nile virus (WNV) RNA genome harbors the characteristic methylated cap structure present at the 5' end of eukaryotic mRNAs. In the present study, we report a detailed study of the binding energetics and thermodynamic parameters involved in the interaction between RNA and the WNV RNA triphosphatase, an enzyme involved in the synthesis of the RNA cap structure. Fluorescence spectroscopy assays revealed that the initial interaction between RNA and the enzyme is characterized by a high enthalpy of association and that the minimal RNA binding site of NS3 is 13 nucleotides. In order to provide insight into the relationship between the enzyme structure and RNA binding, we also correlated the effect of RNA binding on protein structure using both circular dichroism and denaturation studies as structural indicators. Our data indicate that the protein undergoes structural modifications upon RNA binding, although the interaction does not significantly modify the stability of the protein.     

Intramolecular interactions at protein surfaces and their impact on protein function

Protein surfaces play a key role in the biological function of proteins. Consequently, structural features of protein surfaces are the basis for predicting function from structure. A well-established principle of binding by proteins is that ligands must compete with water and other small molecules to form interactions with protein surfaces. A less obvious issue, and the emphasis of this article, is that ligands must also compete with interactions among residues at protein surfaces. Results from structural surveys, a variety of experimental studies and computations suggest that intramolecular interactions are present at protein surfaces and that the energetics of these interactions can change when proteins bind to other molecules.     

Interaction of heparin with annexin V

The energetics and kinetics of the interaction of heparin with the Ca2+ and phospholipid binding protein annexin V, was examined and the minimum oligosaccharide sequence within heparin that binds annexin V was identified. Affinity chromatography studies confirmed the Ca2+ dependence of this binding interaction. Analysis of the data obtained from surface plasmon resonance afforded a Kd of approximately 21 nM for the interaction of annexin V with end-chain immobilized heparin and a Kd of approximately 49 nM for the interaction with end-chain immobilized heparan sulfate. Isothermal titration calorimetry showed the minimum annexin V binding oligosaccharide sequence within heparin corresponds to an octasaccharide sequence. The Kd of a heparin octasaccharide binding to annexin V was approximately 1 microM with a binding stoichiometry of 1:1.     

Energetics of solvent and ligand-induced conformational changes in alpha-lactalbumin

The energetics of structural changes in the holo and apo forms of a-lactalbumin and the transition between their native and denatured states induced by binding Ca2+ and Na+ have been studied by differential scanning and isothermal titration microcalorimetry and circular dichroism spectroscopy under various solvent conditions. Removal of Ca2+ from the protein enhances its sensitivity to pH and ionic conditions due to noncompensated negative charge-charge interactions at the cation binding site, which significantly reduces its overall stability. At neutral pH and low ionic strength, the native structure of apo-alpha-lactalbumin is stable below 14 C and undergoes a conformational change to a native-like molten globule intermediate at temperatures above 25 degrees C. The denaturation of either holo- or apo-alpha-lactalbumin is a highly cooperative process that is characterized by an enthalpy of similar magnitude when calculated at the same temperature. Measured by direct calorimetric titration, the enthalpy of Ca2+-binding to apo-LA at pH 7.5 is -7.1 kJ mol(-1) at 5.0 degrees C. which is essentially invariant to protonation effects. This small enthalpy effect infers that stabilization of alpha-lactalbumin by Ca2+ is primarily an entropy driven process, presumably arising from electrostatic interactions and the hydration effect. In contrast to the binding of calcium, the interaction of sodium with apo-LA does not produce a noticeable heat effect and is characterized by its ionic nature rather than specific binding to the metal-binding site. Characterization of the conformational stability and ligand binding energetics of alpha-lactalbumin as a function of solvent conditions furnishes significant insight regarding the molecular flexibility and regulatory mechanism mediated by this protein.