|
|
||||||||
| ABSTRACT |
|---|
|
|
|---|
H °; 70 ± 20 kJ mol-1), standard entropy change (
S °; 250 ± 70 J mol-1 K-1), and heat capacity change (
Cp; -12 ± 5 kJ mol-1 K-1). These results indicate that the force generation step is an entropy driven, endothermic reaction that accompanies a burial of large surface area. These observations are consistent with the hypothesis that hydrophobic interaction between residues of actin and myosin and between residues of the myosin head underlies the mechanism of force generation.
| INTRODUCTION |
|---|
|
|
|---|
It has been known for some time that tension is maintained at lower energy cost in slow-twitch muscle fibres (STFs) than in fast-twitch muscle fibres (FTFs) when contracting under isometric conditions (Crow & Kushmerick, 1982; Kushmerick et al. 1992; Barclay et al. 1993; Potma et al. 1995). However, STFs work less efficiently than FTFs when shortening with a light load (Curtin & Woledge, 1978; Rome et al. 1988). This functional difference is based on a different distribution of protein isoforms. While actin is common in both fibre types, most contractile proteins in FTFs are of the fast-twitch type, and most contractile proteins in STFs are of the slow-twitch type. It has been generally hypothesized that these molecular differences are the primary cause of a difference in function. One method of characterizing the function is to study elementary steps of the cross-bridge cycle, because a difference in protein isoforms is expected to contribute to a difference in the rate constants of elementary steps. A difference in elementary steps in turn would contribute to a difference in performance of the muscle fibres. Consequently, we decided to correlate fibres of different protein isoforms with a difference in elementary steps of the cross-bridge cycle. In our previous studies (Wang & Kawai, 1996, 1997), we characterized the elementary steps of STFs at 20°C and compared the results with those from FTFs. Ideally, these studies need to be extended to physiological temperature, so that the knowledge obtained becomes more relevant to the contractile mechanisms of mammalian muscles and their diseases.
The temperature study yields yet additional information. It is possible to obtain insights into the molecular mechanisms of contraction from the temperature sensitivity of elementary steps of the cross-bridge cycle. For instance, from a temperature study in the range 5-30°C in rabit psoas (fast-twitch) fibres, Zhao & Kawai (1994) and Murphy et al. (1996) deduced thermodynamic parameters including the enthalpy, entropy and heat capacity changes, and concluded that the molecular force involved in force generation is based on the hydrophobic interaction between amino acid residues of actin and myosin. Similarly, information has been derived from temperature studies on psoas fibres in the low temperature range (Goldman et al. 1987; Ranatunga, 1984, 1989, 1996). However, attempts to study the cross-bridge kinetics at temperatures higher than 30°C have not been satisfactory, for FTFs become less stable at higher temperatures and the cross-bridge kinetics become faster, and hence more difficult to observe. STFs have the advantage over FTFs because they are more stable, particularly at high temperatures, and because they have elementary steps that are 10-30 times slower than those of FTFs (Wang & Kawai, 1996, 1997). With STFs that are more resistant to physical damage and slower enzymatic reactions, we can study the effect of temperature on the elementary steps with better time resolution.
The present study was aimed at finding a difference in thermodynamic parameters between STFs and FTFs, and to correlate the difference to a known structural (molecular) difference. We studied the effects of MgATP and phosphate (Pi) on the exponential processes of tension time course in response to sinusoidal length changes at temperatures ranging from 5 to 37°C. This is a logical extension of our previous work on cross-bridge kinetics of STFs (Wang & Kawai, 1996, 1997) to the physiological temperature range. Our results indicate that (1) the apparent rate constants and the rate constants of elementary steps of the cross-bridge cycle are more temperature sensitive in soleus STFs than in psoas; and (2) the force generation step is the most temperature sensitive step and it accompanies large enthalpy, entropy and heat capacity changes.
| METHODS |
|---|
|
|
|---|
Chemicals and solutions
Chemicals and solutions were prepared as described in Wang & Kawai (1996, 1997). In brief, the relaxing solution contained (mM): 5.03 EGTA, 0.97 CaEGTA (pCa 7.0), 2 MgATP, 5 free ATP, 55 potassium propionate (Kprop), 61 sodium propionate (NaProp), 10 Mops and 10 NaN3. The control activating solution contained (mM): 6 CaEGTA, 0.174 CaProp2 (pCa 4.40), 5.76 MgATP, 1.36 free ATP, 15 creatine phosphate (CP), 8 phosphate (Pi), 53 KProp, 1 NaProp, 10 Mops and 10 NaN3. The rigor solution contained (mM): 8 Pi, 76 NaProp, 103 KProp and 10 Mops. The WS solution was used to remove EGTA and MgATP before calcium activation in the MgATP study, and it contained (mM): 0.84 free ATP, 8 phosphate (Pi), 85.6 KProp, 15 NaProp, 15 CP, 10 Mops and 0.16 mg ml-1 creatine kinase (CK; Boehringer Mannheim GmbH, Germany; 1 mg = 500 U). WP solution was used to remove EGTA and Pi before calcium activation in a Pi study, and it contained (mM): 5.82 MgATP, 1.17 free ATP, 90 KProp, 1 NaProp, 15 CP, 10 Mops and 0.16 mg ml-1 CK. The CK level of all activating solutions was 0.16 mg ml-1 for a temperature of 20°C or less, and it was elevated to 0.32 mg ml-1 for studies carried out at 25°C, to 0.48 mg ml-1 for 30°C, to 0.64 mg ml-1 for 34°C and to 0.8 mg ml-1 for 37°C. The amount of CK we used is sufficient, because the activity of CK also increases with temperature with a Q10 of about 2.2 (Kuby et al. 1954). In all solutions used for experiments, the ionic strength was adjusted to 180 mM with NaProp/KProp, because this is close to physiological ionic strength (Godt & Maughan, 1990), except for the relaxing and rigor solutions that had an ionic strength of 200 mM. The pH of all the solutions was adjusted to 7.00 ± 0.01; the pH of the activating solutions did not change with temperature. EGTA, CaEGTA and Pi were added as neutral potassium salts (their approximate chemical formulae are K2H2EGTA, K2CaEGTA and K1.5H1.5PO4, respectively); MgATP and CP were added as neutral sodium salts (Na2MgATP and Na2CP); free ATP was added as a neutral sodium-potassium salt (approximately Na2K1.7H0.3ATP). Individual concentrations of multivalent ionic species were calculated using our computer program, which assumed multiple equilibria with the following apparent association constants (log values at pH 7.00): CaEGTA 6.28, MgEGTA 1.61, CaATP 3.70, MgATP 4.00, CaCP 1.15 and MgCP 1.30. These values were not corrected for the temperature change, because of a paucity of reliable data obtained under the same condition, and because the temperature may not have a large effect on these association constants. The calculated pCa of all activating solutions was 4.40, and [Mg2+] was 0.5 mM.
Fibre preparations and experimental procedure
Rabbit soleus single slow-twitch fibres (length: 2.4-4.7 mm; diameter: 72-105
m) were obtained as described in our earlier works (Wang & Kawai, 1996, 1997). The use of rabbits conformed with the current 'Guide for the Care and Use of Laboratory Animals' (NIH publication DHSS/USPHS), and was approved by the Institutional Animal Care and Use Committee. Rabbits were killed by intravenous injection of sodium pentobarbital at 150 mg kg-1 in the ear vein. Both ends of fibre preparations were doubly knotted and connected to fish hook-shaped clamps of the length driver and force transducer at 20°C in the relaxing solution. The force transducer was based on an AE801 gauge element (Senso Nor a.s, Horten, Norway). According to the manufacturer's specification, the element was 4.5 % more sensitive at 5°C, and 2.3 % less sensitive at 37°C compared to the sensitivity at 20°C. However, the temperature variation in the element was approximately between 15 and 30°C because of the location of the element, and hence the maximum sensitivity variation was in the range +2.3 to -1.2 % in the experiments reported here. Because this is a small variation, we did not attempt to correct for it. The muscle length (L0) was determined as described in Wang et al. (1999) at the sarcomere length 2.5
m, at which all experiments were performed. This method eliminates the effect of end compliance when present. The baseline record of the complex modulus was collected in the relaxing solution. This was replaced twice by wash solution (WS or WP) to remove EGTA and to add CP and CK before activation. The wash solution was replaced with an activating solution (pCa 4.40). When a steady tension developed, an experimental record of the complex modulus was collected. The fibre was then relaxed with two solution changes. Temperature was changed to the desired experimental temperature at this stage. In the experiments that studied the MgATP or Pi effect, the above sequence was repeated for solutions that contained different concentrations of MgATP or Pi. At the end of each series, temperature was brought back to 20°C, and isometric tension with the control activating solution was tested for reproducibility. Rigor was induced by replacing the control activating solution with the rigor solution (Rig) three times. When the full rigor state developed in about 1 min, the complex modulus data of the rigor state were collected. A representative experimental protocol is shown in Fig. 1.
![]() |
View larger version [in this window] [in a new window] |
|
|
A horizontal bar below the force time course indicates where a computer record was taken and is shown with the record number. The computer record is reflected by length oscillation (0.25 % L0 peak to peak in all frequencies). Both L and F signals were filtered by 3 Hz low-pass filters before being traced by pens, and hence the amplitude at the higher frequency range appears diminished. A single fibre (rabbit soleus STF) was brought to 20°C, a computer record (no. 2) was taken during relaxation, and the fibre was washed twice with the wash solution and activated by the control activating solution (
| ||
The complex modulus Y(f) is defined as the ratio of the stress change to the strain change, and it is a function of frequency (f). Y(f) is fitted to an equation consisting of four terms, each one representing an exponential process (Wang & Kawai, 1996, 1997):
![]() | (1) |
![]() | (2) |
where i =
-1. Each term is called an exponential process, because its time domain expression takes the form of an exponential function. The real part of Y(f) is called the elastic modulus, the imaginary part is the viscous modulus, the absolute value of Y(f) is the dynamic modulus, and the angle defined by Y(f) on the complex plane is the phase shift. At lower temperatures, the complex modulus data were collected in frequencies ranging from 0.04 to 250 Hz, which corresponds to 0.6 ms to 4 s in time domain analysis. At higher temperatures, low frequency points were omitted to reduce the duration of the Ca2+ activation. The amplitude of the length oscillation was fixed to 0.125 % L0 at all frequencies. This length oscillation corresponds to 1.6 nm per half sarcomere. We denote the characteristic frequencies of the respective processes by a, b, c and d, and their magnitudes by A, B, C and D. H is the modulus extrapolated to the zero frequency, and it remained as a small positive value compared to Y
as reported (Kawai & Brandt, 1980) in all experimental conditions. The characteristic frequencies multiplied by 2
are the apparent (= observed) rate constants. Y
defined in eqn (2) is the complex modulus extrapolated to the infinite frequency, often referred to as 'stiffness'. Mathematical transformation demonstrates that Y
corresponds to phase 1 of step analysis, process D corresponds to the fast component of phase 2, process C to the slow component of phase 2, process B to phase 3, and process A to phase 4 (Kawai & Brandt, 1980). The phases of step analysis were previously defined by Huxley & Simmons (1971). Details of the sinusoidal analysis method were published in Kawai & Brandt (1980).
Cross-bridge scheme
The rate constant results were analysed based on a cross-bridge scheme (Scheme 1), which was used for rabbit psoas FTFs (Kawai & Halvorson, 1991; Zhao & Kawai, 1993), soleus STFs (Wang & Kawai, 1996, 1997), and ferret myocardium (Kawai et al. 1993). This scheme is consistent with the effects of MgATP, MgADP and Pi on exponential processes B, C and D. This scheme has a close resemblance to that derived from solution studies of extracted contractile proteins (Taylor, 1979), except that the Pi release step is reversible in muscle fibres.

In Scheme 1, S indicates MgATP, D indicates MgADP and P indicates phosphate (Pi). In step 1a, substrate MgATP binds to the rigor-like cross-bridge state AM (actomyosin) to form a collision complex AM+S, and in step 1b, AM+S isomerizes to form the AM*S state. Cross-bridges then detach at step 2 to form the Det state, which includes weakly attached states (AMS, AMDP) and truly detached states (MS, MDP). The transition between AM+S and AM*S states (step 1b) and the transition between AM*S and Det states (step 2) can be probed by studying the effect of MgATP on the apparent rate constants 2
d and 2
c, respectively. These studies also yield the MgATP binding constant K1a to the AM state. ATP cleavage occurs in step 3. Once the thin filament is activated by Ca2+ binding to troponin C (TnC), weakly attached AMDP state isomerizes to form the strongly attached AM*DP state in step 4. It is at this step that force is generated (Fortune et al. 1991; Kawai & Halvorson, 1991; Dantzig et al. 1992; Wang & Kawai, 1997). The force generation is followed by Pi release (step 5) to form the AM*D state. The force generation and the subsequent Pi release steps can be probed by studying the effect of Pi on the apparent rate constant 2
b. b, c and d are deduced by fitting eqn (1) to the complex modulus data. AM*D slowly isomerizes to the AMD state, which limits the ATP hydrolysis rate.
| RESULTS |
|---|
|
|
|---|
Because about 89 % of the rabbit soleus fibres are slow-twitch type I fibres (STFs) and the remaining 11 % are fast-twitch type IIa fibres (FTFs) (Wang & Kawai, 1996), we studied only fibres which showed a slow-twitch response. Since type I fibres are 10 times slower than type IIa fibres, and the frequency data and the isoforms are perfectly correlated (Kawai & Schachat, 1984), identification of type I fibres can be readily made (Wang & Kawai, 1996).
Temperature effect on tension and stiffness
Figure 1 is a slow pen trace record demonstrating the effect of temperature from 5 to 37°C when the rabbit soleus STF was activated in the control activating solution (pCa 4.40). After nine activations (each lasting for 0.4-3 min) the tension recovery was superb in this soleus fibre and reproduced 98 % of the initial tension (Fig. 1, Rec no. 40 vs. Rec no. 4). Six additional experiments were performed, in which temperature was increased in half of the experiments and decreased in the other half. Tension of the control activation (Tc) at 20°C was 171 ± 7 kN m-2 (mean ± S.E.M., n = 7), and the recovery of the tension of the last control activation was 102 ± 3 % compared to the first control activation. The maximum calcium activated tension, stiffness (Y
), and the tension to stiffness ratio in the control activating solution as functions of the temperature are shown in Fig. 2. Isometric tension increased with temperature in the range 5-30°C, and saturated in the range 30-37°C. Stiffness increased with temperature in the range 5-20°C and saturated in the range 20-37°C. The tension to stiffness ratio increased from 0.90 % L0 to 1.35 % L0 with temperature in the range 5-30°C, and saturated when temperature was further increased from 30 to 37°C.
![]() |
View larger version [in this window] [in a new window] |
|
|
A, isometric tension (
| ||
Temperature effect on the complex modulus and the apparent rate constants
Figure 3 demonstrates the complex modulus data shown in the frequency plots (in A, B, D and E) and in the Nyquist plots (in C and F) during the activation with the control solution and at varying temperatures. The data from 5 to 15°C are shown in Fig. 3A-C, and the data from 20 to 37°C are shown in Fig. 3D-F. A close examination of Fig. 3F demonstrates that the Nyquist plot in the temperature range 20-37°C is typical of skeletal muscle fibres (Kawai & Schachat, 1984; Zhao & Kawai, 1993; Wang & Kawai, 1996) as represented by a work loop, but the Nyquist plots deviate from the typical shape at lower temperatures as shown in Fig. 3C. Figure 3A, B, D and E demonstrate that both the dynamic modulus vs. frequency plot and the phase shift vs. frequency plot shift to the right as the temperature is elevated, indicating that the cross-bridge kinetics became faster when temperature was increased. Equation (1) was then fitted to the data in Fig. 3 (15-37°C) to deduce the parameters of the exponential processes. Note that eqn (1) could not be fitted to some data at 15°C, and to much fewer data at 5-10°C, and hence these data could not be used for further analysis. This difficulty arose because at 5-10°C the centre of the arc that represents processes C and D is displaced below the abscissa (the arc is flattened) in the Nyquist plot (Fig. 3C) indicating the presence of distributed rate constants, and because process B is diminished or absent (Fig. 3C) at 5-10°C. These problems caused unreliable estimation of the rate constants. The theoretical projection for the data corresponding to ˛geł 20°C is given in Fig. 3G, H and I.
![]() |
View larger version [in this window] [in a new window] |
|
|
The data were obtained during control activation as shown in Fig. 1. The data obtained from 5 to 15°C are shown in A-C, and the data from 20 to 37 °C are shown in D-F. Equation (1) was fitted to the data in D-F, and a theoretical projection is given in G-I. The dynamic modulus vs. frequency plot is shown in A, D and G; the phase shift vs. frequency plot is shown in B, E and H; and the viscous modulus vs. elastic modulus plot (Nyquist plot) is shown in C, F and I. The phase shift is positive on phase advance (tension change leads length change), and negative on phase delay (tension change lags behind length change). Each curve represents the average of 7 experiments. Data were obtained from the experiments shown in Fig. 2.
| ||
The effect of temperature on four apparent rate constants is plotted in Fig. 4 for ˛geł 15°C. This figure demonstrates that all the apparent rate constants, 2
a, 2
b, 2
c and 2
d, increased with temperature. The Q10 values of these rate constants in the range 20-30°C are listed in Table 1 and compared with the data obtained from rabbit psoas fibres. This table demonstrates that the Q10 for 2
b is the highest with a value of 7.4, followed by 2
a (5.7), 2
c (3.7) and 2
d (2.1), and they respectively compare to 3.4, 2.5, 1.9 and 0.9 in psoas fibres (Zhao & Kawai, 1994) in the same temperature range. Therefore, it can be concluded that the apparent rate constants in soleus STFs are more temperature sensitive than in psoas (FTFs). The apparent rate constant 2
a was not used for further analysis, because it may not be limited by an elementary step (Wang & Kawai, 1997).
![]() |
View larger version [in this window] [in a new window] |
|
|
The data were obtained from individual experiments similar to that shown in Fig. 3. The rate constants were then averaged for 7 experiments, except for 15°C where the data were from 4 experiments, and shown with S.E.M. error bars.
| ||

The MgATP effect on exponential processes C and D at five temperatures
To examine the temperature sensitivity of the elementary steps surrounding the substrate binding, the effect of MgATP on exponential processes C and D at different temperatures was studied in the range 0.1-5 mM. In Fig. 5A, the apparent rate constant 2
d is plotted against MgATP concentration at temperatures 20, 25, 30, 34 and 37°C. These plots reveal that the MgATP dependence of 2
d is hyperbolic: 2
d increases at low millimolar concentrations, and it approaches saturation by 2-5 mM. Such MgATP dependence can be explained by eqn (3):
![]() | (3) |
where S is the MgATP2- (substrate) concentration, K1a is the association constant of MgATP to cross-bridges, k1b is the rate constant of the ATP-isomerization step, k-1b is the rate constant of its backward reaction, and K1b (= k1b/k-1b) is its equilibrium constant.
![]() |
View larger version [in this window] [in a new window] |
|
|
A, the apparent rate constant 2
| ||
In Fig. 5B, the apparent rate constant 2
c is similarly plotted. These plots reveal that the MgATP dependence of 2
c is also hyperbolic: 2
c increases at low millimolar concentrations, and approaches saturation by 2-5 mM. Such MgATP dependence can be explained by eqn (4):
![]() | (4) |
where k2 is the rate constant of the cross-bridge detachment step 2, k-2 is the rate constant of its reversal step, and K2 = k2/k-2 is the equilibrium constant. Equations (3) and (4) were derived based on a four-state model consisting of AM, AM+S, AM*S, and Det states (Scheme 1), and formulated in our earlier publications (Zhao & Kawai, 1993). Equations (3) and (4) were fitted simultaneously to the data in Fig. 5 to deduce the MgATP association constant, the rate constants and the equilibrium constants. These constants are collectively called 'kinetic constants', and are indicated in Scheme 1. The theoretical projections based on eqns (3) and (4) with best fit parameters are shown by the continuous lines in Fig. 5A and B. The above experiments were carried out in the presence of CP and CK, and hence the ADP concentration is maintained at less than 20
M (Meyer et al. 1985), so that the effect of MgADP can be ignored in eqns (3) and (4).
The kinetic constants were deduced from these fittings and plotted in Fig. 6 as functions of 1000/T, where T is the absolute temperature expressed in kelvins. The Celsius temperature (°C) is also shown at the top of the figure. As seen in Fig. 6, the MgATP binding constant K1a decreased when the temperature was ˛geł 25°C, indicating that MgATP binding becomes weaker at higher temperatures, i.e. this binding reaction is exothermic (
H ° < 0). The average K1a at 20°C was slightly smaller than that at 25°C, but their difference is insignificant. It is also interesting to point out that the MgATP association constant is larger in STFs than in FTFs at all temperatures studied (compare Fig. 6A with Fig. 5A of Zhao & Kawai, 1994). The rate constant of the ATP-isomerization step (k1b) and that of its reversal step (k-1b) increased with a similar Q10 (Table 2), and thus the equilibrium constant of this reaction did not change much with temperature (Fig. 6B). Likewise, the rate constant of the cross-bridge detachment step (k2) and its reversal step (k-2) increased with temperature with a similar Q10 (Table 2), and thus the equilibrium constant K2 did not change much with temperature (Fig. 6B).
![]() |
View larger version [in this window] [in a new window] |
|
|
A, the MgATP association constant (K1a) is plotted against temperature. B, the rate and equilibrium constants of the ATP-isomerization step 1b and those of the cross-bridge detachment step 2 are plotted against temperature. The data were based on the MgATP study on 2
| ||

The phosphate (Pi) effect on exponential process B at five temperatures
To characterize the temperature dependence of the force generation step 4 and Pi release step 5, the effect of Pi on process B was studied at five temperatures. Isometric tension declined as the Pi concentration was increased from 0 to 30 mM as reported in Wang & Kawai (1997) for 20°C, and this trend did not change as the temperature was varied. Stiffness followed the same pattern as isometric tension. In Fig. 7 the apparent rate constant 2
b is plotted against the Pi concentration at temperatures 20, 25, 30, 34 and 37°C. These plots reveal that the Pi dependence is hyperbolic: 2
b increases in the range of 0-8 mM Pi, and it approaches saturation at higher Pi concentrations. Such Pi dependence can be explained by eqn (5):
![]() | (5) |
where
![]() | (6) |
and P represents the Pi concentration. Equation (5) was derived based on a three-state model consisting of AMDP, AM*DP, and AM*D and a rapid equilibrium to the left of the AMDP state as previously described (Kawai & Halvorson, 1991; Dantzig et al. 1992). In Fig. 7 the continuous lines represent theoretical projections based on eqn (5).
was calculated by eqn (6) with K1a, K1b and K2 obtained from the MgATP study and S = 5 mM.
remained in the range 0.44 and 0.57, and it did not change much with temperature. The kinetic constants of the elementary steps are plotted in Fig. 8 as functions of 1000/T. As seen in this figure the rate constant of the force generation step (k4) increased steeply with temperature. It can also be seen that the rate constant of its reversal step (k-4) increased with temperature, but to a lesser extent than k4. Consequently the equilibrium constant K4 increased significantly with temperature (Fig. 8). The Pi association constant (K5) did not change much with temperature (Table 2). Q10 for k4 was 6.7, and Q10 for k-4 was 2.5 in the range 20-30°C. The corresponding Q10 values for rabbit psoas were 3.3 and 1.5, respectively (Table 2). From these results we conclude that the force generation step and its reversal step are more temperature sensitive in soleus STFs than in psoas FTFs.
![]() |
View larger version [in this window] [in a new window] |
|
|
b is plotted against Pi concentration at five temperatures
For these experiments, a fibre was mounted on the apparatus, the control activation was performed at 20°C, and the temperature was changed to the desired experimental temperature. This was followed by a series of activation and relaxation cycles with different Pi concentrations (0, 2, 4, 8, 16, 23 and 30 mM) in increasing or decreasing order. At the end, the temperature was brought back to 20°C, the control activation was repeated, and this was followed by rigor induction. The fibre was then discarded. This procedure was repeated for several STFs. The data are means with S.E.M. error bars. The number of fibres used was 10 (20°C), 7 (25°C), 7 (30°C), 7 (34°C) and 8 (37°C), except there were only 4 experiments for 30 mM Pi at 37°C. The continuous curves are based on eqn (5) and the best fit parameters.
| ||
![]() |
View larger version [in this window] [in a new window] |
|
|
The data were based on the Pi study on 2
| ||
Further analysis of the temperature effect
Once the rate and equilibrium constants are obtained, one can calculate the activation energy, the standard enthalpy change, and the standard entropy change of each reaction step. The activation energy (Ea) and the frequency factor (F) were calculated based on the Arrhenius equation (eqn (7)), and listed in Table 2:
![]() | (7) |
where k is the rate constant of an elementary step, R is the gas constant, and T is absolute temperature. The activation energy was largest (140 kJ mol-1) for the force generation step; it was 80-110 kJ mol-1 for ATP isomerization and cross-bridge detachment steps; and it was smallest (69 kJ mol-1) for the reversal of the force generation step. The activation energy of each elementary step of soleus slow-twitch fibres is 1.1-2.5 times larger than that of psoas FTFs (Table 2).
The standard enthalpy change (
H °), the standard entropy change (
S °), and the standard free energy change (
G°) were calculated by the Van't Hoff equation:
![]() | (8) |
![]() | (9) |
where K is an equilibrium constant. Equation (8) was fitted to the equilibrium constant data in Fig. 6 and 8, and
H ° and
S ° were deduced and listed in Table 3. This table shows that
H ° and
S ° are by far the largest and are positive in step 4 (force generation step), followed by step 2b (ATP-isomerization step), but
H ° and
S ° are close to 0 for the Pi binding and cross-bridge detachment steps.
H ° is a negative number for the ATP binding step.

Because the large positive
H ° and
S ° values imply a formation of hydrophobic interaction (Frank & Evans, 1945; Privalov & Gill, 1988), which accompanies burial of a large surface area, we used the modified Van't Hoff equation (Murphy et al. 1996) and fitted it to K4:
![]() | (10) |
where
Cp is the heat capacity change, and Tr is an arbitrary reference temperature which we chose as 25°C. The derivation of eqn (10) is based on the first-order Taylor expansion of
H ° and
S ° around Tr:
![]() | (11) |
![]() | (12) |
i.e.
H ° is a linear function of the temperature,
S ° is a linear function of lnT, and
Cp is their common coefficient.
Hr° and
Sr° are respective values of
H ° and
S ° at temperature Tr. Figure 9 is a plot of K4 against 1000/T, and the continuous line represents the theoretical projection based on eqn (10). From the fitting, we found
H ° to be 72 ± 20 kJ mol-1,
S ° to be 247 ± 66 kJ mol-1 K-1 and
Cp to be -12.4 ± 5.0 kJ mol-1 K-1 (Table 4).
Cp does not depend on the choice of Tr, but
H ° and
S ° change with temperature (eqns (11) and (12)). The data are best fit values ± 95 % confidence limits based on Student's t test. Thus, the force generation step 4 is an endothermic reaction (
H ° > 0), and it has a large entropy increase. These numbers are comparable to those obtained from rabbit psoas calculated for the same Tr (Table 4).
Cp for the soleus is negative and twice as large as that for the psoas, but the confidence limit is much larger for the soleus than the psoas, and hence
Cp of soleus and psoas may not be different. The larger confidence limit was caused by a larger scatter of the data in soleus fibres. A large negative
Cp would imply burial of a large surface area.

Prediction of isometric tension as a function of temperature
To examine whether the change in the kinetic constants can explain the temperature effect on isometric tension, we deduced the steady-state probability of each cross-bridge state at different temperatures based on the equilibrium constants obtained (Fig. 6 and 8) and eqns (20)-(27) of Kawai & Zhao (1993). The calculation was based on the standard control activating condition (5 mM MgATP, 8 mM Pi). According to this calculation the fraction of all attached cross-bridges (X1a + X1b + X2 + X5 + X6) is 76 % (20°C), 83 % (25°C), 86 % (30°C), 85 % (34°C) and 81 % (37°C). Theoretical tension is then calculated based on eqn (13) assuming that isometric tension can be modelled as a linear combination of the probabilities of cross-bridge states (Kawai & Zhao, 1993):
| Tension = T1aX1a + T1bX1b + T2X2 + T5X5 + T6X6, | (13) |
where the linear coefficient Ti is the value of tension if 100 % of the cross-bridges are in state Xi. The Ti values of soleus STFs were empirically determined (Wang & Kawai, 1997) and are: T1a = T1b = 1.12Tc, T2 = 0.81Tc and T5 = T6 = 1.29Tc, where Tc is isometric tension in the control activating condition. Tc is taken as a constant measured at 20°C, but the equations for Ti do not imply that Ti is a function of temperature. The theoretical tension values based on eqn (13) are entered in Fig. 2A as a dotted line along with isometric tension, assuming that Ti does not change with temperature. As seen in this figure, the theoretical tension matches well with the experimental data, which is consistent with the idea that the number of tension generating cross-bridges increases with temperature, but tension per cross-bridge state remains unchanged.
Adequacy of the ATP regenerating system
To test whether the ATP regenerating system used in the current report is adequate, we studied the effect of CK and CP on 2
b and 2
c at three temperatures, 20, 30 and 37°C. Experiments were carried out at 1 mM MgATP2-, at which 2
c (or 2
b) is most sensitive to a change in the MgATP2- concentration (Fig. 5). 2
c (or 2
b) would decrease if the MgATP2- concentration was decreased, or if the MgADP- concentration was increased (Zhao & Kawai, 1993; Wang & Kawai, 1996). The effect of CK was tested by doubling its level, and the effect of CP was tested by increasing its concentration from 15 mM (standard) to 25 mM. Then the ratio of 2
c (or 2
b) after:before the increase was calculated and summarized in Table 5. If the ATP regenerating system were not adequate, the ratio would become > 1. The fact that the ratio remained < 1 (Table 5) demonstrates that the regenerating system was adequate in our experimental conditions.

| DISCUSSION |
|---|
|
|
|---|
Difference between STFs and FTFs at the molecular level
One purpose of the present study was to correlate the functional difference between STFs and FTFs with the known structural (molecular) difference. Therefore, it is necessary to discuss the difference between STFs and FTFs at the molecular level. Although the two muscle types use the same class of contractile proteins, their difference lies in protein isoforms. In fact, the mechanism of the high adaptability of a specialized muscle function is the result of a cell's ability to generate isoforms. Myosin and the thin filament proteins tropomyosin (Tm), TnT, TnI, and TnC are present in different isoforms in different fibre types, although actin is common. FTFs consist of a fast myosin heavy chain (MHC), three fast myosin light chains (MLCs),
- and
-Tm, and fast TnI, TnT and TnC. STFs consist of a slow MHC, two slow MLCs,
-,
- and
-Tm, and slow TnI, TnT and TnC. Myosin ATPase activity is largely determined by MHC isoforms (Pette & Staron, 1990; Moss et al. 1995; Schiaffino & Reggiani, 1996; Weiss et al. 1999), and hence it is likely that a difference in MHC is the major factor in contributing to the difference in the rate constant of an elementary step. A close examination of the amino acid sequence of MHCs of mammalian striated muscles reveals that the primary sequence is highly conserved among MHC isoforms. However, several regions are significantly different. These include limited regions of the myosin head, the
-helix in the neck region, the hinge between S2 and LMM, and the C-terminal of the MHC (McNally et al. 1989; Moore et al. 1992; Schiaffino & Reggiani, 1996; Weiss et al. 1999). Of particular interest are the two surface loops near the ATP-binding site (loop 1) and actin-binding site (loop 2) in the myosin head. These loops vary in length and amino acid sequence. This variation may determine the ATPase activity and the sliding velocity of the thin filament in vitro, and the maximum velocity of shortening in vivo (Spudich, 1994; Schiaffino & Reggiani, 1996; Weiss et al. 1999). One site with large variability is the junction between the 50 and 20 kDa domains (loop 2) in the myosin head. This site is located at the actin-myosin interface (Rayment et al. 1993a,b), and may contribute to a different ATPase activity (Uyeda et al. 1994).
Tension and stiffness
Isometric tension in the control activating solution of soleus STFs increased when temperature was raised from 5 to 30°C, and plateaued over a further rise to 37°C (Fig. 2A). This observation is consistent with that of Ranatunga (1996) on rabbit psoas and soleus fibres. Stiffness increased similarly to tension but to a lesser extent, and the tension-to-stiffness ratio increased in the low temperature range and saturated in the high temperature range (Fig. 2). These results are generally consistent with previous reports on mammalian muscle fibres (Ranatunga & Wylie, 1983; Stephenson & Williams, 1985; Goldman et al. 1987; Zhao & Kawai, 1994; Ranatunga, 1996; Galler & Hilber, 1998), frog tibialis anterior fibres (Ford et al. 1977) and toad muscles (Stephenson & Williams, 1985). The plateau at the high temperature range in the tension or stiffness plot is not caused by rigor cross-bridges, because the parameter H (eqn (1)) did not increase in the high temperature range. H would increase if non-cycling cross-bridges were formed.
Complex modulus and exponential processes
Figure 3C demonstrates that at lower temperatures such as 5 and 10°C, the Nyquist plot of soleus fibres activated with saturating Ca2+ does not exhibit the work loop as seen at higher temperatures (Fig. 3F; Kawai & Schachat, 1984; Wang & Kawai, 1996). This observation may mean that muscle contracts somewhat differently at lower temperatures.
Figure 4 shows that all the apparent rate constants significantly increased with an increase in temperature, which demonstrates that the rate constants are from typical enzyme reactions. The results are similar to those on rabbit psoas fibres (Zhao & Kawai, 1994), except that a higher temperature sensitivity was observed in soleus STFs than in psoas (Table 1). This result is consistent with a higher temperature sensitivity of the maximum shortening velocity in STFs than in FTFs (Ranatunga, 1984). It is possible that the slower cross-bridge kinetics and higher temperature sensitivity of the apparent rate constants may be correlated to a difference in loop 2 of the MHC of type I and type IIb fibres (Murphy & Spudich, 1999; Weiss et al. 1999).
One advantage in studying soleus STFs is that we can obtain the cross-bridge kinetics at physiological temperature. We found that the cross-bridge Scheme 1 derived at 20°C is well suited to describe the cross-bridge kinetics at the physiological temperature. A comparison of the kinetic constants of soleus STFs at 37°C (Fig. 4, 6 and 8) with those from psoas (Zhao & Kawai, 1994) demonstrates that the cross-bridge kinetics of STFs at the physiological temperature are close to the cross-bridge kinetics of psoas at 15-20°C. One such example is a presence of the work loop in the Nyquist plot in soleus STFs (Fig. 3F) as in psoas fibres.
Nucleotide binding step
It is interesting to point out that the MgATP2- association constant is 2-3 times larger in STFs than in FTFs in the temperature range 20-30°C (compare Fig. 6A with Fig. 5A of Zhao & Kawai, 1994). We infer from this observation that the structure of the nucleotide binding site or its microenvironment is different in the two fibre types. The nucleotide binding site forms a pocket in the N-terminal 25 kDa domain of the myosin head (Rayment et al. 1993b). The pocket is nearly perfectly conserved between type I (rabbit soleus STFs) and type IIb (rabbit psoas) myosin, and hence it may not be the pocket structure that determines the difference in the binding affinity. In terms of the microenvironment, it is known that loop 1 comes close to the pocket (Rayment et al. 1993b), and in fact, there is a significant difference in loop 1 between the two myosin isoforms (Weiss et al. 1999). In particular, the net charge of loop 1 is +2 and more positive in STF myosin (3 K+, 1 R+, 2 E-) than in FTF myosin (4 K+, 3 D-) with the net charge of +1, and thereby STF myosin should attract negatively charged substrate (MgATP2-) twice as much as FTF myosin. This mechanism would result in a higher MgATP association constant in STFs than in FTFs as observed. In fact, it has been suggested that the different length and/or the flexibility of loop 1 may cause the different nucleotide binding affinity and the ATP cleavage rate (Kelly et al. 1993; Rovner et al. 1997; Sweeney et al. 1998; Weiss et al. 1999). The higher nucleotide affinity may also be correlated with the fact that loop 1 is less flexible in ST myosin than in FT myosin (Weiss et al. 1999).
Figure 6A demonstrates that the MgATP2- binding affinity (K1a) for cross-bridges becomes weaker with an increase in temperature (exothermic), as has been the case in FTFs (Zhao & Kawai, 1994). The decrease in K1a may be correlated with the increased isometric tension. Such correlation was demonstrated in a variety of conditions that change isometric tension: 2,3-butanedione monoxime (BDM), [5-(1-(3,4-dimethoxybenzoyl)-1,2,3,4-tetrahydro-6-chinolyl)-6-methyl-3,6-dihydro-2H-1,3,4-thiadiazin-2-one]-53998 (EMD), TnC extraction and temperature change as summarized in Zhao et al. (1996). If the position of loop 1 is coupled with the amount of tension in the myosin head, this mechanism would cause a change in the MgATP2- binding constant.
The rate constants of both cross-bridge detachment step 2 (k2) and its reversal step (k-2) increase similarly with temperature (Fig. 6B). Therefore, the equilibrium constant of the cross-bridge detachment step (K2) does not change much with temperature. This finding is similar to that obtained in psoas (Zhao & Kawai, 1994). As seen in Table 2, the rate constants k2 and k-2 are more temperature sensitive in soleus than in psoas. This finding is in accord with the observation that the apparent rate constant 2
c is more temperature sensitive in soleus STFs than in psoas fibres (Table 1).
Force generation step
Step 4 is the most temperature sensitive step among the elementary steps in the cross-bridge cycle (Table 2 and Fig. 8), represented by a high Q10 (6.7) and a large activation energy (140 kJ mol-1) that corresponds to 2-2.5 times the free energy liberated by ATP hydrolysis. This finding is consistent with our earlier studies on rabbit psoas (Zhao & Kawai, 1994). Step 4 corresponds to the force generation step in muscle fibres (Kawai & Halvorson, 1991; Fortune et al. 1991; Dantzig et al. 1992; Kawai & Zhao, 1993; Wang & Kawai, 1997). The equilibrium constant (K4) of this step increases with temperature in the range 20-30°C, saturates in the range 30-34°C, and decreases for further increases in temperature to 37°C (Fig. 8). The presence of the peak in K4 is predicted by the modified Van't Hoff equation (eqn (10)) and is shown by the continuous line in Fig. 9. The force generation step has a high
H ° indicating that this step is an endothermic reaction and absorbs 72 ± 20 kJ mol-1 of heat at 25°C (Table 4). This observation is consistent with a large force increase with a temperature increase in the range 5-30°C in rabbit soleus fibres (Fig. 2A) and in rabbit psoas fibres (Goldman et al. 1987; Bershitsky & Tsaturyan, 1992; Zhao & Kawai, 1994; Ranatunga, 1996). This is also consistent with solution studies of extracted proteins demonstrating that the actin-myosin binding is an endothermic reaction (Highsmith, 1977). In soleus STFs the force generation step accompanies a large entropy increase corresponding to 247 ± 66 J mol-1 K-1 at 25°C (Table 4). Therefore, it can be concluded that the force generation step is an entropy-driven reaction. This is consistent with the conclusion obtained in psoas fibres (Table 4).
![]() |
View larger version [in this window] [in a new window] |
|
|
The continuous line represents the theoretical curve with best fit parameters. The data are the same as K4 in Fig. 8.
| ||
Our previous studies on psoas that demonstrated large positive
H ° and
S °, and large negative
Cp of the force generation step led us to conclude that the hydrophobic interaction between actin and myosin underlies force generation (Murphy et al. 1996). We found that this mechanism is also true in soleus STFs (Table 4). Water molecules are structured around hydrophobic amino acid residues, and this water structure is removed on the formation of the hydrophobic interaction, which accounts for a large reduction of the heat capacity change (
Cp). Based on the crystallographic structure of myosin and actin it was proposed that force is generated as a result of the stereospecific interaction occurring on myosin as a helix (G516-F542)-loop (P543-T546)-helix (D547-H558) motif, which comes close to actin residues A144-T148 and I341-Q354 (Rayment et al. 1993a; Rayment, 1996). The stereospecific interaction is primarily hydrophobic, and residues P529, M530, I535, M541, F542, P543 of myosin, and residues A144, I341, I345, L349, F352 of actin are involved (Rayment et al. 1993a). Mutational studies were carried out to examine the structure-function relationship of these hydrophobic residues. Miller et al. (1996) demonstrated the importance of I341 of actin in myosin binding, actin-activated S1 ATPase, and actin filament sliding in the in vitro motility assay. Additional hydrophobic interaction between segment R405-K415 of myosin and segment P332-E334 of actin is also suggested (Rayment et al. 1993b). It is interesting to note that these hydrophobic residues are conserved between type I and type IIb fibres (Weiss et al. 1999), and that the
H ° and the
S ° values associated with the force generation step in the two fibre types are similar (Table 4). The formation of these hydrophobic interactions probably corresponds to step 4, because a large enthalpy increase, a large entropy increase, and a large heat capacity decrease are observed with this step.
An additional mechanism to account for the large change in
H °,
S ° and
Cp is that the myosin cleft may close with force generation. Based on the thermodynamic parameters in rabbit psoas fibres, Murphy et al. (1996) calculated the burial of surface area to be 80 nm2, which is larger than the presumed actin-myosin interface, and hence other surface area must be involved. A closure of the cleft between the upper and lower 50 kDa domains is a likely candidate for the burial of a large surface area represented by a large negative
Cp. Another candidate that cannot be ruled out is the conformational change in actin.
As discussed, the equilibrium constant of step 4 (K4) increases with temperature up to 30°C. This will result in a significant increase in the number of cross-bridges in AM*DP and AM*D states. Because both AM*DP and AM*D states bear tension, tension will increase as a result of the temperature increase. The good fit of the theoretical curve (Fig. 2A) to the isometric tension data is consistent with the mechanism that more tension-generating cross-bridges are formed at higher temperature, and that the tension per cross-bridge remains the same when temperature is changed. This mechanism is consistent with earlier studies carried out on rabbit psoas fibres (Zhao & Kawai, 1994) and on single molecules (Kawai et al. 2000; Kawaguchi & Ishiwata, 2000). As seen in Table 2, both the force generation step (k4) and its reversal step (k-4) in soleus STFs are more temperature sensitive than those in psoas. This finding is in accord with the observation that 2
b is more temperature sensitive in soleus than in psoas (Table 1).
Unlike the nucleotide binding step, which is sensitive to temperature, the Pi binding step (K5) is not very sensitive to temperature (Fig. 8). These observations are similar to our earlier observations which showed that K5 was little affected by osmotic compression (Zhao & Kawai 1993), BDM treatment, EMD treatment, or TnC extraction (Zhao et al. 1996, and see reference therein), whereas the nucleotide binding step was significantly affected by them. These results indicate that the Pi release channel of myosin is controlled differently from the nucleotide binding site. This mechanism is consistent with the hypothesis that myosin is a 'back door' enzyme in which Pi is released from the opposite side of the ATP binding site (Yount et al. 1995).
| REFERENCES |
|---|
|
|
|---|
| BARCLAY C. J., CONSTABLE, J. K. & GIBBS, C. L. (1993). Energetics of fast- and slow-twitch muscles of the mouse. Journal of Physiology 472, 61-80 | [Abstract] |
| BERSHITSKY S. Y. & TSATURYAN, A. K. (1992). Tension responses to joule temperature jump in skinned rabbit muscle fibres. Journal of Physiology 447, 425-448 | [Abstract] |
| CROW M. T. & KUSHMERICK, M. J. (1982). Chemical energetics of slow- and fast-twitch muscles of the mouse. Journal of General Physiology 79, 147-166 | [Abstract] |
| CURTIN N. A. & WOLEDGE, R. C. (1978). Energy changes and muscular contraction. Physiological Review 58, 690-761 | |
| DANTZIG J., GOLDMAN, Y., MILLAR, N. C., LACKTIS, J. & HOMSHER, E. (1992). Reversal of the cross-bridge force-generation transition by the photogeneration of phosphate in rabbit psoas muscle fibers. Journal of Physiology 451, 247-278 | [Abstract] |
| FORD L. E., HUXLEY, A. F. & SIMMONS, R. M. (1977). Tension responses to sudden length change in stimulated frog muscle fibres near slack length. Journal of Physiology 269, 441-515 | [Medline] |
| FORTUNE N. S., GEEVES, M. A. & RANATUNGA, K. W. (1991). Tension responses to rapid pressure release in glycerinated rabbit muscle fibers. Proceedings of National Academy of Sciences of the USA 88, 7323-7327 | |
| FRANK H. S. & EVANS, M. W. (1945). Free volume and entropy in condensed system. III. Entropy on binary liquid mixtures; partial molar entropy in dilute solutions; structure and thermodynamics in aqueous electrolytes. Journal of Chemistry and Physics 13, 507-532. | |
| GALLER S. & HILBER, K. (1998). Tensions/stiffness ratio of skinned rat skeletal muscle fibre types at various temperatures. Acta Physiologica Scandinavica 162, 119-126 | [Medline] |
| GOLDMAN Y. E., MCCRAY, J. A. & RANATUNGA, K. W. (1987). Transient tension changes initiated by laser temperature jumps in rabbit psoas muscle fibres. Journal of Physiology 392, 71-95 | [Abstract] |
| HIGHSMITH S. (1977). The effects of temperature and salts on myosin subfragment-1 and F-actin association. Archives of Biochemistry and Biophysics 180, 404-408 | [Medline] |
| HUXLEY A. F. & SIMMONS, R. M. (1971). Proposed mechanism of force generation in striated muscle. Nature 233, 533-538 | [Medline] |
| KAWAI M. & BRANDT, P. W. (1980). Sinusoidal analysis: a high resolution method for correlating biochemical reactions with physiological processes in activated skeletal muscles of rabbit, frog and crayfish. Journal of Muscle Research and Cell Motility 1, 279-303 | [Medline] |
| KAWAI M. & HALVORSON, H. R. (1991). Two step mechanism of phosphate release and the mechanism of force generation in chemically skinned fibers of rabbit psoas muscle. Biophysical Journal 59, 329-342 | [Abstract] |
| KAWAI M., KAWAGUCHI, K., SAITO, M. & ISHIWATA, S. (2000). Temperature change does not affect force between single actin filaments and HMM from rabbit muscles. Biophysical Journal 78, 3112-3119 | [Abstract/Full Text] |
| KAWAI M., SAEKI, Y. & ZHAO, Y. (1993). Crossbridge scheme and the kinetic constants of elementary steps deduced from chemically skinned papillary and trabecular muscles of the ferret. Circulation Research 73, 35-50 | [Abstract] |
| KAWAI M. & SCHACHAT, F. H. (1984). Differences in the transient response of fast and slow skeletal muscle fibers: correlations between complex modulus and myosin light chain. Biophysical Journal 45, 1145-1151 | [Abstract] |
| KAWAI M. & ZHAO, Y. (1993). Cross-bridge scheme and force per cross-bridge state in skinned rabbit psoas muscle fibers. Biophysical Journal 65, 638-651 | [Abstract] |
| KAWAGUCHI K. & ISHIWATA, S. (2000). Temperature dependence of force, velocity, and processivity of single kinesin molecules. Biochemical and Biophysical Research Communications 272, 895-899 | [Medline] |
| KELLEY C. A., TAKAHASHI, M., YU, J. H. & ADELSTEIN, R. S. (1993). An insert of seven amino acids confers functional differences between smooth muscle myosins from the intestines and vasculature. Journal of Biological Chemistry 268, 12828-12854 | |
| KUSHMERICK M. J., MEYER, R. A. & BROWN, T. R. (1992). Regulation of oxygen consumption in fast- and slow-twitch muscle. American Journal of Physiology 263, C598-606 | [Medline] |
MCNALLY E. M., KRAFT, R. & BRAVO-ZEHNDER, M. H. (1989). Full-length rat and cardiac myosin heavy chain sequences. Journal of Molecular Biology 120, 665-671 |
|
| MEYER R. A., BROWN, T. R. & KUSHMERICK, M. J. (1985). Phosphorus nuclear magnetic resonance of fast- and slow-twitch muscle. American Journal of Physiology 248, C279-287 | [Medline] |
| MILLER C. J., DOYLE, T. C., BOBKOVA, E., BOSTEIN, D. & REISLER, E. (1996). Mutational analysis of the role of hydrophobic residues in the 338-348 helix on actin in actomyosin interactions. Biochemistry 35, 3670-3676 | |
| MOORE L. A., TIDYMAN, W. E., ARRIZUBIETA, M. J. & BANDMAN, E. (1992). Gene conversions within the skeletal myosin multigene family. Journal of Molecular Biology 223, 383-387 | [Medline] |
| MOSS R. L., DIFFEE, G. M. & GREASER, M. L. (1995). Contractile properties of skeletal muscle fibers in relation to myofibrillar protein isoforms. Reviews of Physiology, Biochemistry and Pharmacology 126, 1-63 | |
| MURPHY C. T. & SPUDICH, J. A. (1999). The sequence of the myosin 50-20K loop affects myosin's affinity for actin throughout the actin-myosin ATPase cycle and its maximum ATPase activity. Biochemistry 38, 3785-3792 | |
| MURPHY K. P., ZHAO, Y. & KAWAI, M. (1996). Molecular forces involved in force generation during skeletal muscle contraction. Journal of Experimental Biology 199, 2565-2571 | |
| PETTE D. & STARON, R. S. (1990). Cellular and molecular diversities of mammalian skeletal muscle fibers. Reviews of Physiology, Biochemistry and Pharmacology 116, 2-47 | |
| POTMA E. J., VAN GRAAS, I. A. & STIENEN, G. J. M. (1995). Influence of inorganic phosphate and pH on ATP utilization in fast and slow skeletal muscle fibers. Biophysical Journal 69, 2580-2589. | [Abstract] |
| PRIVALOV P. L. & GILL, S. J. (1988). Stability of protein structure and hydrophobic interaction. Advanced Protein Chemistry 39, 191-234 | |
| RANATUNGA K. W. (1984). The force-velocity relation of rat fast- and slow-twitch muscles examined at different temperatures. Journal of Physiology 351, 517-529 | [Abstract] |
| RANATUNGA K. W. (1989). Temperature sensitivity of isometric tension in glycerinated mammalian muscle fibres. In Muscle and Motility, ed. MARECHAL, G. & CARRARO, U., pp. 271-276. Intercept, Andover, UK | |
| RANATUNGA K. W. (1996). Endothermic force generation in fast and slow mammalian (rabbit) muscle fibers. Biophysical Journal 71, 1905-1913 | |
| RANATUNGA K. W. & WYLIE, S. R. (1983). Temperature-dependent transitions in isometric contractions of rat muscle. Journal of Physiology 339, 87-95 | [Abstract] |
| RAYMENT I. (1996). The structural basis of the myosin ATPase activity. Journal of Biological Chemistry 271, 15850-15853 | [Full Text] |
| RAYMENT I., HOLDEN, H. M., WHITTAKER, M., YOHN, C. B., LORENZ, M., HOLMES, K. C. & MILLIGAN, R. A. (1993a). Structure of the actin-myosin complex and its implications for muscle contraction. Science 261, 58-65 | [Medline] |
| RAYMENT I., RYPNIEWSKI, W. R., SCHMIDT-B€SE, K., SMITH, R., TOMCHICK, D. R., BENNING, M. M., WINKELMANN, D. A., WESENBERG, G. & HOLDEN, H. M. (1993b). Three-dimensional structure of myosin subfragment-1: a molecular motor. Science 261, 50-58 | [Medline] |
| ROME L. C., FUNKE, R. P., ALEXANDER, R. M., LUTZ, G., ALDRIDGE, H., SCOTT, F. & FREADMAN, M. (1988). Why animals have different muscle fibre types. Nature 335, 824-827. | [Medline] |
| ROVNER A. S., FREYZON, Y. & TRYBUS, K. M. (1997). An insert in the motor domain determines the functional properties of expressed smooth muscle myosin isoforms. Journal of Muscle Research and Cell Motility 18, 103-110 | [Medline] |
| SCHIAFFINO S. & REGGIANI, C. (1996). Molecular diversity of myofibrillar proteins: gene regulation and functional significance. Physiological Reviews 76, 371-424 | [Medline] |
| SPUDICH J. A. (1994). How molecular motors work. Nature 372, 515-518. | [Medline] |
| STEPHENSON D. G. & WILLIAMS, D. A. (1985). Temperature-dependent calcium sensitivity changes in skinned muscle fibres of rat and toad. Journal of Physiology 360, 1-12 | [Abstract] |
| SWEENEY H. L., ROSENFELD, S. S., BROWN, F., FAUST, L., SMITH, J., XIN, J., STEIN, L. A. & SELLERS, J. R. (1998). Kinetic tuning of myosin via a flexible loop adjacent to the nucleotide binding pocket. Journal of Biological Chemistry 273, 6262-6270 | [Abstract/Full Text] |
| TAYLOR E. W. (1979). Mechanism of actomyosin ATPase and the problem of muscle contraction. CRC Critical Reviews in Biochemistry 6, 103-164 | [Medline] |
| UYEDA T. Q., RUPPEL, K. M. & SPUDICH, J. A. (1994). Enzymatic activities correlate with chimeric substitution at the actin-binding face of myosin. Nature 368, 567-569 | [Medline] |
| WANG G., DING, W. & KAWAI, M. (1999). Does thin filament compliance diminish the cross-bridge kinetics? A study in rabbit psoas fibers. Biophysical Journal 76, 978-984 | [Abstract/Full Text] |
| WANG G. & KAWAI, M. (1996). Effects of MgATP and MgADP on the cross-bridge kinetics of rabbit soleus slow-twitch muscle fibers. Biophysical Journal 71, 1450-1461 | [Abstract] |
| WANG G. & KAWAI, M. (1997). Force generation and phosphate release steps in skinned rabbit soleus slow-twitch muscle fibers. Biophysical Journal 73, 878-894 | [Abstract] |
| WEISS A., SCHIAFFINO, S. & LEINWAND, L. A. (1999). Comparative sequence analysis of the complete human sarcomeric myosin heavy chain family: implications for functional diversity. Journal of Molecular Biology 290, 61-75 | [Medline] |
| YOUNT R. G., LAWSON, D. & RAYMENT, I. (1995). Is myosin a Ňback doorÓ enzyme? Biophysical Journal 68, suppl. 4, 44s-47s. | |
| ZHAO Y. & KAWAI, M. (1993). The effect of lattice spacing change on cross-bridge kinetics in chemically skinned rabbit psoas muscle fibers. II. Elementary steps affected by the spacing change. Biophysical Journal 64, 197-210 | [Abstract] |
| ZHAO Y. & KAWAI, M. (1994). Kinetics and thermodynamic studies of the cross-bridge cycle in rabbit psoas muscle fibers. Biophysical Journal 67, 1655-1668 | [Abstract] |
| ZHAO Y., SWAMY, G. P. M., HUMPHRIES, K. A. & KAWAI, M. (1996). The effect of partial extraction of troponin C on the elementary steps of the cross-bridge cycle in rabbit psoas muscle fibers. Biophysical Journal 71, 2759-2773 | [Abstract] |
Acknowledgements
We are grateful to Drs Kuda W. Ranatunga and Stefan Galler for their critical reading of the manuscript and useful suggestions. We would like to thank Ms Meghan L. Stevenson for her excellent technical assistance, and Dr Wei Ding for his help with computer graphics. This work was supported by grants IBN 93-18120, IBN 96-03858 and IBN 98-14441 from the National Science Foundation, and the Grant-in-Aid 99-50437N from the American Heart Association National Center. G.W. would like to dedicate this paper to Hui Yu and Zoe Q. Wang.
Corresponding author
M. Kawai: Department of Anatomy and Cell Biology, University of Iowa, Iowa City, IA 52242, USA.
Email: masataka-kawai{at}uiowa.edu
This article has been cited by other articles:
![]() |
J. M. Robinson, H. C. Cheung, and W. Dong The Cardiac Ca2+-Sensitive Regulatory Switch, a System in Dynamic Equilibrium Biophys. J., November 15, 2008; 95(10): 4772 - 4789. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. Rice, F. Wang, D. M. Bers, and P. P. de Tombe Approximate Model of Cooperative Activation and Crossbridge Cycling in Cardiac Muscle Using Ordinary Dif |