Force and Stress along Simulated Dissociation Pathways of Cucurbituril-Guest Systems

The field of host-guest chemistry provides computationally tractable yet informative model systems for biomolecular recognition. We applied molecular dynamics simulations to study the forces and mechanical stresses associated with forced dissociation of aqueous cucurbituril-guest complexes with high binding affinities. First, the unbinding transitions were modeled with constant velocity pulling (steered dynamics) and a soft spring constant, to model atomic force microscopy (AFM) experiments. The computed length-force profiles yield rupture forces in good agreement with available measurements. We also used steered dynamics with high spring constants to generate paths characterized by a tight control over the specified pulling distance; these paths were then equilibrated via umbrella sampling simulations and used to compute time-averaged mechanical stresses along the dissociation pathways. The stress calculations proved to be informative regarding the key interactions determining the length-force profiles and rupture forces. In particular, the unbinding transition of one complex is found to be a stepwise process, which is initially dominated by electrostatic interactions between the guest's ammoniums and the host's carbonyl groups, and subsequently limited by the extraction of the guest's bulky bicyclooctane moiety; the latter step requires some bond stretching at the cucurbituril's extraction portal. Conversely, the dissociation of a second complex with a more slender guest is mainly driven by successive electrostatic interactions between the different guest's ammoniums and the host's carbonyl groups. The calculations also provide information on the origins of thermodynamic irreversibilities in these forced dissociation processes.     

Dynamic force spectroscopy of the specific interaction between the PDZ domain and its recognition peptides

To characterize the molecular basis of specific interactions of PDZ proteins, dynamic force spectroscopy (DFS) for the PDZ protein Tax-interacting protein-1 (TIP-1) and its recognition peptide (PDZ-pep) derived from beta-catenin was performed using an atomic force microscope (AFM), together with measurement of thermodynamic and kinetic parameters using surface plasmon resonance (SPR). The unbinding force of this pair was measured under different conditions of AFM tip-retraction velocity. The relationship between the unbinding force and the logarithmic force-loading rate, that is, the dynamic force spectrum, exhibited two different rate regimes, for each of which the forces increased linearly with the force-loading rate. On the basis of the theoretical treatment of the Bell-Evans model, the positions of two different activation barriers in the reaction coordinate and dissociation rate constants in each barrier were evaluated from slopes and x-intercepts of the two linear regimes (first barrier: 0.04 nm and 1.10 x 10 s(-1); second barrier: 0.21 nm and 2.77 x 10(-2) s(-1), respectively). Although two-step unbinding kinetics between TIP-1 and PDZ-pep was suggested from the DFS analysis, SPR results showed single-step dissociation kinetics with a rate constant of 2.89 x 10(-1) s(-1). Different shapes of the free energy profile of the unbinding process were deduced from each result of DFS and SPR. The reason for such topographic differences in the energy landscape is discussed in relation to the differences in the pathways of forced unbinding and spontaneous dissociation.     

Comparison of antibody--antigen interactions on collagen measured by conventional immunological techniques and atomic force microscopy

We have developed a means of using atomic force microscopy (AFM) to repeatedly localize a small area of interest (4 x 4 microm(2)) within a 0.5-cm(2) area on a heterogeneous sample, to obtain and localize high-resolution images and force measurements on nonideal samples (i.e., samples that better reflect actual biological systems, not prepared on atomically flat surfaces). We demonstrate the repeated localization and measurement of unbinding forces associated with antibody--antigen (ab--ag) interactions, by applying AFM in air and in liquid to visualize and measure polyclonal ab--ag interactions, using chicken collagen as a model system. We demonstrate that molecular interactions, in the form of ab--ag complexes, can be visualized by AFM when secondary antibodies are conjugated to 20-nm colloidal gold particles. We then compare those results with established immunological techniques, to demonstrate broader application of AFM technology to other systems. Data from AFM studies are compared with results obtained using immunological methods traditionally employed to investigate ab--ag interactions, including enzyme-linked immunosorbent assay, immunoblotting, and in situ immunofluorescence. Finally, using functionalized AFM tips with a flexible tether [poly(ethylene glycol) 800] to which a derivatized antibody was attached, we analyzed force curve data to measure the unbinding force of collagen antibody from its antigen, obtaining a value of approximately 90 +/- 40 pN with a MatLab code written to automate the analyses of force curves obtained in force--volume mode. The methodology we developed for embedded collagen sections can be readily applied to the investigation of other receptor--ligand interactions.     

Antigen binding forces of single antilysozyme Fv fragments explored by atomic force microscopy

We used atomic force microscopy (AFM) to explore the antigen binding forces of individual Fv fragments of antilysozyme antibodies (Fv). To detect single molecular recognition events, genetically engineered histidine-tagged Fv fragments were coupled onto AFM tips modified with mixed self-assembled monolayers (SAMs) of nitrilotriacetic acid- and tri(ethylene glycol)-terminated alkanethiols while lysozyme (Lyso) was covalently immobilized onto mixed SAMs of carboxyl- and hydroxyl-terminated alkanethiols. The quality of the functionalization procedure was validated using X-ray photoelectron spectroscopy (surface chemical composition), AFM imaging (surface morphology in aqueous solution), and surface plasmon resonance (SPR, specific binding in aqueous solution). AFM force-distance curves recorded at a loading rate of 5000 pN/s between Fv- and Lyso-modified surfaces yielded a distribution of unbinding forces composed of integer multiples of an elementary force quantum of approximately 50 pN that we attribute to the rupture of a single antibody-antigen pair. Injection of a solution containing free Lyso caused a dramatic reduction of adhesion probability, indicating that the measured 50 pN unbinding forces are due to the specific antibody-antigen interaction. To investigate the dynamics of the interaction, force-distance curves were recorded at various loading rates. Plots of unbinding force vs log(loading rate) revealed two distinct linear regimes with ascending slopes, indicating multiple barriers were present in the energy landscape. The kinetic off-rate constant of dissociation (k(off) approximately = 1 x 10(-3) s(-1)) obtained by extrapolating the data of the low-strength regime to zero force was in the range of the k(off) estimated by SPR.     

Energy Landscape of Aptamer/Protein Complexes Studied by Single‐Molecule Force Spectroscopy

Aptamers are single‐stranded nucleic acid molecules selected in vitro to bind to a variety of target molecules. Aptamers bound to proteins are emerging as a new class of molecules that rival commonly used antibodies in both therapeutic and diagnostic applications. With the increasing application of aptamers as molecular probes for protein recognition, it is important to understand the molecular mechanism of aptamer–protein interaction. Recently, we developed a method of using atomic force microscopy (AFM) to study the single‐molecule rupture force of aptamer/protein complexes. In this work, we investigate further the unbinding dynamics of aptamer/protein complexes and their dissociation‐energy landscape by AFM. The dependence of single‐molecule force on the AFM loading rate was plotted for three aptamer/protein complexes and their dissociation rate constants, and other parameters characterizing their dissociation pathways were obtained. Furthermore, the single‐molecule force spectra of three aptamer/protein complexes were compared to those of the corresponding antibody/protein complexes in the same loading‐rate range. The results revealed two activation barriers and one intermediate state in the unbinding process of aptamer/protein complexes, which is different from the energy landscape of antibody/protein complexes. The results provide new information for the study of aptamer–protein interaction at the molecular level.     

Dynamic simulations of adhesion and friction in chemical force microscopy

A hybrid molecular simulation approach has been applied to investigate dynamic adhesion and friction between a chemical force microscope (CFM) tip and a substrate, both modified by self-assembled monolayers (SAMs) with hydrophobic methyl (CH(3)) or hydrophilic hydroxyl (OH) terminal groups. The method combines a dynamic model for the CFM tip-cantilever system and a molecular dynamics (MD) relaxation technique for SAMs on Au(111) at room temperature. The hybrid simulation method allows one to simulate force-distance curves (or adhesion) and friction loops (or friction coefficient) in the CFM on the experimental time scale for the first time. The simulation results also provide valuable molecular information at the interface that is not accessible in CFM experiments, such as the actual tip position with respect to the cantilever support position, molecular and hydrogen-bonding structures at the interface, and load distributions among different molecular chains (or single-molecule forces). Results show that the adhesion force and friction coefficient for the OH/OH contact pair are much larger than those for the CH(3)/CH(3) pair due to the formation of hydrogen bonds. During the retraction of a CFM tip from a surface, the CFM tip is away from the sample surface slightly while the spring undergoes dramatic elongation in the normal direction before rupture occurs. Single-molecule forces are distributed unevenly at the contact area. Surface energies calculated for functionalized surfaces compare well with those determined by experiments.     

A versatile "multiple fishhooks" approach for the study of ligand-receptor interactions using single-molecule atomic force microscopy

Despite the powerfulness of atomic force microscopy (AFM)-based single-molecule force spectroscopy in the study of ligand-receptor interactions, complicated cantilever functionalization and data interpretation have often been a great hurdle for its widespread application. Here, we present a much simplified experimental scheme by using a "multiple fishhooks" approach. In this strategy, multiple ligands are labeled on a single polymer chain, which forms complexes with receptors anchored on the substrate surface. Therefore, multiple single-bond rupture events can be captured in the same force-extension curves, similar to the widely used polyprotein approach. This method also allows nonsingle-molecule events and nonspecific interactions between cantilever and surface to be readily excluded from real data pool and greatly increases the quality and quantity of single-molecule data. The biggest advantage of our approach over the previously reported one is the choice of a naturally occurring polysaccharide, hyaluronan, the conformation of which in solution can be fine-tuned by pH, as the polymer backbone of the "multiple fishhooks" handle. Furthermore, our approach greatly simplifies the chemical synthesis of the polymer handle, allowing bioactive molecules to be easily one-step labeled on the handles in aqueous solution. We validate this strategy using the widely studied streptavidin-biotin system, and our single-molecule AFM results are in good agreement with previously reported ones. We anticipate that this novel strategy can be used as a versatile tool to study other complex and challenging ligand-receptor interactions.     

Single-molecule force spectroscopy study of interaction between transforming growth factor beta1 and its receptor in living cells

Transforming growth factor beta1 (TGF-beta1) regulates many important cellular processes such as cell proliferation, differentiation, and apoptosis, etc. Its signaling is initiated by binding to and bringing together TGF-beta type II receptor (TbetaRII) and type I receptor (TbetaRI). However, it is not fully understood how the TGF-beta1 ligand-receptor interaction occurs in living cells and what is the molecular mechanism of the signaling complex TGF-beta1/TbetaRII/TbetaRI formation. In this study, we have investigated the interaction between TGF-beta1 and its receptors in living cells with single-molecule force spectroscopy for the first time. By positioning TGF-beta1-modified atomic force microscope (AFM) tips on the cells expressing fluorescent protein tagged TGF-beta receptors, the living-cell force measurement was realized with a combined fluorescence microscope and AFM. We found that coexpression of TbetaRI with TbetaRII enhanced the binding force of TGF-beta1 with its receptors, whereas the expressed TbetaRI itself exhibited no binding affinity to TGF-beta1. Moreover, the unbinding dynamics of TGF-beta1/TbetaRII and TGF-beta1/TbetaRI/TbetaRII were investigated with dynamic force spectroscopy under different AFM loading rates. The dissociation rate constants of TGF-beta1 with its receptors as well as other parameters characterizing their dissociation pathways were obtained. The results suggested a more stable binding of TGF-beta1 with the receptor after TbetaRI is recruited and the important contribution of TbetaRI to the signaling complex formation during TGF-beta1 signaling.     

Unbinding of nicotine from the acetylcholine binding protein: steered molecular dynamics simulations

The ligand binding/unbinding process is critical to our understanding of the pharmacology of both the nicotinic acetylcholine receptor (nAChR) and the acetylcholine binding protein (AChBP). Steered molecular dynamics simulations were performed to learn about the unbinding process of the full agonist nicotine. Three different pulling models were designed to investigate the possible binding/unbinding pathways: radial and tangent models, and also a mixed model. Of the three, the tangent pulling model finally failed to dissociate nicotine from the ligand binding pocket. The efficiency of the pulling force profiles was superior, and the opening of the C-loop was smaller in the mixed pulling model than that in the radial model. The most favorable pathway for the cholinergic agonist nicotine to enter or leave the binding pocket is through the principal binding side, following a curvilinear track. Noticeably, it has been seen that the unbinding of the nicotine is concomitant with a global rotation of the protein-ligand complex which could be caused by the interactions of the ligand with protein at the tangent direction.     

Hydrophobic attraction as revealed by AFM force measurements and molecular dynamics simulation

Spherical calcium dioleate particles ( approximately 10 mum in diameter) were used as AFM (atomic force microscope) probes to measure interaction forces of the collector colloid with calcite and fluorite surfaces. The attractive AFM force between the calcium dioleate sphere and the fluorite surface is strong and has a longer range than the DLVO (Derjaguin-Landau-Verwey-Overbeek) prediction. The AFM force between the calcium dioleate sphere and the mineral surfaces does not agree with the DLVO prediction. Consideration of non-DLVO forces, including the attractive hydrophobic force and the repulsive hydration force, was necessary to explain the experimental results. The non-DLVO interactions considered were justified by the different interfacial water structures at calcite- and fluorite-water interfaces as revealed by the numerical computation experiments with molecular dynamics simulation.     

Adhesion Mechanisms of the Contact Interface of TiO(2) Nanoparticles in Films and Aggregates

Fundamental knowledge about the mechanisms of adhesion between oxide particles with diameters of few nanometers is impeded by the difficulties associated with direct measurements of contact forces at such a small size scale. Here we develop a strategy based on AFM force spectroscopy combined with all-atom molecular dynamics simulations to quantify and explain the nature of the contact forces between 10 nm small TiO(2) nanoparticles. The method is based on the statistical analysis of the force peaks measured in repeated approaching/retracting loops of an AFM cantilever into a film of nanoparticle agglomerates and relies on the in-situ imaging of the film stretching behavior in an AFM/TEM setup. Sliding and rolling events first lead to local rearrangements in the film structure when subjected to tensile load, prior to its final rupture caused by the reversible detaching of individual nanoparticles. The associated contact force of about 2.5 nN is in quantitative agreement with the results of molecular dynamics simulations of the particle-particle detachment. We reveal that the contact forces are dominated by the structure of water layers adsorbed on the particles' surfaces at ambient conditions. This leads to nonmonotonous force-displacement curves that can be explained only in part by classical capillary effects and highlights the importance of considering explicitly the molecular nature of the adsorbates.     

Mass transfer limitations at crystallizing interfaces in an atomic force microscopy fluid cell: a finite element analysis

Although atomic force microscopy (AFM) has emerged as the preeminent experimental tool for real-time in situ measurements of crystal growth processes in solution, relatively little is known about the mass transfer limitations that may impact these measurements. We present a continuum analysis of flow and mass transfer in an atomic force microscope fluid cell during crystal growth, using data acquired from calcium oxalate monohydrate (COM) crystal growth measurements as a comparison. Steady-state flows and solute concentration fields are computed using a three-dimensional, finite element method implemented on a parallel supercomputer. Steady-state flow results are compared with flow visualization experiments to validate the model. Computations of the flow field demonstrate how nonlinear momentum transport alters the spatial structure of the flow with increasing flow volume, altering mass transport conditions near the AFM cantilever and tip. The simulations demonstrate that the combination of solute depletion from crystal growth and mass transfer resistance lowers the solute concentration in the region between the tip and the crystal compared with the solute concentration at the inlet of the AFM cell. For example, using experimentally measured growth rates for COM, the solute concentration in this region is 3.1% lower than the inlet value because the solute consumed by crystal growth beneath the AFM tip cannot be replenished fully due to mass transport limitations. The simulations also reveal that increasing the flow rate through the cell does not affect this difference significantly because of the inherent shielding by the AFM tip in proximity with the crystal surface. Models such as the one presented here, used in conjunction with AFM measurements, promise more precise interpretations of measurement data.     

Analysis of the unfolding process of green fluorescent protein by molecular dynamics simulation

Molecular dynamics simulation of the enforced stretching of circularly permuted green fluorescent protein (cpGFP) was performed to observe the detailed process of unfolding of beta-sheets in cpGFP and to clarify the structural change arising from the force. The simulation using the generalized Born method with original force field parameters enabled us to observe the unfolding process of the entire region of the protein and to clarify atom motion of the individual domain during the stretching. The force required for the stretching of cpGFP was estimated from the differential of the computed potential energy. A prominent rise in force appeared three times during the stretching. The amplitude and the position of these three peaks were consistent with the observation in atomic force microscopy (AFM) experiments. Further, the movements of atoms involved in each peak were shown to be closely related to the dissociation of hydrogen bonds. Additional simulations for the unfolding process of titin and spectrin also gave satisfactory interpretation of the results of previous AFM experiments. The difference in the enforced stretching process between cpGFP and wild-type GFP was further discussed through the MD simulation.     

Force measurements on cell repellant and cell adhesive alginic acid coated surfaces

The atomic force microscope (AFM) was used to acquire force versus distance curves between the cantilever tip and samples bearing a surface overlayer of covalently linked alginic acid. The alginic acid coating resists cell-adhesion in in vitro experiments involving a normal and a tumor cell line. However, the surface becomes cell adhesive when alginic acid coated samples are subjected to glow discharge treatment. Force curves show in both cases the typical features resulting from the interaction between the cantilever tip and a hydrophilic, compressible polymer overlayer, suggesting that in both cases a diffuse interface with water exists. Following some recent findings on oligoethyleneglycol-terminated self-assembled-monolayers, it is suggested that conformational and molecular aspects of hydrophilic surface layers, rather than steric repulsion effects, could play a significant role in the mechanism that controls resistance to bio-adhesion.     

Atomic-scale friction on diamond: a comparison of different sliding directions on (001) and (111) surfaces using MD and AFM

Atomic force microscopy (AFM) experiments and molecular dynamics (MD) simulations were conducted to examine single-asperity friction as a function of load, surface orientation, and sliding direction on individual crystalline grains of diamond in the wearless regime. Experimental and simulation conditions were designed to correspond as closely as state-of-the-art techniques allow. Both hydrogen-terminated diamond (111)(1 x 1)-H and the dimer row-reconstructed diamond (001)(2 x 1)-H surfaces were examined. The MD simulations used H-terminated diamond tips with both flat- and curved-end geometries, and the AFM experiments used two spherical, hydrogenated amorphous carbon tips. The AFM measurements showed higher adhesion and friction forces for (001) vs (111) surfaces. However, the increased friction forces can be entirely attributed to increased contact area induced by higher adhesion. Thus, no difference in the intrinsic resistance to friction (i.e., in the interfacial shear strength) is observed. Similarly, the MD results show no significant difference in friction between the two diamond surfaces, except for the specific case of sliding at high pressures along the dimer row direction on the (001) surface. The origin of this effect is discussed. The experimentally observed dependence of friction on load fits closely with the continuum Maugis-Dugdale model for contact area, consistent with the occurrence of single-asperity interfacial friction (friction proportional to contact area with a constant shear strength). In contrast, the simulations showed a nearly linear dependence of the friction on load. This difference may arise from the limits of applicability of continuum mechanics at small scales, because the contact areas in the MD simulations are significantly smaller than the AFM experiments. Regardless of scale, both the AFM and MD results show that nanoscale tribological behavior deviates dramatically from the established macroscopic behavior of diamond, which is highly dependent on orientation.     

Quantitative measurement of friction between single microspheres by friction force microscopy

The sliding friction between single silica microspheres was examined by applying friction force microscopy to probe the interaction between spherical silica particles glued to a tipless atomic force microscopy (AFM) cantilever and another particle glued to a glass slide. A three-dimensional model handling the complex contact geometry between spherical particles was established to compute friction and normal forces at the sliding interface from measured deflections of the AFM cantilever. Results obtained at different loads show a linear relationship between friction and normal force, with a friction coefficient of 0.4 between silica spheres. Friction in this system occurs at multi-asperity contacts. The results show that the macroscopic friction law of Amontons can be used to model the friction behavior of micrometer-sized granular matter. For plasma-treated silica particles, increased friction as well as wear could be observed during sliding.     

AFM study of the behavior of polystyrene and glass particles during the electrodeposition of copper

In this paper, the behavior of polystyrene and glass particles on a copper electrode during the electrodeposition of copper was studied using an atomic force microscope (AFM). Polystyrene or glass particles glued to the tip of the AFM cantilever were kept in contact with the surface of the electrode. The surface forces between the polystyrene or glass particle and the copper electrode were measured before, during, and after electrodeposition. These experiments revealed that glass particles do not make contact with the electrode, probably due to the repulsive hydration force. Polystyrene particles, on the other hand, make contact with the electrode, due to the attractive hydrophobic force. The AFM experiments were correlated with sedimentation co-deposition experiments of polystyrene and glass particles with copper. It was found that 80% of the polystyrene particles added to the plating solution incorporated with copper, while only 0.25% of the glass particles co-deposited under the same conditions.     

CPMD Simulation of a Bimolecular Chemical Reaction: Nucleophilic Attack of a Disulfide Bond under Mechanical Stress

Previous single‐molecule atomic force microscopy (AFM) experiments showed a change in the reactivity of a bimolecular substitution reaction with a definite force acting on a protein containing disulfide bonds. Using Car–Parrinello molecular dynamics (CPMD) simulations, we analyse the relevant reaction pathways for the breaking of a disulfide bond in the presence of nucleophiles.     

分子间相互作用力的直接测量

结合我们近期的研究工作,着重介绍如何将分子组装与单分子力谱相结合,从单分子水平直接研究分子间相互作用力,包括π-π相互作用、多价作用及嵌入作用.在实验中,将一个相互作用单元通过高分子间隔基共价连接到原子力显微镜针尖上,并将另一相互作用单元共价修饰到基底上,通过压电陶瓷管的移动获得力与拉伸长度的曲线.高分子柔性间隔基团的引入既可用来作为判别单链拉伸的"内标",又可避免非特异相互作用对待测的特异相互作用的影响.研究表明,结合静态和动态力学谱,不仅能够实现分子间相互作用力的直接测量,而且还可获得解离速率和相互作用的距离等参数.     

Direct Identification of Protein–Protein Interactions by Single‐Molecule Force Spectroscopy

Single‐molecule force spectroscopy based on atomic force microscopy (AFM‐SMFS) has allowed the measurement of the intermolecular forces involved in protein‐protein interactions at the molecular level. While intramolecular interactions are routinely identified directly by the use of polyprotein fingerprinting, there is a lack of a general method to directly identify single‐molecule intermolecular unbinding events. Here, we have developed an internally controlled strategy to measure protein–protein interactions by AFM‐SMFS that allows the direct identification of dissociation force peaks while ensuring single‐molecule conditions. Single‐molecule identification is assured by polyprotein fingerprinting while the intermolecular interaction is reported by a characteristic increase in contour length released after bond rupture. The latter is due to the exposure to force of a third protein that covalently connects the interacting pair. We demonstrate this strategy with a cohesin–dockerin interaction.