Force-Regulated Spontaneous Conformational Changes of Integrins α5β1 and αVβ3

Integrins are cell surface nanosized receptors crucial for cell motility and mechanosensing of the extracellular environment, which are often targeted for the development of biomaterials and nanomedicines. As a key feature of integrins, their activity, structure and behavior are highly mechanosensitive, which are regulated by mechanical forces down to pico-Newton scale. Using single-molecule biomechanical approaches, we compared the force-modulated ectodomain bending/unbending conformational changes of two integrin species, α5β1 and αVβ3. It was found that the conformation of integrin α5β1 is determined by a threshold head-to-tail tension. By comparison, integrin αVβ3 exhibits bistability even without force and can spontaneously transition between the bent and extended conformations with an apparent transition time under a wide range of forces. Molecular dynamics simulations observed almost concurrent disruption of ∼2 hydrogen bonds during integrin α5β1 unbending, but consecutive disruption of ∼7 hydrogen bonds during integrin αVβ3 unbending. Accordingly, we constructed a canonical energy landscape for integrin α5β1 with a single energy well that traps the integrin in the bent state until sufficient force tilts the energy landscape to allow the conformational transition. In contrast, the energy landscape of integrin αVβ3 conformational changes was constructed with hexa-stable intermediate states and intermediate energy barriers that segregate the conformational change process into multiple small steps. Our study elucidates the different biomechanical inner workings of integrins α5β1 and αVβ3 at the submolecular level, helps understand their mechanosignaling processes and how their respective functions are facilitated by their distinctive mechanosensitivities, and provides useful design principles for the engineering of protein-based biomechanical nanomachines.

I ntegrins are a family of heterodimeric transmembrane molecules on the surface of nearly all cells.−5 In this context, it becomes crucial to understand the mechanosensitivity of integrins, because it not only mediates how cells interact with the (patho)physiological environment, but also critically affects cells' compatibility and interaction with the biophysical properties of therapeutic tools and agents.For instance, the elasticity of nanoparticles has been shown to affect their in vivo localization and therapeutic efficacy. 6,7Of the 24 integrin species currently known, integrins α 5 β 1 and α V β 3 are used by a variety of cells to bind the extracellular matrix (ECM) and form focal adhesion.−10 The molecular basis of such functional distinction is unclear, which was vaguely suggested to be related to the structural differences in α 5 β 1 and α V β 3 ectodomains. 9,10orce modulates the properties and functions of certain proteins by inducing conformational changes, such as coiling/ uncoiling, zipping/unzipping, and folding/unfolding.We previously showed that cell surface α L β 2 and α V β 3 integrins undergo force-modulated conformational changes, such that force facilitates unbending but suppresses bending, shifting the conformational equilibrium toward extension. 11,12From a mechanical perspective, it is surprising that integrins can spontaneously bend against a wide range of forces.It is intuitive that a head-to-tail tension can facilitate a bent integrin to unbend regardless of which conformation is more stable prior to force application, because force can tilt the energy landscape and allow the extended conformation to become more stable, if this is not already the case in the absence of force.However, even if the bent conformation is more stable, its spontaneous bending against a tensile force is still counterintuitive when the tension is as high as 10−20 pN.The mechanical work done by the integrin to bend back under a linearly increasing force is comparable to the free energy of biotin−avidin binding (∼35 k B T), 13 one of the strongest noncovalent interactions, and much greater than the environmental thermal agitation (0.5 k B T).
Our previous experiments on force-modulated integrin bending and unbending were performed on living cells, 11,12 where cell activity can regulate integrin conformational changes biologically.It seems natural to hypothesize that it is the cell that provides a "deactive energy" to bend the integrin against force.However, it is difficult to envision how this presumably cell-provided energy is converted into mechanical work to power the bending of the extended integrin, which occurs distally from the cell surface.To test this hypothesis, we used two force spectroscopic techniques to perform singlemolecule experiments on purified integrins α 5 β 1 and α V β 3 .While both integrins were able to undergo spontaneous bending and unbending under force, their conformational changes exhibited distinctive mechanical and kinetic properties, which may be related to their distinctive structures. 14pecifically, the conformation of integrin α 5 β 1 was mostly bent in Ca 2+ and mostly extended in Mn 2+ , suggesting a canonical energy landscape with a deep energy well that traps the integrin in the bent state in Ca 2+ and the extended state in Mn 2+ .Force could tilt the energy landscape to shift the system in Ca 2+ to bistability such that integrin α 5 β 1 would abruptly transition back-and-forth between the bent and extended states.In contrast, integrin α V β 3 might take either bent or extended conformation in both Ca 2+ /Mg 2+ and Mn 2+ conditions, and can undergo spontaneous bending and unbending with slow kinetics under a wide range of tensile forces similar to cell surface α V β 3 integrins, 11 falsifying our "biological energy" hypothesis and suggesting a physical mechanism.To explain the unusual behaviors of integrin α V β 3 , we developed a multistate conformational energy landscape for this integrin, which was supported by molecular dynamics simulations and could fit our experimental data well.The different mechanosensitivities of integrins α 5 β 1 and α V β 3 likely underlie their distinct biological functions in cell mechanosensing, which should help guide the development of more human compatible nanotherapeutics.The finding that integrin α V β 3 can spontaneously bend and unbend under a wide range of mechanical forces without cell environment or the supply of external energy provides inspirational design concepts for protein-based biomechanical nanomachines.

RESULTS/DISCUSSION
Directly Observing Single Integrin α 5 β 1 Unbending and Bending.Using the atomic force microscopy (AFM), we tested whether mechanical force could induce conformational changes of integrin α 5 β 1 independent of cell regulation.Recombinant integrin α 5 β 1 ectodomain with a human IgG Fc tag at the tail (α 5 β 1 -Fc) was captured on a polystyrene surface (Figure 1A), and driven to touch the fibronectin module III domains 7−10 (abbreviated as FN, containing both the RGD sequence and synergistic site 15 ) adsorbed on a cantilever tip to allow for bond formation.As has been confirmed in our previous work, the binding events measured by this experimental setup were predominantly mediated by specific interactions of α 5 β 1 -Fc (and trα 5 β 1 -Fc used below) with FN 16 .A tensile force was loaded on each integrin α 5 β 1 − FN bond, which was ramped by retracting the polystyrene surface until reaching 20 pN, and then unloaded to 0 pN at the same rate (Figure 1B).Inspection of the force vs time traces often reveals a clearly visible kink in the middle of both the loading and unloading phases, where the slope of the curve suddenly drops from positive to zero or even negative in the loading phase and abruptly jumps from negative to zero or even positive in the unloading phase (Figure 1B).These kinks are clear indications of protein conformational changes such as unfolding−refolding. 17 To identify the origin of these conformational changes, we first replaced α 5 β 1 -Fc with a truncated construct that contains only integrin α 5 β 1 headpiece (trα 5 β 1 -Fc; Figure 1A).In all four cation conditions�2 mM Ca 2+ (Ca 2+ ), 1 mM Ca 2+ plus 1 mM Mg 2+ (Ca 2+ /Mg 2+ ), 1 mM Mg 2+ plus 1 mM EGTA (Mg 2+ / EGTA) and 2 mM Mn 2+ (Mn 2+ )�that favor different integrin conformations, the conformational changes seen in the fulllength α 5 β 1 -Fc were no longer observed (Figure 1C), indicating that these conformational changes are from the α 5 β 1 ectodomain but not the recombinant Fc tail, GG7 or FN, and require the integrin α 5 β 1 tailpiece.Second, such conformational changes occurred progressively less frequently as the cation composition changed to those that activate integrins more and more potently, resulting in a frequency hierarchy of Ca 2+ > Ca 2+ /Mg 2+ > Mg 2+ /EGTA > Mn 2+ (Figure 1C).This suggests the observed structural lengthening/shortening events to be those of integrin unbending/bending, which explains the frequency hierarchy: the activating cation conditions facilitate more integrins to adopt the extended conformation, leaving less integrins in the bent conformation capable of unbending.Third, in all cation conditions, the structural lengthening in the loading phase was almost always ensued by a structural shortening in the unloading phase.By plotting force against the molecular extension of the bond, we found that the molecular complex fully recovers to its original length at the end of the loading−unloading cycle with no hysteresis (Figure 1D) and the loading and unloading phases largely overlap (Figure 1E), suggesting that the conformational changes are highly reversible and ruling out the alternative interpretation that they represent irreversible structure denaturation.Fourthly, the slope of the force−extension curve, which represents molecular stiffness, was seen to always increase after a structural lengthening (Figure 1E,G,H) and decrease after a structural shortening.The curve segments before and after the structural extension were both well-fitted by a linear model, while the worm-like chain (WLC) model did not render better results (Supp.Figure 1A).Furthermore, the slope of the force− extension curve remained constant in the absence of a structural change (Supp.Figure 1B), ruling out the possibility that the molecular stiffening was caused by the WLC nonlinear response.Together, these results indicate that α 5 β 1 becomes stiffer after the structural lengthening and softer after the structural shortening.This is consistent with our previous observations on multiple other integrin species that integrins are stiffer in their extended conformation than in the bent conformation. 11,12,18Using molecular stiffness as a signature readout of integrin conformation, we found that the value of the post-extension stiffness is comparable to integrins in Mn 2+ showing no structural changes (Figure 1F,H,I), which agrees with our hypothesis that integrins would be nearly unable to unbend or bend in this cation condition, because Mn 2+ has already activated most of the integrins to the extended conformation.Finally, the molecular extension change due to structural lengthening centers around 7.5 nm (Figure 1J), which is comparable to the head-to-tail length increase of an unbending integrin α 5 β 1 characterized by theoretical modeling. 19The broad distribution of the extension change was most likely due to the intermolecular variation in the headpiece/ tailpiece angle of both the bent and extended conformations in different α 5 β 1 molecules. 20,21o further validate our discovery, a recombinant integrin α 5 β 1 ectodomain fused with a polyhistidine tag at the Cterminus (α 5 β 1 -Poly-His) was tested on another force spectroscopy technique, Biomembrane Force Probe (BFP).The BFP setup consisted of a micropipette-aspirated human red blood cell (RBC) with a streptavidin (SA) and FN cofunctionalized probe bead attached to its apex to serve as an ultrasensitive force transducer (Figure 1K, left).A bead coated with α 5 β 1 -Poly-His was aspirated by an opposing micropipette (Figure 1K, right) and driven to repeatedly contact the probe bead to induce integrin α 5 β 1 −FN bond formation.The specificity of the bonds was confirmed with the fact that addition of the mAb BMC5 eliminated most of the binding events (Figure 1L).Similar to our observation in the AFM assay, a "kink" was observed in the force vs time traces when the bonds were ramped by force (Figure 1M), which was observed much more frequently in Ca 2+ condition than in Mn 2+ (Figure 1L).The molecular extension change induced by unbending centers at 8.8 nm (Figure 1N), which was slightly longer than the value collected by AFM.This was likely due to the softer force transducer of BFP that does not favor the observation of small molecular extensions in the ramping phase, which could cause a bias in the event detection.Together, these results indicate that real-time integrin α 5 β 1 unbending and bending events could be observed using our force spectroscopy approaches.

ACS Nano
Directly Observing Single Integrin α V β 3 Bending and Unbending.The biophysical characteristics of force-modulated integrin α 5 β 1 unbending and rebending are different from those previously characterized for integrins α L β 2 and α V β 3 on the cell surface.The force range within which structural change events could be observed is quite narrow for integrin α 5 β 1 (<10 pN, Figure 1B,E; also see Figure 5A−D below) but much wider for α L β 2 and α V β 3 (up to 40 pN). 11,12On the other hand, kinetics rates of the conformational changes can be quantified by two parameters: time-to-switch (t 0± ) is the waiting time required for the conformational switch to occur, while switching time (t sw± ) is the time taken for the conformational switch from the start to finish. 11,12Using these definitions, we found that the kinetics are much more rapid for integrin α 5 β 1 than cell surface integrins α L β 2 and α V β 3 (e.g., t 0− at 5 pN is ∼0.02 s vs 2−3 s). 11,12A hypothetical explanation for the different (un)bending behaviors observed here and previously might be the absence of cell regulation for integrin α 5 β 1 .To test this hypothesis, we studied the conformational changes of integrin α V β 3 ectodomain bound to FN, so as to allow the direct comparison of integrins α 5 β 1 and α V β 3 as purified proteins.The AFM approach was first applied.Despite that control experiments confirmed the detection of integrin α V β 3 −FN specific binding, no "kink", i.e., sudden slope change in the force-time curve, was observed over hundreds of force loading−unloading events; clamping the integrin α V β 3 −FN bonds under a constant force or applying cyclic forces did not yield any kink type of conformational changes either.Cell surface integrin α V β 3 (un)bending events were previously observed using BFP 11 .We reasoned that the observation of integrin α V β 3 conformational changes might not be favored by the stiff AFM cantilever but favored by the soft BFP force sensor (spring constants ∼3 vs ∼0.3 pN/nm), because a 10-nm head-to-tail length change caused by integrin α V β 3 bending would result in a ∼30 pN force increase in the AFM, which would severely inhibit bending, but only ∼3 pN force increase in the BFP, which would not.We thus used the BFP for testing, wherein the probe bead was again cofunctionalized with SA and FN, while the target bead was coated with recombinant α V β 3 protein.Binding specificity was confirmed by an integrin α V β 3 -blocking mAb, LM609, which abolished most of the binding events (Figure 2A).
After the FN coating was titrated on probe beads to lower the adhesion frequency to 20%, a necessary condition for most binding events to be mediated by single bonds, 22 integrin α V β 3 was then interrogated under both Ca 2+ /Mg 2+ and Mn 2+ conditions using distance-clamping assay: 11 the integrin α V β 3 −FN bond was first pulled to a certain force level, and the target bead was then clamped at the position until bond dissociation.The bonds could sustain a wide range of forces, with lifetimes much longer in Mn 2+ than in Ca 2+ /Mg 2+ , consistent with the activating role of Mn 2+ (Figure 2B).Unlike the purified integrin α 5 β 1 −FN interaction that forms catch-slip bonds not only in Mn 2+ , but also in Ca 2+ /Mg 2+ and Mg 2+ / EGTA, 16 the purified integrin α V β 3 −FN interaction formed a catch-slip bond in Mn 2+ but a slip-only bond in Ca 2+ /Mg 2+ .This slip-only bond indicates the limited effect of a sustained force to strengthen integrin α V β 3 bonding to FN, which agrees with our previously reported weak integrin α V β 3 −FN catch-slip bond on the cell surface. 11n the clamping phase of some lifetime measurements, we observed integrin α V β 3 unbending or bending events, respectively signified by a concurrent decrease in the mean force and force fluctuation or a concurrent increase in the mean force and force fluctuation (Figure 2C−F; Supp.Table 1), which are clearly distinguishable from formation of an additional bond (signified by an increase in the force and a decrease in thermal fluctuation) and dissociation of a bond from a multibond adhesion (signified by a decrease in the force and an increase in thermal fluctuation). 11,12In most cases, instead of successive back-and-forth transitions, only a single conformational change event could be observed in the distance-clamp cycle, which is likely due to the limitation of integrin α V β 3 −FN bond lifetimes that were too short to provide a long enough observation window to overcome the slow kinetics of integrin α V β 3 conformational changes (shown below).Unlike purified integrin α 5 β 1 and consistent with cell surface integrin α V β 3 , the conformational changes of purified integrin α V β 3 occurred under a wide range of forces (Figure 3A,B) with relatively long time-to-switch (t 0+ and t 0− respectively for unbending and bending) and switching time (t sw+ and t sw− respectively for unbending and bending) 11 (cf.Figure 2C−F).Such unusually slow kinetics ruled out the alternative possibility that these conformational changes were protein domain unfolding/refolding events, which are generally abrupt (e.g., talin 23 ) due to the involved local secondary structure being relatively simple.Replacing FN on the probe beads with LIBS-2, a mAb that binds the α V β 3 βTD domain at its tailpiece, 24 abolished the above signature signals for integrin conformational changes (Supp.Table 1) despite the long lifetimes (Figure 2B), further ruling out the alternative possibility that the putative bending/unbending events are due to multiple bond rupture/formation or instrumental drift.Adding high-concentration LIBS-2 to the solution, which stabilizes β 3 integrins in the extended conformation, 25 also eliminated all bending events (Supp.Table 1).Interestingly, LIBS-2 treatment did not significantly alter the integrin α V β 3 − FN bond type and lifetimes in Ca 2+ /Mg 2+ (Figure 2B), indicating that integrin extension and catch-slip bond formation are decoupled.During (un)bending, the change in the RBC elongation (Figure 2D,F, cyan shaded areas) is equal to the change in the integrin head-to-tail length. 11,12These length changes of both unbending and bending events follow a single-Gaussian distribution (Figure 3A,B) with an indistinguishable average value of ∼13 nm in both Ca 2+ /Mg 2+ and Mn 2+ conditions (Figure 3C), agreeing with our previous observations of cell surface integrin α V β 3 bending/unbending events and with our MD simulation results on integrin α V β 3 unbending. 11,26The length changes observed here on integrin α V β 3 are much longer than α 5 β 1 (Figure 1J,N), which is primarily due to the difference between the two integrin species.Bent integrin α V β 3 adopts a highly compact structure with a headpiece-tailpiece angle of ∼40°(refs 27, 28; also shown in Figure 6E below).In contrast, the bent conformation of integrin α 5 β 1 is less tight where the headpiece-tailpiece angle reaches 71°−93°(ref 21; also shown in Figure 6A below), therefore shortening the traveling distance of its headpiece during conformational changes.
Furthermore, the stiffness of the integrin α V β 3 −FN complex is lower before unbending than before bending (Figure 3D), consistent with the signature integrin stiffening upon unbending. 11,12,18Since the stiffness depends only on the conformation but not the cation condition (Figure 3D), we pooled data from both cation conditions together to examine the stiffness distributions for the bent and extended integrins, finding their respective means and standard deviations of 0.55 ± 0.20 and 0.73 ± 0.24 pN/nm (Figure 3E), comparable to the values previously measured from cell surface α V β 3 (ref 11).Moreover, we plotted the histograms of additional stiffness measurements from each cation condition, regardless of whether integrin (un)bending events were observable, and fitted each by a dual-Gaussian distribution using 0.55 and 0.73 pN/nm as the two means to calculate the proportions of integrins in the bent and extended states.We found that, of those α V β 3 integrins that formed bonds, 62.8% were in the extended conformation in Mn 2+ , but 24.9% were in the extended conformation in Ca 2+ /Mg 2+ (Figure 3F,G), consistent with the activating role of Mn 2+ .More importantly, these results confirm the previous observation that integrin α V β 3 , unlike α 5 β 1 , is already bistable under zero force. 28verall, the data confirm that purified integrin α V β 3 protein can spontaneously transition between the bent and extended conformations under a wide range of forces in the absence of cellular regulation or biological energy supply.
Integrin-mediated mechanosignaling was conventionally believed to require either integrins to cluster, so as to trigger rearrangement of cytoskeletal structure, 29,30 or alternatively, prior inside-out signaling to unbend the integrin for activation and ligand binding ("switch-blade" model 31 ), and/or to activate intracellular scaffold proteins (e.g., talin in "molecular clutch" model 32 ) for signal transduction.Our findings on integrins α 5 β 1 and α V β 3 , together with previous echoing works, 11,12 indicated that bent integrins can also bind to ligands and that integrin unbending can be solely modulated by mechanical force.These shreds of evidence suggest an additional mechanism that allows a single inactive integrin to initiate outside-in mechanosignaling without prior inside-out signaling, wherein the unbending conformational change propagates intracellularly to induce integrin tailpiece separation, 33,34 integrin cluster rearrangement 14,35 and/or the association of cytoplasmic proteins. 14tegrin α V β 3 Showed No Cyclic Mechanical Reinforcement Effect.Integrin spontaneous unbending and bending respectively decrease and increase its ligand binding force (Figure 2C,E), which may help strengthen the bonds through a mechanism called "cyclic mechanical reinforcement" (CMR), where a cyclic force applied to a receptor−ligand bond greatly prolongs its lifetime.CMR was initially observed with integrin α 5 β 1 −FN bonds, 36 but later also observed with actin−actin bonds. 37To test the CMR effect on integrin α V β 3 −FN bonds, we first used AFM as did previously on integrin α 5 β 1 −FN bonds. 36Once a bond was detected, two types of cyclic forces were applied: 1) one loading−unloading cycle that first peaks at 20 pN and then drops to and is held at 5 pN (Figure 4A); and 2) cyclic forces with zero, one, two or three complete loading−unloading cycles followed by ramping to and being clamped at a peak force of 10 pN (Figure 4C).Unexpectedly, neither type of cyclic forces prolonged α V β 3 − FN lifetimes, showing a lack of CMR effect (Figure 4B,D).
We also repeated the above experiments using BFP with integrin α V β 3 -expressing platelets as the target (Figure 4E).Inhibitory mAbs 10E5 and P1D6 were added to respectively block α IIb β 3 and α 5 β 1 , two other FN-binding integrins on platelets, to ensure sole interaction of integrin α V β 3 with FN 18 .The second type of force loading−unloading cycles was applied to integrin α V β 3 −FN bonds followed by ramping to and clamping at 10 pN.Despite the presence of cell environment, the bond lifetime of integrin α V β 3 with FN was still not prolonged by cyclic forces (Figure 4F).
Distinctive force-dependent kinetics of integrins α 5 β 1 and α V β 3 conformational changes.The distinctive biophysical behaviors in the conformational changes of α 5 β 1 and α V β 3 integrins prompted us to analyze and compare the kinetics of their bending and unbending conformational changes as characterized by switching time (t SW± ) and timeto-switch (t 0± ).We employed AFM to pull integrin α 5 β 1 slowly (∼1 nm/s) after performing CMR to strengthen its bond with FN, which prolonged the time for observation of repetitive unbending-bending cycles in a single binding event 36 (Figure 5A), allowing us to collect ensembles of measurements for kinetic analysis.The bending and unbending processes were too fast to measure t SW± values (always beyond the temporal resolution of 1 ms of our AFM instrument) and the individual t 0± values were highly fluctuating (Figure 5A, inset).Nevertheless, the average ⟨t 0+ ⟩ decreased exponentially, and ⟨t 0− ⟩ increased exponentially, with increasing force f (Figure 5B), behaving as a typical slip bond and catch bond, respectively. 38e model the force-dependent ⟨t 0+ ⟩ and ⟨t 0− ⟩ using the Bell equation 39 (eq 1a) and its "catch bond counterpart" (eq 1b): where k B the Boltzmann constant, T is absolute temperature, ⟨t 0± | f=0 ⟩ are the respective values of ⟨t 0± ⟩ at zero force, and Δx ± respectively represent the distances the top of the energy barrier to the bottoms of the energy wells of the bent (Δx + ) and extended (Δx − ) conformations in the energy landscape at zero force.Directly fitting eqs 1a and 1b to the respective ⟨t 0+ ⟩ and ⟨t 0− ⟩ data in Figure 5B yielded excellent agreement and returned ⟨t 0+ | f=0 ⟩= 5.5 ± 2.9 s, Δx + = 1.6 ± 0.16 nm, ⟨t 0− | f=0 ⟩ = 0.004 ± 0.0007 s, and Δx − = 2.4 ± 0.25 nm.Since the reciprocal average time-to-bending and reciprocal average time-to-unbending are the kinetic rates of bending and unbending, respectively, we can calculate the bending equilibrium coefficient as a function of force by taking the ratio of eq 1b to eq 1a, which yields where Δx = Δx + + Δx − .Let f 1/2 be the force at which ⟨t 0− ⟩/ ⟨t 0+ ⟩ = 1, i.e., the force at which the time-to-unbending ⟨t 0+ | f=fd 1/2 ⟩ equals to the time-to-bending ⟨t 0− | f=fd 1/2 ⟩.For α 5 β 1 , Δx = 4.0 ± 0.3 nm, and The definition of [⟨t 0− ⟩/⟨t 0+ ⟩] f=fd 1/2 = 1 predicts that near f 1/2 , integrin α 5 β 1 has an equal chance of residing in the bent and extended states.The value of ⟨t 0+ | f=fd 1/2 ⟩ = ⟨t 0− | f=fd 1/2 ⟩= 0.076 s predicts that the integrin transitions rapidly back-and-forth between these two states.Such consecutive back-and-forth events with brief intermittent durations were indeed observed, but occurred at comparable frequencies only in a narrow force range (6−9 pN, Figure 5A).The probability of time during which the integrin stays in the extended state can be derived from eq 2: Here f 1/2 is defined by the same formula but interpreted as the force at which the integrin has a 50−50 chance of staying in either the bent or extended state.We plotted the measured fraction of extension times (points) and the fitting of eq 3 (curve) to the data (Figure 5C), which showed excellent agreement and returned a slightly larger Δx = 4.40 ± 0.06 nm and a slightly smaller f 1/2 = 6.04 ± 0.01 pN.The consistency between the values obtained by fitting eqs 1a and 1b to the data in Figure 5B and those by fitting eq 3 to the data in Figure 5C supports the quality of our data, the appropriateness of our model, and the robustness of the model parameters.
Across the f 1/2 threshold, force quickly transitioned the integrin from the bent to extended conformation: as force increased from 4.3 to 10.5 pN, the dominant (>95%) population of integrin molecules rapidly changed from the bent to the extended conformation, which increased the population ratio of extended over bent integrins ⟨t 0− ⟩/⟨t 0+ ⟩ by 400-fold, corresponding to an average force sensitivity of >60fold/pN (Figure 5D).Such a high force-sensitivity is due to the relatively large Δx value and agrees with a previous theoretical study inferring that integrin α 5 β 1 unbending is ultrasensitive to force, 19 reflecting nearly "digital" modulation of force on α 5 β 1 conformation.
Compared with the conformational change kinetics of integrin α 5 β 1 , which were rapid and strongly force-dependent, the kinetics of integrin α V β 3 conformational changes were slow and weakly force-dependent.Such characteristics were revealed by using the same approaches as above to analyze the counterpart data for integrin α V β 3 , which occurred over a much broader range of forces (Figure 5E−H).Unlike integrin α 5 β 1 whose switching times t SW± were too brief to measure (Figure 5A), hence mimicking a digital on/off switch, the counterpart values for α V β 3 were long enough to be measurable, exhibiting the characteristic of a more gradual transition.Their ⟨t SW± ⟩ (Figure 5E) and ⟨t 0± ⟩ (Figure 5F) displayed similar trends.Compared to the Ca 2+ /Mg 2+ cation condition, activating the integrin with Mn 2+ resulted in slightly shorter ⟨t sw+ ⟩ and ⟨t 0+ ⟩ and longer ⟨t sw− ⟩ and ⟨t 0− ⟩ (Figure 5E,F), consistent with the known coupling between integrin extension and activation. 28ike integrin α 5 β 1 , increasing force decreased ⟨t 0+ ⟩ and ⟨t sw+ ⟩ but increased ⟨t 0− ⟩ and ⟨t sw− ⟩ of integrin α V β 3 (Figure 5E,F).Quantitatively, however, the response of kinetics to force was very different.Fitting eq 3 to the data in Figure 5G returned much larger f 1/2 values (22.67 ± 0.09 and 17.76 ± 0.14 pN in Mn 2+ and Ca 2+ /Mg 2+ , respectively) and much smaller Δx values (0.5 ± 0.1 and 0.7 ± 0.2 nm in Mn 2+ and Ca 2+ /Mg 2+ , respectively).These values predict that integrin α V β 3 can undergo bending and unbending at a much higher force level and under a much broader range of forces, agreeing with our experimental observations.The much weaker force-dependency of integrin α V β 3 bending/unbending kinetics can be seen in Figure 5H: within the force range of 4.0−28.0pN where sufficient events were collected for statistical analysis, the population ratio ⟨t 0− ⟩/⟨t 0+ ⟩ of extended over bent α V β 3 only increased by 8.8-fold (a force sensitivity of ∼0.37-fold/pN) in Mn 2+ and by 13-fold (a force sensitivity of ∼0.54-fold/pN) in Ca 2+ /Mg 2+ , revealing a >100-fold greater resistance to force modulation than integrin α 5 β 1 .We also reanalyzed our previously published data of force-dependent integrin α V β 3 bending/unbending conformational changes on cell surface, 11 finding f 1/2 and Δx values similar to cell-free integrin α V β 3 (Supp.Table 2), indicating that these conformational changes are mainly modulated by force but not the cell environment.Together, these results demonstrated distinctive mechanisms of force modulation on integrins α 5 β 1 and α V β 3 (un)bending: a "digital" modulation for α 5 β 1 and an "analogous" modulation for α V β 3 .
Interestingly, the distinctive mechanosensitivities of integrins α 5 β 1 and α V β 3 support their respective mechanosignaling roles in focal adhesion.The "digital" unbending of integrin α 5 β 1 by force allows the cell to quickly sense extracellular stretching above a threshold, and initiate integrin α 5 β 1 recruitment and clustering to form strong adhesion. 8Furthermore, around the threshold force (7.4 pN), integrin α 5 β 1 quickly switches backand-forth between the bent and extended conformations (>10 Hz), which could trigger fast oscillation in binding force magnitude, and therefore the strong CMR effect 36 of integrin α 5 β 1 to reinforce adhesion.On the other hand, the "analogous" modulation gradually shifts the conformational equilibrium of integrin α V β 3 over a wide force range.This enables each integrin α V β 3 molecule to act as a "ruler" for the cell to "measure" the local extracellular stretching force and matrix rigidity.As a result, when expressed on the same cell, the two integrin species can cooperate to allow the cell to both quickly adhere to the substrate and sense substrate stiffness.−43 Explaining the Distinctive Switching Times of Integrins α 5 β 1 and α V β 3 Conformational Changes by MD Simulations.The orders of magnitude longer t sw± of integrin α V β 3 than α 5 β 1 is intriguing.To explain this difference, we hypothesize that α V β 3 conformational changes may involve a random sequence of formation/disruption of hydrogen bonds (H-bonds) that does not occur for integrin α 5 β 1 (refs  44, 45), resulting in a slower and more complex submolecular process for integrin α V β 3 than α 5 β 1 .To test this hypothesis, we performed steered molecular dynamics (SMD) simulations on integrins α 5 β 1 (PDB code 7NXD) and α V β 3 (PDB code 3IJE) by applying external pulling forces to their headpieces.Unbending of both integrins was accompanied by the disruption of H-bonds, with larger numbers in integrin α V β 3 than α 5 β 1 (Supp.Figure 2).To acquire more quantitative information while minimizing the artifact introduced by fast force loading, we further performed free molecular dynamics (MD) simulations with integrins α 5 β 1 and α V β 3 at their bent conformation without loading or restraint.In addition, we obtained from the above SMD simulations 3−4 intermediate structures with different head-to-tail lengths, from 14 to 18 nm for α 5 β 1 (Figure 6A) and from 6 to 18 nm for α V β 3 (Figure 6E) and carried out MD simulations on these structures with the molecular length restrained.H-bonds were observed to constantly form and break between the headpiece and tailpiece of both integrins in their respective bent conformations, but the time-averaged numbers differed greatly: ∼2 and ∼7, respectively, for integrins α 5 β 1 and α V β 3 (Figure 6C,G).As the length of integrin α 5 β 1 increased, its ∼2 H-bonds were rapidly disrupted during the initial phase of unbending at an average rate of ∼1 bond/nm of extension (Figure 6A−C).In contrast, the ∼7 H-bonds in integrin α V β 3 were disrupted much slower which occurred across the whole course of unbending (∼0.4 bond/nm of extension) (Figure 6E−G), requiring nearly an order of magnitude longer extension to break all the H-bonds than integrin α 5 β 1 .Interestingly, in integrin α 5 β 1 5 out of the 7 most frequently formed H-bonds were in the integrin knee region (Figure 6D).This contrasts with integrin α V β 3 where the most frequently formed 8 H-bonds were spatially equally distributed along the headpiece-tailpiece interface, and only 2 of them were in the knee region (Figure 6H).Among them, the H-bond most proximal to the integrin knee region (R8-E522) was not disrupted until the integrin reached full extension, whereas H-bonds distal to the knee region (e.g., R633-D393) were disrupted as soon as integrin α V β 3 started to unbend (Figure 6H).These results indicate a direct correlation between the H-bonds' distance to the knee and the chronological sequence of their disruption.The above observations help explain the distinctive (un)bending dynamics of the two integrins studied here and provide the rationale for the energy landscapes and transition models below.
Constructing Energy Landscapes and Transition Models for Integrins α 5 β 1 and α V β 3 Conformational Changes.We wished to construct the corresponding energy landscapes and transition kinetic models for integrins α 5 β 1 and α V β 3 , using the parameters listed in Suppl.Table 2. Noting that Δx and k T ln are the respective differences in the reaction coordinates and energies of the bottoms of the two energy wells for the bent and extended states at zero force, we first built an energy landscape for integrin α 5 β 1 (Figure 7A).Without force, integrin α 5 β 1 dominately stays in the bent conformation (Supp.Video 1).Force tilts the energy landscape such that the energy difference vanishes at f 1/2 , i.e., ln 0 . Thus, around f 1/2 integrin α 5 β 1 switches back and forth between the bent and extended conformations indefinitely using the energy from thermal agitations to hop over the force-tilted energy barrier separating the two states (Figure 7A, Supp.Video 2).With the force further increased, the energy well of the extended conformation is further deepened and the integrin α 5 β 1 mainly stays in the extended conformation (Supp.Video 3).However, such an energy landscape may not be appropriate for integrin α V β 3 , although its force-dependent ⟨t 0± ⟩ data (Figure 5F) could still be fitted by the Bell model. 39This is because the above energy landscape with a single energy barrier corresponds to kinetics in which the integrin stays in one stable conformation for a period of time until rapidly transitioning to the other state−the top of the energy barrier corresponds to the transition state across which the molecule should spend virtually no time jumping.This agrees with the rapid transitions between the bent and extended integrin α 5 β 1 , but contradicts with the much slower conformational changes of integrin α V β 3 (long ⟨t sw± ⟩) observed in our experiments (Figure 5E).Remembering our MD simulation where ∼7 Hbonds holding the integrin in the bent conformation were sequentially disrupted over a long distance traversed by the integrin α V β 3 headpiece during its unbending (Figure 6G; Supp. Figure 3E−I), we reason that the gradual formation and disruption of H-bonds must involve energy release and absorption, respectively, such that each H-bond would create an energy barrier in the energy landscape along the pathway of conformational change.Between the sequential disruptions of two successive H-bonds, integrin α V β 3 would stay for some time in an energy well separated by the two energy barriers, i.e., a metastable state with intermediate energy.This is in sharp contrast to integrin α 5 β 1 because the ∼2 H-bonds between the headpiece and tailpiece of integrin α 5 β 1 were disrupted nearly simultaneously over a much shorter distance over its unbending course (Figure 6C; Supp. Figure 3A−D), likely allowing their corresponding energy barriers to merge into one, which enables us to model its energy landscape by that depicted in Figure 7A.
We thus constructed an energy landscape model of integrin α V β 3 with 7 energy barriers serially distributed between the bent and extended states, thereby creating 8 conformational states (one bent, 6 intermediate and one extended) (Figure 7B).In our cell-free system, the only energy source that drives the conformational transitions is microscopic thermal agitations from the macroscopically thermodynamically equilibrated environment.Since the purified protein may have no mechanism to regulate the directional tendency of conformational changes, the integrin that resides in any intermediate state could transition bidirectionally toward either bending or unbending regardless of the previous direction of its immediate past transition, i.e., the molecule may reversibly transition back-and-forth between any two adjacent states before jumping over the last energy barrier to one of the observed stable states (eq 4 and Supp.Equation 5), giving rise to the slow bending and unbending dynamics observed in our experiment.The switching time ⟨t sw± ⟩ is thus broken down into the unmeasurable times for the integrin to hop over the intermediate energy barriers and the measurable times for it to park in the intermediate energy wells before transitioning over to the next energy barrier.For the sake of simplicity, we further assumed that all energy barriers in integrin α V β 3 between the intermediate states were identical in shape and evenly distributed between the bent and extended conformations, hence having identical transition rates between any two adjacent intermediate states: k − and k + .The respective rates of transition from the bent or extended state to their adjacent intermediate states were designated as k + Bent and k − Extended respectively.By treating the stochastic conformational changes as a Markov process in a finite state space, including bent, intermediate, and extended states, we built a master equation: , where S is the vector of probabilities of the molecule to assume any of the states.T is a [N + 2]-by-[N + 2] matrix of transition rates, in which N = 6 is the number of intermediate states (Supp.Equation 6).Using the probability vector solved from the master equation, we express the average time-totransition ⟨t 0± ⟩ and switch time ⟨t sw± ⟩ in terms of the kinetic rates (see Supp.Methods for details): (1 ) where . Assuming that the transition between every two adjacent states follows the Bell model, 39 all transition rates in eq 4 are regulated by force: where k| f=0 is the value of k under zero force.|Δx ± | is the distance from the bottom of the energy well of any intermediate state to the bottom of its adjacent energy well in the energy landscape that takes the positive sign for unbending and the negative sign for bending.Using this model, we fitted the experimental ⟨t sw± ⟩ and ⟨t 0± ⟩ vs force data simultaneously for both cation conditions (Figure 5E,F,H), showing good agreement.Fitting returned two sets of best-fit parameters, one for each cation condition, which allowed us to evaluate the parameters of the energy landscape, including differences between neighboring states: /k + n ) (n = 0 (Bent), 1, 2, ... 6) ("n + 1 7" represents "Extended" state), and plot the energy landscape of α V β 3 conformational changes (Supp.Table 3, Figure 7B).As a sanity check, for both Ca 2+ /Mg 2+ and Mn 2+ cases we calculate the sum of these parameters, finding To further validate our model, we used Monte Carlo simulations to perform "mock runs" based on this energy landscape, which was able to recreate integrin spontaneous bending and unbending conformational changes over time (Supp.Figure 4A,B, Supp.Videos 4−6).Our multistate model predicts that integrin may jump back-and-forth between adjacent states.Indeed, we observed that integrins occasionally paused in the middle of a bending process and reversed the course to unbend in both Monte Carlo simulations (Supp.Figure 4A, Supp.Video 4) and BFP experiments (e.g., Supp. Figure 4C).These results validate our proposed energy landscape of integrin α V β 3 conformational changes.With force applied to integrin α V β 3 , the energy landscape is tilted and the integrin is shifted toward the extended state (Supp.Video 5, 6)  Although it is still not clear how integrin α V β 3 conformational changes can persist under force conditions that are energetically unfavorable, our model seems to suggest a facilitating mechanism: the sequential formation and disruption of H-bonds serve as "stairs" for integrin α V β 3 to temporarily "rest" as it moves up-and down-stairs, so that the energy differential required in each "step" is reduced.Meanwhile, we would like to point out that the model still has limitations.For instance, it assumes a one-dimensional reaction coordinate and that all energy barriers are identical in shape and evenly distributed, and considers only H-bonds but not other types of noncovalent interactions (e.g., salt bridge, hydrophobic interaction) and covalent bond reactions (e.g., thiol−disulfide exchange), which may be addressed in future studies.

CONCLUSIONS
Force-modulated integrin bending and unbending conformational changes have previously been observed on cell surfaces. 11,12Here, we provide real-time single-molecule experimental data to show that purified integrin ectodomains are capable of undergoing force-modulated bending and unbending conformational changes independent of the cellular environment.Our results reveal very different biophysical characteristics for the two focal adhesion integrins: α 5 β 1 and α V β 3 .The conformational changes of integrin α V β 3 are more gradually modulated by force in an "analogous" fashion as opposed to the "digital" fashion seen in the integrin α 5 β 1 case.It is reasonable to speculate that differences in mechanosensitivity generally exist across all integrin species, which directly affects how different integrins interpret and react to the biomechanical environment.Accordingly, different biomechanical features (e.g., elasticity, viscosity, and surface fluidity) probably should be adopted when designing therapeutical nanoparticles and nanomaterials that target distinctive integrin species to achieve optimal accommodation and avoid undesired cell mechanosignaling.
Among the many macromolecular systems that were found to possess the capability of force-modulated reversible conformational transitions, 46−51 integrin α V β 3 appears to be the only one identified so far that is capable of slow-kinetic sizable spontaneous conformational changes under a wide range of force without external energy source.We cannot help to speculate that more molecules with a similar attribute exist and await to be discovered.Studying these mechanosensitive structures will help us understand how they accumulate and convert small-scale thermal energy into the work required for large-scale molecular conformational changes against force.A nanoscopic module that can fulfill such a task should be of potential use in biomaterial-based nanorobots for certain movement tasks (e.g., "switch" and "hinge" movement). 52In this context, the present work provided not only an actual example but also a critical concept and useful design principles for the engineering of protein biomechanical machines in the field of bionanotechnology. 53
AFM Setup, Preparation, and Experiment.Our AFM was built and calibrated in-house. 16A Petri-dish was directly mounted onto a piezo (P-363, Physik Instrumente, Karlsrube Germany), which was controlled by a computer program (Labview, National Instruments) with a subnanometer spatial resolution through capacitive sensor feedback.A laser (Oz Optics, Ontario, Canada) was focused on the back of the cantilever (TM microscopes, Sunnyvale, CA) end, and deflected onto a photodiode (Hamamatsu, Bridgewater, NJ) to allow the cantilever deflection to be converted to force based on the cantilever spring constant. 56To engage the integrins with FN, cantilever tips were incubated with 10−20 μg/mL FN overnight at 4 °C, rinsed, and incubated in Tris-buffered saline (50 mM Tris-Cl, 150 mM NaCl, pH 7.5) containing 1% bovine serine albumin (BSA) for 15 min at room temperature to block nonspecific binding. 16For integrin coating, anti-Penta-His antibody was adsorbed on the Petridish, rinsed, and then incubated with 10 μg/mL α V β 3 -Hexa-His, or GG-7 was adsorbed on the Petri-dish, rinsed, and incubated with 10 μg/mL α 5 β 1 -Fc or trα 5 β 1 -Fc for 30 min.Control experiments have been performed in a previous work (using the same instrumental setup and molecular systems) 16 and in the present work, which ensured that the detected binding events were mostly mediated by specific binding between the integrins and FN, while nonspecific binding events were negligible. 16ome of the AFM experiment procedures have been described previously. 16,36Briefly, the Petri-dish was added with a buffer of the desired cation composition.The piezo brought the Petri-dish to contact the cantilever tip, retracted slightly and held the Petri-dish close to the tip for 0.5 s to allow bond formation, and then retracted it at a speed of 200 nm/s.The presence of an adhesion event was reflected by a positive force signal in the force-time curves.The coating of the Petri-dish was titrated to keep adhesion infrequent (<20%), a necessary condition for most of the adhesion events (>89%) to be mediated by single bonds. 22For force-induced unbending and rebending measurements, the Petri-dish was driven at a constant speed (200 nm/s) to load the bond to ∼20 pN and retract at the same speed to unload the bond.The (un)bending events were identified and parameters measured from the force-time traces (cf. Figure 1B).For CMR measurements, the Petri-dish was driven to move cyclically so the integrin−FN bond underwent force loading and unloading and then held at a preset force (cf. Figure 4A,C). 36ifetime was measured from the instant when the force reached the desired level to the instant of bond dissociation.The collected lifetime data were categorized into bins of successive force ranges, and averaged within each force bin to plot the lifetime curve.For forceramp after a cyclic loading−unloading cycle with a high peak force, the piezo was retracted at a very low speed (1 nm/s) to allow observation of repetitive unbending and bending events over a prolonged period until bond rupture.
RBC and Glass Bead Preparation.Human blood (8−10 μL) was obtained from finger prick following a protocol approved by the Institutional Review Board of Georgia Institute of Technology (protocol number H12354) and The University of Texas Medical Branch (protocol number 22-0015).RBCs were isolated and biotinylated by incubating with Biotin-PEG3500-NHS solution. 11he biotinylated RBCs were then incubated with nystatin, which would swell the RBCs to near spherical shapes.
The procedure for bead functionalization has been described. 57riefly, after thiolation, glass beads were incubated with streptavidinmaleimide, anti-Penta-His antibody cross-linked with MAL-PEG3500-NHS, or LIBS-2 cross-linked with MAL-PEG3500-NHS overnight.Streptavidin-coated beads were incubated with biotinylated FN solution for 2 h.Anti-Penta-His antibody coated beads were incubated with α V β 3 -Hexa-His or α 5 β 1 -Poly-His solution for 3 h.LIBS-2 coated beads were used without further incubation.All beads after incubation were washed with and resuspended in phosphate buffer (27.6 g/L NaH 2 PO 4 •H 2 O, 28.4 g/L Na 2 HPO 4 ).
Platelet Isolation.The procedure for collecting human venous blood was approved by the Institutional Review Board of the Georgia Institute of Technology (protocol number H12354).Blood was collected from healthy volunteers into tubes containing anticoagulant and activation-suppressing agents, and centrifuged at 200g for 15 min to isolate platelet rich plasma, which was centrifuged at 900g for another 10 min to isolate the platelet pellet.The platelet pellet was resuspended in a platelet washing buffer (4.3 mM K 2 HPO 4 , 4.3 mM Na 2 HPO 4 , 24.3 mM NaH 2 PO 4 , 113 mM NaCl, 5.5 mM D-glucose, 10 mM theophylline, 20 U/mL clexane, 0.01 U/mL apyrase, 1% BSA, pH 6.5), rested for 15 min, and centrifuged again.Finally, the platelet pellet was resuspended into a HEPES-Tyrode buffer (134 mM NaCl, 12 mM NaHCO 3 , 2.9 mM KCl, 0.34 mM sodium phosphate monobasic, 5 mM HEPES, and 5 mM glucose, 0.02 U/mL apyrase, 1% BSA, pH 7.4) ready for experiments.
BFP Setup, Preparation, and Experiment.Our BFP apparatus has been described previously. 12,57A chamber mounted on an inverted microscope (Nikon TiE, Nikon) was filled with an experimental buffer supplemented with 1% BSA to block nonspecific binding and cations (1 mM Ca 2+ /Mg 2+ or 2 mM Mn 2+ ).A biotinylated RBC was aspirated by a micropipette to act as a force transducer (Figures 2A and 3E, left), the spring constant of which was set to 0.5 pN/nm when assessing integrin α 5 β 1 , and to 0.25 or 0.3 pN/nm when assessing integrin α V β 3.
11 A probe bead bearing FN or LIBS-2 was attached to the apex of the RBC via streptavidin−biotin interaction.An integrin α V β 3 -functionalized bead or a platelet was aspirated by an opposing micropipette (Figures 2A and 3E, right) as the target, and driven by a piezoelectric translator (Physical Instrument) to repeatedly touch with the probe bead and retract.The probe bead's position was tracked by a high-speed camera.
The BFP measurement procedures for bond lifetime, (un)bending, and CMR are similar to those for AFM experiments, wherein a tensile force signal indicated an adhesion event between the probe bead and the target.FN coating on the probe bead was titrated to maintain infrequent adhesion (<20%). 22For integrin α 5 β 1 experiments, all adhesion bonds were ramped until they broke.For integrin α V β 3 experiments, upon the detection of an adhesion event, the target pipet was held at a desired position (reflected by the initial clamping force) to wait for the bond to dissociate.
Molecular Stiffness Measurement.As previously described, 11,12 force vs time data from AFM and BFP experiments were transformed to "force vs. extension" data (cf.from Figure 1B to 1E).The tensile force portion of the "force vs. extension" data was fitted by a line and the slope was taken as the stiffness of the integrin-FN complex.The value mainly reflects the integrin stiffness as the contribution from FN is negligible.
Molecular Dynamics (MD) Simulations of Integrins α 5 β 1 and α V β 3 .The ectodomain crystal structure of integrins α 5 β 1 (PDB code 7NXD) 21 and α V β 3 (PDB code 3IJE) 58 was used to perform the MD simulation with GROMACS. 59The TIP3P model was used to depict water molecules.Na + and Cl − were added to neutralize the system and maintain the physiological salt condition (150 mM).The CHARM36 force field 60 was used to describe the interactions of the protein and the solvent.CHARMM Additive All-Atom Force Field 61 was used to describe the sugar.Simulations began with minimizing the energy of the protein using steep decent methods, and then the system temperature was raised from 3 to 300 K in an annealing simulation with controlled volume within 500 ps, followed by another 500 ps simulation in NVT ensemble.Afterward, a 1-ns simulation was performed in an NPT ensemble at 300 K and 1 atm.The temperature and pressure were controlled by a V-rescale thermostat and Parrinello−Rahman barostat, respectively. 62In the annealing, NVT, and NPT simulations, the positions of the heavy atoms of the integrin were restrained.
In the MD simulations of the partially and fully extended integrin structures acquired from the above SMD simulations, the C-terminal Cα atom of β tail was restrained, and the clamping force was applied to the same group of atoms as in the SMD simulation, with the pulling speed set to 0. On the other hand, for MD simulation of the bent α 5 β 1 and α V β 3 integrins, no restraint was applied.The numbers of H-bonds in all the above bent, partially extended, and fully extended integrin structures were analyzed with a threshold distance of 0.3 nm and a donor−acceptor angle of 20°.
Statistical Analysis.Statistical significance was assessed by unpaired or paired, two-tailed Student's t test or one-way ANOVA.

Figure 1 .
Figure 1.Observing and characterizing force-modulated integrin α 5 β 1 unbending and bending.A. Superimposition of AFM experimental setups for integrins α 5 β 1 and α V β 3 .Recombinant integrin α 5 β 1 , truncated integrin α 5 β 1 , and integrin α V β 3 were respectively immobilized on a polystyrene surface using mAbs GG7 (anti-Fc) or anti-Hexa-Histidine.Here and in all following figures, integrins α 5 β 1 and α V β 3 are respectively colored by magenta and yellow.B. A representative AFM force vs time trace of a loading−unloading cycle on an integrin α 5 β 1 − FN bond.Two "kinks", one in the loading and the other in the unloading phase, respectively represent integrin unbending and bending.C. Mean ± standard error of the probability of observing structural changes in integrin α 5 β 1 or trα 5 β 1 in force loading−unloading processes in different metal ion conditions.D. Distribution of the difference of the integrin α 5 β 1 −FN complex molecular length before and after a full force loading−unloading cycle.E,F.Representative force vs extension curves of loading (red) and unloading (blue) in Ca 2+ (E) and Mn 2+ (F).The loading and unloading traces were linearly fitted (black dashed lines) to evaluate molecular stiffness, which shows in (E) two distictive stiffness values (Slope 1 , Slope 2 ) exist for the bent and extended integrin α 5 β 1 , respectively.G-I.Histograms of the integrin α 5 β 1 −FN complex stiffness before (G) and after (H) unbending in Ca 2+ , and with no visible structural change in Mn 2+ (I), and their respective Gaussian distribution fits (mean and standard deviation (SD) annotated).J. Histogram of AFM-measured integrin α 5 β 1 head-to-tail molecular extension change due to unbending in Ca 2+ .K. BFP photomicrograph.L. Mean ± standard error of adhesion frequency (lef t) and α 5 β 1 unbending probability (right) in Ca 2+ and Mn 2+ in BFP assay.The integrin α 5 β 1 -blocking monoclonal antibody (mAb), BMC5 eliminated most adhesion in Ca 2+ .M. Representative BFP force vs time trace of a force ramp cycle on an integrin α 5 β 1 −FN bond.An unbending event is highlighted in the red circle.N. Histogram of BFP-measured α 5 β 1 head-to-tail molecular extension change due to unbending in Ca 2+ .

Figure 2 .
Figure 2. The observation of force-regulated integrin α V β 3 unbending and bending by BFP. A. The adhesion frequency of integrin α V β 3 −FN binding in Mn 2+ and Ca 2+ /Mg 2+ conditions.The addition of mAb LM609 blocked most of the adhesion events in Ca 2+ /Mg 2+ .B. Mean ± s.e.m. of lifetime vs force of single integrin α V β 3 −FN bonds in indicated conditions or integrin α V β 3 − LIBS-2 bonds.C−F.Representative BFP force vs time (C,E) and displacement vs time (D,F) traces respectively showing an integrin unbending (C,D) and bending (E,F) event in the position-clamp phase, along with cartoons depicting different integrin α V β 3 conformations before and after (un)bending.Panels D and F are respectively converted from Panels C and E, where BFP displacement is calculated as Force/k RBC (RBC spring constant).The data (points) is smoothened using the Savitzky-Golay method (curves) to obtain a higher force resolution.Inserts in panels D and F: detailed views of the conformational changes within the cyan-shaded windows that convert the BFP displacement to the integrin α V β 3 extension change, with standard deviations of the signals, σ, indicated as a measure of thermal fluctuation before and after the (un)bending.Definitions of time-to-switch and switching time are indicated.

Figure 3 .
Figure 3. Characterization of force-regulated integrin α V β 3 unbending and bending by BFP.A,B.Scatter plots, histograms (bars) and Gaussian fits (curves) of integrin α V β 3 extension changes due to unbending (lef t) and bending (right) in Mn 2+ (A) and Ca 2+ /Mg 2+ (B).The two solid dots in (A) respectively correspond to the representative unbending and bending events depicted in Figures 2D,E and 2F,G.C. Data (points) and the median and 5−95 percentiles (box and whisker) of integrin α V β 3 extension changes due to unbending and bending in Mn 2+ and Ca 2+ /Mg 2+ .D. Data (points) and the median and 5−95 percentiles (box and whisker) of the integrin α V β 3 −FN molecular stiffness before unbending events and before bending events.N.S. = not significant; ** p < 0.01; **** p < 0.0001, assessed by one-way ANOVA.E. Fitting the integrin α V β 3 −FN molecular stiffness before unbending and before bending with Gaussian distribution to respectively acquire the average molecular stiffness associated with bent and extended integrins.F,G.Fitting the integrin α V β 3 / FN molecular stiffness in Mn 2+ (F) and Ca 2+ /Mg 2+ (G) with dual-Gaussian distribution to calculate the proportions of BFP-detected integrins in bent and extended conformations.The means of the two Gaussian distributions, respectively associated with bent and extended integrin α V β 3 , were derived from (E).

Figure 4 .
Figure 4. Measuring cyclic mechanical reinforcement (CMR) of integrin α V β 3 using AFM and BFP.In both systems the ligand coating was titrated to reach infrequent adhesion (∼20%), a necessary condition for most adhesion events to be mediated by single bonds.A. A representative AFM force vs time trace showing a CMR with one loading−unloading cycle with a ∼20 pN peak force followed by bond lifetime measurement at ∼5 pN, which was used to generate the data in the right group of panel (B).The cartoons indicated how the cantilever would be bent in different segments of the data curve.B. Data (points) and the median and 5-95 percentiles (box and whisker) of integrin α V β 3 −FN bond lifetimes measured after a singlecycle CMR (red, exemplified in panel (A)) or without CMR (black, exemplified in panel (C), top).C. Two representative AFM force vs time traces showing integrin α V β 3 lifetime measurements of a bond with 0.5 (top) and 1.5 (bottom) loading−unloading cycle before clamping at the peak force, which were used to generate the data in the first two groups in (D).D. Data (points) and the median and 5−95 percentiles (box and whisker) of integrin α V β 3 −FN bond lifetimes measured after the indicated numbers of CMR cycles.E. BFP photomicrograph showing the experiment setup used to generate the data in (F), where a platelet aspirated by an opposing micropipette acted as the target.F. Data (points) and the median and 5−95 percentiles (box and whisker) of platelet integrin α V β 3 −FN bond lifetimes measured after the indicated numbers of CMR cycles using the BFP shown in (E).

Figure 5 .
Figure 5. Force-modulated integrin α 5 β 1 and α V β 3 bending and unbending kinetics.A. A representative force vs time trace of applying slow ramping force on an integrin α 5 β 1 −FN bond after a single CMR cycle, which was measured by AFM in Ca 2+ to exemplify reversible and consecutive unbending−bending events of integrin α 5 β 1 .Insert: zoom-in of the curve segment showing repeated bending-unbending events in a narrow force range near ∼7 pN.B. Semilog plots of mean ± s.e.m., integrin α 5 β 1 time-to-unbending t 0+ (square) and time-to-bending t 0− (triangle) vs force data and their fits by the Bell model (curves).The two fitting curves intersect at 7.4 ± 0.6 pN and 0.076 ± 0.017 s (arrow).C. Cumulative histogram of integrin α 5 β 1 unbending force distribution.The distribution was fitted by a theoretical model to derive the parameters of the energy landscape.The equation of the model and the derived parameters were denoted.D. Plot of ⟨t 0− ⟩ to ⟨t 0+ ⟩ ratio of integrin α 5 β 1 conformational changes vs force, calculated based on experimental data (point) and the model fitting in panel (C) (curve).E,F.Semilog plots of mean ± s.e.m., integrin α V β 3 unbending time t sw+ (E) or time-to-unbending t 0 + (F) (hollow triangle and hollow inverted triangle) and bending time t sw− (E) or time-to-bending t 0− (F) (hollow square and hollow circle) vs force data measured in the indicated cation conditions, and their theoretical fits by the multistate model described in the text.The R 2 values of the fittings are 0.95 and 0.96 for Mn 2+ and Ca 2+ /Mg 2+ conditions, respectively.Solid dots: mean ± s.e.m. t 0± and t sw± vs force of cell surface integrin α V β 3 unbending (light magenta) and bending (light cyan) events in Ca 2+ /Mg 2+ .G. Cumulative histogram of integrin α V β 3 unbending force distribution with theoretical model fitting.H. Plots of ⟨t 0− ⟩ to ⟨t 0+ ⟩ ratio of integrin α V β 3 vs force measured under indicated cation conditions and their model fits.

Figure 6 .
Figure 6.MD simulations of integrins α 5 β 1 and α V β 3 unbending conformational change.A,E.Snapshots of representative integrin α 5 β 1 (A) and α V β 3 (E) conformations (bent, 2 or 3 intermediate, and extended) observed from the MD simulations in which the most commonly observed H-bonds are indicated by their donor and acceptor residues.The red number on top of each panel represents the average number of H-bonds observed in 5 independent MD simulations.B,F.Cartoons depicting the average numbers and locations of H-bonds in relation to the head-to-tail distances during integrins α 5 β 1 (B) and α V β 3 (F) unbending.C,G.Change of average number of H-bonds (mean ± s.e.m., from 5 independent runs of MD simulations) between integrin headpiece and tailpiece during integrin α 5 β 1 (C) and α V β 3 (G) unbending.Linear fitting was applied to the first two points in (C) and all points in (G) to estimate the speed of H-bond breakage, as reflected by the slope.D,H.Average occupancy of the most frequently formed 7 H-bonds in a bent integrin α 5 β 1 (D) or the most frequently formed 8 Hbonds in a bent α V β 3 (H) when the integrin unbents to certain head-to-tail distances.Amino acids in integrin α and β chains are shown in red and black, respectively.

Figure 7 .
Figure 7. Energy landscapes of integrins α 5 β 1 and α V β 3 bending and unbending conformational changes.(A) Energy landscapes of integrin α 5 β 1 ectodomain conformation under zero force (dark curve) and f 1/2 (light curve) based on the experimental and model-fit parameters (Figure 5B).(B) Energy landscape of integrin α V β 3 ectodomain conformation in Ca 2+ /Mg 2+ under zero force based on the experimental and model-fit parameters (Figure 5E,F).Energy wells corresponding to the bent (magenta), intermediate (green), and extended (cyan) states are marked by different colors.

m o o o o o o o o o o o o o o o o o o o o o o n o o o o o o o o o o o o o o o o o o o o
corresponding to f 1/2 = ΔG/Δx = 6.14 and 2.64 pN, respectively, which are comparable to the values listed in Supp.Table2, validating that the serial energy barrier model is equivalent to the single energy barrier model in terms of both energetics and force-dependency.