Single-molecule analysis reveals that regulatory light chains fine-tune skeletal myosin II function

Myosin II is the main force-generating motor during muscle contraction. Myosin II exists as different isoforms that are involved in diverse physiological functions. One outstanding question is whether the myosin heavy chain (MHC) isoforms alone account for these distinct physiological properties. Unique sets of essential and regulatory light chains (RLCs) are known to assemble with specific MHCs, raising the intriguing possibility that light chains contribute to specialized myosin functions. Here, we asked whether different RLCs contribute to this functional diversification. To this end, we generated chimeric motors by reconstituting the MHC fast isoform (MyHC-IId) and slow isoform (MHC-I) with different light-chain variants. As a result of the RLC swapping, actin filament sliding velocity increased by ∼10-fold for the slow myosin and decreased by >3-fold for the fast myosin. Results from ensemble molecule solution kinetics and single-molecule optical trapping measurements provided in-depth insights into altered chemo-mechanical properties of the myosin motors that affect the sliding speed. Notably, we found that the mechanical output of both slow and fast myosins is sensitive to the RLC isoform. We therefore propose that RLCs are crucial for fine-tuning the myosin function.

Myosin motors drive diverse motile processes, ranging from intracellular cargo transport and cell division to muscle contraction and the whole cell movement. Striated muscle myosin, such as skeletal or cardiac myosin II, responsible for generating the force during muscle contraction is a hexameric motor composed of two heavy chains (myosin heavy chains (MHCs) 6 ) and four light chains (LCs) (1). Each heavy chain contains a globular motor domain, an ␣-helical light-chain binding domain, and a long coiled-coil rod. The rod parts from different myosin molecules self-assemble to form bipolar thick filaments. The ϳ9-nm-long ␣-helix, also termed as the "lever arm," serves as a link between the motor domain and the rod part of heavy chains. Each MHC has an essential (ELC) and regulatory light chain (RLC) wrapped around the lever arm (see Fig. 1A). The main motor domain has ATP-and actin-binding sites. During the chemo-mechanical coupling, small conformational change in the catalytic domain linked with the ATP hydrolysis gets amplified as a large movement of the lever arm (2,3). The light chains are reported to be essential to maintain the rigidity of the lever arm as force is transduced during the power stroke (4).
Muscle myosin II exists as different isoforms, three "fast" myosin heavy chain (MHC) isoforms (i.e. MHC-IIa, -IId/IIx, and -IIb) and one "slow" MHC-I isoform (5). The MHC-I isoform or ␤ cardiac isoform is primarily expressed in slow skeletal, aerobic muscle fibers and in heart ventricle (6). The muscle expressing "fast" MHCs expresses myosin light-chain (MLC) isoforms (i.e. different ratios of ELC isoforms, MLC1f and MLC3f, and RLC isoform, MLC2f). The "slow" MHC-I isoform is typically equipped with the ELC isoform, MLC1s, and the RLC isoform, MLC2s/MLC2v. RLC or MLC2 (ϳ19 kDa), encoded by the MYL2/Myl2 genes, is a member of the EF-hand superfamily of Ca 2ϩ -binding proteins with a helix-loop-helix structural motif. MLC2 noncovalently binds to the distal end of the lever arm at the junction between the S1 and S2 regions of myosin. Whereas the C-terminal domain of RLC binds in the region between Glu-808 and Val-826 of myosin lever arm, the N-terminal domain of RLC surrounds the region between amino acids Asn-825 and Leu-842 (2). Apart from its structural role to provide rigidity, RLCs are known to regulate cardiac, smooth, and skeletal muscle contraction via phosphorylation. For example, genetic loss-offunction studies in mice revealed an essential function of RLCs in cardiac contraction (7). Compelling reports demonstrated further that RLC phosphorylation regulates cardiac muscle contraction by increasing the number of cross-bridges available to generate the force (8 -10). Reduction in phosphorylation levels of MLC2v are critically linked with human cardiomyopathies (11). Phosphorylation of MLC2 in smooth muscle myosin II determines the active state of the motor required for muscle contraction, thus making RLC a key regulator. However, for striated muscle myosin, the phosphorylation status of RLC is not correlated with the activation of the cross-bridge cycle, but rather Ca 2ϩ binding to troponin activating the thin filaments is the prerequisite to initiate the muscle contraction. MLC2 phosphorylation is known to increase the steady-state force in permeabilized muscle fibers by increasing the number of crossbridges interacting with actin filaments at submaximal levels of Ca 2ϩ activation, without affecting maximum shortening velocity (12,13). Thus, RLCs exert differential effects in a tissuespecific manner. Collectively, MLC-2 is thus considered to play a mechanical role by stabilizing the lever arm and a regulatory role through phosphorylation, but not to affect the ATPase kinetics.
Muscle fiber studies indicated that the cross-bridge kinetics underlying the force transients are determined mainly by the MHC isoforms (14). While the MHCs as the main determinant of the mechanical output is well-agreed on, MLC's primary role in skeletal myosin II is considered to be rather structural. The correlation of the muscle fiber-specific expression and assembly of the light chains with specific MHCs, however, might imply an additional role for the regulatory light chains. Two major observations-1) the distinct sets of light chains associating with different isoforms of myosin heavy chain and 2) single point mutations in the ventricular RLC causing heart disorder in humans, such as hypertrophic and dilated cardiomyopathies-directed us to investigate the importance of RLC in mechanical performance of the myosins.
This investigation is also relevant to understand the mechanical performance of the hybrid/chimeric motors (i.e. combination of specific myosin heavy chains with different light chains) found under different physiological conditions. In some cases, heterogeneity in MHC and LC isoform expression in single muscle fibers has been documented in vertebrate skeletal mus-cles. For instance, extensive analysis of the SDS-polyacrylamide gels from single soleus muscle fibers revealed that in a majority of fibers, the slow heavy chain and corresponding light chains were expressed. In some soleus single fibers, however, fast lightchain isoforms were expressed, which potentially may form a complex with the slow heavy chains (15). In human fast single fibers, co-existence of both fast (MLC2f) and slow (MLC2s) light chains was observed (16). Moreover, analyses of single fiber segments revealed coexpression of fast-and slow-myosin components nonuniformly distributed along the fiber length (17)(18)(19). During the developmental phase, the concentration of such hybrid/chimeric motors is reported to be increased (19). Thus, a mixture of fast and slow motor components (i.e. heavy and light chains) indeed exists in single muscle fibers; however, the mechanical properties of such hybrid motors remain unexplored.
Here, we examined the myosin chimeras where the endogenous regulatory light chain was replaced with a different RLC variant to probe whether, how, and which biochemical and mechanical properties are affected as a result of different trimeric formation(s) of myosin II motor complex (i.e. MHC, ELC, and RLC). We found that the motor properties of a given myosin isoform were sensitive to the type of RLCs. A chimeric motor based on fast myosin heavy chain translocated the actin filaments with reduced speed; conversely the slow myosin improved the gliding speed, simply by replacing native RLC with a different RLC variant. Steady-state and transient ATPase kinetics measurements allowed us to relate the altered motor properties to changes in the overall ATP turnover time and rates of product release. Single-molecule analyses based on optical trapping measurements revealed pronounced alterations in the duration of the strong actin-bound states, size of the power stroke, and stiffness of motors. These studies unraveled the critical role of different RLCs in muscle myosin motor function.

Reconstitution of chimeric myosin II motors and analysis of actin filament gliding
We investigated the role of the RLCs in modulating the motor function of skeletal myosin II isoform found in the fastand slow-twitch muscle fibers. We isolated myosin II motors from rabbit fast-twitch (musculus psoas) and slow-twitch (musculus soleus) muscle fibers. We inspected the native fast and slow myosin extracted from respective muscles in the SDSpolyacrylamide gels to ensure the purity (Fig. S1). We generated single-headed subfragment 1 (S1) and reconstituted the chimeric motors based on fast (MHC-IId) and slow (MHC-I) myosin II. In its native form, MHC-IId, henceforth referred as WT fast (WT-S1f) myosin S1, associates with the 2 light chains, i.e., ELC, MLC1f/MLC3f, and the RLC, MLC2B. MHC-I or WT slow (WT-S1s) myosin S1 is typically equipped with MLC1v and MLC2v.
As shown in Fig. 2A Figure 1. Myosin S1 structure and ATPase cycle. A, crystal structure of scallop myosin II-subfragment-1 (Protein Data Bank code 1SR6) (62). The myosin heavy and light chains are shown as a ribbon diagram; S1 is indicated with the key parts of the motor domain. B, kinetic scheme for actomyosin ATPase. Myosin II with different nucleotide states is illustrated: A, actin; M, myosin; T, ATP; D, ADP; Pi, P i . Strong and weak interaction states of myosin with actin are labeled in black or green, respectively. Rate constants for forward and reverse reactions are indicated as k ϩn and k Ϫn , respectively.

Regulation of myosin II via light chain
ric motors based on the fast and slow myosins. The cgmRLC was previously reconstituted with the chicken skeletal myosin II and was shown to retain the high gliding velocity of the motor (20). For smooth muscle myosin (SMM) lower gliding speeds from 0.2 to ϳ1 m/s were reported under different motility conditions (21)(22)(23). SMM displayed kinetics properties similar to those of slow skeletal soleus myosin (6), with basal Mg 2ϩ ATPase activity of 0.05 s Ϫ1 , acto-S1 ATPase V max of 0.7 s Ϫ1 , ADP dissociation of 15 s Ϫ1 (at 20°C), and a second order rate constant for ATP binding of 0.5 M Ϫ1 s Ϫ1 (24). Despite these SMM properties, when it was assembled with chicken skeletal muscle myosin, the actin filament gliding speed of the myosin subfragment 1 was increased by ϳ4-fold, similar to the native full-length skeletal myosin (20).
Therefore, we used it as a control to test its effect on the fast and slow skeletal myosin under our experimental conditions. As depicted in the Fig. 2A, we replaced the endogenous RLCs of WT myosins and generated six chimeric forms. The fast myosin chimeras comprised, 1) S1f-MLC2B (i.e., fast MHC with skeletal RLC), 2) S1f-MLC2v (i.e. fast MHC with ventricular RLC), and 3) S1f-cgmRLC (i.e. fast MHC with smooth mucle RLC). The slow myosin chimeras included, 1) S1s-MLC2v (i.e. slow MHC with ventricular RLC), 2) S1s-cgmRLC (i.e. slow MHC with smooth mucle RLC), and 3) S1s-MLC2B (i.e. slow MHC with skeletal RLC). We employed human variants of the RLCs because of the nearly identical protein sequences as shown in Fig. S2. Human MLC2B (MLC2B) and rabbit MLC2f (rMLC2f) share 98% sequence similarity. Likewise, human MLC2v is 98% similar to the rabbit slow RLC isoform. However, the fast isoforms MLC2B and rMLC2f has an 86% sequence similarity with the slow MLC2v. The cgmRLC shares low sequence similarity with both fast MLC2B and slow MLC2v (i.e. only 67%).
To understand the main effect of RLCs in these different combinations of chimeric motors, we focused on the mechanical parameters at the ensemble level using an in vitro actin filament gliding assay. We characterized two WT and six chimeric motors in in vitro actin filament gliding assays performed at 2 mM ATP concentrations (and 50 -100 g/ml surface density of motors). The motors were immobilized on a nitrocellulose-coated surface. As shown in Fig. 2B and Table 2, WT-S1f and WT-S1s showed Ͼ20-fold difference in the gliding speed (i.e. 1.25 and 0.047 m/s, respectively) (see Movie S1). Very interestingly, we found gliding velocities on both fast MHC-IId and slow MHC-I motor surface to be sensitive to the variant of RLCs bound to the lever arm. The S1f-MLC2v combination resulted in a decrease in the actin filament gliding velocity by ϳ3-fold, whereas S1f-cgmRLC and S1f-MLC2B retained the high actin filament gliding speed similar to the WT-S1f. On the contrary, the S1s-cgmRLC chimeric motor showed an ϳ10fold increase in the gliding speed of the slow S1s motor from 0.047 to 0.43 m/s. The S1s-MLC2B motor showed a marginal increase in the gliding speed compared with WT S1s. Interestingly, as shown in Movie S2 and Fig. 2B, fast and slow myosindriven actin filaments were indistinguishable when fast MHC was combined with MLC2v (S1f-MLC2v) and slow MHC with cgmRLC (S1s-cgmRLC).
These experiments showed that the motile activity of skeletal muscle myosin II is critically influenced by the regulatory lightchain isoform. Replacement of the native RLCs from fast and slow myosin II S1 with the RLC variants, MLC2B and MLC2v, respectively, showed comparable sliding velocities ( Fig. 2B and Table 2). These results ruled out the alteration of motor properties due to the exchange procedure and further emphasized RLC-specific effects on myosin function.
The altered sliding velocities indicated that RLC interferes with the actomyosin (AM) chemo-mechanical coupling. Actin filament sliding velocities are governed by the net rate of actomyosin cross-bridge cycling (25, 26). The parameters that influence the sliding speed (V) include (i) the lifetime of the strong A)

Regulation of myosin II via light chain
AM-bound state (t on ) and (ii) the stroke size (d) of myosin (i.e. V ϭ d/t on ). The duration of AM-bound and -unbound states during the ATPase cycle described by duty ratio of the myosin is the fraction of the total ATPase cycle time myosin spends bound to actin. The duty ratio can change either by a change in the ADP release rate from the AM complex or the rate of weakto-strong bound transition of AM (i.e. A.M. D.P i to A.M.D state) (see the ATPase scheme in Fig. 1). One other parameter that must be considered in ensemble molecule sliding velocity measurements is the surface density of the myosin molecules. The skeletal myosin II, a low-duty ratio motor (Ͻ0.05), however, showed low dependence of velocities on the motor density (27,28). Therefore, we presumed that the velocities most likely changed due to altered t on and/or d.
For native S1f, while some reports showed an isomerization step before P i release as a rate-limiting step, the others suggested P i release as the slowest step in the AM ATPase cycle (29,30) and ADP release to be fast (Ͼ500 s Ϫ1 ). For the WT S1s, however, ADP was shown to have Ͼ10-fold higher affinity and thus slower ADP release during the cross-bridge cycle, ϳ20 s Ϫ1 (6). Therefore, the detachment rate of the AM complex at the post-powerstroke state is limited by either the rate of ADP release or ATP binding to the rigor AM complex for these slow and fast myosin isoforms.

Ensemble molecule steady-state ATPase and transient kinetics
To investigate the cause of the altered sliding velocities, we measured actin-activated ATPase activities and rates of product release to determine potential rate-limiting steps in the ATPase cycle (Fig. 3). Consistent with previous studies (6), we obtained comparable basal ATPase activities between WT-S1f (0.04 s Ϫ1 ) and S1s (0.05 s Ϫ1 ) and almost 100-fold (3.5 s Ϫ1 ) and 10-fold (0.5 s Ϫ1 ) increase in the actin-activated ATPase activities, respectively ( Fig. 3 (A and B) and Table 1). S1f-cgmRLC displayed comparable actin-activated ATPase activity as the WT construct (Table 1), whereas S1f-MLC2v showed reduced actin activation, as reported previously (31) ( Table 1). For the slow myosin chimeras, we observed pronounced differences in the basal ATPase activity by almost a factor of 10 between S1s (0.05 s Ϫ1 ) and S1s-cgmRLC (0.55 s Ϫ1 ). The high basal ATPase activity of S1s-cgmRLC was unexpected but was reproducibly observed with different preparations. Such comparably high basal activities have been observed in mutant myosins (32,33) and also unconventional myosins, myosin IX (34) and myosin VIIa (35). The S1s-MLC2B construct showed a similar basal activity as the WT S1s (0.04 s Ϫ1 ). The ATPase activity increased for all constructs with increasing actin concentrations, yielding values for activation at 60 M F-actin of 0.14, 0.62, and 0.09 s Ϫ1 for S1s, S1s-cgmRLC, and S1s-MLC2B, respectively. Notably, the efficiency by which actin activated the ATPase (k cat /K app ) was changed 3-5-fold among the motors (Table 1). In additional transient kinetic experiments, we measured the rates of ADP and P i release as major determining factors of ATPase and sliding velocity. Generally, ADP release from actin-bound S1f is a very fast process with an estimated rate constant Ͼ500 s Ϫ1 (6). As seen in Fig. 3C, ADP-release kinetics was too fast to be resolved or not accompanied by signal change. The amplitude change was limited by the dead time (0.5 ms) of the instrument, providing only estimates for the rates of ADP release from actomyosin of Ͼ500 s Ϫ1 for all constructs (Table 1). Notably, ADP release from actin-bound S1s could be monitored due to its slower release kinetics and high ADP affinity (6). We observed two-step ADP-release kinetics for S1s, S1s-cgmRLC, and S1s-MLC2B, yielding rate constants of kЈ ϪAD ϭ 54, 61, and 60 s Ϫ1 for the fast phase and kЈ ϪAD ϭ 3.9, 4.9, and 2.9 s Ϫ1 for the slow phase, respectively (Fig. 3D). We observed P i release as the ratelimiting step of the ATPase reaction for all measured constructs (Fig. 3, E and F). This is reflected in the similarity between P i release rates and actin-activated ATPase rate at the 20 M actin concentration ( Fig. 3 (A and B) and Table 1). Thus, the RLCs appear to dictate the overall cycling time, mainly by affecting ADP release and modulating the P i release kinetics.

Single-molecule analysis
To further probe the change in velocities with chimeric motors, we employed a single-molecule analysis method to gain precise insights into biochemical properties, such as the duration of AM strong bound states, and biophysical properties, such as stroke size (d) and stiffness of native and chimeric motor forms.

Single chimeric motor molecule's kinetic properties
We employed three-bead trapping assays for in-depth analysis of the biochemical and mechanical properties that led to either improved or decreased actin filament gliding speed driven by chimeric motors. Optical trap experimental setup has been described in detail previously (36). Briefly, the actin filament is embedded and stretched between the two optically trapped beads forming a dumbbell, and the motor is immobilized on a glass bead attached on the chamber surface (Fig. 4A). The positions of the two trapped beads are precisely monitored using a quadrant photodiode and acquired.
Single-headed WT or chimeric motor molecules were analyzed for their intermittent interaction with the actin filaments as shown in Fig. 4B at 10 M ATP concentration. We determined the lifetimes of the motor interaction with actin. The reduced Brownian noise in the data traces was characteristic of the AM-bound states, t on , whereas free dumbbell noise, t off , represented myosin-detached states. It is well-known that actin accelerates the P i release, and P i release from the myosin active site is closely associated with the powerstroke generation. However, there is no clear consensus on whether P i is released before, during, or after the powerstroke. In our experimental set up, the initiation of actomyosin association seen as the reduction in the noise amplitude of the bead coincides with the powerstroke generation.
Thus, t on comprised the post-powerstroke states with ADP present in the active site (AM.ADP) plus the time the active site is free of nucleotide (rigor state) until a new ATP molecule binds and initiates the rapid detachment of the myosin head from the actin filament. This argument is valid if the P i is released before or during the powerstroke. Some recent studies have supported the model that the powerstroke precedes the P i release from the active site (37,38). However, to consider the AM.ADP.P i state as a "strong bound state"-contributing to the lifetime of an attached event-the duration of this state should be sufficiently long after the powerstroke (at least 2 ms

Regulation of myosin II via light chain
(i.e. the detection limit with our current experimental setup)). We cannot distinguish between strong bound AM.ADP.P i and AM.ADP states.
The AM-detached state, t off , included the ATP hydrolysis and M.ADP.P i state prior to association with actin (see the ATPase scheme in Fig. 1). As the ADP release from the active Plots of ATP turnover rates as a function of F-actin concentration for fast (A) and slow (B) myosins. Data were fitted to hyperbolas using the Michaelis-Menten formula. The catalytic efficiency, k cat /K app , was obtained from the initial slope of the hyperbola. C and D, ADP-release kinetics from actomyosin. The rate of ADP dissociation from the actomyosin complex was measured by rapidly mixing actomyosin-mantADP complex with excess ADP. The decrease in mant fluorescence followed single exponential kinetics in the case of the fast myosins (C). For the slow myosin constructs, the fluorescence transients were best described by two exponentials with a fast rate constant (k ϪAD fast ) and a slow rate constant (k ϪAD slow ) indicative of two-step ADP-release kinetics (D). ATP turnover experiments of actomyosin constructs used coumarin-labeled phosphate-binding protein (MDCC-PBP; Thermo Fisher Scientific) as a sensor for the detection of liberated P i from fast (E) and slow (F) myosins. The transients follow single exponentials with a corresponding rate constant for P i release (k ϪPi ). All kinetic parameters are summarized in Table 1.

Table 1
Summary of kinetic parameters k basal , ATPase activity in the absence of actin; v max , maximum actin-activated ATPase activity; K app , concentration of actin at one-half v max . v max and K app were estimated from hyperbolic fits. v max (60 M), actin-activated ATPase rate at 60 M F-actin; k cat /K app , catalytic efficiency; k ϪAD fast and k ϪAD slow , rate constants of the ADP release from acto-myosin; k ϪPi , rate constant of P i release at 20 M F-actin concentration.

Regulation of myosin II via light chain
site is in the range of Ͼ500 s Ϫ1 , we believe that for the fast myosins, primarily an AM-bound rigor state is acquired. At saturating ATP concentrations, the duration of the rigor state becomes too short and is essentially undetectable for fast myosins. Therefore, for comparison among different native and chimeric motors, we measured lifetimes, t on , at 10 M ATP concentrations. Fig. 4B shows the comparison of 7-s traces obtained with WT-S1f and WT-S1s, containing several distinct single AMbinding events at 10 M ATP. The obvious difference in the two traces with the two types of motors was the short versus long lifetime of AM-bound events. Consistent with previous observations, S1f showed shorter t on as compared with S1s (39). We estimated the average lifetime of AM-bound states by fitting the t on events from several molecules with a single exponential decay function (Fig. 4, C-J). Fig. 4K and Table 2 display the average lifetimes for WT and the chimeric motors. Duration of AM bound states for WT-S1f, S1f-cgmRLC, and S1f-MLC2B showed no significant difference. However, the S1f-MLC2v motor showed a Ͼ2-fold decrease in t on. For slow myosin, compared with WT-S1s, S1s-MLC2v retained the extended t on . Interestingly, with S1s-cgmRLC, the t on was significantly reduced by ϳ4-fold, and by ϳ3-fold for S1s-MLC2B. From the reciprocal of average time constants, we estimated the AM detachment rate constants and thereby the second-order rate constants of ATP binding (k ATP ) for fast myosins, k ATP ϭ 2.2 M Ϫ1 s Ϫ1 for WT-S1f, 3.26 M Ϫ1 s Ϫ1 for S1f-cgmRLC, 2.2 M Ϫ1 s Ϫ1 for S1f-MLC2B, and 5.9 M Ϫ1 s Ϫ1 for S1f-MLC2v. Our estimation of k ATP for WT-S1f is consistent with previous studies (40,41). Lifetime (ms) S1s-MLC2v S1s-MLC2B S1f-MLC2B S1f-cgmRLC WT-S1f

Regulation of myosin II via light chain
For all of the chimeric constructs except for one (i.e. S1f-MLC2v), the increase in the duration of the AM-bound state was associated with the decrease in the actin filament sliding velocity. The chimeric motors, S1f-cgmRLC and S1f-MLC2B, showing t on (or k ATP ) comparable with WT-S1f, also displayed similar gliding velocities (see Fig. 2B). A similar tendency was observed for WT-S1s and S1s-MLC2v.

Mechanical properties
Stroke size-Previously, single-molecule analysis of motors showed that the lever arm length affects the amplitude of displacement (i.e. stroke size, d) (42). As an accessory protein associated with the lever arm, whether light chain influenced stroke size was not explored.
We therefore investigated whether different variants of the light chains influenced the amplitude of stroke size of the chimeric motors. We determined the average stroke size from the shift of histogram method, introduced by Molloy et al. (43). As shown in Fig. 5 (A-I) and Table 2, we noted a significant difference in the mean displacement for fast and slow motors (i.e. from 4.04 Ϯ 0.41 to 6.11 Ϯ 0.16 nm, respectively). Thus, consistent with previous reports for rat fast and slow skeletal myosin II-S1 (39), the average working stroke for WT-S1f was smaller than for WT-S1s. We found a moderate yet significant increase in the stroke size for chimeric construct S1f-MLC2B as compared with its WT counterpart. For slow myosin heavy chain, however, the S1s-MLC2v and S1s-MLC2B combination resulted in a significant decrease of the average displacement amplitude, with S1s-MLC2B displaying merely 1.86-nm stroke size.
Stiffness-The light chains are reported to be responsible for the structural stability of the lever arm. The ϳ9-nm-long helical structure of lever arm that links the main motor core to the dimerization region and continues as a coiled-coil region that forms the basis of the thick filament backbone may contribute to the compliance of the motor.
One assumption is that variants of RLCs provide different stiffness to the myosin head. Consequently, different myosin head stiffness may result in a different mechanical strain. The kinetics of some steps of the ATPase cross-bridge cycle is assumed to respond to the change of strain under load (44,45).

Regulation of myosin II via light chain
The strain dependence of ATPase kinetics is well-reported for the motors, such as myosin V and smooth muscle myosin (24,46). Could the binding of MLC2v compared with MLC2B change, for example, the stiffness of myosin S1f and resultantly alter the AM cross-bridge cycle? To investigate this notion of whether the different RLCs influenced the stiffness of the chimeric motors, we acquired the AM interaction records by applying positive feedback on the laser-trapped beads as described by Steffen et al. (36,47). Accordingly, the signal/noise ratio (bound versus free dumbbell) was improved by increasing the noise amplitude of the free dumbbell, using analog positive feedback in AC mode. The trapped bead position analog signal from the quadrant detector was directed to the electro-optic deflectors, controlling the trap position. As the actomyosin binding event detection is more difficult for the motors, especially with low stiffness, this procedure allows the event detections and stiffness measurements by increasing the amplitude of Brownian noise by ϳ2-fold, while maintaining the trap stiffness in the y and z directions. Using a bead variance-covariance method as described (48,49), the stiffness of individual motors was calculated, as shown in Fig. 6. We found significant difference in the stiffness for WT-S1f versus WT-S1s (i.e. 1.12 Ϯ 0.03 and 0.437 Ϯ 0.08 pN/nm, respectively) (Fig. 6E). The stiffness values are compa-rable with previous optical trapping measurements from rat fast and slow myosin S1 (39, 50) and single fast and slow muscle fiber measurements for rabbit myosins (50). For the fast myosin, replacing the light chains with two different variants (MLC2B and MLC2v) did not alter the stiffness significantly, as shown in Fig. 6 and Table 2. Similarly, for slow myosin chimeras (S1s-cgmRLC and S1s-MLC2v), there was slight increase in the average stiffness, however not significantly different from the WT form. The S1s-MLC2B motor, however, showed the highest stiffness of 1.62 Ϯ 0.27 pN/nm.
To measure the stiffness of the WT-S1s and S1s-MLC2v chimera, the experiments were performed at higher ATP concentration (i.e. at 50 M ATP) to avoid the occasional inclusion of contaminant fast myosin in the measurements. Although the soleus muscle primarily contains (about 95%) slow myosin II isoform, the presence of ϳ5% fast myosin isoform cannot be avoided. In ensemble molecule gliding assays, 5% contamination of motors is unlikely to influence the actin filament speed as measured previously for mixtures of different ratios of myosin isoforms (22,23). Accordingly, more than 20% contamination of motors was required to cause alteration of gliding speed. In the single-molecule assays, however, this factor needs to be cautiously controlled. At 50 M ATP concentration, the AM interaction events for the fast myosin are too short for detection Variance was calculated for a rolling window of 20 ms, at our sampling rate of 10 kHz. For the example shown here, the variance-Hidden-Markov method yielded a combined trap stiffness of 0.078 pN/nm and a myosin head stiffness of 0.85 pN/nm. E, box plot with overlapped data points, the stiffness measured for WT-S1f, WT-S1s, and chimeric motors. Each rhombus represents the measured stiffness from individual myosins. Average stiffness with S.D. was as follows: for WT-S1f, Ϫ1.12 Ϯ 0.03; for S1f-MLC2B, Ϫ1.24 Ϯ 0.09; for S1f-MLC2v, Ϫ1 Ϯ 0.07; for WT-S1s, Ϫ0.437 Ϯ 0.08; for S1s-cgmRLC, Ϫ0.585 Ϯ 0.024; for S1s-MLC2B, Ϫ1.62 Ϯ 0.27; for S1s-MLC2v, Ϫ0.533 Ϯ 0.11. Unpaired t test was used to calculate the statistical significance. WT-S1f (n ϭ 20) and WT-S1s (n ϭ 13) displayed statistically significant difference with p Ͻ 0.0001. No statistically significant difference between WT-S1f and S1f-MLC2B (N ϭ 12, p ϭ 0.72) or WT-S1f and S1f-MLC2v (N ϭ 15, p ϭ 0.128) was found. For WT-S1s and S1s-cgmRLC (N ϭ 10), marginal but statistically significant difference was found (p ϭ 0.0493). For WT-S1s and S1s-MLC2B, N ϭ 6 and p Ͻ 0.0001. For WT-S1s and S1s-MLC2v, n ϭ 11 and p ϭ 0.4735. N ϭ number of individual motor molecules.

Regulation of myosin II via light chain
with our setup (Ͻ2-ms AM lifetime events could not be observed). This way we excluded the measurements from contaminant fast myosins. For S1s-cgmRLC and S1sMLC2B, stiffness was measured at 10 M ATP concentration, as these motors showed ϳ3-fold shorter t on as compared with the WT-S1s and S1s-MLC2v, and ϳ2-fold longer average t on compared with the S1f chimeras. Altogether, the mechanics measurement suggests that the light chains may cause substantial alteration of the amplitude of powerstroke size and the motor stiffness as seen for S1s-MLC2B. For the other chimeras, the parameters were either moderately altered or comparable with the WT forms (Table 2).

Discussion
For over 2 decades, the issue of myosin heterogeneity in the muscle fiber has been discussed in the literature, and particular emphasis has been put on the mixed content of myosin isoforms. However, until now, it was difficult to systematically examine the mechanical performance of myosins associating with different types of RLCs. Here, by reconstituting a homogeneous population of such well-defined hybrid/chimeric motors and analyzing the kinetic and mechanical features of the constructs using single-molecule and ensemble measurements, we were able to unravel the influence of individual RLCs on the mechanical performance of slow-and fast-muscle myosins.
By swapping the RLCs between fast-and slow-myosin heavychain isoforms, we showed that MLC2 naturally associating with fast myosin motor and smooth muscle myosin enhanced the originally "slow" motor's actin filament velocity, whereas the MLC2 that naturally assembles with the slow soleus or ␤-cardiac myosin slowed down the native fast myosin. It appeared that these specific MLC2 variants can increase or decrease the velocity, depending on the myosin heavy-chain isoform they are attached to. The changes in different kinetic and mechanical parameters of chimeric motors relative to the WT motor forms are listed in Table 1 and 2. In all but two chimeric motors (i.e. S1s-MLC2B and S1f-MLC2v), the t on that represents the duration of strong bound states correlated with the increase or decrease in the gliding velocity (Table 1 and 2). Although the trapping measurements were done at subsaturating ATP concentration, the relative changes in the t on indicated changes in AM detachment rates among WT and chimeric motors. Compared with its native forms, for chimeric motors, shorter and longer t on relating to higher and lower rates of actomyosin dissociation (1/t on ), consequently resulted in faster or slower actin filament speed, respectively. The powerstroke size was comparable between WT-S1f and S1f chimeras and for WT-S1s and S1s chimeras, except for S1s-MLCB. S1s-MLCB displayed t on similar to that of S1s-CgmRLC, which displayed Ͼ5-fold higher gliding speed. However, the other differences between the two chimeras included a 3-fold shorter powerstroke size and slower rate of weak-strong transition (due to slow P i release; Fig. 3F and Table 1) that may be responsible for the slower gliding speed for S1s-MLCB. Light chains are known to provide the rigidity to the myosin lever arm. It appears that this RLC isoform (MLC2B), when combined with the slow myosin heavy chain (S1s), influences the measured stiffness of the actomyosin cross-bridge. The decrease in the displacement (␦) of S1s-MLC2B is intriguing. One possibility could be that the association of the MLC2B hampers the complete swing of the lever arm to produce the displacement comparable with that of the native motors.
Typically, ADP release (strong bound state) rates are known to determine the gliding velocity for slow myosin. The solution kinetics measurement showed two phases of ADP release with marginal differences in the slow phase of ADP release rates for all slow myosin chimeras (Table 1). However, significant difference (up to ϳ3-fold) in the t on was observed in trapping kinetics measurements under subsaturating ATP concentration, and the gliding speeds varied by up to ϳ10-fold among WT and chimeric motors (S1s versus S1s-cgmRLC). A main difference between solution kinetics and the single-molecule trapping experiment is the load experienced by the motor. Thus, the observed differences in the parameters derived from the two methods could be explained by high sensitivity of the motors to the load as in trapping assays (even at the low applied values) and, consequently, load-sensitive ADP release influencing the mechanical output of the motors. In the solution kinetics measurement, the individual myosin molecules remain unaffected by the binding or release of other myosin heads. However, in the motility assay, several myosin molecules working as an ensemble to drive the actin filament movement are not completely independent but are mechanically linked through the actin filament (i.e. the motors may experience low levels of load).
For the chimera, S1f-MLC2v, the reduction in actin filament gliding speed could be associated with the increase in binding affinity to actin, a higher affinity for ADP, and slower ADP release rate, as shown previously (31). These changes may influence the actin-bound fraction of motors at any given time in the ensemble molecule gliding assay. If the bound motors are not coordinated, due to slower ADP release rates, the increased number of attached versus detached motors may generate a drag on actin filaments as compared with the WT counterpart. This may reduce the gliding velocity in the ensemble molecule assay. Consistent with increased duration of the actomyosin strong-binding state or a higher duty ratio for the S1f-MLC2v motor, we found that the dimeric chimera of fast heavy meromyosin, HMMf-MLC2v, even displayed moderately processive behavior along actin filaments in single-molecule studies (31).
We identified a previously unrecognized aspect of the RLC's task (i.e. the effect on AM cross-bridge cycling by influencing the ATPase kinetics). An obvious question raised by our findings is how diverse RLCs attribute such a response on myosin function. This is a very intriguing yet quite complex question, which would require structural investigations and perhaps additional computational approaches to dissect possible longrange allosteric communication between the RLC and motor domain.
The different isoforms of RLC are similar in structure, as predicted from their aa sequences, in particular in the divalent cation-binding sites and proximal phosphorylatable serines in striated muscles (10). They have relatively low sequence similarity, with the light chains MLC2B/MLC2f and MLC2v 86%

Regulation of myosin II via light chain
similar and the cgmRLC exhibiting only 67% sequence similarity with the fast MLC2B/MLC2f isoforms.
RLC is categorized as an EF-hand superfamily member, harboring helix-loop-helix structure and a cation (Ca 2ϩ or Mg 2ϩ ) binding site at the N terminus. Fig. S2 depicts the conserved eight predicted helices, which form four EF-hand domains. The binding of Ca 2ϩ to the RLC is reported to induce a conformational change to a more open state whereby the overall helical content of the molecule increases, which subsequently positively affects the cross-bridge kinetics (51). In our measurements, we have no Ca 2ϩ but have Mg 2ϩ in the reaction buffers. Additionally, the aa sequences for the helix 1 and 2 regions are 100% identical in fast and slow RLC, whereas the loop that serves as the binding site for divalent cation has a 3-aa difference. With the current knowledge, however, it is challenging to foresee which specific feature(s) might define the particular modulating characteristics of RLCs. Detailed structural and functional investigations would be required to address the diverse impact of the regulatory light chains.
Our findings elucidate the specificity of RLCs toward distinct MHCs. In addition to understanding the fundamental role, precise details of the modulatory role of the RLC isoform are of particular clinical relevance, as single point mutations in the human myosin heavy (␤-MHC) and regulatory light chains (MLC2v) have been linked to causing familial hypertrophic or dilated cardiomyopathy. Researchers have gained a plethora of valuable information not only from studies on human muscle tissue, but also from various animal models (e.g. mice and rabbits) and recombinant motor proteins about the motor dysfunction as a result of single point mutations. These studies often constituted motor components either from mixed isoforms (e.g. human mutant MLC2 with mice MHC and MLC1) or the mutations generated in different MHC isoform background (e.g. mouse ␣-MHC or Dictyostelium discoideum myosin II). Some studies were performed without light-chain binding domains or motors with only MLC1 and not MLC2.
Recent micro-mechanical characterization of human induced-pluripotent stem cell-derived cardiomyocyte myofibrils showed different (increased rates) kinetics of active force generation and relaxation from that of the adult human ventricular myofibrils. One possible explanation was that besides other protein isoform differences, the ␤-MHC expression did not fully match with the respective light-chain expressions (i.e. MLC1v and MLC2v); rather, atrial (MLC1a and MLC2a) light chains dominated over ventricle-specific light chains (52). Altogether, a striking revelation from our studies is the point that the complete "holoenzyme" with the respective heavy and light chains is critical to understanding the motor function and therefore should be taken into account in structure-function analyses.
Overall, our findings offer highly compelling evidence of the role of RLCs in modulating myosin motor function by influencing the actomyosin ATPase kinetics. The functional differences among the myosin II family members appear to be only partly defined by the myosin heavy chain alone. Here, we show that variations in light-chain composition crucially determine the mechanical output of a specific myosin II isoform that is commonly found in different muscle types and various tissues. The importance of these findings extends beyond basic understanding of the myosin function. The study opens new perspectives to investigate the physiological relevance of light chainmediated fine-tuning of myosin function in the context of muscle contraction and additionally for other classes of myosin involved in transport processes.

Generation of single-headed myosin motors from native myosin II
Full-length rabbit fast skeletal muscle myosin II (MHC-IId) and MHC-I were isolated from skinned fibers of musculus psoas and musculus soleus muscle, as described previously (53,54). The myosin II was digested with papain to generate singleheaded subfragment 1 (myosin S1), respectively, as described earlier (55,56). S1 was stored at Ϫ80°C in 50% glycerol in a buffer (5 mM sodium phosphate, 10 mM potassium acetate, 4 mM magnesium acetate, and 2 mM DTT at pH 7.0). In the reconstitution experiments, freshly prepared S1 was employed for further downstream processing for light-chain exchange.
To collect the muscle tissues, the rabbits were euthanized as per the guidelines from German animal protection act §4 (Tötung zu wissenschaftlichen Zwecken/sacrifice for scientific purposes). The animal was registered under reference number 2016/122, obtained from Medical School Hannover central animal facility. All of the procedures were carried out in accordance with relevant guidelines and regulations from Hannover Medical School. No experiments were performed on live animals prior to the sacrifice.

Regulatory light-chain expression
The plasmid vector (EX-T0572-B09) containing full-length human cardiac, slow light-chain 2 insert (i.e. MYL2v) and vector (EX-D0356-B09) containing human, fast-skeletal muscle myosin regulatory light-chain (MYLPF) sequence were purchased from Genecoepia (Rockville, MD).The plasmid vector containing biotin-dependent transcarboxylase (BDTC) cgm-RLC was kindly provided by Dr. Atsuko Iwane to Walter Steffen. The details of BDTC-cgmRLC fusion protein are provided by Iwane et al. (20). As mentioned by Iwane et al., a recombinant fusion protein, BDTC-cgmRLC was generated by fusing BDTC with chicken gizzard muscle myosin RLC sequence at the N terminus. BDTC is a 123-amino acid peptide sequence of the 1.3S subunit of Propionibacterium shermanii transcarboxylase. BDTC was post-translationally biotinylated by biotin holoenzyme ligase in Escherichia coli in vivo. Biotin was added to the bacterial culture in growth medium after the induction of the bacterial cells with isopropyl ␤-D-1-thiogalactopyranoside. The molecular weight of the BDTC-cgmRLC protein is ϳ34 kDa.
The plasmids (EX-T0572-B09 and EX-D0356-B09) were transformed in Rosetta TM competent cells (EMD Millipore, Burlington, MA) for expression of the light chains, MLC2v and MLC2B, respectively. The isolation and purification of the protein were as described previously (20,56). Note that bacterially expressed regulatory light chains, BDTC-cgmRLC, MLC2B (or MLC2F), and MLC2v, were without post-translational modifications, such as phosphorylation.

Regulation of myosin II via light chain
Regulatory light-chain exchange RLC exchange with S1 was performed as described previously (20,31,(56)(57)(58), with some modifications. Myosins, S1f (psoas fast myosin subfragment 1) or S1s (soleus slow myosin subfragment 1), were incubated with an ϳ10-fold molar excess of MLC2v/MLC2B/cgmRLC for 30 min at 30°C in a buffer (Mg 2ϩ -free condition) containing 50 mM HEPES (pH 7.6), 0.5 M NaCl, 10 mM EDTA, 10 mM DTT. The exchange reaction was stopped by the addition of 12 mM Mg 2ϩ to the reaction mixture and incubated on ice for 30 min. To remove the excess RLCs from the reconstituted myosin and to obtain the active motors, myosin was further incubated with high concentration of actin filaments for 30 min on ice. The mixture was centrifuged over 20% sucrose cushion buffer at 70,000 rpm for 30 min. The pellet was once washed with buffer (5 mM sodium phosphate, 10 mM sodium acetate, and 4 mM magnesium acetate, 2 mM DTT, pH 7.0) and then resuspended in the buffer with added 10 mM ATP to release active myosin from actin filaments. The actin filaments were separated from myosin solution by centrifugation. The supernatant containing myosin motors was mixed with 50% glycerol, flash-frozen in liquid nitrogen, and stored at Ϫ80°C and used for further experiments as required. The probes were run on the 15% polyacrylamide gels and stained with Coomassie stain (Phastgel Blue R-350 from GE Healthcare). The gel images were acquired, and either ImageQuant or the ImageJ program was used for the densitometric analysis to determine the RLC exchange efficiency for the chimeric myosins. To measure the intensities of the RLC bands, the region of interest (including both the native and exchanged RLC bands) was selected and plotted as an intensity histogram. The area under the curve for each respective band was measured for the intensity values. The difference in the molecular weight of the protein was taken into account while calculating the exchange efficiency. About Ն85% RLC exchange was observed for several preparations of S1 as shown in Fig. S1. The myosin was used for motility and single-molecule assays. For the solution kinetics experiments, the myosin chimeras were prepared without the final actin co-sedimentation step to avoid ATP in the myosin solution. Instead, after the RLC exchange reaction was completed, the myosin was loaded on spin columns with high-molecular mass cut-off (100 kDa) to separate the free RLCs from the assembled myosin complexes. The removal of free RLC was ensured by multiple dilutions with buffer and passage through the spin columns. The myosin in buffer (5 mM sodium phosphate, 10 mM sodium acetate, 4 mM magnesium acetate, 2 mM DTT, pH 7.0) containing 3% sucrose was finally flash-frozen and stored at Ϫ80°C.

Preparation of actin filaments
Actin filaments were generated by incubating rabbit G-actin in polymerization buffer containing 5 mM sodium phosphate, 10 mM sodium acetate, and 4 mM magnesium acetate, supplemented with protease inhibitor overnight at 4°C. An equimolar concentration of fluorescent phalloidin was added to fluorescently mark the actin filaments. To obtain long biotinylated actin filaments (Ն20 m) for optical trapping experiments, additionally, 1 mM DTT and 1 mM ATP was added to the polymerization mixture as described by Steffen et al. (36).

Active myosin heads
Prior to use in in vitro motility assays or optical trapping measurements, the myosin was further purified to discard any inactive motors by co-sedimentation with concentrated F-actin. Actin-myosin complex was dissociated by the addition of 4 mM Mg⅐ATP. Whereas the active heads released from actin in the presence of ATP, the inactive myosin remained bound and separated by centrifugation at 70,000 rpm for 30 min at 4°C. The supernatant containing the enzymatically active motors was further supplemented with 2 mM DTT and protease inhibitor and used in the functional assays. This procedure to remove the inactive motor heads was routinely followed prior to the main experiments.

Steady-state ATPase and transient kinetic experiments
Actin-activated ATPase was measured with a BioTek Synergy 4 multiplate reader (BioTek) using the NADH coupled assay. A solution of myosin (0.2 M), F-actin, NADH (0.4 mM), lactate dehydrogenase (0.02 mg/ml), pyruvate kinase (0.05 mg/ml), and phosphoenolpyruvate (0.5 mM) in buffer containing 25 mM HEPES, pH 7.3, 25 mM KCl, 5 mM MgCl 2 , 0.5 mM DTT was mixed with ATP (2 mM) to start the reaction. ATPase rates at different concentrations of F-actin were obtained from the slopes of the corresponding linear fits of the time-dependent absorbance change at 340 nm. Rates were plotted against the F-actin concentration, and data were fitted according to Michaelis-Menten kinetics. Transient kinetic experiments were carried out with a Hi-Tech Scientific SF-61 DX double mixing stopped-flow system (0.5-ms dead time) at 20°C. Experimental buffer contained 20 mM MOPS, pH 7.0, 25 mM KCl, 1 mM DTT, 5 mM MgCl 2 . For ADP-release measurements, a solution of 2 M AM preincubated with 20 M mantADP was rapidly mixed with 1.6 mM ADP. The reaction was followed through a KV389 cutoff filter by monitoring mantADP fluorescence changes induced by FRET through tryptophan excitation at 296 nm. P i release was measured under single-or multipleturnover conditions with ATP as substrate as described before (59).

In vitro motility assay
An in vitro motility assay was performed with monomeric S1 motors with native or exchanged RLC, by immobilization of the motors on a nitrocellulose-coated surface. The assay was described previously in more detail (53). Briefly, myosin was incubated for 5 min on a nitrocellulose-coated surface, followed by wash and surface blocking with 1 mg/ml BSA in assay buffer (25 mM imidazole hydrochloride, pH 7.2, 25 mM NaCl, 4 mM MgCl 2 , 1 mM EGTA, and 2 mM DTT). To block the inactive or damaged myosin motor heads, 0.25 M short, unlabeled F-actin was injected in the flow cell and incubated for 1 min. 2 mM ATP was introduced in the chamber to release the actin filaments and to make the active motor heads accessible. ATP was washed out with assay buffer, tetramethylrhodamine-labeled F-actin was incubated for 1 min and washed to remove excess filaments, and finally the chamber was infused with buffer con-

Regulation of myosin II via light chain
taining 2 mM MgATP and antibleach system (18 g/ml catalase, 0.1 mg/ml glucose oxidase, 10 mg/ml D-glucose, and 10 mM DTT) to initiate F-actin motility. Images were acquired with a time resolution of 200 ms (i.e. 5 frames/s) using a custom-made objective-type total internal reflection fluorescence microscope. Actin filament gliding speed was analyzed with the Manual Tracking plugin from ImageJ.

Three-bead assay with optical tweezers
For the assay, flow cells with ϳ15-l chamber volumes were assembled using a coverslip with nitrocellulose-coated beads. Glass microspheres (1-1.5 m) suspended in 0.1% nitrocellulose in amyl acetate were applied to 18 ϫ 18-mm coverslips. All of the dilutions of biotin-actin filaments and myosin S1 were made in reaction buffer containing 25 mM KCl, 25 mM Hepes (pH 7.6), 4 mM MgCl 2 , and 1% ␤-mercaptoethanol. For the experiment, the chamber was prepared as follows. 1) Flow cells were first incubated with 1 g/ml myosin S1 for 1 min. 2) This was followed by a wash with 1 mg/ml BSA and further incubation for 1 min to block the surface. 3) Finally, the reaction mixture, containing 0.8 -1-m neutravidin-coated polystyrene beads and 1-2 nM biotinylated actin filaments (60), was flowed in with 10 M ATP, ATP-regenerating system (0.01 mM creatine phosphate and 0.01 unit creating kinase), and deoxygenating system (1 mg/ml catalase, 1 mg/ml glucose oxidase, 2 mg/ml glucose, and 1% ␤-mercaptoethanol). As shown in the schematic view in Fig. 4A, prestretched, biotinylated actin filament was held between the two optically trapped beads via a neutravidin-biotin link forming a dumbbell, such that the actin filament formed a connection between two beads. The approximate distance between the two trapped beads was 10 m. Low-compliance links between the trapped beads and the filament were adjusted such that the ratio of the position variance during free and bound periods was 5-10, as described by Smith et al. (49). All of the experiments were performed at room temperature of ϳ22°C. WT and chimeric motor S1 were immobilized directly on a nitrocellulose-coated surface to have comparable conditions for the single molecules. The actin dumbbell was brought in contact with the myosin bead, and the AM interaction events were registered as a reduction in free Brownian noise of the two beads. The two bead positions were precisely detected with a quadrant photodetector, recorded, and analyzed.
Data traces were collected at a sampling rate of 10,000 Hz and filtered at 5000 Hz. To improve the time resolution and detect short-lived AM binding events, a high-frequency triangular wave of ϳ600 Hz was applied to one of the traps as described by Veigel et al. (61). Matlab routines were employed to evaluate data records for AM interaction lifetime, t on , and stroke size of motors.

Single-myosin molecule interaction with actin filaments in optical trapping measurements
To ensure that each data record was derived from an intermittent interaction between single myosin head and actin, myosin density on the bead surface was adjusted by diluting the myosin solution. In our measurements, typically, one bead was found to interact with dumbbell among 8 -10 beads scanned for the presence of motor on the bead. This measure minimized the likelihood of multiple molecules simultaneously interacting with the actin filaments. From the binomial probability distribution, in our measurements, chances of the presence of two motors per bead are Ͻ1%, and chances of one molecule per bead are ϳ9%. From a total of 126 beads we analyzed for AMbinding events, statistically, 0.4 beads may have two motors. This estimation, however, does not take into account the motors that are inaccessible to actin filament due to the positional limitation on the bead. Therefore, the probability of multiple motors interacting simultaneously with an actin filament is even lower.

Statistical analysis
Binomial probability distribution analysis was used to estimate the chances of more than one molecule interacting with actin filaments in optical trap measurements. Unpaired t test was used to calculate the statistical differences in the gliding velocities, single-motor stroke size, and stiffness of the native versus chimeric motors. The nonparametric Mann-Whitney U test was used to calculate the statistical differences in the duration of AM interaction events, t on , for the WT and chimeric motors. The statistical tests used to calculate statistical significance are included throughout in relevant sections.
All of the experimental protocols were performed in accordance with guidelines for good scientific practices and approved by Hannover Medical School.

Data availability
All data are available within the article.