Quantitative considerations of the consequences of an interplay between ligand binding and reversible adsorption of a macromolecular solute

Explicit expressions are derived which determine the equilibrium composition of mixtures comprising a multivalent, insoluble matrix, a multivalent, macromolecular solute (acceptor) and a univalent ligand. With three-reactant mixtures of this type a range of combinations of interactions is possible wherein the ligand interacts with either the acceptor or the matrix, in either event perturbing the acceptor-matrix equilibria. Theory encompassing this range of possibilities is written in terms of a single site-binding constant for each type of interaction to account, in general terms, for both multiple binding and crosslinking effects. These explicit thermodynamic relationships are discussed, with the use of reported findings on several biological systems, in two frameworks. First, it is established that the theory is applicable to the quantitative interpretation of affinity chromatography experiments designed to elucidate the thermodynamic interaction parameters governing the various types of interacting system. Second, it is emphasized that the relationships are also relevant to metabolite-induced changes in the subcellular distribution of macromolecular species.     

Interpretation of Scatchard plots for aggregating receptor systems.

Aggregation of cell surface receptors, with each other or with other membrane proteins, occurs in a variety of experimental systems. The list of systems where receptor aggregation appears to be important in understanding ligand binding and cellular responses is growing rapidly. In this paper we explore the interpretation of equilibrium binding data for aggregating receptor systems. The Scatchard plot is a widely used tool for analyzing equilibrium binding data. The shape of the Scatchard plot is often interpreted in terms of multiple noninteracting receptor populations. Such an analysis does not provide a framework for investigating the role of receptor aggregation and will be misleading if there is a relation between receptor aggregation and ligand binding. We present a general model for the equilibrium binding of a ligand with any number of aggregating receptor populations and derive theoretical expressions for observable Scatchard plot features. These can be used to test particular models and estimate model parameters. We develop particular models and apply the general results in the cases of six aggregating receptor systems where ligand binding and receptor aggregation are related: cross-linking of monovalent cell surface proteins by monoclonal antibodies, cross-linking of cell surface antibodies by bivalent ligand, antibody-induced co-cross-linking of cell surface antibodies and Fc gamma receptors, ligand-enhanced aggregation of identical epidermal growth factor receptors, aggregation of heterologous receptors for interleukin 2 to form a high-affinity receptor, and association of receptors, including those for interleukins 5 and 6, with nonbinding accessory proteins that influence receptor affinity or effector function.     

Evaluation of equilibrium constants by affinity chromatography

Theoretical expressions are derived for affinity chromatography of systems comprising an acceptor A with one binding site for attachment to a functional group X on the column matrix and one site for interaction with a small ligand B that specifically affects its elution. From a general relationship covering all possible interactions between A, B and X simpler expressions are derived for affinity systems in which only two equilibria operate. Methods are suggested whereby these simpler systems may be characterized in terms of the two pertinent equilibrium constants and the concentration of matrix-bound constituent. The means by which the theory may be adapted to affinity chromatography of acceptors with multiple binding sites for ligand is also illustrated. Results of partition experiments on the Sephadex G-100-lysozyme-d-glucose system in acetate-chloride buffer (I=0.17m), pH5.4, are used to demonstrate the feasibility of evaluating quantitatively affinity-chromatography interactions. Values of 30m(-1) and 1.2x10(6)m(-1) are obtained for the equilibrium constants for the reactions of lysozyme with glucose and Sephadex respectively, there being only an occasional binding site in the polysaccharide matrix (approximately 1 in 10(5) glucose residues). In a second experimental study the phytohaemagglutinin from Ricinus communis is subjected to frontal chromatography on Sepharose 4B in the presence of different concentrations of d-galactose, the results illustrating some of the difficulties and limitations that are likely to be encountered in quantitative studies of affinity-chromatographic systems.     

Stoichiometry of interaction between interferon gamma and its receptor.

The biological response of interferon gamma is mediated by binding to a specific cell-surface receptor. We investigated the stoichiometry of this binding using soluble receptors produced in prokaryotic and eukaryotic expression systems comprising the extracellular ligand-binding domain of the native protein. The ligand-receptor complexes were analyzed by cross-linking, chromatography, analytical ultracentrifugation and laser-light scattering. Cross-linking and chromatography showed that the stoichiometry of the interaction between ligand and receptor depends on the molar ratios of the two components mixed. All approaches confirmed that mixtures of ligand-receptor complexes are formed with one interferon-gamma dimer bound by one or two receptors. The soluble receptor produced in Escherichia coli mainly showed a ligand/receptor stoichiometry of 1:1, while the receptors produced in eukaryotic cells showed a stoichiometry of binding of 1:2. This apparent discrepancy is most likely due to the conformational heterogeneity of the Escherichia-coli-derived protein.     

Binding equations for interacting systems comprising multivalent acceptor and bivalent ligand: application to antigen-antibody systems

Consideration is given to the reversible interaction of a bivalent ligand, B, with a multivalent acceptor, A (possessing f reactive sites) which leads to the formation of a series of complexes, A_iB_j, comprising networks of alternating acceptor and ligand molecules. A binding equation is derived on the basis of a site association constant, k, defined in terms of reacted site probability functions. This equation, which relates the binding function, r (the moles of ligand bound per mole of acceptor) to the concentration of unbound ligand, m_b, is used to show that plots of r vs. 2km_B constructed with fixed but different values of $\text{k}\text{m}{}_{\mathrm{A}}$ intersect at the point ($\text{m}{}_{\mathrm{B}}\text{=}\text{1}\text{2k}\text{, r =}\text{f}\text{2}$) where the extent of reaction and the concentrations of those complexes for which $\text{j}\text{i}\text{=}\text{f}\text{2}$ attain maximal values. Corresponding Scatchard plots are shown by numerical example to be non-linear, their second derivative being positive for all r. It follows that such deviations from linearity cannot be taken alone as evidence for site heterogeneity in cross-linking systems. The binding equation obtained directly is shown to be identical with that obtained with f = 2 by summation procedures involving the general expression for concentrations of complexes, m_AiBj, formulated in terms of appropriate statistical factors. In this way, previous findings on precipitation and gel formation in cross-linking systems are correlated with the present development of binding theory.     

Multiscale Simulations Examining Glycan Shield Effects on Drug Binding to Influenza Neuraminidase

Influenza neuraminidase is an important drug target. Glycans are present on neuraminidase and are generally considered to inhibit antibody binding via their glycan shield. In this work, we studied the effect of glycans on the binding kinetics of antiviral drugs to the influenza neuraminidase. We created all-atom in silico systems of influenza neuraminidase with experimentally derived glycoprofiles consisting of four systems with different glycan conformations and one system without glycans. Using Brownian dynamics simulations, we observe a two- to eightfold decrease in the rate of ligand binding to the primary binding site of neuraminidase due to the presence of glycans. These glycans are capable of covering much of the surface area of neuraminidase, and the ligand binding inhibition is derived from glycans sterically occluding the primary binding site on a neighboring monomer. Our work also indicates that drugs preferentially bind to the primary binding site (i.e., the active site) over the secondary binding site, and we propose a binding mechanism illustrating this. These results help illuminate the complex interplay between glycans and ligand binding on the influenza membrane protein neuraminidase.     

Molecular insights into the binding selectivity of a synthetic ligand DAAG to Fc fragment of IgG

Affinity chromatography with synthetic ligands has been focused as the potential alternative to protein A‐based chromatography for antibody capture because of its comparable selectivity and efficiency. Better understanding on the molecular interactions between synthetic ligand and antibody is crucial for improving and designing novel ligands. In this work, the molecular interaction mechanism between Fc fragment of IgG and a synthetic ligand (DAAG) was studied with molecular docking and dynamics simulation. The docking results on the consensus binding site (CBS) indicated that DAAG could bind to the CBS with the favorable orientation like a tripod for the top‐ranked binding complexes. The ligand‐Fc fragment complexes were then tested by molecular dynamics simulation at neutral condition (pH 7.0) for 10 ns. The results indicated that the binding of DAAG on the CBS of Fc fragment was achieved by the multimodal interactions, combining the hydrophobic interaction, electrostatic interaction, hydrogen bond, and so on. It was also found that multiple secondary interactions endowed DAAG with an excellent selectivity to Fc fragment. In addition, molecular dynamics simulation conducted at acidic condition (pH 3.0) showed that the departure of DAAG ligand from the surface of Fc fragment was the result of reduced interaction energies. The binding modes between DAAG and CBS not only shed light on the molecular mechanisms of DAAG for antibody purification but also provide useful information for the improvement of ligand design. Copyright © 2014 John Wiley & Sons, Ltd.     

A novel approach to local similarity of protein binding sites substantially improves computational drug design results

We present a novel notion of binding site local similarity based on the analysis of complete protein environments of ligand fragments. Comparison of a query protein binding site (target) against the 3D structure of another protein (analog) in complex with a ligand enables ligand fragments from the analog complex to be transferred to positions in the target site, so that the complete protein environments of the fragment and its image are similar. The revealed environments are similarity regions and the fragments transferred to the target site are considered as binding patterns. The set of such binding patterns derived from a database of analog complexes forms a cloud-like structure (fragment cloud), which is a powerful tool for computational drug design. It has been shown on independent test sets that the combined use of a traditional energy-based score together with the cloud-based score responsible for the quality of embedding of a ligand into the fragment cloud improves the self-docking and screening results dramatically. The usage of a fragment cloud as a source of positioned molecular fragments fitting the binding protein environment has been validated by reproduction of experimental ligand optimization results.     

Cooperation effects in binding of large ligands to DNA. II. Contact interactions between adsorbed ligands

Cooperative effects arising upon binding of biologically active ligands to DNA are considered. Equations are derived which enable one to describe the binding of two different ligands to DNA. We also consider the case when ligand can form two type of DNA complexes. The cooperative binding of the ligand in the vicinity of saturation level of binding can be described with a good accuracy by equation derived for the non-cooperative adsorption of the same ligand with some effective binding constant Keff. It is shown that cooperative effects arising upon binding of proteins and other ligands to DNA can be divided into two groups depending on the symmetry of interactions between the bound ligand molecules. In particular, if such interactions favor the formation of dimeric ligand species on the DNA, Keff approximately a1/2, where a is the ligand-ligand interaction constant. If cooperative interactions favor the formation of aggregates of unrestricted size, then Keff approximately aL+Y, where L is the size of the binding site for the ligand on DNA.     

Molecular docking study and development of an empirical binding free energy model for phosphodiesterase 4 inhibitors

In the present work, several computational methodologies were combined to develop a model for the prediction of PDE4B inhibitors' activity. The adequacy of applying the ligand docking approach, keeping the enzyme rigid, to the study of a series of PDE4 inhibitors was confirmed by a previous molecular dynamics analysis of the complete enzyme. An exhaustive docking procedure was performed to identify the most probable binding modes of the ligands to the enzyme, including the active site metal ions and the surrounding structural water molecules. The enzyme-inhibitor interaction enthalpies, refined by using the semiempirical molecular orbital approach, were combined with calculated solvation free energies and entropy considerations in an empirical free energy model that enabled the calculation of binding free energies that correlated very well with experimentally derived binding free energies. Our results indicate that both the inclusion of the structural water molecules close to the ions in the binding site and the use of a free energy model with a quadratic dependency on the ligand free energy of solvation are important aspects to be considered for molecular docking investigations involving the PDE4 enzyme family.     

A consensus-binding structure for adenine at the atomic level permits searching for the ligand site in a wide spectrum of adenine-containing complexes

Attempts to derive structural features of ligand-binding sites have traditionally involved seeking commonalities at the residue level. Recently, structural studies have turned to atomic interactions of small molecular fragments to extract common binding-site properties. Here, we explore the use of larger ligand elements to derive a consensus binding structure for the ligand as a whole. We superimposed multiple molecular structures from a nonredundant set of adenosine-5'-triphosphate (ATP) protein complexes, using the adenine moiety as template. Clustered binding-site atoms of compatible atomic classes forming attractive contacts with the adenine probe were extracted. A set of atomic clusters characterizing the adenine binding pocket was then derived. Among the clusters are three vertices representing the interactions of adenine atom N6 with its protein-binding niche. These vertices, together with atom C6 of the purine ring system, complete the set of four vertices for the pyramid-like structure of the N6 anchor atom. Also, the sequence relationship for the adenine-binding loop interacting with the C2-N6 end of the conjugated ring system is expanded to include a third hydrophilic cluster interacting with atom N1. A search procedure involving interatomic distances between cluster centers was formulated and applied to seek putative binding sites in test cases. The results show that a consensus network of clusters, based on an adenine probe and an ATP-complexed training set of proteins, is sufficient to recognize the experimental cavity for adenine in a wide spectrum of ligand-protein complexes.     

Cross-linking reconsidered: binding and cross-linking fields and the cellular response.

We analyze a model for the reversible cross-linking of cell surface receptors by a collection of bivalent ligands with different affinities for the receptor as would be found in a polyclonal anti-receptor serum. We assume that the amount of cross-linking determines, via a monotonic function, the rate at which cells become activated and divide. In addition to the density of receptors on the cell surface, two quantities, the binding field and the cross-linking field, are needed to characterize the cross-linking curve, i.e., the equilibrium concentration of cross-linked receptors plotted as a function of the total ligand site concentration. The binding field is the sum of all ligand site concentrations weighted by their respective binding affinities, and the cross-linking field is the sum of all ligand site concentrations weighted by the product of their respective binding and cross-linking affinity and the total receptor density. Assuming that the cross-linking affinity decreases if the binding affinity decreases, we find that the height of the cross-linking curve decreases, its width narrows, and its center shifts to higher ligand site concentrations as the affinities decrease. Moreover, when we consider cross-linking-induced proliferation, we find that there is a minimum cross-linking affinity that must be surpassed before a clone can expand. We also show that under many circumstances a polyclonal antiserum would be more likely than a monoclonal antibody to lead to cross-linking-induced proliferation.     

Interdomain interactions in the tumor suppressor discs large regulate binding to the synaptic protein GukHolder

Multidomain scaffolding proteins are central components of many signaling pathways and are commonly found at membrane specializations. Here we have shown that multiple interdomain interactions in the scaffold Discs Large (Dlg) regulate binding to the synaptic protein GukHolder (GukH). GukH binds the Src homology 3 (SH3) and guanylate kinase-like (GK) protein interaction domains of Dlg, whereas an intramolecular interaction between the two domains inhibits association with GukH. Regulation occurs through a PDZ domain adjacent to the SH3 that allows GukH to interact with the composite SH3-GK binding site, but PDZ ligands inhibit GukH binding such that Dlg forms mutually exclusive PDZ ligand and GukH cellular complexes. The PDZ-SH3-GK module is a common feature of membrane associate guanylate kinase scaffolds such as Dlg, and these results indicate that its supramodular architecture leads to regulation of Dlg complexes.     

EQUIL: Simulation and data analysis of binding reactions with arbitrary chemical models

We have developed an algorithm for simulation and analysis of arbitrary chemical systems in equilibrium, with emphasis on ligand binding reactions. The program EQUIL can treat reactions involving multiple ligands, multiple binding sites, ternary complex models, allosteric effectors, competitive and noncompetitive binding, conformational changes, cooperativity, and generally any scheme that can be represented as a set of chemical equations. EQUIL is based on a general thermodynamic model of chemical equilibria; it does not involve nonlinear transformation of experimental data, but it does require the user to define the model of interaction between ligands and receptors by writing down the appropriate chemical reactions. EQUIL contains features of particular importance to ligand binding experiments: variable binding capacities, nonspecific binding, and the ability to simultaneously analyze data from different types of experiments. Furthermore, the simulation feature of EQUIL allows the user to investigate the feasibility of experiments that could possibly distinguish between different reaction models. We illustrate the use of this program on personal computers to analyze and simulate simple and complicated interactions between ligands and receptors.     

Merging the binding sites of aldose and aldehyde reductase for detection of inhibitor selectivity-determining features

Inhibition of human aldose reductase (ALR2) evolved as a promising therapeutic concept to prevent late complications of diabetes. As well as appropriate affinity and bioavailability, putative inhibitors should possess a high level of selectivity for ALR2 over the related aldehyde reductase (ALR1). We investigated the selectivity-determining features by gradually mapping the residues deviating between the binding pockets of ALR1 and ALR2 into the ALR2 binding pocket. The resulting mutational constructs of ALR2 (eight point mutations and one double mutant) were probed for their influence towards ligand selectivity by X-ray structure analysis of the corresponding complexes and isothermal titration calorimetry (ITC). The binding properties of these mutants were evaluated using a ligand set of zopolrestat, a related uracil derivative, IDD388, IDD393, sorbinil, fidarestat and tolrestat. Our study revealed induced-fit adaptations within the mutated binding site as an essential prerequisite for ligand accommodation related to the selectivity discrimination of the ligands. However, our study also highlights the limits of the present understanding of protein-ligand interactions. Interestingly, binding site mutations not involved in any direct interaction to the ligands in various cases show significant effects towards their binding thermodynamics. Furthermore, our results suggest the binding site residues deviating between ALR1 and ALR2 influence ligand affinity in a complex interplay, presumably involving changes of dynamic properties and differences of the solvation/desolvation balance upon ligand binding.     

The ligand binding site of the platelet integrin receptor GPIIb-IIIa is proximal to the second calcium binding domain of its alpha subunit

The extreme carboxyl-terminal amino acid sequence of the gamma chain of fibrinogen is involved in the binding of this adhesive protein to the platelet integrin glycoprotein (GP) IIb-IIIa, and synthetic peptides corresponding to this region inhibit fibrinogen as well as fibronectin and von Willebrand factor binding to platelets. A chemical cross-linking approach was used to characterize the interaction of a 16-amino acid fibrinogen gamma chain peptide with platelets and to localize the site of its binding to GPIIb-IIIa. This peptide became specifically cross-linked to GPIIb, and platelet stimulation selectively enhanced its cross-linking to this alpha subunit. The cross-linking reaction was specifically inhibited by fibrinogen and an Arg-Gly-Asp peptide but not by an unrelated protein or a substituted peptide. Utilizing a combination of immunochemical mapping, enzymatic and chemical digestions, and amino acid sequencing, the cross-linking site of the gamma chain peptide in GPIIb was localized to a stretch of 21 amino acids. The identified region, GPIIb 294-314, contains the second putative calcium binding domain within GPIIb. The primary structure of this region is highly conserved among alpha subunits of other integrin adhesion receptors. These results identify a discrete region of GPIIb that resides in close proximity to a ligand binding site within GPIIb-IIIa. The homologous region may be involved in the functions of other integrin receptors.