Energy landscape of chelated uranyl: antibody interactions by dynamic force spectroscopy

We used dynamic force spectroscopy (DFS) to explore the energy landscape of interactions between a chelated uranyl compound and a monoclonal antibody raised against the uranyl-dicarboxy-phenanthroline complex. We estimated the potential energy barrier widths and the relevant thermodynamic rate constants along the dissociation coordinate. Using atomic force microscopy, four different experimental setups with or without the uranyl ion in the chelate ligand, we have distinguished specific and nonspecific binding in the binding affinity of the uranyl compound to the antibody. The force loading rates for our system were measured from 15 to 26,400 pN/s. The results showed two regimes in the plot of the most probable unbinding force versus the logarithm of the loading rate, revealing the presence of two (at least) activation barriers. Analyses of DFS suggest parallel multivalent binding present in either regime. We have also built a molecular model for the variable fragment of the antibody and used computational graphics to dock the chelated uranyl ion into the binding pocket. The structural analysis led us to hypothesize that the two regimes originate from two interaction modes: the first one corresponds to an energy barrier with a very narrow width of 0.5 +/- 0.2 A, inferring dissociation of the uranyl ion from its first coordination shell (Asp residue); the second one with a broader energy barrier width (3.9 +/- 0.3 A) infers the entire chelate compound dissociated from the antibody. Our study highlights the sensitivity of DFS experiments to dissect protein-metal compound interactions.     

Cooperative adhesion of ligand-receptor bonds

Cooperative (simultaneous) breakage of multiple adhesive bonds has been proposed as a mechanism for enhanced binding strength between adhesion molecules on apposing cell surfaces. In this report, we used the atomic force microscopy (AFM) to study how changes in binding affinity and separation rate of force-induced ligand-receptor dissociation affect binding cooperativity. The AFM force measurements were carried out using (strept)avidin-functionalized cantilever tips and biotinylated agarose beads under conditions where multiple (strept)avidin-biotin linkages were formed following surface contact. At slow surface separation of the AFM cantilever from the bead's surface, the (strept)avidin-biotin linkages appeared to rupture sequentially. Increasing the separation rate from 210 to 1950 nm/s led to a linear increase in the average rupture force. Moreover, force histograms revealed a quantized force distribution that shifted toward higher values with increasing separation rate. In measurements of streptavidin-iminobiotin adhesion, the force distribution also shifted toward higher values when the buffer was adjusted to a higher pH to raise the binding affinity. Together, these results demonstrate that the cooperativity of ligand-receptor bonds is significantly enhanced by increases in surface separation rate and/or binding affinity.     

Force spectroscopy of barnase-barstar single molecule interaction

Results of the single molecule force spectroscopy study of specific interactions between ribonuclease barnase and its inhibitor barstar are presented. Experimental data obtained for the force loading rate ranging 2-70 nN/s are well approximated by a single straight line, from which the dissociation barrier of the width of 0.12 nm and height of 0.75-0.85 × 10(-19)J can be inferred. The measured value of specific interaction does not depend on the NaCl concentration. This apparently contradicts the well-known dependence of the binding energy of this pair on the salt concentration, but such a "contradiction" is explained by the insensitivity of the force spectroscopy data to the relatively long-range electrostatic interaction. The latter essentially contributes to the value of barnase-barstar binding energy revealed by biochemical measurements, and it is exactly this electrostatic interaction which is influenced by the salt concentration.     

Stable peptide binding to MHC class II molecule is rapid and is determined by a receptive conformation shaped by prior association with low affinity peptides

Formation of stable class II MHC/peptide complex involves conformational changes and proceeds via an intermediate. Although this intermediate complex forms and dissociates in minutes, its conversion to a stable complex is a very slow process, taking up to a few days to reach completion. Here, we investigate the different steps of this binding and demonstrate that the conformational changes necessary to generate a receptive molecule is the rate-determining slow step in the process, while formation of the stable MHC/peptide complex is very rapid. With HLA-DR1 as our model class II molecule, we first used low affinity variants of hemagglutinin peptide (HA306-318), which lack the principal anchor, to shape the conformation of the MHC and then studied the kinetics of stable binding of HA306-318 to such an induced conformation. We found that the apparent association rate of HA306-318 is equivalent to the dissociation rate of the low affinity peptide. A 4- to 18-fold enhancement in the binding rates of HA306-318 was observed depending on the dissociation rates of the low affinity peptides. These results establish that 1) formation of stable MHC/peptide complexes is very rapid and 2) prior binding of low affinity peptide induces a receptive conformation in MHC for efficient stable peptide binding. Furthermore, in the absence of any free peptide, this receptive molecule rapidly reverts to slow binding behavior toward the subsequently offered peptide. These results have important implications for the roles of low affinity MHC/peptide complexes in Ag presentation.     

Metal ion binding and coordination geometry for wild type and mutants of metallo-beta -lactamase from Bacillus cereus 569/H/9 (BcII): a combined thermodynamic, kinetic, and spectroscopic approach

One high affinity (nm) and one low affinity (microM) macroscopic dissociation constant for the binding of metal ions were found for the wild-type metallo-beta-lactamase from Bacillus cereus as well as six single-site mutants in which all ligands in the two metal binding sites were altered. Surprisingly, the mutations did not cause a specific alteration of the affinity of metal ions for the sole modified binding site as determined by extended x-ray absorption fine structure (EXAFS) and perturbed angular correlation of gamma-rays spectroscopy, respectively. Also UV-visible absorption spectra for the mono-cobalt enzymes clearly contain contributions from both metal sites. The observations of the very similar microscopic dissociation constants of both binding sites in contrast to the significantly differing macroscopic dissociation constants inevitably led to the conclusion that binding to the two metal sites exhibits negative cooperativity. The slow association rates for forming the binuclear enzyme determined by stopped-flow fluorescence measurements suggested that fast metal exchange between the two sites for the mononuclear enzyme hinders the binding of a second metal ion. EXAFS spectroscopy of the mono- and di-zinc wild type enzymes and two di-zinc mutants provide a definition of the metal ion environments, which is compared with the available x-ray crystallographic data.     

The analysis of peptide affinity and its binding kinetics to DR1DW1 major histocompatibility complex protein.

The connection between experimentally measured values of ED50 (concentration of added peptide required to bind half of the protein), which characterize peptide-protein binding and the equilibrium dissociation constant of peptide-protein complex Kd (affinity) is considered. It is shown and confirmed by experimental studies that in certain cases, as a result of the absence of equilibrium in the system, the value of Kd could be much less than the experimental value of ED50, but not equal to that as commonly assumed. This is especially applicable to the formation of peptide-MHC complexes with low dissociation rates (strong binding), which may require longer time-intervals to reach equilibrium. Thus the search of the good binding peptides based on finding ones with the smallest measured values' of ED50 may result in missing the best binders with the lowest values of dissociation constant (highest affinity). To analyze the problem we considered the formal chemical kinetics of peptide-protein binding. Experimental studies of peptide binding was performed to obtain the parameters of the kinetic model. According to the predictions of the model, it was confirmed that peptide binding occurs through the preceding step, which is either a release of an endogenous peptide or some conformational change of the molecule. The half decay time for this process was determined to be approximately 3 h. Based on the model developed, a new effective method for determination of the dissociation rates of peptide-MHC complexes and the equilibrium dissociation constants Kd was proposed, which implies the comparison of binding levels (ED50) at different instants of time. This method works especially well for the peptide-MHC complexes with relatively slow dissociation rates (stable complexes), for which the direct off-rate measurements as well as obtaining equilibrium binding data to determine Kd are highly time consuming and not very reliable.     

Muscarinic binding to mouse brain receptor sites

In mouse brain the binding of [³H]-Atropine to the muscarinic receptor seems to be a simple mass-action determined process as gauged both by approach to equilibrium kinetics and binding at equilibrium. In contrast, using isotopic dilution technique, dissociation measurements indicate the existence of two receptor-ligand complexes. It would appear that association and dissociation rates of binding of the muscarinic antagonists atropine, scopolamine, N-methyl-4-piperidyl benzilate (4NMPB) and 3-quinuclidinyl benzilate (QNB) decrease with increasing affinity based on comparisons of kinetic binding data. The differences between the association rate constants are small whereas those between the dissociation rate constants differ markedly. This kinetic behavior is similar to the well-known time profile of antimuscarinic activity in isolated tissues. These phenomena are discussed in terms of possible isomerization of the receptor-ligand complex, as has been proposed recently for [³H]-scopolamine and [³H]-4NMPB binding.     

Investigating the binding behaviour of two avidin‐based testosterone binders using molecular recognition force spectroscopy

Molecular recognition force spectroscopy, a biosensing atomic force microscopy technique allows to characterise the dissociation of ligand–receptor complexes at the molecular level. Here, we used molecular recognition force spectroscopy to study the binding capability of recently developed testosterone binders. The two avidin‐based proteins called sbAvd‐1 and sbAvd‐2 are expected to bind both testosterone and biotin but differ in their binding behaviour towards these ligands. To explore the ligand binding and dissociation energy landscape of these proteins, we tethered biotin or testosterone to the atomic force microscopy probe while the testosterone‐binding protein was immobilized on the surface. Repeated formation and rupture of the ligand–receptor complex at different pulling velocities allowed determination of the loading rate dependence of the complex‐rupturing force. In this way, we obtained the molecular dissociation rate (k_off) and energy landscape distances (x_β) of the four possible complexes: sbAvd‐1‐biotin, sbAvd‐1‐testosterone, sbAvd‐2‐biotin and sbAvd‐2‐testosterone. It was found that the kinetic off‐rates for both proteins and both ligands are similar. In contrast, the x_β values, as well as the probability of complex formations, varied considerably. In addition, competitive binding experiments with biotin and testosterone in solution differ significantly for the two testosterone‐binding proteins, implying a decreased cross‐reactivity of sbAvd‐2. Unravelling the binding behaviour of the investigated testosterone‐binding proteins is expected to improve their usability for possible sensing applications. Copyright © 2014 John Wiley & Sons, Ltd.     

Affinity profiling using the peptide microarray technology: a case study

Little about the reliability of measurements obtained using synthetic peptide microarrays is known. We report results from a study on the quantitative reliability of microarrays manufactured by robot-supported immobilization of presynthesized peptides for different microarray platforms. Technological precision is assessed for inter- and intra-device readout comparisons. Correlations between measured signals and known dissociation constants using a phenomenological model derived from the mass action law are discussed. Special emphasis is on discussing the pitfalls of high-throughput affinity measurements. We show that the quantitative determination of binding affinities is prone to be biased toward a mean affinity of around 10(-7)M, while the classification of peptides into either "binders" or "nonbinders" provides very high prediction accuracy. The experimental requirements needed to obtain reliable binding affinity predictions are discussed.     

Energy landscape roughness of the streptavidin-biotin interaction

Molecular interactions between receptors and ligands can be characterized by their free energy landscape. In its simplest representation, the energy landscape is described by a barrier of certain height and width that determines the dissociation rate of the complex, as well as its dynamic strength. Some interactions, however, require a more complex landscape with additional barriers and roughness along the reaction coordinate. This roughness slows down the dissociation kinetics of the interaction and contributes to its dynamic strength. The streptavidin-biotin complex has been extensively studied due to its remarkably low dissociation kinetics. However, single molecule measurements from independent experiments showed scattered and disparate results. In this work, the energy landscape roughness of the streptavidin-biotin interaction was estimated to be in the range of 5-8kBT using dynamic force spectroscopy (DFS) measurements at three different temperatures. These results can be used to explain both its slow dissociation kinetics and the discrepancies in the reported force measurements.     

Site-selective agonist binding to the nicotinic acetylcholine receptor from Torpedo californica

Fluorescent energy transfer measurements of dansyl-C6-choline binding to the nicotinic acetylcholine receptor (AChR) from Torpedo californica were used to determine binding characteristics of the alpha gamma and alpha delta binding sites. Equilibrium binding measurements show that the alpha gamma site has a lower fluorescence than the alpha delta site; the emission difference is due to differences in the intrinsic fluorescence of the bound fluorophores rather than differences in energy transfer at the two sites. Stopped-flow fluorescence kinetics showed that dissociation of dansyl-C6-choline from the AChR in the desensitized conformation occurs 5-10-fold faster from the alpha gamma site than from the alpha delta site. The dissociation rates are robust for distinct protein preparations, in the presence of noncompetitive antagonists, and over a broad range of ionic strengths. Equilibrium fluorescent binding measurements show that dansyl-C6-choline binds with higher affinity to the alpha delta site (K = 3 nM) than to the alpha gamma site (K = 9 nM) when the AChR is desensitized. Similar affinity differences were observed for acetylcholine itself. The distinct dissociation rates permit the extent of desensitization to be measured at each site during the time course of binding. This sequential mixing method of measuring the desensitized state population at each agonist site can be applied to study the mechanism of AChR activation and subsequent desensitization in detail.     

3(17)beta-Hydroxysteroid Dehydrogenase of Pseudomonas testosteroni. Ligand binding properties.

The binding of NAD and NADH to electrophoretically pure 3(17)beta-hydroxysteroid dehydrogenase of Pseudomonas testosteroni was determined by Fluorescence spectroscopy and gel filtration. Four moles of cofactor are bound/mol of tetrameric enzyme; the binding sites are equivalent and independent. The dissociation constants for NAD and NADH are 16 and 0.25 micronM, respectively. As measured by gel filtration in the absence of cofactor, 0.4 mol of estradiol-17 beta is bound/mol of tetrameric enzyme. Data obtained from isotope exchange at equilibrium indicate that the binding of the cofactor to the enzyme is favored over the binding of steroid, although each may bind in the absence of the other. The rates of cofactor dissociation from the ternary complexes are slower than the rates of steroid dissociation; cofactor dissociation is probably the rate-limiting step. Cofactor analogs modified in the pyridine moiety are cosubstrates, whereas modified adenine derivatives are not. The enzyme also utilized as substrate a number of potential steroid affinity labels; no enzyme inactivation by these compounds was observed.     

From affinity selection to kinetic selection in Germinal Centre modelling

Affinity maturation is an evolutionary process by which the affinity of antibodies (Abs) against specific antigens (Ags) increases through rounds of B-cell proliferation, somatic hypermutation, and positive selection in germinal centres (GC). The positive selection of B cells depends on affinity, but the underlying mechanisms of affinity discrimination and affinity-based selection are not well understood. It has been suggested that selection in GC depends on both rapid binding of B-cell receptors (BcRs) to Ags which is kinetically favourable and tight binding of BcRs to Ags, which is thermodynamically favourable; however, it has not been shown whether a selection bias for kinetic properties is present in the GC. To investigate the GC selection bias towards rapid and tight binding, we developed an agent-based model of GC and compared the evolution of founder B cells with initially identical low affinities but with different association/dissociation rates for Ag presented by follicular dendritic cells in three Ag collection mechanisms. We compared an Ag collection mechanism based on association/dissociation rates of B-cell interaction with presented Ag, which includes a probabilistic rupture of bonds between the B-cell and Ag (Scenario-1) with a reference scenario based on an affinity-based Ag collection mechanism (Scenario-0). Simulations showed that the mechanism of Ag collection affects the GC dynamics and the GC outputs concerning fast/slow (un)binding of B cells to FDC-presented Ags. In particular, clones with lower dissociation rates outcompete clones with higher association rates in Scenario-1, while remaining B cells from clones with higher association rates reach higher affinities. Accordingly, plasma cell and memory B cell populations were biased towards B-cell clones with lower dissociation rates. Without such probabilistic ruptures during the Ag extraction process (Scenario-2), the selective advantage for clones with very low dissociation rates diminished, and the affinity maturation level of all clones decreased to the reference level.     

Diversification of specificity after maturation of the antibody response to the HIV gp41 epitope ELDKWA

During maturing antibody responses the increase in affinity for target antigens is achieved by genetic diversification of antibody genes followed by selection for improved binding. The effect this process has on the specificity of antibody for variants of the antigen is not well-defined, despite the potential role of antibody diversification in generating enhanced protection against pathogen escape mutants, or novel specificities after vaccination. To investigate this, a library of single amino-acid substitution epitope variants has been screened with serum obtained at different time-points after immunization of mice with the HIV gp41 peptide epitope ELDKWA. The serum IgG response is shown to mature and increase affinity for ELDKWA, and the titre and affinity of IgG against most epitope variants tested increases. Furthermore there is a bias towards high affinity serum IgG binding to variant epitopes with conservative substitutions, although underlying this trend there is also significant binding to many epitopes with non-conservative substitutions. Thus, maturation of the antibody response to a single epitope results in a broadening of the high-affinity response toward variant epitopes. This implies that many pathogen epitope escape variants that could manifest as single amino-acid substitutions would not emerge by escaping immune surveillance.     

Shotgun alanine scanning shows that growth hormone can bind productively to its receptor through a drastically minimized interface

The high affinity binding site (Site1) of the human growth hormone (hGH) binds to its cognate receptor (hGHR) via a concave surface patch containing about 35 residues. Using 167 sequences from a shotgun alanine scanning analysis of Site1, we have determined that over half of these residues can be simultaneously changed to an alanine or a non-isosteric amino acid while still retaining a high affinity interaction. Among these hGH variants the distribution of the mutation is highly variable throughout the interface, although helix 4 is more conserved than the other binding elements. Kinetic and thermodynamic analyses were performed on 11 representative hGH Site1 variants that contained 14-20 mutations. Generally, the tightest binding variants showed similar associated rate constants (k(on)) as the wild-type (wt) hormone, indicating that their binding proceeds through a similar transition state intermediate. However, calorimetric analyses indicate very different thermodynamic partitioning: wt-hGH binding exhibits favorable enthalpy and entropy contributions, whereas the variants display highly favorable enthalpy and highly unfavorable entropy contributions. The heat capacities (DeltaCp) on binding measured for wt-hGH and its variants are significantly larger than normally seen for typical protein-protein interactions, suggesting large conformational or solvation effects. The multiple Site1 mutations are shown to indirectly affect binding of the second receptor at Site2 through an allosteric mechanism. We show that the stability of the ternary hormone-receptor complex reflects the affinity of the Site2 binding and is surprisingly exempt from changes in Site1 affinity, directly demonstrating that dissociation of the active signaling complex is a stepwise process.     

Kinetics of calcium dissociation from its high-affinity transport sites on sarcoplasmic reticulum ATPase.

We investigated the kinetics of calcium dissociation from its high-affinity transport sites on sarcoplasmic reticulum Ca2(+)-ATPase by combining fast filtration with stopped-flow fluorescence measurements. At pH 6 and 20 degrees C, in the absence of potassium and in the presence of 20 mM MgCl2, isotopic exchange of bound calcium exhibited biphasic kinetics, with two phases of equal amplitude, regardless of the initial extent of binding site saturation. The rapidly exchangeable site, whose occupancy by calcium controlled the rate constant of the slow phase, had an apparent affinity for calcium of about 3-6 microM. A similar high affinity was also deduced from measurements of the calcium dependence of the rate constant for ATPase fluorescence changes. This affinity was higher than the overall affinity for calcium deduced from the equilibrium binding measurements (dissociation constant of 15-20 microM); this was consistent with the occurrence of cooperativity (Hill coefficient of 1.6-1.8). The drop in intrinsic fluorescence observed upon chelation of calcium was always slightly faster than the dissociation of calcium itself, although the rates for both this drop in fluorescence and calcium dissociation varied slightly from one preparation to the other. This fluorescence drop was therefore mainly due to dissociation of the bound ions, not to slow transconformation of the ATPase. Dissociation of the two bound calcium ions in a medium containing EGTA exhibited monophasic kinetics in the presence of a calcium ionophore, with a rate constant about half that of the fast phase of isotopic exchange. This particular pattern was observed over a wide range of experimental conditions, including the presence of KCl, dimethyl sulfoxide, 4-nonylphenol, or a nucleotide analogue, at pH 6 or 7, and at various temperatures. The kinetics of calcium dissociation under the above various conditions were not correlated with the ATPase affinity for calcium deduced from equilibrium measurements under the same conditions. These results are consistent with sequential dissociation of calcium from a narrow binding pocket inside which a single calcium ion can move fairly easily. Escape of calcium might be controlled by a structural compartment acting as a gate.     

Revisiting the phosphotyrosine binding pocket of Fyn SH2 domain led to the identification of novel SH2 superbinders.

Protein engineering through directed evolution is an effective way to obtain proteins with novel functions with the potential applications as tools for diagnosis or therapeutics. Many natural proteins have undergone directed evolution in vitro in the test tubes in the laboratories worldwide, resulting in the numerous protein variants with novel or enhanced functions. we constructed here an SH2 variant library by randomizing 8 variable residues in its phosphotyrosine (pTyr) binding pocket. Selection of this library by a pTyr peptide led to the identification of SH2 variants with enhanced affinities measured by EC50. Fluorescent polarization was then applied to quantify the binding affinities of the newly identified SH2 variants. As a result, three SH2 variants, named V3, V13 and V24, have comparable binding affinities with the previously identified SH2 triple‐mutant superbinder. Biolayer Interferometry assay was employed to disclose the kinetics of the binding of these SH2 superbinders to the phosphotyrosine peptide. The results indicated that all the SH2 superbinders have two‐orders increase of the dissociation rate when binding the pTyr peptide while there was no significant change in their associate rates. Intriguingly, though binding the pTyr peptide with comparable affinity with other SH2 superbinders, the V3 does not bind to the sTyr peptide. However, variant V13 and V24 have cross‐reactivity with both pTyr and sTyr peptides. The newly identified superbinders could be utilized as tools for the identification of pTyr‐containing proteins from tissues under different physiological or pathophysiological conditions and may have the potential in the therapeutics.     

Direct and Model Free Calculation of Force-Dependent Dissociation Rates from Force Spectroscopic Data

Force spectroscopy allows testing the free energy landscapes of molecular interactions. Usually, the dependency of the most probable rupture force on the force rate or the shape of the rupture force histogram is fitted with different models that contain approximations and basic assumptions. We present a simple and model free approach to extract the force-dependent dissociation rates directly from the force curve data. Simulations show that the dissociation rates at any force are given directly by the ratio of the number of detected rupture events to the time this force was acting on the bond. To calculate these total times of acting forces, all force curve data points of all curves measured are taken into account, which significantly increases the amount of information which is considered for data analysis compared to other methods. Moreover, by providing force-dependent dissociation rates this method allows direct testing and validating of any energy landscape model.     

Monitored whole gene in vitro evolution of an anti-hRaf-1 affibody molecule towards increased binding affinity

The use of library technologies for the generation of affinity proteins often includes an affinity maturation step, based on the construction of secondary libraries from which second generation variants with improved affinities are selected. Here, we describe for the first time the affinity maturation of affibody molecules based on step-wise in vitro molecular evolution, involving cycles of error-prone PCR (epPCR) amplification for the introduction of diversity over the entire 58-residue three-helix bundle structure and ribosome display (RD) for the selection of improved variants. The model affibody molecule for the process was Z(RAF322), binding with a 1.9μm equilibrium dissociation constant (K(D)) to human Raf-1 (hRaf-1), a protein kinase of central importance in the MAPK/ERK proliferation pathway. The molecular evolution process was followed on both gene and protein levels via DNA sequencing and a biosensor-based binding analysis of pools of selected variants. After two cycles of diversification and selection, a significant increase in binding response of selected pools was seen. DNA sequencing showed that a dominant alanine to valine substitution had been effectively enriched, and was found in 83% of all selected clones, either alone or in combination with other enriched substitutions. The evolution procedure resulted in variants showing up to 26-fold increases in affinity to the hRaf-1 target. Noteworthy, for the two variants showing the highest affinities, substitutions were also found in affibody framework positions, corresponding to regions of the protein domain not addressed by traditional affibody molecule affinity maturation strategies. Interestingly, thermal melting point (T(m)) analyses showed that an increased affinity could be associated with both higher and lower T(m) values. All investigated variants showed excellent refolding properties and selective binding to hRaf-1, as analysed using a multiplexed bead-based binding assay, making them potentially valuable affinity reagents for cell biology studies.