Interaction of protein with a self-associating ligand. Deviation from a hyperbolic binding curve and the appearance of apparent co-operativity in the Scatchard plot

We derived a general formula to analyze a binding system in which a ligand self-associates, in terms of experimentally determinable quantities, i.e. r, the average number of bound ligands per protein molecule, and Lft, the total free ligand concentration, which are expressed as a ligand monomer unit. The limiting behaviors of the Scatchard plot (r/Lft vs r plot), that is, the intercepts on the r-axis and the r/Lft-axis, and the limiting slopes, are generally given. Three models that may be encountered are considered in detail. Numerical examples are also presented to illustrate how the self-association of a ligand affects the binding curves. The ligand self-association alone can cause deviation of the profile of the binding curve (r vs Lft plot) from a hyperbola, resulting in a nonlinear Scatchard plot. Therefore, analysis of the binding data without consideration of ligand self-association may lead to erroneous conclusions as to the numbers and classes of binding sites, co-operativity among the sites and binding parameter values.     

A simple and rapid orientation procedure for the estimation of affinities (1/Kd) and maximum binding capacities (Bmax) of ligand binding to two receptor subpopulations

This paper describes a simple and rapid procedure for the estimation of specific parameters (dissociation constants, Kd and maximum binding capacities, Bmax) of ligand binding to two receptor subpopulations. This procedure provides, in a few minutes, the investigator, performing the actual binding studies, the necessary information about receptor heterogeneity, enabling the investigator to plan further experiments. The procedure is based on the graphical comparison of experimental binding data (ligand binding to one or two receptor subpopulations) with the theoretical values of ligand binding to one receptor population at four levels). The values of Kd and Bmax for high- and low-affinity receptors are derived from 4 horizontal deviations of experimental data from a theoretical data plot at these levels by their comparison with tabulated deviations. The correctness of the estimated parameters can be confirmed by the comparison of experimental data with those simulated on the basis of applying the values of Kd and Bmax found in the formula for ligand binding to two receptor subpopulations. The practical applicability of this procedure was demonstrated both on simulated and experimental data, and confirmed by the well known computer programs for evaluating receptor heterogeneity, namely "LIGAND" and "Affinity spectra".     

Ligand-receptor interaction rates in the presence of convective mass transport.

The rate of binding of a ligand to receptors on the cell surface can be diffusion limited. We analyze the kinetics of binding, diffusion-limited in a stationary liquid, in the presence of convective mass transport. We derive a formula that expresses the reaction kinetics in terms of the mass transfer coefficient. A moderately transport-limited kinetics is not readily recognizable from the shape of the binding curve and may lead to erroneous estimates of the rate coefficients. We apply our results to practically important cases: a cell suspension in a stirred volume of liquid and a confluent cell colony under a laminar stream. Using typical numbers characterizing the ligand-receptor interactions, we show that stirring and perfusion can be important factors determining the reaction rates. With the confluent colony, the early reaction kinetics requires a different treatment, and we provide it for the case of low receptor occupancy. We show that, even with a fast perfusion, a cell monolayer can transiently generate a zone of depletion of the ligand, and that would affect the early stages of the reaction. Our results are expressed in a simple analytical form and can be used for the design and interpretation of experimental data.     

A general method of analysis of ligand-macromolecule equilibria using a spectroscopic signal from the ligand to monitor binding. Application to Escherichia coli single-strand binding protein-nucleic acid interactions

We describe a general method for the analysis of ligand-macromolecule binding equilibria for cases in which the interaction is monitored by a change in a signal originating from the ligand. This method allows the absolute determination of the average degree of ligand binding per macromolecule without any assumptions concerning the number of modes or states for ligand binding or the relationship between the fractional signal change and the fraction of bound ligand. Although this method is generally applicable to any type of signal, we discuss the details of the method as it applies to the analysis of binding data monitored by a change in fluorescence of a ligand upon binding to a nucleic acid. We apply the analysis to the equilibrium binding of Escherichia coli single-strand binding (SSB) protein to single-stranded nucleic acids, which is monitored by the quenching of the intrinsic tryptophan fluorescence of the SSB protein. With this method, one can quantitatively determine the relationship between the fractional signal change of the ligand and the fraction of bound ligand, LB/LT, and rigorously test whether the signal change is directly proportional to LB/LT. For E. coli SSB protein binding to single-stranded nucleic acids in its (SSB)65 binding mode [Lohman, T. M., & Overman, L. B. (1985) J. Biol. Chem. 260, 3594; Chrysogelos, S., & Griffith, J. (1982) Proc. Natl. Acad. Sci. U.S.A. 79, 5803], we show that the fractional quenching of the SSB fluorescence is equal to the fraction of bound SSB.     

DynaBiS: A hierarchical sampling algorithm to identify flexible binding sites for large ligands and peptides

Knowing the ligand or peptide binding site in proteins is highly important to guide drug discovery, but experimental elucidation of the binding site is difficult. Therefore, various computational approaches have been developed to identify potential binding sites in protein structures. However, protein and ligand flexibility are often neglected in these methods due to efficiency considerations despite the recognition that protein–ligand interactions can be strongly affected by mutual structural adaptations. This is particularly true if the binding site is unknown, as the screening will typically be performed based on an unbound protein structure. Herein we present DynaBiS, a hierarchical sampling algorithm to identify flexible binding sites for a target ligand with explicit consideration of protein and ligand flexibility, inspired by our previously presented flexible docking algorithm DynaDock. DynaBiS applies soft-core potentials between the ligand and the protein, thereby allowing a certain protein–ligand overlap resulting in efficient sampling of conformational adaptation effects. We evaluated DynaBiS and other commonly used binding site identification algorithms against a diverse evaluation set consisting of 26 proteins featuring peptide as well as small ligand binding sites. We show that DynaBiS outperforms the other evaluated methods for the identification of protein binding sites for large and highly flexible ligands such as peptides, both with a holo or apo structure used as input.     

Mutational analysis of the ligand binding domain of the low density lipoprotein receptor

The ligand binding domain of the low density lipoprotein (LDL) receptor contains seven imperfect repeats of a 40-amino acid cysteine-rich sequence. Each repeat contains clustered negative charges that have been postulated as ligand-binding sites. The adjacent region of the protein, the growth factor homology region, contains three cysteine-rich repeats (A-C) whose sequence differs from those in the ligand binding domain. To dissect the contribution of these different cysteine-rich repeats to ligand binding, we used oligonucleotide-directed mutagenesis to alter expressible cDNAs for the human LDL receptor which were then introduced into monkey COS cells by transfection. We measured the ability of the mutant receptors to bind LDL, which contains a single protein ligand for the receptor (apoB-100), and beta-migrating very low density lipoprotein (beta-VLDL), which contains apoB-100 plus multiple copies of another ligand (apoE). The results show that repeat 1 is not required for binding of either ligand. Repeats 2 plus 3 and repeats 6 plus 7 are required for maximal binding of LDL, but not beta-VLDL. Repeat 5 is required for binding of both ligands. Repeat A in the growth factor homology region is required for binding of LDL, but not beta-VLDL. Repeat B is not required for ligand binding. These results support a model for the LDL receptor in which various repeats play additive roles in ligand binding, each repeat making a separate contribution to the binding event.     

Theoretical Analysis of the Footprinting Experiment

Abstract

In the footprinting experiment, an end-radiolabeled DNA restriction fragment is subjected to digest by an endonuclease in the presence and absence of a ligand which alters the endonuclease cleavage rate at sites of ligand-DNA contact. The location of these sites, and the strength of the ligand binding, are then deduced from the measured concentrations of the different oligonucleotides produced by the digest. We analyze the experiment in terms of coupled kinetic equations which take into account the cutting rates of endonuclease for sites with ligand present and absent, and the rates of binding and dissociation of the ligand to a site. As long as the ligand concentration remains essentially constant (which occurs, for example, if digest is terminated early enough to assure that all fragments result from single cuts by the endonuclease), the oligonucleotide concentrations reflect only the ligand binding equilibrium constant (ratio of rate constants) and the cutting rates in the presence and absence of ligand. We also show how the measured oligonucleotide concentrations (from, e.g. an autoradiogram) can be used to deduce the ligand equilibrium binding constants for the various sites on the polymer.     

Ligand binding properties of the human erythropoietin receptor extracellular domain expressed in Escherichia coli.

We developed an assay to directly measure the ligand binding properties of the cloned human erythropoietin receptor (EpoR). The cDNA encoding the extracellular domain of the human EpoR was amplified by polymerase chain reaction and ligated into the prokaryotic expression vector pGEX3X. Synthesis in Escherichia coli was induced and a soluble glutathione S-transferase fusion protein, EREx, was purified by erythropoietin affinity chromatography. Purified EREx was bound to GSH agarose beads and used in a solid phase ligand binding assay. Specific binding of 125I-erythropoietin to EREx beads was demonstrated. A single affinity class (Kd = 1.5 nM) of the binding site was evident on Scatchard analysis. The Kd of this site is quantitatively equivalent to that of the "low" affinity cellular binding site. Kinetic analysis of ligand binding to EREx revealed both the on and off rates to be rapid, with t1/2 of 60 and 40 s, respectively. EREx ligand binding exhibits no obvious metal ion dependence or cross-competition by other hemopoietins. Antibodies to EREx block the binding of erythropoietin to the cellular EpoR. We conclude that the 66-kDa EpoR protein is capable of specific ligand binding and that no covalent modifications or associated molecules are required for this interaction. We speculate that the "high" affinity cellular binding site (Kd less than 0.2 nM) results from the interaction of the EpoR with another molecule, either additional EpoR or associated subunits, that decreases the ligand off rate.     

19F NMR studies of tryptophan/serum albumin binding

19F NMR provides direct measures of the Trp binding avidity of 'fatty acid free' bovine serum albumin when D- and L-6-fluorotryptophan are used as the probes. Both a high and low affinity binding site are present. The addition of octanoate either displaces the ligand from both sites or greatly decreases the affinity such that little binding occurs at 2 mM levels. In the case of L-6-fluorotryptophan separate signals are observed for the high and low affinity binding sites and titrations with competing ligands can be used to establish the relative affinities of ligands at the high affinity site. Binding at this site appears to be hydrophobic and shape specific with L-Phe being a very poor ligand (K(D)[L-Phe]/K(D)[L-Trp]=800) while both GHKalphaNal and GHKW displace L-6-fluorotryptophan from this site. In tripeptides of the general formula GHK[ epsilon NH(CH(2))(n)(CO)W], affinity increases with tether length and binding at the low affinity site is restored. This NMR assay appears well-suited for the discovery of selective binding agents in this and other biorecognition phenomena.     

Mutational analysis of the human KDEL receptor: distinct structural requirements for Golgi retention, ligand binding and retrograde transport.

The KDEL receptor is a seven-transmembrane-domain protein that is responsible for the retrieval of endoplasmic reticulum (ER) proteins from the Golgi complex. It is a temporary resident of the Golgi apparatus: upon binding a KDEL-containing ligand, it moves to the ER, where the ligand is released. We have expressed mutant forms of the human receptor in COS cells and examined their intracellular locations and ligand-binding capacities. We show that ligand binding is dependent on charged residues within the transmembrane domains. Surprisingly, retrograde transport of occupied receptor is unaffected by most mutations in the cytoplasmic loops, but is critically dependent upon an aspartic acid residue in the seventh transmembrane domain. Retention in the Golgi apparatus requires neither ligand binding nor this aspartate residue, and thus is independent of receptor recycling. We suggest that movement of the receptor is controlled by conformational changes and intermolecular interactions within the membrane bilayer.     

Validation of an Allosteric Binding Site of Src Kinase Identified by Unbiased Ligand Binding Simulations

Allostery plays a primary role in regulating protein activity, making it an important mechanism in human disease and drug discovery. Identifying allosteric regulatory sites to explore their biological significance and therapeutic potential is invaluable to drug discovery; however, identification remains a challenge. Allosteric sites are often “cryptic” without clear geometric or chemical features. Since allosteric regulatory sites are often less conserved in protein kinases than the orthosteric ATP binding site, allosteric ligands are commonly more specific than ATP competitive inhibitors. We present a generalizable computational protocol to predict allosteric ligand binding sites based on unbiased ligand binding simulation trajectories. We demonstrate the feasibility of this protocol by revisiting our previously published ligand binding simulations using the first identified viral proto-oncogene, Src kinase, as a model system. The binding paths for kinase inhibitor PP1 uncovered three metastable intermediate states before binding the high-affinity ATP-binding pocket, revealing two previously known allosteric sites and one novel site. Herein, we validate the novel site using a combination of virtual screening and experimental assays to identify a V-type allosteric small-molecule inhibitor that targets this novel site with specificity for Src over closely related kinases. This study provides a proof-of-concept for employing unbiased ligand binding simulations to identify cryptic allosteric binding sites and is widely applicable to other protein–ligand systems.     

Identification of fetal liver tyrosine kinase 3 (flt3) ligand domain required for receptor binding and function using naturally occurring ligand isoforms

We used a comparative approach to identify the fetal liver tyrosine kinase 3 (flt3) ligand structure required for binding and function. Two conserved bovine flt3 ligand isoforms, which differ in a defined region within the extracellular domain, were identified and shown to be uniformly transcribed in individuals with diverse MHC haplotypes. Notably, at the amino acid level, the extracellular domain of the bovine flt3 ligand isoform 1 is 81 and 72% identical with the extracellular domains of the human and murine flt3 ligands, respectively, whereas isoform-2 has a deletion within this domain. Bovine flt3 ligand isoform 1, but not 2, bound the human flt3 receptor and stimulated murine pro B cells transfected with the murine flt3 receptor. This retention of binding and function allowed definition of key residues by identifying sequences conserved among species. We have shown that a highly conserved, 18 aa sequence within the flt3 ligand extracellular domain is required for flt3 receptor binding and function. However, a peptide representing this sequence is insufficient for receptor binding as demonstrated by its failure to inhibit the bovine flt3 ligand isoform 1 binding to the human flt3 receptor. The requirement for flanking structure was confirmed by testing bovine flt3 ligand isoform 1 constructs truncated at specific residues outside the 18 aa sequence. Overall, the flt3 ligand structure required for function is markedly similar to that of the related hemopoietic growth factors, CSF-1 and steel factor. This definition of the required flt3 ligand structure will facilitate development of agonists to enhance dendritic cell recruitment for vaccines and immunotherapy.     

The role of molecular size in ligand efficiency

Ligand efficiency is a simple metric for assessing whether a ligand derives its potency from optimal fit with the protein target or simply by virtue of making many contacts. Comparison of protein-ligand binding affinities for over 8000 ligands with 28 protein targets shows conclusively that the average ligand binding affinities are not linear with molecular size. It is therefore important to scale ligand efficiencies by the size of the ligand, particularly where small ligands (e.g., fragments) are involved. We propose a simple 'fit quality' metric that removes this dependence.     

Theoretical analysis of the footprinting experiment

In the footprinting experiment, an end-radiolabeled DNA restriction fragment is subjected to digest by an endonuclease in the presence and absence of a ligand which alters the endonuclease cleavage rate at sites of ligand-DNA contact. The location of these sites, and the strength of the ligand binding, are then deduced from the measured concentrations of the different oligonucleotides produced by the digest. We analyze the experiment in terms of coupled kinetic equations which take into account the cutting rates of endonuclease for sites with ligand present and absent, and the rates of binding and dissociation of the ligand to a site. As long as the ligand concentration remains essentially constant (which occurs, for example, if digest is terminated early enough to assure that all fragments result from single cuts by the endonuclease), the oligonucleotide concentrations reflect only the ligand binding equilibrium constant (ratio of rate constants) and the cutting rates in the presence and absence of ligand. We also show how the measured oligonucleotide concentrations (from, e.g. an autoradiogram) can be used to deduce the ligand equilibrium binding constants for the various sites on the polymer.     

Competitive binding assays for high-affinity binders in the presence of endogenous ligands: application to biotin-binding proteins

Endogenous ligands complicate radioligand-binding assays of high-affinity binding proteins by obscuring binding sites or by diluting the labeled ligand. We have developed a mathematical model for such systems where radioligand and endogenous ligand are structurally identical. Data which relate radioligand binding at equilibrium as a function of sample volume can be plotted such that the concentrations of endogenous ligand and binder are graphically determined; however, a more precise determination may be done by nonlinear regression with the aid of a microcomputer. The method is demonstrated for the assay of biotin-binding proteins in the presence of a range of endogenous biotin concentrations below and above that required to saturate the binding sites.     

Stability of Ligand-induced Protein Conformation Influences Affinity in Maltose-binding Protein

Our understanding of what determines ligand affinity of proteins is poor, even with high-resolution structures available. Both the non-covalent ligand–protein interactions and the relative free energies of available conformations contribute to the affinity of a protein for a ligand. Distant, non-binding site residues can influence the ligand affinity by altering the free energy difference between a ligand-free and ligand-bound conformation. Our hypothesis is that when different ligands induce distinct ligand-bound conformations, it should be possible to tweak their affinities by changing the free energies of the available conformations. We tested this idea for the maltose-binding protein (MBP) from Escherichia coli. We used single-molecule Förster resonance energy transfer (smFRET) to distinguish several unique ligand-bound conformations of MBP. We engineered mutations, distant from the binding site, to affect the stabilities of different ligand-bound conformations. We show that ligand affinity can indeed be altered in a conformation-dependent manner. Our studies provide a framework for the tuning of ligand affinity, apart from modifying binding site residues.     

Stabilization of the FK506 binding protein by ligand binding.

Although the rotamase activity of the FK506 binding protein is inhibited by ligand binding, it is hypothesized that the ligand/protein complex itself may be responsible for the immunosuppressive effects of FK506. We have therefore examined the structure of the FK506 binding protein in the presence of an analog of FK506 (FK520) by a combination of fluorescence, CD, FTIR and calorimetry. While only small changes in the overall structure of the protein may be induced by ligand, a large change in thermal stability of the binding protein is observed.