Volume of Hsp90 ligand binding and the unfolding phase diagram as a function of pressure and temperature

Volume changes that accompany protein unfolding and ligand binding are important but largely neglected thermodynamic parameters that may facilitate rational drug design. Here, we determined the volume of lead compound ICPD47 binding to an anticancer target, heat shock protein 90 N-terminal domain, using a pressure shift assay (PressureFluor). The ligand exhibited a stabilizing effect on the protein by increasing its melting pressure and temperature. The Gibbs free energy of unfolding depends on the absence or presence of ligand and has an elliptical shape. Ellipse size increases upon addition of the strongly binding ligand, which stabilizes the protein. The three-dimensional (3D) ellipsoidal surface of the Gibbs free energy of unfolding was calculated with increasing ligand concentrations. The negative volume of ligand binding was relatively large and significantly exceeded the volume of protein unfolding. The pressure shift assay technique could be used to determine the volume changes associated with both protein unfolding as well as ligand binding to protein.     

Identification of two transactivation domains in the mouse oestrogen receptor.

We have identified two discrete transactivation domains within the mouse oestrogen receptor whose relative activities vary according to the target promoter. One domain lies within the N-terminal region and is active in the absence of oestradiol. The second domain is contained within the C-terminal portion of the protein and depends upon oestrogen binding for its activity. The location and oestrogen dependence of this domain has been confirmed using chimaeric receptors containing the Lex A DNA binding domain. Although transactivation by the C-terminal domain is dependent upon ligand binding the analysis of receptor deletion mutants has demonstrated that these two functions are not entirely coincident.     

Unusual shapes of absorption isotherms in three-component systems

It has been shown that the shape of Scatchard isotherms upon competitive binding of two ligands to the same binding site in the three-component ligand 1-ligand 2-DNA system depends crucially on the binding constant values. The binding isotherm of ligand 2 in the presence of the competitive ligand 1 turns back (has a bow-like form) when the binding constant of the first ligand is larger than the binding constant of the second one.     

Identification of constitutive and steroid-dependent transactivation domains in the mouse oestrogen receptor

We have identified two transactivation domains in the mouse oestrogen receptor whose activities depend on the target promoter. The major domain is contained within the C-terminal portion of the protein and depends upon oestrogen binding for its activity. The location and oestrogen dependence of this domain has been confirmed using chimaeric receptors containing the Lex A DNA binding domain. Although transactivation by the C-terminal domain is dependent upon ligand binding the analysis of receptor deletion mutants has demonstrated that these two functions are not entirely coincident. The second transactivation domain lies within the N-terminal region and is active in the absence of oestradiol. The differences in oestrogen requirement for the activity of the two transactivation domains may account for the partial agonist activity of certain antihormones.     

Mechanisms for activation and antagonism of an AMPA-sensitive glutamate receptor: crystal structures of the GluR2 ligand binding core

Crystal structures of the GluR2 ligand binding core (S1S2) have been determined in the apo state and in the presence of the antagonist DNQX, the partial agonist kainate, and the full agonists AMPA and glutamate. The domains of the S1S2 ligand binding core are expanded in the apo state and contract upon ligand binding with the extent of domain separation decreasing in the order of apo > DNQX > kainate > glutamate approximately equal to AMPA. These results suggest that agonist-induced domain closure gates the transmembrane channel and the extent of receptor activation depends upon the degree of domain closure. AMPA and glutamate also promote a 180 degrees flip of a trans peptide bond in the ligand binding site. The crystal packing of the ligand binding cores suggests modes for subunit-subunit contact in the intact receptor and mechanisms by which allosteric effectors modulate receptor activity.     

Order changes within receptor systems upon ligand binding: receptor tightening/oligomerisation and the interpretation of binding parameters

Recent hydrogen-deuterium exchange experiments have highlighted tightening and loosening of protein structures upon ligand binding, with changes in bonding (DeltaH) and order (DeltaS) which contribute to the overall thermodynamics of ligand binding. Tightening and loosening show that ligand binding respectively stabilises or destabilises the internal structure of the protein, i.e. it shows positive or negative cooperativity between ligand binding and the receptor structure. In the case of membrane-bound receptors, such as G protein-coupled receptors (GPCRs) and ligand gated ion channel receptors (LGICRs), most binding studies have focussed on association/dissociation constants. Where these have been broken down into enthalpic and entropic contributions, the phenomenon of "thermodynamic discrimination" between antagonists and agonists has often been noted; e.g. for a receptor where agonist binding is predominantly enthalpy driven, antagonist binding is predominantly entropy driven and vice versa. These data have not previously been considered in terms of the tightening, or loosening, of receptor structures that respectively occurs upon positively, or negatively, cooperative binding of ligand. Nor have they been considered in light of the homo- and hetero-oligomerisation of GPCRs and the possibility of ligand-induced changes in oligomerisation. Here, we argue that analysis of the DeltaH and DeltaS of ligand binding may give useful information on ligand-induced changes in membrane-bound receptor oligomers, relevant to the differing effects of agonists and antagonists.     

Mathematical analysis of ligand-induced monomerization and dimerization in the monomer dimer equilibrium of proteins. Application to D-amino acid oxidase.

We have mathematically analyzed ligand-induced monomerization and dimerization in a protein monomer-dimer equilibrium system, in which the monomer has one and the dimer two binding sites. These dimer sites have the same binding constants for the first ligand but may cooperatively interact when one of them is occupied by a ligand molecule. In this system, the apparent dimerization constant and the apparent molecular weight are functions of free ligand concentration, and depend on the intrinsic binding constants of the ligand molecule to the monomer and the dimer. The behavior of these functions is classified into 17 cases according to the values of the three intrinsic binding constants, and some calculated examples are shown graphically for selected parameters. The theory was also applied to D-amino acid oxidase [EC 1.4.3.3], a flavoprotein, and the pH dependence of the apparent dimerization constant and the apparent molecular weight in the presence of ligand, p-aminobenzoate, were studied theoretically using parameters obtained in our previous experiments (5).     

Structural and thermodynamic studies of binding saturated fatty acids to bovine β-lactoglobulin

Lactoglobulin is a globular milk protein for which physiological function has not been clarified. Due to its binding properties lactoglobulin might serve as a carrier for bioactive molecules. Binding of 12-, 14-, 16- and 18-carbon saturated fatty acids to bovine β-lactoglobulin has been characterised by isothermal titration calorimetry and X-ray crystallography as a part of systematic studies of lactoglobulin complexes with ligands of biological importance. The thermodynamic parameters have been determined for lauric, myristic and palmitic acid complexes revealing systematic decrease of enthalpic and increase of entropic component of ΔG with elongation of aliphatic chain. In all crystal structures determined with resolution 1.9-2.1?, single fatty acid molecule was found in the β-barrel in extended conformation with individual pattern of interactions. Location of a fatty acid in the binding site depends on the length of aliphatic chain and influences polar interactions between protein and ligand. Systematic changes of entropic component indicate important role of water in binding process.     

Pre-existing soft modes of motion uniquely defined by native contact topology facilitate ligand binding to proteins

Modeling protein flexibility constitutes a major challenge in accurate prediction of protein-ligand and protein-protein interactions in docking simulations. The lack of a reliable method for predicting the conformational changes relevant to substrate binding prevents the productive application of computational docking to proteins that undergo large structural rearrangements. Here, we examine how coarse-grained normal mode analysis has been advantageously applied to modeling protein flexibility associated with ligand binding. First, we highlight recent studies that have shown that there is a close agreement between the large-scale collective motions of proteins predicted by elastic network models and the structural changes experimentally observed upon ligand binding. Then, we discuss studies that have exploited the predicted soft modes in docking simulations. Two general strategies are noted: pregeneration of conformational ensembles that are then utilized as input for standard fixed-backbone docking and protein structure deformation along normal modes concurrent to docking. These studies show that the structural changes apparently "induced" upon ligand binding occur selectively along the soft modes accessible to the protein prior to ligand binding. They further suggest that proteins offer suitable means of accommodating/facilitating the recognition and binding of their ligand, presumably acquired by evolutionary selection of the suitable three-dimensional structure.     

Calculations of DNA Condensation Caused by Ligand Adsorption

A method for calculating the curves of DNA transition from linear to condensed state upon binding of condensing ligands has been developed. The character of the transition and ligand concentration necessary for condensation have been shown to be governed by the length of DNA molecule, energy and stoichiometry parameters of the DNA–ligand complex (equilibrium constant between linear and condensed form in the absence of ligands, constants for ligand binding to linear and condensed forms, the number of base pairs covered by one ligand, etc.). The results of the calculations indicate that a slight difference in the free energies of these DNA states (less than 6 cal/mol(bp) for a DNA of 500 bp) is sufficient for the existence of a stable linear state in the absence of ligands (in free DNA) and the formation of stable condensed state upon complexation.     

Two-dimensional differential scanning calorimetry: simultaneous resolution of intrinsic protein structural energetics and ligand binding interactions by global linkage analysis.

A general theoretical development for the design and analysis of two-dimensional thermal stability surfaces of proteins is presented. The surfaces are generated from multiple excess heat capacity profiles ( vs T) obtained at varying concentrations of an interacting ligand. The energetics of both the intrinsic protein stability and the protein-ligand interaction are simultaneously resolved by employing statistical thermodynamic models in global linkage analysis. This formalism allows resolution of the intrinsic protein folding-unfolding parameters (enthalpy, entropy, and heat capacity changes) as well as the ligand interaction parameters (binding stoichiometry, enthalpy, entropy, and heat capacity changes). The theory has been applied to the case of ribonuclease A and its interaction with cytidine-2'-monophosphate. The accuracy of the thermodynamic parameters obtained by this approach compares within error with those parameters that can be obtained by direct measurements.     

The forward rate of binding of surface-tethered reactants: effect of relative motion between two surfaces

The reaction of molecules confined to two dimensions is of interest in cell adhesion, specifically for the reaction between cell surface receptors and substrate-bound ligand. We have developed a model to describe the overall rate of reaction of species that are bound to surfaces under relative motion, such that the Peclet number is order one or greater. The encounter rate between reactive species is calculated from solution of the two-dimensional convection-diffusion equation. The probability that each encounter will lead to binding depends on the intrinsic rate of reaction and the encounter duration. The encounter duration is obtained from the theory of first passage times. We find that the binding rate increases with relative velocity between the two surfaces, then reaches a plateau. This plateau indicates that the increase in the encounter rate is counterbalanced by the decrease in the encounter duration as the relative velocity increases. The binding rate is fully described by two dimensionless parameters, the Peclet number and the Damk?hler number. We use this model to explain data from the cell adhesion literature by incorporating these rate laws into "adhesive dynamics" simulations to model the binding of a cell to a surface under flow. Leukocytes are known to display a "shear threshold effect" when binding selectin-coated surfaces under shear flow, defined as an increase in bind rate with shear; this effect, as calculated here, is due to an increase in collisions between receptor and ligand with increasing shear. The model can be used to explain other published data on the effect of wall shear rate on the binding of cells to surfaces, specifically the mild decrease in binding within a fixed area with increasing shear rate.     

Calculation of DNA condensation, caused by adsorption of ligands

A method for calculating the curves of DNA transition from linear to condensed state upon binding of condensing ligands has been developed. The character of the transition and ligand concentration necessary for condensation have been shown to be governed by the length of DNA molecule, energy and stoichiometry parameters of DNA-ligand complex (equilibrium constant between linear and condensed form in the absence of ligands, constants for ligand binding to linear and condensed forms, the number of base pairs covered by one ligand, etc.). The results of the calculations indicate that only slight difference in the free energies of these states in free DNA (less than 6 cal/mole(bp) for DNA of 500 bp long) is sufficient for the existence of stable linear state in the absence of ligands (in free DNA) and the formation of stable condensed state upon complexation.     

Ligand-induced monomerization of <Emphasis Type="Italic">Allochromatium vinosum</Emphasis> cytochrome <Emphasis Type="Italic">c</Emphasis>′ studied using native mass spectrometry and fluorescence resonance energy transfer

Cytochrome c' from Allochromatium vinosum is an attractive model protein to study ligand-induced conformational changes. This homodimeric protein dissociates into monomers upon binding of NO, CO or CN(-) to the iron of its covalently attached heme group. While ligand binding to the heme has been well characterized using a variety of spectroscopic techniques, direct monitoring of the subsequent monomerization has not been reported previously. Here we have explored two biophysical techniques to simultaneously monitor ligand binding and monomerization. Native mass spectrometry allowed the detection of the dimeric and monomeric forms of cytochrome c' and even showed the presence of a CO-bound monomer. The kinetics of the ligand-induced monomerization were found to be significantly enhanced in the gas phase compared with the kinetics in solution, however. Ligand binding to the heme and the dissociation of the dimer in solution were also studied using energy transfer from a fluorescent probe to both heme groups of the protein. Comparison of ligand binding kinetics as observed with UV-vis spectroscopy with changes in fluorescence suggested that binding of one CO molecule per dimer could be sufficient for monomerization.     

Selected-fit versus induced-fit protein binding: kinetic differences and mutational analysis

The binding of a ligand molecule to a protein is often accompanied by conformational changes of the protein. A central question is whether the ligand induces the conformational change (induced-fit), or rather selects and stabilizes a complementary conformation from a pre-existing equilibrium of ground and excited states of the protein (selected-fit). We consider here the binding kinetics in a simple four-state model of ligand-protein binding. In this model, the protein has two conformations, which can both bind the ligand. The first conformation is the ground state of the protein when the ligand is off, and the second conformation is the ground state when the ligand is bound. The induced-fit mechanism corresponds to ligand binding in the unbound ground state, and the selected-fit mechanism to ligand binding in the excited state. We find a simple, characteristic difference between the on- and off-rates in the two mechanisms if the conformational relaxation into the ground states is fast. In the case of selected-fit binding, the on-rate depends on the conformational equilibrium constant, whereas the off-rate is independent. In the case of induced-fit binding, in contrast, the off-rate depends on the conformational equilibrium, while the on-rate is independent. Whether a protein binds a ligand via selected-fit or induced-fit thus may be revealed by mutations far from the protein's binding pocket, or other "perturbations" that only affect the conformational equilibrium. In the case of selected-fit, such mutations will only change the on-rate, and in the case of induced-fit, only the off-rate.     

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.     

A dynamic dialysis method for studying protein-ligand binding using chromatographic theory

A new dynamic dialysis method has been developed for studying protein-ligand binding phenomena. The method depends on analysis of the elution pattern of ligand in a single dialyzing process where the ligand concentration in the sample compartment changes greatly with time. The dialyzer is composed of a long, narrow chamber (the sample compartment) between two sheets of semipermeable membrane and two outside chambers (the sink compartment) connected as a single path. Eluting buffer flows in the sink compartment to exchange the ligand with the solution in the sample compartment. Therefore, the ligand concentration gradient in the sink compartment is in the longitudinal direction. The mathematical expressions to analyze the experimental data were derived from a modified theory of chromatography. Examination of the binding of sulfanilamide to bovine serum albumin using this method shows that these equations are valid for use in studying protein-ligand binding.     

Kinetics of ligand binding and quaternary conformational change in the homodimeric hemoglobin from Scapharca inaequivalvis

The kinetics of the reaction with oxygen and carbon monoxide of the homodimeric hemoglobin from the bivalve mollusc Scapharca inaequivalvis has been extensively investigated by flash and dye-laser photolysis, temperature jump relaxation, and stopped flow methods. The results indicate that cooperativity in ligand binding, already observed for oxygen at equilibrium, finds its kinetic counterpart in a large decrease of the oxygen dissociation velocity in the second step of the binding reaction. In the case of carbon monoxide, cooperativity is clearly evident in the increase of the combination velocity constant as the reaction proceeds. Therefore, the ligand-binding kinetics of this dimeric hemoglobin shows the characteristic features of the corresponding reactions of tetrameric hemoglobins. Analysis of the data in terms of the allosteric model proposed by Monod et al. (Monod, J., Wyman, J., and Changeux, J. P. (1965) J. Mol. Biol. 12, 88-118) has shown that the values of the allosteric parameters cannot be fixed uniquely for a dimeric hemoglobin. The rapid changes in absorbance observed at the isosbestic points of unliganded and liganded hemoglobin following laser photolysis provided a value of 7 X 10(4) S-1 at 20 degrees C for the rate of the ligand-free quarternary conformational change, postulated on the basis of cooperative ligand binding. Comparison of the rapid absorbance changes observed during ligand rebinding in this hemoglobin with those observed in tuna hemoglobin indicate that, at full photolysis, binding to the T state is followed by further binding and conversion to the liganded R state; at partial photolysis, population of the liganded T state occurs immediately and is followed by a decay to the liganded R state upon further ligand binding. These new results, in conjunction with previous equilibrium data on the same system, show unequivocally that the presence of two different types of chain is not an absolute prerequisite for cooperativity in hemoglobins, contrary to currently accepted ideas.     

Decoupling of melting domains in immobilized ribonuclease A

Ribonuclease A has been immobilized on silica beads through glutaraldeyde-mediated chemical coupling in order to improve the stability of the protein against thermal denaturation. The thermodynamic and binding properties of the immobilized enzyme have been studied and compared with those of the free enzyme. The parameters describing the binding of the inhibitor 3′ -CMP (K_a and ΔH) as monitored by spectrophotometry and calorimetry were not significantly affected after immobilization. Conversely both the stability and unfolding mechanism drastically changed. Thermodynamic analysis of the DSC data suggests that uncoupling of protein domains has occurred as a consequence of the immobilization. The two state approximation of the protein unfolding process is not longer valid for the immobilized RNase. Protein stability strongly depends on the hydrophobicity properties of the support surface as well as on the presence of the inhibitor and pH. For example, after immobilization on a highly hydrophobic surface, the enzyme is partially in the unfolded state. The binding of a ligand is able to reorganize the protein structure into a native-like conformation. The refolding rates are different for the two protein domains and vary as a function of pH and presence of the inhibitor 3′-CMP. © 1994 Wiley-Liss, Inc.     

Simultaneous recognition of nucleobase and sites of DNA damage: effect of tethered cation on the binding affinity

BACKGROUND: The 3,5-diamino-N-(3-aminopropyl)-6-chloropyrazine-2-carboxamide (DCPC-NH(2)) has been synthesized and characterized by Mass and (1)H NMR. The selective binding of the ligand to thymine (T) target base is investigated by the melting temperature (T(m)) and fluorescence measurements. METHODS: Thermal denaturation study of DNA duplex containing T target base revealed the DeltaT(m) of 5.1 degrees C, while least influence was observed for other target bases. The fluorescence of the ligand DCPC-NH(2) is quenched only upon adding the DNA containing T target base. RESULTS: The binding constant for the interaction of the ligand to T target base containing DNA duplex was determined to be 4.7 (+/-0.3)x10(6) M(-1). The tethered cation in the ligand is found to enhance the binding constant. The ligand binds to both a target nucleotide and an AP site on the complimentary strand for the target strand in a DNA duplex. GENERAL SIGNIFICANCE: Interestingly, the electronic behavior of the ligand depends on the bases flanking the AP site. Its fluorescence is quenched with guanine flanking bases, while it is enhanced with DNA duplex containing T bases flanking an AP site. Finally, the binding modes were visualized by molecular modeling.