## Force and kinetics of fast and slow muscle myosin determined with a synthetic sarcomere-like nanomachine

Valentina Buonfiglio 1 , Irene Pertici 1 , Matteo Marcello 1 , Ilaria Morotti 1 , Marco Caremani 1 , Massimo Reconditi 1 , Marco Linari 1 , Duccio Fanelli 2* , Vincenzo Lombardi 1* , Pasquale Bianco 1

1

PhysioLab, University of Florence, Sesto Fiorentino (FI), Italy 2 Department of Physics and Astronomy, University of Florence, Sesto Fiorentino (FI), Italy

## Abstract

Myosin II is the muscle molecular motor that works in two bipolar arrays in each thick filament of the striated (skeletal and cardiac) muscle, converting the chemical energy into steady force and shortening by cyclic ATP-driven interactions with the nearby actin filaments. Different isoforms of the myosin motor in the skeletal muscles account for the different functional requirements of the slow muscles (primarily responsible for the posture) and fast muscles (responsible for voluntary movements). To clarify the molecular basis of the differences, here the isoform-dependent mechanokinetic parameters underpinning the force of slow and fast muscles are defined with a unidimensional synthetic nanomachine powered by pure myosin isoforms from either slow or fast rabbit skeletal muscle. Data fitting with a stochastic model provides a self-consistent estimate of all the mechanokinetic properties of the motor ensemble including the motor force, the fraction of actin-attached motors and the rate of transition through the attachment-detachment cycle. The achievements in this paper set the stage for any future study on the emergent mechanokinetic properties of an ensemble of myosin molecules either engineered or purified from mutant animal models or human biopsies.

Keywords : Myosin motor; muscle myosin isoforms; synthetic myosin-based nanomachine

## Introduction

Steady force and shortening in the half-sarcomere, the functional unit of the muscle cell, are due to ATP-driven cyclic interactions of the subfragment 1 (S1, the head) of the heavy meromyosin fragment (HMM) of the myosin II molecule extending from the thick filament with the actin monomers on the nearby overlapping thin filament (Figure 1). In each interaction the free energy of the splitting of MgATP to MgADP and inorganic phosphate (Pi) in the head is associated to a structural working stroke consisting in a tilting between the motor domain firmly attached to actin (red) and the light chain binding domain (or lever arm, violet) connected to the myosin filament backbone through the subfragment 2 (S2, the tail, green). In isometric contraction, lever arm tilt raises the force exerted by the half-sarcomere, increasing the strain of all the elastic elements represented for simplicity in Figure 1, state (b) by the bending of the lever arm.

When the load is lower than the maximum steady force exerted under isometric conditions ( T 0 , that is conventionally expressed as force per cross-sectional area of the contractile material, kN m -2 ), lever arm tilting results in relative filament sliding with a reduced strain in the elastic components

∗ Corresponding authors E-mail addresses: vincenzo.lombardi@unifi.it, duccio.fanelli@unifi.it

Schematic diagram of the chemo-mechanical cycle of the myosin motor during its interaction with the actin filament.

Fig. 1. The HMM fragment of the myosin molecule is a dimer with each monomer made by the subfragment 1 ( S1 or head containing the motor domain (red) and the light chain domain (the lever arm, violet)) and the subfragment 2 ( S2 or tail, green) extending from the myosin filament backbone (blue). For simplicity, only one S1 and S2 are represented here. The myosin ADP Pi · · complex attaches to actin (white circles) (a) , forming the cross-bridge, which triggers tilting of the lever arm and Pi release with generation of force and actin filament sliding. If the mechanical load is high, it opposes filament sliding and tilting of the lever arm causes increase in strain in the system, represented here by the distortion of the lever arm (b) . If the load is low (c) tilting of the lever arm causes actin filaments sliding (yellow arrow), keeping the strain low. ADP release from and ATP binding to the motor domain cause myosin detachment from actin. ADP release is slower at high load, ( b ) → ( d ) , and becomes faster at lower load ( c ) → ( d ) . Hydrolysis of ATP in the detached head and reversal of the lever arm tilting (recovery stroke, ( d ) → ( e ) ) completes the cross-bridge cycle. The absence of ATP causes the cycle to stop before detachment so that all motors stay attached to actin (rigor).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-a6dd7635-e7ba-4cb2-90f8-1790b69a34d4.jpg)

(Figure 1, state (c)). Cyclic asynchronous interactions of myosin motors with the actin filament account for the generation of steady force and shortening. The shortening velocity V is inversely proportional to the force T (force-velocity relation, T -V 1 ). At physiological concentrations of ATP, ADP release is the rate-limiting step for motor detachment from actin (step ( ) b / c ( ) → ( d )). The rate of ADP release is conformation-dependent, increasing during steady shortening when motors at the end of the working stroke would become negatively strained. This explains the increased rate of energy liberation ˙ E (and the underlying ATP hydrolysis rate, ϕ ) when the load is reduced and the shortening speed is increased 1-8 . Faster detachment of negatively strained (compressed) motors prevents the ones at the end of their working stroke to oppose positively strained motors, a requirement for the maximisation of the power and efficiency of an array of motors working in parallel. The performance of different types of skeletal muscles depends on the myosin II isoform expressed in the muscle. Slow muscles, which are primarily involved in maintenance of posture and are characterised by the dominant presence of the isoform 1 of Myosin Heavy Chain (MHC -1 isoform), exhibit lower shortening speed at any given load, thus develop lower power and consume ATP at a lower rate than fast muscles which are involved in movement and are characterised by the dominant presence of the isoforms MHC -2 A , -2 B or -2 X isoforms . 9 Strikingly, the functional

difference between slow and fast isoforms is due to a difference of only 20% in the amino-acid composition.

During an isometric contraction the power is zero, so that the rate of energy consumption accounting for the steady force T 0 (denoted ˙ E 0 ) corresponds to the rate of heat production ( ˙ H 0 ) 1 . ˙ E 0 has been found ∼ 5-fold larger in fast muscles than in slow muscles 10-13 . The underlying rate of ATP hydrolysis at T 0 can be obtained from ˙ E 0 by dividing it by the energy liberated per molecule of ATP hydrolysed (∆ G ATP = 60kJmol -1 in mammalian muscle according to 14 ). In this way the energetic cost of the isometric force in the intact muscle can be compared with that in the demembranated fibres, in which the rate of energy liberation is determined by measuring the rate of ATP hydrolysis. A further normalisation for the concentration of myosin motors in the mammalian muscle (0 18 mM, . 11 ) gives the rate of ATP hydrolysed per myosin motor ( ϕ ). In demembranated fibres of rat fast muscle 15,16 , rabbit 17,18 and human muscle 19 at 12 ◦ C, ϕ is 5-fold (or more) larger than in slow muscle, in agreement with muscle measurements (Supplementary Table 1). In both fast muscles 10-13 and fast demembranated muscle fibres 15-18 T 0 is either similar or at max 1 5-fold . larger than in slow muscles and muscle fibres. Thus the tension cost of the isometric contraction ˙ E /T 0 0 results to be systematically larger in fast muscles by on average 5-fold (with a minimum of 3-fold). The justification for the elevated tension cost of the fast muscle can be only partly found in the intrinsic larger actin-activated myosin ATPase in solution, which for the fast myosin is twice that of the slow myosin 20 .

The bulk of data characterising the energetics of slow and fast muscles at cell and tissue levels, first of all the ∼ 5-fold larger isometric tension cost, leaves open the question of the underlying molecular mechanism. Inferring the definition of the molecular mechanism from cell and tissue is complicated by the scaling factors related to the structural organisation of the molecular motors in the three-dimensional lattice, the co-presence of different isoforms in the same muscle and even in the same muscle fibre and the possible confounding contribution of the other sarcomeric (cytoskeleton and regulatory) proteins. Even assuming that the tension cost is solely related to intrinsic properties of the motor isoform, the question remains about the role played by the differences in the mechanokinetic properties of the motor, as the force developed in a single motor interaction or the fraction of the ATP-ase cycle time each motor spends attached (the duty ratio) while working in situ in the half-sarcomere of the striated muscle. In vitro , the definition of the emergent properties of the half-sarcomere, which cannot be studied with single molecule mechanics 21,22 , became recently accessible by exploiting a unidimensional synthetic nanomachine powered by myosin motors purified from the skeletal muscle 23 . The nanomachine allows the performance of the half-sarcomere, the generation of steady force and shortening, to be reproduced by an ensemble of pure myosin isoforms interacting with the actin filament without the confounding effects of other sarcomeric proteins and higher hierarchical levels of organisation of the muscle. In the nanomachine 8 HMM fragments extending from the functionalised surface of a micropipette carried on a three-way nanopositioner acting as a length transducer interact with an actin filament attached with the correct polarity to a bead trapped by a Dual Laser Optical Tweezers (DLOT) acting as a force transducer (Figure 2). In solution with physiological ATP concentration, in which the two motors of each dimer act independently 23 , myosin motors, after entering in contact with the actin filament, establish continuous interactions underpinning force development to a steady maximum value (equivalent to the force generated by the muscle in isometric contraction). In the original design 23 the system was operated either in position clamp (achieved using as feedback signal the position of the nanopositioner carrying the motor array ( x ), red branch in Figure 2), to reproduce the isometric contraction, or in force clamp (achieved using as feedback signal the position of the bead in the laser trap ( x bead ), green branch in Figure 2), to reproduce isotonic contraction.

A major limit of the nanomachine working in position clamp was the large trap compliance in series with the motor array, two order of magnitude larger than the native compliance in series with the half-sarcomere. This makes each addiction-subtraction of force by individual motor attachment-detachment to induce substantial sliding undermining the condition of independent force generators of the motors in the native half-sarcomere. Consequently, in position clamp, the kinetics of the attached motors is influenced by the push-pull experienced when actin slides away-toward the bead for the addition-subtraction of the force contribution by a single motor ( 23 , Supplementary Figure 7).

Block diagram of the system for nanomachine mechanics.

Fig. 2. HMMfragments (blue) deposited on the functionalised lateral surface of a pulled micropipette (cyan) are brought to interact with the actin filament (red) attached with the correct polarity (+) via gelsolin (yellow) to the bead trapped in the focus of the DLOT (pink). Force generation produces the movement of the bead away from the focus of the DLOT. The switch selects the feedback signal that, together with the command (black), feeds the summing amplifier Σ that drives the piezo nanopositioner: in position clamp (red) the feedback signal is the position of the nanopositioner x carrying the support for the myosin array; in force clamp (green) the feedback signal is the force ( F , calculated as the product of the stiffness of the trap ( e ) times the change in position of the bead x bead ); in length clamp (blue) the feedback signal is the change in the distance L between the position of the bead and the myosin array support.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-f1bf0b7f-cd79-4493-83f4-a2347b5d0b70.jpg)

In the present experiments the system has been implemented to operate in length clamp (achieved using a feedback signal the difference between the position of the bead and that of the nanopositioner x bead -x = L ), blue branch in Figure 2. In length clamp the sliding between the actin filament and the motor array caused by force generating interactions is eliminated because any movement of the bead is counteracted by the movement of the nanopositioner. In this way the condition of the motors as independent force generators in the array is recovered and the rate of development of the steady isometric force as well as the force fluctuations superimposed on the steady force are direct expression of attachment/force-generation and detachment of the myosin motors. The data collected from either myosin isoform are used to feed a stochastic model providing a self-consistent estimate of all the relevant mechanokinetic parameters of the isometric performance of the motor ensemble: the force of a motor, f 0 , the fraction of actin-attached motors, r , and the rate of transition through the attachment-detachment cycle, ϕ , without assumptions from cell mechanics and solution kinetics as in previous studies 23-25 . The combined experimental and theoretical achievements reported in this paper set the stage for any future study on the emergent mechanokinetic properties of an ensemble of myosin molecules, either engineered or purified from mutant animal models or human biopsies.

Estimating the number of HMMs available for actin interaction from rigor rupture events in ATP-free solution.

Fig. 3. a . 1 . Formation of the rigor bonds between the HMM array and the actin filament. 2 . The motor support is moved away first in the direction ( z ) perpendicular to the plane of the actin-myosin interface and then in the direction ( x ) parallel to the plane, as indicated by the arrow. Panel modified from 23 . b . Force ( F x , lower record) of an ensemble of soleus HMMs in response to the movement of the nanopositioner away from the actin filament in the x direction (upper record, velocity 50 nms -1 ). The small vertical bars indicate the rupture events (force drop complete in less than 50 ms ), the last of which corresponds to complete detachment of the actin filament. c . Records with the same protocol applied to an ensemble of psoas HMM.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-e1ce516e-a1aa-4624-93d8-2fccc0ae3fcb.jpg)

## Results

## Estimate of the number of HMM molecules on the support available for the interaction with the actin filament.

The number of motors on the micropipette surface able to interact with the actin filament ( N ) is initially determined by measuring the number of mechanical rupture events when the motor array is brought to interact with the actin filament in ATP-free solution 23 . Following the formation of rigor bonds between the HMM-coated support and the actin filament (panel 1 in Figure 3 a), the HMM support is moved away from the actin filament, first by 1-2 µ m in the direction orthogonal to the plane of the support, in order to raise a force from the trapped bead to the first bound HMM at an angle greater than 30 ◦ with the plane of the support, and then in the direction parallel to the plane, at constant velocity (50 nm s -1 ), to pull the motors away from the actin filament diagonally. This allows the first bonded HMM to undergo a pulling force that is higher than the axial component shared among the other motors. In this way the myosin-actin bonds brake one at a time and the attached motors cannot bind back once detached. Moreover, following each detachment the force drops because the length of actin filament segment between the bead and the next attached motor is transiently increased. Thus an additional pull is necessary to get to the next rupture event, the occurrence of which will vary in time according to the distance between the two neighbouring motors. With HMMs purified from soleus muscle the number of rupture events per interaction (Figure 3b) attains a saturating value of 7 9 . ± 1 1 ( . n = 8), with [HMM] used to coat the pipette of 0 2 mg L . -1 . A similar saturating value of rupture events, 8 1 . ± 1 4 ( . n = 8), is obtained for the HMM purified from psoas muscle with a [HMM] of 0 1 mg L . -1 (Figure 3c). Notably, similar saturating values of [HMM] and number of rupture events (8 1 . ± 1 2) were found for the psoas . motors in the previous study in which an optical fibre etched to the same diameter was used as support 23 . In 2 mM ATP each head of an HMM dimer works independently and thus the number of available motors is twice the number of HMM ruptures: N = 16 ± 2 and 16 ± 3 for the soleus and psoas respectively.

## Isometric force development by the nanomachine powered by slow and fast myosin motors.

The experiment starts in position clamp, because, for the system to operate in length clamp, it is necessary that first the feedback loop is closed by the establishment of actin-myosin interactions. When an array of motors from the slow soleus muscle is brought to interact with a bead-tailed actin filament in solution with 2 m M ATP (Figure 4a), the establishment of continuous ATP-driven actin-myosin interactions causes the force ( F , blue trace) to rise pulling on the actin filament, which in position clamp (HMM support position x = 0, red trace), slides in the shortening direction (∆ L , black trace, negative for shortening) due to the trap compliance (phase 1). A steady maximum force ( F 0 ) of ∼ 12 pN is attained with a shortening of ∼ 55 nm. The control is switched to length clamp in correspondence of the vertical dashed line separating phase 1 and 2. The switch time is marked by the increase in noise of the force trace as a consequence of the reduction of the compliance in series with the motor system. In fact in length clamp the force change generated in each individual attachment and detachment is no longer dissipated in filament sliding against the large in series trap compliance. A shortening of ∼ 500nm completed within ∼ 700ms is superimposed on the steady isometric force in correspondence of the second vertical dashed line to drop and keep the force at zero (phase 3). When actin filament sliding stops (third vertical dashed line) force starts to redevelop towards F 0 (phase 4) with just a minimum delay, indicating that the motor array was able to cope with the imposed 500 nm shortening maintaining continuous interactions under zero load. The extent of shortening minus the amount taken by the trap compliance, (500 -55 =) 445nm, divided by the time passed from the imposition of the shortening to the start of force redevelopment (0 88 s) gives an estimate of the velocity of unloaded shortening ( . V 0 ) of 0 5 . µ ms -1 . Force redevelopment in length clamp is much faster than the original force rise in position clamp and occurs in truly isometric conditions, as the movement of the bead due to trap compliance is counteracted by a corresponding movement of the nanopositioner in the lengthening direction (red trace, ∼ 55nm), that keeps ∆ L = 0 (black trace). Notably, the force redevelopment following a 500 nm release attains the same F 0 value as that attained during the original rise in position clamp thanks to the architecture of the machine, in which the dimension of the motor array remains constant independent of the amount of reciprocal sliding 23 . F 0 from 33 records shows a Gaussian distribution with centre 10 5 pN (Figure 4 b). . The rate of force redevelopment in length clamp only depends on the attachment/detachment kinetics of myosin motors in isometric conditions. Force redevelopment is roughly exponential, and its time course is quantified by the rise time t r (the time from 10% to 90% of F 0 ). t r estimated on the record (Figure 4c, black) obtained by averaging the blue 6 traces from as many experiments (grey) is 238 ± 13 ms. The time constant of the underlying exponential force rise of the soleus powered nanomachine, τ is ( t r / . 2 2 =) 108 ± 5 ms, and the rate of force development, a is (1 /τ =) 9 3 . ± 0 5 s . -1 . The sequence of events accompanying the interaction of the array of motors purified from psoas muscle with the actin filament is the same as for the soleus motors (Figure 4 d). The Force develops in position clamp (phase 1), while the actin filament slides in the shortening direction due to the trap compliance. A steady isometric force F 0 (15 9 pN), . is attained with a shortening of 70 nm. In the 47 records of the psoas HMM F 0 shows a Gaussian distribution with centre 17 pN (Figure 4e). Following the switch to length clamp, a rapid shortening of ∼ 500 nm is imposed so that the force drops to zero. The shortening in this case is just sufficient to drop the force to zero, given the much faster shortening velocity afforded by the fast motor array, so that, as soon as the actin filament sliding stops (third vertical dashed line), V 0 , calculated by the extent of shortening minus the amount taken by the trap compliance, (500 -70 =) 430 nm, divided by the time passed from the imposition of shortening to the start of force redevelopment (0 22 s), . is 1 95 . µ ms -1 (3 9 times larger than that of slow muscle). . It must be considered, however, that V 0 in this case is somewhat underestimated, as most of the shortening occurs with force greater than zero. Force redevelopment in length clamp (phase 4) occurs with a rate that is not influenced by the trap compliance and thus is the expression of the kinetics responsible for the transition to the steady force F 0 by the fast isoform array. A t r of 77 ± 4 ms is estimated on the record (black in Figure 4f) obtained by averaging the traces from 7 experiments (grey). The corresponding τ (= t r / 2 2) and . a (= 1 /τ ) are 35 0 . ± 1 8 ms and 28 6 . . ± 1 4 s . -1 respectively.

The -3 dB upper frequency characterising the force rise f c = 0 35 . /t r is 4 5 . ± 0 2 Hz. .

Two main points emerge from these measurements on the synthetic machine operating in length

Active force generation by the nanomachine powered by slow (soleus) and fast (psoas) myosin motors.

Fig. 4. a -c . Slow myosin array. a . Force ( F , blue trace), movement of the nanopositioner carrying the motor array ( ∆ x , red trace) and relative sliding between the motor array and the actin filament ( ∆ L , black trace) during the actin myosin interaction. Phase 1 : following the establishment of the contact between the actin filament and myosin motors, the force rises in position clamp to the maximum isometric value ( F 0 ≃ 12 pN ), with the simultaneous sliding of the actin filament by ∼ 55 nm toward the shortening direction to load the trap compliance. Phase 2 : the switch to length clamp (marked by the first vertical line) is followed by the increase in force fluctuations superimposed on F 0 . Phase 3 : force drops to zero in response to a rapid shortening of ∼ 500 nm imposed in length clamp (start marked by the second vertical line) with actin filament sliding under zero force. Phase 4 : following the end of the imposed shortening (marked by the third vertical line) force redevelops in length clamp with the nanopositioner moving by ∼ 55 nm to counteract the trap compliance and keep the filament sliding at zero. b . Frequency distribution of F 0 . Data are plotted in classes of 1pN and fitted with a Gaussian (continuous line) with centre 10 5 pN . and standard deviation σ = 1 8 pN . . c . Time course of force redevelopment after rapid shortening (black trace) averaged from 6 records from as many experiments (grey traces). The red line is the single exponential fit to measure t r (the time elapsed from 10% , horizontal thin dashed line, to 90% , thick horizontal dashed line, of F 0 recovery). d -f . Fast myosin array. d . F , ∆ x and ∆ L , defined and colour coded as in a . Phases 1 -4 as described in panel a . e . Frequency distribution of F 0 plotted in classes of 2 pN and fitted with a Gaussian (continuous line) with centre 17 pN and standard deviation 3 pN . f . Time course of force redevelopment after rapid shortening (black trace) averaged from 7 records from as many experiments (grey traces). The red line is the single exponential fit to measure t r labelled as in c .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-d45b6e1a-04fe-423e-8bfc-7c33629d2205.jpg)

clamp conditions. The first is that the rate of force redevelopment, which only depends on the

attachment/detachment kinetics of myosin motors in isometric conditions, is three times slower in the soleus powered nanomachine than in the psoas powered nanomachine. The second point is that the force fluctuations around the average value displayed by the force record at the steady state are stemming from individual attachment/detachment events. Both pieces of information can be used to feed the stochastic model described in the next section.

## Modelling the mechanical output of the nanomachine powered by the slow and the fast isoform ensembles.

In the stochastic model each motor exists in three possible states (or motor configurations): one detached state and two different force-generating attached states. Fitting the experimental records allows a self-consistent estimate of all the relevant mechanokinetic parameters of the nanomachine including the force exerted by a single myosin motor and the average number of attached motors in the stationary state, without assumptions from cell and solution kinetic studies.

As detailed in the Introduction, the implementation of the length clamp mode allows to recover the condition of the motors as independent force generators in the array. Therefore we consider an ensemble of N independent ATP-fuelled molecular motors interacting with an actin filament in isometric conditions. Each motor can be found in the detached state D , in the attached low force-generating state A 1 , or in the attached high force-generating state A 2 . The corresponding kinetic scheme, which exemplifies the possible transitions between distinct allowed motor states is:

$$D \xrightarrow { k _ { 1 } } A _ { 1 } \xrightarrow { k _ { 2 } } A _ { 2 } \xrightarrow { k _ { 3 } } D$$

The rate constants k j , j ∈ { 1 , -1 2 , , -2 3 , } represent the probability per unit of time for the reaction j to occur, and are expressed in s -1 . The state of the system at any time t is characterised by the vector n ( ) = t ( n D ( ) t , n 1 ( ) t , n 2 ( ) t ) whose entries specify the number of molecular motors in each of the considered configurations. Specifically, n D stands for the number of motors in the state D n , 1 is the number of motors in the state A 1 and n 2 denotes the number of motors in the state A 2 . The system admits the obvious conservation law N = n D + n 1 + n 2 where N stands for the total number of motors in any of the considered states. Accounting for the above relation enables one to employ just two scalar (discrete) entries to photograph the state of the system, namely n ( ) = t ( n 1 ( ) t , n 2 ( ) t ) . Under the Markov hypothesis, the stochastic dynamics of the scrutinised system is ruled by a master equation which sets the evolution of the probability P ( n , t ) of finding the system in the state specified by the vector n at time t . The master equation accounts for the balance of opposing contributions: on the one side the transitions towards the reference state (the associated terms bearing a plus sign). On the other, the transitions from the reference state (terms with a minus). From the master equation one can readily derive the mean field equations that govern the deterministic dynamics for the continuous concentrations y and z of the molecular motors in states A 1 and A 2 , as detailed in Methods. Further, a self-consistent elimination of the variable y can be performed to yield a simpler description of the examined process in terms of the variable z (the derivation is given in Methods, see also Supplementary Figure 1):

$$z ( t ) = \frac { b } { a } \left ( 1 - e ^ { - a t } \right ) = \frac { k _ { 1 } } { k _ { 1 } + G } \frac { k _ { 2 } } { k _ { 2 } + k _ { - 2 } + k _ { 3 } } \left ( 1 - e ^ { - \left [ k _ { - 2 } + k _ { 3 } + \frac { k _ { 2 } \ ( k _ { 1 } - k _ { - 2 } ) } { k _ { 1 } + k _ { - 1 } + k _ { 2 } } \right ] t } \right )$$

where a and b are positive quantities, self-consistently defined by the latter equality and G = ( k -1 ( k -2 + k 3 ) + k k 2 3 ) / k ( 2 + k -2 + k 3 ). We define the duty ratio r as the fraction of attached motors (or the fraction of the ATPase cycle time a motor spends attached). Assuming that for mammalian muscle myosin at temperature T ≃ 24 ◦ C motors are prevalently found in the state A 2 , a straightforward analysis outlined in the Methods yields r ≃ z ∗ = k 1 k 1 + G . Here, z ∗ stands for the equilibrium concentration of A 2 motors. Through parameter a , we have also access to a closed estimate for the characteristic time scale of the exponential evolution of z . Let us notice that a is the inverse of the time constant of the development of the steady force, τ as defined in the experiment, hence a = 2 2 . /t r . According to this simplified scheme, the effective rate of ATP consumption can be estimated as the flux of motors through the cycle per unit time. This equals to the rate of motors in A 2 detaching from the actin, in formula ϕ = z ∗ ( a -b ).

Starting from these premises we can characterise the average force exerted by a small ensemble of molecular motors in isometric conditions. This is obtained by combining the contributions from each individual motor of the collection: motors in the configuration A 1 , each exerting a force f 1 and motors in A 2 , each exerting a force f 2 . In the nanomachine, motors have a random orientation with respect to the actin filament. As a consequence, we assume that the intensity of the force exerted by a motor depends on the binding angle θ , as measured from the correct orientation. Depending on the specific orientation of the molecule, the force progressively decreases up to a minimum value that is 0 1 . f 0 22 . In particular, the force of a single motor can change within a bounded interval: the largest value of the force f 0 is exerted when the motor orientation is correct (corresponding to the in situ orientation). Then, in accordance with 24 (see Supplementary Figure 2) we postulate that the exerted force f 1 is a random variable, uniformly distributed within the interval I 1 = [ -f , f 0 0 ] . Similarly, the force f 2 is randomly extracted from the interval

I 2 = [ f 0 10 , f 0 ] . The mean field average force exerted by the ensemble of myosin motors, at any time t , is ⟨ F t ( ) ⟩ = ⟨ n 1 ( ) t ⟩ ⟨ f 1 ⟩ + ⟨ n 2 ( ) t ⟩ ⟨ f 2 ⟩ = ⟨ n 2 ( ) t ⟩ ⟨ f 2 ⟩ , given that ⟨ f 1 ⟩ = 0 since the interval I 1 is symmetric with respect to zero. In the stationary state, ⟨ F t ( ) ⟩ converges to the asymptotic plateau F 0 , and thus:

$$r f _ { 0 } = \frac { 1 } { N } \frac { 2 0 } { 1 1 } F _ { 0 } \, _ { \mu }$$

where use has been made of the conditions ⟨ f 2 ⟩ = (11 / 20) f 0 and z ∗ ≃ r .

The experimental value of the stationary force exerted by a pool of N motors acting in the state A 2 solely constrains the product of f 0 and r . That is, on deterministic means, we cannot access a direct estimate of the maximum force exerted by an individual motor ( f 0 ) and the associated duty ratio ( r ), but just constraint this latter pair to fall on a hyperbole, set by F 0 . Accounting for the fluctuations superimposed on F 0 , and thus by properly gauging the stochastic component of the dynamics, enables us to resolve the above degeneracy.

To this end we consider the dynamics of the system at finite N , so as to account for the role played by finite size fluctuations. To quantify the contribution as stemming from the intimate graininess of the investigated system, we ought to solve the master equation, focusing in particular on the stationary state probability distribution P st . As discussed in Methods, we are in a position to solve exactly the stochastic model in its stationary state, and thus get a closed expression for P st , as function of the parameters of the model. This knowledge can be used to compute P F ( ) the distribution of the exerted force F (see Methods). Remark that P F ( ) is ultimately shaped by the kinetic constants of the model (namely, k , k 1 -1 , k 2 , k -2 , k 3 ) and also reflects the maximum force f 0 , as applied by individual motors. Recall that N is directly determined (Figure 3). We can therefore construct an inverse procedure to recover information on the underlying parameters, by confronting the predicted distribution of the force P F ( ) to the homologous curve recorded experimentally. In particular we will prove that, by exploiting the information stemming from the fluctuations, it is eventually possible to unambiguously determine both f 0 and r (see Supplementary Table 2, Supplementary Figure 8 and Supplementary Figure 9). The relevant steps that define the envisaged fitting strategy are summarised in the following and made explicit in the Methods:

- 1. The first step amounts to analyse the time evolution of the force in its mean field approximation: the asymptotic force F 0 and the time scale a , as defined above, are extracted via a direct two parameters - fit that exploits expression (2).
- 2. We turn to study the distribution of the fluctuation of the force around the equilibrium value. To this end we make use of P st .
- 3. From P st we extract the N +1 marginal probabilities ρ k , namely the probabilities to find k ≤ N motors in the force-generating configuration A 2 . This is achieved by summing over n 1 = 0 , . . . , N the stationary probability distribution P st .
- 4. We then make use of the marginal probabilities ( ρ , ρ 0 1 , ρ 2 , . . . , ρ N ) to weight the probability distributions Π k ( f ) of the force exerted by a set of k motors. These latter are computed as generalised Irwing-Hall distributions for independent and identically distributed random variables f drawn from the considered interval I 2 . The distribution of the force is hence estimated as P F ( ) = ∑ N k =0 ρ k Π ( k f ).

- 5. For fixed size N (previously estimated by the counting of rupture events in ATP-free solution, see also Supplementary Figure 10 where the possibility to modulate N is accounted for) we adjust the kinetic rate constants k , k 1 -1 , k 2 , k -2 , k 3 , so as to minimise the root mean square distance between the recorded distribution and its analytic estimate. The best fit values are used to compute the parameter b and thus determine the sought estimates for r and f 0 , as well as the rate of motors detaching from the actin, i.e. ϕ = z ∗ ( a -b ).

The above fitting strategy is successfully challenged against synthetic data as discussed in Methods. Then we proceed by applying the validated procedure to the experimental data collected with the nanomachine powered by either HMM isoforms. As mentioned, the number of available molecular motors ( N = 16) estimated from number of ruptures in rigor for both isoforms (Figure 3), is assumed as the reference value in the following, unless otherwise specified. We interpolate the distribution of the fluctuations as recorded experimentally, given the analytical solution obtained above. A representative example of the fitting outcome for the soleus HMM ensemble is reported in the Supplementary Figure 11. The estimated values for f 0 , r and ϕ , for both the psoas and the soleus HMM, are listed in the Supplementary Note 1 (Supplementary Table 3 and Supplementary Table 4 respectively), and their respective average values are reported in Table 1. In Figure 5

Table 1. Average values (mean ± SD obtained by averaging over 6 data records for each isoform) of the three relevant parameters estimated by the stochastic model: the force of a motor, f 0 , the fraction of actin-attached motors, r and the rate of transition through the attachment-detachment cycle, ϕ .

| Estimated Parameters   | fast            | slow            | ratio   |
|------------------------|-----------------|-----------------|---------|
| f 0 (pN)               | 6 . 8 ± 1 . 0   | 2 . 4 ± 0 . 4   | 2 . 8   |
| r                      | 0 . 32 ± 0 . 02 | 0 . 50 ± 0 . 03 | 0 . 64  |
| ϕ (s - 1 )             | 6 . 0 ± 0 . 2   | 2 . 27 ± 0 . 04 | 2 . 6   |

the results of the analysis are plotted in the parameters plane ( f , r 0 ) (symbols and lines refer to different isoforms according to the colour: blue for psoas, red for soleus; different tones identify different experiments). The solid lines highlight the ensemble of distinct - though equivalent solutions ensuing from the average force profile. By accounting for the fluctuations one breaks the degeneracy inherent to the system when analysed in its mean-field version, getting just one pair ( f , r 0 ) (identified by the symbol) compatible with each individual experimental curve.

One can relax the constraint N = 16 obtained from the rigor experiments (Figure 3) and scan the range of N that yields convergence of the optimisation algorithm, for the imposed level of accuracy. The results of the analysis for the soleus isoform is reported in the Supplementary Note 1 (Supplementary Figure 12), where the estimated f 0 is plotted against N . The histogram computed from the collection of fitted parameters can be conceptualised as an indirect imprint of the degree of experimental variability as associated to f 0 and N .

## Estimated motor force f 0 and fraction of attached motors r from the experimental

## data.

Fig. 5. Best fit parameters from the experimental data sets of rabbit soleus HMMs (red symbols) and rabbit psoas HMMs (blue symbols). Mean values and standard deviations are obtained by averaging over 20 independent realisations of the stochastic fitting procedure for each data record. Different tones refers to different experiments. Each solid line represents the hyperbola on which each of the pair ( f , r 0 ) is constrained to be, according to the mean-field analysis.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-80a4d177-b312-4f61-b6cb-9f29cc6678ca.jpg)

## Discussion

We use the one-dimensional synthetic nanomachine described in 23 to define the isometric mechanical output of an array of 16 myosin motors purified from either fast (psoas), or slow (soleus), muscle of the rabbit. To eliminate the large trap compliance and recover the condition for the motors to operate as independent force generators as in the native half-sarcomere, once the interaction is established the system control is switched from position clamp to length clamp. The array of 16 motors in physiological ATP concentration (2 mM) at 24 ◦ C exhibits a steady isometric force that in the fast isoform is 17 pN, and in the slow isoform is 10 5 pN. . The finding that the force exerted by the same number of motors is 1 6-fold larger in the fast isoform disagrees with the . most common finding in muscles and muscle fibres that the isometric force normalised for the cross sectional area of the fibre T 0 , is either similar or at max 1 5-fold larger in the fast isoform . 10-13,15-19 . Notably in skinned fibres from the same rabbit muscles from which the nanomachine motor proteins are purified, T 0 in psoas at 25 ◦ C has been found 317 ± 14kPa, 1 9-fold larger than . T 0 in soleus, 165 ± 12 kPa 26 .

Recording the development of the steady isometric force in length clamp eliminates the contamination of the large trap compliance, showing a roughly exponential time course characterised by the parameter t r that is 238 ms for the slow isoform and 77 ms for the fast isoform. Thus the rise of the force to the maximum steady value takes a 3-fold longer time for the slow isoform than for the fast isoform. How this emergent property of the motor ensemble relates to the corresponding event in situ and how it is affected by the different isoforms has been tested by comparing the nanomachine output with that of Ca 2+ -activated skinned fibres, from the same rabbit muscles from which the motor proteins were purified. According to the sarcomere-level mechanics for skinned fibres developed in our laboratory, the compliance of the attachments of the skinned fibre segment to the transducer levers is negligible (see the Supplementary Note 2, Supplementary Figure 13). Under these conditions, the force redevelopment following a fast shortening able to drop the isometric force to zero is characterised by a t r = 265 ± 15ms in soleus fibres and 62 ± 5ms in psoas fibres. Thus, the time course of force development recorded by the nanomachine in length clamp and its modulation by the two isoforms, are in quite satisfactory agreement with those recorded at the cell level. The corresponding rates of force development ( a ) by the nanomachine are 28 6 s . -1 for the fast isoform and 9 3 s . -1 for the slow isoform. Considering that in length clamp a is direct expression of the sum of the effective rate constant of attachment/force-generation and the effective rate constant of detachment of the myosin motors, we conclude that the interaction kinetics in isometric condition is 3-fold higher in the fast isoform than in the slow isoform. The attachment/detachment kinetics is expected to increase if the load on the motor ensemble is reduced, due to the strain-dependent increase in rate constant of detachment, which underpins the maximum velocity of shortening ( V 0 ) attained under zero load. V 0 estimated by the time taken by the ensemble to redevelop force following a release able to drop the isometric force to zero (Figure 4a,d) is 0 5 . µ ms -1 and 1 95 . µ ms -1 in the slow and fast HMM respectively, showing a V 0 ∼ 4-fold larger in the fast isoform. Thus, the isoform-dependent increase of V 0 is 33% larger than the increase in a and even larger if one considers that V 0 of the fast isoform is underestimated by the proportionally larger fraction of time spent for the force to drop to zero following the release (compare the records in Figure 4a,d). This suggests that the fast isoform exhibits a specifically larger strain dependence of the detachment rate constant.

The rate of development of the isometric steady force and the force fluctuations superimposed on the steady force in length clamp have been exploited to implement a stochastic three-state model which is able to fit the experimental responses, allowing self-consistent estimates of all the relevant mechanokinetic parameters underlying the isometric performance of the motor ensemble: f 0 , the force of a single correctly oriented motor, r , the fraction of attached motors, and ϕ , the rate of transition through the attachment/detachment cycle (Table 1 in the Results). f 0 of the fast isoform (6 8 . ± 1 0 pN) is 2 8-fold larger than . . f 0 of the slow isoform (2 4 . ± 0 4 pN), while the . ensemble force F 0 is only 1 6 times larger (Figure 4b,e). . This is in a great part explained by the different fraction of attached motors r , which in the fast isoform (0 32 . ± 0 02) is 0 64 that of the . . slow isoform (0 50 . ± 0 03). . The corresponding number of attached motors ( Nr ) is ∼ 5 and ∼ 8 for the fast and the slow isoform respectively. The average force of a single randomly oriented motor (0 55 . f 0 ) is 3 7 pN for the fast isoform and 1 3 pN for the slow isoform, from which the predicted . .

ensemble force is (3 7 . × 5 =) 18 5 pN and (1 3 . . × 8 =) 10 4 pN respectively. . These values are in quite good agreement with the observed values: 17 ± 3 pN ( σ ) for the fast isoform and 10 5 . ± 1 8 pN ( . σ ) for the slow isoform.

The model predicts a rate of transition of a motor through the interaction cycle, and thus a frequency of ATP splitting per motor ( ϕ ) 2 6 times higher for the fast isoform (6 0 s . . -1 ) than for the slow isoform (2 3 s . -1 ) (Table 1). The value of ϕ of the slow isoform array is in a remarkably good agreement with that estimated on the slow muscle of mouse and rat (2 3-2 9 s . . -1 , Supplementary Table 1). On the other hand, ϕ for the fast isoform array is less than half of the one estimated in the fast muscle of the same animals (12 4-13 3 s . . -1 , Supplementary Table 1). The same discrepancy for the isoform-dependent increase in ϕ is found between the model prediction and the skinned fibre experiments (Supplementary Table 1, 15-19 ). However, it must be noted that ( ) the absolute values i of ϕ in skinned fibres is 10-fold smaller than the one in the muscle for both slow and fast myosin isoforms 15-18 ; ( ii ) the difference can only in minor part be explained by the different temperature of the experiments (21-27 ◦ C for the muscle and 12 ◦ C for the skinned fibres), taking into account that the Q 10 of ϕ is ≤ 2 5 in either preparation . 11,19,27 . ϕ predicted by the model for the output of the fast isoform nanomachine is 2 5-fold larger than that predicted for the slow isoform, but . still 2-fold smaller than that indicated by the energy rate measured in the fast muscle. Thus the 5-fold larger ϕ of the fast isoform with respect to the slow isoform found in muscle is only partly explained by higher rate constants of transitions through the conventional attachment/force generation and detachment cycle operating in isometric conditions and recorded by the nanomachine force fluctuations. The actin-activated myosin ATPase activity in solution is 2 5 times larger in . fast than in slow muscle 20 , which can be accounted for by a higher rate of ADP release (which is followed by a fast ATP binding and detachment, step (c)-(d) in Figure 1, 28 ) and/or a higher rate of the hydrolysis step (d)-(e), and/or a higher rate of actin attachment (step (a)-(b)). In isometric contraction at physiological ATP concentration, ADP release is the rate-limiting step for detachment and is 10-fold slower in slow myosin than in fast myosin 28 and this may per se explain the finding that during steady isometric force generation the duty ratio of the fast myosin nanomachine is lower than that of the slow myosin nanomachine. However, it must be taken into account that under isometric conditions (or high load) the transitions through the different force-generating states of the motor (step (a)-(b)/(c) in Figure 1, 29,30 ) slow down due to the strain dependence of the transition rate and thus the subsequent conformation-dependent release of ADP also gets slower (5, 31 ). As far as the difference in ϕ between slow and fast myosin ensembles in isometric contraction, the finding that the force of fast myosin is 2 5-fold higher should suggests . that the equilibrium distribution between different force-generating states is shifted toward the end of the working stroke in the fast myosin, in this way explaining a larger flux through the detachment step and thus the reduction in the duty ratio and the increase in ϕ with respect to the slow myosin (Table 1). However, it must be considered that the stiffness of the myosin motor, determined in situ with fast sarcomere-level mechanics applied to skinned fibres from rabbit muscle, is larger in the fast muscle in proportion to the motor force, so that the extent of the force-generating structural change is the same in either fast or slow myosin motor 32 .

In conclusion, the 2 5-fold larger isometric . ϕ of the fast myosin isoform found with the analysis of force fluctuations is accounted for by an intrinsic faster rate of the relevant kinetic steps of the fast myosin isoform which underpins a 2 5-fold larger ATPase rate in solution . 20 . Instead, the 5-fold larger isometric ϕ of the fast isoform reported in the literature Supplementary Table 1, exceeds by a factor of 2 the one recorded by the nanomachine force fluctuations at 25 ◦ C and could be explained by a further kinetic adaptation of fast myosin isoform hypothesising that, also in isometric conditions, a futile faster actin-activated ATPase cycle is present. In terms of the kinetic scheme in Ref. 33 , this cycle implies the working stroke transition to occur in the motor undergoing weak actin-binding interactions and does not imply strong/force-generating attachment unless the load is reduced and the muscle shortens.

A comparison of the parameters estimated in this work with those obtained in previous nanomechanical approaches is possible for the fast isoform purified from rabbit psoas investigated by Yanagida's group 22 through the microneedle manipulation technique. In close-to-isometric conditions, obtained through a stiff microneedle, both the force of the motor (5 9 pN) and the fraction of actin-attached . motors (0 36) estimated in that work are in exceptional good agreement with the values calculated . here from the output of the nanomachine. A peculiar difference that makes our nanomachine

unique is the possibility to define the performances emerging from the array arrangement of the motors in the half-sarcomere, as the force - velocity relation and the maximum power output. The novelty of the present nanomachine application in relation to the previous ones 23-25 , is the interpretation of the output of the motor ensemble and of the isoform-dependent differences on the basis of the mechanokinetic molecular properties of either isoform is self-consistent way without any assumptions from cell mechanics and solution kinetics. The combined experimental and theoretical achievements in this paper set the stage for any future studies on the emergent mechanokinetic properties of the half-sarcomere like arrangement of any myosin motors, either engineered or purified from mutant animal models or human biopsies.

## Methods

## Preparation of Proteins

Adult male rabbits (New Zealand white strain), provided by Envigo, were housed at Centro di servizi per la Stabulazione Animali da Laboratorio (CeSAL, University of Florence), under controlled conditions of temperature (20 ± 1) ◦ C and humidity (55 ± 10)%, and were euthanized by injection of an overdose of sodium pentobarbitone (150 mg kg -1 ) in the marginal ear vein, in accordance with the Italian regulation on animal experimentation (Authorisation 956 / 2015-PR) in compliance with Decreto Legislativo 26 / 2014 and EU directive 2010 / 63. Three rabbits were used for the experiments. HMM fragments of myosin were purified from rabbit soleus and psoas muscles as reported previously in 23,24 . The functionality of the purified motors was always preliminarily checked with in vitro motility assay. Actin was prepared from leg muscles of the rabbits according to 34 , and polymerised F-actin was fluorescently labelled by incubating it overnight at 4 ◦ C with an excess of phalloidin-tetramethyl rhodamine isothiocyanate 35 . For the mechanical measurements in the nanomachine, the correct polarity of the actin filament was pursued by attaching the + end of the filament to a polystyrene bead (3 µ m diameter) (Bead-Tailed Actin, BTA, 36 ) with either the Ca 2+ -sensitive capping protein gelsolin 23 or the Ca 2+ insensitive gelsolin fragment TL 40 (Hypermol, Germany) 24,25 ).

## Mechanical experiments

The mechanical apparatus, described in detail in 23 , is depicted in Figure 2. HMM fragments of myosin were deposited randomly on the lateral surface of a glass pipette pulled to a final diameter of ∼ 3-4 µ m and functionalised with nitrocellulose 1 % ( w/v ). The glass pipette was mounted in the flow chamber carried on a three-way piezoelectric nanopositioner (nano-PDQ375, Mad City Lab, Madison WI, USA) that acts as a displacement transducer, and was brought to interact with a BTA trapped in the focus of a Dual Laser Optical Tweezers (DLOT) that acts as a force transducer. The DLOT system has a dynamic range for both force (0 -200 pN, resolution 0 3 pN) and displacement . (0 -75 µ m, resolution 1 1 nm) adequate for the measuring the output of the nanomachine. . The buffer solutions used for all the experiments are already reported in 23 and contained physiological concentrations of ATP (2 mM) unless differently specified. 0 5 % methylcellulose was added to the . running buffer in order to inhibit the lateral diffusion of F-actin 37 and minimise the probability of loss of acto-myosin interaction. The concentration of HMM from soleus and psoas muscle used for the experiments was defined by the concentration at which the number of rupture events in rigor attained a saturating value.

The mechanical apparatus, as already reported in 23 , can be operated either in position clamp (Figure 2, red branch), when the feedback signal is the position of the nanopositioner carrying the motor array ( x ), or in force clamp (green branch), when the feedback signal is the force ( F ), calculated as the product of the stiffness of the trap ( e ) times the change in position of the bead in the laser trap ( x bead ). Recording of the nanomachine performance in true isometric condition, however, cannot be achieved in position clamp, due to the large trap compliance ( ∼ 4nmpN -1 ), which implies both several tens of nanometres movement to develop the maximum steady force and blunting of the force of individual attachment-detachment events (Supplementary Figure 7 in 23 ). To eliminate the trap compliance the system has been implemented with a length clamp (blue branch in Figure 2), which uses as a feedback signal the change in distance ( L ) between the position of the actin attached bead in the laser trap ( x bead ) and that of the nanopositioner ( x ), so that the movement of the bead with the force change is counteracted by the movement of the nanopositioner. In this way the effective trap compliance is reduced to 0 2 nm pN . -1 .

In length clamp the frequency response of the system is reduced by the propagation time of the mechanical signal through the loop from the force transducer to the nanopositioner, which also includes the array of actin attached myosin motors. The power density spectrum (PDS) of the system, measured with sinusoidal oscillations at different frequencies, changes depending on the selected feedback mode: in position clamp the PDS shows an upper -3dB frequency (or corner frequency f c ) of 59 Hz (Figure 6, red); in length clamp, when the feedback loop is closed with the array of actin-attached myosin motors in rigor, f c decreases to 32 Hz with HMM from fast muscle (violet) and to 17 Hz with HMM from slow muscle (cyan). The mass of the system ( m ) is the same

Power density spectrum of the system.

Fig. 6. Superimposed power density spectrum either in position clamp (red circles interpolated by the red Lorentzian curve), or in length clamp with the array of actin attached motors in rigor from both fast muscle (violet circles and curve) and slow muscle (cyan circles and curve). The upper -3 dB frequency f c is: 59 ± 3 Hz (red), 31 ± 6 Hz (violet) and 17 ± 3 Hz (cyan). The coloured area delimited by thinner lines indicates the confidence limits.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-bdfe5572-9ac9-4251-b916-6402b376b416.jpg)

with either isoform array thus the different corner frequency between the two nanomachines should almost in part depend on the different stiffness of the two arrays in rigor.

The architecture of the machine (with the length of the motor array much shorter than the length of the overlapping actin filament) implies that, for a given HMM concentration, the measured number of rupture events does not significantly change from experiment to experiment, therefore there is no need to normalise the mechanical response obtained in different experiments at physiological [ATP] by the actin-filament length. All the experiments were conducted at room temperature (24 ◦ C).

## Statistics and reproducibility

Data are expressed as mean ± standard deviation unless otherwise stated. The number of replicates is defined in the text and in the figure legends.

## Stochastic model: on the governing master equation

The master equation can be cast in the general form:

̸

$$\frac { \partial P ( n, t ) } { \partial t } = \sum _ { n ^ { \prime } \neq n } \left [ T ( n | n ^ { \prime } ) P ( n ^ { \prime }, t ) - T ( n ^ { \prime } | n ) P ( n, t ) \right ]$$

where T ( n n ′ | ) represent the transition rates from the state n to a new state n ′ , compatible with the former. In the following, to identify the arrival/departure state n ′ we solely highlight the individual component of the vector n that changes due to the considered reaction. The explicit expression for the transition rates that originates from the chemical equations (1) is given in the annexed Supplementary Note 1.

In the Supplementary Note 1 are also presented the details of the numerical simulation of a single stochastic orbit of the considered dynamics, obtained via the celebrated Gillespie algorithm 38,39 . In Supplementary Figure 3 it is shown the time evolution of the (discrete) concentration of molecular motors in each configuration, while in Supplementary Figure 4 the results of the stochastic simulations for the probability distributions of the fractions of motors in the force-generating configurations are compared with the theoretical prediction.

## The deterministic limit

From the master equation one can readily derive the mean field equations that governs the deterministic dynamics for the continuous concentrations of the molecular motors in configurations A 1 and A 2 :

$$\begin{cases} \frac { d y } { d t } & = k _ { 1 } - \left ( k _ { 1 } + k _ { - 1 } + k _ { 2 } \right ) \, y - \left ( k _ { 1 } - k _ { - 2 } \right ) \, z \\ d z & = k _ { 2 } \, y - \left ( k _ { - 2 } + k _ { 3 } \right ) \, z \end{cases} \\ \intertext { z \text{ identifier to aroad fraction of the molecular matrices in state } \, \Delta \, \text{ and } \, \Delta \, \text{ recognize}$$

where y and z identify the averaged fraction of the molecular motors in states A 1 and A 2 , respectively. Equations (4) can be studied at equilibrium, yielding the fixed point solutions:

$$\begin{cases} y ^ { * } = & \left ( \frac { k _ { 1 } } { k _ { 1 } + G } \right ) \frac { k _ { - 2 } + k _ { 3 } } { k _ { 2 } + k _ { - 2 } + k _ { 3 } } \\ \\ z ^ { * } = & \left ( \frac { k _ { 1 } } { k _ { 1 } + G } \right ) \frac { k _ { 2 } } { k _ { 2 } + k _ { - 2 } + k _ { 3 } } \end{cases}$$

where:

$$G = \frac { k _ { - 1 } ( k _ { - 2 } + k _ { 3 } ) + k _ { 2 } k _ { 3 } } { k _ { 2 } + k _ { - 2 } + k _ { 3 } }$$

We define the duty ratio r as the fraction of attached motors (or the fraction of the ATPase cycle time a motor spends attached). It can be computed as:

$$r = y ^ { * } + z ^ { * } = \frac { k _ { 1 } } { k _ { 1 } + G } \.$$

The temporal evolution of the mean field concentrations of motors in different configurations, for a suitable choice of the kinetics parameters is shown in the Supplementary Figure 2, in Supplementary Note 1. A straightforward calculation can be performed to show that z ∗ ≃ r = k 1 k 1 + G and y ∗ ≪ 1, when k -2 /k , k 2 3 /k 2 ≪ 1. In practical terms, under this operating assumption, which for the mammalian muscle myosin under consideration is approached at temperature T ≃ 24 ◦ C, motors are solely found in state A 2 . As discussed earlier, this is the relevant setting for the specific case study at hand.

Equations (4) can be drastically simplified by performing a self-consistent elimination of the variable y . To this end we set dy/dt = 0 in the first of equations (4) to eventually express y as a function of z . This procedure is customarily invoked to carry out the so called adiabatic elimination, which proves correct when there is a clear separation of time scales between co-evolving variables. Although this is not a priori the case for the system at hand, we will postulate the validity of the aforementioned condition and operate with the ensuing approximation that, as we shall prove, will materialise in an accurate interpretative framework. Further details can be found in the Supplementary Note 1 (see Supplementary Figure 1). Plugging the expression for y as a function of z into the second of equations (4) and solving the ensuing differential equation readily yields:

$$\frac { d z } { d t } = \frac { k _ { 1 } \ k _ { 2 } } { k _ { 1 } + k _ { - 1 } + k _ { 2 } } - z \, \left [ k _ { - 2 } + k _ { 3 } + \frac { k _ { 2 } \ ( k _ { 1 } - k _ { - 2 } ) } { k _ { 1 } + k _ { - 1 } + k _ { 2 } } \right ] \equiv b - a z$$

which immediately yields solution (2), as reported in the Results. From equation (2) we can write z ∗ = b/a and this latter condition matches the homologous estimate derived from the original two dimensional model and reported in equations (5).

## Exact solution of the stochastic problem

At first, we remark that the probability distribution P ( n ; t ) ≡ P n , n ( 1 2 ; t ) can be written as a vector P ( ) t of dimension ( N + 1) × ( N + 1). This latter returns the probability at time t , of finding the system in the state characterised by n 1 motors in configuration A 1 and n 2 motors

in configuration A 2 . Here, n 1 and n 2 can in principle assume every integer values in the range [0 , N ], i.e. a total of N +1 values each. Observe however that the populations of motors in the configurations A 1 and A 2 , must satisfy the obvious constraint that reflects the conservation law, i.e. n 1 + n 2 ≤ N . A simple way to express the condition above is to consider that for each possible value of n 1 , n 2 can assume values in the range [0 , N -n 1 ]. This readily implies that the total number of possible states for the system is identically equal to M = ( N +1)( N +2) 2. The number of allowed / states are hence smaller than what anticipated above. Indeed the non trivial entries of P ( ) are t M = ( N +1)( N +2) 2. We will consequently focus on the non trivial elements of vector / P ( ) which t we shall denote as P m ( ) with t m = 1 , . . . , M . For the relevant case of N = 16 molecular motors, instead of ( N +1) × ( N +1) = 289 configurations we only have M = 153 possible states that can be eventually visited by the system, and that we explicitly list in the Supplementary Note 1. As we shall also discuss in the same Supplementary Note 1, the stationary probability distribution P st defines the kernel of a M × M matrix Q and can be hence computed as the eigenvector of Q relative to the null eigenvalue. The entries of the matrix Q can be computed from the transition rates of the underlying master equation, as highlighted in the Supplementary Note 1.

The marginal probability ρ k to find k ≤ N motors in A 2 can be extracted from the stationary probability distribution P st , the stationary solution of the master equation. This is done by summing the elements of P st that refer to the selected k , and that account for all possible partitioning of the remaining N -k motors among configurations D and A 1 . The knowledge of the stationary probabilities ( ρ , ρ 0 1 , ρ 2 , ..., ρ N ) opens up the perspective to calculate the stationary state distribution of the applied force F , an important asset when aiming at a refined fitting scheme which meticulously accounts for the role played by fluctuations.

To work along these lines we shall assume that the contribution to the force (including fluctuations) of the motors in the state A 1 is always negligible. This assumption is motivated by the fact that, for the experimental setting here explored, only a tiny fraction of motors is found to populate state A 1 , at any time t . In the Supplementary Note 1, we will however relax this working assumption so as to provide a rigorous theoretical framework that extends to account for the relevant setting where the population of A 1 motors is instead significant in size.

Let us focus on k ≤ N distinct motors in state A 2 . As postulated earlier, each motor can exert a constant random force f , uniformly spanning the assigned interval I 2 . For each choice of k , we can compute the distribution of the forces Π k ( f ) applied by the selected k motors, by combining independent and identically uniformly distributed random variables drawn for the interval of pertinence I 2 (see the Supplementary Note 1 for further technical information, Supplementary Figure 6 and Supplementary Figure 7). As stated in the Results, functions Π k ( f ) need to be combined together with proper weighting factors that reflect the stationary probability ρ k of having exactly k motors in the force-generating state A 2 , namely P F ( ) = ∑ N k =0 ρ k Π ( k f ).

The relevant steps of the fitting strategy are listed below. We begin by focusing on the average force profile, hence disregard the impact of finite size fluctuations. As mentioned above, the time evolution of the recorded force can be approximated by an effective, two-parameters model (see Supplementary Figure 5 in the Supplementary Note 1). The latter parameters - respectively denoted F 0 and a - can be estimated via a direct fit. Having accessed to preliminary estimated values for the average force at the stationary plateau F 0 and for the rate of isometric force development a , one can then set forth to characterise the other kinetic parameters by analysing the distribution of the fluctuations of the force around the asymptotic plateau. To this end we note that f 0 , one of the unknown of the model, can be written as:

$$f _ { 0 } = \frac { 2 0 } { 1 1 } \frac { F _ { 0 } } { N } \frac { a } { b }$$

where a is constrained to the value determined above while b = k k / k 1 2 ( 1 + k -1 + k 2 ) as defined by equation (8).

Armed with the above knowledge, we can proceed further by comparing the probability density function of the force fluctuations as obtained analytically to the homologous histogram computed from the stochastic simulations. The former is adjusted to the latter by modulating the free parameters k 1 , k -1 , k 2 , k -2 and k 3 , for a fixed choice of N . As discussed in the Supplementary Note 1, testing the method against synthetic data generated in silico enables us to conclude that parameters f 0 and r = k / k 1 ( 1 + G ) can be correctly estimated, following the above fitting scheme

(see Supplementary Table 2). Also the estimated b and a (recomputed from the best fitted values for the kinetic constants) are pretty close to their nominal values as imposed in the simulations. Remarkably ϕ , the rate of motors completing the interacting cycle with the actin, is also correctly recovered. A graphic comparison between estimated and exact parameters (i.e. those employed in the inspected simulations) is also shown in Supplementary Figure 8 of the Supplementary Note 1. Notice however that multiple combinations of the parameters k 1 , k -1 , k 2 , k -2 and k 3 exists that yields the same fitted profile (with almost identical estimates for the relevant quantities f 0 , r , a and b ).

While the kinetics of the scrutinised model cannot be solved unequivocally, we are in a position to accurately determine crucial parameters - as e.g. the maximum force exerted by a single motor and the associated duty ratio of the ensemble - which proved elusive under the deterministic viewpoint, as can be appreciated in Supplementary Figure 9 of the Supplementary note 1.

The above analysis refers to a fixed value of N , the size of the system that we assumed (from the experiment results shown in Figure 3) N = 16. In principle the correct value of N is not a priori known. To overcome this intrinsic limitation, one could repeat the analysis by varying N and recording the parameter estimated as follows the fitting scheme. Here, we will consider the simplified setting where a and b are frozen to the values determined for N = 16 (so that z ∗ stays unchanged). This is implemented by removing two parameters from the pool of quantities to be fitted. Specifically k -1 and k 3 are constrained to match two constitutive relations, that involve k 1 , k 2 and k -2 , in addition to a and b . The parameters to be fitted are hence k 1 , k 2 and k -2 , while k -1 , k 3 and f 0 can be self-consistently determined from the their best fit values. Notice that f 0 is expected to change as a function of N as specified by relation (9). The result of the analysis are reported in Supplementary Figure 10 in the Supplementary Note 1: the fitting procedure converges (with the requested limit of precision) only over a finite range of values of N , centred around the value adopted when performing the simulations. This observation implies in turn that we are in a position to obtain a reasonable estimate for the interval of pertinence of N , as follows the procedure outlined above.

The introduced theoretical framework and the ensuing fitting strategy, thoroughly validated against synthetic data, can be hence applied to the analysis of the experimental data so to yield a selfconsistent estimate of the underlying mechanokinetic parameters. For a detailed validation of the fitting scheme against synthetically generated data refer to the Supplementary Note 1.

Data collection and analysis: A custom built program written in LabVIEW (National Instruments) was used for signal recording. Data analysis was carried out using LabVIEW (National Instruments) and MATLAB (MathWorks) dedicated scripts, and Origin 2018 (OriginLab Corp., Northampton, MA, USA) and Igor Pro 8 (WaveMetrics, Portlan, OR, USA) software.

Data availability: The source data for all figures and tables are provided as Supplementary Data 1. All remaining data will be available from the corresponding authors upon reasonable request.

## Acknowledgments

We thank Gabriella Piazzesi for critical evaluation of this work. We thank the staff of the mechanical workshop of the Department of Physics and Astronomy (University of Florence) for mechanical engineering support. This work was supported by Fondazione Cassa di Risparmio di Firenze (2020.1582 (Italy), the Italian Ministry of Education, Universities and Research (PRIN2020, CUP: B55F22000540001), the European Joint Program on Rare Diseases 2019 (PredACTINg, EJPRD19033), and the Next Generation EU programme [DM 1557 11.10.2022] in the context of the National Recovery and Resilience Plan, Investment PE8-Project Age-It: 'Ageing Well in an Ageing Society', Italy.

## Author contributions

PB, VL and DF designed the research. IP, PB, ML, MC, MM, IM performed the experiments and analysed the data. VB and DF developed the stochastic model. VL, DF, VB, PB and MR wrote the paper or revised it critically for important intellectual content. All authors participated in discussions on this work and approved the final version of the manuscript.

Competing interests: The authors declare no competing interests.

## Supplementary Information

## Supplementary Table 1

Supplementary Table 1. Energetic cost of isometric contraction in slow and fast mammalian skeletal muscle in terms of rate of ATP hydrolysis per myosin motor, ϕ , and ratio of ϕ of fast over ϕ of slow muscle ( R ). Data from the literature (ref, in brackets) except the last row that reports the data from the nanomachine simulation

| ϕ (s - 1 )                  | [ref.]   | fast   | slow   | R      |
|-----------------------------|----------|--------|--------|--------|
| Mouse muscle 21 ◦ C         | [10]     | 12 . 4 | 2 . 47 | 5 . 0  |
| Mouse muscle 25 ◦ C         | [11]     | 13 . 3 | 2 . 95 | 4 . 5  |
| Rat muscle 27 ◦ C           | [13]     | 12 . 5 | 2 . 3  | 5 . 4  |
| Rat skinned fibre 12 ◦ C    | [15]     | 1 . 28 | 0 . 25 | 5 . 12 |
| Rabbit skinned fibre 12 ◦ C | [18]     | 1 . 79 | 0 . 23 | 7 . 78 |
| Human skinned fibre 12 ◦ C  | [19]     | 3 . 22 | 0 . 65 | 4 . 95 |
| Nanomachine 24 ◦ C          |          | 6 . 0  | 2 . 3  | 2 . 6  |

## Supplementary Note 1

## Additional information on the stochastic model

The quantities T ( n n ′ | ), as introduced in the main body of the paper, stand for the transition probabilities from the state n to the state n ′ . In the following, to identify the arrival/departure state n ′ we solely highlight the individual component of the vector n that changes due to the considered reaction. The transition rates as stemming from the chemical equations that define the stochastic model under investigation take the following explicit form:

$$\text{stochastic model under investigation take the following explicit form: } \\ ( 1 ) \ A T T A C H M E N T \quad \ \ T ( n _ { 1 } + 1 | n ) = k _ { 1 } \ \frac { n _ { D } } { N } = k _ { 1 } \left [ 1 - \frac { 1 } {$$

The governing master equation can be hence written in the following explicit form:

$$\text{$\text{governing master equation can be hence written in the following explicit form:} \\ \frac { \partial P ( n, t ) } { \partial t } = \left [ T ( n | n _ { 1 } - 1 ) P ( n _ { 1 } - 1 ; t ) - T ( n _ { 1 } + 1 | n ) P ($$

## Additional information on the adiabatic approximation

As anticipated in the main text, the validity of the adiabatic elimination was postulated even in the absence of a clear time scales separation of the dynamics of the system, and here we test the validity of this approach, for the operating conditions, in reproducing the correct time scale of the dynamics of the concentration for motors in A 2 . In the Supplementary Figure 1 it is shown the evolution of the concentration of the motors in the force generating configuration A 2 , from a initial condition with all the motors detached from the actin, from which it is possible to estimate the time scale on which the stationary state is approached.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-1ed02a5a-5aa8-4a3b-b095-50606ea61981.jpg)

## Supplementary Figure 1. Checking the adiabatic approximation.

In red (solid and dashed lines) the evolution of the concentration of motors in configuration A 1 , in blue (solid and dashed lines) the evolution of the concentration of motors in configuration A 2 . The solid lines represent the integration of the dynamical system without any approximation, the dashed lines correspond to the integration of the dynamical system when the adiabatic approximation is adopted. The approximation allows to reproduce the correct the time scale at which the concentration of motors in A 2 reaches the stationary state.

## Additional information on the structure of the stochastic matrix

As remarked in the main body of the paper for each possible value of n 1 , n 2 can assume values in the range [0 , N -n 1 ]. Hence, the total number of possible states of the system is M = ( N +1)( N +2) 2. / Consider the relevant setting N = 16. We thus have to deal with M = 153 possible states as listed in the following:

$$\text{use} \, \stackrel { \underset { m \in [ 1, 1 5 3 ] } { m \in [ 1, \, n _ { 2 } ) } } { \underset { 1 \, 1 \, 1 } { m \, \stackrel { \underset { m \, \stackrel { \, \, \,$$

| m ∈ [1 ,   | ( n 1 , n 2 )         |                 |
|------------|-----------------------|-----------------|
| 1 2 . .    | (0 , 0) (0 , 1) . . . |      k = 0 |
| . 17       | (0 , 16)              |            |
| 18 . .     | (1 , 0) . .           |   k = 1       |
| . 33       | . (1 , 15)            |               |
| 34 . . .   | (2 , 0) . . .         | . . .           |
| 150        | (14 , 2)              |                 |
|            |                       | }               |
| 151        | (15 , 0)              | k = 15          |
| 152        | (15 , 1)              | }               |
| 153        | (16 , 0)              | k = 16          |

Focus now on the generic element W lm that enters the definition of matrix Q . Assume in particular m to label the reference initial state; l identifies the state that can be eventually reached through the chemical dynamics. Five possible types of transitions exist, organised in k = N +1 blocks, which corresponds to the selected value of n 1 , while n 2 can freely varies within [0 , N -n 1 ]:

$$\text{ecu value on $n_{1},$ where $n_{2}$ can merely $\alpha$} \\ l _ { 1 } & = m ^ { \max } ( k + 1 ) + n _ { 2 } + 1 \\ l _ { 2 } & = m ^ { \max } ( k ) - n _ { 2 } ^ { \max } + n$$

where we denote m max ( k ) the largest possible index as associated to block k , for the selected choice of m n ; max 2 stands the largest values that can eventually take the discrete variable n 2 . Hence:

$$n _ { 2 } ^ { m a x } & = N - n _ { 1 } \\ m ^ { m a x } ( k + 1 ) & = m ^ { m a x } ( k ) + N + 1 - n _ { 1 }$$

The non trivial elements W m,l are hence:

$$( \begin{matrix} 1 \\ 1 \end{matrix} )$$

$$\text{The non trivial elements } W _ { m, l } \text{ are hence} \colon \\ W _ { l, m } & \equiv T ( n _ { 1 } + 1, n _ { 2 } | n _ { 1 }, n _ { 2 } ) = \frac { k _ { 1 } } { N } ( N - n _ { 1 } - n _ { 2 } ) \quad \begin{cases} n _ { 2 } \in [ 0, n _ { 2 } ^ { m a x } - 1 ] \\ n _ { 1 } \in [ 0, N - 1 ] \end{cases} \\ W _ { l, 2, m } & \equiv T ( n _ { 1 } - 1, n _ { 2 } | n _ { 1 }, n _ { 2 } ) = \frac { k _ { - 1 } } { N } n _ { 1 } \quad \begin{cases} n _ { 2 } \in [ 0, n _ { 2 } ^ { m a x } ] \\ n _ { 1 } \in [ 1, N ] \end{cases} \\ W _ { l, 3, m } & \equiv T ( n _ { 1 } - 1, n _ { 2 } + 1 | n _ { 1 }, n _ { 2 } ) = \frac { k _ { 2 } } { N } n _ { 1 } \quad \begin{cases} n _ { 2 } \in [ 0, n _ { 2 } ^ { m a x } ] \\ n _ { 1 } \in [ 1, N ] \end{cases} \\ W _ { l, 4, m } & \equiv T ( n _ { 1 } + 1, n _ { 2 } - 1 | n _ { 1 }, n _ { 2 } ) = \frac { k _ { - 2 } } { N } n _ { 2 } \quad \begin{cases} n _ { 2 } \in [ 1, n _ { 2 } ^ { m a x } ] \\ n _ { 1 } \in [ 0, N - 1 ] \end{cases} \\ W _ { l, 5, m } & \equiv T ( n _ { 1 }, n _ { 2 } - 1 | n _ { 1 }, n _ { 2 } ) = \frac { k _ { 3 } } { N } n _ { 2 } \quad \begin{cases} n _ { 2 } \in [ 1, n _ { 2 } ^ { m a x } ] \\ n _ { 1 } \in [ 0, N ] \end{cases}$$

Given the above structure, it is possible to identify for every choice of m , the corresponding combination of n 1 and n 2 , and associate the m -component of vector P ( ) t to a specific state ( n , n 1 2 ). That is because the m max ( k ) are the partial sums of the finite sequence:

$$m ^ { m a x } ( k ) = \sum _ { i = 0 } ^ { k } N + 1 - i \text{ \quad \ \ } i \in \{ 0, \dots, N \}$$

for k = 0 , . . . , N .

For a given m , we thus identify the index k that matches the relation:

$$m ^ { m a x } ( k ) \geq m$$

and then set:

$$\begin{cases} n _ { 1 } = k & \\ n _ { 2 } = m - m ^ { m a x } ( k - 1 ) - k & \end{cases}$$

The stationary solution of the stochastic dynamics is found as the eigenvector of Q associated to the null eigenvalue.

## Deterministic and stochastic simulations of the underlying model

Consider first the set of differential equations that represents the mean field approximation for the examined system. These equations can be numerically integrated, for a representative choice of the chemical reaction parameters, so as to resolve the temporal evolution of the concentrations of the motors, in configurations D A , 1 and A 2 . As shown in Supplementary Figure 2, after a transient, the concentration of the motors in the various populations approach their stationary state x ∗ , y ∗ , z ∗ .

We now turn to studying the stochastic model that implements the microscopic dynamics described by the chemical equations (1) in the main text. It is possible to numerically simulate a single stochastic orbit of the considered dynamics (working with the same kinetic constants as assumed in the the mean field simulations) via the celebrated Gillespie algorithm. In doing so one eventually obtains the time evolution of the (discrete) concentration of molecular motors in each configuration, i.e. n D ( ) t /N , n 1 ( ) t /N , n 2 ( ) t /N . A typical solution is displayed in Supplementary Figure 3. As expected the stochastic trajectories fluctuate around the corresponding deterministic orbit (solid lines). The observed fluctuations are a material imprint of the inherent discreteness of the simulated system.

The distribution of fluctuations can be numerically accessed from individual stochastic simulation, by averaging over a large set of independent stochastic realisations. In Supplementary Figure 4 the equilibrium distribution of the fluctuations (i.e. the fluctuations displayed around the fixed point, once the initial transient has faded away) is depicted and compared to the analytical solution obtained from the governing master equation, via the procedure discussed above. The agreement is satisfying and testifies on the correctness of the proposed analytical treatment.

We are now in the position to simulate also the temporal series of the force of the ensemble, that is obtained by assigning to each individual motor the force that it is able to exert, based on its configuration (as stipulated by the stochastic dynamics), and following the prescriptions described in the Results. The time series of the force exerted by motors of the populations A 1 and A 2 is displayed in Supplementary Figure 5, for the specific choice of the parameters here operated.

By accessing the temporal evolution of the force, including the equilibrium fluctuations, one can recover key information on the underlying structural and chemical parameters. As shown in the main body of the paper, this corresponding to solving an inverse problem, from the observed force back to the relevant parameters.

Supplementary Figure 2. Temporal behaviour of the mean field concentrations of motors in different states.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-77cf4886-0a8f-4f88-b2ba-7be9b40a6652.jpg)

The solid yellow line is the evolution of motors in the detached configuration D (yellow); the solid red and blue lines are the evolution of motors in the configuration A 1 and A 2 respectively. The dotted lines correspond to the stationary state values x ∗ , y ∗ , z ∗ .

Supplementary Figure 3. Stochastic simulations obtained with the Gillespie algorithm. Temporal behaviour of the concentration of the three populations of motors in each configuration: D (yellow line), A 1 (red line), A 2 (blue line). The black solid lines represents the deterministic evolution of the concentrations.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-ad285c5c-f6e8-4a25-a920-efaa3f84b0a8.jpg)

Supplementary Figure 4. Checking the theory predictions vs. stochastic simulations. Comparison between the probability distributions of the concentrations of motors in the forcegenerating configurations as obtained from the simulated dynamics and the stationary solution of the master equation. Normalised histograms (light coloured lines) refers to the populations of motors in configurations A 1 (in red) and A 2 (in blue), as obtained from the simulated dynamics for a suitable choice of the kinetic parameters. Dark coloured lines (with symbols) stand for the homologous distribution as derived from the stationary solution of the master equation.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-b1fa1c3f-0eaa-45fe-a277-21fb73b34e97.jpg)

Supplementary Figure 5. Force development as obtained by averages of stochastic simulations.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-9eaed640-ca69-4611-af8a-5ac9d4a1b60f.jpg)

The trajectory has been obtained averaging over 100 simulations. The force is measured in pN and it is exerted by a collection of N = 16 molecular motors with a suitable choice of the kinetic parameters.

## Analytical characterisation of the force probability distribution

As discussed in the main body of the paper, the fluctuations of the force around the average value stem from finite size corrections. To estimate the distribution we focus on k ∈ [1 , N ] force-generating distinct motors and postulate that each of them can exert a uniform, randomly selected force f . For each choice of k , one can compute the distribution of the force Π k ( f ) exerted by the k motors. This is a particular case of the more general problem of calculating the probability distribution for the total force exerted by all the n 1 and n 2 molecular motors in the force-generating configurations A 1 and A 2 . In this section we discuss the problem under this general perspective. To this end we denote by f 1 the random force uniformly distributed in the interval I 1 and exerted by the motors in configuration A 1 and by f 2 the one extracted from the interval I 2 , exerted by the motors in the configuration A 2 . As suggested in 40 we are dealing with the problem of finding the probability distribution of the sum of n random variables Y i , for i = 1 , . . . , n each of them uniformly distributed in the interval [ b , c i i ]. In our case we have only two classes of variables Y i : the ones relative to the forces exerted by motors in the configuration A 1 , and those relative to the forces exerted by motors in the configuration A 2 , so that:

$$\sum _ { i = 1 } ^ { n } Y _ { i } = \sum _ { i = 1 } ^ { n _ { 1 } } Y _ { i } + \sum _ { i = n _ { 1 } + 1 } ^ { n } Y _ { i }$$

and:

$$\text{and} \colon \\ Y _ { i } = \begin{cases} f _ { 1 } & \text{and} & \quad [ b _ { i }, c _ { i } ] = [ - f _ { 0 }, f _ { 0 } ] & \text{if} \quad i = 1, \dots, n _ { 1 } \\ f _ { 2 } & \text{and} & \quad [ b _ { i }, c _ { i } ] = \left [ \frac { f _ { 0 } } { 1 0 }, f _ { 0 } \right ] & \text{if} \quad i = n _ { 1 } + 1, \dots, n \\ \dots & \quad. \quad. \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{cases}$$

We can observe that the probability distributions of the sum of the variables Y i is the same as the probability distribution of the sum of the variables X i = Y i -b i , which are defined in the intervals [0 , a i ], where a i = c i -b i . Introduce the sum s :

$$s = \sum _ { i = 1 } ^ { n } X _ { i }$$

with n ≥ 2. Then the distribution of the sum as defined above reads 40 :

$$\Pi _ { n } ( s ; n _ { 1 }, n _ { 2 } ) = \frac { 1 } { ( n - 1 )! } \frac { 1 } { ( a _ { 1 } ) ^ { n _ { 1 } } ( a _ { 2 } ) ^ { n _ { 2 } } } \left [ \ s ^ { n - 1 } + \sum _ { k = 1 } ^ { n } ( - 1 ) ^ { k } \left ( \sum _ { J _ { k } } \left ( s - \sum _ { l = 1 } ^ { k } a _ { j _ { l } } \right ) _ { + } \right ) ^ { n - 1 } \right ] \quad ( 1 5 )$$

where we adopted the notation: ( f ) + = max { 0 , f } .

From this expression it is possible to compute the probability distribution for the sum of the variables of our interest:

$$\sum _ { i = 1 } ^ { n } Y _ { i } = \sum _ { i = 1 } ^ { n } \left ( X _ { i } + b _ { i } \right ) = \sum _ { i = 1 } ^ { n } X _ { i } + \sum _ { i = 1 } ^ { n _ { 1 } } b _ { i } + \sum _ { i = n _ { 1 } + 1 } ^ { n } b _ { i } = s + n _ { 1 } b _ { 1 } + n _ { 2 } b _ { 2 }$$

where b 1 = -f 0 and b 2 = f 0 10 .

If we consider just one class of variables Y i , meaning if we consider only the forces exerted by one of the two force-generating populations of motors, the distributions Π n can be computed as a specific case of the generalisation of the Irwin-Hall distribution, 41 , the uniform sum distribution. These refer to the sum of n random variables x i , each of them defined in the interval [ a, b ] and take the form:

$$\Pi _ { n } ( x ) = \frac { 1 } { b - a } \ g ( y ; n ) & & \text{ with} \colon & & y = \frac { x - n a } { ( b - a ) } & & ( 1 6 )$$

where:

$$g ( y, n ) = \frac { 1 } { 2 ( n - 1 )! } \sum _ { k = 0 } ^ { n } ( - 1 ) ^ { k } \binom { n } { k } ( y - k ) ^ { n - 1 } \, \text{sgn} ( y - k ) \.$$

In order to obtain the correct expression for the probability distribution of the total force exerted by both n 1 motors in the configuration A 1 and n 2 motors in the configuration A 2 at a given time, one has to weight the above expressions with the probability to find the system in the considered condition. This latter follows the stationary solution of the master equation as discussed in the main body of the paper. From P st , we can in particular extract the probability ρ k to find k ≤ N motors in A 2 . This is obtained by summing the elements of P st that refer to a given k , thus accounting for all possible partitioning of the remaining N -k motors among states D and A 1 (see also discussion in the main body of the paper). Labelling these marginal probabilities ( ρ , ρ 0 1 , ρ 2 , ..., ρ N ) one gets:

$$P ( F ) = \rho _ { 1 } \Pi _ { 1 } ( f ) + \rho _ { 2 } \Pi _ { 2 } ( f ) + \cdots + \rho _ { 1 6 } \Pi _ { 1 6 } ( f ) \ = \sum _ { i = 1 } ^ { N } \rho _ { i } \ \Pi _ { i } ( f ) \.$$

where use has been made of the fact that Π 0 = 0.

The obtained profiles are reported in Supplementary Figure 6 for the relevant case N = 16 and for i = 1 2 3 4. , , , In Supplementary Figure 7 the global distribution of fluctuations is depicted for a specific choice of the kinetic parameters of the model.

Supplementary Figure 6. Many body force distributions.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-03801ca4-3c26-4f26-a689-27940db5e3c0.jpg)

The distributions Π ( k f ) are plotted for k = 1 2 3 4 , , , .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-b258c28e-aed4-4319-abf9-8ea6314e6ae9.jpg)

## Supplementary Figure 7. Theoretical distribution of force fluctuations.

Probability density function P F ( ) as resulting from the sum of the Π ( k f ) , for k = 1 , . . . , N weighted with the stationary solution of the master equation for a system of N = 16 . The adopted parameters are: f 0 = 6 pN , k 1 = 30 s -1 , k -1 = 500 s -1 k 2 = 2000 s -1 k -2 = 100 s -1 k 3 = 10 s -1 .

## Parameters estimation

The fitting procedure described in Results and applied to the experimental data set of soleus and psoas HMM has been validated against synthetic data, and the details are presented in this Section. The parameters f 0 and r are correctly estimated, as it follows from inspection of Supplementary Table 2, where the error associated with the estimated value of r is the standard deviation. Also the estimated b and a (recomputed from the best fitted values for the kinetic constants) are pretty close to their nominal values as imposed in the simulations, as well as the value of the rate ϕ . Supplementary Figure 8 shows the graphic comparison between the results of the optimisation procedure for the parameters estimation performed under the assumptions of the effective model, and the values of the parameters adopted in the stochastic simulations of the dynamics (these latter are referred to as true parameters, as they are imposed in the simulation and thus known with infinite precision). The symbols refers to the mean and SD of the duty ratio r , the force of a single motor f 0 , the rate ϕ and the two parameters a and b that describe the dynamics of the system in the deterministic limit. The blue line represents the value of the corresponding parameter adopted in the simulations.

Supplementary Figure 9 refers to the parameters estimation in the plane ( f , r 0 ), where the solutions resulting from the analysis of the force of the ensemble in the deterministic framework are represented by the solid line, while the symbols shows the value of the parameter f 0 and r that can be estimated by taking into account the fluctuations of the force of the ensemble.

The result of the analysis for different system size N are reported in Supplementary Figure 10: the solid line represent f 0 as a function of N , while the symbols refers to the best fit value of f 0 as determined for different choices of N . We remark that the fitting procedure converges (with the requested limit of precision) only over a finite range of values of N , centred around the value adopted when performing the simulations. This observation implies in turn that we are in a position to obtain a reasonable estimate for the interval of pertinence of N , as follows the procedure outlined above.

Supplementary Table 2. Estimated parameters via the inverse scheme fed with simulated data. Errors are below 10 -3 if not explicitly provided.

|                      | F 0 (pN)   | f 0 (pN)   | r               | a (s - 1 )   | b (s - 1 )   | ϕ (s - 1 )   |
|----------------------|------------|------------|-----------------|--------------|--------------|--------------|
| True parameters      | 22 . 9     | 6 . 0      | 0 . 46          | 54 . 7       | 23 . 7       | 13 . 43      |
| Estimated parameters | 22 . 7     | 6 . 1      | 0 . 46 ± 0 . 03 | 54 . 8       | 23 . 5       | 13 . 42      |

Supplementary Figure 8. Results of the optimisation procedure for the parameters estimation.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-08fca6e6-813b-40ce-bf1d-7ef2269a124c.jpg)

The parameters estimation has been performed in the framework of the effective model and the results are identified by the yellow symbols (mean ± SD). Mean and SD for each parameter are computed from 58 independent iterations of the fitting procedure (black dots). The estimated values are compared with the parameters adopted in the stochastic simulations of the dynamics (blue solid line).

Supplementary Figure 9. Parameters estimates in the reference plane ( f , r 0 ) .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-80d7f531-890d-4ff0-a3a2-31b767430eb7.jpg)

The symbol follows from the integrated fitting strategy that accounts for fluctuations. Mean and SD is computed from 58 independent iterations of the fitting procedure (black dots). The solid line is the hyperbole populated with the degenerate mean field, hence deterministic solutions. Remark that the fitted symbol is close but not on top of the hyperbole. The observed deviation is eventually due to the residual population y ∗ that is adequately estimated via the generalised fitting strategy base on the stochastic description. The error is obtained from different replica of the stochastic optimisation algorithm.

Supplementary Figure 10. Checking f 0 against N .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-d034800c-36b6-408f-aa9b-b99770d75f04.jpg)

The parameter f 0 is estimated as a function of N , applying the inverse scheme to the simulated data.

## Parameters estimation for soleus and psoas HMM data

A representative example of the fitting outcome for the soleus HMM ensemble is reported in the Supplementary Figure 11. The histogram is generated from the experimental data series of the isometric force in the stationary state (i.e. at the isometric plateau). The results of the fitting procedure and the parameters estimation performed on the experimental data sets from rabbit psoas and rabbit soleus are listed in the Supplementary Table 3 and Supplementary Table 4 respectively. The result of the analysis for different system size N are reported in Supplementary Figure 10: the symbols refers to the best fit value of f 0 as determined for different choices of N . The shaded region identifies the portion of the parameters plane where the solutions are expected to be found. Specifically, it is assumed to lay in between the two curves:

$$f _ { 0 } N = \frac { 2 0 } { 1 1 } \left \langle F _ { 0 } \pm 2 \Delta F _ { 0 } \right \rangle \left \langle \frac { a } { b } \right \rangle$$

Here the relative error associated with the average value of the quantity ⟨ a/b ⟩ is assumed negligible, as compared to that stemming from the average stationary force, ∆ F 0 . The histogram computed from the collection of fitted parameters (each choice of symbols refers to a different experimental series) can be conceptualised as an indirect imprint of the degree of experimental variability as associated to f 0 and N .

Supplementary Table 3. Estimated parameters for psoas data

| Exp       | F 0 (pN)       | f 0 (pN)      | r               | ϕ (s - 1 )    |
|-----------|----------------|---------------|-----------------|---------------|
| PSO 1     | 16 . 4         | 5 . 6         | 0 . 36          | 6 . 4         |
| PSO 2     | 16 . 1         | 5 . 8         | 0 . 34          | 6 . 1         |
| PSO 3     | 20 . 7         | 7 . 6         | 0 . 32          | 6 . 2         |
| PSO 4     | 17 . 7         | 6 . 9         | 0 . 32          | 5 . 9         |
| PSO 5     | 19 . 7         | 8 . 0         | 0 . 30          | 5 . 8         |
| PSO 6     | 17 . 0         | 7 . 0         | 0 . 30          | 5 . 8         |
| mean ± SD | 17 . 9 ± 1 . 9 | 6 . 8 ± 1 . 0 | 0 . 32 ± 0 . 02 | 6 . 0 ± 0 . 2 |

Supplementary Table 4. Estimated parameters for soleus data

| Exp       | F 0 (pN)       | f 0 (pN)      | r               | ϕ (s - 1 )      |
|-----------|----------------|---------------|-----------------|-----------------|
| SOL 1     | 8 . 7          | 1 . 8         | 0 . 55          | 2 . 25          |
| SOL 2     | 10 . 9         | 2 . 5         | 0 . 51          | 2 . 27          |
| SOL 3     | 9 . 1          | 2 . 2         | 0 . 49          | 2 . 29          |
| SOL 4     | 13 . 4         | 3 . 2         | 0 . 49          | 2 . 26          |
| SOL 5     | 9 . 7          | 2 . 3         | 0 . 48          | 2 . 34          |
| SOL 6     | 9 . 3          | 2 . 4         | 0 . 48          | 2 . 23          |
| mean ± SD | 10 . 2 ± 1 . 7 | 2 . 4 ± 0 . 4 | 0 . 50 ± 0 . 03 | 2 . 27 ± 0 . 04 |

Supplementary Figure 11. Result of the fitting procedure on the force of the soleus HMMs ensemble.

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-fa06cb14-5150-4fba-b358-ba5c3bf60776.jpg)

The histogram of the soleus force at the isometric plateau is fitted against the analytical profile, via a self-consistent optimisation procedure which aims at estimating the unknown kinetic parameters.

Supplementary Figure 12. f 0 vs. N on the force of the soleus HMMs ensemble. The estimated parameter f 0 of the soleus HMMs is plotted as a function of the imposed N ; each choice of symbols refers to a different experimental series. The shaded region is drawn from the theoretical curve that resolve the dependence of f 0 on N .

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-9f8faed5-39f1-4e7e-8171-c28bb4bc8b50.jpg)

## Supplementary Note 2

## Rate of force development in situ

The time course of isometric force redevelopment following a release that drops the isometric force to 0 was determined in Ca 2+ -activated demembranated fibres from the same slow and fast muscles of the rabbit used to extract the myosin isoforms for the nanomachine.

As previously described linari2007stiffness , 26,32 , small bundles dissected from the two muscles were stored in skinning solution containing 50% glycerol at -20 ◦ C for 3-4 weeks and single fibres were prepared just before the experiment. A fibre segment 4-6 mm long was clamped at its extremities by T-clips and mounted between the lever arms of a loudspeaker motor and a capacitance force transducer lombardi1990contractile . To prevent sliding of the ends of the fibre segment inside the clips and minimise the shortening of the activated fibre against the damaged sarcomeres at the ends of the segment during force development, the extremities of the fibre were fixed first with a rigor solution containing glutaraldehyde and then glued to the clips with shellac dissolved in ethanol. Fibres were activated by temperature jump using a solution exchange system as previously described linari2007stiffness . A striation follower huxley1981millisecond allowed nanometre-microsecond resolution recording of length changes in a selected population of sarcomeres.

The composition of the solutions has been reported previously ( 26 , Supplementary Table 1, 25 ◦ C). The increase of interfilamentary distance following cell membrane permeabilisation was reversed by the addition of the osmotic agent Dextran T-500 (4% weight/volume).

The rate of force development was determined on the isometric force rise recorded after superimposing on the isometric contraction of the maximally Ca 2+ -activated fibre (pCa 4 5) a fast ramp shortening . (5-6% of the initial fibre length) able to drop the force to zero (Supplementary Figure 13).

![Image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00373v1-figure-d8c0dd51-c814-4e18-a0f6-d7bc182b8cd1.jpg)

## Supplementary Figure 13. Time course of force redevelopment in fast and slow skinned fibres.

Force redevelopment (lower traces) and corresponding half-sarcomere shortening (upper traces) after a period of unloaded shortening in a skinned fibre from rabbit psoas (black traces) and soleus (red traces) muscles. The vertical line indicates the time at which the force development starts. Force is normalised for the isometric value ( T 0 ) before the imposed large shortening ( ∼ 5% of the fibre length or ∼ 60 nm per hs). Further shortening against end compliance during force redevelopment was 29 ± 6nm per hs ( n = 4 ) and 19 ± 3nm per hs ( n = 11 ) in fast and slow fibres respectively. T 0 was 276 ± 44kPa and 195 ± 26kPa in fibres from psoas and soleus respectively. Temperature, 25 2 . ◦ C . 4% dextran T500 was added to reduce the lateral filament spacing of the relaxed fibre to the value before skinning.

## References

- 1. Hill, A. V. The heat of shortening and the dynamic constants of muscle. Proceedings of the Royal Society of London. Series B-Biological Sciences 126, 136-195 (1938).
- 2. Huxley, A. F. Muscle structure and theories of contraction. Prog. Biophys. Biophys. Chem 7, 255-318 (1957).
- 3. Fenn, W. O. The relation between the work performed and the energy liberated in muscular contraction. The Journal of physiology 58, 373 (1924).
- 4. Woledge, R. C., Curtin, N. A. &amp; Homsher, E. Energetic aspects of muscle contraction. Monographs of the physiological society 41, 1-357 (1985).
- 5. Smith, D. &amp; Geeves, M. Strain-dependent cross-bridge cycle for muscle. Biophysical Journal 69, 524-537 (1995).
- 6. Smith, D. &amp; Geeves, M. Strain-dependent cross-bridge cycle for muscle. II. Steady-state behavior. Biophysical journal 69, 538-552 (1995).
- 7. Smith, D., Geeves, M. A., Sleep, J &amp; Mijailovich, S. M. Towards a unified theory of muscle contraction. I: foundations. Annals of biomedical engineering 36, 1624-1640 (2008).
- 8. Linari, M. et al. Straightening out the elasticity of myosin cross-bridges. Biophysical Journal 118, 994-1002 (2020).
- 9. Schiaffino, S. &amp; Reggiani, C. Fiber types in mammalian skeletal muscles. Physiological reviews 91, 1447-1531 (2011).
- 10. Barclay, C., Constable, J. &amp; Gibbs, C. Energetics of fast-and slow-twitch muscles of the mouse. The Journal of physiology 472, 61-80 (1993).
- 11. Barclay, C., Woledge, R. &amp; Curtin, N. Is the efficiency of mammalian (mouse) skeletal muscle temperature dependent? The Journal of physiology 588, 3819-3831 (2010).
- 12. Wendt, I. &amp; Gibbs, C. Energy production of rat extensor digitorum longus muscle. American Journal of Physiology-Legacy Content 224, 1081-1086 (1973).
- 13. Wendt, I. &amp; Gibbs, C. Energy production of mammalian fast-and slow-twitch muscles during development. American Journal of Physiology-Legacy Content 226, 642-647 (1974).
- 14. Barclay, C. J., Woledge, R. C. &amp; Curtin, N. A. Inferring crossbridge properties from skeletal muscle energetics. Progress in Biophysics and Molecular Biology 102, 53-71 (2010).
- 15. Bottinelli, R, Betto, R, Schiaffino, S &amp; Reggiani, C. Unloaded shortening velocity and myosin heavy chain and alkali light chain isoform composition in rat skeletal muscle fibres. The Journal of physiology 478, 341-349 (1994).
- 16. Reggiani, C et al. Chemo-mechanical energy transduction in relation to myosin isoform composition in skeletal muscle fibres of the rat. The Journal of physiology 502, 449 (1997).
- 17. Potma, E., Van Graas, I. &amp; Stienen, G. Effects of pH on myofibrillar ATPase activity in fast and slow skeletal muscle fibers of the rabbit. Biophysical journal 67, 2404-2410 (1994).
- 18. Reggiani, C., Bottinelli, R. &amp; Stienen, G. J. Sarcomeric myosin isoforms: fine tuning of a molecular motor. Physiology 15, 26-33 (2000).
- 19. He, Z.-H., Bottinelli, R., Pellegrino, M. A., Ferenczi, M. A. &amp; Reggiani, C. ATP consumption and efficiency of human single muscle fibers with different myosin isoform composition. Biophysical journal 79, 945-961 (2000).
- 20. Barany, M, Barany, K, Reckard, T. t. &amp; Volpe, A. Myosin of fast and slow muscles of the rabbit. Archives of Biochemistry and Biophysics 109, 185-191 (1965).
- 21. Finer, J. T., Simmons, R. M. &amp; Spudich, J. A. Single myosin molecule mechanics: piconewton forces and nanometre steps. Nature 368, 113-119 (1994).
- 22. Ishijima, A. et al. Multiple-and single-molecule analysis of the actomyosin motor by nanometerpiconewton manipulation with a microneedle: unitary steps and forces. Biophysical journal 70, 383-400 (1996).

- 23. Pertici, I. et al. A myosin II nanomachine mimicking the striated muscle. Nature communications 9, 1-10 (2018).
- 24. Pertici, I., Bianchi, G., Bongini, L., Lombardi, V. &amp; Bianco, P. A myosin II-based nanomachine devised for the study of Ca2+-dependent mechanisms of muscle regulation. International journal of molecular sciences 21, 7372 (2020).
- 25. Pertici, I. et al. Muscle myosin performance measured with a synthetic nanomachine reveals a class-specific Ca2+-sensitivity of the frog myosin II isoform. The Journal of Physiology 599, 1815-1831 (2021).
- 26. Caremani, M. et al. The force of the myosin motor sets cooperativity in thin filament activation of skeletal muscles. Communications Biology 5, 1266 (2022).
- 27. Stienen, G., Kiers, J., Bottinelli, R &amp; Reggiani, C. Myofibrillar ATPase activity in skinned human skeletal muscle fibres: fibre type and temperature dependence. The Journal of physiology 493, 299-307 (1996).
- 28. Iorga, B., Adamek, N. &amp; Geeves, M. A. The slow skeletal muscle isoform of myosin shows kinetic features common to smooth and non-muscle myosins. Journal of Biological Chemistry 282, 3559-3570 (2007).
- 29. Huxley, A. F. &amp; Simmons, R. M. Proposed mechanism of force generation in striated muscle. Nature 233, 533-538 (1971).
- 30. Piazzesi, G. &amp; Lombardi, V. A cross-bridge model that is able to explain mechanical and energetic properties of shortening muscle. Biophysical journal 68, 1966-1979 (1995).
- 31. Caremani, M., Melli, L., Dolfi, M., Lombardi, V. &amp; Linari, M. Force and number of myosin motors during muscle shortening and the coupling with the release of the ATP hydrolysis products. The Journal of Physiology 593, 3313-3332 (2015).
- 32. Percario, V. et al. Mechanical parameters of the molecular motor myosin II determined in permeabilised fibres from slow and fast skeletal muscles of the rabbit. The Journal of Physiology 596, 1243-1257 (2018).
- 33. Zeng, W. et al. Dynamics of actomyosin interactions in relation to the cross-bridge cycle. Philosophical Transactions of the Royal Society B: Biological Sciences 359, 1843 (2004).
- 34. Pardee, J. D. &amp; Spudich, J. A. Mechanism of K+-induced actin assembly. The Journal of cell biology 93, 648-654 (1982).
- 35. Kron, S. J., Toyoshima, Y. Y., Uyeda, T. Q. &amp; Spudich, J. A. in Methods in enzymology 399-416 (Elsevier, 1991).
- 36. Suzuki, N., Miyata, H., Ishiwata, S. &amp; Kinosita Jr, K. Preparation of bead-tailed actin filaments: estimation of the torque produced by the sliding force in an in vitro motility assay. Biophysical journal 70, 401-408 (1996).
- 37. Uyeda, T. Q., Kron, S. J. &amp; Spudich, J. A. Myosin step size: estimation from slow sliding movement of actin over low densities of heavy meromyosin. Journal of molecular biology 214, 699-710 (1990).
- 38. Gillespie, D. T. A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. Journal of computational physics 22, 403-434 (1976).
- 39. Gillespie, D. T. Exact stochastic simulation of coupled chemical reactions. The journal of physical chemistry 81, 2340-2361 (1977).
- 40. Sadooghi-Alvandi, S., Nematollahi, A., Habibi, R, et al. On the distribution of the sum of independent uniform random variables. Statistical papers 50, 171-175 (2009).
- 41. Hall, P. The distribution of means for samples of size n drawn from a population in which the variate takes values between 0 and 1, all such values being equally probable. Biometrika, 240-245 (1927).
Valentina Buonfiglio Irene Pertici Matteo Marcello Ilaria Morotti Marco Caremani Massimo Reconditi Marco Linari Duccio Fanelli Vincenzo Lombardi Pasquale Bianco
Force and kinetics of fast and slow muscle myosin determined with a synthetic sarcomere-like nanomachine
Myosin II is the muscle molecular motor that works in two bipolar arrays in
each thick filament of the striated (skeletal and cardiac) muscle, converting
the chemical energy into steady force and shortening by cyclic ATP--driven
interactions with the nearby actin filaments. Different isoforms of the myosin
motor in the skeletal muscles account for the different functional requirements
of the slow muscles (primarily responsible for the posture) and fast muscles
(responsible for voluntary movements). To clarify the molecular basis of the
differences, here the isoform--dependent mechanokinetic parameters underpinning
the force of slow and fast muscles are defined with a unidimensional synthetic
nanomachine powered by pure myosin isoforms from either slow or fast rabbit
skeletal muscle. Data fitting with a stochastic model provides a
self--consistent estimate of all the mechanokinetic properties of the motor
ensemble including the motor force, the fraction of actin--attached motors and
the rate of transition through the attachment--detachment cycle. The
achievements in this paper set the stage for any future study on the emergent
mechanokinetic properties of an ensemble of myosin molecules either engineered
or purified from mutant animal models or human biopsies.