| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Cavendish Laboratory, Cambridge CB3 0HE, United Kingdom
Correspondence: Address reprint requests to Andrej Vilfan, J. Stefan Institute, Jamova 39, 1000 Ljubljana, Slovenia. Tel.: 386-1-4773-900; Fax: 386-1-4263269; E-mail: andrej.vilfan{at}ijs.si.
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
In an experiment to determine the isometric transient response, a muscle fiber is held at both ends to prevent it from contracting. The muscle is then suddenly shortened (or stretched) by a fixed amount, and the tension T that it generates is measured. Immediately after the imposed change of length, the tension shifts from the isometric value T0 to a new value, which is termed T1. But shortly afterwards (typically within 2 ms), the tension adjusts to a new value, termed T2. Subsequently, it gradually reverts to the original isometric value T0, and the entire transient response is usually completed in a fraction of a second. It is generally accepted that the initial response T1 corresponds to the mechanical deformation of cross-bridges and provides a direct measure of their elasticity (Huxley and Simmons, 1971
). The interpretation of T2 is rather more controversial. It is often attributed to force generation by the working stroke of bound myosin molecules (Huxley and Simmons, 1971
; Hill, 1974
; Eisenberg et al., 1980
; Huxley and Tideswell, 1996
; Brenner et al., 1995
; Duke, 1999
, 2000
) and this interpretation has recently gained support from x-ray interference techniques applied to shortening fibers (Irving et al., 2000
; Piazzesi et al., 2002b
). But alternative models suggest that the force regeneration might be due, in part, to the rapid binding of new myosin heads to the thin filament (Brenner, 1991
; Howard, 2001
), or that it might involve the activation of the second myosin head (Huxley and Tideswell, 1997
).
In this article, we wish to address a fundamental problem connected with the interpretation of force transients. The present theories are all based on the consideration of a single pair of filaments, i.e., one filament containing myosin molecules, interacting with one actin filament. The dynamics of this filament pair is generalized to that of a whole muscle fiber by assuming that all filament pairs in a fiber behave in exactly the same way. This assumption is certainly justified as long as there are no static or dynamic instabilities in the system. However, the possibility of such instabilities has been known for a long time (Hill, 1974
). Moreover, a stochastic model of the actomyosin cycle, based on the swinging lever-arm hypothesis, has shown that instabilities do arise when values of parameters such as the lever-arm displacement and the cross-bridge elasticity are chosen to provide effective energy transduction (Duke, 1999
, 2000
). Such instabilities would give rise to a region of negative slope in the T2 curve of a single filament pair. Several reasons have been advanced for the absence of any negative slope in the experimentally determined T2 curve. Huxley and Simmons (1971)
argued that the power stroke is subdivided into several small steps, and fixed the step size so that the T2 curve had zero slope for limitingly small changes of length. In the model proposed by Eisenberg et al. (1980)
, the flatness of the T2 curve was explained by a broad distribution of cross-bridge strain after attachment, combined with a specific strain dependence of the transition rates to ensure the proper occupancies of the two bound states. A further explanation involved the compliance of the filaments and the distribution of binding sites on the thin filament in addition to a subdivided power stroke (Huxley and Tideswell, 1996
). Duke (1999)
has suggested that the flat T2 curve of a muscle fiber can arise despite an instability in the dynamics of a single pair of filaments, owing to the symmetry of a sarcomere. We investigate this possibility further in this article.
| METHODS |
|---|
|
|
|---|
We assumed complete mechanical relaxation of the system in each step, i.e., the strain of all elastic elements is equilibrated before the next transition takes place. The structure involving sarcomeres, filaments, and myosin heads was described as a circuit of elements with given resting lengths and compliances. The strain of every cross-bridge was calculated, given the constraint of fixed total length of the system (isometric conditions), or of fixed force acting on the ends of the system (isotonic conditions).
The T2 transients were always determined 2 ms after the stretch/release. All other parameter values are summarized in Table 1.
|
| MOTORS ACTING BETWEEN A SINGLE PAIR OF FILAMENTS |
|---|
|
|
|---|
, associated with ADP release. In the ATP-rich physiological environment, this step is quickly followed by detachment of the head and we therefore skip the transitional state (A3) in our model. The attachment and detachment rates, which are summarized in Table 1, determine the shape of the force-velocity relation, but have little influence on the transient response to length steps, which mainly depend on the power-stroke displacement d and the elastic constant K of the myosin cross-bridge. We assume that the transitions A1
A2 take place on a faster timescale than the detachment and reattachment of heads. Therefore, once a head has bound with strain
in the state A1, its state can be described as a statistical ensemble of the states A1 and A2 with probabilities given by the Boltzmann factors
![]() | (1) |
) =
-P2(
), where
Gstroke is the free energy change associated with phosphate release.
|
(
) (normalized to 1), then the isometric force that they generate is
![]() | (2) |
x is suddenly applied to the pair of filaments, the force will first change to
![]() | (3) |
+
x). As a result the force adjusts to
![]() | (4) |
As a first approximation, we can neglect the distribution of strains and assume
= 0 for all myosin heads, which allows us to calculate the T2 curves analytically. Their shape depends on two dimensionless parameters,
![]() | (5) |
![]() | (6) |
A2). The parmeter
0 measures this energy relative to the thermal energy, while
1 compares it to the magnitude of the chemical free-energy change
Gstroke that accompanies the power-stroke transition. Note that if
1 < 1, the power stroke can occur immediately after the myosin head binds to the thin filament; but if
1 > 1, the conformational change is energetically inhibited initially, and it is only likely to occur once the thin filament has been moved forward by the action of other motors. On grounds of efficiency of energy transduction, we expect the value of
1 to be as high as possible, while ensuring that there is sufficient chemical energy to drive the power stroke, i.e.,
1
1 (Duke, 1999
50% (Kushmerick and Davies, 1969
GATP
20 kBT, which leads to
0
10. Taking into account that not all energy stored in the spring can be converted to mechanical work one obtains a better estimate
0
20 (Duke, 1999
1
1 if most of the energy of hydrolysis is used to power the stroke
Gstroke
GATP). This corresponds to an elastic constant K = 2.5 pN/nm if a power-stroke distance of d = 8 nm is assumed.
The form of the T2 curve depends critically on the values of these dimensionless parameters. If
0 > 2, an interval with a negative slope (a hysteresis) exists (Hill, 1974
). This means that under isotonic (constant load) conditions the system can be bistable. The location of the interval of negative slope depends on the value of
1. If the hysteresis spans the origin, then the state with
x = 0 will be unstable under isotonic conditions. Putting
we see that this occurs if
![]() | (7) |
1
1 the T2 curve always has a negative slope around the stationary point. For example, with
0 = 20, Eq. 7 yields 0.85 <
1 < 1.22. However, it should be borne in mind that this calculation does not take into account the distribution of strains on cross-bridges within the ensemble, and so the actual range might deviate slightly from this estimate. An example of the T2 transient for a group of myosin heads between two firmly clamped filaments is shown in Fig. 2.
|
|
1 = 1.33, dashed line) the difference between this T2 curve and the one in the strictly isometric state is only quantitative. The most important result is that the range of negative slope in T2 remains. On the other hand, in the oscillating case (
1 = 0.8, dot-dashed line) the negative slope is flattened out because of the broader distribution of strains on cross-bridges.
Effect of filament elasticity
The compliance of the thick and thin filaments is also expected to affect the T2 curves. We model it by introducing a linear elasticity in the backbone connecting the heads (which is, to first order, equivalent to an elasticity in the track the motors are running on) (Vilfan et al., 1998
). With
being the linear modulus of the filaments, which has the measured value
= 44,000 pN (Kojima et al., 1994
), the spring constant of a filament segment of length L is Kf = L-1
. In (Vilfan et al., 1998
) a linear two-state model was used to study the effect of filament elasticity on force-velocity relations, showing that they reduce the isometric force if the compliance of a filament segment between two bound motors becomes comparable to the motor spring constant K. Thin filaments are just stiff enough to prevent significant losses due to this effect. However, the filament compliance can have a significantly bigger effect on the transient response. The result of a simulation is shown in Fig. 4. Although the filament compliance almost halves the total stiffness of the system, it leaves the major part of the interval of negative slope unaffected.
|
38 nm (7 monomers). We therefore propose that the binding rate to a particular site is proportional to the Boltzmann factor determined by the spring distortion energy, composed of a longitudinal and a lateral component
![]() | (8) |
= x - (i + 7j)c and
where j denotes the repeat on the helix and i the consecutive number of the actin monomer. The total binding rate at a given position x then reads
![]() | (9) |
These curves have been measured experimentally using S1 myosin heads by Steffen et al. (2001)
and the data were fitted with KA = 15 pN/nm, the value which we use here. The position dependence of the binding rate for these parameters is shown in Fig. 5. A similar distribution has also been measured for myosin V (Veigel et al., 2002
).
|
| MOTORS IN MUSCLE SARCOMERES |
|---|
|
|
|---|
|
|
|
|
|
|
1, the model shows an interesting feature. In the isometric state the number of cross-bridges in the state A2 can be very low, although the generated force per attached cross-bridge is as high as 5 pN, in agreement with single-filament measurements (Kawai et al., 2000
Isotonic response
Another important class of experiments which provides information on the actomyosin interation involves the isotonic transient. Here the applied force is initially set at the value of the stalling load T0, so that the fiber is prevented from contracting. The force is then suddenly changed and held constant at a different value, while the length of the fiber is recorded. Some early experiments showed that after a small step change of load, damped oscillations were imposed on the steady contraction or extension of the fiber (Podolsky, 1960
; Granzier et al., 1990
). Such oscillations are particularly clear in recent experiments on single muscle fibers (Edman and Curtin, 2001
).
One possible cause of this oscillatory response has previously been suggested on the basis of the stochastic model of the actomyosin interaction used in this article. When
1 > 1, the chemical cycles of myosin motors on the same filament can become synchronized at loads close to the stalling force (Duke, 1999
, 2000
). A pair of filaments then slides in a step-wise fashion under isotonic conditions. But during steady shortening, the motors on different filaments within the same muscle fiber operate out of phase, so that there is no macroscopic manifestation of the steps. However, an abrupt change in the load can cause the synchronization of a large fraction of the bound motors, whereupon the steps do become observable (Duke, 1999
). Because the correlation of the motors soon decays, the macroscopic steps fade and a damped oscillation is seen.
A much stronger oscillatory response is seen in the regime where individual pairs of filaments perform oscillations in the near-isometric state (Fig. 3). The synchronized oscillations can then be very pronounced after a small decrease in the load, as shown in Fig. 11 A. On the other hand, no damped oscillations are observed after a larger drop in the load, e.g. to T0/3 (Fig. 11 B), because the individual filament pairs immediately move out of the hysteretic regime. These properties are in agreement with recent experiments on the isotonic response of single muscle fibers carried out by Edman and Curtin (2001)
.
|
| DISCUSSION |
|---|
|
|
|---|
Our model shows that there are two different regimes of the microscopic dynamics in near-isometric conditions. For a range of values of the parameter
1 close to unity, the isometric point falls in the interval where the slope of the T2 curve is negative. In this case individual filament pairs oscillate with small amplitude. In the other regime, where the T2 curve has positive slope at the isometric point, the individual filaments are stationary, apart from stochastic fluctuations. The macroscopic manifestations of these two regimes differ in few respects. Because the oscillations of different filament pairs have different phases, oscillatory motion is not normally observable on the scale of a whole muscle fiber in steady conditions. However, a sudden change of load can synchronize the oscillations and thereby make them visible. The existence of damped oscillations in the isotonic transient response of single muscle fibers (Edman and Curtin, 2001
) therefore argues in favor of the oscillating regime. We note, however, that damped oscillations can also be a manifestation of step-wise shortening (Duke, 1999
), which can exist in both regimes. Because the efficient transduction of energy demands
1
1, which is close to value of this parameter at the boundary of the two regimes, it is possible that both regimes exist depending on conditions such as the myosin isoform, phosphate concentration, pH, ionic strength and temperature. Indeed, measurements by Edman and Curtin (2001)
show a dependence of the oscillation decay on the solution pH and on muscle fatigue. Further experiments, in which conditions are systematically varied, could shed more light on the mechanism of oscillation.
As a final remark, we emphasize that according to our model, the T2 curve of an individual pair of filaments differs from that of a muscle fiber. Recently, assays have been developed to measure the force-velocity relation of a single filament within a half sarcomere (Kawai et al., 2000
). In order to test our predictions, it would worthwhile to develop such high-precision techniques to measure the transient response of a single filament.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Submitted on November 5, 2002;
Revision received April 22, 2003.
| REFERENCES |
|---|
|
|
|---|
Barclay, C. J. 1998. Estimation of cross-bridge stiffness from maximum thermodynamic efficiency. J. Muscle Res. Cell Motil. 19:855864.[Medline]
Barsotti, R. J., and M. A. Ferenczi. 1988. Kinetics of ATP hydrolysis and tension production in skinned cardiac muscle of the guinea pig. J. Biol. Chem. 263:1675016756.
Brenner, B. 1991. Rapid dissociation and reassociation of actomyosin cross-bridges during force generation:a newly observed facet of cross-bridge action in muscle. Proc. Natl. Acad. Sci. USA. 88:1049010494.
Brenner, B., J. M. Chalovich, and L. C. Yu. 1995. Distinct molecular processes associated with isometric force generation and rapid tension recovery after quick release. Biophys. J. 68:106s111s.[Medline]
Chen, Y., and B. Brenner. 1993. On the regeneration of the actin-myosin power stroke in contracting muscle. Proc. Natl. Acad. Sci. USA. 90:51485152.
Cooke, R., and E. Pate. 1985. The effects of ADP and phosphate on the contraction of muscle fibers. Biophys. J. 48:789798.
Duke, T. 2000. Cooperativity of myosin molecules through strain-dependent chemistry. Philos. Trans. R. Soc. Lond. B Biol. Sci. 355:529538.[Medline]
Duke, T. A. J. 1999. Molecular model of muscle contraction. Proc. Natl. Acad. Sci. USA. 96:27702775.
Ebus, J. P., and G. J. Stienen. 1996. ATPase activity and force production in skinned rat cardiac muscle under isometric and dynamic conditions. J. Mol. Cell. Cardiol. 28:17471757.[Medline]
Edman, K. A. P., and N. A. Curtin. 2001. Synchronous oscillations of length and stiffness during loaded shortening of frog muscle fibres. J. Physiol. 534:553563.
Eisenberg, E., T. L. Hill, and Y. Chen. 1980. Cross-bridge model for muscle contraction. Biophys. J. 29:195227.
Ford, L. E., A. F. Huxley, and R. M. Simmons. 1977. Tension responses to sudden length change in stimulated frog muscle fibres near slack length. J. Physiol. 269:441515.
Ford, L. E., A. F. Huxley, and R. M. Simmons. 1985. Tension transients during steady shortening of frog muscle fibres. J. Physiol. 361:131150.
Ford, L. E., A. F. Huxley, and R. M. Simmons. 1986. Tension transients during the rise of tetanic tension in frog muscle fibres. J. Physiol. 372:595609.
Granzier, H. L., A. Mattiazzi, and G. H. Pollack. 1990. Sarcomere dynamics during isotonic velocity transients in single frog muscle fibers. Am. J. Physiol. 259:C266C278.[Medline]
Hill, A. V. 1939. The heat of shortening and dynamic constants of muscle. Proc. R. Soc. London Ser. B. 126:136195.
Hill, T. L. 1974. Theoretical formalism for the sliding filament model of contraction of striated muscle. Part I. Prog. Biophys. Mol. Biol. 28:267340.[Medline]
Howard, J. 2001. Mechanics of Motor Proteins and the Cytoskeleton. Sinauer, Sunderland, MA.
Huxley, A. F. 1957. Muscle structure and theories of contraction. Prog. Biophys. Biophys. Chem. 7:255318.[Medline]
Huxley, A. F. 2000. Cross-bridge action:present views, prospects, and unknowns. J. Biomech. 33:11891195.[Medline]
Huxley, A. F., and R. M. Simmons. 1971. Proposed mechanism of force generation in striated muscle. Nature. 233:533538.[Medline]
Huxley, A. F., and S. Tideswell. 1996. Filament compliance and tension transients in muscle. J. Muscle Res. Cell Motil. 17:507511.[Medline]
Huxley, A. F., and S. Tideswell. 1997. Rapid regeneration of power stroke in contracting muscle by attachment of second myosin head. J. Muscle Res. Cell Motil. 18:111114.[Medline]
Irving, M., G. Piazzesi, L. Lucii, Y. B. Sun, J. J. Harford, I. M. Dobbie, M. A. Ferenczi, M. Reconditi, and V. Lombardi. 2000. Conformation of the myosin motor during force generation in skeletal muscle. Nat. Struct. Biol. 7:482485.[Medline]
Jülicher, F., and J. Prost. 1995. Cooperative molecular motors. Phys. Rev. Lett. 75:26182621.[Medline]
Kawai, M., K. Kawaguchi, M. Saito, and S. Ishiwata. 2000. Temperature change does not affect force between single actin filaments and HMM from rabbit muscles. Biophys. J. 78:31123119.
Kojima, H., A. Ishijima, and T. Yanagida. 1994. Direct measurement of stiffness of single actin filaments with and without tropomyosin by in vitro nanomanipulation. Proc. Natl. Acad. Sci. USA. 91:1296212966.
Kushmerick, M. J., and R. E. Davies. 1969. The chemical energetics of muscle contraction. II. The chemistry, efficiency and power of maximally working sartorius muscles. Proc. R. Soc. Lond. B Biol. Sci. 174:315347.[Medline]
Piazzesi, G., L. Lucii, and V. Lombardi. 2002b. The size and the speed of the working stroke of muscle myosin and its dependence on the force. J. Physiol. 545:145151.
Piazzesi, G., M. Reconditi, M. Linari, L. Lucii, Y. B. Sun, T. Narayanan, P. Boesecke, V. Lombardi, and M. Irving. 2002a. Mechanism of force generation by myosin heads in skeletal muscle. Nature. 415:659662.[Medline]
Podolsky, R. J. 1960. Kinetics of muscular contraction:the approach to the steady state. Nature. 188:666668.[Medline]
Steffen, W., D. Smith, R. Simmons, and J. Sleep. 2001. Mapping the actin filament with myosin. Proc. Natl. Acad. Sci. USA. 98:1494914954.
Tyska, M. J., and D. M. Warshaw. 2002. The myosin power stroke. Cell Motil. Cytoskeleton. 51:115.[Medline]
Veigel, C., M. L. Bartoo, C. S. White, J. S. Sparrow, and J. E. Molloy. 1998. The stiffness of rabit skeletal acotmyosin cross-bridges determined with an optical tweezers transducer. Biophys. J. 75:14241438.
Veigel, C., F. Wang, M. L. Bartoo, J. R. Sellers, and J. E. Molloy. 2002. The gated gait of the processive molecular motor, myosin V. Nat. Cell Biol. 4:5965.[Medline]
Vilfan, A., E. Frey, and F. Schwabl. 1998. Elastically coupled molecular motors. Eur. Phys. J. B. 3:535546.
Vilfan, A., E. Frey, and F. Schwabl. 1999. Force-velocity relations of a two-state cross-bridge model for molecular motors. Europhys. Lett. 45:283289.
This article has been cited by other articles:
![]() |
A. Vilfan Myosin V Passing over Arp2/3 Junctions: Branching Ratio Calculated from the Elastic Lever Arm Model Biophys. J., May 1, 2008; 94(9): 3405 - 3412. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. S. Davis and N. D. Epstein Mechanism of Tension Generation in Muscle: An Analysis of the Forward and Reverse Rate Constants Biophys. J., April 15, 2007; 92(8): 2865 - 2874. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Lan and S. X. Sun Dynamics of Myosin-Driven Skeletal Muscle Contraction: I. Steady-State Force Generation Biophys. J., June 1, 2005; 88(6): 4107 - 4117. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Vilfan Elastic Lever-Arm Model for Myosin V Biophys. J., June 1, 2005; 88(6): 3792 - 3805. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS |