Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC2896253

Formats

Article sections

- Abstract
- 1 Introduction
- 2 The two-state model
- 3 Perturbation
- 4 Tension transients and delayed tension
- 5 The expanded two-state model
- 6 The three-state model
- 7 Three-state model, a simple case analysis (α + α′ β + β′ γ + γ′)
- 8 Steady state
- 9 The mechanisms of tension transients that satisfy the Le Châtelier–Brown Principle
- 10 Effects of MgATP and phosphate (Pi)
- 11 Substrate-binding step with a conformational change
- 12 ADP as a competitive inhibitor of ATP
- 13 Phosphate release step with a conformational change
- 14 Combined cross-bridge model
- 15 Force generation step, isometric tension, and work performance
- 16 Exponential process D and one extra cross-bridge state
- 17 Stiffness
- 18 ATP hydrolysis rate and the rate constant of step 6
- 19 Series compliance, force development, and the slowest step in the cross-bridge cycle
- 20 More realistic conditions
- 21 Cardiac muscle fibers
- 22 Conclusion and future direction
- References

Authors

Related links

J Muscle Res Cell Motil. Author manuscript; available in PMC 2010 July 2.

Published in final edited form as:

Published online 2008 April 19. doi: 10.1007/s10974-008-9131-3

PMCID: PMC2896253

NIHMSID: NIHMS207701

Masataka Kawai, Department of Anatomy and Cell Biology, University of Iowa, Iowa City, IA 52245, USA;

Masataka Kawai: ude.awoiu@iawak-akatasaM; Herbert R. Halvorson: moc.loa@531rovlaHH

The publisher's final edited version of this article is available at J Muscle Res Cell Motil

See other articles in PMC that cite the published article.

Two- and three-state cross-bridge models are considered and examined with respect to their ability to predict three distinct phases of the force transients that occur in response to step change in muscle fiber length. Particular attention is paid to satisfying the Le Châtelier–Brown Principle. This analysis shows that the two-state model can account for phases 1 and 2 of a force transient, but is barely adequate to account for phase 3 (delayed force) unless a stretch results in a sudden increase in the number of cross-bridges in the detached state. The three-state model (**A** → **B** → **C** → **A**) makes it possible to account for all three phases if we assume that the **A** → **B** transition is fast (corresponding to phase 2), the **B** → **C** transition is of intermediate speed (corresponding to phase 3), and the **C** → **A** transition is slow; in such a scenario, states **A** and **C** can support or generate force (high force states) but state **B** cannot (detached, or low-force state). This model involves at least one ratchet mechanism. In this model, force can be generated by either of two transitions: **B** → **A** or **B** → **C**. To determine which of these is the major force-generating step that consumes ATP and transduces energy, we examine the effects of ATP, ADP, and phosphate (Pi) on force transients. In doing so, we demonstrate that the fast transition (phase 2) is associated with the nucleotide-binding step, and that the intermediate-speed transition (phase 3) is associated with the Pi-release step. To account for all the effects of ligands, it is necessary to expand the three-state model into a six-state model that includes three ligand-bound states. The slowest phase of a force transient (phase 4) cannot be explained by any of the models described unless an additional mechanism is introduced. Here we suggest a role of series compliance to account for this phase, and propose a model that correlates the slowest step of the cross-bridge cycle (transition **C** → **A**) to: phase 4 of step analysis, the rate constant *k*_{tr} of the quick-release and restretch experiment, and the rate constant *k*_{act} for force development time course following Ca^{2+} activation.

Numerous contraction models involving multiple cross-bridge states have been proposed to account for isometric tension and tension transients, but many of them are not easy to follow due to their mathematical complexity and the numerical approximations that are used for models. Indeed, the numerical approach makes it especially difficult to grasp the significance and usefulness of a particular model. Moreover, some models address only a part of the cycling scheme, and others rely on assumptions that are not clearly delineated. In an extreme case, too many cross-bridge states are employed and too many parameters are used to fit scanty data, which makes the reader wonder how much of the model is actually supported by the data on fiber studies, how much is an extrapolation from the data obtained from solution studies, and how much is pure speculation. Is a model with more cross-bridge states necessarily better than one with fewer? We believe that the purpose of the model is to account for experimental results, and to predict results for future experimentation. For this reason, we consider simplicity to be one of the key aspects of model building. In this mini-review, we focus on simple cases and attempt to explain some fundamental experimental observations. The mathematical principles used here are not complex and are, in fact, taught in high school and college courses. Hence if there are difficulties, one should consult with a math teacher or look into appropriate textbooks. Alternatively, the authors are always available for consultation. The cases covered here are limited to the isometric, or near-isometric, state of contraction; isotonic contraction is not dealt with because it introduces extra complexity. The mathematical symbols used in the text are summarized in Table 1.

This is the simplest model of all and, hence, easiest for beginners to understand while still of value for experts because it involves fundamental principles. The two-state model is extensively used in muscle physiology and in mechanics (e.g., Huxley 1957; Thorson and White 1969; Huxley and Simmons 1971; Abbot 1973; Brenner 1988; Palmer et al. 2007). When a reactant (**A**) changes into a product (**B**) with the forward rate constant *α* and its reversal rate constant *α*′, this kinetic scheme can be written as:

For application to muscle mechanics, state **A** is assigned to the attached “high force state”, and state **B** to the detached “no force state” (or low force state). The rate of the forward reaction is *αA*, and the rate of the reversal reaction is *α*′*B*, where *A* and *B* represent the concentrations of the respective enzymatic species. Although concentration is usually written as [A] and [B], here the brackets are left out for simplicity and *italics* are used to denote mathematical variables. In contrast, when discussing molecular species we use bold letters such as **A** and **B**. As *A* and *B* are functions of time (*t*), they can be written as *A* = *A*(*t*) and *B* = *B*(*t*). In Scheme 1, because the forward rate is a loss and the reverse rate is a gain for **A**,

$$dA/dt=-\alpha A+{\alpha}^{\prime}B$$

(1)

The reaction in Scheme 1 is called the first-order reaction because it is represented by the first-order differential equation (Eq. 1). Given that **A** and **B** are two forms of the same molecular species,

$$A+B={A}_{\text{T}}$$

(2)

where *A*_{T} is the total concentration and is conserved (i.e., it does not change with time). In the case of concentrated solutions, the word “activity” is used in place of the concentration because not all molecules are chemically active. The consequent elimination of *B* from Eqs. 1 and 2 results in:

$$dA/dt+\lambda A={\alpha}^{\prime}{A}_{\text{T}}$$

(3)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\lambda \equiv \alpha +{\alpha}^{\prime}$$

(4)

Note that indicates definition. Equation 3 takes the form of Eq. A1 in Appendix 1, and it can be solved (Eq. A6) to result in:

$$A(t)={A}_{0}exp(-\lambda t)+{A}_{1}$$

(5)

$$\text{and}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}B(t)=-{A}_{0}exp(-\lambda t)+{B}_{1}$$

(6)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{A}_{1}\equiv ({\alpha}^{\prime}/\lambda ){A}_{\text{T}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\text{and}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{B}_{1}\equiv (\alpha /\lambda ){A}_{\text{T}}$$

(7)

Thus, the first-order reaction of Scheme 1 results in a time course with one exponential process whose rate constant is *λ*, which is the sum of the forward rate constant *α* and the reverse rate constant *α*′ (Eq. 4). *λ* is termed the “apparent” rate constant, indicating that this is the rate constant experimentally observed as in Eq. 5 or 6. In contrast *α* and *α*′ are termed the “intrinsic” or “fundamental” rate constants (Gutfreund 1995). Although the adjective “intrinsic”, “fundamental” or “apparent” is often left out, which of these would apply is usually self evident. *A*_{0} is an integration constant that is determined by the initial condition, and it is termed the amplitude (or magnitude) of the exponential process in Eqs. 5 and 6.

As time approaches infinity (*t* → ∞), Scheme 1 achieves an equilibrium, resulting in *dA*(*t*)/*dt* = 0 (in Eqs. 1 and 3). From Eqs. 5 and 6, it becomes clear that *A*(∞) = *A*_{1} and *B*(∞) = *B*_{1}. Thus, *K* *α*/*α*′ = *B*_{1}/*A*_{1} is the equilibrium constant, and *A*_{1} and *B*_{1} are the equilibrium concentrations of individual molecular species. *K* = *B*_{1}/*A*_{1} is the well-known law of mass action.

A perturbation (*δx*) is experimentally applied after the equilibrium described in Scheme 1 (*A* = *A*_{1}, *B* = *B*_{1}) is reached. The perturbation can be a length change (e.g., Thorson and White 1969; Huxley and Simmons 1971; Abbott 1973; Heinl et al. 1974; Kawai 1978), a temperature change (*e.g*., Bershitsky and Tsaturyan 1992; Davis and Rodgers 1995; Ranatunga 1999), a pressure change (e.g., Fortune et al. 1991), etc., but it must be small and fast. A perturbation can also be a sudden change in the ligand (ATP, ADP, Pi) concentration, such as that used in caged Pi (e.g., Homsher and Millar 1990; Dantzig et al. 1992; Walker et al. 1992; Araujo and Walker 1996), caged ATP (e.g., Goldman et al. 1984; Martin and Barsotti 1994), caged Ca (e.g., Preston et al. 2007), or rapid solution-switch (e.g., Tesi et al. 2000, 2002; Stehle et al. 2002; Piroddi et al. 2003) experiments. The purpose of the perturbation is to probe the rate constant(s) and to construct a kinetic scheme. In Scheme 1, the perturbation is assumed to modify the rate constants *α* and *α*′:

$$\alpha \to \alpha +\delta \alpha $$

(8)

$${\alpha}^{\prime}\to {\alpha}^{\prime}+\delta {\alpha}^{\prime}$$

(9)

$$\text{Therefore},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\lambda \to \lambda +\delta \lambda ,\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\delta \lambda =\delta \alpha +\delta {\alpha}^{\prime}$$

(10)

The changes described in Eqs. 8–10 are assumed to occur simultaneous to perturbation *δx*, with *δα* being the sensitivity of *α*, *δα*′ the sensitivity of *α*′, and *δλ* the sensitivity of *λ*, to perturbation *δx*. Because *δx* is a small quantity, all other *δ*s are also small quantities.

The perturbation (Eqs. 8–10) results in a change in *A* that is time dependent:

$${A}_{1}\to {A}_{1}+\delta A(t)$$

(11)

With the substitution of Eqs. 9, 10 and 11 into Eq. 3 we obtain:

$$d({A}_{1}+\delta A)/dt+(\lambda +\delta \lambda )({A}_{1}+\delta A)={A}_{\text{T}}({\alpha}^{\prime}+\delta {\alpha}^{\prime})$$

(12)

Sidewise subtraction of Eq. 3 from Eq. 12 gives:

$$d\delta A/dt+\lambda \delta A\cong {A}_{\text{T}}\delta {\alpha}^{\prime}-{A}_{1}\delta \lambda $$

(13)

The cross term *δλδA* is dropped from Eq. 13 because this is an extremely small quantity. If preferred, Eq. 13 can be derived directly from Eq. 3 by differentiation. The right-hand side of Eq. 13 is a constant of time, because all variables are time-independent quantities. Note the similarity between Eqs. 3 and 13. By solving Eq. 13 (Appendix 1),

$$\delta A(t)={A}_{2}[1-exp(-\lambda t)]$$

(14)

$$\delta B(t)=-\delta A(t)=-{A}_{2}[1-exp(-\lambda t)]$$

(15)

where *A*_{2} (*A _{T}δα*′ −

An experimenter will observe the same rate constant *λ* whether starting a new reaction as in Scheme 1 (Eq. 5), or applying a small perturbation (of any kind) during equilibrium and measuring the transient as in Eq. 14. Although not shown explicitly, this conclusion holds true for any multiple-state model.

When an isometrically contracting muscle fiber is stretched by *δl* (with *δl* < 0 in the case of release), force rises simultaneous to the stretch because of the elastic link of attached cross-bridges. This is termed phase 1 (Fig. 1A; Huxley and Simmons 1971; Heinl et al. 1974). Force then decays quickly (phase 2), after which it again rises, but more gradually (phase 3 in Fig. 2). Phase 3 is particularly evident in insect indirect flight muscles (Fig. 3A; Pringle 1978), and is called “delayed tension”. To explain the delayed tension, Thorson and White (1969, 1983) used the two-state model, but had to assume that either *α*′ increases (as does in Huxley’s model 1957, and in Julian et al.’s model 1974) or *α* decreases with a stretch (*δl*), in other words, that *δα*′/*δl* > 0 or *δα*/*δl* < 0 (see also Abbott 1973; Abbott and Steiger 1977; Pringle 1978). The perturbation analysis used was exactly the same as in Eqs. 8–14 above. This assumption (*δα*′/*δl* > 0 or *δα*/*δl* < 0), however, does not satisfy the Le Châtelier–Brown Principle, which states that a thermodynamic system tends to balance the effects of any stress inflicted upon it (Kirkwood and Oppenheim 1961). In the case of muscle mechanics, this principle is interpreted in the following way.

A sketch of stretch (at A) and release (at B) experiments during a plateau of active tension in frog semitendinosus fibers. Four phases of tension transients are indicated. Modified from Fig. 1 of Heinl et al. (1974) and reproduced with permission from **...**

A record of the force time course that occurs in response to a step-length increase (at ↑, 1.5 nm/half sarcomere, which is about 0.12%) in rabbit psoas fibers at 5°C during Ca^{2+} activation. Numbers indicate four phases of tension transients. **...**

A record of the force time course that occurs in response to a step-length increase (0.1%) (at A) and the ensuing decrease (at B), in insect fibrillar muscle fibers. Three phases of the tension transient are indicated. Phase 3 is most prominent in this **...**

If a cross-bridge that is generating force () is stretched by *δl*, then the work (*δW*) performed by the forcing apparatus on this cross-bridge is:

$$\delta W=(\text{force})\times (\text{distance})=\phi \delta l$$

(16)

Equation 16 also holds true for release (*δl* < 0), in which case *δW* < 0 and the work is performed by the cross-bridge on the forcing apparatus. is force/cross-bridge ( > 0), and the same as the unitary force used in single molecule experiments (e.g. Molloy et al. 1995; Guilford et al. 1997; Kitamura et al. 1999). It is likely that the unitary force is distributed because the mismatch between actin and myosin repeats in most muscles causes the strain on cross-bridges to be distributed. Thus, it is more accurate to state that represents the averaged unitary force over the distributed strain. From here, there are two possibilities: in case 1, the activation energy *E _{α}* for the reaction

$$\delta {E}_{\alpha}=-\delta W=-\phi \delta l$$

(17)

The rate constant *α* is related to the activation energy by the Arrhenius equation (Gutfreund 1995):

$$\alpha ={\alpha}_{0}exp(-{E}_{\alpha}/{k}_{B}T)$$

(18)

where *k _{B}* is the Boltzmann’s constant, and

$$\delta \alpha =-\frac{{\alpha}_{0}}{{k}_{B}T}exp(-\frac{{E}_{\alpha}}{{k}_{B}T})\delta {E}_{\alpha}=-\frac{\alpha}{{k}_{B}T}\delta {E}_{\alpha}$$

(19)

$$\frac{\delta \alpha}{\alpha}=\frac{\phi}{{k}_{B}T}\delta l$$

(20)

For case 2, the chain of events is
$\delta l\to \delta W\to \delta {E}_{\alpha}^{\prime}\to {\alpha}^{\prime}$ and this changes the reversal rate *α*′*B*. Here, because
$\delta {E}_{\alpha}^{\prime}=+\delta W$, the expression equivalent to Eq. 20 is

$$\frac{\delta {\alpha}^{\prime}}{{\alpha}^{\prime}}=-\frac{\phi}{{k}_{B}T}\delta l$$

(21)

From Eq. 20, *δα*/*δl* = *α*/*k _{B}T* > 0 (it cannot be negative). Similarly, from Eq. 21,

When attached cross-bridges **A** are stretched by *δl*, they follow different kinetics from ones that are not stretched. One way of simplifying this complication is to consider that the stretched cross-bridges are in a special state **G**, as shown in Scheme 2. As time progresses, **G** decays at the rate constant *μ* to become **B**, and is not regenerated. This aspect is in contrast to the two-state model discussed in Section 4, in which the stretched state can be regenerated. In Scheme 2, state **A** vanishes at the time of stretch, but because the **B** ↔ **A** transition takes place, **A** reappears over time.

The following differential equations (Eqs. 22–24) describe the kinetics of Scheme 2:

$$dA(t)/dt=-\alpha A+{\alpha}^{\prime}B$$

(22)

$$dB(t)/dt=\alpha A-{\alpha}^{\prime}B+\mu G$$

(23)

$$dG(t)/dt=-\mu G$$

(24)

$$\text{Evidently}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}A(t)+B(t)+G(t)={A}_{\text{T}}$$

(25)

The initial conditions and the final steady state concentrations are:

$$\begin{array}{l}A(-0)={A}_{1},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}A(+0)=0,\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}A(\infty )\equiv {A}_{1}\\ B(-0)={B}_{1},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}B(+0)={B}_{1},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}B(\infty )\equiv {B}_{1}\\ G(-0)=0,\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}G(+0)={A}_{1},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}G(\infty )\equiv {G}_{1}=0\end{array}$$

where *t* = −0 means immediately before the stretch, and *t*= +0 means immediately after the stretch. *A*_{1} and *B*_{1} are defined in Eq. 7. Eqs. 22–25 can be solved in reverse order, i.e. first yielding *G*(*t*) (Appendix 1), then *B*(*t*), and finally *A*(*t*). Alternatively, the method described in Appendix 2 for three states can be used, because the expanded two-state model is a special case of the three-state model. The results are as follows.

$$A(t)={A}_{1}+\frac{{A}_{1}}{{\alpha}^{\prime}-\delta \alpha}[\delta \alpha exp(-\lambda t)-{\alpha}^{\prime}exp(-\mu t)]$$

(26)

$$B(t)={B}_{1}+\frac{\delta \alpha}{{\alpha}^{\prime}-\delta \alpha}{A}_{1}[-exp(-\lambda t)+exp(-\mu t)]$$

(27)

$$G(t)={A}_{1}exp(-\mu t)$$

(28)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\mu \equiv \alpha +\delta \alpha $$

*λ* is defined in Eq. 4. *λ* > *μ*, because *α*′ > *δα* when *δ* is small. Relating a chemical state to force can be simplified by considering a myofibril of cross sectional area *A _{C}* and a half sarcomere of length

$$\begin{array}{l}{F}_{\text{A}}=\phi {AA}_{c}{l}_{0}{N}_{\text{A}}\\ {F}_{\text{G}}=(\phi +\rho \delta l){GA}_{c}{l}_{0}{N}_{\text{A}}\end{array}$$

is the force/cross-bridge (unitary force) in state **A**, + *ρδl* is the force/cross-bridge in state **G**, *ρ* is the spring constant of the cross-bridge, and *N _{A}* is Avogadro’s number. Therefore, in the half sarcomere

$$\begin{array}{l}F(t)={F}_{G}+{F}_{A}=[\phi (A+G)+G\rho \delta l]{A}_{c}{l}_{0}{N}_{A}\\ ={A}_{1}[\phi -zexp(-\mu t)+wexp(-\lambda t)]{A}_{c}{l}_{0}{N}_{A}\end{array}$$

(29)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}z\equiv \phi \delta \alpha /({\alpha}^{\prime}-\delta \alpha )-\rho \delta l$$

(30)

$$\text{and}\phantom{\rule{0.16667em}{0ex}}w\equiv \phi \delta \alpha /({\alpha}^{\prime}-\delta \alpha )$$

(31)

$$\begin{array}{l}F(-0)=\phi {A}_{1}{A}_{c}{l}_{0}{N}_{A}\\ F(+0)={A}_{1}(\rho \delta l+\phi ){A}_{c}{l}_{0}{N}_{A}\\ F(\infty )=\phi {A}_{1}{A}_{c}{l}_{0}{N}_{A}\end{array}$$

Equations 29, 30, and 31 illustrate that a force transient comprises three phases. The first is a rise in the force simultaneous with the stretch. The second is a rapid decay of force at rate constant *μ*, and reflects the decay of **G**. The outcome of the third phase depends on the polarity of *z*: if *z* > 0, then there is a delayed rise in force as it approaches *A*_{1}*A _{c}l*

According to the above analysis, the length change can be seen to induce an extra state **G**, which decays with a unique rate constant (*μ*). Steiger and Abbott (1981) used this line of modeling on a four-state model (with two attached and two detached states after Lymn and Taylor 1971), to generate one with six states. Although we believe that Scheme 2 and its theory represent one of the most rigorous systems for treating the stretched (or released) cross-bridges, they may not allow for practical analysis of experimental results because the apparent rate constants *λ* and *μ* are similar and, hence, difficult to distinguish from one another. In addition, modeling must include one extra cross-bridge state per attached cross-bridge state, which adds to the complexity of the algebra required when multiple attached states are considered. For these reasons, we do not use this line of approach in the following discussions.

To explain delayed tension more explicitly, it is necessary to introduce the three-state model. Given the extreme usefulness of this model, we discuss it in detail. The three state model has been used frequently among muscle physiologists and mechanics (e.g., Julian et al. 1974; Thorson and White 1983; Woledge et al. 1985; Murase et al. 1986; Dickinson et al. 1997; Campbell et al. 2001). Scheme 3 represents its most general form.

This is a cyclic scheme with the forward reaction consisting of the outer clockwise cycle: **A** → **B** → **C** → **A**, for which the respective rate constants are *α*,*β* and *γ*. The reversal reaction consists of the inner counter-clockwise cycle: **A** → **C** → **B** → **A**, with the respective rate constants *γ*′, *β*′ and *α*′. In Scheme 3, we assume that **A** and **C** are attached states, which generate and/or support force, and that **B** is a detached state, which does not support force. Because the analysis of Scheme 3 is very complex, we will discuss a simple case here. The complete analysis of Scheme 3 can be found in the Appendix 2.

Now let us consider a simple case in which *α* + *α*′ *β* + *β*′ *γ* + *γ* is satisfied. In the following discussion, these rate constants are referred to fast, medium, and slow, respectively. The case *α* + *α*′ *β* + *β*′ < *γ* + *γ*′ is included in this discussion by an exchange: *α*, *α*′ ↔ *γ*, *γ*′. The more realistic case (*α* + *α*′ > *β* + *β*′ > *γ* + *γ*′) is discussed later in Section 20. *X* *Y* indicates that *X* is much larger than *Y*, so that *Y* can be ignored when *X* + *Y* is evaluated. An example is 1000 s^{−1} and 100 s^{−1}. If we focus on a fast time scale, comparable to 1/(*α* + *α*′) (such as ~1 ms in the cross-bridge cycle), then steps 2 (**B** ↔ **C**) and 3 (**C** ↔ **A**) are too slow to occur. Hence they could not be witnessed on the time-scale of observation. In this case, only step 1 is observed, and the problem is reduced to Scheme 1, with the apparent rate constant *λ*_{2}.

$${\lambda}_{2}=\alpha +{\alpha}^{\prime}$$

(32)

The apparent rate constant of a fast step (*λ*_{2}) is the sum of the forward rate constant and the reversal rate constant of this step (Eq. 32), and it is not influenced by either the rate or equilibrium constants of slower steps.

If we focus on a medium-speed time scale, comparable to 1/(*β* + *β*′) (such as ~10 ms in the cross-bridge cycle), then step 1 is very fast and looks like an equilibrium, whereas step 3 is too slow to be witnessed. Therefore, in this scenario, Scheme 3 actually looks like Scheme 4.

where *K _{α}*

$$\frac{d(A+B)}{dt}=-\beta B+{\beta}^{\prime}C$$

(33)

That is, **A** and **B** gain/lose as a group. Because **A**, **B** and **C** are three different states of myosin,

$$A+B+C={A}_{\text{T}}$$

(34)

where *A _{T}* is the total concentration of myosin subfragment one (S1). Because step 1 represents an equilibrium, the mass action law applies:

$$B/A={K}_{\alpha}$$

(35)

By eliminating *B* and *C* from Eqs. 33–35, we obtain

$$dA/dt+{\lambda}_{3}A={\beta}^{\prime}{A}_{\text{T}}/(1+{K}_{\alpha})$$

(36)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{\lambda}_{3}\equiv \frac{{K}_{\alpha}}{1+{K}_{\alpha}}\beta +{\beta}^{\prime}$$

(37)

By solving Eq. 36 (Appendix 1), we arrive at:

$$\begin{array}{l}A(t)={A}_{0}exp(-{\lambda}_{3}t)+{A}_{1},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}B(t)={A}_{0}{K}_{\alpha}exp(-{\lambda}_{3}t)+{B}_{1},\\ C(t)=-{A}_{0}(1+{K}_{\alpha})exp(-{\lambda}_{3}t)+{C}_{1}\end{array}$$

(38)

where *A*_{1} *A*_{T}/[1 + *K _{α}*(1 +

The apparent rate constant of a medium-speed step (*λ*_{3}) is a weighted sum of the intrinsic rate constants *β* and *β*′ of step 2; although *λ*_{3} is influenced by the equilibrium constant of the fast steps, it is not influenced by the rate or equilibrium constants of the slow steps.

Thus, the three-state model (Scheme 3) has 2 exponential processes with 2 apparent rate constants, *λ*_{2} and *λ*_{3} (*λ*_{2} > *λ*_{3}) (generally, an N-state model has (N−1) exponential processes). Each apparent rate constant is the weighted sum (linear combination) of the forward and reversal rate constants of the respective steps (Eqs. 32 and 37). The time courses *A*(*t*), *B*(*t*) and *C*(*t*) take the form described in Eqs. A21–A23 of Appendix 2. Step 3 occurs slowly, which is important for calculation of the turnover rate. In perturbation analysis, the rate constants of the fast steps (1 and 2) can be observed (Eqs. 32 and Eq. 37), whereas the rate constant of the slowest step (3) cannot be observed. However, because of the presence of in-series compliance, the forward rate constant of step 3 can be observed as discussed later in Section 19.

In Scheme 3, *t*→∞ will result in the steady state, in which *dA*(*t*)/*dt* = 0, *dB*(*t*)/*dt* = 0, and *dC*(*t*)/*dt* = 0 are achieved. From these conditions, it follows that the steady-state concentrations of *A*_{1}, *B*_{1}, *C*_{1} (still assuming that *γ*,*γ*′ *α*, *α*′, *β*,*β*′) are:

$${A}_{1}={A}_{\text{T}}{\alpha}^{\prime}{\beta}^{\prime}/M,\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{B}_{1}={A}_{\text{T}}\alpha {\beta}^{\prime}/M,\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{C}_{1}={A}_{\text{T}}\alpha \beta /M$$

(39)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}M\equiv \alpha \beta +\alpha {\beta}^{\prime}+{\alpha}^{\prime}{\beta}^{\prime}$$

The turnover rate *J* is calculated as:

$$J=\gamma {C}_{1}-{\gamma}^{\prime}{A}_{1}={A}_{\text{T}}(\alpha \beta \gamma -{\alpha}^{\prime}{\beta}^{\prime}{\gamma}^{\prime})/M$$

(40)

In the case of the cross-bridge cycle, *J* is the same as the ATP hydrolysis rate. The equilibrium is a special case (*J* = 0) of the steady state. If *αA*_{1} − *α*′*B*_{1} or *βB*_{1} − *β*′*C*_{1} is evaluated in place of Eq. 40, both become 0 because of the approximations we used. In this case, *B*_{1}/*A*_{1} = *α*/*α*′ = *K _{α}*, and

At steady state, the fast steps of the cyclic reaction can be approximated by the equilibrium, whereas the slowest (rate-limiting) step cannot be approximated by the equilibrium.

Based on Scheme 3, it is now possible to explain delayed tension without violating the Le Châtelier–Brown Principle. When the length of active muscle fiber is stretched (*δl* > 0) rapidly, force increases simultaneously because of the combined effect of attached cross-bridge elasticity and series compliance (phase 1, Figs. 1A and and2).2). The stretch simultaneously increases the elastic energy *δW* (mechanical energy) stored in the cross-bridge state **A** (Eq. 16); hence it either reduces the activation energy *E _{α}* (Eq. 17; Fig. 4A) or elevates the activation energy
${E}_{\alpha}^{\prime}$ (Fig. 4B). This increases the rate constant

$$\frac{1}{{B}_{1}}{\left(\frac{\partial {B}_{1}}{\partial \alpha}\right)}_{{\alpha}^{\prime}}=\frac{{\alpha}^{\prime}}{\alpha \lambda}$$

(41)

$$\text{and}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\frac{1}{{B}_{1}}{\left(\frac{\partial {B}_{1}}{\partial {\alpha}^{\prime}}\right)}_{\alpha}=-\frac{1}{\lambda}$$

(42)

In case 1 (Section 4), a side-wise multiplication of Eq. 20 to Eq. 41 yields Eq. 43.

$$\begin{array}{l}\frac{\delta {B}_{1}}{{B}_{1}}=\frac{{\alpha}^{\prime}\phi}{\lambda kT}\delta l\\ =\frac{1\times 2\phantom{\rule{0.16667em}{0ex}}\text{pN}}{2\times 1.38\times {10}^{-23}{\text{JK}}^{-1}\times 293\phantom{\rule{0.16667em}{0ex}}\text{K}}\times 0.75\phantom{\rule{0.16667em}{0ex}}\text{nm}=0.19\end{array}$$

(43)

In case 2 (Section 4), a side-wise multiplication of Eq. 21 by Eq. 42 likewise yields Eq. 43. Here we assumed that unitary force () is 2 pN (range 1–6 pN: Finer et al. 1994; Miyata et al. 1995; Molloy et al. 1995; Ishijima et al. 1996; Guilford et al. 1997), the value of *α* is comparable to that of *α*′, *T* = 293 K (20°C), and the length change is 0.1% stretch. Consequently, *δl* = 0.001 × 1.25 μm × 60% = 0.75 nm, where half sarcomere length is 1.25 μm, at which skinned psoas fiber experiments are usually carried out. 60% is to account for cross-bridge compliance, and the remaining 40% series compliance; the results with respect to the temperature effect on stiffness (rabbit psoas fibers) fit well if the series compliance is 40 ± 5% during standard activation at 20°C (Kawai 2003). The reported series compliance varies significantly, the lowest value being 19% (Bagni et al. 1990) and the highest 69% (Wakabayashi et al. 1994), and mid-range values having been published as well (44% in Higuchi et al. 1995 and 50–60% in Dobbie et al. 1998). Equation 43 suggests that the number of cross-bridges in state **B** increases as much as 19% in response to a stretch as small as 0.1%. The increase becomes about 38% for a 0.2% stretch, and so on. Here it is possible that a “ratchet” mechanism is at work in step 1 (i.e., preferred detachment of unstrained, or perhaps negatively strained, cross-bridges: Huxley 1957; Vale and Oosawa 1990; Smith and Geeves 1995; Smith 1998), increasing the number of cross-bridges in state **B** even further. However, it must be recognized that the ratchet mechanism is against the Le Châtelier-Brown principle, thereby the ratchet has its value.

The next thing that happens is an increase in the slower step 2 (**B** → **C**), because its rate is *βB*, and the number (*B*) of cross-bridges in state **B** is now larger than steady-state (*δB*_{1} > 0: Eq. 43). The stretch also increases the elastic energy in state **C**, which in turn reduces the activation energy
${E}_{\beta}^{\prime}$. Hence *β*′ is also accelerated (or *E _{β}* is elevated and

A mirror image of the above events occurs on release (*δl* < 0) of the active muscle fibers. On release, force drops (phase 1, Figs. 1B and and3B)3B) due to simultaneous release of the combined elasticity of series compliance and the attached cross-bridges. This results in a decrease in the elastic energy stored in state **A** (Eq. 16), which in turn increases *E _{α}* (Eq. 17) and reduces

In all cases, the transitions satisfy the Le Châtelier–Brown Principle. From the above discussion, it is clear that phase 2 represents **A** ↔ **B** interconversion (step 1), with the fast apparent rate constant of *λ*_{2}, and that phase 3 represents **B** ↔ **C** interconversion (step 2), with the medium apparent rate constant *λ*_{3} as indicated in Scheme 3. What it is important to notice here is that both phase 2 of stretch and phase 2 of release arise from the same elementary step of the cross-bridge cycle, because tension transients on stretch and release are almost in mirror image to each other (Figs. 1 and and3).3). Similarly, phase 3 of stretch and phase 3 of release arise from the same elementary step. Due to energetic considerations (Eqs. 16–21), however, the resulting rate constants for stretch and release are asymmetric. This is referred to as “non-linearity” of the rate constants (Huxley and Simmons 1971), and accounts for the fact that tension transients are not in perfect mirror image on stretch and release, as shown in Fig. 1A and 1B. From the above discussions, it can be concluded that tension develops (increases) by two mechanisms: through the fast **B** → **A** transition (phase 2) that occurs on release (counter-clockwise cycle in Scheme 3; Fig. 1B), and through the medium-speed **B** → **C** transition (phase 3) that occurs on stretch (clockwise cycle; Figs. 2 and and3A3A).

The three-state model (Scheme 3) can explain phase 2 and phase 3 of force transients without violating the Le Châtelier–Brown Principle, provided that step 3 is the rate-limiting step, the middle state (**B**) is a low-or no-force state, and there is a ratchet mechanism to increase the **A** → **B** transition or to diminish the **C** → **B** transition in response to stretch. Both states **A** and **C** must bear force, and step 1 must be faster than step 2.

In performing perturbation experiments in muscle fibers and probing the elementary steps of the cross-bridge cycle, it is important to keep the length change small. We typically use a sinusoidal length change with a 0.25% peak-to-peak amplitude, which corresponds to 1.9 nm (or ±0.95 nm) at the cross-bridge level when 40% series compliance is considered. Because this is smaller than the step size (4–11 nm: Molloy et al. 1995; Guilford et al. 1997; Kitamura et al. 1999), the elementary steps can be probed. A 10% length change, by contrast, corresponds to 75 nm at the cross-bridge level, and the cross-bridge would have to cycle many (7 to 19) times to compensate for this large length change. The time course of such an experiment would be limited by the slowest step of the cycle. As is seen in Eq. 43, even a small length change has a strong effect on cross-bridge occupancy, and therefore results in a large change in force.

The next question is whether, in Scheme 3, cross-bridges cycle in the clockwise direction or in the counter-clockwise direction to consume ATP and to perform work. To answer this, we need to look into the effects the substrate (MgATP) and product (phosphate and MgADP) concentrations have on the apparent rate constants.

In rabbit psoas fibers, an increase in the the MgATP concentration was shown to accelerate the rate constants of phases 2 and 3 in sinusoidal length-change experiments (Kawai 1978). In sinusoidal analysis, the apparent rate constant of phase 2 is termed 2*πc*, and that of phase 3 is termed 2*πb* (see Scheme 3). Because the temporal resolving power (signal-to-noise ratio) is greater in sinusoidal length-change experiments than in step length-change experiments, we adapted this technique. According to conclusions 2 and 3 above, the effect of ATP implies that the MgATP binding is associated with a fast reaction step whose rate constant is comparable to 2*πc*. We have also found that an increase in the Pi concentration accelerates the rate constant of phase 3 (2*πb*), but that Pi does not significantly affect the rate constant of phase 2 (2*πc*) (Kawai 1986). In light of conclusion 3, this observation implies that the Pi binding step is coupled to a medium-speed step with a rate constant comparable to 2*πb*. If we fit these results to Scheme 3, then it follows that ATP binding is associated with step 1, whereas Pi binding is associated with step 2. These possibilities are considered quantitatively in the following Sections 11–13.

To indicate a cross-bridge state we now use **X*** _{i}* (i = 0, 1, 2, …). When

The math that corresponds to Scheme 5 is identical to that of Scheme 4 with the following substitutions: **X**_{1} in place of **A**, **X**_{2} in place of **B**, **X**_{3} in place of **C**, *K*_{1}*S* in place *K _{α}*,

$${\lambda}_{2}=\frac{{K}_{1}S}{1+{K}_{1}S}{k}_{2}+{k}_{-2}$$

(44)

Thus, *λ*_{2} increases and saturates hyperbolically with an increase in S. The data of 2*πc vs. S* that were obtained from rabbit psoas fibers were fitted to Eq. 44 to derive three kinetic constants (*K*_{1} = 1.4 mM^{−1}, *k*_{2} = 440 s^{−1}, *k*_{−2} = 100 s^{−1} at 20°C: Kawai and Halvorson 1989) that characterize Scheme 5. Based on caged ATP experiments, Goldman et al. (1984) reported that the second order rate constant of cross-bridge detachment is at least 500 mM^{−1} s^{−1} (rabbit psoas, 20–22°C). This number is equivalent to *K*_{1}*k*_{2} (initial slope of Eq. 44) of Scheme 5, and is in approximate agreement with the value calculated from it (1.4 mM^{−1} × 440 s^{−1} = 620 mM^{−1} s^{−1}). The NPE-caged ATP that was used did not allow faster rate constants to be probed, in part because of the slow rate of photolysis (~100 s^{−1}) (Goldman et al. 1984) and also because of its competition with ATP (Thirlwell et al. 1995). Hence, the saturation phase of Eq. 44 cannot be defined based on caged ATP experiments. Sinusoidal analysis does not have these limitations. Hence Eq. 44 has been utilized extensively to characterize other fasttwitch (Galler et al. 2005) and slow-twitch (Wang and Kawai 1996) skeletal muscle fibers, and porcine (Zhao and Kawai 1996) and bovine (Fujita et al. 2002) myocardium as well as *Drosophila melanogaster* indirect flight muscles (Swank et al. 2006), as summarized in Table 2. Solution studies of S1 from rabbit back muscles (performed at 20–23°C) yielded *K*_{1} = 4.5 mM^{−1}, *k*_{2} = 400 s^{−1}, and *k*_{−2}<0.02 s^{−1} (Bagshaw et al. 1974), and similar studies of actoS1 (performed at 20°C) yielded *K*_{1}*k*_{2} = 1500 mM^{−1} s^{−1} (Millar and Geeves 1988). Similarly, studies of myofibril suspension (at 20°C) yielded *K*_{1}*k*_{2} = 1800 mM^{−1} s^{−1} and *k*_{−2}~0 (Herrmann et al. 1994). Thus, it can be concluded that *K*_{1} and *k*_{2} were in approximate agreement in four very different preparations.

From the above discussion, it is clear that **X**_{1} and **X**_{2} correspond to state **A** in Scheme 3, and that **X**_{3} corresponds to state **B**. This assignment cannot be transposed because that would require that the slow rate constant, *λ*_{3}(= 2*πb*), decreases as *S* increases, yet 2*πb* was found to increase concomitant with an increase in tHe MgATP concentration (Kawai 1978; Kawai and Zhao 1993). This assignment demonstrates that a cross-bridge cycles in the clockwise direction in Scheme 3. With regard to the ratchet mechanism discussed in Section 9, we have experimental evidence that *K*_{1} (ATP binding) becomes larger as isometric tension is decreased by various experimental manipulations (Zhao et al. 1996), indicating that cross-bridges detach more readily as load is reduced. This mechanism is consistent with the Fenn effect (Fenn 1923).

It has been shown that an increase in the MgADP concentration results in a decrease in the rate constants of phase 2 and phase 3 (Kawai 1986; Kawai and Halvorson 1989). This result can be interpreted as the effect of a competitive inhibitor (Scheme 6).

where **D** is MgADP, *K*_{0} is the association constant for **D** to **X**_{1}, and **X**_{0} is the ADP-bound state. It is known that myosin S1 has only one nucleotide-binding site (Rayment et al. 1993), and that it can bind either MgATP or MgADP. In Scheme 6, *K*_{0} is written underneath the reaction, because the association is the reversal reaction, from **X**_{1} to **X**_{0}.

Because **X**_{0}, **X**_{1} and **X**_{2} lose and gain as a group in Scheme 6, Eq. 45 can be set up.

$$d({X}_{0}+{X}_{1}+{X}_{2})/dt=-{k}_{2}{X}_{2}+{k}_{-2}{X}_{3}$$

(45)

From the mass action law,

$${K}_{0}={X}_{0}/{X}_{1}D\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\text{and}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{K}_{1}={X}_{2}/{X}_{1}S$$

(46)

where *D* is the MgADP concentration. Because *X _{i}* is the probability,

$${X}_{0}+{X}_{1}+{X}_{2}+{X}_{3}={X}_{T}$$

(47)

where *X _{T}* is a constant of time while we make the observation at the high speed. By eliminating

$$d{X}_{3}/dt+{\lambda}_{2}{X}_{3}=\epsilon {k}_{2}$$

(48)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\epsilon \equiv {K}_{1}S/(1+{K}_{0}D+{K}_{1}S)$$

(49)

$$\text{and}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{\lambda}_{2}\equiv \epsilon {k}_{2}+{k}_{-2}=\frac{{K}_{1}S}{1+{K}_{0}D+{K}_{1}S}{k}_{2}+{k}_{-2}$$

(50)

From Eq. 48 (Appendix 1), we get:

$$\begin{array}{l}{X}_{3}(t)={X}_{30}exp(-{\lambda}_{2}t)+{X}_{31},\\ {X}_{2}(t)=-\epsilon {X}_{30}exp(-{\lambda}_{2}t)+{X}_{21},\phantom{\rule{0.38889em}{0ex}}\text{etc}.\end{array}$$

(51)

where *X*_{30} is an integration constant and *X*_{31} is the steady state probability of X_{3}:

$${X}_{31}\equiv \epsilon {k}_{2}/{\lambda}_{2}=\frac{{K}_{1}S{K}_{2}}{1+{K}_{0}D+{K}_{1}S(1+{K}_{2})}\equiv \sigma $$

(52)

Other steady state probabilities are:

$${X}_{21}=\sigma /{K}_{2},\phantom{\rule{0.38889em}{0ex}}{X}_{11}=\sigma /{K}_{1}S{K}_{2},\phantom{\rule{0.38889em}{0ex}}{X}_{01}=\sigma {K}_{0}D/{K}_{1}S{K}_{2}$$

Therefore, Scheme 6 includes an exponential process with an apparent rate constant *λ*_{2}, as defined by Eq. 50. *λ*_{2} decreases hyperbolically as *D* is increased. Fitting the data from the ADP study in rabbit psoas fibers at 20°C to Eq. 50 led us to deduce that *K*_{0} = 2.8 mM^{−1} (Kawai and Halvorson 1989) (Table 2). In comparison, solution studies of S1 performed at 20–23°C yielded *K*_{0} = 3.7 mM^{−1} (Bagshaw et al. 1974). *K*_{0} is least important among all kinetic constants, however, because *D* is minimal (0.01–0.02 mM) in the presence of CP/CK (Meyers et al. 1985). Hence *X*_{0}< 0.01 and few cross-bridges are distributed in the AMD state.

The above discussion makes it clear that **X**_{0}, **X**_{1} and **X**_{2} correspond to state **A** in Scheme 3, and that**X**_{3} corresponds to state **B** in Scheme 3. Specifically, **X**_{1} corresponds to AM (A = actin. M = myosin), to which no nucleotide is bound. **X**_{2} corresponds to AM*ATP, and **X**_{0} to AM.ADP. **X**_{3} corresponds to both the weakly attached AM.ATP and the detached M.ATP. Recently, Palmer et al. (2007) showed a cross-bridge model in which 2*πc* corresponds to the detachment step (**X**_{2} ↔ **X**_{3} interconversion), based on the stochastic attachment and detachment cycle of cross-bridges.

Assuming that a cross-bridge state **X**_{4} makes a conformational change to generate state **X**_{5}, then the phosphate (**P**) is released to result in state **X**_{6} (Scheme 7):

If this reaction is viewed in the opposite direction, it follows that **P** forms a collision complex with **X**_{6} to result in **X**_{5}, and that **X**_{5} makes a slow conformational change to result in **X**_{4}. Mathematically this situation is exactly the same as that described in Scheme 5, with the following substitutions (but not correlations): **X**_{6} in place of **X**_{1}, **X**_{5} in place of **X**_{2}, **X**_{4} in place of **X**_{3}, **P** in place of **S**, *K*_{5} in place of *K*_{1}, *k*_{−4} in place of *k*_{2}, and *k*_{4} in place of *k*_{−2}. Therefore, Eq. 44 allows us to derive the apparent rate constant of Scheme 7 as:

$${\lambda}_{3}={k}_{4}+\frac{{K}_{5}P}{1+{K}_{5}P}{k}_{-4}$$

(53)

where *P* is the Pi concentration, and *K*_{5} is the association constant for **P** to **X**_{6}. In Scheme 7, *K*_{5} is written underneath the reaction, because the Pi association is the reversal reaction, from **X**_{6} to **X**_{5}.

The above discussion makes it clear that **X**_{4} corresponds to state **B**, and **X**_{5} and **X**_{6} to state **C**, of Scheme 3. This assignment cannot be transposed, because from that it would follow that an increase in the Pi concentration results in an increase in force (state **C**), which is opposite to what was observed (Cooke and Pate 1985; Kawai 1986). A transposition of this assignment would also contradict to the previous assignment of the ATP binding step (cross-bridges cycle in the clockwise direction in Scheme 3), because Pi must be released from the cycle. Consequently, **X**_{4} corresponds to the detached M.ADP.Pi state and the weakly attached AM.ADP.Pi state; **X**_{5} corresponds to the strongly attached AM*ADP.Pi state, and **X**_{6} corresponds to the AM*ADP state. It must be pointed out here that the AM*ADP state (**X**_{6}) is completely distinct from the AM.ADP state (**X**_{0}), which is consistent to solution studies (Sleep and Hutton 1980).

When examining the Pi effect, there is also a fast reaction to the left of Scheme 7 (step 1 in Scheme 3) that cannot be ignored, based on conclusion 3 above. The analysis is carried out in the same way as in Scheme 4, to result in the apparent rate constant *λ*_{3} (Kawai and Halvorson 1991):

$${\lambda}_{3}=\sigma {k}_{4}+\frac{{K}_{5}P}{1+{K}_{5}P}{k}_{-4}$$

(54)

*σ* is defined in Eq. 52 (0<*σ*<1). Thus, the apparent rate constant *λ*_{3} (2*πb*) increases and saturates with an increase in either *S* or *P*, and it decreases with an increase in *D*. As seen in Eq. 53 or Eq. 54, *λ*_{3} is a hyperbolic function of *P*, and fitting the data to this hyperbolic function makes it possible to deduce three kinetic constants (*k*_{4}, *k*_{−4}, *K*_{5}). *σ* of Eq. 52 is obtained from the ATP study, where the *K*_{0}*D* term can be dropped because *D* is very small (0.01–0.02 mM) in the presence of CP/CK. Pi is thought to be released via a narrow deep cleft (back door) in S1 (Yount et al. 1995); hence it is not likely that ADP competitively binds to the Pi binding site. An increase in the ADP concentration decreases the rate constant 2*πb* (*λ*_{3}), but this is via a reduction of *σ* in Eqs. 52 and 54.

Equations 53 or 54 has been used extensively to analyze cross-bridge kinetics in rabbit skeletal muscle fibers (Fortune et al. 1991; Kawai and Halvorson 1991; Dantzig et al. 1992; Millar and Homsher 1992; Walker et al. 1992; Wang and Kawai 1997; Galler et al. 2005). With sinusoidal analysis, *k*_{4} = 56 s^{−1}, *k*_{−4} = 129 s^{−1}, and *K*_{5} = 0.069 mM^{−1} were deduced for rabbit psoas fibers at 20°C (Kawai and Halvorson 1991). These values compare with those from caged Pi experiments (*k*_{4} = 79.2 s^{−1}, *k*_{−4} = 114.7 s^{−1}, *K*_{5} = 0.27 mM^{−1}; Dantzig et al. 1992, 20°C) as well as pressure release experiments (*k*_{4} ~ 15 s^{−1}, *k*_{−4} ~ 36 s^{−1}, *K*_{5} = 0.255 mM^{−1} at 12°C: deduced from Fortune et al. 1991; see note b of Table 2 legend), and are summarized in Table 2. If the data of the pressure-release experiment are extrapolated to 20°C using the Q_{10} values reported for *k*_{4} (6.8), *k*_{−4} (1.6), and *K*_{5} (1.1) (Zhao and Kawai 1994), then *k*_{4} = 69 s^{−1}, *k*_{−4} = 52 s^{−1}, and *K*_{5} = 0.275 mM^{−1} will result. If the data of the caged Pi experiment by Walker et al. (1992) (Table 2) are extrapolated to 20°C using the same Q_{10} values, then *k*_{4} = 70 s^{−1}, *k*_{−4} = 145 s^{−1}, and *K*_{5} = 0.172 mM^{−1}. As seen with these numbers, approximate agreement was achieved based on three different techniques from four laboratories. *K*_{5} deduced from sinusoidal analysis is somewhat (2.5–4 ×) smaller, primarily because Fortune et al. (1991), Dantzig et al. (1992), and Walker et al. (1992) fitted the rate constant results to Eq. 53, whereas Kawai and Halvorson (1991) fitted the Pi dependence data of both rate constant (Eq. 54) and amplitude simultaneously. In slow-twitch fibers of rabbit soleus, Millar and Homsher (1992) found: *k*_{4} = 1.96 s^{−1} and the second order Pi binding constant *K*_{5}*k*_{−4} = 1.994 mM^{−1} s^{−1} when using caged Pi at 20°C. These numbers compare well to those of Wang and Kawai (1997) in Table 2: *k*_{4} = 5.7 s^{−1}, *k*_{−4} = 4.5 s^{−1}, *K*_{5} = 0.18 mM^{−1}, hence *K*_{5}*k*_{−4} = 0.81 mM^{−1} s^{−1}. *k*_{4} is usually larger in our case, because of the inclusion of *σ* (Eq. 54) that accounts for fast equilibria to the left of Scheme 7. Because Pi released by the photolysis of caged Pi is limited to 2–3 mM, it is difficult to probe the saturation phase using this technique (*K*_{5} = 0.18 mM^{−1} indicates that half saturation occurs at 5.5 mM Pi). Although Eq. 54 demonstrates a sensitivity to *S* through *σ*, this equation has not been used to deduce *K*_{1} or *K*_{2} because of uncertainties about the ATP cleavage step (step 3), which skinned-fiber studies have not been able to resolve (see Section 14).

When the effects of ATP, ADP and Pi are considered, the only model consistent to the Le Châtelier–Brown Principle involves assigning phase 2 to steps surrounding the ATP binding step, and assigning phase 3 to steps surrounding the Pi release step. Consequently, cross-bridges cycle in the clockwise direction in Scheme 3, consuming ATP, generating force, and releasing Pi and ADP. Therefore, the transient force increase in phase 2 that was seen by Huxley and Simmons (1971) following step-length release is not the result of energy transduction (that requires ATP), but of a transient accumulation of the number of cross-bridges in the **A** state (which includes both the AM and AM*ATP species). Phase 2 was not correctly assigned because Huxley and Simmons (1971) used intact fibers, and hence did not examine the effects of ATP or Pi concentration.

Schemes 6 and and77 are combined to summarize the complete cross-bridge Scheme 8.

Note that the steps are renumbered. These are called “elementary steps” or “fundamental steps” (Gutfreund 1995). In this scheme, states **X**_{3} and **X**_{4} are equated and represented as **X**_{34} for simplicity. In reality, **X**_{3} ↔ **X**_{4} (step 3) corresponds to the ATP cleavage step (Lymn and Taylor 1971; Bagshaw et al. 1974;Xu et al. 2003), which is difficult to resolve in skinned fiber studies as they depend on strongly attached cross-bridges. Therefore, cross-bridge states before and after the cleavage of ATP are merged together and referred to as the **Det** (detached) state. The **Det** state includes both the weakly attached states (AM.ATP,AM.ADP.Pi) and the truly detached states (M.ATP, M.ADP.Pi). Because weakly attached states are based on ionic interactions (Brenner et al. 1982), at 200 mMionic strength the majority of **Det** consists of truly detached states. An accounting of step 3 is one remaining problem that should be solved by skinned fiber experiments in the coming years. What we know already is that step 3 is not very much slower than 2*πb* (*λ*_{3}), because if it were slower, then the rate constant 2*πb* would not be sensitive to the ATP concentration. Therefore, step 3 cannot be rate-limiting for the entire cross-bridge cycle as once suggested (Stein et al. 1985) based on four-state model of Lymn and Taylor (1971). Based on solution studies, it is known that in fast-twitch muscles at 20°C, *k*_{3} and *k*_{−3} are on the order of 30–120 s^{−1} (Bagshaw et al. 1974; Xu et al. 2003). The fact that this is similar to the rate constants of step 4 complicates the problem.

Step 6 is the slowest step in the cycle, and can be considered the “rate-limiting” step. If *k*_{6} + *k*_{−6} is similar to or faster than *k*_{4} + *k*_{−4}, then it follows that an increase in [MgATP] reduces 2*πb* and an increase in [MgADP] increases 2*πb*, both of which are contrary to the experimental evidence (Kawai 1978, 1986). Therefore, *k*_{6} + *k*_{−6} < *k*_{4} + *k*_{−4}, which makes step 6 the slowest in the cross-bridge cycle. The product of step 6 can be either **X**_{0} (as shown in Scheme 8) or can be converted directly to **X**_{1} by ADP release. Experiments in solution carried out by Sleep and Hutton (1980) indicated that *K*_{6}~50 (i.e., *k*^{−}_{6} *k*_{6}), and that reversal of the **X**_{0} state to the **X**_{6} state is possible (shown in a dashed arrow in Scheme 8). That is, **X**_{0} is likely to be an intermediate product of the hydrolysis pathway, and a part of the cross-bridge cycle.

An examination of Scheme 8 reveals that collision complex formation is followed by a conformational change in three places: (1) step 1 (ATP binding) to step 2, (2) the reversal of step 5 (Pi binding) to step 4, and (3) the reversal of step 0 (ADP binding) to step 6. The conformational change presumably accompanies a change in the force on a cross-bridge and secures ligand binding. This is consistent with the suggestion from Huxley (1980), that chemical reaction (collision complex formation) and mechanical reaction (conformational change) alternate in the cross-bridge cycle. Because these sequential reactions are difficult to resolve, they are frequently merged together. Therefore, the “Pi release step” may include both steps 4 and 5. In fact, the strongly attached AM*ADP.Pi state was not recognized in solution studies. Hence, in solution systems the conformational change (step 4) must have been merged together and referred to the “Pi release step”. Similarly, the “ADP release step” may include both steps 6 and 0. It has been often stated that the ADP release step is slow, but what that really means is that the preceding conformational change (step 6) is slow. The ADP release step itself (step 0) is a reversal of collision complex formation and, hence, cannot be slow. Similarly, “ATP binding” may include steps 1 and 2.

For the steady-state calculations, step 0 through step 5 can be approximated by equilibrium, and *X _{i}* can be calculated from five equilibrium constants (

Solution studies have shown that the Pi release step accompanies a large reduction in free energy (Δ*G*°), amounting to ~1/2 the energy generated by ATP hydrolysis (White and Taylor 1976). Hence, this step is almost irreversible (Sleep and Hutton 1980) (note: Δ*G*° = −*RT* log *K*, where *R* is the gas constant and *K* is the equilibrium constant), requiring 10–100 M Pi for reversal (Taylor 1979). This is why White and Taylor (1976) suggested that the Pi release step is associated with the force generation step. Skinned-fiber studies, in contrast, demonstrated that even at mM levels, Pi can diminish isometric tension (Cooke and Pate 1985; Kawai 1986), which indicates that the Pi release step is reversible. This reversibility has been demonstrated through ^{18}O exchange experiments (Webb et al. 1986), and by equilibrium constant measurements showing that *K*_{4} = 0.43 (step 4), and 1/*K*_{5}*P* 1.8 (step 5 at 8 mM Pi) based on values listed in Table 2 for rabbit psoas fibers (Kawai and Halvorson 1991). Therefore,
$\mathrm{\Delta}{G}_{4}^{\circ}=+2.0\phantom{\rule{0.16667em}{0ex}}\text{kJ}/\text{mole}$ and
$\mathrm{\Delta}{G}_{5}^{\circ}=-1.5\phantom{\rule{0.16667em}{0ex}}\text{kJ}/\text{mole}$, and their absolute values are less than *RT* (= 2.44 kJ/mole), hence there is not much change in the free energy in steps 4 or 5. It is likely that in skinned fibers, the free energy of ATP hydrolysis is stored in the contractile apparatus as elastic (mechanical) energy, which can be sensed as force, and can be readily reversed to chemical energy as the Pi concentration is elevated in the mM range. In the solution system, in contrast, a great deal of the liberated free energy is lost as heat. Hence, such a reaction cannot be reversed.

If we focus on one half sarcomere, isometric tension in this structure can be represented by Eq. 55 (Kawai and Zhao 1993). Because all cross-bridges are arranged mechanically in parallel (Fig. 6), overall tension is their sum of the tension of each.

A mechanical equivalent of one half sarcomere. Cross-bridges are arranged in parallel, and this assembly is connected to a series compliance (q). The parallel stiffness (*S*_{P}) is added to the entire structure. Although each element is written as elastic, **...**

$$(\text{Tension})={T}_{0}{X}_{0}+{T}_{1}{X}_{1}+{T}_{2}{X}_{2}+{T}_{5}{X}_{5}+{T}_{6}{X}_{6}$$

(55)

where *T _{i}* is tension when all cross-bridges are in state

Force generated (or supported) by each cross-bridge state. □ Rabbit psoas fibers (replotted from Kawai and Zhao 1993). ■ Rabbit soleus slow twitch fibers (plotted from data published by Wang and Kawai 1997). The data are normalized to **...**

**X**_{5} (AM*ADP.Pi) is a collision complex between **P** and **X**_{6} (AM*ADP). Tension supported (or generated) by **X**_{5} and **X**_{6} must be the same, i.e., *T*_{5} = *T*_{6} because the conformation of a protein (or proteins) cannot change instantly as the collision complex is formed. If the conformation is the same, the tension must be the same. Because *T*_{34} = 0, force must be generated at step 4, with the transition of **X**_{34} to **X**_{5}. Experimental evidence in support of this insight comes from an evaluation of force as the function of the Pi concentration (Fig. 5) (Kawai and Halvorson 1991; Kawai and Zhao 1993; Wang and Kawai 1997; Ranatunga 1999; Tesi et al. 2002); force does not suddenly decrease as Pi is bound to cross-bridges (Dantzig et al. 1992; Takagi et al. 2004). If additional force is generated upon ADP release, as reported in the case of smooth muscles (Rosenfeld et al. 2001), then it follows that *T*_{6}<*T*_{0}, *T*_{0}<*T*_{1}, or *T*_{6}<*T*_{1}. However, there is no evidence for this in skeletal muscle fiber studies as shown in Fig. 5.

An examination of Fig. 5 shows that forces supported by AM.ADP and AM are not significantly different (*T*_{0}≈*T*_{1}). This observation demonstrates that AM.ADP is indeed the collision complex. However, force supported by AM*ATP is significantly less than that supported by AM, indicating that a conformational change may occur by the time AM*ATP is formed. This supposition is strengthened by the finding of yet another state AM^{†}ATP(**X**_{1}* _{b}*) between the AM and AM*ATP states, as discussed in Section 16 below. In this case, the tension data as a function of

In parallel to force generation at step 4 and the transduction of chemical energy into elastic (mechanical) energy, the series compliance (including cross-bridge compliance) is strained. This must be the point of the cross-bridge cycle at which a swing of the lever arm (Dominguez et al. 1998; Geeves and Holmes 1999; Cooke 2005) takes place. Once force is generated, the cross-bridge is stabilized in this conformation by Pi release to result in **X**_{6}. The contractile system is now ready to perform external work. For this, it is convenient that the next step (step 6) is slow, as the muscle has to transfer both the mechanical energy and the momentum to the load before/while lifting it. The momentum is calculated according to Eq. 56 and related to the force application *F*(*t*) and its duration Δ*t*:

$$(\text{momentum})=(\text{mass})\times (\text{velocity})={\int}_{\mathrm{\Delta}t}F(t)dt={F}_{av}\mathrm{\Delta}t$$

(56)

where *F _{av}* is the average force. Equation 56 represents a quantity known as the “force-time integral” among muscle physiologists. This equation shows that the momentum increases with increased duration (Δ

In some of our publications on rabbit psoas fibers (e.g., Kawai and Zhao 1993; Zhao and Kawai 1994) as well as soleus slow-twitch fibers (Wang and Kawai 1996, 1997) and on myocardium (Kawai et al. 1993; Zhao and Kawai 1996), an extra exponential process (termed process D) that occurs at higher frequency than process C (2*πd*>2*π c*) was described. Similarly, in the case of phase 2 of step analysis, Abbott and Steiger (1977), Coupland et al. (2001), and Davis et al. (2002) resolved two exponential processes in rabbit psoas fibers, and Ford et al. (1977) recognized four exponential processes in frog tibialis anterior fibers. Because of the ATP dependence of process D, we ascribed this process to an extra step that occurs immediately after ATP binding to the myosin head as in Scheme 9:

where **X**_{1}* _{b}* is the collision complex, step 1b is a conformational change, and step 2 is yet another conformational change.

Stiffness depends on the number of attached cross-bridges as well as on series compliance (Huxley et al. 1994; Wakabayashi et al. 1994; Kojima et al. 1994; Higuchi et al. 1995) and parallel stiffness (Fig. 6). Because these elements are intermixed and distributed, the formal mathematical treatment becomes extremely complex, as was shown by Mijailovich et al. (1996). The particular patterns in which cross-bridges are arrayed and linked to thin filaments add another degree of complexity with respect to determining the relationship between the force of individual cross-bridges and measured stiffness or tension (Tanner et al. 2007). For practical applications, however, we have previously proposed a simple model in which a parallel array of attached cross-bridges “y” is connected to a series compliance *q* (Fig. 6) (Wang and Kawai 1997; Kawai 2003), and the entire structure is connected to parallel stiffness *S _{P}*. Both series compliance and the parallel stiffness can have viscous components at the same time. Hence, in general,

$$(\text{Stiffness})={Y}_{\text{a}}{X}_{\text{a}}/(1+{Y}_{\text{a}}{X}_{\text{a}}q)+{S}_{\text{P}}$$

(57)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{X}_{\text{a}}\equiv {X}_{0}+{X}_{1}+{X}_{2}+{X}_{5}+{X}_{6}=1-{X}_{34}$$

and *Y*_{a} is the stiffness contribution from the cross-bridges when all cross-bridges are attached (*X*_{a} = 1). The contribution of the parallel stiffness can be removed by subtraction of the viscoelasticity of a muscle fiber when all cross-bridges are detached (when is open in Fig. 6). Classically, *S _{P}* would be the viscoelasticity of relaxed muscle fibers. However, it may be better to measure this quantity on fibers that are relaxed by using 40 mM 2,3- butanedione monoxime (BDM) (e.g., Fujita et al. 2002) in the presence of Ca

The formula for calculation of the ATP hydrolysis rate can be derived from Eq. 40 (Zhao and Kawai 1994):

$$J=(\text{ATPase})={A}_{\text{T}}({k}_{6}{X}_{6}-{k}_{-6}{X}_{0})\cong {A}_{\text{T}}{k}_{6}{X}_{6}$$

(58)

In practice, the *k*_{−6} term is dropped as shown, both because this reaction is highly unidirectional (*K*_{6}~50, hence *k*_{6} *k*_{−6}; Sleep and Hutton 1980) and due to the fact that *X*_{0} is very small (<0.01) under normal experimental conditions in which the ADP concentration is 0.01–0.02 mM in the presence of CP/CK. *A _{T}* is the total myosin S1 concentration, and it is 0.17–0.2 mM in psoas fibers (Yates and Greaser 1983; Tikunov et al. 2000), 0.13–0.14 mM in soleus fibers (Tikunov et al. 2000), and 0.15 mM in myocardium (Barsotti and Ferenczi 1988).

As discussed in Section 17, series compliance makes an important contribution on muscle mechanics (Huxley et al. 1994; Wakabayashi et al. 1994; Kojima et al. 1994; Higuchi et al. 1995; Wang et al. 1999; Kawai 2003). This is particularly evident at the onset of contraction, when force develops from a zero load to the maximum load. In the following, we define that *η* is the step size, *ν* is the number of steps that take place in 1 sec, and *ρ*′ is the stiffness of one half sarcomere, which includes the contribution from both series compliance and cross-bridges (Eq. 57). For simplicity, we assume that *S _{P}* = 0. Under isometric conditions, the increase in force (

$$dF={\rho}^{\prime}\phantom{\rule{0.16667em}{0ex}}\eta v\phantom{\rule{0.16667em}{0ex}}dt$$

(59)

In Eq. 59, *ην* is the distance in 1 sec to pull the thin filament toward the M-line (hence *ην* is the velocity of filament elongation), which at the same time stretches series compliance and other compliant elements of cross-bridges to result in a force increase by *ρ*′ (*ην*). Let us assume that developed force (*F*) affects *η* in the following way:

$$\eta ={\eta}_{0}(1-F/{F}_{0})$$

(60)

where *F*_{0} is the maximum isometric force. That is, a step of the cross-bridge becomes increasingly more difficult as force develops, and as discussed by Piazzesi et al. (2007). Equation 60 can be actually derived from the force–velocity relationship (see e.g., Woledge et al. 1985; Debold et al. 2005) by a linear extrapolation of the velocity at *F*_{0}. An alternative to Eq. 60 can be generated if we assume that:

$$v={v}_{0}(1-F/{F}_{0})$$

(61)

Equation 61 is an expression of the Fenn effect (Fenn 1923), i.e. the rate of stepping becomes increasingly more difficult as force develops. From Eq. 59 and Eq. 60 (or Eq. 61), we get:

$$dF/dt+{\lambda}_{4}F={\lambda}_{4}{F}_{0}$$

(62)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{\lambda}_{4}\equiv {\rho}^{\prime}{\eta}_{0}{v}_{0}/{F}_{0}$$

(63)

$$\text{Therefore},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{F}_{0}/{\rho}^{\prime}={\eta}_{0}{v}_{0}/{\lambda}_{4}$$

(64)

Thus, from Eq. 62 (Appendix 1) we arrive at:

$$F(t)={F}_{0}[1-exp(-{\lambda}_{4}t)]$$

(65)

This system has one exponential process with an apparent rate constant of *λ*_{4} and amplitude *F*_{0}. Note that *F*(0) = 0 is assumed. Equation 65 would be the form of a force time course in which a fiber is suddenly activated by Ca^{2+}, or in which the length of the active fiber is suddenly released to zero load and restretched to its original length (resulting in force redevelopment like that seen by *k _{tr}* measurement). The exponential function (Eq. 65) may be a good approximation for describing the time course data. In fact, available data are consistent with the exponential function (e.g. Regnier et al. 1995; Stehle et al. 2002; Piroddi et al. 2003).

Now, let us evaluate the right side of Eq. 64. With single myofibrils, Piroddi et al. (2003) observed *k*_{tr} = 8.0 s^{−1} and *k*_{act} = 7.9 s^{−1} in rabbit psoas at 15°C. These are the values of *λ*_{4}. 1/*ν*_{0} is the average time a cross-bridge spends for one cycle, and is related to the elementary steps by:

$$\frac{1}{{v}_{0}}=\frac{1}{{\lambda}_{2}}+\frac{1}{{\lambda}_{3}}+\frac{1}{{k}_{6}+{k}_{-6}}$$

(66)

In Eq. 66, each term on the right-hand side corresponds to a time constant: 1/*λ*_{2} is the time constant of step 2, which is the average time a cross-bridge spends to complete step 2; similarly, 1/*λ*_{3} to is the average time spent to complete step 4, and 1/(*k*_{6} + *k*_{−6}) is the average time spent to complete step 6. If *k*_{6} *λ*_{3} (note that: *k*_{−6} *k*_{6}, *λ*_{3}<*λ*_{2}), then from Eq. 66,

$${v}_{0}\cong {k}_{6}$$

(67)

*k*_{6} 16–18 s^{−1} was reported in rabbit psoas fibers at 15–20°C (Zhao and Kawai 1994); hence *ν*_{0}~17 s^{−1}. *η*_{0} = 5.3 nm was reported by Kitamura et al. (1999) for the step size (see also Molloy et al. 1995; Guilford et al. 1997). Thus, the right side of Eq. 64 is:

$${\eta}_{0}{v}_{0}/{\lambda}_{4}=5.3\phantom{\rule{0.16667em}{0ex}}\text{nm}\times 17\phantom{\rule{0.16667em}{0ex}}{\text{s}}^{-1}/7.95\phantom{\rule{0.16667em}{0ex}}{\text{s}}^{-1}=11\phantom{\rule{0.16667em}{0ex}}\text{nm}$$

The left side of Eq. 64 (*F*_{0}/*ρ*′) is the same as the (force):(stiffness) ratio, or the instantaneous length release needed to abolish full tension. The ratio (active tension):(elastic modulus) was measured to be T:Y_{∞} = 1.14% in rabbit psoas fibers at 15°C (Zhao and Kawai 1994). Therefore,

$$(\text{Force})/(\text{stiffness})=1.25\phantom{\rule{0.16667em}{0ex}}\mu \text{m}\times 0.0114=14\phantom{\rule{0.16667em}{0ex}}\text{nm}$$

where 1.25 μm is the half sarcomere length of the experiment. Thus, approximate agreement can be reached for both sides of Eq. 64, which may validate this equation and the theory behind it. Because *k*_{tr} is slowest of all rate constants measured in the cross-bridge cycle, this parameter has been thought to be limited by the slowest step in the cross-bridge cycle (Gordon et al. 2000). The above analysis provides support for this hypothesis.

The fact that *k*_{tr} is affected by Pi (Regnier et al. 1995; Tesi et al. 2000) is presumably related to the approximation used in Eq. 67, which can be improved by inclusion of the Pi effect (*λ*_{3} is affected by the Pi concentration: Eq. 54). Assuming that the 1/*λ*_{2} term can be ignored in Eq. 66 (because *λ*_{2} *k*_{6}), Eq. 66 can be rewritten by using Eq. 54 to result in Eq. 68.

$${v}_{0}={v}_{00}+\frac{\psi {K}_{5}P}{1+\xi {K}_{5}P}$$

(68)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}{v}_{00}\equiv \frac{\sigma {k}_{4}{k}_{6}}{\sigma {k}_{4}+{k}_{6}}={k}_{6}-\frac{{k}_{6}^{2}}{\sigma {k}_{4}+{k}_{6}}$$

(69)

$$\text{and}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\psi \equiv \frac{{k}_{-4}{{k}_{6}}^{2}}{{(\sigma {k}_{4}+{k}_{6})}^{2}},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\xi \equiv 1+\frac{{k}_{-4}}{\sigma {k}_{4}+{k}_{6}}$$

(70)

Equation 68 is a hyperbolically increasing and saturating function with respect to *P*, which is similar to Eqs. 53 and 54. Thus, an increase in *P* results in the increase in *λ*_{4} via Eq. 63 and as reported (Regnier et al. 1995; Tesi et al. 2000). Similarly, an increase in *S* results in the increase in *λ*_{4} via *S* → *σ* → *λ*_{3} → *ν*_{0} → *λ*_{4} (see Eqs. 52, 54, 66, 63, respectively).

Several colleagues in the field have suggested that exponential process A of sinusoidal analysis (and hence phase 4 of step analysis) may represent the slowest step of the cross-bridge cycle (step 3 in Scheme 3, and step 6 in Scheme 8). Until recently, we have discounted this possibility because the three-state model, as represented in Scheme 3, can have only two exponential processes (see Appendix 2 for an analysis of a general case) and, in principle, the rate constants of the slowest step cannot be detected by the perturbation analysis method. The same is true for Scheme 8. However, the introduction of series compliance changes this outlook significantly. As shown above, series compliance may introduce an extra exponential process with a slow rate constant, such as 2*πa* of process A, which is present in all three subtypes of fast-twitch muscle fibers examined (Galler et al. 2005). In contrast, the amplitude of process A(phase 4) is very small or even absent in the cases of insect indirect flight muscles (Fig. 3) (Thorson and White 1969; Abbott 1973; Pringle 1978; Marcussen and Kawai 1990; Swank et al. 2006), myocardium (Saeki et al. 1978, 1991; Zhao and Kawai 1996; Wannenburg et al. 2000; Fujita et al. 2002), and covalently but partially cross-linked rabbit psoas fibers (Tawada and Kawai 1990). In slow-twitch fibers, process A is present but its amplitude is small (Kawai and Schachat 1984; Wang and Kawai 1996, 1997). These observations suggest that the sarcomere structure is more rigid in these muscle preparations than in fast-twitch skeletal muscle fibers, minimizing the amplitude of process A. In fact, in insect muscle fibers, an extra C-filament is present to secure the connection between the thick filament and the Z-line, and to stabilize the sarcomere structure (White 1983). In cardiac muscles, there is a possibility that myosin binding protein C (MyBP-C) reaches the thin filament (Spirito et al. 1997; Squire et al. 2003), which may give extra stability to sarcomeres. In partially cross-linked rabbit psoas fibers, the extra covalent linkage (18% of cross-bridges) was created to stabilize the preparation, causing process A to disappear (Tawada and Kawai 1990).

Let us perform a calculation to examine the possibility that process A (phase 4) represents the slowest step of the cross-bridge cycle. Because 2*πa* = 6 s^{−1} was observed in rabbit psoas fibers at 20°C (Zhao and Kawai 1994), the right side of Eq. 64 becomes:

$${\eta}_{0}{v}_{0}/{\lambda}_{4}=5.3\phantom{\rule{0.16667em}{0ex}}\text{nm}\times 17\phantom{\rule{0.16667em}{0ex}}{\text{s}}^{-1}/6\phantom{\rule{0.16667em}{0ex}}{\text{s}}^{-1}=15\phantom{\rule{0.16667em}{0ex}}\text{nm}$$

Because T:Y_{∞} = 1.3% was observed in the same study (20°C), the left side of Eq. 64 becomes:

$$(\text{Force})/(\text{stiffness})=1.25\phantom{\rule{0.16667em}{0ex}}\mu \text{m}\times 0.013=16\phantom{\rule{0.16667em}{0ex}}\text{nm}$$

Thus, the values are consistent with process A being the slowest step of the cross-bridge cycle. This insight is, however, limited to preparations in which process A (phase 4) can be observed. In any case, this is a subject that will be exciting to pursue in future studies.

The apparent rate constant of tension development (*k*_{act}, *k*_{tr}) or of process A of sinusoidal analysis (phase 4 of step analysis) appears to be proportionate to the rate constant of the slowest step of the cross-bridge cycle (*k*_{6}), with the proportionality constant *ρ*′*η*_{0}/*F*_{0} (Eq. 63). This proportionality exists because of the presence of in-series compliance in sarcomeres.

In the discussion in Sections 7–18, the condition (*α* + *α*′ *β* + *β*′ *γ* + *γ*′) was used for the sake of simplicity. However, what would happen if this condition were not met, and if the conditions were merely *α* + *α*′ > *β* + *β*′ > *γ* + *γ*′ as suggested from the kinetic constants listed in Table 2? The answer to this question is simply that the approximation becomes less accurate as these numbers get closer to each other, and that the analysis may become qualitative rather than quantitative. That is, there will be an increasing contribution of steps 4–5 to the apparent rate constant *λ*_{2} of phase 2, with the Pi concentration consequently having increasingly more influence on phase 2. The conclusions derived then become qualitative, and their usefulness depends on one’s own expectations. They include the assignment of phase 2 and its rate constant *λ*_{2} to step 2 of Scheme 8, and the assignment of phase 3 and its rate constant *λ*_{3} to step 4 of the same scheme. If, however, one is not satisfied with this approximation, then an exact derivation can be undertaken according to Appendix 2. Alternatively, one could use the method used to analyze the data from experiments in myocardium (Section 21), as this does not rely on approximation.

In cardiac muscle fibers (myocardium), the apparent rate constants *λ*_{2} (2*π c*) and *λ*_{3} (2*π b*) are closer to one another than their counterparts in skeletal muscle fibers are, and thus the approximation approach is not as accurate as in skeletal fibers. In fact, we have found that the Pi concentration affects the apparent rate constant 2*πc* significantly. In this case, the sum (Eq. 71) and the product (Eq. 72) of the two apparent rate constants can be used to deduce the kinetic constants of the elementary steps. If we take this approach, there is no need for approximation except that we still have to assume that *k*_{6} is much smaller () than other rate constants. Under these conditions,

$${\lambda}_{2}+{\lambda}_{3}=\epsilon {k}_{2}+{k}_{-2}+{k}_{4}+\zeta {k}_{-4}$$

(71)

$${\lambda}_{2}{\lambda}_{3}=\epsilon {k}_{2}{k}_{4}+\epsilon \zeta {k}_{2}{k}_{-4}+\zeta {k}_{-2}{k}_{-4}$$

(72)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\zeta \equiv {K}_{5}P/(1+{K}_{5}P)$$

(73)

*ε* is defined in Eq. 49 [*ε* = *K*_{1}*S*/(1 + *K*_{0}*D* + *K*_{1}*S*)]. For the derivation of Eqs. 71 and 72, see Eqs. 14 and 15 of Kawai and Halvorson (1991). Both Eqs. 71 and 72 are hyperbolic functions of *S*, *P* and *D*, as are the previously discussed Eqs. 50 and 54. By fitting the data (effects of *S*, *P* and *D*) to these equations, we can deduce all rate constants (*k*_{2}, *k*_{−2}, *k*_{4}, *k*_{−4}) and association constants (*K*_{0}, *K*_{1}, *K*_{5}). This method enabled us to resolve the kinetic constants of the elementary steps in ferret myocardium (Kawai et al. 1993) and bovine myocardium (Lu et al. 2005).

Martin and Barsotti (1994) used caged ATP on guinea pig heart muscle, and deduced the second order ATP binding constant (our *K*_{1}*k*_{2}) to be ~39 mM^{−1} s^{−1} at 21°C. This number compares with the values we have calculated from Table 2: ferret myocardium 48 mM^{−1} s^{−1} (2°C, Kawai et al. 1993), porcine myocardium 138 mM^{−1} s^{−1} (20°C, Zhao and Kawai 1996), and bovine myocardium 240 mM^{−1} s^{−1} (25°C, Fujita et al. 2002). Araujo and Walker (1996) used caged Pi on rat ventricular myocytes, and found the second order Pi binding constant (our *K*_{5}*k*_{−4}) to be 3.1 mM^{−1} s^{−1} (15°C). This number compares well with our values calculated from Table 2 (same references): ferret myocardium 6.4 mM^{−1} s^{−1} (20°C), porcine myocardium 1.1 mM^{−1} s^{−1} (20°C), and bovine myocardium 1.8 mM^{−1} s^{−1} (25°C). Thus, even when very different techniques and different cardiac preparations are used, we are able to achieve approximate agreement.

As shown in this mini review, the three-state model is adequate to account for the three phases of force transients that are generated in response to a step length change. The six-state model is an extension of the three-state model that includes ligand (ATP,ADP, Pi)-bound states. In order to correlate each phase of a transient to individual elementary step of the cross-bridge cycle correctly, it is necessary to examine the effects of these ligand concentrations on the apparent rate constants and to fit their results to the model. Therefore, any future cross-bridge models must be able to predict the effects of these ligands, in particular those of ATP and Pi. When studying force transients, it is best not to make a large increase in force, because this may cause a delay in the time course to stretch series elastic elements. Since multiple cross-bridge cycles are needed to stretch series elastic elements, the time course may be limited by the slowest step (rate-limiting step) of the cross-bridge cycle (Section 19).

The authors are indebted to Dr. David W. Maughan of the University of Vermont (USA) for his critical reading of the manuscript and his constructive comments and useful suggestions. The authors are also grateful to Dr. Christine Blaumueller for her critical reading of the manuscript and creative suggestions. This work was supported by a grant from NIH HL70041 to MK. The contents of this work are solely the responsibility of the authors and do not necessarily represent the official view of the awarding organization.

The solution of the differential equation (A1) can be found in the following way.

$$dY(t)/dt+\lambda Y(t)=h$$

(A1)

Now we substitute *Y*(*t*) with *Z*(*t*) so that

$$dY(t)/dt=[-\lambda Z(t)+dZ(t)/dt]exp(-\lambda t)$$

(A3)

By substituting Eqs. A2 and A3 to Eq. A1, we arrive at

$$dZ(t)/dt=0$$

(A4)

By integrating Eq. A4, we get

$$Y(t)={Y}_{0}exp(-\lambda t)+h/\lambda $$

(A6)

Therefore, Eq. A1 has an exponential process with the rate constant *λ*, amplitude *Y*_{0}, and the steady state value *h/λ* (Eq. A6). *Y*_{0} is determined by the initial conditions.

In Scheme 3, the following differential equations can be set up.

$$dA/dt=-(\alpha +{\gamma}^{\prime})A+{\alpha}^{\prime}B+\gamma \text{C}$$

(A7)

$$dB/dt=\alpha A-(\beta +{\alpha}^{\prime})B+{\beta}^{\prime}\text{C}$$

(A8)

$$dC/dt={\gamma}^{\prime}A+\beta B-(\gamma +{\beta}^{\prime})C$$

(A9)

Eqs. A7–A9 can be written in matrix form.

$$\frac{dX}{dt}=-HX$$

(A10)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}X(t)\equiv \left(\begin{array}{c}A(t)\\ B(t)\\ C(t)\end{array}\right)$$

(A11)

$$\text{And}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}H\equiv \left(\begin{array}{ccc}(\alpha +{\gamma}^{\prime})& -{\alpha}^{\prime}& -\gamma \\ -\alpha & (\beta +{\alpha}^{\prime})& -{\beta}^{\prime}\\ -{\gamma}^{\prime}& -\beta & (\gamma +{\beta}^{\prime})\end{array}\right)$$

(A12)

$$\text{Also},\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}A=B+C={A}_{\text{T}}$$

(A13)

Note the similarity of Eqs. A1 and A10. Note also that |*H*| = det(*H*) = 0, which is consistent with Eq. A13. The eigen values (*λ*) of matrix *H* can be found by determining the roots of Eq. A14.

$$\begin{array}{rl}\mid H-I\lambda \mid & =\mid \begin{array}{ccc}(\alpha +{\gamma}^{\prime})-\lambda & -{\alpha}^{\prime}& -\gamma \\ -\alpha & (\beta +{\alpha}^{\prime})-\lambda & -{\beta}^{\prime}\\ -{\gamma}^{\prime}& -\beta & (\gamma +{\beta}^{\prime})-\lambda \end{array}\mid \\ & =0\end{array}$$

(A14)

where *I* is the identity matrix. From Eq. A14,

$${\lambda}^{3}-2Q{\lambda}^{2}+M\lambda =0$$

(A15)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}Q\equiv (\alpha +{\alpha}^{\prime}+\beta +{\beta}^{\prime}+\gamma +{\gamma}^{\prime})/2$$

and

$$M\equiv \alpha \beta +\beta \gamma +\gamma \alpha +\alpha {\beta}^{\prime}+\beta {\gamma}^{\prime}+\gamma {\alpha}^{\prime}+{\alpha}^{\prime}{\beta}^{\prime}+{\beta}^{\prime}{\gamma}^{\prime}+{\gamma}^{\prime}{\alpha}^{\prime}$$

(A16)

Equation A15 has three roots *λ*_{1}, *λ*_{2}, and *λ*_{3}.

$${\lambda}_{1}=0$$

$${\lambda}_{2}=Q+\sqrt{R}$$

(A17)

$${\lambda}_{3}=Q-\sqrt{R}=\frac{M}{Q+\sqrt{R}}$$

(A18)

$$\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}R\equiv {Q}^{2}-M$$

where *R* *Q*^{2} − *M*

With Eqs. A17 and A18, Eq. A10 can be solved to result in:

$$X(t)=UE(t)$$

(A19)

$$\begin{array}{l}\text{where}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}U\equiv \left(\begin{array}{ccc}{A}_{1}& {U}_{12}& {U}_{13}\\ {B}_{1}& {U}_{22}& {U}_{23}\\ {C}_{1}& {U}_{32}& {U}_{33}\end{array}\right),\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\phantom{\rule{0.38889em}{0ex}}\text{and}\\ E(t)\equiv \left(\begin{array}{c}1\\ exp(-{\lambda}_{2}t)\\ exp(-{\lambda}_{3}t)\end{array}\right)\end{array}$$

(A20)

The elements *U _{i}*

$${U}^{-1}HU=\left(\begin{array}{ccc}{\lambda}_{1}& 0& 0\\ 0& {\lambda}_{2}& 0\\ 0& 0& {\lambda}_{3}\end{array}\right)$$

From Eqs. A19 and A20, the individual solutions are:

$$A(t)={A}_{1}+{U}_{12}exp(-{\lambda}_{2}t)+{U}_{13}exp(-{\lambda}_{3}t)$$

(A21)

$$B(t)={B}_{1}+{U}_{22}exp(-{\lambda}_{2}t)+{U}_{23}exp(-{\lambda}_{3}t)$$

(A22)

$$C(t)={C}_{1}+{U}_{32}exp(-{\lambda}_{2}t)+{U}_{33}exp(-{\lambda}_{3}t)$$

(A23)

Therefore, according to Eq. 34, a force transient has two exponential processes that correspond to phases 2 and 3, with the rate constants λ_{2} and λ_{3}, respectively.

If R<0, then λ_{2}, λ_{3}, *U*_{i2} and *U*_{i3} are complex numbers, and λ_{3} = λ_{2}* (Eqs. A17 and A18), and *U*_{i3} = *U*_{i2}*, where * indicates the complex conjugate. In this case, Eq. A19 can be rearranged to result:

$$X(t)=\left(\begin{array}{ccc}{A}_{1}& 2\mathfrak{R}({U}_{12})& 2\mathfrak{J}({U}_{12})\\ {B}_{1}& 2\mathfrak{R}({U}_{22})& 2\mathfrak{J}({U}_{22})\\ {C}_{1}& 2\mathfrak{R}({U}_{32})& 2\mathfrak{J}({U}_{32})\end{array}\right)\phantom{\rule{0.16667em}{0ex}}\left(\begin{array}{c}1\\ exp(-Qt)cos\sqrt{-R}\xb7t\\ exp(-Qt)sin\sqrt{-R}\xb7t\end{array}\right)$$

(A24)

where (*U _{ij}* refers to the real, and (

By setting *dA*/*dt* = 0 and *dB*/*dt* = 0 in Eqs. A7 and A8, and with the constraints of Eq. A13, we can solve the steady-state concentrations of *A*_{1}, *B*_{1} and *C*_{1}. These three values constitute an eigen vector (column vector) which belongs to the eigen value λ_{1} *M* is defined in Eq. A16.

$${A}_{1}={A}_{T}(\beta \gamma +\gamma {\alpha}^{\prime}+{\alpha}^{\prime}{\beta}^{\prime})/M$$

(A25)

$${B}_{1}={A}_{T}(\gamma \alpha +\alpha {\beta}^{\prime}+{\beta}^{\prime}{\gamma}^{\prime})/M$$

(A26)

$${C}_{1}={A}_{\text{T}}(\alpha \beta +\beta {\gamma}^{\prime}+{\gamma}^{\prime}{\alpha}^{\prime})/M$$

(A27)

The turnover (ATP hydrolysis) rate is:

$$\begin{array}{l}J=\alpha {A}_{1}-{\alpha}^{\prime}{B}_{1}=\beta {B}_{1}-{\beta}^{\prime}{C}_{1}=\gamma {C}_{1}-{\gamma}^{\prime}{A}_{1}\\ ={A}_{T}(\alpha \beta \gamma -{\alpha}^{\prime}{\beta}^{\prime}{\gamma}^{\prime})/M\end{array}$$

(A28)

A similar analysis based on Scheme 8 is found in the Appendix of Kawai (2003).

Masataka Kawai, Department of Anatomy and Cell Biology, University of Iowa, Iowa City, IA 52245, USA.

Herbert R. Halvorson, 4661 Devonshire Rd., Detroit, MI 48224, USA.

- Abbott RH. An interpretation of the effect of fiber length and calcium on the mechanical properties of insect flight muscle. Cold Spring Harb Symp Quant Biol. 1973;37:647–654.
- Abbott RH, Steiger GJ. Temperature and amplitude dependence of tension transients in glycerinated skeletal and insect fibrillar muscle. J Physiol. 1977;266:13–42. [PubMed]
- Araujo A, Walker JW. Phosphate release and force generation in cardiac myocytes investigated with caged phosphate and caged calcium. Biophys J. 1996;70:2316–2326. [PubMed]
- Bagni MA, Cecchi G, Colomo F, Poggesi C. Tension and stiffness of frog muscle fibres at full filament overlap. J Muscle Res Cell Motil. 1990;11:371–377. [PubMed]
- Bagshaw CR, Eccleston JF, Eckstein F, Goody RS, Gutfreund H, Trentham DR. The magnesium ion-dependent adenosine triphosphatase of myosin. Two-step processes of adenosine triphosphate association and adenosine diphosphate dissociation. Biochem J. 1974;141:351–364. [PubMed]
- Barsotti RJ, Ferenczi MA. Kinetics of ATP hydrolysis and tension production in skinned cardiac muscle of the guinea pig. J Biol Chem. 1988;263:16750–16756. [PubMed]
- Bershitsky SY, Tsaturyan AK. Tension responses to joule temperature jump in skinned rabbit muscle fibres. J Physiol. 1992;447:425–448. [PubMed]
- Brenner B. Effect of Ca
^{2+}on cross-bridge turnover kinetics in skinned single rabbit psoas fibers: implications for regulation of muscle contraction. Proc Natl Acad Sci USA. 1988;85:3265–3269. [PubMed] - Brenner B, Schoenberg M, Chalovich JM, Greene LE, Eisenberg E. Evidence for cross-bridge attachment in relaxed muscle at low ionic strength. Proc Natl Acad Sci USA. 1982;79:7288–7291. [PubMed]
- Burton K, Simmons RM, Sleep J. Kinetics of force recovery following length changes in active skinned single fibres from rabbit psoas muscle. J Physiol. 2006;573(2):305–328. [PubMed]
- Campbell KB, Razumova MV, Kirkpatrick RD, Slinker BK. Nonlinear myofilament regulatory processes affect frequency-dependent muscle fiber stiffness. Biophys J. 2001;81:2278–2296. [PubMed]
- Cooke R. Milestone in physiology: the sliding filament model: 1972–2004. J Gen Physiol. 2005;123:643–656. [PMC free article] [PubMed]
- Cooke R, Pate E. The effects of ADP and phosphate on the contraction of muscle fibers. Biophys J. 1985;48:789–798. [PubMed]
- Coupland ME, Puchert E, Ranatunga KW. Temperature dependence of active tension in mammalian (rabbit psoas) muscle fibres: effect of inorganic phosphate. J Physiol. 2001;536:879–891. [PubMed]
- Dantzig J, Goldman Y, Millar NC, Lacktis J, Homsher E. Reversal of the cross-bridge force-generating transition by the photogeneration of phosphate in rabbit psoas muscle fibers. J Physiol. 1992;451:247–278. [PubMed]
- Davis JS, Rodgers ME. Force generation and temperature-jump and length-jump tension transients in muscle fibers. Biophys J. 1995;68:2032–2040. [PubMed]
- Davis JS, Satorius CL, Epstein ND. Kinetic effects of myosin regulatory light chain phosphorylation on skeletal muscle contraction. Biophys J. 2002;83:359–370. [PubMed]
- Dickinson MH, Hyatt CJ, Lehmann F-O, Moore JR, Reedy MC, Simcox A, Tohtong R, Vigoreaux JO, Yamashita H, Maughan DW. Phosphorylation-dependent power output of transgenic flies: an integrated study. Biophys J. 1997;73:3122–3134. [PubMed]
- Dobbie I, Linari M, Piazzesi G, Reconditi M, Koubassova N, Ferenczi M, Lombardi V, Irving M. Elastic bending and active tilting of myosin heads during muscle contraction. Nature. 1998;396:383–387. [PubMed]
- Dominguez R, Freyzon Y, Trybus KM, Cohen C. Crystal structure of a vertebrate smooth muscle myosin motor domain and its complex with the essential light chain: visualisation of the pre-power stroke state. Cell. 1998;94:559–571. [PubMed]
- Fenn WO. A quantitative comparison between the energy liberated and the work performed by the isolated sartorius muscle of the frog. J Physiol. 1923;58:175–203. [PubMed]
- Finer JT, Simmons RM, Spudich JA. Single myosin mechanics: piconewton forces and nanometre steps. Nature (Lond) 1994;368:113–119. [PubMed]
- Ford LE, Huxley AF, Simmons RM. Tension responses to sudden length change in stimulated frog muscle fibres near slack length. J Physiol. 1977;269:441–515. [PubMed]
- Fortune NS, Geeves MA, Ranatunga KW. Tension responses to rapid pressure release in glycerinated rabbit muscle fibers. Proc Natl Acad Sci USA. 1991;88:7323–7327. [PubMed]
- Fujita H, Sasaki D, Ishiwata S, Kawai M. Elementary steps of the cross-bridge cycle in bovine myocardium with and without regulatory proteins. Biophys J. 2002;82:915–928. [PubMed]
- Galler S, Wang BG, Kawai M. Elementary steps of the cross-bridge cycle in fast-twitch fiber types from rabbit skeletal muscles. Biophys J. 2005;89:3248–3260. [PubMed]
- Geeves MA, Holmes KC. Structural mechanism of muscle contraction. Annu Rev Biochem. 1999;68:687–728. [PubMed]
- Geeves MA, Holmes KC. The molecular mechanism of muscle contraction. Adv Protein Chem. 2005;71:161–193. [PubMed]
- Goldman YE, Hibberd MG, Trentham DR. Relaxation of rabbit psoas muscle fibres from rigor by photochemical generation of adenosine-5′-triphosphate. J Physiol. 1984;354:577–604. [PubMed]
- Gordon AM, Homsher E, Regnier M. Regulation of contraction in striated muscle. Physiol Rev. 2000;80:853–924. [PubMed]
- Guilford WH, Dupuis DE, Kennedy G, Wu J, Patlak JB, Warshaw DM. Smooth muscle and skeletal muscle myosins produce similar unitary forces and displacements in the laser trap. Biophys J. 1997;72:1006–1021. [PubMed]
- Gutfreund H. Receptors, transmitters and catalysts. Cambridge University Press; 1995. Kinetics for the life sciences.
- Heinl P, Kuhn HJ, Rüegg JC. Tension responses to quick length changes of glycerinated skeletal muscle fibres from the frog and tortoise. J Physiol. 1974;237:243–258. [PubMed]
- Herrmann C, Lionne C, Travers F, Barman T. Correlation of ActoS1, myofibrillar, and muscle fiber ATPases. Biochemistry. 1994;33:4148–4154. [PubMed]
- Higuchi H, Yanagida T, Goldman YE. Compliance of thin filaments in skinned fibers of rabbit skeletal muscle. Biophys J. 1995;69:1000–1010. [PubMed]
- Homsher W, Millar NC. Caged compounds and striated muscle contraction. Annu Rev Physiol. 1990;52:875–896. [PubMed]
- Huxley AF. Muscle structure and theories of contraction. Prog Biophys Biophys Chem. 1957;7:255–318. [PubMed]
- Huxley AF. Muscular contraction. J Physiol. 1974;243:1–43. [PubMed]
- Huxley AF. Reflections on muscle. Princeton University Press; Princeton: 1980.
- Huxley AF, Simmons RM. Proposed mechanism of force generation in striated muscle. Nature. 1971;233:533–538. [PubMed]
- Huxley HE, Stewart A, Sosa H, Irving T. X-ray diffraction measurements of the extensibility of actin and myosin filaments in contracting muscle. Biophys J. 1994;67:2411–2421. [PubMed]
- Ishijima A, Kojima H, Higuchi H, Harada Y, Funatsu T, Yanagida T. Multiple- and single-molecule analysis of the actomyosin motor by nanometer-piconewton manipulation with a microneedle: unitary steps and force. Biophys J. 1996;70:383–400. [PubMed]
- Julian FJ, Sollins KR, Sollins MR. A model for the transient and steady-state mechanical behavior of contracting muscle. Biophys J. 1974;14:546–562. [PubMed]
- Kawai M. Head rotation or dissociation? A study of exponential rate processes in chemically skinned rabbit muscle fibers when MgATP concentration is changed. Biophys J. 1978;22:97–103. [PubMed]
- Kawai M. The role of orthophosphate in crossbridge kinetics in chemically skinned rabbit psoas fibers as detected with sinusoidal and step length alterations. J Muscle Res Cell Motil. 1986;7:421–434. [PubMed]
- Kawai M. What do we learn by studying the temperature effect on isometric tension and tension transients in mammalian striated muscle fibres? J Muscle Res Cell Motil. 2003;24:127–138. [PubMed]
- Kawai M, Brandt PW. Sinusoidal analysis: a high resolution method for correlating biochemical reactions with physiological processes in activated skeletal muscles of rabbit, frog and crayfish. J Muscle Res Cell Motil. 1980;1:279–303. [PubMed]
- Kawai M, Halvorson HR. Role of MgATP and MgADP in the crossbridge kinetics in chemically skinned rabbit psoas fibers. Study of a fast exponential process (C) Biophys J. 1989;55:595–603. [PubMed]
- Kawai M, Halvorson HR. Two step mechanism of phosphate release and the mechanism of force generation in chemically skinned rabbit psoas muscle. Biophys J. 1991;59:329–342. [PubMed]
- Kawai M, Schachat FH. Differences in the transient response to fast and slow skeletal muscle fibers: correlations between complex modulus and myosin light chains. Biophys J. 1984;45:1145–1151. [PubMed]
- Kawai M, Zhao Y. Cross-bridge scheme and force per cross-bridge state in skinned rabbit psoas muscle fibers. Biophys J. 1993;65:638–651. [PubMed]
- Kawai M, Saeki Y, Zhao Y. Cross-bridge scheme and the kinetic constants of elementary steps deduced from chemically skinned papillary and trabecular muscles of the ferret. Circ Res. 1993;73:35–50. [PubMed]
- Kirkwood JG, Oppenheim I. Chemical thermodynamics. McGraw-Hill Book Co. Inc; New York: 1961.
- Kitamura K, Tokunaga M, Iwane AH, Yanagida T. A single myosin head moves along an actin filament with regular steps of 5.3 nanometres. Nature. 1999;397(6715):129–134. [PubMed]
- Kojima H, Ishijima A, Yanagida T. Direct measurement of stiffness of single actin filaments with and without tropomyosin by in vitro nanomanipulation. Proc Natl Acad Sci USA. 1994;91:12962–12966. [PubMed]
- Labeit D, Watanabe K, Witt C, Fujita H, Wu Y, Lahmers S, Funck T, Labeit S, Granzier H. Calcium-dependent molecular spring elements in the giant protein titin. Proc Natl Acad Sci. 2003;100:13716–13721. [PubMed]
- Lu X, Bryant MK, Bryan KE, Rubenstein PA, Kawai M. Role of the N-terminal negative charge of actin in cross-bridge kinetics and force generation in reconstituted bovine myocardium. J Physiol. 2005;564(1):65–82. [PubMed]
- Lymn RW, Taylor EW. Mechanism of adenosine triphosphate hydrolysis by actomyosin. Biochemistry. 1971;10:4617–4624. [PubMed]
- Marcussen BL, Kawai M. Role of MgATP and inorganic phosphate ions in crossbridge kinetics in insect (
*Lethocerus colossicus*) flight muscle. In: Sperelakis N, Wood JD, editors. Frontiers in smooth muscle research. Wiley-Liss; New York: 1990. pp. 805–813. [PubMed] - Martin H, Barsotti RJ. Relaxation from rigor of skinned trabeculae of the guinea pig induced by laser photolysis of caged ATP. Biophys J. 1994;66:1115–1128. [PubMed]
- Meyer RA, Brown TR, Kushmerick MJ. Phosphorus nuclear magnetic resonance of fast- and slow-twitch muscle. Am J Physiol. 1985;248:C279–C287. [PubMed]
- Mijailovich SM, Fredberg JJ, Butler JP. On the theory of muscle contraction: filament extensibility and the development of isometric force and stiffness. Biophys J. 1996;71:1475–1484. [PubMed]
- Millar NC, Geeves MA. Protein fluorescence changes associated with ATP and adenosine 5′-[γ-thio]triphosphate binding to skeletal muscle myosin subfragment 1 and actomyosin subfragment 1. Biochem J. 1988;249:735–743. [PubMed]
- Millar NC, Homsher E. Kinetics of force generation and phosphate release in skinned rabbit soleus muscle fibers. Am J Physiol. 1992;252 [PubMed]

Cell Physiol. 31:C1239–C1245. [PubMed] - Miyata H, Yoshikawa H, Hakozaki H, Suzuki N, Furuno T, Ikegami A, Kinosita K, Jr, Nishizaka T, Ishiwata S. Mechanical measurements of single actomyosin motor force. Biophys J. 1995;68:286S–290S. [PubMed]
- Molloy JE, Burns JE, Kendrick-Jones J, Tregear RT, White DCS. Movement and force produced by a single myosin head. Nature (Lond) 1995;378:209–212. [PubMed]
- Murase M, Tanaka H, Nishiyama K, Shimizu H. A three-state model for oscillation in muscle: sinusoidal analysis. J Muscle Res Cell Motil. 1986;7:2–10. [PubMed]
- Palmer BM, Suzuki T, Wang Y, Barnes WD, Miller MS, Maughan DW. Two-state model of acto-myosin attachment-detachment predicts C-process of sinusoidal analysis. Biophys J. 2007;93:760–769. [PubMed]
- Piazzesi G, Reconditi M, Linari M, Lucii L, Bianco P, Brunello E, Decostre V, Stewart A, Gore DB, Irving TC, Irving M, Lombardi V. Skeletal muscle performance determined by modulation of number of myosin motors rather than motor force or stroke size. Cell. 2007;131:784–795. [PubMed]
- Piroddi N, Tesi C, Pellegrino MA, Tobacman LS, Homsher E, Poggesi C. Contractile effects of the exchange of cardiac troponin for fast skeletal troponin in rabbit psoas single myofibrils. J Physiol. 2003;552.3:917–931. [PubMed]
- Preston LC, Lipscomb S, Robinson P, Mogensen J, McKenna WJ, Watkins H, Ashley CC, Redwood CS. Functional effects of the DCM mutant Gly159Asp troponin C in skinned muscle fibres. Pflugers Arch. 2007;453:771–776. [PubMed]
- Pringle JWS. Stretch activation of muscle: function and mechanism. Proc R Soc Lond B. 1978;201:107–130. [PubMed]
- Ranatunga KW. Effects of inorganic phosphate on endothermic force generation in muscle. Proc Biol Sci. 1999;266:1381–1385. [PMC free article] [PubMed]
- Rayment I, Rypniewski WR, Schmidt-Bäse K, Smith R, Tomchick DR, Benning MM, Winkelmann DA, Wesenberg G, Holden HM. Three-dimensional structure of myosin subfragment-1: a molecular motor. Science. 1993;261:50–58. [PubMed]
- Regnier M, Morris C, Homsher E. Regulation of the cross-bridge transition from a weakly to strongly bound state in skinned rabbit muscle fibers. Am J Physiol. 1995;269(6 Pt 1):C1532–C1539. [PubMed]
- Rosenfeld SS, Xing J, Whitaker M, Cheung HC, Brown F, Wells A, Milligan RA, Sweeney HL. Kinetic and spectroscopic evidence for three actomyosin: ADP states in smooth muscle. J Biol Chem. 2001;275:25418–25426. [PubMed]
- Saeki Y, Kawai M, Zhao Y. Comparison of crossbridge dynamics between intact and skinned myocardium from ferret right ventricles. Circ Res. 1991;68:772–781. [PubMed]
- Saeki Y, Sagawa K, Suga H. Dynamic stiffness of cat heart muscle in Ba
^{2+}-induced contracture. Circ Res. 1978;42:324–333. [PubMed] - Sleep JA, Hutton RL. Exchange between inorganic phosphate and adenosine 5′-triphosphate in the medium by actomyosin subfragment 1. Biochemistry. 1980;19:1276–1283. [PubMed]
- Smith DA. A strain-dependent ratchet model for [phosphate]- and [ATP]-dependent muscle contraction. J Muscle Res Cell Motil. 1998;19:189–211. [PubMed]
- Smith DA, Geeves MA. Strain-dependent cross-bridge cycle for muscle. Biophys J. 1995;69:524–537. [PubMed]
- Spirito P, Seidman CE, McKenna WJ, Maron BJ. The management of hypertrophic cardiomyopathy. N Engl J Med. 1997;336:775–785. [PubMed]
- Squire JM, Luther PK, Knupp C. Structural evidence for the interaction of C-protein (MyBP-C) with actin and sequence identification of a possible actin-binding domain. J Mol Biol. 2003;331:713–724. [PubMed]
- Stehle R, Krüger M, Scherer P, Brixius K, Schwinger RH, Pfitzer G. Isometric force kinetics upon rapid activation and relaxation of mouse, guinea pig and human heart muscle studied on the subcellular myofibrillar level. Basic Res Cardiol. 2002;97(Suppl 1):I127–I135. [PubMed]
- Steiger GJ, Abbott RH. Biochemical interpretation of tension transients produced by a four-state mechanical model. J Muscle Res Cell Motil. 1981;2:245–260. [PubMed]
- Stein LA, Greene LE, Chock PB, Eisenberg E. Rate-limiting step in the actomyosin adenosinetriphosphatase cycle: studies with myosin subfragment 1 cross-linked to actin. Biochemistry. 1985;24:1357–1363. [PubMed]
- Swank DM, Vishnudas VK, Maughan DW. An exceptionally fast actomyosin reaction powers insect flight muscle. Proc Natl Acad Sci USA. 2006;103:17543–17547. [PubMed]
- Takagi Y, Shuman H, Goldman YE. Coupling between phosphate release and force generation in muscle actomyosin. Philos Trans R Soc Lond B Biol Sci. 2004;359:1913–1920. [PMC free article] [PubMed]
- Tanner BC, Daniel TL, Regnier M. Sarcomere lattice geometry influences cooperative myosin binding in muscle. PLoS Comput Biol. 2007;3:1195–1211. [PMC free article] [PubMed]
- Tawada K, Kawai M. Covalent cross-linking of single muscle fibers from rabbit psoas increases oscillatory power. Biophys J. 1990;57:643–647. [PubMed]
- Taylor EW. Mechanism of actomyosin ATPase and the problem of muscle contraction. CRC Crit Rev Biochem. 1979;6:103–164. [PubMed]
- Tesi C, Colomo F, Nencini S, Piroddi N, Poggesi C. The effect of inorganic phosphate on force generation in single myofibrils from rabbit skeletal muscle. Biophys J. 2000;78:3081–3092. [PubMed]
- Tesi C, Colomo F, Piroddi N, Poggesi C. Characterization of the cross-bridge force-generating step using inorganic phosphate and BDM in myofibrils from rabbit skeletal muscles. J Physiol. 2002;541(1):189–199. [PubMed]
- Thirlwell H, Sleep JA, Ferenczi MA. Inhibition of unloaded shortening velocity in permeabilized muscle fibres by caged ATP compounds. J Muscle Res Cell Motil. 1995;16:131–137. [PubMed]
- Thorson J, White DC. Distributed representations for actinmyosin interaction in the oscillatory contraction of muscle. Biophys J. 1969;9:360–390. [PubMed]
- Thorson J, White DC. Role of cross-bridge distortion in the small-signal mechanical dynamics of insect and rabbit striated muscle. J Physiol. 1983;343:59–84. [PubMed]
- Tikunov BA, Sweeney HL, Rome LC. Quantitative electrophoretic analysis of myosin heavy chains in single muscle fibers. J Appl Physiol. 2000;90:1927–1935. [PubMed]
- Vale RD, Oosawa F. Protein motors and Maxwell’s demons: does mechanochemical transduction involve a thermal ratchet? Adv Biophys. 1990;26:97–134. [PubMed]
- Wakabayashi K, Sugimoto Y, Tanaka H, Ueno Y, Takezawa Y, Amemiya Y. X-ray diffraction evidence for the extensibility of actin and myosin filaments during muscle contraction. Biophys J. 1994;67:2422–2435. [PubMed]
- Walker JW, Lu Z, Moss RL. Effects of Ca
^{2+}on the kinetics of phosphate release in skeletal muscle. J Biol Chem. 1992;267:2459–2466. [PubMed] - Wang G, Kawai M. Effects of MgATP and MgADP on the cross-bridge kinetics of rabbit soleus slow-twitch muscle fibers. Biophys J. 1996;71:1450–1461. [PubMed]
- Wang G, Kawai M. Force generation and phosphate release steps in skinned rabbit soleus slow-twitch muscle fibers. Biophys J. 1997;73:878–894. [PubMed]
- Wang G, Ding W, Kawai M. Does thin filament compliance diminish the cross-bridge kinetics? A study in rabbit psoas fibers. Biophys J. 1999;76:978–984. [PubMed]
- Wannenburg T, Heijne GH, Geerdink JH, Van Den Dool HW, Janssen PM, De Tombe PP. Cross-bridge kinetics in rat myocardium: effect of sarcomere length and calcium activation. Am J Physiol. 2000;279:H779–H790. [PubMed]
- Webb MR, Hibberd MG, Goldman YE, Trentham DR. 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. 1986;261:15557–15564. [PubMed]
- White DC. The elasticity of relaxed insect fibrillar flight muscle. J Physiol. 1983;343:31–57. [PubMed]
- White HD, Taylor EW. Energetics and mechanism of actomyosin adenosine triphosphatase. Biochemistry. 1976;15:5818–5826. [PubMed]
- Woledge RC, Curtin NA, Homsher E. Energetic aspects of muscle contraction (Monographs of the Physiological Society, No. 41) Academic Press; New York: 1985. [PubMed]
- Xu S, Offer G, Gu J, White HD, Yu LC. Temperature and ligand dependence of conformation and helical order in myosin filaments. Biochemistry. 2003;42:390–401. [PubMed]
- Yates LD, Greaser ML. Quantitative determination of myosin and actin in rabbit skeletal muscle. J Mol Biol. 1983;168:123–141. [PubMed]
- Yount RG, Lawson D, Rayment I. Is myosin a “back door” enzyme? Biophys J. 1995;68:44S–47S. [PubMed]
- Zhao Y, Kawai M. The effect of the lattice spacing change on cross-bridge kinetics in chemically skinned rabbit psoas muscle fibers. II. Elementary steps affected by the spacing change. Biophys J. 1993;64:197–210. [PubMed]
- Zhao Y, Kawai M. Kinetic and thermodynamic studies of the cross-bridge cycle in rabbit psoas muscle fibers. Biophys J. 1994;67:1655–1668. [PubMed]
- Zhao Y, Kawai M. Inotropic agent EMD 53998 weakens nucleotide and phosphate binding to cross bridges in porcine myocardium. Am J Physiol. 1996;271 [PubMed]

Heart Circ Physiol. 40:H1394–H1406. - Zhao Y, Swamy PMG, Humphries KA, Kawai M. The effect of partial extraction of troponin C on the elementary steps of the cross-bridge cycle in rabbit psoas fibers. Biophys J. 1996;71:2759–2773. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |