Evaluation of temperature effect on the interaction between β-lactoglobulin and anti-β-lactoglobulin antibody by atomic force microscopy

Molecular recognition such as antigen-antibody interaction is characterized by the parameters of kinetics and the energy landscape. Examinations of molecules involved in the interaction at different temperatures using atomic force microscopy (AFM) can provide information on not only the effects of temperature on the unbinding force between a molecule of interest and a complementary molecule but also the parameters of kinetics and the energy landscape for dissociation of the molecular complex. We investigated the effect of temperature on the dissociation process of the complex of β-lactoglobulin and anti-bovine β-lactoglobulin IgG polyclonal antibody using AFM. Measurements of the unbinding forces between β-lactoglobulin and the antibody were performed at 25, 35, and 45 °C. The following results were obtained in our present study: (i) The unbinding forces decreased as temperature increased, suggesting that the binding force between β-lactoglobulin and the antibody includes the force originating from temperature-dependent interactions (e.g., hydrogen bonding). (ii) At each temperature, the unbinding force exhibited two linear regimes in the force spectra, indicating that the dissociation process of the β-lactoglobulin-antibody complex passes at least two energy barriers from the bound state to the dissociated state. (iii) The dissociation rates at zero force and the position of energy barriers increased as temperature increased. (iv) The heights of the two energy barriers in the reaction coordinates were 49.7 k(B)T and 14.5 k(B)T. (v) The values of roughness of the barriers were ca. 6.1 k(B)T and 3.2 k(B)T. Overall, the present study using AFM revealed more information about the β-lactoglobulin-antibody interaction than studies using conventional bulk measurement such as surface plasmon resonance.     

Energy landscape of streptavidin-biotin complexes measured by atomic force microscopy

The dissociation of ligand and receptor involves multiple transitions between intermediate states formed during the unbinding process. In this paper, we explored the energy landscape of the streptavidin-biotin interaction by using the atomic force microscope (AFM) to measure the unbinding dynamics of individual ligand-receptor complexes. The rupture force of the streptavidin-biotin bond increased more than 2-fold over a range of loading rates between 100 and 5000 pN/s. Moreover, the force measurements showed two regimes of loading in the streptavidin-biotin force spectrum, revealing the presence of two activation barriers in the unbinding process. Parallel experiments carried out with a streptavidin mutant (W120F) were used to investigate the molecular determinants of the activation barriers. From these experiments, we attributed the outer activation barrier in the energy landscape to the molecular interaction of the '3-4' loop of streptavidin that closes behind biotin.     

Dynamic force measurements of avidin-biotin and streptavdin-biotin interactions using AFM

Using atomic force microscopy (AFM) we performed dynamic force measurements of the adhesive forces in two model systems: avidin-biotin and streptavidin-biotin. In our experiments we used glutaraldehyde for immobilization of (strept)avidin on the tip and biotin on the sample surface. Such interface layers are more rigid than those usually reported in the literature for AFM studies, when (strept)avidin is coupled with biotinylated bovine albumin and biotin with agarose polymers. We determined the dependence of the rupture forces of avidin-biotin and streptavidin-biotin bonds in the range 300-9600 pN/s. The slope of a semilogarithmic plot of this relation changes at about 1700 pN/s. The existence of two different regimes indicates the presence of two activation barriers of these complexes during the dissociation process. The dissociation rates and activation energy barriers, calculated from the Bell model, for the avidin-biotin and streptavidin-biotin interactions are similar to each other for loading rates > 1700 pN/s but they are different from each other for loading rates < 1700 pN/s. In the latter case, the dissociation rates show a higher stability of the avidin-biotin complex than the streptavidin-biotin complex due to a larger outer activation barrier of 0.8 k(B)T. The bond-rupture force is about 20 pN higher for the avidin-biotin pair than for the streptavidin-biotin pair for loading rates < 1700 pN/s. These two experimental observations are in agreement with the known structural differences between the biotin binding pocket of avidin and of streptavidin.     

Dynamic strength of molecular adhesion bonds.

In biology, molecular linkages at, within, and beneath cell interfaces arise mainly from weak noncovalent interactions. These bonds will fail under any level of pulling force if held for sufficient time. Thus, when tested with ultrasensitive force probes, we expect cohesive material strength and strength of adhesion at interfaces to be time- and loading rate-dependent properties. To examine what can be learned from measurements of bond strength, we have extended Kramers' theory for reaction kinetics in liquids to bond dissociation under force and tested the predictions by smart Monte Carlo (Brownian dynamics) simulations of bond rupture. By definition, bond strength is the force that produces the most frequent failure in repeated tests of breakage, i.e., the peak in the distribution of rupture forces. As verified by the simulations, theory shows that bond strength progresses through three dynamic regimes of loading rate. First, bond strength emerges at a critical rate of loading (> or = 0) at which spontaneous dissociation is just frequent enough to keep the distribution peak at zero force. In the slow-loading regime immediately above the critical rate, strength grows as a weak power of loading rate and reflects initial coupling of force to the bonding potential. At higher rates, there is crossover to a fast regime in which strength continues to increase as the logarithm of the loading rate over many decades independent of the type of attraction. Finally, at ultrafast loading rates approaching the domain of molecular dynamics simulations, the bonding potential is quickly overwhelmed by the rapidly increasing force, so that only naked frictional drag on the structure remains to retard separation. Hence, to expose the energy landscape that governs bond strength, molecular adhesion forces must be examined over an enormous span of time scales. However, a significant gap exists between the time domain of force measurements in the laboratory and the extremely fast scale of molecular motions. Using results from a simulation of biotin-avidin bonds (Izrailev, S., S. Stepaniants, M. Balsera, Y. Oono, and K. Schulten. 1997. Molecular dynamics study of unbinding of the avidin-biotin complex. Biophys. J., this issue), we describe how Brownian dynamics can help bridge the gap between molecular dynamics and probe tests.     

Direct measurement of protein energy landscape roughness

The energy landscape of proteins is thought to have an intricate, corrugated structure. Such roughness should have important consequences on the folding and binding kinetics of proteins, as well as on their equilibrium fluctuations. So far, no direct measurement of protein energy landscape roughness has been made. Here, we combined a recent theory with single-molecule dynamic force spectroscopy experiments to extract the overall energy scale of roughness epsilon for a complex consisting of the small GTPase Ran and the nuclear transport receptor importin-beta. The results gave epsilon > 5k(B)T, indicating a bumpy energy surface, which is consistent with the ability of importin-beta to accommodate multiple conformations and to interact with different, structurally distinct ligands.     

Force-based analysis of multidimensional energy landscapes: application of dynamic force spectroscopy and steered molecular dynamics simulations to an antibody fragment-peptide complex

Multidimensional energy landscapes are an intrinsic property of proteins and define their dynamic behavior as well as their response to external stimuli. In order to explore the energy landscape and its implications on the dynamic function of proteins dynamic force spectroscopy and steered molecular dynamics (SMD) simulations have proved to be important tools. In this study, these techniques have been employed to analyze the influence of the direction of the probing forces on the complex of an antibody fragment with its peptide antigen. Using an atomic force microscope, experiments were performed where the attachment points of the 12 amino acid long peptide antigen were varied. These measurements yielded clearly distinguishable basal dissociation rates and potential widths, proving that the direction of the applied force determines the unbinding pathway. Complementary atomistic SMD simulations were performed, which also show that the unbinding pathways of the system are dependent on the pulling direction. However, the main barrier to be crossed was independent of the pulling direction and is represented by a backbone hydrogen bond between Gly^H-H40 of the antibody fragment and Glu^Oε-6_peptide of the peptide. For each pulling direction, the observed barriers can be correlated with the rupture of specific interactions, which stabilize the bound complex. Furthermore, although the SMD simulations were performed at loading rates exceeding the experimental rates by orders of magnitude due to computational limitations, a detailed comparison of the barriers that were overcome in the SMD simulations with the data obtained from the atomic force microscope unbinding experiments show excellent agreement.     

Reaction pathway and free energy profile determined for specific recognition of oligosaccharide moiety of carboxypeptidase Y

The interaction of mannose specific lectin (from Lens culinaris, LcL) with the carbohydrate moiety of carboxypeptidase Y (CaY) was studied using both atomic force microscope (AFM) and quartz crystal microbalance with dissipation monitoring (QCM-D). The AFM enables to determine the positions of energy barriers present in the energy landscape of the single complex undergoing dissociation. The QCM-D measurements allow the estimation of the quantitative parameters characterizing the kinetics of the studied molecular interaction (namely the association and dissociation rate constants and the association constant). The use of both methods not only delivers the complementary characterization of kinetic and thermodynamic parameters but also permits to investigate the mechanism of the binding and unbinding of the molecules. The results for LcL were compared with those obtained for concanavalin A i.e. lectin, which interacts with the carbohydrate moiety on a similar way.     

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.     

Force spectroscopy of LFA-1 and its ligands, ICAM-1 and ICAM-2

Single-molecule measurements of the interaction of leukocyte function-associated antigen-1 (LFA-1), expressed on Jurkat T cells, with intercellular adhesion molecules-1 and -2 (ICAM-1 and ICAM-2) were conducted using atomic force microscopy (AFM). The force spectra (i.e., unbinding force versus loading rate) of both the LFA-1/ICAM-1 and LFA-1/ICAM-2 interactions were acquired at a loading rate range covering 3 orders of magnitude (50-60,000 pN/s) and revealed a fast loading regime and a slow loading regime. This indicates that the dissociation of both complexes involves overcoming a steep inner and a wide outer activation barrier. LFA-1 binding to ICAM-1 and ICAM-2 was strengthened in the slow loading regime by the addition of Mg(2+). Differences in the dynamic strength of the LFA-1/ICAM-1 and LFA-1/ICAM-2 interactions can be attributed to the presence of wider barriers in the ICAM-2 complex, making it more responsive to a pulling force than the ICAM-1 complex.     

LFA-1 binding destabilizes the JAM-A homophilic interaction during leukocyte transmigration

Leukocyte transendothelial migration into inflamed areas is regulated by the integrity of endothelial cell junctions and is stabilized by adhesion molecules including junctional adhesion molecule-A (JAM-A). JAM-A has been shown to participate in homophilic interactions with itself and in heterophilic interactions with leukocyte function-associated antigen-1 (LFA-1) via its first and second immunoglobulin domains, respectively. Using competitive binding assays in conjunction with atomic force microscopy adhesion measurements, we provide compelling evidence that the second domain of JAM-A stabilizes the homophilic interaction because its deletion suppresses the dynamic strength of the JAM-A homophilic interaction. Moreover, binding of the LFA-1 inserted domain to the second domain of JAM-A reduces the dynamic strength of the JAM-A homophilic interaction to the level measured with the JAM-A domain 2 deletion mutant. This finding suggests that LFA-1 binding cancels the stabilizing effects of the second immunoglobulin domain of JAM-A. Finally, our atomic force microscopy measurements reveal that the interaction of JAM-A with LFA-1 is stronger than the JAM-A homophilic interaction. Taken together, these results suggest that LFA-1 binding to JAM-A destabilizes the JAM-A homophilic interaction. In turn, the greater strength of the LFA-1/JAM-A complex permits it to support the tension needed to disrupt the JAM-A homophilic interaction, thus allowing transendothelial migration to proceed.     

Force spectroscopy of the leukocyte function-associated antigen-1/intercellular adhesion molecule-1 interaction

Interactions between leukocyte function-associated antigen-1 (LFA-1) with its cognate ligand, intercellular adhesion molecule-1 (ICAM-1) play a crucial role in leukocyte adhesion. Because the cell and its adhesive components are subject to external perturbation from the surrounding flow of blood, it is important to understand the binding properties of the LFA-1/ICAM-1 interaction in both steady state and in the presence of an external pulling force. Here we report on atomic force microscopy (AFM) measurements of the unbinding of LFA-1 from ICAM-1. The single molecule measurements revealed the energy landscape corresponding to the dissociation of the LFA-1/ICAM-1 complex and provided the basis for defining the energetic determinants of the complex at equilibrium and under the influence of an external force. The AFM force measurements were performed in an experimental system consisting of an LFA-1-expressing T cell hybridoma, 3A9, attached to the end of the AFM cantilever and an apposing surface expressing ICAM-1. In measurements covering three orders of magnitude change in force loading rate, the LFA-1/ICAM-1 force spectrum (i.e., unbinding force versus loading rate) revealed a fast and a slow loading regime that characterized a steep inner activation barrier and a wide outer activation barrier, respectively. The addition of Mg(2+), a cofactor that stabilizes the LFA-1/ICAM-1 interaction, elevated the unbinding force of the complex in the slow loading regime. In contrast, the presence of EDTA suppressed the inner barrier of the LFA-1/ICAM-1 complex. These results suggest that the equilibrium dissociation constant of the LFA-1/ICAM-1 interaction is regulated by the energetics of the outer activation barrier of the complex, while the ability of the complex to resist a pulling force is determined by the divalent cation-dependent inner activation barrier.     

The relationship between ligand-binding thermodynamics and protein-ligand interaction forces measured by atomic force microscopy.

The interaction forces between biotin and a set of streptavidin site-directed mutants with altered biotin-binding equilibrium and activation thermodynamics have been measured by atomic force microscopy. The AFM technique readily discriminates differences in interaction force between the site-directed (Trp to Phe or Ala) mutants. The interaction force is poorly correlated with both the equilibrium free energy of biotin binding and the activation free energy barrier to dissociation of the biotin-streptavidin complex. The interaction force is generally well correlated with the equilibrium biotin-binding enthalpy as well as the enthalpic activation barrier, but in the one mutant where these two parameters are altered in opposite directions, the interaction force is clearly correlated with the activation enthalpy of dissociation. These results suggest that the AFM force measurements directly probe the enthalpic activation barrier to ligand dissociation.     

Two-dimensional kinetics of inter-connexin interactions from single-molecule force spectroscopy

Gap junction channels are intercellular channels that form by docking the extracellular loops of connexin protein subunits. While the structure and function of gap junctions as intercellular channels have been characterized using different techniques, the physics of the inter-connexin interaction remain unknown. Moreover, as far as we know, the capacity of gap junction channels to work as adhesion complexes supporting pulling forces has not yet been quantitatively addressed. We report the first quantitative characterization of the kinetics and binding strength of the interaction of a short peptide mimicking extracellular loop 2 of Cx26 with membrane-reconstituted Cx26, combining the imaging and force spectroscopy capabilities of atomic force microscopy. The fast dissociation rate inferred a dynamic bond, while the slow association rate reflected the reduced flexibility and small size of extracellular loops. Our results propose the gap junction channel as an adhesion complex that associates slowly and dissociates fast at low force but is able to support important pulling forces in its native, hexameric form.     

Extending Bell's model: how force transducer stiffness alters measured unbinding forces and kinetics of molecular complexes

Forced unbinding of complementary macromolecules such as ligand-receptor complexes can reveal energetic and kinetic details governing physiological processes ranging from cellular adhesion to drug metabolism. Although molecular-level experiments have enabled sampling of individual ligand-receptor complex dissociation events, disparities in measured unbinding force F_R among these methods lead to marked variation in inferred binding energetics and kinetics at equilibrium. These discrepancies are documented for even the ubiquitous ligand-receptor pair, biotin-streptavidin. We investigated these disparities and examined atomic-level unbinding trajectories via steered molecular dynamics simulations, as well as via molecular force spectroscopy experiments on biotin-streptavidin. In addition to the well-known loading rate dependence of F_R predicted by Bell's model, we find that experimentally accessible parameters such as the effective stiffness of the force transducer k can significantly perturb the energy landscape and the apparent unbinding force of the complex for sufficiently stiff force transducers. Additionally, at least 20% variation in unbinding force can be attributed to minute differences in initial atomic positions among energetically and structurally comparable complexes. For force transducers typical of molecular force spectroscopy experiments and atomistic simulations, this energy barrier perturbation results in extrapolated energetic and kinetic parameters of the complex that depend strongly on k. We present a model that explicitly includes the effect of k on apparent unbinding force of the ligand-receptor complex, and demonstrate that this correction enables prediction of unbinding distances and dissociation rates that are decoupled from the stiffness of actual or simulated molecular linkers.     

Direct force measurements between siRNA and chitosan molecules using force spectroscopy

Information on the interaction strength between small interfering RNA (siRNA) and chitosan can contribute to the understanding of the formation and stability of chitosan/siRNA nanoparticles used as siRNA delivery systems for gene silencing. In this study, we utilize atomic force microscopy to obtain force spectroscopy results of the interaction strengths between siRNA and chitosan measured in physiological phosphate buffered saline buffer at different pH. The force measurements revealed that the adhesive interactions decreased in force strength and force frequency as the pH was increased from 4.1 to 6.1, 7.4, and 9.5, exhibiting distinct multimodal distributions of the interaction forces between siRNA and chitosan molecules at acidic pH and only negligible adhesive forces were observed at neutral or high pH. The strong pH dependence of siRNA-chitosan interactions can provide a convincing rationale for siRNA/chitosan complex formation and nanoparticle stability under low acidic conditions. These findings demonstrate that the use of force spectroscopy for the adhesive force measurements allows an evaluation of the complexing ability between siRNA and chitosan that can be utilized to predict nanoparticle stability.     

Probing effects of pH change on dynamic response of Claudin-2 mediated adhesion using single molecule force spectroscopy

Claudins belong to a large family of transmembrane proteins that localize at tight junctions (TJs) where they play a central role in regulating paracellular transport of solutes and nutrients across epithelial monolayers. Their ability to regulate the paracellular pathway is highly influenced by changes in extracellular pH. However, the effect of changes in pH on the strength and kinetics of claudin mediated adhesion is poorly understood. Using atomic force microscopy, we characterized the kinetic properties of homophilic trans-interactions between full length recombinant GST tagged Claudin-2 (Cldn2) under different pH conditions. In measurements covering three orders of magnitude change in force loading rate of 10²–10⁴ pN/s, the Cldn2/Cldn2 force spectrum (i.e., unbinding force versus loading rate) revealed a fast and a slow loading regime that characterized a steep inner activation barrier and a wide outer activation barrier throughout pH range of 4.5–8. Comparing to the neutral condition (pH 6.9), differences in the inner energy barriers for the dissociation of Cldn2/Cldn2 mediated interactions at acidic and alkaline environments were found to be < 0.65 k_BT, which is much lower than the outer dissociation energy barrier (> 1.37 k_BT). The relatively stable interaction of Cldn2/Cldn2 in neutral environment suggests that electrostatic interactions may contribute to the overall adhesion strength of Cldn2 interactions. Our results provide an insight into the changes in the inter-molecular forces and adhesion kinetics of Cldn2 mediated interactions in acidic, neutral and alkaline environments.     

Model energy landscapes and the force-induced dissociation of ligand-receptor bonds

We discuss models for the force-induced dissociation of a ligand-receptor bond, occurring in the context of cell adhesion or single molecule unbinding force measurements. We consider a bond with a structured energy landscape which is modeled by a network of force dependent transition rates between intermediate states. The behavior of a model with only one intermediate state and a model describing a molecular zipper is studied. We calculate the bond lifetime as a function of an applied force and unbinding forces under an increasing applied load and determine the relationship between both quantities. The dissociation via an intermediate state can lead to distinct functional relations of the bond lifetime on force. One possibility is the occurrence of three force regimes where the lifetime of the bond is determined by different transitions within the energy landscape. This case can be related to recent experimental observations of the force-induced dissociation of single avidin-biotin bonds.     

Temperature Dependence of Unbinding Forces between Complementary DNA Strands

Force probe techniques such as atomic force microscopy can directly measure the force required to rupture single biological ligand receptor bonds. Such forces are related to the energy landscape of these weak, noncovalent biological interactions. We report unbinding force measurements between complementary strands of DNA as a function of temperature. Our measurements emphasize the entropic contributions to the energy landscape of the bond.     

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.     

Dissecting streptavidin-biotin interaction with a laminar flow chamber

A laminar flow chamber was used to study single molecule interactions between biotinylated surfaces and streptavidin-coated spheres subjected to a hydrodynamic drag lower than a piconewton. Spheres were tracked with 20 ms and 40 nm resolution. They displayed multiple arrests lasting between a few tens of milliseconds and several minutes or more. Analysis of about 500,000 positions revealed that streptavidin-biotin interaction was multiphasic: transient bound states displayed a rupture frequency of 5.3 s(-1) and a rate of transition toward a more stable configuration of 1.3 s(-1). These parameters did not display any significant change when the force exerted on bonds varied between 3.5 and 11 pN. However, the apparent rate of streptavidin-biotin association exhibited about 10-fold decrease when the wall shear rate was increased from 7 to 22 s(-1), which supports the existence of an energy barrier opposing the formation of the transient binding state. It is concluded that a laminar flow chamber can yield new and useful information on the formation of molecular bonds, and especially on the structure of the external part of the energy landscape of ligand-receptor complexes.