Multiple pathways of the actin-myosin cycle in energy transduction and the release of orthophosphate in muscle

In the striated muscle, the molecular motor myosin II functions in two bipolar arrays in each thick filament, converting chemical energy into steady force and shortening by cyclic ATP-driven interactions with nearby actin filaments. The fundamental steps in energy transduction are the working stroke, an inter-domain tilting of the lever arm about the actin-attached catalytic domain, generating up to ∼5 pN force or ∼10 nm of filament sliding, and the release of the ATP hydrolysis product orthophosphate (Pi) from the nucleotide-binding site, which is associated with a large free energy release. The two events are not simultaneous, as first demonstrated by the force response to a stepwise change in [Pi] (the Pi transient), showing the saturation kinetics characteristic of a two-step reaction. However, while high-resolution crystal structures of the myosin motor suggest that Pi release precedes the working stroke, in vitro functional studies indicate that it follows the working stroke. High-resolution sarcomere-level mechanics applied to single muscle fibers, allowing myosin motor synchronization by step perturbations in length or load, revealed that the kinetics of the working stroke is independent of [Pi] and depends only on the load. Moreover, this approach highlights the need for two unconventional pathways of the chemo-mechanical cycle: an early detachment of the force-generating motors and the possibility for attached motors to slip to the next actin monomer farther from the sarcomere center during shortening. Transient and steady-state responses to stepwise changes in load or [Pi] can be fitted with a structurally and biochemically explicit model in which the Pi release step is orthogonal to the progression of the working stroke. Model simulations indicate that the rate of Pi release depends on motor conformation, which resolves longstanding unanswered questions such as the dependence of Pi transient kinetics on the final level of [Pi] under any load and clarifies the issue of the relative timing between the working stroke and Pi release: at high loads, Pi release precedes the execution of the working stroke, while at low loads, the working stroke state transitions are fast enough to occur with Pi still bound to the catalytic site.

Introduction

Force and shortening are generated in the sarcomere, the structural unit of striated muscle, by cyclical, ATP-driven interactions between myosin motors, which extend in two bipolar arrays from the center of the thick filament, and the overlapping thin actin-containing filaments that originate from the sarcomere extremities (Figure 1).

Figure 1

Figure 1. Cartoon representing the sarcomere and its shortening induced by the working stroke of myosin motors. Reciprocal sliding between the actin filaments (white), originating from the Z-line at the sarcomere extremities, and the myosin filaments (black), bound to the M-line at the sarcomere center and carrying two bipolar arrays of myosin motors (represented by their motor domains or sub-fragment 1 S1, red; all actin-attached for simplicity), is powered by the working stroke of the actin-attached myosin motor (red to green oval). The working stroke consists of inter-domain tilting of the lever arm about the catalytic domain, which remains firmly attached to the actin site (gray).

In each interaction, the free energy of the hydrolysis of ATP in the catalytic site of the myosin motor is converted into mechanical work through a structural change (the working stroke) that, according to the crystallographic model (Rayment et al., 1993a; Rayment et al., 1993b), consists of an inter-domain tilting of the lever arm about the catalytic domain firmly attached to the actin site. The working stroke is associated with the sequential release of the hydrolysis products orthophosphate (Pi) and ADP. Calorimetric and protein solution studies have shown that a large enthalpy change is related to Pi release from the quaternary complex actin–myosin-ADP-Pi (AMADPPi), suggesting that the Pi dissociation step is the chemical transition associated with the mechanical working stroke (White and Taylor, 1976; Hibberd et al., 1985; Hibberd and Trentham, 1986).

Photo-liberation of caged Pi in skinned muscle fibers showed that the force during an isometric contraction is reduced by a stepwise increase in [Pi] with a rate that increases asymptotically with [Pi] (Dantzig et al., 1992; Millar and Homsher, 1992; Walker et al., 1992; Homsher et al., 1997), implying that, according to the minimum five-step reaction of Scheme 1, the Pi release (step 4) follows the force-generating transition from a weak AMADPPi to a strong AM'ADPPi state (step 3).

Scheme 1

Scheme 1. 5 step reaction scheme from Dantzig et al., 1992.

In Reaction Scheme 1, which summarizes the results from solution protein studies and caged compound experiments in skinned fibers, ATP binding to the actomyosin (AM) complex, after the release of the hydrolysis products and the execution of the working stroke, promotes rapid dissociation of myosin from actin (step 1), which is followed by the recovery stroke (the reversal of the working stroke) and ATP hydrolysis (step 2). The detached myosin with the hydrolysis products is in rapid equilibrium with the weakly bound AMADPPi state. The closure of the actin-binding cleft of the myosin catalytic domain (Geeves and Holmes, 2005) forms the strongly attached, stiffness-generating AM'ADPPi state, which is allosterically associated with the structural change that leads to tilting of the lever arm responsible for force generation and filament sliding (step 3). Thus, according to the simplified Scheme 1, step 3 is the combination of two processes, namely, the formation of a strongly bound motor and the generation of force, with kinetics limited by the much slower attachment process. This definition of step 3 is in agreement with the fast kinetics of force generation, as estimated from the quick force recovery following a step release (Huxley and Simmons, 1971; Ford et al., 1977), and it is supported by the evidence that the force increases in proportion to the number of myosin motors during isometric force development by a muscle fiber (Brunello et al., 2006; Caremani et al., 2008). Pi is released without any further contribution to force (step 4, Dantzig et al., 1992). Thus, the sum of the fractional occupancies of the AM'ADPPi and AM'ADP states constitutes the fraction (f) of actin-attached motors that contribute equally to the force and stiffness. The subsequent ADP release (step 5) occurs at a rate that, in the isometric contraction, is low (Nyitrai and Geeves, 2004; Sleep et al., 2005; West et al., 2005). According to Scheme 1, the development of isometric force is rate-limited by both the ATP hydrolysis step (step 2) and the attachment of motors (step 3), while the steady-state flux through the whole cycle (the rate of ATP hydrolysis) is limited by the rate of the ADP release (step 5). ADP release is conformation-dependent and becomes relatively fast only when attached motors preferentially populate the final state of the working stroke, such as during shortening at low load (Nyitrai and Geeves, 2004; Caremani et al., 2015). Thus, in vitro and in situ kinetic studies support Scheme 1 as the unique path through which myosin–actin interaction in muscle hydrolyzes ATP and produces force and power.

This conclusion, however, is questioned on the basis of several pieces of challenging evidence, which are discussed in detail in the next sections. First, most crystallographic models of the myosin motor domain (sub-fragment 1, S1) (Smith and Rayment, 1996; Houdusse et al., 2000; Geeves and Holmes, 2005; Llinas et al., 2015) support the alternative view that Pi is released prior to the working stroke (see also Robert-Paganin et al., 2020; Moretto et al., 2022; Mansson et al., 2023; Rassier and Mansson, 2025). Second, the finding that an increase in [Pi] reduces the isometric force of the Ca²⁺-activated fibers more than the ATPase rate (Cooke and Pate, 1985; Caremani et al., 2008) suggests that an increase in [Pi] increases the probability of an alternative pathway of the chemo-mechanical cycle (Linari et al., 2010; Debold et al., 2013), with premature detachment of motors from a strongly bound force-generating state, followed by the release of the hydrolysis products. Third, the attachment–detachment kinetics that fits the rate of isometric force generation is too slow in relation to the maximum power developed during steady shortening, which suggests that during shortening, an attached motor could complete its ATPase cycle by slipping to the next actin monomer farther from the sarcomere center (Lombardi et al., 1992; Piazzesi and Lombardi, 1995; Caremani et al., 2013; Pertici et al., 2023).

Relation between the release of Pi and the working stroke

Crystallographic studies (Smith and Rayment, 1996; Houdusse et al., 2000; Geeves and Holmes, 2005; Llinas et al., 2015) have the unique ability to provide details of the inter-domain movements within the motor domain (S1, Figure 1) that associate the closure of the actin-binding cleft, which is responsible for the formation of the strongly-bound motor, to the catalytic site and then to the converter that controls the orientation of the lever arm (Figure 2). At the catalytic site, modifications in the orientation of the relevant structural elements (P-loop, switch 1, and switch 2) that join the subdomains promote both a decrease in affinity for orthophosphate and a chain of changes in neighboring residues, culminating in a tilt in the converter.

Figure 2

Figure 2. Schematic drawing of the subdomains within the myosin catalytic domain and their reciprocal structural changes accompanying the working stroke transition. The four subdomains identified by different colors are the N-terminal (black), the upper 50-kDa subdomain (blue), the lower 50-kDa subdomain (gray), and the converter (green). Two (yellow and red) of the structural elements joining the subdomains are visible. The α-helix constituting the lever arm is cyan. The structures are obtained from the scallop myosin: (A) S1 complexed with VO₄^3-, representing the pre-working stroke state. (B) Nucleotide-free S1, representing the end-of-stroke state. The vertical arrows indicate the approximate direction of the actin filament axis. The orientation of the two structures is defined by conserving the orientation of the lower 50-kDa relative to actin. The release of Pi is indicated with the dashed arrow. The exit direction is an open question, and more recent evidence (Llinas et al., 2015) indicates that it is opposite (backdoor) to what is shown here. The direction of the movement of the subdomains in the transition between the two states is shown by the short continuous arrows. The ∼70° change in the orientation of the converter is turned into ∼11 nm axial movement by the length of the lever arm (modified from Figure 2B by Houdusse et al. (2000). Copyright (2000) National Academy of Sciences, United States).

In contrast to this view, Förster resonance energy transfer (FRET) which, unlike crystallography, allows for the direct recording of structural dynamics within the myosin motor, indicates that the working stroke precedes Pi release (Muretta et al., 2015). Similar conclusions have been reached with single-molecule mechanics using a three-bead assay system (Takagi et al., 2004), which were later confirmed with an enhanced time resolution obtained by eliminating the effect of the trap compliance by imposing resistive loads from the very beginning of the interaction (Capitanio et al., 2012; Woody et al., 2019). A critical issue in all the above in vitro mechanical studies is that the loading condition under which the myosin–actin interaction occurs is far from reproducing the physiological load range, while the phenomena under investigation, the working stroke and the release of Pi, have individual unknown load sensitivities as far as their extent and kinetics are concerned. In particular, it must be considered that FRET measurements were achieved in unloaded conditions, and thus, the fractional occupancy of attached motors derived from this study is biased to the state at the end of the working stroke. On the other hand single-molecule measurements showed a working stroke rate that increases and a lifetime of actin-attached motors that decreases with the increase of the resisting load, opposite to the relations expected from the bulk of in situ mechanical and energetic data (Hill, 1938; Huxley and Simmons, 1971; Reconditi et al., 2004; Pertici et al., 2023).

In our experience, the most reliable approach for clarifying the issue of the relative timing and load-dependence of the working stroke and the Pi release is the sarcomere-level mechanics with nanometer-microsecond resolution achieved by integrating fast force transducers (Huxley and Lombardi, 1980) with the striation follower (Huxley et al., 1981), which was first developed for intact fibers isolated from frog skeletal muscle (Lombardi and Piazzesi, 1990) and then applied to single skinned fibers of the mammalian skeletal muscle (Linari et al., 2007; Caremani et al., 2013) (Supplementary Figure 1A). The striation follower is an optoelectronic apparatus that enables the measurement and control of half-sarcomere (hs) length changes in a selected population of sarcomeres, with sub-nanometer precision and 2-μs time resolution. As shown in Figure 3A and Supplementary Figure 1B, a stepwise reduction of half-sarcomere length by a few nanometers (complete in ∼100 μs) superimposed on the isometric contraction elicits a force transient composed of a force change simultaneous with the length step (phase 1) to a value T₁, which depends on elastic properties of the half-sarcomere, followed by a force recovery due to the active properties of the myosin motor. The earliest component is a rapid force recovery (phase 2), which represents the mechanical manifestation of the working stroke of attached motors synchronized by the step release. For releases larger than 6 nm per hs, the force recovered in phase 2 is only a fraction of the force decrease during the elastic response and is followed by a pause in force recovery (phase 3) that precedes the final recovery to the original isometric force T₀ (phase 4). Phase 3 results from the synchronous attainment of the end of the working stroke by attached motors and their accelerated detachment, whereas phase 4 is explained by the isometric attachment–detachment kinetics. Phase 2 recovery holds the information used for developing a model of the working stroke (Huxley and Simmons, 1971; Ford et al., 1977): the maximum extent of filament sliding for which attached motors can maintain force (estimated by the abscissa intercept of the T₂ relation) was ∼11 nm (Figure 3C), which is in agreement with the crystallographic model of the working stroke proposed 20 years later (Rayment et al., 1993a). The speed of phase 2 force recovery, r₂, measured as the reciprocal of the time to recover from T₁ to T₁ + 0.63*(T₂-T₁), is larger for release than for stretch and increases with the release size (Figure 3D). This dependence of r₂ on the step size and direction excludes the response of a passive viscoelasticity and holds the constraints for the kinetic model of force generation (Huxley and Simmons, 1971).

Figure 3

Figure 3. Force and velocity transients. (A) Early components of force transients (middle trace in each panel) elicited by a step release (upper trace) superimposed on the steady force of the isometric contraction (T₀) of a single fiber from frog skeletal muscle. Release amplitude per half-sarcomere (hs), 2.3 nm (upper panel) and 6.5 nm (lower panel). The lower horizontal line is the baseline for force. (B) Isotonic velocity transients following steps in force superimposed on the isometric force (T₀) of a single fiber from frog skeletal muscle. Upper trace, stepwise decrease to a force T = 0.5 T₀; middle trace, force baseline; lower traces, change in half-sarcomere length (L) corresponding to steps to the forces indicated by the figures on the right. The figures on the left of the 0.1 T₀ trace show the phases of the transient. Reproduced from Figure 1 of Piazzesi et al. (2002) with permission. (C) Comparison of L₁ (filled circles) and L₂ (filled diamonds) relations from velocity transients with T₁ (open circles) and T₂ (open diamonds) relations from force transients. Reproduced from Figure 4 of Piazzesi et al. (2002) with permission. (D) Rate of early force recovery (r₂) versus step amplitude. The line is obtained by fitting the r₂ data points with an exponential. Reproduced from Figure 3C of Piazzesi and Lombardi (1995). (E) Relation between the initial shortening velocity of phase 2 (V₂) and force (T/T₀). The line is an exponential fit to the data. Reproduced from Figure 3A of Piazzesi et al. (2002).

Figure 4

Figure 4. Relevant parameters of the kinetic model of the working stroke. (A) Free energy diagram of the myosin motor states (M₁–M₅) as a function of the relative position (x) between the myosin motor and the actin monomer. x is set to 0 for the minimum free energy of an attached motor in the state M₁. The free energy minima of the various states are separated by a distance z (the step size) = 2.75 nm so that the working stroke of 11 nm is preserved with a number of structural transitions = 4. The thick line marks the axial distribution of motors during the isometric contraction. (B) Functions expressing the x-dependence of the forward (continuous lines) and backward (dashed lines) transitions. The figures indicate the order of the transition. (C) Fractional occupancy of the various states as a function of x during isometric contraction. (D) Sketch of the conformation of the five states (M₁–M₅) of the myosin motor characterized by orientations of the lever arm separated by angles corresponding to an axial displacement of 2.75 nm of the tip of the lever under zero force. The lower panel shows the fractional occupancy of the different states in isometric contraction. (E) Simulated velocity transients following force decreases to the values indicated by the figures close to the traces. (F) Filled circles and continuous line: simulated V₂-T/T₀ relation; open circles: experimental data from Figure 3E. (G) Filled symbols and lines: simulated relations between r₂ and the size of the step release in the absence of filament compliance (triangles) and in the presence of 1.7 nm/T₀ (squares), 2 nm/T₀ (filled circles), and 2.3 nm/T₀ (diamonds) filament compliance. Open circles are experimental data from Figure 3D. Panels A–C, reproduced from Figure A1 of Piazzesi et al. (2014); panel D, reproduced from Figure 4 of Fusi et al. (2014); panels E and F, reproduced from Figure 4 of Piazzesi et al. (2014); panel G, reproduced from Figure 5D of Piazzesi et al. (2014).

However, the presence of a significant compliance in the myofilaments, functionally in series with that of the motor array in each half-sarcomere (Huxley et al., 1994; Wakabayashi et al., 1994), complicates the use of the rate of quick force recovery r₂ as the constraint for modeling the working stroke kinetics (Linari et al., 2009; Piazzesi et al., 2014). The rationale is illustrated with a simplified mechanical model of the half-sarcomere (Supplementary Figure 2), in which the distributed compliance of actin and myosin filaments is represented by an elastic element in series with the array of myosin motors. Under these conditions, force recovery occurs while motors move in the shortening direction to increase stress in the filaments, and its rate r₂ decreases in proportion to filament compliance. The availability of a capacitance force transducer with a resonant frequency of 50 kHz (Huxley and Lombardi, 1980) and, thus, an adequately fast force clamp enabled an alternative protocol in which the velocity transient was recorded following stepwise decreases in force superimposed on the isometric steady force (Figure 3B; Piazzesi et al., 2002). The elastic shortening in phase 1 (L₁) simultaneous with the force step is followed by a rapid phase-2 shortening (due to the synchronization of the working stroke in the attached motors) with an extent (L₂ - L₁, where L₂ is the hs shortening attained at the end of phase 2) that is larger and attained earlier at lower force. L₂ (the sum of the elastic strain and the working stroke) for a decrease in force to zero coincides with the abscissa intercept of the T₂ relation (Figure 3C). Phase-2 shortening is followed by a pause (phase 3, due to synchronized detachment of motors) and then a steady shortening (phase 4) at a lower velocity that is characteristic of the detachment–attachment kinetics underpinning the force–velocity relation of the contracting muscle. The phases that follow the elastic (phase 1) response occur under force clamp and, thus, are not affected by the compliance of myofilaments. Under this condition, the dependence on force of the speed V₂ estimated by the tangent to the initial part of phase-2 shortening (Figure 3B) directly measures the load dependence of the speed of the working stroke (Figure 3E), providing the constraints for a mechanical–kinetic model of the myosin working stroke in situ (Figures 4A–C).

For simplicity, the model considers only the mechanical and kinetic properties of the attached myosin motors that determine phase 2 of the responses to length and force steps; thus, the kinetic scheme is based on Huxley–Simmons 1971 model (Huxley and Simmons, 1971) and does not include further steps of the actin–myosin interaction, such as the attachment and detachment of motors. In this study, we provide a synthetic description of the model, which is described in detail by Piazzesi et al. (2014). The stiffness of the myosin motor is 3 pN/nm, and this sets the slope of the parabolas representing the free energy profile of attached states and, thus, the size (2.75 nm) and number (4) of the force-generating steps necessary to fit the transition and equilibrium kinetics of the working stroke. The rate functions of the transitions between the five structural states marking the progression of the working stroke (Figure 4D) were selected by the model simulation to achieve the best fit of the phase-2 velocity transient and its dependence on load (Figures 4E, F). Then, the same kinetic scheme was applied to simulate the relation between the rate of phase-2 force transient (r₂) and the size of the step release (Figure 4G, filled triangles). The simulated relation was shifted upward with respect to the observed relation (open circles), with the simulated r₂ already four times higher for 1.5 nm release. The difference between the observed and simulated relations was the consequence of neglecting the depressant effect of filament compliance on the rate of force generation by the motor array in the model. Assuming that the filament compliance was the only adjustable parameter, it was found that the same kinetic scheme can fit the observed relation when an equivalent filament compliance of ∼2 nm/T₀ per hs is assumed (filled circles), a value in fairly good agreement with the filament compliance estimated with mechanical and X-ray diffraction experiments on the same preparation (Brunello et al., 2014).

The simplified kinetic scheme does not imply detachment/attachment; thus, the simulation of the velocity transients ends with the equilibrium distribution of attached states at each force, which underpins working stroke amplitude that reduces with the load, mostly due to the reduced contribution of the recoil of elastic elements (Figure 4E). The reduction in the observed working stroke amplitude with increase in the load is larger (Figure 3B) because the completion of the working stroke is progressively truncated by the ensuing synchronized detachment in phase 3.

These conclusions were straightforwardly confirmed in their structural counterpart, achieved by using the same mechanical protocol (Figure 5A) in combination with time-resolved X-ray diffraction from synchrotron light (Supplementary Figure 1) (Reconditi et al., 2004).

Figure 5

Figure 5. Motion of myosin motors during the velocity transient. (A) Force change normalized by isometric force T₀ and length change in nm/hs. (B) Axial X-ray intensity distribution in the region of the M₃ reflection recorded at the BioCAT beamline of the APS Synchrotron (Argonne, United States). Colors denote X-ray exposure periods, as shown in (A) brown, L₀; orange, L_2s; magenta, L_2e; and cyan, L₃. (C) Structural organization of myosin motors in the muscle sarcomere at each X-ray exposure, as predicted by the M3 interference fine structure in (B). (A–C) From Figure 1 of Reconditi et al. (2004). (D) Superimposed half-sarcomere shortenings (left panel) following force decreases to three different values (identified by the color code in the inset) and intensity ratio (R_M3, right panel) of the higher over the lower angle component of the M3 intensity distribution. From Figure 2A, B of Reconditi et al. (2004).

X-ray diffraction by itself lacks the phase information required to directly define the change in shape of the diffracting unit, but the size and direction of the movement can be recovered with sub-nanometer precision by exploiting the X-ray interference between the two bipolar arrays of motors in each sarcomere. Following the decrease in force, the changes in the interference fine structure of the so-called M3 reflection originating from the 14.5-nm axial periodicity of myosin motors (Figure 5B) indicate the movement of attached motors toward the center of the sarcomere in phase 2, coherent with the execution of the working stroke, and their movement away from the center at the onset of phase 3, associated with motor detachment at the end of the working stroke followed by re-attachment farther from the sarcomere center (Figure 5C). The structural data (Figure 5D, right panel) match the striation follower signal (left panel), indicating a larger and faster working stroke with decreasing load.

A reliable method for in situ recording of the working stroke is the prerequisite for investigating its sensitivity to [Pi]. By applying the striation follower technology and the fast force clamp protocols to skinned fibers from rabbit psoas (Caremani et al., 2013), it has been shown that increase in [Pi] to 10 mM (i) reduces the isometric force by 35%, without a significant effect on the unloaded shortening velocity (V₀) (see also Cooke et al., 1988), (ii) does not affect the rate constant of the working stroke, as measured by the speed of phase-2 shortening, and (iii) shortens phase-3 pause and accelerates the subsequent transition to the final steady shortening velocity (Figure 6).

Figure 6

Figure 6. Effect of [P_i] on the isotonic velocity transient. (A) Superimposed isotonic velocity transients following a force step to 0.5 T₀ in control (black) and 10 mM P_i (green). (B) Superimposed phase 2 of the velocity transients in control (black) and in 10 mM P_i (green), which were obtained after subtracting phase-1 and phase-3 shortenings. (C) Relation between the shortening accounted for by the working stroke (L_T) and the force expressed in relative units T/T₀ in control (black) and 10 mM P_i (green). (D) Dependence of the rate constant of the working stroke (r₂) on T/T₀ in control (black) and 10 mM P_i (green). (E) Dependence of the rate of phase 3 (r₃) on T/T₀ in control (black) and 10 mM P_i (green). Reproduced from Figures 2, 3 of Caremani et al. (2013) with permission.

Most of the literature, in contrast to the above conclusion, reports a direct effect of Pi on the force generation process. This finding can be explained considering that those studies lack the resolution to discriminate between motor attachment and force generation (for instance, the process 2πb identified with sinusoidal analysis by Kawai and Halvorson, 1991). A study on the effect of [Pi] on the rate of the working stroke, which was determined in the same preparation using the quick force recovery elicited by length steps (Ranatunga et al., 2002), produced contradictory results: increasing [Pi] did not affect the rate of quick force recovery following a step release, whereas it increased the rate following a step stretch, mainly at the expense of its latter component. These results, however, can be justified by taking into account that in the quick force recovery following a step stretch, there is a contamination of the reversal of the working stroke by rapid detachment–attachment. With double-step protocols in single frog muscle fibers, it has been demonstrated that (i) the quick force recovery from a step release up to 6 nm per hs is complete within 2 ms and is almost fully accounted for by the synchronous execution of the working stroke because the detachment–attachment responsible for the rapid repriming of the working stroke is three times slower (Lombardi et al., 1992); (ii) the recovery from a stretch of 4 nm (a size comparable to that of Ranatunga et al., 2002), instead, takes ∼20 ms, and only 1/3 of it is explained by the synchronous reversal of the working stroke while the remaining part is due to the ensuing detachment–attachment that, in the force transient following a stretch, merges with the reversal within phase 2 (Piazzesi et al., 1997). Thus, the sensitivity to [Pi] of the later part of the recovery from a stretch is explained by the Pi-dependent acceleration of detachment in phase 3 of the velocity transient (Caremani et al., 2013).

In summary, the analysis of the velocity transients shows that the rate of the working stroke is independent of [Pi] and solely depends on the load and that the increase in [Pi] accelerates the phase-3 detachment–attachment (Caremani et al., 2013). These conclusions constrain the modeling of the chemically and structurally explicit cycle implemented from Reaction Scheme 1 (Reaction Scheme 4 presented later) to simulate both transient (phases 2 and 3) and steady-state (phase 4) shortening responses following the decrease in force below T₀ (Caremani et al., 2013), which is in agreement with the original idea of Huxley and Simmons (1971) that the progression through the working stroke is controlled by rate functions that solely depend on the strain under which the transition from a lower to a higher force-generating state occurs. At the same time, the release of Pi from the catalytic site of the myosin motor can occur at any stage of the structural transitions. In this way, the release of Pi is an orthogonal process to the working stroke, and its rate does not directly depend on the motor strain, which is compelling for a chemical step. At the same time, as detailed in the section describing the properties of Reaction Scheme 4, the rate of Pi release can be made conformation-dependent, in which case it increases as the motor progresses through the working stroke.

Evidence for premature motor detachment from a strongly bound force-generating state induced by an increase in [Pi]

Addition of Pi to the solution induces a decrease in the Ca²⁺-activated steady isometric force (T₀) of a skinned fiber (Figure 7A) (Cooke and Pate, 1985; Kawai and Halvorson, 1991; Dantzig et al., 1992; Caremani et al., 2008) and, to a lesser extent, a decrease in the isometric ATPase rate (Figure 7B) (Webb et al., 1986; Bowater and Sleep, 1988; Cooke et al., 1988; Potma et al., 1995; Potma and Stienen, 1996). In 10 mM Pi, T₀ is reduced to ½ the value in control solution (with 0.7 mM–1 mM of contaminating Pi Pate and Cooke, 1989a; Millar and Homsher, 1990), while the ATPase rate is reduced by less than 20%.

Figure 7

Figure 7. Effects of [Pi] on the isometric force, ATPase rate, and rate of force change during initial force increase and following a Pi jump. (A) Pi-dependence of isometric force (relative to the force in the control solution, 1 mM Pi). Filled circles, data from Figures 1F, 6A of Caremani et al. (2008). Model predictions: blue line, Scheme 1; red line, Scheme 2. (B) Pi-dependence of the ATPase rate: filled circles, pooled data from Bowater and Sleep (1988), Potma et al. (1995) and Potma and Stienen (1996). Model predictions: blue line, Scheme 1; red line, Scheme 2. Values relative to those in the control solution. (C) Pi-dependences of k_TR (the rate of force redevelopment from a large release bringing the force to 0) and k_Pi (the rate of force decrease following a Pi jump). Filled circles, observed k_TR relation from Figure 6C of Caremani et al. (2008). Model prediction: blue continuous line, Scheme 1; red continuous line, Scheme 2. Open squares, observed k_Pi relation from Figure 6 of Dantzig et al. (1992). Model prediction: blue dashed line, Scheme 1; red dashed line, Scheme 2. Modified from Linari et al. (2010).

In terms of a conventional chemo-mechanical cycle in which Pi release occurs with the working stroke (Pate and Cooke, 1989a; Pate and Cooke, 1989b), the apparent contradiction was explained by hypothesizing that an increase in [Pi] shifts attached motors toward a state generating lower force without any effects on the ATP turnover. However, this idea was contradicted by evidence that the Pi-dependent reduction in force is accompanied by a proportional reduction in the number of attached motors (Caremani et al., 2008), a finding anticipated by work using sinusoidal analysis (Kawai and Halvorson, 1991), even though the conclusion in that case appeared questionable because of both the limited time resolution and the absence of analysis of the effects of filament compliance. A Pi-dependent reduction in the number of attached motors in isometric contraction without a comparable reduction in the ATP turnover indicates that the flux of myosin motors through the attachment/force-generating step (step 3 in Reaction Scheme 1) is larger than the flux through the conventional cycle estimated through steps 5 and 1, demonstrating that the conventional model cannot explain the reduced effect of [Pi] on the ATPase rate with respect to that on the force. The problem is solved by assuming a branched pathway that allows the myosin in the AM'ADPPi force-generating state to detach from actin before the release of ADP and Pi and then quickly release the hydrolysis products through an unconventional cycle (see Reaction Scheme 2; Linari et al., 2010). This conclusion was later supported by the results and modeling of the effect of Pi on the in vitro mechanics of a mini-ensemble of myosin motors (Debold et al., 2013). Notably, to be effective, the release of hydrolysis products from the AM'ADPPi state must be assumed irreversible; so, for the mass action, an increase in Pi promotes only the flux through the AMADPPi state.

Scheme 2

Scheme 2. 5 step reaction scheme integrated by a branched pathway accounting for the detachment from a strong bound AM'ADPPi state followed by release of hydrolysis products.

Notably, the branched pathway introduced in Reaction Scheme 2, while enhancing the ATP consumption at any [Pi] fitting the observed relation (red line in Figure 7B), maintains the property of Reaction Scheme 1. It fits the dependence on [Pi] of either the rate of isometric force development (k_TR, Figure 7C; filled symbols, experiment; continuous blue line, Reaction Scheme 1; continuous red line Reaction Scheme 2) or the rate of force reduction elicited by a stepwise increase in [Pi] produced by photo-liberation of caged Pi during the isometric contraction (k_Pi, Figure 7C; open symbols experiment; dashed blue line Reaction Scheme 1; dashed red line Reaction Scheme 2) equally well. This is because both the schemes assume that force generation and Pi release occur as a two-step process, with the first step (force generation, corresponding to step 3 in Reaction Scheme 1) and the second step (the release of Pi, corresponding to step 4 in Reaction Scheme 1) being sufficiently fast (>500/s) to be considered in rapid equilibrium, with a K ∼10 mM (Dantzig et al., 1992; Homsher et al., 1997).

Pros and cons of Reaction Scheme 1 from Pi-transient experiments

The force decrease in response to an increase in [Pi] (Pi transient) imposed on Ca²⁺-activated fibers exhibits an amplitude A_Pi and a rate k_Pi that provide the fundamental constraints of the kinetics of the underlying reactions. Dantzig et al. (1992) found that the observed k_Pi can be restricted to a section of the whole cycle involving the two-step reaction:

$k_{Pi}=k_a+k_b*\frac{\left[\text{Pi}\right]}{\left(K_c+\left[Pi\right]\right)},(1)$

where k_a and k_b are the forward and backward rate constants of step 3 in Reaction Scheme 1 and K_c is the equilibrium constant of Pi release. The process 2πb of sinusoidal analysis (Kawai and Halvorson, 1991) follows the same two-step reaction kinetics. Notably, the rate of isometric force development following a large release that reduces the force to 0, k_TR, which involves the whole attachment–detachment cycle, shows a Pi-dependence with saturation kinetics (Figure 7C, filled circles) that holds the footprint of the two-step reaction of Equation 1.

To test the strain-dependence of the rate constants of step 3, the analysis of the Pi transient was extended to steady shortening at different velocities, thus reducing the load on the attached myosin motors in this way (Figure 8; Homsher et al., 1997). It was shown that A_Pi reduces with shortening velocity (Figure 8A), while k_Pi increases in proportion with the shortening velocity (Figure 8B) and, at a relatively slow shortening velocity, shows the same hyperbolic dependence on [Pi] as in the isometric contraction (Figure 8C). The increase in k_Pi and the reduction in A_Pi with increasing shortening velocity could be explained by introducing a strain-dependence of the rate constants of step 3 in the Reaction Scheme 1, as proposed by Dantzig et al. (1992). This model, however, was not able to fit the dependence of k_Pi on [Pi] during shortening, even at the low velocity used in Figure 8C (filled circles). As demonstrated with the model simulation reported in the final section, the explanation for the failure, which is common to all models with Pi release in series with the working stroke (Homsher et al., 1997; Ranatunga, 1999; Smith, 2014; Stehle and Tesi, 2017; Offer and Ranatunga, 2020), can be found in the assumption that Pi leaves the active site from a unique motor conformation either before or after the working stroke and, thus, with a unique value of its rate constant. Instead, a model where the working stroke and the Pi release are orthogonal processes (Caremani et al., 2013; Caremani et al., 2015) can fit the data in Figure 8C as it allows the Pi release rate to increase with the change in motor conformation during its progression through the working stroke.

Figure 8

Figure 8. Pi transients in isometric contraction and during isovelocity shortening. (A) Force responses to Pi jump (points) imposed in isometric contraction and during shortening at the velocities (ML/s means fiber length/s) indicated on the traces. Records were interpolated by the sum of an exponential and a straight line (continuous line). (B) Relation of k_Pi versus shortening velocity. (C) Relation between k_Pi and final Pi concentration in isometric contraction (open circles) and during shortening at 0.1 fiber length/s (filled circles). Dashed lines interpolated to data. Continuous lines in (B–C) are the predictions of the model reported in Scheme 1 (reproduced from Figures 3, 4, 8 in Homsher et al. (1997) with permission).

Matching the mechanics and energetics of shortening muscle requires a myosin motor to interact with two actin sites during one ATP hydrolysis cycle

The first evidence of a loose coupling between the mechanical and biochemical cycle of actin–myosin interactions during shortening was the finding that the working stroke on the attached myosin motors could be reprimed much faster than expected from the rate of the ATP hydrolysis (Lombardi et al., 1992). With double-step releases, each <5 nm per half sarcomere, superimposed on the isometric contraction (Figure 9), it was demonstrated that, if delivered within 1–2 ms following the first release, the second step release elicits a quick recovery that corresponds to that elicited by the sum of the two releases. Meanwhile, the quick recovery from a second step release imposed progressively later, increases and recovers the value expected from a single-step release with a time constant of ∼6 ms (at 4 °C), which is at least four times faster than expected from the kinetics of the attachment step, underpinning the rate of isometric force development from 0 (time constant ∼26 ms; Figure 3A, B in Linari et al. (2015)).

Figure 9

Figure 9. Force transient elicited by a step release imposed at different times following a conditioning step. (A) Superimposed records of force responses (lower traces) to (a) a single step release of 7 nm per half-sarcomere (upper traces) and to (b) a step of 2 nm delivered 2 ms after a conditioning step of 5 nm. (B) Superimposed records of force responses to test steps of 2 nm delivered 2, 4, 8, and 15 ms after the conditioning step of 5 nm. Single muscle fiber from Rana esculenta tetanically stimulated at 2.1 μm sarcomere length and 4 °C. Half-sarcomere length changes measured in a selected fiber segment measuring ∼1 mm in length. From Figure 1 in Lombardi et al. (1992).

Structural evidence of the fast repriming of the working stroke was obtained shortly afterward using X-ray diffraction (Irving et al., 1992), exploiting changes in the intensity of the third-order myosin-based meridional reflection to track the changes in the orientation of the lever arm of the myosin motor. In isometric contraction, the motor lever is oriented nearly perpendicular to the filament axis. During the quick force recovery after the step release, it tilts toward the center of the sarcomere, marking the synchronous execution of the working stroke, and then, within the next 15 ms, it recovers the isometric conformation in parallel with the recovery of the ability to generate a second working stroke, which is much faster than expected from the rate of the isometric force development.

The problem that the rate of motor attachment responsible for the transition from 0 to T₀ in isometric conditions is too low to account for the rate of the working stroke regeneration during shortening is further exacerbated by the finding that the kinetics of the attachment step underpinning the development of isometric force is too low to account for the maximum power developed during isotonic contractions against loads of approximately 1/3 T₀ (Pertici et al., 2023). On the other hand, the increase in the rate constant for the motor attachment required to fit the maximum power would imply an increase in the rate of force development to values larger than those observed. A way to keep the rate of isometric force development as low as required by that observed, preserving the power of the shortening contraction, is to assume a specific geometrical hindrance in the isometric condition consequent to the mismatch between actin and myosin periodicities (Marcucci et al., 2021). In this case, however, the fit of the observed power produced during shortening implies that the ATP hydrolysis rate (or the rate of energy liberation) increases by 8–9 times with respect to the isometric rate. This is a major drawback of the model of Marcucci et al. (2021) because the observed rate of energy liberation during an isotonic contraction developing maximum power is much lower, ranging from 2 (mammalian muscle) to 4 (frog muscle) times larger the isometric energy rate (Pertici et al., 2023, and references therein).

As shown by the model of Marcucci et al. (2021), any model that attempts to simultaneously simulate the limited rate of isometric force development and the power during isotonic shortening by increasing the apparent rate of attachment during shortening under the assumption that more motors become available with shortening (Chen and Brenner, 1993) finds its limit in an exaggerated increase in the rate of energy consumption for the maximum power relative to the isometric rate. This same limitation also applies to the Huxley and Tideswell (1997) model, in which the first attached motor of a myosin dimer during shortening facilitates the attachment of the partner motor. Eventually, the hypothesis that the reduced rate of isometric force development can be explained by mechanosensing-based thick filament activation (Linari et al., 2015), masking the effect of a rate of attachment that is high enough to account for the power developed during isotonic contraction, is contradicted on the energetic basis. In fact, in this case, the underlying ATP hydrolysis rate during steady isometric contraction would be fivefold higher than that observed (Marcucci and Reggiani, 2016).

The Reaction Scheme 3 (Figure 10) (Pertici et al., 2018) has the characteristics of a model that reconciles both the rapid regeneration of the working stroke determined in the transient regime (Lombardi et al., 1992) and the maximum power determined in the steady-state regime with the constraints that (i) the rate of isometric force development is four times slower than the rate of the working stroke regeneration following a step release and (ii) the rate of energy liberation in the isotonic contraction at the maximum power is no more than four times larger than that in the isometric contraction.

Figure 10

Figure 10. Simplified model with myosin slipping (Reaction Scheme 3). Kinetic scheme with three states of the myosin motor (blue): D, detached state; A1 and A2, low- and high-force states, respectively, attached to an actin monomer (brown). During shortening, the motor attached in the A2 state can slip to the next actin monomer farther from the center of the sarcomere (red) within the same ATPase cycle. The probability of a second slipping to the pink monomer is limited to 1/10 of that of the first slipping. From Figure 4 of Pertici et al. (2018).

The demonstration, summarized in Figure 11 and Table 1, is given for the performance of the frog skeletal muscle at 4 °C–5 °C (Pertici et al., 2023), but, as shown in the same paper, it is equally valid for the fast mammalian skeletal muscle.

Figure 11

Figure 11. Model simulation of the performance of the half-thick filament of the frog skeletal muscle at 4.6 °C. (A) Force–velocity (F–V) relations. The force per half-thick filament (F, pN) is calculated from the force per cross-sectional area (T, with T₀ = 165 kN m^-2) taking a density of filament of 5.87 10¹⁴ m^-2 (Piazzesi et al., 2018). Black circles and lines, experimental F–V relation; green circles and line, simulation with Reaction Scheme 3a, which is the integral version of Reaction Scheme 3 in Figure 10; violet circles and line, simulation with the suppression of slip transition (Reaction Scheme 3b); brown circles and line, simulation with the suppression of slip transition and ad hoc increase in the rate constants of attachment/detachment (Reaction Scheme 3c). (B) Power–force (P–F) relations calculated from data in A (same color code as in A). (C) Dependence of the rate of ATP hydrolysis per motor (φ) on V, according to the simulation with the model identified by the same color code as in (A). Arrows indicate the values of φ at V for P_max (φ_Pmax). In Reaction Scheme 3c, without slipping as in Reaction Scheme 3b, the increase in the attachment/detachment rate constants able to predict the experimental P_max implies an increase in φ_Pmax to ∼10 times the isometric value (the ordinate intercept). (D) Redevelopment of isometric force following a 5% rapid shortening able to decrease the isometric force F₀ to 0. The rise time (from 10% to 90% of full change) is indicated next to the record. Black, observed (from Figure 3A green in Linari et al., 2015); simulated rises are identified by the color code as in (A) green (Reaction Scheme 3a), violet (Reaction Scheme 3b), and brown (Reaction Scheme 3c). Modified from Pertici et al. (2023).

Table 1

Table 1. Simulations of the relevant mechanical and energetic parameters of the half-thick filament with Reaction Scheme 3a (integral version of the model shown in Figure 10), Reaction Scheme 3b (slipping suppressed), and Reaction Scheme 3c (slipping surrogated by ad hoc increase of the attachment–detachment kinetics). F₀, steady isometric force per half-thick filament; φ₀, steady isometric ATPase rate per myosin motor; V₀, unloaded shortening velocity; P_max, maximum power; φ_Pmax, ATPase rate at P_max; t_r, rise time of isometric force development. Reaction Scheme 3a implies the possibility of the attached motor in state A2 to slip to the next Z-ward actin site. Reaction Scheme 3b retains the same kinetic scheme as Reaction Scheme 3a but excludes the possibility of slipping. In Reaction Scheme 3c, the attachment and detachment rate constants are increased to fit the observed F–V relation and P–F relation as shown in Figure 10. Data are presented as the mean ± SEM from at least 12 simulations. Table from Pertici et al. (2023).

The three-state mechano-kinetic model of the actin–myosin interaction in Figure 10 has already been described in detail (Pertici et al., 2018; Pertici et al., 2020). The way in which the most relevant energetic features of the model are constrained by literature data is summarized here. In isometric contraction, the rate-limiting step in the cycle is detachment from the high-force generating state A2: the rate of ATP splitting per myosin motor (φ) is the minimum under the isometric conditions (φ₀). During steady shortening, φ increases as the rate of motor detachment increases due to the accelerated execution of the working stroke. φ for the maximum power (φ_Pmax) is higher than φ₀ by a factor of four (Barclay et al., 2010). Under this condition, the curvature of the F–V relation and the resulting maximum power can be reproduced only by assuming the integral version of Reaction Scheme 3 (3a), in which, during shortening, the attached myosin motors can rapidly regenerate the working stroke during the same ATPase cycle by slipping to the next actin monomer farther from the center of the sarcomere (Lombardi et al., 1992; Piazzesi and Lombardi, 1995) (step “slip” in Figure 10) and undergoing A1′–A2′ state equilibration according to the strain dependency of step-2 kinetics (Huxley and Simmons, 1971). Detachment from either A1′ or A2’ (step 3′) implies ATP hydrolysis.

All the relevant emergent properties of the array of motors in the half-thick filament (htf) of fast skeletal muscle are reproduced by Reaction Scheme 3a. In Figure 11, the F–V relation (A), where F is the force per htf, is shown, and the P–F relations (B), simulated with Reaction Scheme 3a (green lines), are superimposed on the experimental data (black lines). The corresponding simulated φ–V relation, where φ is calculated by the flux through step 1, is shown in Figure 11C (green line). The relevant parameters are shown in the first line of Table 1. Notably, the observed P_max is simulated with φ (φ_Pmax), which is four times larger than the value for the isometric force F₀ (φ₀).

The limits of the conventional mechano-chemical cycle in which the ATP hydrolysis is completed within a single actin–myosin interaction are summarized in the simulations of Reaction Scheme 3b and c, where the slipping transition is removed. If no other adjustments are introduced (Reaction Scheme 3b), F₀ and φ₀ are not substantially affected, while P_max is reduced to less than ½ (violet line in Figure 11B) without a marked change in φ_Pmax (violet arrows in Figure 11C; Table 1). With Reaction Scheme 3c, the mechanical performance of fast skeletal muscle (brown lines in Figures 11A, B) is recovered in the absence of slipping with ad hoc increases in the relevant rate constants for attachment and detachment (for the details of the kinetic assumptions, see Pertici et al., 2023). The consequences are that while φ₀ can be maintained in substantial agreement with the experimental value, φ_Pmax (brown arrow in Figure 11C; Table 1) increases by ∼10-fold with respect to φ₀.

The mechanical parameter that also in isometric conditions is sensitive to the kinetic differences between Reaction Scheme 3a and c is the rate at which force reaches the steady-state isometric value F₀. This is expected from the consideration that the time course of force increase depends on the rate constants governing attachment (step 1) and force generation (step 2) and, specifically, in this model, on step 1, which is the rate limiting step in the process. To avoid the influence of the time required for thick filament activation on initial force development (Linari et al., 2015), the selected reference parameter is the force redevelopment following a 5% rapid shortening, which can reduce and maintain the isometric force at 0 (Figure 11D, black trace). The rise time t_r (the time from 10% to 90% of F₀), estimated from the exponential fit to the recorded trace, is 52 ms, and it is perfectly reproduced by the simulation with Reaction Scheme 3a (Figure 11D, green trace, and first row in Table 1). In contrast, t_r is reduced to less than ½ in the simulation with Reaction Scheme 3c (Figure 11D, brown trace, and third row in Table 1), demonstrating that Reaction Scheme 3c fails to predict the observed low rate of the attachment-force development. Notably, t_r simulated by Reaction Scheme 3b (Figure 11D, violet trace, and second row in Table 1) fits the observed value. This is because Reaction Scheme 3b shares the same kinetic scheme as Reaction Scheme 3a, apart from the slipping transition, which is effective only for the performance of the shortening muscle.

Chemically and structurally explicit model where Pi release is orthogonal to the working stroke transition

Reaction Scheme 4 (Figure 12) (Caremani et al., 2013; Caremani et al., 2015) integrates the properties of Reaction Schemes 1–3, described in the previous sections, into a chemically and structurally explicit model that retains the predictive properties of the previous schemes and adds the possibility to both simulate the velocity transients elicited by step perturbations in force and their modulation by increase in [Pi] and identify the biochemical states that are involved in unconventional pathways such as early detachment from a force-generating state or, during shortening, slippage to the next actin monomer away from the sarcomere center.

Figure 12

Figure 12. Structurally and chemically explicit cycle of the myosin motor integrating the properties of the three preceding schemes (Reaction Scheme 4). Transitions shown in black indicate the cycle undergone by myosin motors (M) that interact with only one actin monomer (A). Transitions shown in red indicate the cycle undergone by myosin motors slipping to the next actin monomer 5.5 nm away from the center of the sarcomere (A′). For simplicity, this legend is mostly limited to black transitions. M₁–M₄ represent the structural states defining the progression of the working stroke in a given biochemical state of the attached myosin motor. Binding of an ATP molecule induces dissociation of the myosin motor from actin (step 1, 1′), which is followed by the recovery of the conformation of the motor at the beginning of the working stroke and the hydrolysis step (step 2). The MADPPi–AMADPPi state represents both the detached motor just after the hydrolysis step and recovery stroke and the weakly actin-bound motor with the hydrolysis products still bound. Strong binding of M to an actin monomer A (step 3) implies the formation of the first of four different force-generating states, AM₁ADPPi, which, without significant delay, undergoes the structural transition, leading to the strained conformation responsible for the isometric force. The working stroke, in any of the biochemical states, AMADPPi, AMADP, and AM, implies three subsequent force-generating steps (M₁→ M₂, M₂→ M₃, and M₃→ M₄) controlled by strain-dependent rate constants (Huxley and Simmons, 1971). Both biochemical events in the attached motor, release of Pi (step 4, 4′) and release of ADP (step 5, 5′), can occur in any of the four structural states. Unconventional pathways are represented by the possibility of the force-generating AMADPPi state to detach (step 6, 6′) and rapidly release the hydrolysis products and rebind ATP (step 7, 7′) and the possibility that during shortening, the motor attached to the first actin monomer (A) slips to the next Z-ward actin monomer (A′, steps 8 and 9). This possibility is specifically relevant to the AMADP state (step 9; see the rate functions listed in Supplementary Table 1). Reproduced from Figure 1 in Caremani et al. (2015) with permission.

The following is a summary of the Reaction Scheme 4. The two-step reaction responsible for force generation and release of Pi in Reaction Scheme 1 (steps 3 and 4, respectively), predicting the hyperbolic dependence of k_Pi (but also of k_TR) on [Pi] (Figure 7C), is reproduced in Reaction Scheme 4 (Figure 12) by steps 3 and 4 (and 4′). The early detachment from the force-generating AM'ADPPi state and the subsequent rapid release of the hydrolysis products in Reaction Scheme 2 (steps 6 and 7, respectively), which are required to explain the reduced effect of the increase of [Pi] on the ATPase rate (Figure 7B), are reproduced in Reaction Scheme 4 by steps 6 (and 6′) and 7 (and 7′), respectively. The slipping of an actin-attached motor during shortening to the next actin monomer farther from the sarcomere centre within the same ATP hydrolysis cycle (step “slip” in Reaction Scheme 3a, Figure 11), required for a rapid regeneration of the working stroke (Figure 9) and for the observed maximum power (Figure 11B) with a limited increase in the rate of energy liberation (Figure 11C) during steady shortening, is reproduced in Reaction Scheme 4 by steps 8 (slipping from the AMADPPi states) and 9 (slipping from the AMADP states).

The biochemically different force-generating states of Reaction Scheme 1 that are significantly populated in physiological conditions (AMADPPi and AMADP) in Reaction Scheme 4 exist in several conformations depending on the progression through the working stroke. Strong binding of M to an actin monomer A (step 3), caused by the closure of the actin-binding cleft (Geeves and Holmes, 2005), implies the formation of the first, AM₁ADPP_i, of several states generating progressively higher force. The kinetic, mechanical, and energetic features constraining the simulation of the working stroke have already been described in detail for the intact frog muscle fiber (Figure 4). In that case, the stiffness of the motor was 3 pN/nm so that the size of the force-generating transition, corresponding to the progressive swinging of the lever arm, was set to 2.75 nm (Piazzesi et al., 2014), and the number of transitions for an 11-nm working stroke was four (Figure 4D). In the case of the mammalian-skinned fiber, the smaller stiffness of the motor (≤1.7 pN/nm; Linari et al., 2007) allowed the assumption of a larger size of the transition (3.1 nm) and only three force-generating transitions, namely, M₁→ M₂, M₂→ M₃, and M₃→ M₄. According to the principle of the nearest-neighbor interaction, attachment to the actin monomer (with approximately 5.5-nm axial separation between monomers along a strand of the actin helix) occurs for a range of x from −2.75 nm to 2.75 nm, where x is the relative axial position between the motor and A, and is 0 for the position of the center of distribution of attachments of the motors in the M₁ state (Figure 13A). AM₁ADPPi undergoes the first transitions leading to the strained conformation responsible for the isometric force (Figure 13A).

Figure 13

Figure 13. x-distribution of the attached states in isometric contraction and during shortening and simulation of the force–velocity and power–force relations with Reaction Scheme 4. (A-C) x-distribution of the attached states in isometric contraction and during shortening. Dashed lines, AMADPPi states; continuous lines, AMADP states. Progression of states through the working stroke identified by the progression from the thinnest to the thickest line, as indicated in the inset. (A) Distribution at T₀. R_m, the ratio between AMADP and AMADPPi motors, is ∼7. (B). Distribution during isotonic shortening at V = 1,000 nm/s per hs, R_m is ∼2. (C). Distribution at V₀, R_m is ∼1. (D) T–V relations. Circles, experimental data; lines, simulated relations with Reaction Scheme 4 with (green) and without (violet) steps 8 and 9. (E) P–T relations calculated from T–V relations in (D). Circles and lines, same codes as in (D).

Both the release of P_i (step 4) and the release of ADP (step 5) can occur from any of the four attached states, that is P_i release and and ADP release are orthogonal to the working stroke state transitions. The most important consequences for the aims of this report are as follows: first, the working stroke kinetics is independent of the concentration of Pi, as shown from the data in Figure 6D; second, the kinetics of Pi release can be made to depend on the conformation of the motor, increasing with the progression of the motor through the working stroke from M₁ to M₄. Notably, this kinetic feature is excluded in all the reaction schemes in which Pi-release coincides (Pate and Cooke, 1989b) or is in series (Dantzig et al., 1992; Homsher et al., 1997; Ranatunga, 1999; Smith, 2014; Stehle and Tesi, 2017; Offer and Ranatunga, 2020) with the force-generating transition. The rate functions of the reaction scheme are given in Supplementary Table 1. Further details on the calculation procedure can be found in Caremani et al. (2013) and Caremani et al. (2015).

Reaction Scheme 4 (Figure 12) retains all the properties of the three preceding schemes, and most of the mechanical and energetic outputs, observed for mammalian skinned fibers both without and with added P_i, are reported by Caremani et al. (2013). Therefore, its predictive power is demonstrated in detail here, mainly for the transient mechanical properties not simulated in the previous sections. In addition, with Reaction Scheme 4, as already demonstrated with the simplified Reaction Scheme 3 (Figures 11A, B), suppressing the possibility for an attached motor to slip to the next actin monomer farther from sarcomere centre during shortening increases the curvature of the T–V relation (Figure 13D, violet line) and reduces the maximum power (E, violet line) with respect to the experimental data (circles).

Again, as shown in the Reaction Scheme 4, the rate of transition from 0 to the steady-state isometric force is not affected by the suppression of the slipping possibility (Figure 14A; black trace, experimental record; green and violet lines, force rise simulated with Reaction Scheme 4 with and without steps 8–9, respectively).

Figure 14

Figure 14. Simulation of transient responses and the effect of Pi in skinned fibers from rabbit psoas with Reaction Scheme 4. (A) Superimposed time course of transition from 0 to T₀ in the isometric contraction. Black trace, experiment; green trace, simulation with the complete scheme; violet trace, simulation after suppression of slipping possibility (steps 8 and 9). (B) Superimposed velocity transients elicited by a force drop to 0.5 T₀. Black trace, experiment; green and violet traces, simulation with and without slipping, respectively. (C) Dependence on T/T₀ of the rate constant of the working stroke (r₂) elicited by a stepwise decrease in force. Black circles, experimental data from Figure 5D without (filled) and with 10 mM added Pi (open); green circles, simulation without (filled) and with 10 mM added P_i (open, masked by the filled circles). (D) Simulation of the effect of 10 mM Pi on the velocity transient elicited by a force drop to 0.5 T₀ without (green) and with (orange) 10 mM added Pi. (E) Dependence on [Pi] of k_TR, the rate of the transition from 0 to T₀ in the isometric contraction. Circles, experimental k_TR–Pi relation from Figure 7C; line, simulated relation. (F) Dependence on the shortening velocity of the ratio of the rate of ATP hydrolysis per motor with 10 mM added Pi over that without (R). (G) Relation between the rate of the Pi transient, k_Pi, and the shortening velocity. Dots, experimental data from Figure 8B; line, model simulation. (H) Relations between k_Pi and final [Pi] in isometric contraction (filled circles) and during shortening at 0.1 fiber length/s (open circles) from Figure 8C. Continuous and dashed lines, simulations fitted to the data in isometric condition and during shortening, respectively. Green, simulations with Reaction Scheme 4; orange, simulation with the same reaction scheme as green but with the assumption of the same value (100 s^-1) for the rate constant of the Pi release (k₄) for the four conformations.

Reaction Scheme 4 is unique in allowing simulation of the force and velocity transients elicited in response to a step in length or load that synchronizes the execution of the working stroke and its dependence on mechanical and biochemical conditions. As already discussed, the presence of a significant compliance in the myofilaments, functionally in series with the motor array in each half-sarcomere, complicates the use of the rate of quick force recovery as a constraint for modeling the working stroke kinetics (Piazzesi et al., 2014). Instead, phase-2 shortening elicited by a stepwise decrease in force (Figures 3B, 5D) directly provides information on the load dependence of the rate and the size of the working stroke, and the following phase-3 pause provides information on the kinetics of the events leading to working stroke regeneration. Reaction Scheme 4 can satisfactorily fit the isotonic velocity transient (Figure 14B: black, experiment; green, model), providing a direct demonstration of the role of slipping in the rapid regeneration of the working stroke, which is the condition for the maximization of the power during the more physiological steady shortening. Following the removal of steps 8–9, the phase-3 duration increases (Figure 14B, violet), delaying the transition to the steady shortening velocity of phase 4, with consequences on the steady-state (physiological) performance, such as the increase in the curvature of the force–velocity relation and the reduction in the power output (violet lines in Figures 13D, E). As already discussed for the simplified Reaction Scheme 3, this drawback, intrinsic to any conventional chemo-mechanical cycle, cannot be remedied by increasing the rate constant for the attachment step, because of the limit imposed to the kinetics of this step by the observed rate of isometric force development (Figures 11D, 14A).

A consequence of the idea that the release of Pi is orthogonal to the working stroke transitions is that only the kinetics of the working stroke is strain-dependent. Instead, Pi-release kinetics is sensitive to the structural state of the attached motor (conformation-dependent) and, thus, indirectly connected to the mechanical conditions by the speed of the progression through the different conformations (M₁–M₄). Under these conditions, the model satisfies the constraints that the speed of the working stroke (the rate of phase-2 rapid shortening, r₂) elicited by a stepwise force decrease (i) increases with the size of the step (Figure 14C: filled symbols: black, experiment; green, model) and (ii) is not affected by the increase in [Pi] (Figure 14C; open symbols: black, experiment; green, model).

With smaller steps, i.e., at higher loads, the working stroke is slower (Figures 3B, 5D), which increases the probability that the other reactions taking over in the phase-3 pause (early detachment with the hydrolysis products in the catalytic site (step 6) and slipping to the next actin monomer (steps 8 and 9)) truncate the working stroke and lead to its regeneration. An increase in [Pi] alone also truncates the working stroke and accelerates phase-3 pause for the same relative decrease in force (Figures 6A–C, E). Notably, this effect is accompanied by an increase in both k_Pi and k_TR (Figure 7C), revealing that the involved kinetic steps are not restricted to the unconventional early detachment but also include the conventional cycle.

Reaction Scheme 4 can discriminate two different mechanisms leading to the accelerated working stroke regeneration during shortening: (i) at low Pi, the operating mechanism is based on the possibility that, for critical negative values of x, the force-generating motors, mainly those in the AMADP state, slip to the next Z-ward actin monomer (step 9), reducing phase-3 pause (Figure 14B, green) and, in this way, increasing the sliding distance per one ATPase cycle (see Figure 7A in Caremani et al., 2013), without any effect on k_TR (Figure 14A); (ii) the operating mechanism induced by the increase in Pi (Figure 14D: green, no added Pi; orange, 10 mM Pi) is based on the acceleration of the whole attachment-force generation and detachment cycle, as revealed by a corresponding increase in k_TR (Figure 14E; circles: experiment, line: model simulation), as a consequence of the increase in the second-order rate constant of the reversal of Pi release. Under these conditions, the model predicts that, in isometric contraction and at high loads, an increase in [Pi] increases the rate of ATP hydrolysis per motor (Figure 14F), which is the tension cost.

Conventional models in which attachment-force generation is strain-dependent can predict the finding from Pi jump experiments that the rate of the Pi transient is larger during isovelocity shortening than under isometric conditions (Figure 8A) and increases with an increase in shortening velocity (Figure 8B). However, those models cannot predict that during isovelocity shortening, the rate of Pi transient increases with an increase in final [Pi], as in the isometric contraction (Figure 8C). This is because in all those models, the Pi release is a step in series with the working stroke and thus can be defined by only one rate constant. In Reaction Scheme 4, Pi release is orthogonal to the working stroke, and its kinetics depends on the effects that the progression of the working stroke exerts on the catalytic site. Assuming that the rate of Pi release (k₄) increases with the transition of motor conformation from M₁ to M₄ (Figure 15, see also Supplementary Material) is the sole condition that allows the rate of Pi transient to depend on either the velocity in isovelocity contractions (Figure 14G) or the final [Pi] both in isometric contraction (Figure 14H, green continuous line) and during shortening (green dashed line), as observed (circles, filled during isometric contraction, and open during shortening). If k₄ is given a convenient constant value (100 s^-1), the simulation still predicts k_Pi to depend on the final [Pi] in isometric contraction (orange continuous line) but not during shortening (orange dashed line).

Figure 15

Figure 15. Strain-dependence of the rate constants of the working stroke transitions and Pi release. Forward transitions, continuous lines; backward transitions, dashed lines. Black, working stroke transitions; line thickness increases with the progression of the transition from the first (M₁–M_2, the thinnest line) to the third (M₃–M₄, the thickest line). Magenta, Pi release step; line thickness increases with the progression of motor conformation from M₁ (the thinnest line) to M₄ (the thickest line). For simplicity, the transitions considered are those concerning the myosin motor while attached to the first actin monomer A (and the slip to the next actin monomer A′, 5.5 nm away from the sarcomere center can be deduced by shifting the rate functions for A leftward by 5.5 nm).

Contribution of Reaction Scheme 4 in clarifying unsolved questions of the coupling between the myosin working stroke and Pi release

Kinetic aspects

First, the structurally and biochemically explicit Reaction Scheme 4 allows clarifying the longstanding issue of the relative timing between the myosin working stroke and Pi release in fast skeletal muscle. According to this scheme, the release of Pi can occur at any stage of the working stroke with a conformation-dependent kinetics, implying that the rate of Pi release increases with the transition of motors from the M₁ to M₄ states (Figure 15, magenta lines). Reaction Scheme 4, first published in 2013 (Caremani et al., 2013), resolves all contradictions inherent to conventional models in which Pi release occurs in series with the working stroke, either before (Smith, 2014; Offer and Ranatunga, 2020) or after the stroke (Kawai and Halvorson, 1991; Dantzig et al., 1992; Homsher et al., 1997). Both hypotheses have the limitation that they do not explain (i) why the increase in [Pi] decreases the number of attached/force generating motors much more than the ATPase rate and (ii) why Pi transient kinetics depends on [Pi] both in the isometric contraction and during shortening. Moreover, the hypothesis that Pi is released before the working stroke has the further contradiction that the working stroke kinetics depends on the concentration of Pi (see Supplementary Figure 5 in Moretto et al., 2022).

The powerful tool offered by a structurally explicit model is the possibility to test the consequences on the muscle performance of the idea that the Pi release rate is conformation-dependent. This exquisitely emerges from the analysis of the fractional occupancies of the two biochemically distinct attached states in relation to the load of the contraction (Figure 13): at T₀, the distribution of the attached motors is around x = 0 and is biased toward the early structural states (mostly M₁ and M₂, Figure 13A; the thinnest and second thinnest lines) by the strain-dependent kinetics of the working stroke transitions (Figure 15, black lines; see Caremani et al., 2015, while the Pi release kinetics (magenta lines), although relatively slow, favors the AMADP state). The ratio between the AMADP and AMADPPi states (R_m) is 6.8. During isotonic shortening at V = 1 μm/s per hs (close to that for the maximum power in the skinned rabbit fiber) (Figure 13B), the motor distribution spreads toward negative values of x, favoring the transition toward structural states ahead in the working stroke (M₃ and M₄). Because of the conformation-dependent increase in the rate of Pi release, the result is a progressive shift from the AMADPPi state to the AMADP state of the M₃ and M₄ motors. At the same time, the Pi release becomes the rate-limiting step for the short-lived M₁ and M₂ motors. The result is that the fractional occupancy of the AM₁ADPPi state (the thinnest dashed line) and the AM₂ADPPi state (the second thinnest dashed line) is higher than the corresponding AM₁ADP (the thinnest continuous line) and AM₂ADP (the second thinnest continuous line) states, respectively, while the proportion is inverted for the M₃ and M₄ motors (the second thickest and thickest lines, respectively). Altogether, R_m decreases to 2.1. R_m further reduces to 1.1 for the maximum shortening velocity (V₀ = 3.5 μm/s per hs, Figure 13C). This is because V₀ implies a further shift of the distribution toward negative x and, thus, even faster strain-dependent working stroke transitions. Consequently, Pi release becomes the rate-limiting step for all motor conformations, which are mostly in the AMADPPi state, except M₄, for which the fractional occupancy of the AMADPPi state (the thickest dashed line) progressively decreases for more negative x with respect to that of the AMADP state (the thickest continuous line). This is because M₄, representing the final working stroke conformation, becomes a progressively more long-lived state as it attains more negative x. At V₀, the distribution of attached motors extends well beyond x = −10 nm, and an M₄ motor requires more than 3 ms to reach values beyond −10 nm, a time three times longer than the time constant for the release of Pi.

The inhibitory effect of the increase in [Pi] on the number of force-generating motors is much higher than that on the ATPase rate and is explained in Reaction Scheme 4 with the presence of an unconventional short pathway for the completion of the ATPase cycle, consisting of the early detachment of the motors from the force-generating AMADPPi states (step 6 in Figure 12), followed by rapid and irreversible release of the hydrolysis products and binding of a new ATP (step 7). Under this condition, the flux through step 6 is intrinsically sensitive to any increase in [Pi] through the mass action exerted by the second-order rate constant of the reversal of step 4. This pathway reflects its effects on the Pi transient kinetics, weakening the interpretation of Pi transients with the two-step reaction expressed by Equation 1.

The role of the two-step reaction is further undermined by considering the refined study of the Pi transient kinetics, which was made possible by the use of rabbit psoas myofibrils (Tesi et al., 2000). Unlike skinned fibers, in which a stepwise increase in [Pi] can be imposed by photo-liberation of caged Pi, in myofibrils, stepwise changes in [Pi] in either direction can be achieved by a rapid exchange system made by a double-barreled pipette. It was found that both the rate of the decrease in force elicited by a stepwise increase in [Pi] (k_+Pi) and its hyperbolic dependence on the final [Pi] were similar to those observed with Pi jumps in skinned fibers. Instead, the rate of the force rise elicited by a stepwise reduction in [Pi] (k_-Pi) was 2–3 times lower than that of the force reduction elicited by a stepwise increase to the same final [Pi]. In addition, the relation expressing k_-Pi dependence on the final [Pi] was shifted downward with respect to the k_+Pi-Pi relation, becoming similar to the k_TR–Pi relation characterizing the force redevelopment from 0 following a release superimposed on the isometric force. Two questions emerge from these results, both challenging the solidity of the interpretation of the Pi transient with a two-step reaction scheme. The first concerns the asymmetry between k_+Pi and k_-Pi, which is unexpected on the basis of the two-step reaction; the second concerns the finding that k_-Pi exhibits the same kinetics as k_TR and, thus, overcomes the limits of the two-step reaction, encompassing the whole attachment–detachment cycle. An explanation for the finding that k_+Pi > k_-Pi (and k_TR) was that, in analogy with the rapid phase of force relaxation following the decrease in activating Ca²⁺, the decrease in force in the isometric condition, whatever the reason is, it is accompanied by the development of sarcomere inhomogeneity and formation of weak sarcomeres that speed-up force decrease (Stehle, 2017). However, this argument was contradicted by evidence from Tesi et al. (2000), showing that the force transient following Ca²⁺ or ADP jumps in either direction did not show any direction-dependent kinetics. Moreover, while k_-Pi and k_TR were sensitive to the concentration of activating Ca²⁺, as expected for a reaction scheme implying de novo motor attachment (Brenner, 1988), k_+Pi was not. These results (i) exclude the possibility to explain the Pi transient following an increase in [Pi] simply with an increase in the flux through the reversal of the attachment step (step 3 in Figure 12)—a possibility that is also excluded by the evidence that an increase in [Pi] decreases the force and the number of attached motors much more than the ATP turnover (Webb et al., 1986; Bowater and Sleep, 1988; Cooke et al., 1988; Potma et al., 1995; Potma and Stienen, 1996; Caremani et al., 2008)—and (ii) suggest that the increase in [Pi] promotes the flux through a detachment path different from the reversal of attachment and committed for ATP hydrolysis, as the alternative cycle through step 6. If this unconventional pathway is irreversible, the increase in [Pi] by mass action would promote the increase in the flux through both the conventional (reversal of step 3) and the unconventional pathways (steps 6 and 7), leading to a faster Pi transient. Moreover, the finding that k_-Pi is smaller than k_+Pi and similar to k_TR could be explained according to Reaction Scheme 4 by considering that a reduction in [Pi] would promote only de novo attachment through the conventional pathway. According to the ¹⁸O exchange experiments, the Pi release step is reversible in skinned fibers but not in actomyosin S1 in solution (Bowater and Sleep, 1988). The same irreversibility could characterize the release of Pi from the unconstrained MADPPi state generated by the detachment of the force generating AMADPPi motors (Figure 12). In conclusion, for fast skeletal muscle, any alternative model (Stehle, 2017; Mansson, 2025) based on the assumptions that k_Pi+, k_Pi-, and k_TR are similar and the asymmetry in the Pi transient arises from an increase in sarcomere inhomogeneity with force decrease is not sustained by experimental evidence. Instead, the idea that the asymmetry of the Pi transient kinetics in fast skeletal muscle is because of an early detachment of force-generating motors with the hydrolysis products, introduced to explain the limited effect of the increase in [Pi] on the rate of ATP hydrolysis (Linari et al., 2010), is further supported by the behavior of these parameters in slow (soleus) skeletal muscle. It was found that at temperatures ≤15 °C, the increase in [Pi] has a similar inhibitory power of isometric force and ATPase rate (Potma et al., 1995; Kerrick and Xu, 2004). In this case, one expects that the probability of a flow through the unconventional cycle remains low independent of Pi. Accordingly, in the same preparation, it was found that there is no asymmetry in the Pi transient kinetics, and both k_+Pi and k_-Pi are similar to k_TR (Millar and Homsher, 1992; Stehle and Tesi, 2017).

Structural aspects

Crystallographic, EM, and transient biochemical studies of the decade across the end of the last century (Rayment et al., 1993a; Zeng et al., 2004; Geeves and Holmes, 2005) provided the classical synthetic description of the relevant molecular aspects of actin attachment and force generation mentioned in the Introduction section (Figure 2). The transition of the interactions between myosin S1 and actin from weak to strong, with the closure of the cleft between the upper and lower 50 kDa domains, is the trigger for the change in a joining structure (switch 1) that induces both the release of Pi from the active site and the rotation of the converter-tilting of the lever arm (Figure 2). This scenario has been challenged by a more recent crystallographic study on non-muscle myosin VI integrated with site-specific mutations and kinetic measurements (consisting of rapid freezing of the crystal soaked in high [Pi] for different times) (Llinas et al., 2015). It is shown that Pi leaves the active site before the complete closure of the cleft (recorded by quenching of actin pyrene label; De La Cruz et al., 1999) because the early changes promoted by the formation of the actin–myosin interface trigger changes of another switch (switch 2) that opens the way for Pi to leave the active site from the backdoor. Thus, an intermediate state (Pi release state, PiR) precedes the cleft closure, switch 1 movement, and the execution of the working stroke. The new scenario includes the presence of possible interactions of Pi within the tunnel that leads to the backdoor, which would delay the appearance of Pi in solution with respect to the time it leaves the active site. Therefore, even if the working stroke is promoted by Pi leaving the active site, Pi may appear in solution after the end of the working stroke, especially at a low load, which explains the results of the FRET experiment (Muretta et al., 2015).

How do the mechanical and kinetic features of Reaction Scheme 4 integrate with this new structural scenario linking Pi release and the working stroke? The first question concerns the identification of the functional properties of the PiR state. A hypothesis promoted by the new crystallographic data (Llinas et al., 2015; Robert-Paganin et al., 2020; Moretto et al., 2022; Mansson, 2025) is that the PiR state represents the first attached state that triggers further cleft closure and the working stroke, while Pi progression through the tunnel still allows the reversal of Pi release. Taking into account the finding that the force and number of attached motors change in proportion during the increase in isometric force (Brunello et al., 2006; Caremani et al., 2008), it seems excluded that this state is a strongly attached, stiffness-generating state unless it coincides with the AM₁ADPPi state in Figure 12, which is in rapid equilibrium with the higher force-generating states. On the other hand, if it is a weakly bound state, it coincides with the AMADPPi state, upstream with respect to step 3, which we cannot distinguish from the MADPPi state.

It must be noted that identifying the PiR state with the first attached state and, thus, with the AM₁ADPPi state implies a contradiction with the results of sarcomere-level mechanical experiments on the effect of temperature on the force-generating transition of attached motors (Decostre et al., 2005; Linari et al., 2007). In those experiments, it was demonstrated that increasing the temperature potentiates the force developed by a motor during isometric contraction, at the expense of the size of the working stroke elicited by a stepwise decrease in force to zero superimposed on the steady isometric contraction. This indicates that a myosin motor in situ exploits the same working stroke mechanism for both the limited movement, underpinning the increase in strain and force in isometric conditions, and the maximum movement, underpinning the sliding under low load.

The second question concerns whether the conformation dependence of the Pi release rate can be structurally explained by the interactions of Pi within the tunnel after its release from the active site. This would be the case if the interaction kinetics and, thus, the time taken to appear in solution depended on the degree of progression of the working stroke.

In a first attempt to exploit the Llinas et al. crystallographic model for a kinetic “multistep Pi-release” model (Moretto et al., 2022), it is assumed that the working stroke may occur starting not only from the PiR state but also from any of the Pi states following the PiR state or from the AMADP state. In this respect, the model supports the idea proposed by Caremani et al. (2013) that the working stroke and Pi release are independent processes, with the simulated working stroke appearing unaffected by Pi. However, the simulation of the working stroke under unloaded conditions is not obtained by perturbing the equilibrium distribution of motors at the steady state of isometric contraction but by starting from a pre-working stroke state. Under this condition, the output of a “multistep Pi-release” model should consist of a mixture of working strokes equal to the number of states from which the working stroke occurs, and if they are separated by a time-consuming and Pi-sensitive process, such as the progression of Pi through the tunnel, the strokes cannot occur in parallel, and the resulting movement would be temporally dispersed and Pi-sensitive. Presently, the Llinas et al. (2015) crystallographic model is still untested for its functionally sound aspects in relation to the multifaceted mechanical and kinetic aspects of the in situ actin–myosin chemo-mechanical coupling recollected in this report, such as (i) the load dependence of the rate of the working stroke, (ii) the conformation dependence of the rate of Pi release accounting for the Pi transient kinetics, (iii) the much larger inhibitory effect of the increase of [Pi] on the number of attached motors than on the ATPase rate in contractions at high load, and (iv) the limited increase in the rate of energy consumption to achieve the observed power production during shortening.

Author contributions

MC: Writing – original draft, Writing – review and editing, Software. IP: Writing – review and editing, Software. IM: Writing – review and editing. PB: Writing – review and editing, Software. MR: Writing – review and editing, Software. GP: Writing – review and editing, Software. VL: Writing – review and editing, Writing – original draft. ML: Writing – original draft, Writing – review and editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Italian Ministry of University and Research (MUR) in the context of the National Recovery and Resilience Plan (NRRP), Investment PE8—Project Age-It: “Aging Well in an Aging Society”, PRIN2022-PNRR (P2022XPT32, CUP:B53D23033280001) and PRIN2022 (P2022A4XZA2, CUP:B53D23033280001).

Acknowledgments

The authors thank Mario Dolfi and Lorenzo Bongini for their support in the development of the model algorithms.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2025.1664568/full#supplementary-material

References

Barclay C. J., Woledge R. C., Curtin N. A. (2010). Inferring crossbridge properties from skeletal muscle energetics. Prog. Biophys. Mol. Biol. 102 (1), 53–71. doi:10.1016/j.pbiomolbio.2009.10.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Bowater R., Sleep J. (1988). Demembranated muscle fibers catalyze a more rapid exchange between phosphate and adenosine triphosphate than actomyosin subfragment 1. Biochemistry 27 (14), 5314–5323. doi:10.1021/bi00414a055

PubMed Abstract | CrossRef Full Text | Google Scholar

Brenner B. (1988). Effect of Ca2+ on cross-bridge turnover kinetics in skinned single rabbit psoas fibers: implications for regulation of muscle contraction. Proc. Natl. Acad. Sci. U. S. A. 85 (9), 3265–3269. doi:10.1073/pnas.85.9.3265

PubMed Abstract | CrossRef Full Text | Google Scholar

Brunello E., Bianco P., Piazzesi G., Linari M., Reconditi M., Panine P., et al. (2006). Structural changes in the myosin filament and cross-bridges during active force development in single intact frog muscle fibres: stiffness and X-ray diffraction measurements. J. Physiol. 577 (Pt 3), 971–984. doi:10.1113/jphysiol.2006.115394

PubMed Abstract | CrossRef Full Text | Google Scholar

Brunello E., Caremani M., Melli L., Linari M., Fernandez-Martinez M., Narayanan T., et al. (2014). The contributions of filaments and cross-bridges to sarcomere compliance in skeletal muscle. J. Physiol. 592 (17), 3881–3899. doi:10.1113/jphysiol.2014.276196

PubMed Abstract | CrossRef Full Text | Google Scholar

Capitanio M., Canepari M., Maffei M., Beneventi D., Monico C., Vanzi F., et al. (2012). Ultrafast force-clamp spectroscopy of single molecules reveals load dependence of myosin working stroke. Nat. Methods 9 (10), 1013–1019. doi:10.1038/nmeth.2152

PubMed Abstract | CrossRef Full Text | Google Scholar

Caremani M., Dantzig J., Goldman Y. E., Lombardi V., Linari M. (2008). Effect of inorganic phosphate on the force and number of myosin cross-bridges during the isometric contraction of permeabilized muscle fibers from rabbit psoas. Biophys. J. 95 (12), 5798–5808. doi:10.1529/biophysj.108.130435

PubMed Abstract | CrossRef Full Text | Google Scholar

Caremani M., Melli L., Dolfi M., Lombardi V., Linari M. (2013). The working stroke of the myosin II motor in muscle is not tightly coupled to release of orthophosphate from its active site. J. Physiol. 591 (20), 5187–5205. doi:10.1113/jphysiol.2013.257410

PubMed Abstract | CrossRef Full Text | Google Scholar

Caremani M., Melli L., Dolfi M., Lombardi V., Linari M. (2015). Force and number of myosin motors during muscle shortening and the coupling with the release of the ATP hydrolysis products. J. Physiol. 593 (15), 3313–3332. doi:10.1113/JP270265

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen Y. D., Brenner B. (1993). On the regeneration of the actin-myosin power stroke in contracting muscle. Proc. Natl. Acad. Sci. U. S. A. 90 (11), 5148–5152. doi:10.1073/pnas.90.11.5148

PubMed Abstract | CrossRef Full Text | Google Scholar

Cooke R., Pate E. (1985). The effects of ADP and phosphate on the contraction of muscle fibers. Biophys. J. 48 (5), 789–798. doi:10.1016/S0006-3495(85)83837-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Cooke R., Franks K., Luciani G. B., Pate E. (1988). The inhibition of rabbit skeletal muscle contraction by hydrogen ions and phosphate. J. Physiol. 395, 77–97. doi:10.1113/jphysiol.1988.sp016909

PubMed Abstract | CrossRef Full Text | Google Scholar

Dantzig J. A., Goldman Y. E., Millar N. C., Lacktis J., Homsher E. (1992). Reversal of the cross-bridge force-generating transition by photogeneration of phosphate in rabbit psoas muscle fibres. J. Physiol. 451, 247–278. doi:10.1113/jphysiol.1992.sp019163

PubMed Abstract | CrossRef Full Text | Google Scholar

De La Cruz E. M., Wells A. L., Rosenfeld S. S., Ostap E. M., Sweeney H. L. (1999). The kinetic mechanism of myosin V. Proc. Natl. Acad. Sci. U. S. A. 96 (24), 13726–13731. doi:10.1073/pnas.96.24.13726

PubMed Abstract | CrossRef Full Text | Google Scholar

Debold E. P., Walcott S., Woodward M., Turner M. A. (2013). Direct observation of phosphate inhibiting the force-generating capacity of a miniensemble of myosin molecules. Biophys. J. 105 (10), 2374–2384. doi:10.1016/j.bpj.2013.09.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Decostre V., Bianco P., Lombardi V., Piazzesi G. (2005). Effect of temperature on the working stroke of muscle myosin. Proc. Natl. Acad. Sci. U. S. A. 102 (39), 13927–13932. doi:10.1073/pnas.0506795102

PubMed Abstract | CrossRef Full Text | Google Scholar

Ford L. E., Huxley A. F., Simmons R. M. (1977). Tension responses to sudden length change in stimulated frog muscle fibres near slack length. J. Physiol. 269 (2), 441–515. doi:10.1113/jphysiol.1977.sp011911

PubMed Abstract | CrossRef Full Text | Google Scholar

Fusi L., Brunello E., Reconditi M., Piazzesi G., Lombardi V. (2014). The non-linear elasticity of the muscle sarcomere and the compliance of myosin motors. J. Physiol. 592 (5), 1109–1118. doi:10.1113/jphysiol.2013.265983

PubMed Abstract | CrossRef Full Text | Google Scholar

Geeves M. A., Holmes K. C. (2005). The molecular mechanism of muscle contraction. Adv. Protein Chem. 71, 161–193. doi:10.1016/S0065-3233(04)71005-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Hibberd M. G., Trentham D. R. (1986). Relationships between chemical and mechanical events during muscular contraction. Annu. Rev. Biophys. Biophys. Chem. 15, 119–161. doi:10.1146/annurev.bb.15.060186.001003

PubMed Abstract | CrossRef Full Text | Google Scholar

Hibberd M. G., Dantzig J. A., Trentham D. R., Goldman Y. E. (1985). Phosphate release and force generation in skeletal muscle fibers. Science 228 (4705), 1317–1319. doi:10.1126/science.3159090

PubMed Abstract | CrossRef Full Text | Google Scholar

Hill A. V. (1938). The heat of shortening and the dynamic constants of muscle. Proc. R. Soc. Lond B Biol. Sci. 126, 136–195. doi:10.1098/rspb.1938.0050

CrossRef Full Text | Google Scholar

Homsher E., Lacktis J., Regnier M. (1997). Strain-dependent modulation of phosphate transients in rabbit skeletal muscle fibers. Biophys. J. 72 (4), 1780–1791. doi:10.1016/S0006-3495(97)78824-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Houdusse A., Szent-Gyorgyi A. G., Cohen C. (2000). Three conformational states of scallop myosin S1. Proc. Natl. Acad. Sci. U. S. A. 97 (21), 11238–11243. doi:10.1073/pnas.200376897

PubMed Abstract | CrossRef Full Text | Google Scholar

Huxley A. F., Lombardi V. (1980). A sensitive force transducer with resonant frequency 50 kHz. J. Physiol. 305, 15–16P.

Google Scholar

Huxley A. F., Simmons R. M. (1971). Proposed mechanism of force generation in striated muscle. Nature 233 (5321), 533–538. doi:10.1038/233533a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Huxley A. F., Tideswell S. (1997). Rapid regeneration of power stroke in contracting muscle by attachment of second myosin head. J. Muscle Res. Cell Motil. 18 (1), 111–114. doi:10.1023/a:1018641218961

PubMed Abstract | CrossRef Full Text | Google Scholar

Huxley A. F., Lombardi V., Peachey L. D. (1981). A system for fast recording of longitudinal displacement of a striated muscle fibre. J. Physiol. 317, 12P–13P.

Google Scholar

Huxley H. E., Stewart A., Sosa H., Irving T. (1994). X-ray diffraction measurements of the extensibility of actin and myosin filaments in contracting muscle. Biophys. J. 67 (6), 2411–2421. doi:10.1016/S0006-3495(94)80728-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Irving M., Lombardi V., Piazzesi G., Ferenczi M. A. (1992). Myosin head movements are synchronous with the elementary force-generating process in muscle. Nature 357 (6374), 156–158. doi:10.1038/357156a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kawai M., Halvorson H. R. (1991). Two step mechanism of phosphate release and the mechanism of force generation in chemically skinned fibers of rabbit psoas muscle. Biophys. J. 59 (2), 329–342. doi:10.1016/S0006-3495(91)82227-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kerrick W. G., Xu Y. (2004). Inorganic phosphate affects the pCa-force relationship more than the pCa-ATPase by increasing the rate of dissociation of force generating cross-bridges in skinned fibers from both EDL and soleus muscles of the rat. J. Muscle Res. Cell Motil. 25 (2), 107–117. doi:10.1023/b:jure.0000035841.04314.16

PubMed Abstract | CrossRef Full Text | Google Scholar

Linari M., Caremani M., Piperio C., Brandt P., Lombardi V. (2007). Stiffness and fraction of myosin motors responsible for active force in permeabilized muscle fibers from rabbit psoas. Biophys. J. 92 (7), 2476–2490. doi:10.1529/biophysj.106.099549

PubMed Abstract | CrossRef Full Text | Google Scholar

Linari M., Piazzesi G., Lombardi V. (2009). The effect of myofilament compliance on kinetics of force generation by myosin motors in muscle. Biophys. J. 96 (2), 583–592. doi:10.1016/j.bpj.2008.09.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Linari M., Caremani M., Lombardi V. (2010). A kinetic model that explains the effect of inorganic phosphate on the mechanics and energetics of isometric contraction of fast skeletal muscle. Proc. Biol. Sci. 277 (1678), 19–27. doi:10.1098/rspb.2009.1498

PubMed Abstract | CrossRef Full Text | Google Scholar

Linari M., Brunello E., Reconditi M., Fusi L., Caremani M., Narayanan T., et al. (2015). Force generation by skeletal muscle is controlled by mechanosensing in myosin filaments. Nature 528 (7581), 276–279. doi:10.1038/nature15727

PubMed Abstract | CrossRef Full Text | Google Scholar

Llinas P., Isabet T., Song L., Ropars V., Zong B., Benisty H., et al. (2015). How actin initiates the motor activity of myosin. Dev. Cell 33 (4), 401–412. doi:10.1016/j.devcel.2015.03.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Lombardi V., Piazzesi G. (1990). The contractile response during steady lengthening of stimulated frog muscle fibres. J. Physiol. 431, 141–171. doi:10.1113/jphysiol.1990.sp018324

PubMed Abstract | CrossRef Full Text | Google Scholar

Lombardi V., Piazzesi G., Linari M. (1992). Rapid regeneration of the actin-myosin power stroke in contracting muscle. Nature 355 (6361), 638–641. doi:10.1038/355638a0

PubMed Abstract | CrossRef Full Text | Google Scholar

Mansson A. (2025). Theoretical treatment of tension transients in muscle following sudden changes in orthophosphate concentration: implications for energy transduction. J. Muscle Res. Cell Motil. doi:10.1007/s10974-025-09698-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Mansson A., Usaj M., Moretto L., Matusovsky O., Velayuthan L. P., Friedman R., et al. (2023). New paradigms in actomyosin energy transduction: critical evaluation of non-traditional models for orthophosphate release. Bioessays 45 (9), e2300040. doi:10.1002/bies.202300040

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcucci L., Reggiani C. (2016). Mechanosensing in myosin filament solves a 60 years old conflict in skeletal muscle modeling between high power output and slow rise in tension. Front. Physiol. 7, 427. doi:10.3389/fphys.2016.00427

PubMed Abstract | CrossRef Full Text | Google Scholar

Marcucci L., Fukunaga H., Yanagida T., Iwaki M. (2021). The synergic role of actomyosin architecture and biased detachment in muscle energetics: insights in cross bridge mechanism beyond the lever-arm swing. Int. J. Mol. Sci. 22 (13), 7037. doi:10.3390/ijms22137037

PubMed Abstract | CrossRef Full Text | Google Scholar

Millar N. C., Homsher E. (1990). The effect of phosphate and calcium on force generation in glycerinated rabbit skeletal muscle fibers. A steady-state and transient kinetic study. J. Biol. Chem. 265 (33), 20234–20240. doi:10.1016/s0021-9258(17)30494-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Millar N. C., Homsher E. (1992). Kinetics of force generation and phosphate release in skinned rabbit soleus muscle fibers. Am. J. Physiol. 262 (5 Pt 1), C1239–C1245. doi:10.1152/ajpcell.1992.262.5.C1239

PubMed Abstract | CrossRef Full Text | Google Scholar

Moretto L., Usaj M., Matusovsky O., Rassier D. E., Friedman R., Mansson A. (2022). Multistep orthophosphate release tunes actomyosin energy transduction. Nat. Commun. 13 (1), 4575. doi:10.1038/s41467-022-32110-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Muretta J. M., Rohde J. A., Johnsrud D. O., Cornea S., Thomas D. D. (2015). Direct real-time detection of the structural and biochemical events in the myosin power stroke. Proc. Natl. Acad. Sci. U. S. A. 112 (46), 14272–14277. doi:10.1073/pnas.1514859112

PubMed Abstract | CrossRef Full Text | Google Scholar

Nyitrai M., Geeves M. A. (2004). Adenosine diphosphate and strain sensitivity in myosin motors. Philos. Trans. R. Soc. Lond B Biol. Sci. 359 (1452), 1867–1877. doi:10.1098/rstb.2004.1560

PubMed Abstract | CrossRef Full Text | Google Scholar

Offer G., Ranatunga K. W. (2020). The location and rate of the phosphate release step in the muscle cross-bridge cycle. Biophys. J. 119 (8), 1501–1512. doi:10.1016/j.bpj.2020.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Pate E., Cooke R. (1989a). Addition of phosphate to active muscle-fibers probes actomyosin states within the powerstroke. Pflugers Archiv-European J. Physiology 414 (1), 73–81. doi:10.1007/Bf00585629

PubMed Abstract | CrossRef Full Text | Google Scholar

Pate E., Cooke R. (1989b). A model of crossbridge action: the effects of ATP, ADP and Pi. J. Muscle Res. Cell Motil. 10 (3), 181–196. doi:10.1007/BF01739809

PubMed Abstract | CrossRef Full Text | Google Scholar

Pertici I., Bongini L., Melli L., Bianchi G., Salvi L., Falorsi G., et al. (2018). A myosin II nanomachine mimicking the striated muscle. Nat. Commun. 9 (1), 3532. doi:10.1038/s41467-018-06073-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Pertici I., Bianchi G., Bongini L., Lombardi V., Bianco P. (2020). A myosin II-Based nanomachine devised for the study of Ca(2+)-Dependent mechanisms of muscle regulation. Int. J. Mol. Sci. 21 (19), 7372. doi:10.3390/ijms21197372

PubMed Abstract | CrossRef Full Text | Google Scholar

Pertici I., Bongini L., Caremani M., Reconditi M., Linari M., Piazzesi G., et al. (2023). Matching mechanics and energetics of muscle contraction suggests unconventional chemomechanical coupling during the actin-myosin interaction. Int. J. Mol. Sci. 24 (15), 12324. doi:10.3390/ijms241512324

PubMed Abstract | CrossRef Full Text | Google Scholar

Piazzesi G., Lombardi V. (1995). A cross-bridge model that is able to explain mechanical and energetic properties of shortening muscle. Biophys. J. 68 (5), 1966–1979. doi:10.1016/S0006-3495(95)80374-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Piazzesi G., Linari M., Reconditi M., Vanzi F., Lombardi V. (1997). Cross-bridge detachment and attachment following a step stretch imposed on active single frog muscle fibres. J. Physiol. 498 (Pt 1), 3–15. doi:10.1113/jphysiol.1997.sp021837

PubMed Abstract | CrossRef Full Text | Google Scholar

Piazzesi G., Lucii L., Lombardi V. (2002). The size and the speed of the working stroke of muscle myosin and its dependence on the force. J. Physiol. 545 (Pt 1), 145–151. doi:10.1113/jphysiol.2002.028969

PubMed Abstract | CrossRef Full Text | Google Scholar

Piazzesi G., Dolfi M., Brunello E., Fusi L., Reconditi M., Bianco P., et al. (2014). The myofilament elasticity and its effect on kinetics of force generation by the myosin motor. Arch. Biochem. Biophys. 552-553, 108–116. doi:10.1016/j.abb.2014.02.017

PubMed Abstract | CrossRef Full Text | Google Scholar

Piazzesi G., Caremani M., Linari M., Reconditi M., Lombardi V. (2018). Thick filament mechano-sensing in skeletal and cardiac muscles: a common mechanism able to adapt the energetic cost of the contraction to the task. Front. Physiol. 9, 736. doi:10.3389/fphys.2018.00736

PubMed Abstract | CrossRef Full Text | Google Scholar

Potma E. J., Stienen G. J. (1996). Increase in ATP consumption during shortening in skinned fibres from rabbit psoas muscle: effects of inorganic phosphate. J. Physiol. 496 (Pt 1), 1–12. doi:10.1113/jphysiol.1996.sp021660

PubMed Abstract | CrossRef Full Text | Google Scholar

Potma E. J., van Graas I. A., Stienen G. J. (1995). Influence of inorganic phosphate and pH on ATP utilization in fast and slow skeletal muscle fibers. Biophys. J. 69 (6), 2580–2589. doi:10.1016/S0006-3495(95)80129-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Ranatunga K. W. (1999). Effects of inorganic phosphate on endothermic force generation in muscle. Proc. Biol. Sci. 266 (1426), 1381–1385. doi:10.1098/rspb.1999.0791

PubMed Abstract | CrossRef Full Text | Google Scholar

Ranatunga K. W., Coupland M. E., Mutungi G. (2002). An asymmetry in the phosphate dependence of tension transients induced by length perturbation in Mammalian (rabbit psoas) muscle fibres. J. Physiol. 542 (Pt 3), 899–910. doi:10.1113/jphysiol.2002.019471

PubMed Abstract | CrossRef Full Text | Google Scholar

Rassier D. E., Mansson A. (2025). Mechanisms of myosin II force generation: insights from novel experimental techniques and approaches. Physiol. Rev. 105 (1), 1–93. doi:10.1152/physrev.00014.2023

PubMed Abstract | CrossRef Full Text | Google Scholar

Rayment I., Holden H. M., Whittaker M., Yohn C. B., Lorenz M., Holmes K. C., et al. (1993a). Structure of the actin-myosin complex and its implications for muscle contraction. Science 261 (5117), 58–65. doi:10.1126/science.8316858

PubMed Abstract | CrossRef Full Text | Google Scholar

Rayment I., Rypniewski W. R., Schmidt-Base K., Smith R., Tomchick D. R., Benning M. M., et al. (1993b). Three-dimensional structure of myosin subfragment-1: a molecular motor. Science 261 (5117), 50–58. doi:10.1126/science.8316857

PubMed Abstract | CrossRef Full Text | Google Scholar

Reconditi M., Linari M., Lucii L., Stewart A., Sun Y. B., Boesecke P., et al. (2004). The myosin motor in muscle generates a smaller and slower working stroke at higher load. Nature 428 (6982), 578–581. doi:10.1038/nature02380

PubMed Abstract | CrossRef Full Text | Google Scholar

Robert-Paganin J., Pylypenko O., Kikuti C., Sweeney H. L., Houdusse A. (2020). Force generation by myosin motors: a structural perspective. Chem. Rev. 120 (1), 5–35. doi:10.1021/acs.chemrev.9b00264

PubMed Abstract | CrossRef Full Text | Google Scholar

Sleep J., Irving M., Burton K. (2005). The ATP hydrolysis and phosphate release steps control the time course of force development in rabbit skeletal muscle. J. Physiol. 563 (Pt 3), 671–687. doi:10.1113/jphysiol.2004.078873

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith D. A. (2014). A new mechanokinetic model for muscle contraction, where force and movement are triggered by phosphate release. J. Muscle Res. Cell Motil. 35 (5), 295–306. doi:10.1007/s10974-014-9391-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith C. A., Rayment I. (1996). X-ray structure of the magnesium(II).ADP.vanadate complex of the Dictyostelium discoideum myosin motor domain to 1.9 A resolution. Biochemistry 35 (17), 5404–5417. doi:10.1021/bi952633+

PubMed Abstract | CrossRef Full Text | Google Scholar

Stehle R. (2017). Force responses and sarcomere dynamics of cardiac myofibrils induced by rapid changes in [Pi]. Biophys. J. 112 (2), 356–367. doi:10.1016/j.bpj.2016.11.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Stehle R., Tesi C. (2017). Kinetic coupling of phosphate release, force generation and rate-limiting steps in the cross-bridge cycle. J. Muscle Res. Cell Motil. 38 (3-4), 275–289. doi:10.1007/s10974-017-9482-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Takagi Y., Shuman H., Goldman Y. E. (2004). Coupling between phosphate release and force generation in muscle actomyosin. Philos. Trans. R. Soc. Lond B Biol. Sci. 359 (1452), 1913–1920. doi:10.1098/rstb.2004.1561

PubMed Abstract | CrossRef Full Text | Google Scholar

Tesi C., Colomo F., Nencini S., Piroddi N., Poggesi C. (2000). The effect of inorganic phosphate on force generation in single myofibrils from rabbit skeletal muscle. Biophys. J. 78 (6), 3081–3092. doi:10.1016/S0006-3495(00)76845-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Wakabayashi K., Sugimoto Y., Tanaka H., Ueno Y., Takezawa Y., Amemiya Y. (1994). X-ray diffraction evidence for the extensibility of actin and myosin filaments during muscle contraction. Biophys. J. 67 (6), 2422–2435. doi:10.1016/S0006-3495(94)80729-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Walker J. W., Lu Z., Moss R. L. (1992). Effects of Ca2+ on the kinetics of phosphate release in skeletal muscle. J. Biol. Chem. 267 (4), 2459–2466. doi:10.1016/s0021-9258(18)45901-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Webb M. R., Hibberd M. G., Goldman Y. E., Trentham D. R. (1986). Oxygen exchange between Pi in the medium and water during ATP hydrolysis mediated by skinned fibers from rabbit skeletal muscle. Evidence for Pi binding to a force-generating state. J. Biol. Chem. 261 (33), 15557–15564. doi:10.1016/s0021-9258(18)66750-x

PubMed Abstract | CrossRef Full Text | Google Scholar

West T. G., Ferenczi M. A., Woledge R. C., Curtin N. A. (2005). Influence of ionic strength on the time course of force development and phosphate release by dogfish muscle fibres. J. Physiol. 567 (Pt 3), 989–1000. doi:10.1113/jphysiol.2005.087106

PubMed Abstract | CrossRef Full Text | Google Scholar

White H. D., Taylor E. W. (1976). Energetics and mechanism of actomyosin adenosine triphosphatase. Biochemistry 15 (26), 5818–5826. doi:10.1021/bi00671a020

PubMed Abstract | CrossRef Full Text | Google Scholar

Woody M. S., Winkelmann D. A., Capitanio M., Ostap E. M., Goldman Y. E. (2019). Single molecule mechanics resolves the earliest events in force generation by cardiac myosin. Elife 8, e49266. doi:10.7554/eLife.49266

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeng W., Conibear P. B., Dickens J. L., Cowie R. A., Wakelin S., Malnasi-Csizmadia A., et al. (2004). Dynamics of actomyosin interactions in relation to the cross-bridge cycle. Philos. Trans. R. Soc. Lond B Biol. Sci. 359 (1452), 1843–1855. doi:10.1098/rstb.2004.1527

PubMed Abstract | CrossRef Full Text | Google Scholar