Calorimetric determination of cooperative interactions in high affinity binding processes

It is demonstrated that isothermal titration calorimetry can be used to determine cooperative interaction energetics even for extremely tight binding processes in which the binding affinity constants are beyond the limits of experimental determination. The approach is based on the capability of calorimetry to measure the apparent binding enthalpy at any degree of ligand saturation. When calorimetric measurements are performed under conditions of total association at partial saturation, the dependence of the apparent binding enthalpy on the degree of saturation is a function only of the cooperative binding interactions. The method developed in this paper allows an independent estimation of cooperative energetic parameters without the need to simultaneously estimate or precisely know the value of the association constants. Since total ligand association at partial saturation is achieved only at macromolecular concentrations much larger than the dissociation constants, the method is especially suited for high and very high affinity processes. Biological associations in this category include fundamental cellular processes like cell surface receptor binding or protein-DNA interactions.     

Elucidating slow binding kinetics of a protein without observable bound resonances by longitudinal relaxation NMR spectroscopy

We developed a new method to elucidate the binding kinetics k_on and k_off, and the dissociation constant K_D (=k_off/k_on), of protein-protein interactions without observable bound resonances of the protein of interest due to high molecular weight in a complex with a large target protein. In our method, k_on and k_off rates are calculated from the analysis of longitudinal relaxation rates of free resonances measured for multiple samples containing different concentration ratios of ¹⁵N-labeled protein and substoichiometric amounts of the target protein. The method is applicable to interactions that cannot be analyzed by relaxation dispersion spectroscopy due to slow interactions on millisecond to second timescale and/or minimal conformational (chemical shift) change upon binding. We applied the method to binding of the B1 domain of protein G (GB1) to immunoglobulin G, and derived the binding kinetics despite the absence of observable bound GB1 resonances.     

Farseer-NMR: automatic treatment,analysis and plotting of large,multi-variable NMR data

Molecular dynamics play a significant role in how molecules perform their function. A critical method that provides information on dynamics, at the atomic level, is NMR-based relaxation dispersion (RD) experiments. RD experiments have been utilized for understanding multiple biological processes occurring at micro-to-millisecond time, such as enzyme catalysis, molecular recognition, ligand binding and protein folding. Here, we applied the recently developed high-power RD concept to the Carr–Purcell–Meiboom–Gill sequence (extreme CPMG; E-CPMG) for the simultaneous detection of fast and slow dynamics. Using a fast folding protein, gpW, we have shown that previously inaccessible kinetics can be accessed with the improved precision and efficiency of the measurement by using this experiment.     

A rapid assay for affinity and kinetics of molecular interactions with nucleic acids

The Differential Radial Capillary Action of Ligand Assay (DRaCALA) allows detection of protein interactions with low-molecular weight ligands based on separation of the protein-ligand complex by differential capillary action. Here, we present an application of DRaCALA to the study of nucleic acid-protein interactions using the Escherichia coli cyclic AMP receptor protein (CRP). CRP bound in DRaCALA specifically to (32)P-labeled oligonucleotides containing the consensus CRP binding site, but not to oligonucleotides with point mutations known to abrogate binding. Affinity and kinetic studies using DRaCALA yielded a dissociation constant and dissociation rate similar to previously reported values. Because DRaCALA is not subject to ligand size restrictions, whole plasmids with a single CRP-binding site were used as probes, yielding similar results. DNA can also function as an easily labeled carrier molecule for a conjugated ligand. Sequestration of biotinylated nucleic acids by streptavidin allowed nucleic acids to take the place of the protein as the immobile binding partner. Therefore, any molecular interactions involving nucleic acids can be tested. We demonstrate this principle utilizing a bacterial riboswitch that binds cyclic-di-guanosine monophosphate. DRaCALA is a flexible and complementary approach to other biochemical methods for rapid and accurate measurements of affinity and kinetics at near-equilibrium conditions.     

Transient Raman study of horseradish peroxidase. Evidence for a lack of extensive heme pocket relaxation subsequent to carbon monoxide photolysis

Time-resolved resonance Raman spectroscopy is a valuable tool for the study of the dynamics of heme-protein interactions. In particular, the photolysis of a ligand by short laser pulses allows for the examination of the dynamic evolution of heme-protein interactions subsequent to ligand dissociation. To date, such studies have been confined largely to hemoglobins and myoglobins. Here we present the results of the first transient Raman study of a peroxidase. Resonance Raman spectra of horseradish peroxidase were obtained within 10 ns of ligand (CO) photolysis at a variety of pH values. We find that there is only minimal relaxation of the heme pocket of horseradish peroxidase in response to ligand photolysis. This relaxation is pH-dependent and most probably involves the heme vinyl substituents. Such behavior is in sharp contrast to the transient behavior of most hemoglobins and beef heart cytochrome oxidase.     

Accuracy of bound peptide structures determined by exchange transferred nuclear Overhauser data: A simulation study

The exchange-transferred NOE method to determine the three-dimensional structure of peptides bound to proteins, or other macromolecular systems, is becoming increasingly important in drug design efforts and for large or multicomponent assemblies, such as membrane receptors, where structural analysis of the full system is intractable. The exchange-transferred nuclear Overhauser effect spectroscopy (etNOESY) method allows the determination of the bound-state conformation of the peptide from the intra-molecular NOE interactions between ligand protons. Because only ligand–ligand NOEs are generally observable, the etNOESY method is restricted to fewer NOEs per residue than direct protein structure determination. In addition, the averaging of relaxation rates between free and bound states affects the measured cross-peak intensities, and possibly the accuracy of distance estimates. Accordingly, the study reported here was conducted to examine the conditions required to define a reliable structure. The program CORONA was used to simulate etNOE data using a rate-matrix including magnetic relaxation and exchange rates for two peptide–protein complexes derived from the reference complex of cAMP-dependent protein kinase ligated with a 24-residue inhibitor peptide. The results indicate that reasonably accurate peptide structures can be determined with relatively few NOE interactions when the interactions occur between non-neighboring residues. The reliability of the structural result is suggested from the pattern of NOE interactions. A structure with an accuracy of approximately 1.3 Å rms difference for the main-chain atoms can be obtained when etNOE interactions between non-neighboring residues occur over the length of the peptide. The global precision is higher (approximately 0.9 Å rms difference) but is not correlated to global accuracy. A local definition of precision along the backbone appears to be a good indicator of the local accuracy.     

Proton magnetic relaxation of aspartate transcarbamylase - succinate complexes.

Nuclear magnetic relaxation methods were used to investigate the interaction of the inhibitor succinate with aspartate transcarbamylase from Escherichia coli. Over the pH range 7 to 9, the dissociation constant for succinate remains less than the inhibitor concentration used for most of this work (0.05 M). As a result, the enzyme predominantly exists in a single "gross" conformational state. Succinate binding to this enzyme state (generally known as the R form) parallels the behavior seen previously with the isolated catalytic subunit (Beard, C. B., and Schmidt, P.G. (1973) Biochemistry 12, 2255-2264). The pH and temperature dependence of succinate proton relaxation rates, 1/T2 - 1/T1, in the presence of carbamyl phosphate, is interpreted in terms of a binding mechanism involving three forms of the enzyme, differing by their states of protonation. The least protonated form of the enzyme does not interact with succinate, the singly protonated species binds succinate to form a rapidly dissociating complex, and the doubly protonated species undergoes a conformational isomerization upon succinate binding, yielding a slow exchange complex. Relaxation data provide sufficient information to determine pKa values of 7.2 and 8.9 for two ionizing groups, as well as the dissociation constant for succinate in the fast exchange complex, Kd =1.6 X 10(-2) M. Rate constants for the forward and reverse steps of the isomerization, 1.3 X 10(3) s-1 and 33 s-1, respectively, indicate a significantly slower reverse rate from that obtained in the earlier NMR study of the isolated catalytic subunit. In experiments where the succinate concentration was varied, the relaxation rates showed sigmoidal binding of that ligand to the fast exchange complex above pH 9.1, (a) indicating cooperative binding of succinate, and (b) suggesting that above pH 9.1, the system cannot be characterized by a single dissociation constant, ionization constant, or relaxation effect. CTP and ATP were tested for their ability to affect succinate binding to the fast exchange complex. Heterotropic interactions were observed for CTP but not for ATP. Addition of low concentrations of the transition state analog N-(phosphonacetyl)-L-aspartate to the enzyme-carbamyl phosphate-succinate complex sharply decreased the relaxation rate, indicating that the measurements are sensitive only to succinate bound specifically to the active site.     

Direct observation of protein-ligand interaction kinetics

Internal dynamics on the micro- to millisecond time scale have a strong influence on the affinity and specificity with which a protein binds ligands. This time scale is accessible through relaxation dispersion measurements using NMR. By studying the dynamics of a protein with different concentrations of a ligand, one can determine the dynamic effects induced by the ligand. Here we have studied slow internal dynamics of the N-terminal src homology 2 domain of phosphatidylinositide 3-kinase to probe the role of individual residues for the interaction with a tyrosine-phosphorylated binding sequence from polyoma middle T antigen. While slow dynamic motion was restricted to a few residues in the free SH2 and in the SH2 complex, motion was significantly enhanced by adding even small amounts of ligand. Kinetic rates induced by ligand binding varied between 300 and 2000 s(-1). High rates reflected direct interactions with the ligand or rearrangements caused by ligand binding. Large differences in rates were observed for residues adjacent in the primary sequence reflecting their individual roles in ligand interaction. However, rates were similar for residues involved in the same side chain interactions, reflecting concerted motions during ligand binding. For a subset of residues, exchange must involve structural intermediates which play a crucial role in high-affinity ligand binding. This analysis supports a new view of the dynamics of individual sites of a protein during ligand interaction.     

Interactions of glucagon and glucagon analogs with isolated canine hepatocytes

We have used glucagon and nine glucagon analogs to investigate the interactions of these ligands with glucagon-binding sites present on isolated canine hepatocytes. Curves reflecting the inhibition of 125I-labeled glucagon or 125I-labeled analog binding to cells by the 10 peptides spanned, overall, a 10(6)-fold range of hormone concentration, were consistent with hormone binding to two classes of binding sites in each case, and fell into two groups, one of which contained curves that were considerably more shallow than the other. Only conditions that emphasized prior binding to low affinity sites resulted in the rapid and extensive dissociation of receptor-bound ligand from isolated cells. Finally, all 10 peptides exhibited a concentration-dependent inhibition of the incorporation of [14C]fructose into hepatocyte glycogen that correlated best with dissociation constants for high affinity rather than for low affinity binding. We conclude that (a) the association of ligand with the high and low affinity glucagon-binding sites of isolated canine hepatocytes is a characteristic of analogs modified at diverse sites throughout the peptide hormone, (b) the different rates of dissociation of ligand from the two populations of binding sites most probably account for the biphasic dissociation of ligand from isolated cells and for the different affinities of the two receptor populations for ligand, and (c) the activity of glucagon and glucagon analogs to inhibit the incorporation of fructose into hepatocyte glycogen arises from the association of ligand with high affinity binding sites.     

Accurate determination of the binding free energy for KcsA-charybdotoxin complex from the potential of mean force calculations with restraints

Free energy calculations for protein-ligand dissociation have been tested and validated for small ligands (50 atoms or less), but there has been a paucity of studies for larger, peptide-size ligands due to computational limitations. Previously we have studied the energetics of dissociation in a potassium channel-charybdotoxin complex by using umbrella sampling molecular-dynamics simulations, and established the need for carefully chosen coordinates and restraints to maintain the physiological ligand conformation. Here we address the ligand integrity problem further by constructing additional potential of mean forces for dissociation of charybdotoxin using restraints. We show that the large discrepancies in binding free energy arising from simulation artifacts can be avoided by using appropriate restraints on the ligand, which enables determination of the binding free energy within the chemical accuracy. We make several suggestions for optimal choices of harmonic potential parameters and restraints to be used in binding studies of large ligands.     

Consequences of ligand bivalency in interactions involving particulate receptors: equilibrium and kinetic studies with Sephadex-concanavalin A, butylagarose-phosphorylase b, and Fc receptor-IgG dimer interactions as model systems

Theoretical consideration is given to the interaction of a bivalent ligand with particulate receptor sites, not only from the viewpoint of quantitatively describing the binding behavior but also from that of the kinetics of ligand release upon infinite dilution of a receptor-ligand mixture. In the latter regard, a general expression is derived that describes the time dependence of the amount of ligand bound as a function of two rate constants for the stepwise dissociation of cross-linked ligand-receptor complex and a thermodynamic parameter expressing the initial ratio of singly linked to doubly linked ligand-receptor complexes. An experimental study of the interaction between Sephadex and concanavalin A is then used to illustrate application of this recommended theoretical approach for characterizing the binding behavior and dissociation kinetics of a bivalent ligand for a system in which all ligand-receptor interactions may be described by a single intrinsic association constant. Published results on the interaction of phosphorylase b with butylagarose are also shown to comply with this simplest model of the bivalent ligand hypothesis; but those for the interaction between immunoglobulin G (IgG) dimers and Fc receptors require modification of the model by incorporation of different intrinsic association constants for the successive binding of receptor sites to the bivalent ligand. These results emphasize the need to consider ligand bivalency as a potential phenomenon in studies of interactions between protein ligands and particulate receptors and illustrate procedures by which the effects of ligand bivalency may be identified and characterized.     

Semiquantitative and quantitative analysis of protein–DNA interactions using steady-state measurements in surface plasmon resonance competition experiments

One method commonly used to characterize protein–DNA interactions is surface plasmon resonance (SPR). In a typical SPR experiment, chip-bound DNA is exposed to increasing concentrations of protein; the resulting binding data are used to calculate a dissociation constant for the interaction. However, in cases in which knowledge of the specificity of the interaction is required, a large set of DNA variants has to be tested; this is time consuming and costly, in part because of the requirement for multiple SPR chips. We have developed a new protocol that uses steady-state binding levels in SPR competition experiments to determine protein-binding dissociation constants for a set of DNA variants. This approach is rapid and straightforward and requires the use of only a single SPR chip. Additionally, in contrast to other methods, our approach does not require prior knowledge of parameters such as on or off rates, using an estimate of the wild-type interaction as the sole input. Utilizing relative steady-state responses, our protocol also allows for the rapid, reliable, and simultaneous determination of protein-binding dissociation constants of a large series of DNA mutants in a single experiment in a semiquantitative fashion. We compare our approach to existing methods, highlighting specific advantages as well as limitations.     

Beta-adrenergic receptors: evidence for negative cooperativity.

The specific binding of [³H] (?)alprenolol to sites in frog erythrocyte membranes provides a tool for directly assessing ligand binding to adenylatecyclase coupled β-adrenergic receptors. Hill Plots of such binding data yield slopes (n_H=“Hill Coefficients”) less than 1.0, suggesting that negatively cooperative interactions among the β-adrenergic receptors may occur. The existence of such negative cooperativity was confirmed by a direct kinetic method. The dissociation of receptor bound [³H] (?)alprenolol was studied under two conditions: 1) with dilution of the ligand-receptor complex sufficient to prevent rebinding of the dissociated tracer and 2) with this same dilution in the presence of excess unlabeled (?)alprenolol. If the sites are independent, the dissociation rates must be the same in both cases. However, the presence of (?)alprenolol increases the rate of [³H] (?)alprenolol dissociation, indicating that negatively cooperative interactions among the β-adrenergic receptor binding sites do occur.     

Quantitative affinity electrophoresis of RNA–small molecule interactions by cross-linking the ligand to acrylamide

We show that the affinity electrophoresis analysis of RNA–small molecule interactions can be made quantifiable by cross-linking the ligand to the gel matrix. Using an RNA–aminoglycoside model system to verify our method, we attached an acryloyl chloride molecule to the aminoglycosides paromomycin and neomycin B to synthesize an acrylamide–aminoglycoside monomer. This molecule was then used as a component in gel polymerization for affinity electrophoresis, covalently attaching an aminoglycoside molecule to the gel matrix. To test RNA binding to the cross-linked aminoglycosides, we used the aminoglycoside binding RNA molecule derived from thymidylate synthase messenger RNA (mRNA) that contains a C–C mismatch. Binding is indicated by the difference in RNA mobility between gels with cross-linked ligand, with ligand embedded during polymerization, and with no ligand present. Critically, the predicted straight line relationship between the reciprocal of the relative migration of the RNA and the ligand concentration is obtained when using cross-linked aminoglycosides, whereas a straight line is not obtained using embedded aminoglycosides. Average apparent dissociation constants are determined from the slope of the line from these plots. This method allows an easy quantitative comparison between different nucleic acid molecules for a small molecule ligand.     

Solid-phase assay for determination of binding parameters of ligand-protein complexes with high dissociation rates

The binding parameters, the affinity constant (Ka) and binding capacity (Q), of a protein possessing ligand-protein complexes with a high dissociation rate (Sex Steroid Binding protein from Bufo arenarum) were determined using a solid-phase method. The technique is based upon the adsorption of the steroid-protein complex to DEAE-cellulose. This method was compared with a nonequilibrium method (charcoal adsorption of free ligand), and the latter resulted in underestimation of both binding parameters, Ka and Q. The solid-phase method reported here is appropriate to determine the binding parameters of proteins with high dissociation rates because the results are independent of the complex half-time. The method also possesses advantages compared to other equilibrium assays such as dialysis or steady-state electrophoresis. With minor modifications, it may be useful to characterize different proteins, particularly those possessing ligand-protein complexes with very high dissociation rates.     

Protein-ligand binding detected using ultrafiltration Raman difference spectroscopy

A new ultra-filtration-Raman-difference (UFRD) method facilitates the tag-free screening and quantitation of protein-ligand binding constants. The method relies on drop-coating-deposition-Raman (DCDR) combined with ultrafiltration and difference spectroscopy. Ultrafiltration is used to remove free (unbound) ligands from pre-equilibrated protein/ligand solutions. Difference DCDR spectroscopy is used to detect binding-induced vibrational spectral changes obtained from proteins with and without a bound ligand. The capabilities of the UFRD method are demonstrated using the binding of 2,4-dinitrophenol (DNP) to transthyretin (TTR), as well as preliminary measurements in several other systems. The UFRD results clearly reveal DNP spectral features induced by binding to TTR and confirm that only a 1:1 complex is formed even under 10-fold excess DNP conditions. The UFRD method is shown to be most useful when applied to strongly Raman active ligands (such as aromatic compounds). Weakly Raman-active ligands (such as sugars) are typically not compatible with UFRD detection (unless they produce a sufficiently large binding-induced change in protein secondary structure). Theoretical predictions suggest that UFRD may be used to screen binding events with a dissociation constant cut-off of the order of 10 microM, and perhaps also to quantify dissociation constants in the 100 nM to 100 microM range.     

Quantitative evaluation of solution equilibrium binding interactions by affinity partitioning: application to specific and nonspecific protein-heparin interactions

A variation of the quantitative affinity chromatography (QAC) method of Winzor, Chaiken, and co-workers for the analysis of protein-ligand interactions has been developed and used to characterize sequence-specific and nonspecific protein-heparin interactions relevant to blood coagulation. The method allows quantitation of the binding of two components, A and B, from the competitive effect of one component, B, on the partitioning of the other component, A, between an immobilized acceptor phase and solution phase at equilibrium. Under the conditions employed, the differences in total A concentrations yielding an equivalent degree of saturation of the immobilized acceptor in the absence and presence of B defines the concentration of A bound to B in solution, thereby enabling conventional Scatchard or nonlinear least-squares analysis of the A-B equilibrium interaction. Like the QAC method, quantitation of the competitor interaction does not depend on the nature of the affinity matrix interaction, which need only be described empirically. The additional advantage of the difference method is that only the total rather than the free competitor ligand concentration need be known. The method requires that the partitioning component A be univalent, but allows for multivalency in the competitor, B, and can in principle be used to study binding interactions involving nonidentical, interacting, or nonspecific overlapping sites. Both the binding constant and the stoichiometry for the specific antithrombin-heparin interaction as well as the apparent binding constant for the nonspecific thrombin-heparin interaction at low thrombin binding densities obtained using this technique were in excellent agreement with values determined using spectroscopic probes.     

Two simple methods for quantifying low-affinity dye-substrate binding

Binding with low-affinity ligands, such as histological dyes, can be difficult to quantitate owing to the dissociation of bound ligand with washing or the retention of nonspecifically bound ligand because of incomplete washing. The present report describes two simple, rapid methods of discriminating bound from free ligand without the need for washing steps. One method is based on the spectral changes induced in a dye ligand, Congo red, on binding to the "receptor" insulin fibrils. This method discriminates spectrophotometrically between bound and free ligand without requiring any physical separation of the two forms. No radioactive ligands are necessary, and, by using disposable cuvettes, the entire binding assay can be done in a single container without the need for transfers. The second method employs a non-traditional filtration approach that avoids the need for a washing step by measuring the decrease in concentration of the dye ligand in the filtrate rather than by applying the usual approach of measuring the absolute amount of ligand bound to the precipitated "receptor." Both methods show saturation of binding sites and give similar values for the KD and Bmax.     

Mechanisms of ligand binding by monoclonal anti-fluorescyl antibodies

Binding of fluorescyl ligand by five IgG anti-fluorescyl hybridoma proteins (4-4-20, 6-10-6, 20-4-4, 20-19-=1, 20-20-3) was examined. Relative reduction in fluorescence of bound fluorescein, deuterium oxide (D2O)-induced enhancement of fluorescence, and the effects of pH on binding kinetics were measured for each clone. Individual hybridoma proteins (all of which bind fluorescein with relatively high affinity) exhibited significant differences in the relative contribution of various forces (hydrophobicity, hydrogen bonding, ionic interactions) to binding and hence, affinity. The extent of such variations in binding mechanisms among monoclonal antibodies binding the same hapten is indicative of the extreme functional diversity of active sites. In addition, ligand binding by clone 20-20-3 was examined in greater detail. ABsorption spectra of ligand bound by purified intact antibody, Fab fragments, and reassociated heavy and light chains indicated that protonation of the fluorescyl ligand by a residue within the active site contributed significantly to the binding free energy. Comparative dissociation rates of fluorescein and a structural analog, rhodamine 110, were used to quantitatively substantiate the contribution of this interaction. Association and dissociation rate studies with fluorescein and antibody indicated that: 1) the active site appeared to undergo a conformational change upon ligand binding, and 2) neither intact disulfides nor intersite cooperativity affected the dissociation rate of bound ligand. Observed mechanisms of ligand binding are discussed in terms of proposed mechanisms of antibody affinity maturation and diversity.