https://www.nature.com/articles/s41467-024-52749-w Skip to main content Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript. Advertisement Advertisement Nature Communications * View all journals * Search * Log in * Explore content * About the journal * Publish with us * Sign up for alerts * RSS feed 1. nature 2. nature communications 3. articles 4. article Engineering tunable catch bonds with DNA Download PDF Download PDF * Article * Open access * Published: 12 October 2024 Engineering tunable catch bonds with DNA * Micah Yang ORCID: orcid.org/0009-0006-9822-2892^1, * David t. R. Bakker ORCID: orcid.org/0009-0005-5284-6205^1 & * Isaac T. S. Li ORCID: orcid.org/0000-0002-8450-5326^1 Nature Communications volume 15, Article number: 8828 (2024) Cite this article * 3332 Accesses * 1 Citations * 74 Altmetric * Metrics details Subjects * DNA nanostructures * Molecular conformation * Nanostructures * Single-molecule biophysics Abstract Unlike most adhesive bonds, biological catch bonds strengthen with increased tension. This characteristic is essential to specific receptor-ligand interactions, underpinning biological adhesion dynamics, cell communication, and mechanosensing. While artificial catch bonds have been conceived, the tunability of their catch behaviour is limited. Here, we present the fish-hook, a rationally designed DNA catch bond that can be finely adjusted to a wide range of catch behaviours. We develop models to design these DNA structures and experimentally validate different catch behaviours by single-molecule force spectroscopy. The fish-hook architecture supports a vast sequence-dependent behaviour space, making it a valuable tool for reprogramming biological interactions and engineering force-strengthening materials. Similar content being viewed by others [41557_2024] De novo DNA-based catch bonds Article Open access 24 June 2024 [43246_2023] A catch bond mechanism with looped adhesive tethers for self-strengthening materials Article Open access 14 August 2023 [41467_2024] Engineering an artificial catch bond using mechanical anisotropy Article Open access 08 April 2024 Introduction Cells communicate mechanically by applying tension across cell-cell and cell-environment adhesions^1,2,3,4,5. The transduction of specific signals across these adhesive junctions relies on the force-dependent behaviour of the receptors and ligands. Many biochemical pathways activate as a response to external force, and the degree of this activation can vary depending on the magnitude, direction, loading period, and loading rate of mechanical strain. For example, T-cell receptor (TCR) specificity between agonists and antagonists is enhanced 4-fold when tension is applied^6, focal adhesion maturation is mediated by sustained high forces and fast loading rates^7, cell polarization during migration is determined by direction-dependent tension between vinculin and actin^8, and E. coli avoids being flushed out by the high shear stress of urinal and intestinal tracts while retaining mobility under low shear stress^9. In all these cases, measuring the dissociation constant (K[d]) or lifetime in a zero-force environment (t[0]) is not enough to see the full scale of biologically relevant behaviours. In a purely hypothetical space, there are three adhesion responses to increasing force--the slip bond, which weakens (force-dependent lifetime (t[F]) decreases); the catch bond, which strengthens (t[F] increases); and the ideal bond, which is insensitive to force (t[F] remains constant). Given that all covalent and hydrogen bonds are slip bonds individually^10, one might assume that adhesive interactions that involve a number of these slip bonds in different geometries--as is the case with receptor-ligand interactions - would also result in slip bonds, albeit with a higher t[0]. However, this is not always the case. In biology, proteins demonstrate a colourful array of multi-phasic t[F], including biphasic catch-slip^8,9,11,12, 13,14,15,16,17,18,19 and triphasic slip-ideal-slip^20 or slip-catch-slip^13,21 bonds. The counter-intuitive catch behaviour is vital to many cellular functions and biological behaviours. For instance, the catch bond behaviour in TCR interactions is essential in tuning their specificity, making them more sensitive to antigens presented on sessile, tension-producing cells than the unbound ones. Cells that function in environments with fluid shear, such as E. coli in the intestine and neutrophils in the bloodstream^22, rely exclusively on catch bonds to facilitate rolling adhesion motility across a large range of external forces. Cell-matrix adhesion and mechanosensing involve focal adhesion proteins like integrins and vinculin, both of which exhibit catch bond character: Integrins enhance their substrate binding through a force-dependent conformation change and vinculin binds to actin in a direction-dependent catch bond^3,8,13. Similarly, cadherins and downstream mechanotransduction proteins such as catenins require catch behaviour to regulate cell-cell communication and maintain tissue integrity^11,23. Overall, catch bonds play a pivotal role in cell adhesion, motility, mechanosensing, communication, and organization across a wide spectrum of biological contexts. Biological systems have evolved catch bond behaviours as adaptive solutions to many challenges, employing a variety of mechanisms^9,24, 25,26,27,28,29,30. Given their biological importance, designs for artificial catch bonds have been put forward in an effort to expand our molecular toolbox^10,31,32,33,34,35,36,37. To create a framework for thinking about artificial catch bond designs, two very similar base-principle models were conceptualized in parallel - one from a synthetic chemistry perspective^10, and one from a biomolecular perspective^31. While these architectures were purely theoretical, they cemented a series of guiding principles which subsequent work^32 ,33,35 relied upon: i.e. the interactions must have a weak state S [w], a strong state S[s], and a force-dependent interconversion between the two that biases the S[w] at low force and S[s] at high force^12,24,38. These designs pull directly from the two-state, two-pathway model, one of the simplest theories that fit experimental data (the other being the one-state, two-pathway model) among many ways to model biological catch behaviour^12,15,39,40,41. Small-molecule mechanophores exhibiting catch bond behaviour were created by using force-dependent molecular geometry to change the impact of work on the scission point^36,37. While applicable to creating self-protecting materials, breaking covalent bonds occurs at forces on the nano-Newton scale. In contrast, molecular forces in biology occur on the pico-Newton (pN) scale. Using similar principles, a protein-based artificial catch bond was created^34. A protein-ligand complex (Dockerin G:Cohesin E) with two native unbinding pathways was pulled from different non-native anchor points, of which one combination produced catch behaviour. While more biologically relevant, the design methodology would only work under specific circumstances: the initial protein-ligand complex must have two existing unbinding pathways, with no predictability of whether different pulling geometry will produce catch behaviour. Around the same time as this work, a DNA-based catch bond was also created using a cryptic binding site hidden in a force-unfoldable DNA hairpin called the barrier^35. While the authors showed that the barrier length can alter the catch behaviour, they stated that the tunability of the design is relatively limited due to high sensitivity to the barrier length. Given these advancements in catch bond design, it is clear that a tunable catch bond architecture that enables a variety of catch behaviours within a biologically relevant range is needed. In this work, we present the fish-hook architecture, a highly tunable DNA catch bond. We show the conceptual design, probe the catch behaviour with single-molecule force spectroscopy (SMFS), and demonstrate its wide range of tunability both experimentally and through analytical simulation. Results Conceptual design of a DNA catch bond To create a catch bond, our strategy is to design a DNA construct with two dissociation pathways of distinct mechanical strengths. By biasing towards the weak pathway at low force and the strong pathway at high force, we can effectively create the catch behaviour. To do so, we designed a force-switchable structure with two rupturable domains, nicknamed the hook and the jaw (Fig. 1a). When force is initially applied at the two points, both the hook and jaw are in the unzipping geometry and compete to unbind. If the hook opens first, the closed jaw locks it in the weak unzipping geometry. On the other hand, if the jaw opens first, the hook converts to the stronger shearing geometry. We have accomplished our goal if we can predictably control which of the two pathways is active by force. Our previous work demonstrated that the unzipping of two DNA placed in series can be preferentially controlled by force and sequence (Supplementary Fig. 1)^42. Thus, we designed the jaw and hook sequences such that the t[F] of jaw unzipping is higher than hook unzipping at low force and vice versa at high force, intersecting each other at a crossover force (F[c]) (Fig. 1b). The differential unzipping t[F] ensures that the hook is more likely to open at F < F [c]. In contrast, the jaw is more likely to open at F > F[c], effectively creating a force-dependent bias towards the weak and strong pathways around F[c] (Fig. 1c). Given that the ultimate hook dissociation in the weak and strong pathways are unzipping and shearing, respectively, the switching from the weak pathway to the strong pathway as force increases will effectively create a slip-catch-slip behaviour (Fig. 1d). Fig. 1: Schematics of DNA catch bond design and expected behaviour. figure 1 a Schematic of a two-state mechanism that allows the DNA construct to dissociate via two pathways directed by force above or below a crossover force (F[c]). In the weak pathway, the closed jaw locks the hook in the unzipping geometry. In the strong pathway, the open jaw switches the hook to the shearing geometry. b The lifetimes (t[F]) of the hook unzipping and jaw opening are designed to intersect, enabling force-dependent pathway selection. c Below F[c], the construct follows the weak pathway, as the jaw remains closed, and thus the hook unzips. Above F[c], the jaw opens before the hook, steering the construct into the strong pathway, thereby increasing the proportion of constructs with an open jaw (P[jaw-open]) as force escalates. d As the force increases, the transition from the hook unzipping to shearing creates a slip-catch-slip behaviour. e Slip bonds are characterized by a unimodal rupture force probability density function (PDF) at a specific loading rate. f In contrast, catch bonds exhibit a bimodal rupture force PDF (catch-slip-catch shown here). Here, the high rupture force population (red) corresponds to high-force slip behaviour and the low rupture force population (green) arises from the catch behaviour. For comparison, a slip-only PDF is represented by a grey dashed line. The catch behaviour shifts subpopulations from high-force to low-force ruptures (grey arrows). g Additionally, the fraction of high rupture force in the bimodal distribution increases with force ramp speed or loading rate. Full size image Relating our system to the two-state, two-pathway model, the weak, closed-jaw conformation is S[w], while the strong, open-jaw conformation is S[s]. The conversion rates between S[w] and S[s] are the unbinding and rebinding rates of the jaw hairpin, which are slow due to the large loop. This slow conversion rate is characteristic of the two-state model and is necessary to explain the catch-slip behaviour of FimH/mannose^24, vinculin/F-actin^8, and cadherins in the X-dimer configuration^11. At faster conversions, the two-state model collapses to the one-state model^12, which assumes a single bound state with two competing unbinding pathways^12,15,39,41. The survival probabilities predicted from the one-state model fit catch-slip data collected from P- and L-selectins/PSGL-1^12,15 and myosin/actin^43. This general framework for building an artificial catch bond aligns with the guiding principles described earlier^10,31 ,32,33. The force-dependent dissociation of interactions can be experimentally examined by SMFS. The rupture force distributions in a force ramp experiment provide the mechanical signatures of the bond characteristics. For a slip bond, the dissociation probability monotonically increases with force over time until the cumulative dissociation probability approaches 1, leading to a unimodal rupture force distribution (Fig. 1e). However, for a slip-catch-slip bond, the local minimum in t[F] creates a local maximum in dissociation probability, resulting in an additional peak in the rupture force distribution. This bimodal distribution is the first signature of a catch bond (Fig. 1f). Additionally, the force ramp speed controls the population ratio of the bimodal rupture force distribution, with higher ramp speed biasing towards more high-force ruptures (Fig. 1g). This force-dependent ratio of the two rupture force populations is a second signature of a catch bond. Characterization of the DNA catch bond by SMFS To implement the above strategy, we created a fish and a hook, two DNA pieces that act as receptors and ligands in this artificial catch bond. In our dual trap optical tweezers experiment, the hook was bound to one bead and the fish to the other (Fig. 2a). With this configuration, we began an approach-dwell-retract fishing pattern where the fish bead would approach the hook, dwell for a few seconds, and then retract. The fish contains a closed jaw and a single-stranded region complementary to the hook (Fig. 2a, b - regrettably, the jaw does not bite the hook). In order to show the two signatures of a catch bond (Fig. 1f, g), we designed the hook and jaw sequences by considering the ramp speed limits of our optical tweezers (max 2000 nm s^-1). We used an analytical model to find t[F] for sequences of various lengths and CG contents, and then used Monte Carlo simulations that emulate our pulling experiments to find a set of sequence parameters that would show a clear shift towards high-force rupture population at our available ramp speeds. Further details about the model and simulation can be found in the Methods section. The construct we selected has a jaw length of 18 base pairs (bp), of which 8 are C/G (jaw 8/18); and a hook length of 9 bp, of which all 9 are C/G (hook 9/9). (Fig. 2b) This construct is predicted via the analytical model to have the slip-catch-slip t[F] shown in Fig. 2d with an F[c] of 12.6 pN. To ensure the jaw can open fully, the loop must be at least as long as the two jaw strands combined; in our case, our jaw was 18 bp long, so the loop must be at least 36 nucleotides (nt) (Fig. 2b). We designed the poly-T section of the loop to be 42 nt long as an extra precaution. Fig. 2: Experimental setup and catch bond sequence design. figure 2 a Schematic representation of the dual-trap optical tweezers experiment. The hook is attached to bead A, and the fish is attached to bead B via DNA handles of 2633 and 272 bp, respectively. Both handles are attached to the beads by a biotin-streptavidin linkage on the 5' end of each handle (top inset). Bead B undergoes an oscillatory motion with a dwell period before retraction (bottom inset) to increase the fishing success rate (binding between the hook and the fish). b The construct sequence. A 9-mer polyethylene glycol (PEG) linker attaches the hook to its DNA handle to minimize nonspecific base interactions with the fish and increase flexibility. c Schematic showing the weak and strong pathways. d The expected t[F] of the selected construct sequence based on our analytical model, showing the overall expected slip-catch-slip behaviour resulting from transitioning from hook unzipping to hook shearing around the intersection of the jaw opening and the hook unzipping t[F] at F[c] = 13.6 pN. The anticipated force-extension curves differ for the weak (e) and strong (f) pathways, as predicted by the extensible worm-like chain model (XWLC). In the weak pathway, we expect a singular low-force rupture event indicative of hook unzipping (e). In the strong pathway, we expect a jaw-opening transition preceding a high-force rupture event corresponding to hook shearing (f). Full size image Once the hook successfully catches the fish, it will unbind via one of the two pathways (Fig. 2c). In the weak pathway (F < F[c]), the hook unzips first, leaving a low rupture force signature characteristic of DNA unzipping (Fig. 2e). Conversely, in the strong pathway (F > F[c]), the jaw opens before the hook unbinds, creating an unfolding event with a change in contour length (DL[c]) corresponding to the loop size. As the force continues to increase, the hook, now in the shearing geometry, will eventually rupture at a much higher force (Fig. 2f). The two pathways can be clearly distinguished by examining their force-extension signatures. If the system undergoes the weak pathway, the fish remains folded in the original conformation with its jaw closed. If the system undergoes the strong pathway, the jaw will remain open until the hook unbinds but will refold and revert to its initial jaw-closed conformation. Hence, the fish-hook catch interaction can be indefinitely repeated. In the fishing experiment, we pulled every tether until it broke to ensure both a single tether interaction and a new fish-hook pair formed each time. The force-extension curves showed two distinct behaviours: a population of curves that break at low forces (10-15 pN), corresponding to hook unzipping in the weak pathway (Fig. 3a, blue traces), and another population of curves that break at high forces (30-50 pN), corresponding to hook shearing in the strong pathway (Fig. 3a, red traces). The latter is always accompanied by an additional low force (10-15 pN) unfolding event prior to rupture with a contour length difference (DL[c]) of 20.4 +- 0.3 nm, which is in excellent agreement with the theoretical DL[c] of 20.2 nm (Fig. 3b) of the jaw-opening event in the strong pathway. Additionally, the force-extension curves of the weak pathway align with the pre-rupture portion of the strong pathway, indicating both are in the jaw-closed conformation. This evidence demonstrates that the fish-hook system exhibits the two designed pathways. We used an approach-dwell-retract fishing pattern to maximize the likelihood of catching a fish on any given approach (Fig. 2a, inset). The dwell time of ~8 s allowed the 9 bp hook to find a fish reliably while also allowing any fish with open jaws from the previous pull to close their jaws, preventing any history-dependent behaviour. We tested multiple beads in multiple samples for each experimental condition. Fig. 3: Experimental and simulation results from the optical tweezers. figure 3 a Experimental force-extension curves in chronological order from a medium ramp experiment, showing non-history dependent switching between strong (red) and weak (blue) rupture events for the hook9/ 9-jaw 8/18 construct shown in Fig. 2b. Traces end when the hook unbinds. Fishing attempts that do not yield single tethers are not shown. b The same force-extension curves in (A), but with theoretical worm-like chain curves overlaid for the two possible states (dashed lines). The distribution of DL[c] (inset) from all medium pulls with a jaw opening event, as compared with the theoretical DL[c] of 20.20 nm. c Rupture force distributions at three pulling rates compared to the Monte-Carlo simulation. The slow pulling rate has a higher proportion of hook unzipping, while the fastest pulling rate has a higher proportion of hook shearing. d Rupture force plotted against force loading rate for the three ramps, calculated at the moment of rupture. Due to the long DNA handles, each pulling speed (nm s^-1) has a range of loading rates (pN s^-1) obtained from the worm-like chain model. e The proportion of open-jawed constructs as a function of pulling rate compared to the simulation. Two additional catch bonds were tested in addition to the one shown in the rest of Figs. 2 and 3 (h[9/9]j[8/18]), one with an 11/11 CG content hook (h [11/11]j[8/18]; yellow; n = 101, 178, 130 for slow, medium, and fast ramp speed) and one with a 10/18 CG content jaw (h[9/9]j[10/18]; green; n = 253, 490, 187 for slow, medium, and fast ramp speed). Error bars are the 95% confidence interval of the proportion; n = 1000 for each simulated point. Full size image Next, we examined the rupture force distribution of the fish-hook catch bond at fast (2067 nm s^-1), medium (206.7 nm s^-1), and slow (51.7 nm s^-1) ramp speeds (Fig. 3c). In all three experimental sets, the fish-hook system exhibits the characteristic bimodal distribution of a catch bond. By checking for the presence of a jaw-opening event, we are able to determine the dissociation pathway and colour code them in the histograms, where the lower force peak is dominated by the weak hook unzipping pathway, while the high force peak is dominated by the strong hook shearing pathway. Our simulation of the system demonstrates quantitatively similar bimodal distribution as the experiment, taking into account the non-linear elasticity of the DNA handles in the optical tweezers experiment. To examine this further, we plotted rupture force as a function of the loading rate at the time of rupture, colour-coded by the identity of the pathway (Fig. 3d). As expected, the rupture forces from both pathways increase with ramp speed and loading rate. The simulated results under the same conditions overlayed well with the experimental results. The scatter distribution from each ramp speed is curved due to the non-linear elasticity of the DNA handle used in the optical tweezer experiment, where the loading rate (pN s^-1) scales with the spring constant (pN nm^-1) at constant ramp speeds (nm s^-1). Lastly, we examined the fraction of rupture events undergoing the strong and the weak pathways at different ramp speeds. The fraction of rupture events from the strong pathway increased from 0.33 at the low ramp speed to 0.74 at the fast ramp speed (Fig. 3e, pink). This, again is in good agreement with our simulation results indicating that higher ramp speed leads to a greater fraction of rupture events from the strong pathway (Fig. 3e), demonstrating our fish-hook system exhibits the designed catch bond behaviour. Our experimental results cover approximately 2 orders of magnitude in ramp speed, which are at the limit of our instrument (high ramp speed) and reasonable experimental time scale (low ramp speed). We additionally tested 2 other fish-hook catch bonds, one with a longer hook (Fig. 3e, yellow, CG content of 11 bp, length of 11 bp (11/11 hook)) and one with a more CG-rich jaw (Fig. 3e, green, CG content of 10 bp, length of 18 bp, (10/18 jaw)). The longer hook resulted in a higher proportion of jaw-opening at the same speeds. This makes sense because a longer hook is a stronger hook, making the jaw weaker by comparison and allowing it to open more often. The opposite is true of creating a more CG-rich jaw: if the jaw is stronger, the hook unzips more often, shifting the proportion of jaw opening events down at the same speeds. As with the original fish-hook, these trends are mirrored by the model. Design of fish-hook catch bonds in biological ranges The design space of the fish-hook catch bond is vast. Using the analytical model, we searched all hook/jaw sequence combinations with varying CG contents and lengths from 7-30 bp for catch behaviour. For a combination to have catch behaviour, the following conditions must be met: First, the log(t[F]) curves of the hook and jaw must intersect. The intersection creates the force-dependent state switching behaviour and dictates the force range (around F[c]) at which it happens. Second, the log(t[F]) curves of the hook and jaw must intersect in the correct order: the t[F] of the hook must be lower than the jaw at lower force and higher than the jaw at higher force. In other words, the slope of log(t[F]) for the jaw must be steeper than the hook. In the opposite case, the hook would shear at low force and unzip at high force, leading to a more exaggerated slip behaviour rather than slip-catch-slip. The slope of log(t[F]) in hairpin unzipping is almost entirely dependent on the sequence length, rather than CG content^42, thus the jaw should be longer than the hook. Third, there must exist a positively-sloped (catch) region in the overall log(t[F]). Even if the hook and jaw intersect correctly, the overall t[F] may not have a catch region for two reasons: Firstly, when the slopes of the hook and jaw are very similar, the state switching from hook unzipping to shearing occurs too gently for log(t[F]) to be positively sloped. Secondly, as F[c] approaches 0, the difference between the shearing and unzipping log(t [F]) of the hook is less prominent and disappears at F = 0, thus reducing both the height and width of the catch region. Within these bounds, we found 96645 combinations satisfying the first two conditions; of those, 54354 had a positively-sloped catch region (Supplementary Figs. 4, 5); of those, 3493 were within a biologically relevant range (we define this as t[F] in the catch region between 0.01-100 s and the force marking the start of catch region (F[start]) between 0-10 pN). We found t[F] located at F[start] and F[end] (Fig. 4d) for the catch region of each construct (Fig. 4a) and calculated the force range (DF) and t[F] range (Dlog(t[F])) that characterize the catch region of each construct (Fig. 4b). Both graphs (Fig. 4a, b ) are colour-coded by the 2D colormap shown in the inset of Fig. 4b. Fig. 4: Fish-hook catch bond design space. figure 4 a Catch regions of fish-hook catch bonds within the biologically relevant window of 0.01