Force-Dependent Binding Constants

Life is an emergent property of transient interactions between biomolecules and other organic and inorganic molecules that somehow leads to harmony and order. Measurement and quantitation of these biological interactions are of value to scientists and are major goals of biochemistry, as affinities provide insight into biological processes. In an organism, these interactions occur in the context of forces and the need for a consideration of binding affinities in the context of a changing mechanical landscape necessitates a new way to consider the biochemistry of protein–protein interactions. In the past few decades, the field of mechanobiology has exploded, as both the appreciation of, and the technical advances required to facilitate the study of, how forces impact biological processes have become evident. The aim of this review is to introduce the concept of force dependence of biomolecular interactions and the requirement to be able to measure force-dependent binding constants. The focus of this discussion will be on the mechanotransduction that occurs at the integrin-mediated adhesions with the extracellular matrix and the major mechanosensors talin and vinculin. However, the approaches that the cell uses to sense and respond to forces can be applied to other systems, and this therefore provides a general discussion of the force dependence of biomolecule interactions.

T he development of life has evolved in the context of physical forces acting on biological systems. From individual molecules to organelles to cells, tissues, and organs, every part of every organism is exposed to and experiences forces. These forces, generated or experienced, impact every aspect of physiology; 1−3 every cell interprets "classical" signaling pathways (growth factors, hormones, etc.) in the context of its physical environment. 5−8 The purpose of this review is to provide a brief introduction to the concept of force-dependent binding constants, and we will introduce the study of how forces can impact biomolecular interactions. This review will be divided into three sections. The first section will be a general discussion of biomolecular interactions and their importance in biological processes; in particular, this section will focus on the protein interactions involved in mechanotransduction leading to the appreciation that many of these protein interactions have a force-dependent component. In the second section, we will discuss the modifications to the theory of binding constants required to enable force dependence to be considered. Finally, in the third section, we will discuss some of the novel approaches that are emerging and/or required to enable force-dependent binding constants to be measured.

■ BIOMOLECULAR INTERACTIONS
For the purposes of this review, we focus on interactions with proteins and DNA, although the concepts and principles can be applied to other systems. Interactions between proteins and other proteins, DNA, lipid membranes, inorganic metal ions, etc., are mediated by compatible interacting residues and surfaces, i.e., a surface on the substrate protein that has the optimal shape to recognize its ligand (a moiety that forms a complex with that biomolecule to serve a biologically relevant purpose). The ligand can be any biomolecule or non-organic molecule that interacts with the biomolecule in a meaningful way. Interacting surfaces that have important biological functions tend to be highly conserved through evolution, to preserve and maintain the interaction.
The binding affinity of an interaction describes the strength of the binding between a target molecule and its ligand. This binding affinity is usually reported as the equilibrium dissociation constant, K d . This quantity is defined as the ratio between the off-rate, k off , typically in units of s −1 , and the on-rate, k on , typically in units of M −1 s −1 . Therefore, the dissociation constant, K k k d off on = , has a dimension of concentration typically expressed in molar concentration. K d can also be expressed by the equilibrium ratio of the fractions of the bound target (α on ) and that of the unbound target (α off = 1 − α on ) molecules, affinity K d with a single-molecule level of accuracy. Such quantification of K d is mainly achieved either by measuring the on-(k on ) and off-rates (k off ) through the equation and Δg 0 is the free energy difference between the unbound and bound states, which is related to the dissociation constant by the equation K d = ce −βΔg 0 . Measurement and quantitation of biological interactions are of value to scientists and are major goals of biochemistry, as affinities provide insight into biological processes. In drug discovery, measuring affinities is important to aid the design of drugs with a higher affinity for the target and thus a higher efficacy. If there are two potential ligands available to bind to a single target molecule, at similar concentrations, then the relative affinities of each for the biomolecules will dictate which binds preferentially. If one has a higher affinity (lower K d ), then it will bind preferentially. Modulation of the affinities for the two ligands (via post-translational modification, including alteration of conformation via mechanical force) can alter the complexes that form.

■ FORCES IN BIOLOGY
Life that lives under the sea has forces acting on it that are very different from those on land or the forces on those that take to the skies, and changes in mechanical signaling that arise from these different physical environments enable the stunningly beautiful diversity of creatures. Strikingly, despite this incredible diversity, the adhesive structures holding the cells in place, via contact with neighboring cells and the surrounding extracellular matrix, are largely made of the same building blocks. The appreciation of forces in biology has been the subject of many excellent reviews, 10−13 and these all provide excellent accounts to which we refer the reader. The focus here is how these forces can be sensed by the cell and how they can alter signaling outcomes.

■ FORCES ON BIOMOLECULES
Many forces exist in cells, arising from collisions, flow (both retrograde flow of proteins inside the cell and flow of blood past proteins on the surface of cells), force generation machinery and motor proteins (myosins, kinesins, etc.), and forces exerted from the outside world, gravity, pressure, friction, etc. 14−16 All of these forces are sensed by mechanosensors in the cell and used to control cell behavior. 17−19 As these forces are ubiquitous, it seems safe to assume that many biomolecules in the cell will experience forces and as such that the binding constants of interactions involving these molecules will have a force-dependent component.
Actomyosin Contraction. Motor proteins harness the energy from ATP hydrolysis to generate mechanical energy that drives conformational changes that act on other cellular structures (reviewed in refs 16 and 20). In the case of myosin, its interaction against actin filaments enables motion, with the myosin either moving along the filament or pulling the actin filament toward it. 21,22 These forces generated by myosin motors pulling on the actin filament is dubbed actomyosin contraction and is a major method the cell utilizes to generate and maintain forces. Much of the discussion herein will focus on the mechanisms with which the cell responds to this internally generated force.
For a force to act on a biomolecule, it requires the force generation machinery to couple to the protein. This coupling can be either direct, if the protein directly binds to the actin filament, or indirect, if the actomyosin pulls on another protein that acts on the protein. If the biomolecule is bound just to the force generator, then this will pull the biomolecule toward the force. Trafficking of proteins can occur in this manner; for instance, MyosinX can couple to cargo proteins and drag them along actin tracks to the tips of filopodia. 23,24 In this scenario, the biomolecule is dragged along (and will experience the forces associated with drag).
However, if a protein couples to the force generation machinery but is also tethered to a second less mobile system, then the forces "pull" on the tethered protein and the forces are exerted on all proteins constituting the force-transmission molecular linkage. An example of this tethered system is seen for talin, bound to the integrin−ECM complexes at the plasma membrane 18,25 (Figure 1). Here, when actomyosin contractions act on talin, as the talin is tethered, the forces are exerted on the length of the molecule, and as these forces exceed force thresholds of stability for any talin domains, they can trigger domain unfolding (see the next section for a discussion of the consequences of domain unfolding).
This action of actomyosin contraction pulling on a tethered protein is not exclusive to adhesion; any protein, DNA, or membrane that experiences forces has the potential to undergo changes in its shape and thus its function. Therefore, any consideration of the mechanobiology of biomolecules needs to consider where the force is originating, how it is acting on the biomolecule, and how the biomolecule will react, which will depend on its mobility (can it move, or is it tethered?).
Structural Mechanobiology. One exciting, though not altogether unsurprising, aspect of the mechanosensitive events identified to date is that knowledge of the structural basis of the interaction provides an atomic basis of the mechanosensitive mechanism of the interaction. For instance, in the case of talin it is possible to identify the exact amino acids that render the molecule mechanically sensitive. A striking example of this is the R3 domain of talin; 25−27 this domain has reduced mechanical stability due to a cluster of four threonine residues in the central hydrophobic core, which due to their more polar nature destabilize the domain. Using this precise structural information, it is possible to modulate mechanosensing via targeted point mutations that alter mechanosensing. In the case of R3, modification of the four threonine residues to isoleucine and valine residues (a so-called "IVVI mutant") results in stabilization of the R3 domain, 25−27 shifting its mechanical stability to a higher level (i.e., the modified domain unfolds at a higher force at the same force loading rate). Therefore, a comprehensive appreciation of mechanobiology at the atomic level requires atomic-resolution structural biology.

■ TYPES OF FORCE-DEPENDENT BINDING
There are many ways that force can impact binding, and these effects can quickly stack to give diverse responses. Here we offer a non-exhaustive summary of some of the common force dependencies. Protein−protein interactions can change their affinity by orders of magnitude under different force constraints, and as a result, biochemistry done in bulk solution in vitro captures only part of the picture and lacks consideration of the mechanical regulation of the interaction. As a consequence, it is necessary to consider the force dependence of biomolecular interactions.
The focus of this review is on equilibrium binding constants under mechanical force, K d (F). Discussion of the K d (F) will facilitate the description of forces impacting the affinity of binary interactions, autoinhibition, and the interplay of these factors. This will enable the description of two common forcedependent processes that regulate mechanosensitivity through talin, namely, exposure of cryptic binding sites, whereby force exposes hidden binding sites, a type of autoinhibition, and disruption of binding sites, where the binding site is accessible to the ligand, but force results in the domain unfolding and destroying the binding site. The talin rod contains 13 rod domains, R1−R13, which all combine both of these processes together, to create a series of mechanochemical switches 25 ( Figure 1).
Mechanochemical Switches. Tension-sensitive conformational change is a very rapid way for proteins to respond to mechanical force. In many ways, this can be considered as a post-translational modification; force alters the conformation of a protein, and if this change in shape elicits a change in biochemical function, then force can be converted into biological signals. If the conformational change is reversible, that is when force is released the domain reverts back to its low-force condition, then this provides incredible plasticity and can enable rapid and dynamic changes in signaling outcomes.
The theoretical basis of these processes will be discussed in the next section, but to illustrate this concept of forcedependent binding constants, consider the interactions of a talin rod domain (here the example is R3) with three different ligands.
Conventional Protein Interactions. Here, two proteins interact in a "classical" manner, whereby one ligand binds to a binding site on a folded domain (i.e., RIAM binding to the folded talin rod domain R3, as demonstrated in the left panel in Figure 1). In the absence of force, the talin domains are folded, so the interaction can occur. In this scenario, the level of binding is highest at low force and so the K d is lowest. If the domain is unfolded by mechanical force, then the binding surface on the domain is destroyed, so force, above a certain threshold, drives a sharp decrease in the level of binding. Beyond this point, additional force has a weaker effect on the affinity of the interaction as the interaction is already destroyed.
Cryptic Binding Sites That in the Absence of Force Are Inaccessible to Binding. In conventional protein interactions, force results in the loss of ligands binding to folded rod domains; however, eight of the 13 talin rod domains contain vinculin binding sites (VBSs), amphipathic helices in which the vinculin binding epitope is buried inside the domain. At low force, the affinity of the VBS interaction with vinculin is weak (it has a high K d ) as the binding site is not accessible to vinculin. Here, the K d profile with force is different, at low force, the K d is high as binding is not possible, at forces above the unfolding threshold, then the K d is low [exposed VBS bind tightly to the vinculin head, with a nanomolar K d (although as we will see the vinculin itself is also regulated by forces acting on vinculin)]. However, the force dependence is complicated as at high forces the VBS helix can unfold and lead to a loss of vinculin binding.
Cryptic Binding Sites That Are Exposed Only When All Secondary Structure Is Destroyed. While no such ligand has yet been identified for talin, there is also the third scenario in which a fully extended talin polypeptide creates linear epitopes that bind ligands. Here, the affinity will be the highest at high force when the domain is completely unfolded.
Autoinhibition. Autoinhibition is not always considered to be mechanosensitive, but if proteins that are regulated by autoinhibition form part of the mechanical linkages in the cell, then the affinity of the protein for its ligands is directly correlated to the force on the system. The cryptic nature of VBS in talin is also a type of autoinhibition, 28 and talin is further autoinhibited by the molecule folding up into a tight globular compact structure. 29−31 Vinculin is also regulated by a head−tail interaction, where the binding sites for talin and actin are rendered cryptic. 32,33 All of these layers of autoinhibition have a strong force-dependent component, as the affinity of autoinhibition is controlled by whether the protein is under force (if the autoinhibitory domains are held Biochemistry Perspective apart by mechanical force, then the protein is maintained in an activated state). 27,34,35 What is emerging is that autoinhibition of proteins in force-transmitting pathways represents a major force-dependent mechanism to enable mechanotransduction.

■ FORCE-DEPENDENT KINETIC CHANGES IN BINDING
Another key force dependence of interactions arises from scenarios in which forces alter the kinetics of the interactions. As the binding affinity is determined by the ratio of the dissociation rate to the association rate, the influence of force on the binding affinity must be through force-mediated changes in the dissociation and association rates. While a detailed discussion of these is beyond the scope of this review, it is worth mentioning these effects as they impact the mechanical functioning of force-dependent interactions. The force-dependent dissociation rate is of particular interest, as it defines the average lifetime of force-bearing molecular complexes once they are formed. The two most well characterized occurrences of force-dependent dissociation kinetics affecting biomolecular interactions are the slip bond and the catch bond.
Slip Bonds. A slip bond refers to the phenomenon in which the rate of dissociation of a molecular complex increases as the applied force increases. It indicates that a protein interaction is weakened under force, as the force pulls the two interacting components apart.
Catch Bonds. A catch bond refers to an anti-intuitive phenomenon in which the molecular lifetime of an interaction increases as force applied to the interacting molecules increases, 36−39 which plays critical roles in cell−matrix and cell−cell adhesions. 40−42 Interestingly, catch bond kinetics can have a geometric component, where the interaction exhibits catch bond behavior only when force is exerted in a particular pulling geometry. This leads to directionally asymmetric catch bonds as are seen for the interaction between vinculin and Factin. 43 Directional asymmetry can be rationalized at the atomic level by looking at the geometry of the force vectors on the interacting proteins ( Figure 2). Here, the two extremes are "unzipping" and "shearing" geometries. 44−47 The unzipping geometry typically exhibits slip bond kinetics, while in contrast, the shearing geometry often exhibits catch bond kinetics. At the same force, the dissociation rate is often faster for the unzipping force geometry than for the shearing force geometry. We refer the readers to our recent publications 44,46 for the physical principles underlying the effects of pulling geometry on the force-dependent dissociation kinetics.
Intriguingly, many force-bearing protein−protein interactions, such as the integrin−talin connection, the talin−vinculin connection, and the vinculin−actin connection, are under shearing force geometry. 43,48−50 These protein−protein interactions form the interfaces in various force-transmission supramolecular linkages in cells, to enable mechanosensing. Perhaps these linkages evolved to achieve high mechanical stability for their functions through making shearing−force connections. In contrast, both force geometries occur frequently in force-dependent unfolding of protein domains. For example, when force is exerted through the N-and Ctermini of proteins, Ig domains and α-helix bundles consisting of an odd number of β-strands/α-helices are typically subjected to shearing force geometry, while domains with an even number of β-strands/α-helices are under the unzipping force geometry, as illustrated in Figure 2B.
While important to our understanding of mechanobiology, these two force-dependent kinetic phenomena will be the subject of a subsequent review and will not be considered further here.
The theoretical description of the dissociation constant defined in the previous section can be extended to include the effects of force dependence.
Two-State Binary Interactions. The force-dependent affinity of binary interactions has been discussed previously for simple two-state interactions, 51 where a molecule can exist in either an unbound state or a bound state. When force is applied to this target molecule, each state is associated with a force-induced free energy, which is additional to the free energy change associated with its ligand. Therefore, the applied force can cause an additional change in the free energy from the bound to the unbound states, Δg(F) = Δg 0 + Δϕ(F), where Δg 0 is the binding energy (i.e., the free energy cost of unbinding) at zero force and Δϕ(F) is the force-dependent conformational free energy difference between the unbound and bound states. On the basis of the Boltzmann distribution of the states and the definition of the dissociation constant, e Therefore, the force-dependent affinity of such simple twostate binary interactions is solely determined by Δϕ(F), which can be calculated by the equation is the extension difference between the unbound state and the bound state of the target molecule at the same applied force. 52−54 Over more than a decade of single-molecule manipulation studies, the force−extension curves of many interesting molecules, such as double-stranded DNA (dsDNA), single-stranded DNA (ssDNA), folded protein domains, and unfolded protein peptide chains, have Biochemistry Perspective been investigated. Therefore, the K d (F) of two-state binary interactions can be considered well understood. Figure 3 shows three test-case examples of force-dependent dissociation constants calculated on the basis of the well-characterized force−extension curves of DNA molecules. In each case, a ligand binding causes extension changes of the target DNA molecule at the same applied force ( Figure 3B,E,H), which leads to force-dependent binding affinity K d (F) ( Figure  3C,F,I). In each of these scenarios, the strength of the interaction behaves markedly differently when force is exerted.
Binary Interactions Involving Autoinhibition. Many proteins can adopt multiple conformational states, where the binding sites in these proteins are exposed at different levels of force. Autoinhibition is where the native state (i.e., the conformational state with the lowest energy) results in suppression of the accessibility of the binding site. In this scenario, binding can be discussed on the basis of a three-state model, which involves two unbound states (closed state "off,1" and exposed unbound state "off,2") and one bound state (exposed bound state "on"). In the absence of force, on the basis of the Boltzmann distribution of the three states and the definition of the dissociation constant, is the zeroforce dissociation constant of the ligand molecule binding to the target molecule and K d,o 0 is the zero-force dissociation constant of the exposed binding site in a constitutively open conformation of the target molecule. μ c = g off,2 − g off,1 is the autoinhibition energy, which is the chemical potential energy difference between the open conformation and the closed conformation of the target molecule. This reveals that the dissociation constant approximately increases exponentially as the autoinhibition energy increases.
Autoinhibition can be relieved via a number of mechanisms that reduce the value of μ c , including biochemical processes such as phosphorylation or binding of an activating molecule. 55−58 For force-bearing mechanosensing proteins, mechanical stretching provides another possible means of releasing autoinhibition, which has not been extensively studied in the field. Relief of autoinhibition by mutation is one way to study these processes, as this shifts the autoinhibition dynamics toward a more open conformation, This provides an effective way to study protein dynamics that would normally be observed for the wild-type protein only when it is under force.
Here, we provide a succinct discussion of force-dependent release of autoinhibition and its impact on molecular interactions.
A binary interaction involving autoinhibition can still be understood on the basis of the aforementioned three-state model. The only difference from the zero-force binding case is that each state now contains an additional force-dependent conformational free energy. A similar derivation based on where Δϕ 1,2 (F) = ϕ off,1 (F) − ϕ off,2 (F) and Δϕ on,2 (F) = ϕ on (F) − ϕ off,2 (F). Test case 4: ssDNA Binding to an Autoinhibited Region in a dsDNA Hairpin. We demonstrate the application of this equation using a "simple" model system of autoinhibition. Here, the annealing of a short single-stranded DNA (ssDNA) oligo, acting as the "ligand", binds to a 10-nucleotide complementary region that is buried, cryptic, inside a 20-bp double-stranded DNA (dsDNA) hairpin ( Figure 4A). The dsDNA hairpin can exist in two distinct unbound states: a closed hairpin state (state "off,1" in Figure 4A) and an open unzipped state (state "off,2" in Figure 4A). Considering the bound state (state "on" in Figure 4A) where the ligand ssDNA binds to the complementary region of the dsDNA hairpin, there are in total three states that need to be considered. In this case, eq 1 can be directly applied to calculate the forcedependent dissociation constant of ssDNA binding, where the autoinhibition energy, μ c , is the base pairing energy in the hairpin. Under physiological conditions, a base pair energy is in the range of 1−4 k B T depending on the nearest-neighbor dinucleotide sequences. 59,60 Assuming an average 2 k B T per base pair, for a 20-bp DNA hairpin, μ c is around 40 k B T, which completely inhibits binding of the ligand ssDNA. However, such strong autoinhibition can be easily released by forces. On the basis of μ c = 40 k B T and the force−extension curves of ssDNA and dsDNA, it is found that forces of ∼15 pN can decrease the dissociation constant (i.e., increase the binding affinity) by more than 10 15 fold ( Figure 4B).
Interestingly, the predicted K d (F) has a biphasic force dependence, which can be divided into two regions: a monotonically decreasing function at forces below 15 pN due to the force-dependent release of autoinhibition, maximal binding affinity at 15 pN where the binding region is no longer autoinhibited, and a monotonically increasing function at forces above 15 pN due to force-dependent destabilization of the force-bearing DNA duplex. 52 Even in this simple model system, the force dependence on the K d is complex.
Test Case 5: Vinculin D1 Domain Binding to Talin. Another important example of autoinhibition affecting the K d (F) is the binding of the vinculin D1 domain (Vd1) to a vinculin binding site (VBS) buried in a talin rod α-helical bundle. In the absence of force, the VBS is cryptic in the folded rod domains. In contrast to a dsDNA hairpin that can exist in only two distinct unbound states, an α-helical bundle can exist in three distinct unbound states: an autoinhibited folded state (state "off,1" in Figure 4C), an unfolded state in which the VBS exists in an α-helical conformation (state "off,2" in Figure  4C), and an unfolded state in which the VBS becomes an unstructured peptide polymer (state "off,3" in Figure 4C). Considering the bound state (state "on" in Figure 4C) where Vd1 binds to the α-helical conformation of VBS, there are in total four states that need to be considered. Denoting ε as the chemical potential energy between the unstructured and αhelical conformations of the VBS (broadly equivalent to the stability of one α-helix) and K d,o 0 as the zero-force dissociation constant of Vd1 for the exposed α-helical conformation of VBS, on the basis of a similar analysis of the force-dependent energies of the four states, we can show that where Δϕ 1,2 (F) = ϕ off,1 (F) − ϕ off,2 (F) and Δϕ 3,2 (F) = ϕ off,3 (F) − ϕ off,2 (F), which can be computed on the basis of the force− extension curves of the states. Using a fixed μ c value of 11 k B T, the fold change of K d (F) relative to K d,o 0 can be calculated for several values of ε ( Figure 4D). At these parameter values, a force of ∼5 pN can decrease the dissociation constant (i.e., increase the binding affinity) of Vd1 for the talin rod α-helical bundle by >10000 fold. In other words, the high autoinhibition energy μ c of 11 k B T that limits vinculin binding in the absence of force can be released by a small force of ∼5 pN.
The predicted K d (F) exhibits an overall biphasic profile, which can be divided into three regions. A monotonically decreasing function at forces of <5 pN resulted from forcedependent release of autoinhibition energy (μ c = 11 k B T), an almost force-independent basin region, followed by a monotonically increasing function due to force-dependent destabilization of the α-helical conformation and thus the Biochemistry Perspective bound complex. 27,34,61 In contrast to the sharp switch from the decreasing profile of K d (F) to the increasing profile of K d (F) at ∼15 pN in the case of ssDNA binding to an autoinhibited region in a dsDNA hairpin ( Figure 4B), the switch is much less sharp in the case of Vd1 binding to a VBS-containing domain ( Figure 4D). Here, after the domain unfolds and autoinhibition is relieved, the VBS binding helix is exposed and has maximal affinity for Vd1 all the while the VBS helix is folded. The more stable the VBS helix, the greater the force range over which binding affinity is maximal. As a result, maximal affinity is present over a range of forces, as seen by the forceindependent basin region in Figure 4D with the force range of the basin determined by the stability, ε, of the VBS α-helix.
Mutually Exclusive Binary Interactions Involving Autoinhibition. The previous examples describe the scenarios of the complex force dependence on K d (F) where the applied force can drastically increase the binding affinity between the target molecule and its ligand by releasing the autoinhibition of the target molecule, thus exposing the binding site(s) for its ligand, and where the applied force can also then decrease the binding affinity by destabilizing the conformation of the bound state.
In the case of talin rod domains, these force-dependent components can be multiplexed, to create force-dependent switching of binding partners. The force-dependent switching of binding partners on talin rod domain 3 (R3) provides an example of this (Figure 1). At low forces, the Rap1 effector RIAM binds the α-helical bundle form of R3 where the VBSs are cryptic (Figure 1 and top panel in Figure 5A), whereas at high forces, the autoinhibition of VBSs in R3 is released and exposure of the VBSs significantly increases the binding affinity for Vd1 (Figure 1, Figure 4C,D, and bottom panel in Figure  5A). As such, force drives a change in binding partners.
To be specific, the force-dependent binding constants K d (F) of R3−RIAM interaction and R3−Vd1 interaction can be derived on the basis of a similar analysis of the force-dependent energies of all of the states involved. Via analysis of the four states involved in the R3−RIAM interaction (illustrated in the top panel in Figure 5A) where K d 0 is the zero-force dissociation constant of RIAM for the α-helical bundle form of R3, Δϕ 2,1 (F) = ϕ off,2 (F) − ϕ off,1 (F), and Δϕ 3,1 (F) = ϕ off,3 (F) − ϕ off,1 (F).
With regard to the R3−Vd1 interaction, there are two VBSs in talin R3, and in such a scenario, a complete description needs to consider the multiple bound states and the effect of volume exclusion. For the sake of simplicity of demonstrating the idea of force-dependent switching of binding partners, here we consider only binding to one VBS in R3. As such, there are one bound state and three unbound states (illustrated in the bottom panel of Figure 5A), and force-dependent dissociation constant K d (F) of the R3−Vd1 interaction can be directly calculated by eq 2.
In eqs 2 and 3, μ c is the chemical potential energy between the α-helical bundle form of R3 and its extended α-helix chain, and ε is the chemical potential energy between the unstructured and α-helical conformations of one α-helix.
Using fixed values for μ c of 11 k B T and ε of 5 k B T for the purpose of demonstration, the force-dependent switching between R3−RIAM and R3−Vd1 interactions can be shown in Figure 5B.

■ MEASUREMENT OF K d (F) IN EXPERIMENTS
To further explore the impact of forces on the affinity of binary interactions, direct measurement of K d (F) in experiments will provide a straightforward understanding of the force dependence of binding affinity. This section is devoted to providing a brief discussion of the measurement of K d (F) in experiments, which includes (i) the measurement of zero-force dissociation constant K d (0) = K d 0 by bulk technology and (ii) the measurement of dissociation constant under force K d (F) by single-molecule manipulation technology.
Bulk Technology for Measuring the K d (0). There are many methods for biochemically measuring zero-force dissociation constants, which are generally based on two approaches.
The first way to quantify the K d (0) is based on the measurement of the bound (α on ) or unbound fraction (α off ) of target molecules through K c

Biochemistry
Perspective molecules in the corresponding bands. Such methods require an additional assumption that the bound and unbound fractions of the target molecules remain fixed during the gel shift assay. The other commonly applied way to quantify the binding affinity is based on the measurement of the association (k on ) and dissociation rates (k off ) through K k k d off on = . The surface plasmon resonance (SPR) technology is a representative example, which detects binding of ligands to target molecules tethered on a gold surface based on the binding-induced shift in the resonant oscillation of conduction electrons ( Figure  6A). 64,65 In typical experiments ( Figure 6B), SPR measures the evolution of the resonance signal after flowing a ligandcontaining solution until it reaches equilibrium. On the basis of the most commonly used Langmuir model, 66,67 this time evolution follows the single-exponential relation R(Δt a ) = R eq [1 − e −(ck on +k off )Δt a ], where R eq is the resonance signal when the binding and unbinding of ligands reach equilibrium, c is the ligand concentration, and Δt a is the time duration after flowing in the ligand-containing solution. After the removal of the ligand from the solution, the bound ligands dissociate, resulting in an SPR signal time evolution: R(Δt d ) = R eq e −k off Δt d , where Δt d is the duration after the removal of the ligand from the solution. By fitting the association and dissociation SPR data with the two equations, one can obtain the values of association and dissociation rates and thus K d .
Besides the EMSA and the SPR assay, a number of other methods have been developed to quantify the binding affinity of molecules based on either measurement of the bound and unbound fractions of target molecules or measurement of the kinetic rates, which include (but are not limited to) fluorescence, fluorescence polarization (FP), nuclear magnetic resonance (NMR), isothermal titration calorimetry (ITC), stopped flow kinetics, etc., and these are discussed elsewhere. 68−70 However, these biochemical assays are not applicable for the study of the force dependence of molecular interactions, because they do not apply mechanical constraints to molecules.
Single-Molecule Technology for Measuring K d (F). As mentioned above, single-molecule technologies can also be used to determine the dissociation constant K d either by quantifying the equilibrium binding probability p on through the  74−77 For the measurement of K d (F), single-molecule mechanical manipulation is necessary to apply force to a target molecule and measure its force-dependent interactions with the ligand. Therefore, the following section will focus on the measurement of K d (F) using single-molecule mechanical manipulation technology.
Single-molecule mechanical manipulation technologies 78 can be categorized into two groups on the basis of the types of mechanical constraints they apply to a molecule. In one group represented by optical tweezers (OT) 79 and atomic force microscopy (AFM), 80 an external Hookean spring is attached to one end of the tethered molecule and its distance R from the other end of the molecule is controlled (Figure 7, top panel).
In the other group, represented by magnetic tweezers (MT), 81 centrifuge tweezers, 82 and acoustic tweezers, 83 an external force is applied to a bead attached to the tethered molecule and the level of force is controlled (Figure 7, bottom panel). Through a force-clamping feedback control, OT and AFM can also apply an external force control to molecules. 84, 85 To measure the force-dependent dissociation constant K d (F), it is most convenient to measure the interaction across a range of constant forces.
Because the measurement is performed on a single target molecule, it is desirable to record repetitive binding and unbinding events over a long duration of measurement. This imposes strong requirements on the stability of the instrument over long durations. Magnetic tweezers can make measurements over a time course of hours to days with negligible spatial and force drifts. 86,87 The force-dependent binding to, and unbinding from, the mechanically manipulated target

Biochemistry
Perspective molecule can in principle be detected on the basis of integration of single-molecule fluorescence imaging with single-molecule manipulation. 88 However, such integration has a drawback of photobleaching, which impairs the long duration measurement for interactions with slow kinetics.
Hence, it is important to develop label-free measurement approaches for detection and quantitation of single-molecular interactions under force. In a force−constraint single-molecule manipulation experiment, the target molecule is tethered between a coverslip surface at one end and a bead at the other end. 89,90 A well-controlled stretching force, which can be calibrated at a sub-piconewton resolution, is applied to the target molecule through the bead; the bead height from the surface can be measured at nanometer resolution. Hence, it can detect molecular extension changes of a few nanometers. Utilizing this spatial resolution, label-free measurement of single-molecule interaction under force can be based on detecting (i) binding-induced deformation of the target molecule ( Figure 8A,B) or (ii) binding-induced delay to a structural transition ( Figure 8C,D). These two label-free detection approaches, discussed in the next section, have been applied in a number of recent studies to detect DNA− protein 76,77,91,92 and protein−protein interactions. 27,34 More recently, these detection approaches have been further developed to quantify the dissociation constant K d (F). 76,77 Detection Based on Target Molecule Deformation. Ligand binding can be probed on the basis of the changes in the endto-end extension resulting from binding-induced structural transition or deformation of the binding site ( Figure 8A,B). This approach is suitable for ligands that induce a detectable change in the extension of the tethered molecule. An example of this is the interaction between Vd1 and the mechanically exposed VBS in talin R3 27 and α-catenin. 34 At forces of >15 pN, the exposed VBS exists in a randomly coiled peptide conformation. Binding of Vd1 induces the formation of the αhelical conformation of the VBS, resulting in a detectable stepwise extension decrease of 2−3 nm depending on the applied force. Similarly, when Vd1 dissociates, it will be accompanied by a 2−3 nm stepwise extension increase ( Figure  8A,B). This 2−3 nm step provides a visible and quantifiable readout of ligand binding and unbinding. From the time trace of such two-state stepwise extension fluctuation, the association and dissociation rates can be determined by obtaining the dwell time in bound and unbound states; therefore, the forcedependent dissociation constant can be determined by the Detection Based on Delayed Structural Transitions. A bound ligand on a target molecule can also be detected if binding results in a delayed structural transition. This can be (i) delayed refolding if the ligand is bound on an unfolded structure via an exposed binding site ( Figure 8C) and (ii) delayed unfolding if the ligand is bound on a folded structure ( Figure 8D).
An example of delayed protein refolding due to ligand binding is the interaction between Vd1 and VBS in talin rod and α-catenin domains. 27,34,35 Previous studies from our group have shown that Vd1 bound on the mechanically exposed VBS in the domains can keep the domains in the unfolded conformation for a longer duration after force is released, compared with in the absence of Vd1. This resulting longer extension ( Figure 8C) can be detected and quantified.
A ligand bound on a folded target molecule can be detected similarly ( Figure 8D). As the structural unfolding of the target molecule can happen only after the ligand dissociates, ligand binding can stabilize the target, which often results in a slower unfolding transition. Thus, if the ligand results in a detectable delay in the unfolding transition after jumping to a higher A ligand bound on a mechanically exposed binding site at a higher force causes a delay in refolding after a force jump to a lower force. (D) A ligand bound on a folded structure at a lower force causes a delay in unfolding after a force jump to a higher force.

Biochemistry
Perspective force, it becomes a readout on whether the target molecule is bound by a ligand or not right before the force jump. While delayed unfolding has mainly been applied to quantify protein−DNA interactions, 77,91,92 the same principle can be utilized to detect ligand binding to a protein domain (e.g., RIAM binding to talin R3) and thereby quantify the forcedependent dissociation constant.
On the basis of these delayed structural transitions, the bound and unbound states of the target molecule can be determined, which enables the determination of the equilibrium binding probability p on (F) using a force jump assay. In such an assay ( Figure 8C,D), many cycles of force jump between two force levels are used: (i) a binding force (F) at which the target binding site is stable for a duration of

■ DISCUSSION
In this short review, we have sought to provide a brief discussion of some of the force-dependent considerations at the heart of mechanobiology and outline some of the strategies for measuring and studying them. The requirement to consider the force dependence of binding constants in biology necessitates development of the existing mathematical descriptions of binding constants to include mechanical descriptors. In addition, we discuss the novel experimental approaches required to measure them. The examples of force-dependent binding events illustrated in Figures 3−5 highlight some of the diverse and ingenious ways that biological systems sense, and respond to, forces. Even in these simplified in vitro systems, it is evident that the ways in which binding constants in each scenario are affected by force are complex and often biphasic. The consequence of these complex force dependencies is that the binding affinities between two ligands can change, either increasing or decreasing, by ≥1000 fold dynamically over a physiological force range. This creates incredible complexity in these mechanical linkages, with the same components assembling differently in different force environments.
Multivalent Interaction. Our review of the force-dependent affinity has been developed on the basis of single-site binary interactions. In many cases, such as antibody−antigen interactions, however, multivalent interactions play a crucial role in biological functions. In such cases, the binding strength of multivalent interactions cannot be formulated on the basis of a two-state model. The binding strength or the functional activity of such multivalent interactions is often termed "avidity". Such multivalent interactions are also implicated in mechanosensing reactions. For example, talin contains 11 VBS. 35 When multiple VBSs in talin are mechanically exposed for binding to vinculin, the functional activity of the forcedependent talin−vinculin interaction is expected to be further boosted.
Multiplexing Force-Dependent Factors. A hint about the immense amount of information that can be encoded by such mechanosensitive complexes is evident when you consider these principles in a simple two-component system and the multiplexing of the multiple force-dependent contributions that can emerge. Take, for example, the mechanosensitive interactions between talin and vinculin at the core of the focal adhesion. There is the well-documented force dependence of the exposure of cryptic VBS from within the core of the folded talin bundles discussed above. Above a certain threshold, talin bundles unfold exposing cryptic VBS and allowing vinculin to bind. The description in Figure 4 presents a single VBS binding a single vinculin. However, talin contains 11 VBS, and the mechanical response of talin is complex, 35 with diverse force thresholds governing exposure of each VBS. This setup is further impacted by additional forcedependent considerations.
For instance, both talin and vinculin are further regulated by autoinhibition, and in both cases, the affinity of the autoinhibition (which is a head−tail interaction) is reduced when the protein is under force; here, the effective binding constant of the autoinhibition is massively increased as the two interacting domains (the head and the tail) are physically held apart from each other ( Figure 1). As such, layers upon layers of inhibition exist on these molecules all exhibiting force dependence. 93 The forces exerted on talin change constantly, and if the force decreases, then the talin rod domain and the head tail autoinhibition will have a higher effective affinity (they will no longer be held apart from each other) and can be trying to revert back to their closed inactive conformations.
There is also considerable hysteresis on talin domain refolding. 25,35 A domain that unfolds at 15 pN will not refold when the force drops to <15 pN; instead, it requires forces of <3 pN to refold on a reasonable time scale. When refolding of the talin domains occurs, it will affect the exposure and thus the affinity of the talin−vinculin interactions. With up to 11 VBS in talin, many force linkages can be coupled.
However, additional force dependence is introduced by the connectivity of vinculin with the actin filaments. If vinculin engages talin via its head domain and couples to actin via its tail domain, then this will exert force on the talin−vinculin interaction, which will affect the affinity, and as the geometry of this force on the VBS−vinculin interaction is likely to be shearing force geometry, it means that force will strengthen this interaction (with catch bond kinetics). In addition, tethering vinculin to talin and actin also restricts vinculin autoinhibition, so the effective K d of vinculin autoinhibition is also greatly increased, which will further enhance the VBS− vinculin interaction.
This description is complex but still does not include the potential force dependence of unfolding of vinculin domains, the strength of the talin linkages with the integrin and actin, or the recruitment, or displacement, of factors as a result of domain unfolding that might enhance or decrease contractility. All of these will further augment the mechanical connections. Therefore, even within this simplified description of the system, there are almost endless possibilities for diverse outcomes. When you layer on the myriad of other mechanoeffectors and regulators that assemble on this hyperplastic framework, it becomes apparent that there is a huge capacity in these linkages to encode vast amounts of data. 25 Biochemistry Perspective Future Perspective. This review is centered around the molecular interactions involved in mechanosensing, with a focus on how force applied to a molecule may influence the affinity with its binding partners. In addition, the review also briefly discusses how force applied to interacting partners may affect the lifetime of the complex.
What has been missed in these discussions is how force applied to one molecule affects the binding and unbinding rates of its binding partners in solution. Because force can drastically change the affinity of the interaction, it should have a significant influence on the binding rate, the unbinding rate, or both. Compared to force-dependent affinity, the forcedependent kinetics of binary interactions have not been as extensively studied and remain less well understood. A deeper understanding of the force-dependent interaction kinetics is required because in cells many interactions do not reach equilibrium. The force-dependent reaction rates will provide crucial insights into how non-equilibrium molecular interactions under mechanical force in cells can be understood.
Another interesting topic, which is not included here, is how other types of mechanical constraints may affect interactions. Cytoskeleton filaments, such as actin filaments, are often subjected not only to tensile force but also to rotational constraint. 94 The latter will result in torque in the filament, which is transmitted to the proteins, such as formin or other actin capping proteins, linked to the end of the actin filaments. 95 Another example is that DNA in many topologically isolated chromatin domains is supercoiled, which also results in torque applied to the DNA. 96,97 How torque may affect the affinity and the kinetics of the molecular interactions is currently poorly understood, which should be another interesting future direction.
The rapid advances in the field of mechanobiology are making great strides in advancing our understanding of these complex mechanosensing signaling systems.

■ MODEL AND PARAMETERS OF FORCE−EXTENSION CURVES OF TARGET MOLECULES
Two-State Binary Interactions. Single-stranded DNA (ssDNA) is described by the worm-like chain (WLC) polymer model that contains two parameters, bending persistence length A ssDNA and contour length L = n nt l nt , where n nt is the number of nucleotides. ssDNA has a persistence length A ssDNA of ∼0.7 nm and a contour length per nucleotide l nt of ∼0.7 nm. 98  Double-stranded DNA (dsDNA) is described by the WLC model where the dsDNA has a persistence length A dsDNA of ∼50 nm and a contour length per base pair l bp of ∼0.34 nm. 100 Similarly, the force−extension curve of dsDNA x dsDNA (F) can be obtained by solving the inverse function of the Marko− Siggia formula.
Stiffening protein−dsDNA and bending protein−dsDNA are both described by the WLC model where the only difference compared with the naked dsDNA is the bending persistence length. In panels E and F of Figure 3, the bending persistence length of stiffening protein−dsDNA (A sf−dsDNA ) is taken to be ∼174 nm, which is adapted from the case of H-NS binding. 4 In panels H and I of Figure 3, the bending persistence length of bending protein−dsDNA (A bd−dsDNA ) is taken to be ∼30 nm, which is adapted from the case of IHF binding. 9 Binary Interactions Involving Autoinhibition. The dsDNA hairpin is described by the rigid body with a size L hp of ∼2 nm, which is the width of dsDNA. 101 Its force−extension curve is ( ) The α-helical bundle is described by the rigid body with a size L bd of ∼3.4 nm, which is the N-terminus−C-terminus distance obtained from talin R10 (Protein Data Bank entry 2KVP). 102 The α-helical chain is described by the freely jointed chain (FJC) polymer model where contour length L ahc = n ah L ah , where n ah is the number of α-helices and L ah is the size of each α-helix. In Figures 4D and 5B, fixed values of n ah and L ah are taken to be 4 and 4.8 nm, respectively. 102 Its force−extension curve is ( ) The α-helix form of VBS is described by the rigid body where size L ah,VBS = N VBS l ah,aa , where N VBS = 25 is the number of amino acids in VBS 103 and l ah,aa = 0.15 nm is the extension per amino acid in the form of an α-helix. 104,105 The unstructured peptide chain is described by the WLC model where the peptide chain has a persistence length A pt of ∼0.8 nm and a contour length per amino acid l aa of ∼0.38 nm. 106 Similarly, the force−extension curve of peptide chain x pt (F) can be obtained by solving the inverse function of the Marko−Siggia formula.