Increasing evidence of mechanical force as a functional regulator in smooth muscle myosin light chain kinase.

Mechanosensitive proteins are key players in cytoskeletal remodeling, muscle contraction, cell migration and differentiation processes. Smooth muscle myosin light chain kinase (smMLCK) is a member of a diverse group of serine/threonine kinases that feature cytoskeletal association. Its catalytic activity is triggered by a conformational change upon Ca2+/calmodulin (Ca2+/CaM) binding. Due to its significant homology with the force-activated titin kinase, smMLCK is suspected to be also regulatable by mechanical stress. In this study, a CaM-independent activation mechanism for smMLCK by mechanical release of the inhibitory elements is investigated via high throughput AFM single-molecule force spectroscopy. The characteristic pattern of transitions between different smMLCK states and their variations in the presence of different substrates and ligands are presented. Interaction between kinase domain and regulatory light chain (RLC) substrate is identified in the absence of CaM, indicating restored substrate-binding capability due to mechanically induced removal of the auto-inhibitory regulatory region.

Introduction

All cells need to withstand as well as actively generate forces during division, differentiation or for their differentiated function. These mechanically governed adaptive processes require the translation of mechanical signals into physicochemical signals inside the cell that trigger an appropriate biological response (Shivashankar et al., 2015). A number of mechanosensors have evolved for these purposes, including mechanosensitive ion channels. The cytoskeleton, comprising microtubules, intermediate filaments and contractile actin-myosin filaments, plays a key role in maintaining cell shape against internal and external forces, but also emerges as a main hub for mechanosignaling. The regulation of actomyosin contraction differs significantly between the calcium-activated striated muscles and myosin-light chain phosphorylation in smooth and non-muscle cells. Four myosin light chain kinases (MLCK) exist, transcribed from the MYLK1 to 4 genes (reviewed in Chang et al., 2016).

Smooth muscle myosin light chain kinase (smMLCK) is a ubiquitously expressed serine/threonine kinase that particularly contributes to the regulation of smooth muscle contraction (Gallagher et al., 1997; Hong et al., 2011). SmMLCK phosphorylates the regulatory light chain (RLC) of smooth muscle myosin, which in turn triggers ATPase activity of the myosin heads. This activation results in the myosin power stroke – the fundamental process of muscle contraction (Dillon et al., 1981). Different MLCK homologues exist for skeletal and cardiac muscle tissue, but only in smooth muscle cells the contraction is directly linked to the activity of smMLCK (Zhang et al., 2010). Other regulatory proteins such as troponin are not present in smooth muscle cells. In contrast to skeletal or cardiac MLCK, smMLCK is not restricted to muscle tissue but is expressed in almost all mammalian cells. Further studies have revealed the importance of smMLCK in additional cellular pathways besides muscle contraction such as platelet aggregation, exocytosis, and cell migration (Hashimoto et al., 1994; Kumakura et al., 1994; Chen et al., 2014).

The most important regulator of smMLCK activity is Ca2+-loaded calmodulin (Ca2+/CaM). Without external activation, smMLCK’s catalytic turnover is suppressed by a pseudosubstrate mechanism, i.e. its catalytic core is auto-inhibited by a regulatory element that mimics the RLC substrate sequence (Pearson et al., 1988). Binding of Ca2+/CaM, however, initiates a conformational change in this regulatory region and thus removes auto-inhibition. Release of the regulatory sequence facilitates RLC binding and catalytic activity of smMLCK (Bagchi et al., 1992).

Several isoforms are encoded by the smMLCK gene (MYLK1) by alternative initiation sites and differential splicing (Lazar and Garcia, 1999), containing an N-terminal extension in addition to the core catalytic kinase domain. The core sequence is highly conserved in its domains and in their respective order. It comprises immunoglobulin-like domains, a fibronectin-like domain, a proline-rich presumably elastic region as well as the catalytic kinase domain. By homology with other members of the MLCK family (Chang et al., 2016), the kinase domain is composed of a smaller N-terminal lobe that contains the ATP binding site and a larger C-terminal lobe that is responsible for substrate recognition (Gallagher et al., 1997). The latter is auto-inhibited with respect to RLC binding by the regulatory element with adjacent CaM binding region. The active site is located at the interface between the two lobes. So far, the molecular structure of smMLCK and how it conveys this regulation have not been fully elucidated. However, a high degree of structural similarity to titin kinase (TK) and twitchin kinase – two prominent serine/threonine kinases from the giant muscle protein titin and the titin-like protein twitchin found in invertebrate muscles – as well as to the recently solved MYLK4 can be assumed (Gautel, 2011; Chang et al., 2016). The terminal domains of smMLCK are also highly conserved and form binding regions to F-actin (N-terminal) and myosin (C-terminal) (Sellers and Pato, 1984; Hong et al., 2009; Gautel, 2011); longer isoforms (MLCK-210) also seem to interact with other cytoskeletal components through their N-terminal domains (Kudryashov et al., 2004). Specific binding to myosin presumably increases the affinity between the kinase domain and RLC due to local proximity (Silver et al., 1997). The actin-binding domain allows smMLCK to associate along actin filaments and thus enhances its phosphorylation rate in smooth muscle (Hong et al., 2015). The fact that these binding sites are located at the termini of the molecule suggests that smMLCK might connect simultaneously with both myosin and actin, and is theoretically capable of bridging thick and thin filaments in smooth muscle due to extensible linker regions in the proline rich repeat segment (Mabuchi et al., 2010). This cytoskeletal association of smMLCK could hence significantly contribute to stiffness and passive tension of smooth muscle, or to responses in external stress. An intriguing interplay exists between smMLCK and mechanical forces in some tissues; for example, repeated contractile activation leads to increased contractility in airway smooth muscle (Fairbank et al., 2008). Such mechanosensitive conformational modulation might be comparable to the role of TK in striated muscle cells. Titin bridges the thick and thin filament systems in the sarcomere in an analogous manner to smMLCK. It may thus act as a muscle mechanosensor, signaling exposure to mechanical tension in the contracting and mechanically stressed sarcomere. TK is also intrasterically regulated by a pseudosubstrate mechanism, but its auto-inhibition is understood to be released upon mechanical stress rather than allosterically by an effector molecule. Single-molecule studies as well as molecular dynamics simulations established that partial unfolding of TK upon external force results in a controlled release of the regulatory segment (Puchner et al., 2008; Gräter et al., 2005). This process forms an enzymatically active intermediate capable of ATP binding and subsequent substrate turnover. Ca2+/CaM affinity has also been detected for TK, but its binding shows only a stimulating effect rather than a full activation of its turnover (Mayans et al., 1998). Due to these striking similarities between smMLCK and TK both in structure and in their actin-myosin association in the muscle, the existence of a comparable Ca2+/CaM-independent regulation mechanism for smMLCK is plausible (Chang et al., 2016).

While the established Ca2+/CaM activation mechanism is the most prominent and best-understood for the activation of smMLCK, other activation/regulatory factors are likely to exist, some of which are already identified (Stull et al., 1993; Pfitzer, 2001). In this study, the mechanical response of smMLCK was probed via single-molecule force spectroscopy with an atomic force microscope (AFM) to understand how this cytoskeletal kinase is conformationally modulated by external forces. We investigated the effects of the presence of different ligands such as ATP, Ca2+/CaM or RLC peptide substrate or a combination of these on smMLCK’s pathway through its different conformational states until fully unfolded. Ligand interactions and the resulting changes to the smMLCK energy potential give further insights into its intrasterical regulatory mechanism and whether substrate binding can be enabled by the application of mechanical force.

Results

Force response of smMLCK during AFM-based force spectroscopy

The investigated molecular construct is a truncated version of the smooth muscle isoform encoded in the smMLCK gene (Lazar and Garcia, 1999) lacking the proline rich region and the N-terminal actin-binding domain. The remaining sequence comprises the kinase domain, its neighboring fibronectin-like domain (Fn3) and two N-terminal Ig-like domains (Ig, Ig) as well as the C-terminal Ig-like myosin-binding domain called telokin (Ig) (Figure 1A). For force spectroscopy experiments, the smMLCK is specifically tethered by an N-terminal Strep-tag II via an AFM cantilever tip that is functionalized with a monovalent variant of Strep-Tactin (Baumann et al., 2016). Force is applied in physiological pulling geometry with the PEG spacers minimizing unspecific protein-surface interactions.

Figure 1.

Overview of the experimental configuration for applying controlled mechanical stress to smMLCK.

(A) Schematic illustration of the investigated smMLCK construct. It consists of the kinase domain surrounded by several Ig-like domains (Ig) and a fibronectin-like domain (Fn3). Possible substrate interactions are indicated (ATP, Ca2+/CaM and RLC). RLC interaction is prevented by the auto-inhibitory pseudosubstrate sequence that is released upon Ca2+/CaM binding. For covalent attachment onto the surface, the construct harbors a C-terminal ybbR-tag. (B) Representative force-distance curves (red, blue) depicting the characteristic transitions of the kinase through different conformational states (S, S, S) and subsequent unfolding of the adjacent Fn3 and Ig-like domains. Whereas most force-distance curves rupture before or after Fn3 unfolding (as shown in red) due to comparable rupture forces of Fn3 and the employed handle system, the blue curve illustrates a descriptive example with additional unfolding of Fn3 and Ig-like domains depicting the further force-distance pattern given by the construct. Structural interpretation and assignment of the detected force-distance pattern is schematically depicted above the curve. (C) Contour length transformation of 99 unfolding events with respective contour length increments. L and L are released at the transition of the kinase domain from conformational state S to S and S to S respectively. The contour lengths of Fn3 and Ig-like domains are additionally depicted.

DOI: http://dx.doi.org/10.7554/eLife.26473.002

DOI: http://dx.doi.org/10.7554/eLife.26473.003

Contour length plot of 99 unfolding events of MLCK with 0 mM ATP present, aligned as described in the data analysis section.

When the molecular construct is stretched with a constant speed, the investigated protein passes through a characteristic sequence of conformational states. These protein states represent semistable folding intermediates on the guided way through the protein's energy potential starting from a fully folded native structure. The conformations might also correspond to partially unfolded structures which are stable enough to form functional states. The transitions between these states are marked by distinct drops in the force-distance curves. In all cases reported here these transitions are accompanied by a lengthening of the molecular construct, which is manifested in a characteristic release of hidden contour length. The recorded force-distance response of smMLCK reveals a strict hierarchy in mechanical stability of the individual domains (Figure 1B). At low forces, the kinase domain unfolds by passing an intermediate conformational state S indicated by two force peaks at the end of state S at 30 nm and of S at 60 nm. Both transitions appear at forces of approximately 30 pN, measured at a retraction speed of 800 nm/s. At state S the kinase domain is fully unfolded. This characteristic sequence is followed by Fn3 unfolding, which is assigned to a single rupture event at around 100 pN in accordance with the respective contour length and previous data on domain strength (Rief et al., 1998; Li et al., 2005). Finally, the three Ig-like domains in the construct denature at forces of around 200–250 pN (Rief et al., 1998). Significant elongation of the surrounding linker length due to preceding kinase unfolding does not alter rupture forces of these latter domains in our study and is therefore unlikely to influence the hierarchy of rupture forces in the probed structure or its characteristic sequence of unfolding. Since the employed tethering complex (Strep-tag II:Strep-Tactin) substantially overlaps in its rupture force regime with the unfolding force of Fn3, most force distance traces end either before or after the Fn3 domain because the tethering complex ruptured before the Ig-like domains could unfold. A transformation of 99 force-distance curves into a contour length histogram (Figure 1C, Figure 1—source data 1) indicates that the rarely unfolded Ig-like domains and Fn3 contribute a contour length of 30 nm each. These increments are consistent with results obtained in other force spectroscopy measurements (Rief et al., 1998). The kinase domain is assigned with respect to its contour length to a rather clear peak at 61.79.6 nm (L released at transition S→S) and a less distinct broad ridge in the range of 25–40 nm (L released at transition S→S). The relatively wide contour length distribution of the broad ridge may be caused by the low force of this event and the consequently poor transformation of this segment: further aspects which might contribute to this poorly defined transformation are discussed in the context of the ATP-binding measurements below. For an estimation of the released contour length at S→S transition, individual force-distance curves with high rupture forces at this point were fitted with the Worm-Like Chain (WLC) model yielding an approximate length increment of around 30 nm. The measured contour length of approximately 92 nm for the full kinase domain is in good agreement with a simple approximation assuming 0.365 nm per amino acid (255 aa 0.365 nm = 93.1 nm) (Dietz and Rief, 2004). It indicates a fully unfolded kinase domain at S. The two-step unfolding behavior is consistent with the bi-lobed structure of the kinase’s catalytic core and suggests assignment of S→S transition to the unfolding of the smaller kinase lobe and of S→S to the larger lobe.

Conformational changes of smMLCK upon ATP binding

In contrast to TK, smMLCK’s capability of binding ATP is not inhibited through its pseudosubstrate mechanism. Independent from Ca2+/CaM activation, ATP interacts with a K of around 10 µM (Kennelly et al., 1992). In this study, ATP binding to smMLCK and the corresponding effects on its structure were identified by changes in the characteristic transition pattern through the kinase’s different conformational states during AFM-based force spectroscopy. ATP was added in buffer solution with a final concentration of 3 mM. A heatmap in force-distance space of 560 aligned and overlaid unfolding curves emphasizes ATP-induced changes: the smaller kinase domain lobe (S→S transition) is significantly stabilized upon ATP binding (Figure 2A). Since the order of released contour lengths L and L – associated with the transitions S→S and S→S – is not altered by this enhanced stability, the S→S transition seems to remain structurally shielded from S→S. The increased force signal of S→S allows precise extraction of the small lobe increment L in contour-length space. The determined 30.89.8 nm add to a total length of 92.67.5 nm for the full kinase domain with L being unchanged with a contour length of 61.88.9 nm (Figure 2B). Other contour length increments in the overall construct including Fn3 and Ig-like domains were unaffected by the interaction of kinase domain and ATP. Unfolding traces in absence and presence of ATP were realized within one experiment that is, same cantilever and sample surface. The measured forces are therefore directly comparable without uncertainties in the force signal that could originate from deviations in AFM spring constant calibration. The histograms for the unfolding force of the S→S transition with and without ATP indicate an increase of the most probable rupture force of about 30 pN (Figure 2C). Due to saturated binding conditions using 3 mM ATP, the obtained force histogram most likely represents exclusively the ATP-bound conformation and not a mixed population of ATP-bound and ATP-free states of the tethered smMLCK molecules. The peak forces for S→S and Fn3 were found to be unaffected by addition of ATP (Figure 2—figure supplement 1). The fact that Fn3 remains unchanged in its unfolding properties independent of ligand addition allows for relative comparison of measurements from different cantilevers. In the following, forces from different experiments are normalized according to the most probable peak force of Fn3 unfolding. Binding of ADP as well as the non-hydrolysable ATP analogue adenylyl-imidodiphosphate (AMP-PNP) could analogously be detected for smMLCK by a stabilization of the S→S transition in the tethered construct (Figure 2—figure supplement 2). Stabilization, however, appears to be weaker in comparison to the interaction with ATP.

Figure 2.

Structural effects of ATP binding on smMLCK’s characteristic sequence of conformational states.

(A) Stabilization of the S→S transition upon ligand binding. For better illustration, a heatmap of 560 aligned curves is depicted. (B) Contour-length transformation of the respective events in the presence of ATP. L is associated to the contour length released at the transition from state S to S. (C) Statistical evaluation of S→S stabilization via force histograms fitted with the Bell-Evans model. An increase in the most probable transition force of about 30 pN is observed upon ATP addition. Both data sets were recorded within one experiment with the same cantilever.

DOI: http://dx.doi.org/10.7554/eLife.26473.004

DOI: http://dx.doi.org/10.7554/eLife.26473.005

DOI: http://dx.doi.org/10.7554/eLife.26473.006

DOI: http://dx.doi.org/10.7554/eLife.26473.007

All data was collected within one experiment and absolute force values can directly be compared. Whereas S→S is stabilized by ATP as described in the main part of the manuscript, S→S and Fn3 appear not to be significantly changed by substrate interaction. The force histograms of the Fn3 domain are used in other experiments for normalizing forces to the same value.

DOI: http://dx.doi.org/10.7554/eLife.26473.008

Same as Figure 2—source data 2.

DOI: http://dx.doi.org/10.7554/eLife.26473.009

Same as Figure 2—source data 2.

DOI: http://dx.doi.org/10.7554/eLife.26473.010

DOI: http://dx.doi.org/10.7554/eLife.26473.011

DOI: http://dx.doi.org/10.7554/eLife.26473.012

DOI: http://dx.doi.org/10.7554/eLife.26473.013

DOI: http://dx.doi.org/10.7554/eLife.26473.014

DOI: http://dx.doi.org/10.7554/eLife.26473.015

DOI: http://dx.doi.org/10.7554/eLife.26473.016

DOI: http://dx.doi.org/10.7554/eLife.26473.017

DOI: http://dx.doi.org/10.7554/eLife.26473.018

DOI: http://dx.doi.org/10.7554/eLife.26473.019

Actual histogram of Fn3 forces is not shown here and is just used for normalization.

DOI: http://dx.doi.org/10.7554/eLife.26473.020

DOI: http://dx.doi.org/10.7554/eLife.26473.021

Actual histogram of Fn3 forces is not shown here and is just used for normalization.

DOI: http://dx.doi.org/10.7554/eLife.26473.022

Figure 2—figure supplement 1.

Effects of ATP or Ca2+/CaM addition on the peak forces for the respective transitions S→S and S→S and for the Fn3 unfolding force.

All data was collected within one experiment and absolute force values can directly be compared. Whereas S→S is stabilized by ATP as described in the main part of the manuscript, S→S and Fn3 appear not to be significantly changed by substrate interaction. The force histograms of the Fn3 domain are used in other experiments for normalizing forces to the same value.

DOI: http://dx.doi.org/10.7554/eLife.26473.008

Same as Figure 2—source data 2.

DOI: http://dx.doi.org/10.7554/eLife.26473.009

Same as Figure 2—source data 2.

DOI: http://dx.doi.org/10.7554/eLife.26473.010

DOI: http://dx.doi.org/10.7554/eLife.26473.011

DOI: http://dx.doi.org/10.7554/eLife.26473.012

DOI: http://dx.doi.org/10.7554/eLife.26473.013

DOI: http://dx.doi.org/10.7554/eLife.26473.014

DOI: http://dx.doi.org/10.7554/eLife.26473.015

DOI: http://dx.doi.org/10.7554/eLife.26473.016

DOI: http://dx.doi.org/10.7554/eLife.26473.017

Figure 2—figure supplement 2.

Stabilization of the S→S transition upon ADP or AMP-PNP binding.

DOI: http://dx.doi.org/10.7554/eLife.26473.018

DOI: http://dx.doi.org/10.7554/eLife.26473.019

Actual histogram of Fn3 forces is not shown here and is just used for normalization.

DOI: http://dx.doi.org/10.7554/eLife.26473.020

DOI: http://dx.doi.org/10.7554/eLife.26473.021

Actual histogram of Fn3 forces is not shown here and is just used for normalization.

DOI: http://dx.doi.org/10.7554/eLife.26473.022

Contour length plot of 560 unfolding events of MLCK in the presence of 3 mM ATP, aligned as described in the data analysis section.

Force histogram of S→S transition in the presence of 0 mM ATP.

Force histogram of S→S transition in the presence of 3 mM ATP.

Force histogram of S→S transition in the presence of 3 mM ATP, 25 µM CaM, 2 mM Ca2+.

Force histogram of Fn3 unfolding in the presence of 0 mM ATP.

Force histogram of Fn3 unfolding in the presence of 3 mM ATP.

Force histogram of Fn3 unfolding in the presence of 3 mM ATP, 25 µM CaM, 2 mM Ca2+.

Force histogram of S→S transition in the presence of 3 mM AMP-PNP, 30 µM CaM, 3 mM Ca2+, 280 µM RLC.

Force histogram of Fn3 unfolding in the presence of 3 mM AMP-PNP, 30 µM CaM, 3 mM Ca2+, 280 µM RLC is used for normalizing forces to the same value.

Force histogram data of S→S transition in the presence of 4 mM ADP.

Force histogram of the Fn3 unfolding in the presence of 4 mM ADP is used for normalizing forces to the same value.

Energy barrier attenuation upon Ca2+/CaM interaction

Importantly, an additional peculiarity of the kinase unfolding becomes evident from the ATP-stabilized case: an unusual stretching behavior of the S state at AFM distances of about 30–40 nm (Figure 2A). The force response in this region clearly deviates from a typical WLC behavior exhibiting an unusual kink. Apparently, a small but distinct energy barrier - a mechanical barrier stabilizing a specific protein conformation - has to be overcome at low forces (15 pN) before the kinase domain acquires conformational state S. The energy barrier is attributed to an additional transition S→S preceding the formerly described transitions S→S and S→S (Figure 3A). S→S represents a small, uncoupled conformational change with a barely resolvable contour length release. This additional energy barrier might also explain the ambiguous contour length transformations for the S→S transition in the measurements in the absence of ATP (Figure 1C). Due to the small forces of transitions S→S and S→S they are hard to distinguish in the unfolding pattern and their transformations result in a broadened overlap of their respective distributions.

Figure 3.

Structural effects of Ca2+/CaM binding on smMLCK’s characteristic sequence of conformational states.

(A) Attenuated S→S transition in the characteristic force-distance pattern of the smMLCK construct due to conformational changes upon Ca2+/CaM binding. This effect is emphasized by a heatmap comparison of several hundred overlaid force-distance curves. Both data sets were collected within one measurement. (B) Structural model interpretation. The S→S transition is assigned to a force-induced rearrangement in the kinase domain that correlates with the conformational changes induced by Ca2+/CaM binding – the release of the inhibitory pseudosubstrate.

DOI: http://dx.doi.org/10.7554/eLife.26473.023

The atypical stretching behavior indicating a transition from state S to S is still observable in the presence of Ca2+: only in combination with CaM is the barrier not detected in the unfolding pattern.

DOI: http://dx.doi.org/10.7554/eLife.26473.024

The existence of this energy barrier becomes particularly relevant as the presence of Ca2+/CaM affects this initial feature. Namely, addition of Ca2+/CaM to the measurement buffer leads to an absence of this barrier otherwise observed in the kinase force response (Figure 3A). Ca2+/CaM is understood to be a trigger for structural rearrangements in the kinase domain that enable RLC binding. The detected energy barrier could represent this rearrangement mechanism. We thus propose the following hypothesis: the release of the pseudosubstrate is already realized upon initial Ca2+/CaM binding to smMLCK and is therefore not observed. If not activated by Ca2+/CaM, however, this regulatory modulation of the smMLCK structure appears as distinct part of the transition pathway indicated by the S→S transition (Figure 3B). Addition of Ca2+ without CaM has no effect on the observed feature (Figure 3—figure supplement 1). Since the energy barrier precedes the complete process of guiding smMLCK through different conformations, the regulatory fragment appears to be released before the intact kinase domain gets denatured at all and therefore becomes inactive. This mechanism could therefore imply the existence of intermediate smMLCK conformations (S, S) with an intact active site and RLC binding capability independent of Ca2+/CaM induced by mechanical forces comparable to the regulation of TK.

Figure 3—figure supplement 1.

Missing effects by addition of Ca2+ without CaM.

The atypical stretching behavior indicating a transition from state S to S is still observable in the presence of Ca2+: only in combination with CaM is the barrier not detected in the unfolding pattern.

DOI: http://dx.doi.org/10.7554/eLife.26473.024

RLC peptide interaction

Based on these considerations, we designed experiments to test if conformational changes in the unfolding pattern can also be detected in the presence of RLC substrate. A minimal MLCK substrate peptide was used as RLC substitute that interacts only with the catalytic core of the kinase domain (Knighton et al., 1991). It consists of a truncated and slightly modified smooth muscle myosin light chain (MLC 11–23, P14A, Q15A) from chicken gizzard with the amino acid sequence KKRAARATSDVFA. It is effectively phosphorylated by smMLCK with a turnover rate of K = 7.5 µM (for chicken gizzard smMLCK) (Kemp and Pearson, 1985). Assuming an effectively blocked RLC binding site without activation by Ca2+/CaM, RLC addition should in principle not influence the structure of the smMLCK kinase domain due to the lack of interaction – unless there is binding by an active intermediate during the process of unfolding.

The characteristic patterns of transitions through different smMLCK states in the presence of RLC peptide are plotted via heatmaps of smMLCK force-distance curves in presence (n = 1013) and absence (n = 864) of ATP in the buffer solution in Figure 4A. The heatmaps clearly reveal a significant stabilization of the S state independent of ATP. This structural change clearly indicates interaction with the substrate and is thus interpreted as RLC binding. Quantitative evaluation of S→S transition force histograms was performed via kernel-density functions, since the recorded forces show deviations from a single-bond Bell-Evans distribution (Figure 4B). This might partly be based on the overlap of different populations e.g. RLC-bound and RLC-free kinase domain, but deviations from the typical Bell-Evans model are also observed for S→S in measurements without RLC substrate. Other factors such as interactions between the catalytic domain and adjacent Fn3-like domain or the pseudosubstrate sequence could be examples that already promote different transition pathways and thus yield this atypical behavior. Independent of shape, however, the recorded histograms reveal a clear shift in most probable rupture force in S→S transition by about 30 pN.

Figure 4.

Structural effects of RLC peptide binding on smMLCK’s characteristic sequence of conformational states.

(A) Qualitative observation of an increased mechanical stability in the large kinase lobe illustrated by the higher forces in the S→S transition. The effect is emphasized by heatmaps of aligned force-distance curves obtained under different substrate conditions. The stabilizing effect is detected independently of the presence of ATP. (B) Quantitative evaluation of the increased S→S transition force. The force histograms were approximated with a kernel-density function for extracting the most probable rupture force. It reveals a significant shift of about 30 pN due to the stable interaction of the RLC peptide with the catalytic core. Since this binding is stated to be prevented by an auto-inhibition process according to the conventional view of smMLCK activation, the experimental observation hints at an additional path of kinase regulation modulated by force.

DOI: http://dx.doi.org/10.7554/eLife.26473.025

DOI: http://dx.doi.org/10.7554/eLife.26473.026

DOI: http://dx.doi.org/10.7554/eLife.26473.027

DOI: http://dx.doi.org/10.7554/eLife.26473.028

Western blot probing RLC phospho-serine 19 following smMLCK in vitro kinase assay time-course shows only phosphorylation in presence of Ca2+/CaM.

DOI: http://dx.doi.org/10.7554/eLife.26473.029

Force histogram of S→S transition in the presence of 3 mM ATP and 280 µM RLC.

Force histogram of S→S transition in the presence of 280 µM RLC.

Discussion

Binding of RLC peptide is directly observed in our measurements in the absence of Ca2+/CaM as it clearly alters the characteristic transition pattern of single smMLCK molecules by stabilizing the S→S transition. This observation is in conflict with the blocking of the RLC binding site without Ca2+/CaM activation as proposed by the established pseudosubstrate inhibition mechanism (Bagchi et al., 1992). The experimental results could therefore represent a first indicator for a force-driven activation of the catalytic pathway of smMLCK (Figure 5). In the living organism, we assume that smMLCK bridges actin and myosin filaments in smooth muscle tissue and thus allows for a stretched conformation that additionally promotes kinase activity. By relative movements in the cytoskeleton, sufficient mechanical stress may be created within the protein regulating its catalytic activity besides Ca2+/CaM binding, as proposed in this study. Even for mechanical forces below 15 pN - the force that is identified to be necessary for activation in our experiments - steady physiological stress could likely remove autoinhibition if (semi-)permanently applied as the kinase is spanned between two filaments instead of being probed in a constant speed single-molecule force spectroscopy experiment. Despite a K in the µM range between smMLCK and myosin or actin (Hong et al., 2015), long-term association to the cytoskeleton is assumed for smMLCK especially due to several additional binding domains to both filaments in its long isoform (Kudryashov et al., 2004). Necessary forces and lifetimes of the spanned conformation could therefore in principle be reached in a physiological context to allow for the proposed force-driven activation. The in vivo existence and relevance of this regulation pathway, however, has to be examined in further studies. Our conclusions are drawn from single-molecule force spectroscopy measurements where we specifically tethered smMLCK and forced it through several conformations. The presented AFM-based approach represents a sensitive means of detecting ligand binding on the single-molecule level but can provide only limited temporal information about the interaction. From the barrier pattern, it is not directly discernible if the substrate binds in the course of the pulling (as assumed for the RLC) or if it is bound right from the beginning of measurements (as for ATP binding). In order to exclude pre-binding of the substrate, basal binding activity of RLC in absence of Ca2+/CaM activation was tested for the used construct via isothermal titration calorimetry. Due to the presumably very low binding affinities, interaction kinetics could not be detected at moderate concentrations to distinguish between basal and mechanically induced binding. The investigated construct does not show Ca2+/CaM-independent enzymatic activity, suggesting that RLC does not bind to smMLCK in the absence of Ca2+/CaM under zero tension (Figure 4—figure supplement 1). While beyond the scope of this study, hybrid approaches will be required to obtain a complete and conclusive picture of the binding mechanism and its physiological relevance, ideally complemented by structural information of the auto-inhibited state. In particular, combined force and fluorescence spectroscopy and molecular dynamics simulations will aid in this. Ultimately, we need to develop an assay that directly visualizes substrate turnover upon force-activation of these enzymes, as in-situ unbinding forces from actin and myosin filaments (which are yet unknown) will interplay with the mechanically induced intramolecular conformational changes. Direct measurements of force-induced substrate binding and activity, that is, RLC phosphorylation, will then lead to full comprehension of this alternative activation path for smMLCK. To this end, the powerful combination of single-molecule force spectroscopy and fluorescence spectroscopy in nanoapertures can provide the basis for in vitro force-activation assays (Heucke et al., 2013), to complement the presented findings that RLC substrate binds smMLCK under force in the absence of Ca2+/CaM.

Figure 5.

Structural interpretation of the stabilized S state upon RLC interaction.

Mechnical stress forces the construct into a conformational state S equivalent to the state reached by Ca2+/CaM binding. By release of the pseudosubstrate sequence the conformational state is capable of RLC binding which is detected by a significant stabilization of the S state. Both initially different activation pathways eventually result in the same sequence of conformational states, with the only difference being the presence or absence of bound Ca2+/CaM, depicted in light grey in the S state.

DOI: http://dx.doi.org/10.7554/eLife.26473.030

Figure 4—figure supplement 1.

Ca2+/CaM-dependent RLC phosphorylation of the investigated smMLCK construct.

Western blot probing RLC phospho-serine 19 following smMLCK in vitro kinase assay time-course shows only phosphorylation in presence of Ca2+/CaM.

DOI: http://dx.doi.org/10.7554/eLife.26473.029

Materials and methods

Expressed construct

The smMLCK construct used for this study is an 858 amino acid protein (808 aa, from 1097 to 1904 in human smooth muscle isoform 1 myosin light chain kinase, accession number NP_444253.3) incorporating Strep-Tag II (WSHPQFEK) at the N-terminus and a ybbR-tag (DSLEFIASKLA) and hexa-histidine tag at the C-terminus. The cDNA encoding smMLCK was cloned into the EcoRI and XhoI sites of a modified pENTR11 (Invitrogen) baculovirus shuttle plasmid region surrounded by attL1 and attL2 recombination sites (5’- attL1 - TCG AAG GAG ATA GAA CCA ATT CTC TAA GGA AAT ACT TAA CCA TGG CTA GCT GGA GCC ACC CGC AGT TCG AAA AAG GCG CCG AGA CCG CGG TCC CGA ATT CG - smMLCK - CCC TCG AGC GGT TCC GGT GGT GAC TCC CTG GAG TTC ATC GCT TCC AAG CTG GCT TCA GGC CTG AGA GGA TCG CAT CAC CAT CAC CAT CAC TAA GAT CCG TCG AGA TAT CTA G - attL2 - 3’). After generation of recombinant virus using a BaculoDirectBaculovirus Expression System (Invitrogen), production of smMLCK was carried out in suspension cultures of Spodoptera frugiperda Sf9 (Sf9 cells in Sf-900II SFM, Invitrogen). Cells were routinely maintained at 28°C, 100 rpm in the concentration range from 1  106 cells/ml to 8  106 cells/ml. 1 l culture of Sf9 cells at a concentration of 2.5  106 cells/ml was infected with smMLCK recombinant baculovirus of third generation (P3) and left growing for 3 days. Infected cells were then pelleted at 1000  g for 10 min. They then were resuspended in an ice-cold buffer containing 20 mM Tris (pH 8.0), 100 mM NaCl, 40 mM Imidazole, 14 mM 2-mercaptoethanol (buffer A) supplemented with cOmpleteEDTA-free Protease Inhibitor Cocktail (Roche) as per manufacturer’s recommendations. Cells were lysed by passing the mixture three times through a 0.8  40 mm syringe needle and treated with DNase I at a final concentration of 25 µg/ml for 10 min at 4°C. After centrifugation at 4000  g, supernatant was loaded onto a 2 ml His-Trap crude column (GE Healthcare Life Sciences) pre-equilibrated with buffer A. The column was then washed with 20 column volumes of buffer A and smMLCK protein was eluted with a step of buffer A containing 250 mM Imidazole. Peak fraction was purified further on HiLoad 26/600 Superdex 200 (GE Healthcare Life Sciences) equilibrated with 20 mM Tris (pH 8.0), 100 mM NaCl, 1 mM DTT. Peak fractions were pooled together and analysed on SDS PAGE for purity.

Enzymatic activity assay

To assess CaM-dependent and independent enzymatic activity of smMLCK, the construct used in AFM experiments was incubated with human RLC (NP_291024.1) in the presence and absence of Ca2+/CaM and the phosphorylation of RLC serine 19 was probed.

10 nM smMLCK was mixed with 10 µM RLC, 500 µM ATP in 20 mM HEPES (pH 7.5), 50 mM NaCl, 10 mM MgCl, 1 mM DTT, 500 µM ATP in the presence of either 40 nM CaM/1 mM CaCl or 1 mM EGTA and incubated at 20°C. Samples were taken at various intervals with the reaction quenched by addition of SDS-PAGE loading buffer.

Samples containing 300 ng RLC were run on SDS-PAGE alongside BioRad Precision Plus protein marker, transferred onto nitrocellulose membrane, incubated in 5% milk at room temperature for 45 min and probed with Cell Signaling Technology antibody 3671 against phospo-myosin light chain 2 (Ser19) in 5% milk for 2 hr at room temperature. The membrane was washed with low-salt buffer (10 mM Tris (pH 7.4), 150 mM NaCl, 0.1% Tween-20) three times and then incubated with horseradish peroxidase conjugated anti-mouse IgG (Dako P0260) in 5% milk for 45 min. After washing three times in low-salt buffer the membrane was stained using the enhanced chemiluminescence method.

Sample preparation

Glass coverslip and AFM cantilever were identically passivated for unspecific interactions using heterobifunctional succinimide-PEG-maleimide spacers (Rapp Polymere, Tübingen, Germany) with a molecular weight of 5000 Da (Celik and Moy, 2012). The succinimide group is attached via (3-aminopropyl)-dimethyl-ethoxysilane (APDMES, Karlsruhe, Germany) to the respective surface. The reactive maleimide group covalently conjugates to accessible thiol groups on applied bio-molecules used for specific immobilization. On the coverslip, this reaction is employed for covalent attachment of coenzyme A. The smMLCK construct harbors a C-terminal ybbR-tag (Wong et al., 2008) that site-selectively reacts to coenzyme A catalysed by the Sfp-synthase system (Yin et al., 2005, 2006). The investigated construct additionally contains an N-terminal Strep-tag II that is pulled via a monovalent variant of Strep-Tactin in the force spectroscopy experiments (Baumann et al., 2016). The monovalent Strep-Tactin is engineered to contain a unique cysteine residue on its single functional subunit, which is utilized for specific immobilization onto maleimide-PEG functionalized cantilevers (Biolever Mini, Olympus, Tokyo, Japan) (Zimmermann et al., 2010). For the AFM experiments, 40 mM HEPES (pH 7.2) with 2 mM MgCl and 1 mM DTT was used as measurement buffer.

Force spectroscopy experiments

AFM force spectroscopy data was acquired on a custom-built AFM operated in closed loop mode by a MFP3D controller (Asylum Research, Santa Barbara, CA, USA). Software for the automated control of AFM head and xy-piezos was programmed in IgorPro6 (Wavemetrics, OR, USA). Strep-Tactin coated Biolever Mini cantilevers were briefly brought in contact with the sample surface and then retracted at 800 nm/s. After each recorded force-distance curve, the surface was horizontally moved in steps of 100 nm distance for iteratively probing a new position. The cantilever spring constant was calibrated using the equipartition theorem method (Hutter and Bechhoefer, 1993). Typically, datasets contain 30000 force-distance curves and the addition of the substrate was performed in the course of the experiment. If data was not collected within one experiment, but was directly compared via rupture force histograms, the recorded force values were normalized according to the most probable rupture force of the Fn3 domain. For this normalization, only those force-distance curves in a dataset were regarded that detached after unfolding of the Fn3 domain.

Data analysis

Force heatmaps were assembled from all curves showing the characteristic unfolding pattern of the stressed smMLCK construct. Respective force spectroscopy data was aligned in force-distance space and transformed to a heatmap based on raw data points with 750 bins per axis. Data evaluation was carried out in Python 2.7. The rupture forces were evaluated from the AFM force-distance curves utilizing a quantum mechanically corrected WLC model (Hugel et al., 2005). The AFM distance was corrected for cantilever bending. 20 nm force baseline after the last rupture event – typically representing detachment – was used for determining zero force in the transformation of deflection signal to force values. Force-extension data was transformed into contour length space via an inverse worm-like chain model assuming a persistence length of 0.4 nm and a thermal energy of 4.1 pN nm (Jobst et al., 2013). On the transformed data a Gaussian kernel density estimate is applied with a bandwidth of 1 nm. Data set alignments in contour length space are created by the following process: the full set of transformed force-distance curves is aligned to a random curve from this data set according to least residual in cross correlation. This process results in a first superimposition which is used as a template in a second iteration of this process. Again, all contour-length transformed curves are aligned to a template curve but this time to the one formed by the first iteration. This two-step approach diminishes biasing effects given by the choice of the random curve used for initial alignment. Contour lengths of the individual domains are determined by a Gaussian fit of each determined peak and subtraction of the respective fitted means. The error of an increment is given by the standard deviations of both peaks defining the individual increment. Force-distance curves were denoised with Total Variation Denoising in order to detect rupture events as significant drops in force. For the force histograms, detected peaks in the typical region of respective domain unfolding were assigned to S→S, S→S or Fn3. All assignments were manually double-checked - especially to reassure small rupture events close to the noise level or to delete erroneously assigned peaks. In some curves no distinct rupture event could be detected for a specific domain due to too small forces (below noise level of about 10–15 pN) and was not included to the evaluation. Rupture forces of respective domains in the unfolding pattern were binned to histograms and fitted with the Bell-Evans model yielding the most probable rupture force (Bell, 1978; Evans and Ritchie, 1997). In the case of the S→S rupture event, kernel-density estimates with a bandwidth of 1 pN were applied to the data and used for extracting the most probable rupture force.

Amino acid sequence: SII-smMLCK-ybbR-His

MASWSHPQFEKGAETAVPNSAPAFKQKLQDVHVAEGKKLLLQCQVSSDPPATIIWTLNGKTLKTTKFII

LSQEGSLCSVSIEKALPEDRGLYKCVAKNDAGQAECSCQVTVDDAPASENTKAPEMKSRRPKSSLPPVLG

TESDATVKKKPAPKTPPKAAMPPQIIQFPEDQKVRAGESVELFGKVTGTQPITCTWMKFRKQIQESEHMK

VENSENGSKLTILAARQEHCGCYTLLVENKLGSRQAQVNLTVVDKPDPPAGTPCASDIRSSSLTLSWYG

SSYDGGSAVQSYSIEIWDSANKTWKELATCRSTSFNVQDLLPDHEYKFRVRAINVYGTSEPSQESE

LTTVGEKPEEPKDEVEVSDDDEKEPEVDYRTVTINTEQKVSDFYDIEERLGSGKFGQVFRLVEKKTRK

VWAGKFFKAYSAKEKENIRQEISIMNCLHHPKLVQCVDAFEEKANIVMVLEIVSGGELFERIIDEDFELTE

RECIKYMRQISEGVEYIHKQGIVHLDLKPENIMCVNKTGTRIKLIDFGLARRLENAGSLKVLFGTPEFVAPEVI

NYEPIGYATDMWSIGVICYILVSGLSPFMGDNDNETLANVTSATWDFDDEAFDEISDDAKDFISNLLKKD

MKNRLDCTQCLQHPWLMKDTKNMEAKKLSKDRMKKYMARRKWQKTGNAVRAIGRLSSMAMISGLSG

RKSSTGSPTSPLNAEKLESEEDVSQAFLEAVAEEKPHVKPYFSKTIRDLEVVEGSAARFDCKIEGYPDPE

VVWFKDDQSIRESRHFQIDYDEDGNCSLIISDVCGDDDAKYTCKAVNSLGEATCTAELIVETMEEPSSG

SGGDSLEFIASKLASGLRGSHHHHHH

In the interests of transparency, eLife includes the editorial decision letter and accompanying author responses. A lightly edited version of the letter sent to the authors after peer review is shown, indicating the most substantive concerns; minor comments are not usually included.

Thank you for submitting your article "Mechanical Force as a Functional Regulator in Smooth Muscle Myosin Light Chain Kinase" for consideration by eLife. Your article has been favorably evaluated by Arup Chakraborty (Senior Editor) and three reviewers, one of whom is a member of our Board of Reviewing Editors. The following individual involved in review of your submission has agreed to reveal his identity: James R Sellers (Reviewer #3).

The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.

Summary:

The smooth muscle myosin kinase (smMLCK) is a catalytically active serine/threonine kinase which is triggered by a conformational change upon Ca2+/calmodulin (Ca2+/CaM) binding. smMLCK shows a structural homology with the force-activated titin kinase and thus may be expected also to be regulated by mechanical pulling stress. In this paper, the authors studied the mechanical response of smMLCK applying AFM-based single-molecule force spectroscopy (SMFS) to understand how this cytoskeletal kinase is conformationally modulated by external forces. It was particularly investigated whether the presence of different ligands including ATP, Ca2+/CaM or RLC peptide substrate or a combination of these affect the mechanical unfolding intermediates along the unfolding pathway of smMLCK. To reveal statistically solid insights these mechanisms have been investigated via high throughput SMFS. The characteristic pattern of transitions between different smMLCK states and their variations in the presence of different substrates and ligands are presented. The results show a very interesting interaction between kinase domain and regulatory light chain (RLC) substrate in the absence of CaM, which indicates a restored substrate-binding capability due to mechanically induced removal of the auto-inhibitory regulatory region. The data also suggest that the conformational change occurs in the presence of the substrate, suggesting that the substrate cannot bind prior to the said conformational change. In addition, the increase in unfolding force shows that the substrate does bind after the conformational change. The AFM-based single-molecule force spectroscopy experiments are conducted at outstanding quality, the experimental data is carefully analyzed and interpreted. Suitable controls have been performed and support the experimental conclusions very convincingly. The paper is written very well.

Essential revisions:

The evidence for mechanical regulation of kinase activity is indirect and thus the title may be a bit presumptuous. We realize that this is a very difficult thing to actually prove and the necessary technical advances are brought up in the Discussion. Based on this, we wonder whether the title oversells the manuscript.

Some of the speculation about possible mechanoregulation of the kinase stems from early biochemical studies that showed smooth muscle MLCK could bind myosin filaments via its C-terminus and actin filaments via the N-terminus. However, the K for both of these bindings were in the μm range and the recent results from the Cremo and Baker labs (Hong et al., 2015) show that the kinase diffuses along actin filaments. Based on this how would the authors imagine that the kinase could be subjected to sufficient force to pull out the autoinhibited domain rather than just dissociating from actin filaments or myosin filaments if the two systems were sliding relative to one another? Perhaps catch bonds are involved?

Reviewer #1:

In this work, they present evidence that the inhibitory domain of the kinase undergoes a tension-induced conformational change that is very similar to the change that occurs upon calmodulin binding, which is a known activator. Because the N and C termini are known to bind actin filaments and xxx, respectively, it is proposed that tension may function as an additional regulator mechanism for kinase activation. The data also suggest that the conformational change occurs in the presence of the substrate, suggesting that the substrate cannot bind prior to the said conformational change. In addition, the increase in unfolding force shows that the substrate does bind after the conformational change. Overall, I find the work interesting and of high quality, and there is a good reason to suspect that force can activate the enzyme activity.

Reviewer #2:

The smooth muscle myosin kinase (smMLCK) is a catalytically active serine/therosine kinase which is triggered by a conformational change upon Ca2+/calmodulin (Ca2+/CaM) binding. smMLCK shows a structural homology with the force-activated titin kinase and thus may be expected also to be regulated by mechanical pulling stress. In this paper, the authors studied the mechanical response of smMLCK applying AFM-based single-molecule force spectroscopy (SMFS) to understand how this cytoskeletal kinase is conformationally modulated by external forces. It was particularly investigated whether the presence of different ligands including ATP, Ca2+/CaM or RLC peptide substrate or a combination of these affect the mechanical unfolding intermediates along the unfolding pathway of smMLCK. To reveal statistically solid insights these mechanisms have been investigated via high throughput SMFS. The characteristic pattern of transitions between different smMLCK states and their variations in the presence of different substrates and ligands are presented. The results show a very interesting interaction between kinase domain and regulatory light chain (RLC) substrate in the absence of CaM, which indicates a restored substrate-binding capability due to mechanically induced removal of the auto-inhibitory regulatory region. The AFM-based SMFS experiments are conducted at outstanding quality, the experimental data is carefully analyzed and interpreted. Suitable controls have been performed and support the experimental conclusions very convincingly. The paper is written very well. I strongly recommend publication.

Reviewer #3:

The manuscript by Baumann et al. uses AFM to unfold smooth muscle myosin light chain kinase. From the data they infer that force may be a regulator of kinase function. The data show that various domains of the kinase unfold sequentially in a reproducible manner and the authors are able to assign unfolding events to particular domains. They show that the kinase domain unfolds first followed by the FN3 domain and then the Ig domains unfold last. The evidence for mechanical regulation of kinase activity is very indirect and thus the title may be a bit presumptuous. The main data to suggests this behavior is that the S2 state is stabilized in the presence of an RLC substrate peptide.

The study is of considerable interest and is well executed.

Regarding my comment about the data not really showing mechanical activation of the enzyme, I realize that this is a very difficult thing to actually prove and the necessary technical advances are brought up in the Discussion. Based on this, I do wonder whether the title oversells the manuscript.

Some of the speculation about possible mechanoregulation of the kinase stems from early biochemical studies that showed smooth muscle MLCK could bind myosin filaments via its C-terminus and actin filaments via the N-terminus. However, the K for both of these bindings were in the μm range and the recent results from the Cremo and Baker labs (Hong et al., 2015) show that the kinase diffuses along actin filaments. Based on this how would the authors imagine that the kinase could be subjected to sufficient force to pull out the autoinhibited domain rather than just dissociating from actin filaments or myosin filaments if the two systems were sliding relative to one another?

Also in the cartoon in Figure 5, the authors show the while the RLC stabilizes the lower portion of the kinase domain, the upper portion (depicted as the ATP binding site) is already unfolded. This does not seem to be a mechanism for the mechanical regulation of activity.

The legend to Figure 5 does not describe the significance of the color change of the triangular motif (CaM). Does the light gray color mean that this is apo-calmodulin?

Essential revisions:

The evidence for mechanical regulation of kinase activity is indirect and thus the title may be a bit presumptuous. We realize that this is a very difficult thing to actually prove and the necessary technical advances are brought up in the Discussion. Based on this, we wonder whether the title oversells the manuscript.

The shown results give strong experimental indication for mechanical force being a functional regulator in smMLCK comparable to the established Ca2+/CaM activation pathway. However, we agree that those evidences are drawn rather indirectly by means of AFM force spectroscopy data. Although the data is conclusive and in our opinion convincing, a direct proof of the force activation cannot be given so far. To make this clear to the reader, we would therefore like to suggest the following title:

Increasing Evidence of Mechanical Force as a Functional Regulator in Smooth Muscle Myosin Light Chain Kinase (see title).

Some of the speculation about possible mechanoregulation of the kinase stems from early biochemical studies that showed smooth muscle MLCK could bind myosin filaments via its C-terminus and actin filaments via the N-terminus. However, the K for both of these bindings were in the μm range and the recent results from the Cremo and Baker labs (Hong et al., 2015) show that the kinase diffuses along actin filaments. Based on this how would the authors imagine that the kinase could be subjected to sufficient force to pull out the autoinhibited domain rather than just dissociating from actin filaments or myosin filaments if the two systems were sliding relative to one another? Perhaps catch bonds are involved?

Hong et al. indeed show sliding of actin on a myosin-coated surface with MLCK having µM affinity to myosin and actin. However, the measured affinity does not necessarily implicate low binding forces as could be shown in the past by other force spectroscopy studies. Other factors such as pulling geometry or force trajectory through the receptor-ligand system could be identified as much more dominant parameters. Best example is the Strep-Tag II:Strep-Tactin handle system employed in this study itself: despite a similar K in the µM range, we could show that rupture forces easily suffice to not only activate but also unfold the kinase and its side domains.

A crucial factor, however, is given by the timescales of the binding partners represented in the dissociation rate. The system has to form a stable linkage lasting long enough to bridge mechanical stress and therefore allow activation of the kinase domain before it dissociates again. For a system diffusing along the filaments, however, long-term association can certainly be assumed. Additional binding domains in the long isoforms of MLCK additionally increase this cytoskeletal association as shown by Kudryashov et al. in 2004 (Myosin light chain kinase (210 kDa) is a potential cytoskeleton integrator through its unique N-terminal domain).

Our hypothesis is that relative movement of actin to myosin in the muscle sets smMLCK under tension. Even for physiological forces below values of 15 pN as shown in our AFM force spectroscopy data, activation can be expected if the molecule is kept in a stretched conformation for sufficient time (comparable to a force-clamp configuration). A catch-bond mechanism might be possible but is not essential for our interpretation.

To address this we included a paragraph in the discussion section of the manuscript.