ARTICLE

Nonlinear summation of force in cat soleus muscle results primarily from stretch of the common-elastic elements

Abstract

The complex connective tissue structure of muscle and tendon suggests that forces from two parts of a muscle may not summate linearly. This study measured the nonlinear summation of force (Fnl) in whole cat soleus during isometric and ramp movements. In six anesthetized cats, the soleus was attached to a servomechanism to control muscle length and record force. The ventral roots were divided into two bundles, each innervating about half the soleus; thus the two parts could be stimulated alone or together. In all experiments, Fnl was small (<6% of maximum tetanic tension). Peak Fnl occurred during changes in muscle force, either as a result of imposed muscle movement or the onset or offset of a stimulus train. The data were fit to a model in which both parts of the muscle were assumed to stretch to a common elasticity. The servomechanism was programmed to compensate for reduced stretch of the common elasticity during partial compared with whole muscle activation. These compensatory movements showed how the model could account for most, but not all, of Fnl.

control of a muscle depends in part on the recruitment and derecruitment of motor units. Because of the complex connective tissue matrix within a muscle (14, 20, 26, 33, 35), the effective mechanical properties of newly recruited motor units may depend on how much of the muscle is already active and to what degree these active fibers strain neighboring passive fibers and associated connective tissue and tendon. Little is known about the mechanical interaction of muscle fibers, and a quantitative measure of this interaction is essential for understanding how whole muscle properties arise from its constituent motor units.

Fibers within a muscle may mechanically interact in several ways. Even with the simplest possible model of muscle, that is, independent parallel fibers attached to an external tendon, activation of one fiber will stretch the tendon and affect the length and velocity of neighboring fibers (see Fig.1A). Recent research has shown that the connective tissue structure of muscle is far more complex (Fig. 1B) (26, 33, 35). Proteins within the muscle fibers are capable of transmitting lateral forces between myofibrils and to the extracellular connective tissue matrix pervasive throughout the muscle. The links between neighboring fibers may prevent sarcomere length heterogeneity observed in isolated fibers, thus preventing damage and altering muscle force (8, 21). In some muscles, fibers do not run the length of the muscle but rather rely on this extracellular matrix and other muscle fibers to transmit force (18, 24). Thus neighboring fibers are not mechanically independent. Additional sources of interaction between muscle fibers occur in the aponeurosis and tendon, where cross-links likely exist between tendon fibrils (33). Furthermore, in pennate muscle, when part of the muscle contracts, the angle of pennation of the remaining muscle fibers is likely to be altered (10, 36). It is not presently known which of these factors are functionally important.

Fig. 1.

Fig. 1.Muscle and tendon models. A: simplified model of the common elasticity. Part A and part B refer to the muscle fibers in groups A and B, respectively. The elasticity ofpart A is attributed to the cross bridges within the muscle fibers, as well as elastic links within the fibers and between fibers and tendon, that act independently from part B. The common elasticity, K, refers to any elastic component that is stretched by the fibers in both parts A and B. B: detailed model with elastic links between fibers and tendons.


One measure of the mechanical interaction between parts of a muscle, and the measure that will be used in this study, is the change of force produced by one part when another part of the muscle is active. This is commonly referred to as nonlinear summation. In this paper, nonlinear summation of force (Fnl) will be defined as

Fnl(t)=FAB(t)FA(t)FB(t)Equation 1
FAB(t) is the force measured when both parts of the muscle are activated together. FA(t) and FB(t) are the force waveforms when part A and part B, respectively, are stimulated alone. Although the measure of Fnl is only an indirect measure of the strain throughout a muscle, it has the important advantage in that force is the key variable of interest in a systems' approach to motor control studies. Thus Fnl is a good functional representation of the total mechanical effects between parts of a muscle.

Nonlinear summation has been previously studied between parts of whole muscle during isometric (fixed end) contractions (4, 15) and between individual motor units during both isometric and slow ramp movements (7, 9, 27). In whole muscle studies, forces summed linearly during static isometric contractions. Linearity is surprising, considering the many ways fibers within a muscle can mechanically interact. This study expands on those studies by examining nonlinear summation during dynamic movements. In contrast to whole muscle, motor unit summation has been shown to be quite nonlinear (7, 9, 27). However, single motor units represent a small portion of a muscle; thus sources of nonlinearity, such as compression of passive muscle, may be insignificant when larger parts of the muscle are active. For this reason, it seemed important to study Fnl using larger active portions of the muscle.

The purpose of this study was to measure the nonlinear interaction between two parts of cat soleus during static and dynamic conditions. The soleus was experimentally divided into two large pseudo-motor units by splitting the ventral roots into two bundles. Fnl was measured during equal-duration isometric tetani to each part of the muscle, different-duration tetani to each part of the muscle, and during stretches and releases of the muscle. Finally, the lumped parameter model of Fig. 1A was fit to the data to determine to what extent the model of Fig. 1A can account for the nonlinearity. Each part of the muscle was assumed to have its own independent-series elastic element (representing the elasticity of the cross bridges plus the independent portion of the tendon), which in turn connect to a common-series elasticity. The physical location of the common elasticity is not restricted to the external tendon but rather represents any compliance stretched by both parts of the muscle. The model was shown to provide a good enough fit so more complex models may not be needed for most applications.

METHODS

The data were obtained from six cats (male and female). All surgical and experimental procedures conformed with the policies of Northwestern University and the National Institutes of Health. The surgical procedures are only briefly summarized here because a detailed account was recently presented (30). The cats were anesthetized with isoflurane during the surgical procedures and switched to pentobarbital sodium for data collection. The left hind legs were partially denervated and mounted in a rigid frame. The nerve and blood supplies to the soleus were preserved. The calcaneus was cut and attached to a servomechanism (custom device with a compliance of 0.01 mm/N), allowing the soleus to be moved by computer while simultaneously measuring muscle force. Before the tendon was cut, the foot was manually dorsiflexed and a thread was tied to the tendons of both the soleus and medial gastrocnemius. After the soleus was attached to the puller, realignment of the threads was defined as maximum physiological length (0 mm). Lo is defined in this paper as the length at which peak isometric tetanic tension (Po) occurs. (This corresponds to Loin Eq. 2.) The ventral roots were exposed via laminectomy and divided into two bundles, each part innervating roughly half of the soleus, and placed on separate hook electrodes so that each part could be stimulated independently. Stimulus isolation units were used, and the roots were stimulated using suprathreshold pulses 0.1 ms in duration and 2–10 V in amplitude. Stimulus trains were 100 Hz. The muscle force and length signals were generally sampled at 1 kHz, although for some quick stretches and releases sampling was at 10 kHz. Passive tension was always measured and subtracted from active tension.

Nonlinear summation was measured under three experimental conditions:1) isometric contractions with equal-duration tetani to both parts of the muscle, 2) isometric contractions with unequal-duration tetani, simulating recruitment of new motor units to an already partially active muscle, and 3) quick stretches and releases of the muscle, simulating a changing external load.Experiment 1 was selected because of its simplicity.Experiments 2 and 3 were selected because they were expected to produce the largest nonlinear summation. All three experiments share a common protocol. First, the whole muscle was stimulated by activating the ventral roots for part A andpart B simultaneously. Next, the ventral roots forpart A and part B were stimulated individually, and muscle force was measured. The length of the muscle-tendon complex was the same in each of the three trials. Nonlinear summation error, Fnl, was defined by Eq. 1. It was essential that there was no cross-stimulation of the ventral roots. This possibility was ruled out because tetanic tension from the two parts of the muscle summed linearly to 0.5% when the muscle was stimulated near maximum physiological length, even when the stimulating voltage was increased by a factor of two. It is also critical that the muscle did not fatigue during the times when parts A,B, and AB were measured. This was not a problem in cat soleus because the muscle is resistant to fatigue and because 1-min rest intervals were allowed between trials. In every group of trials (A, B, and AB), at least one waveform was remeasured to test for fatigue. In some experiments, the measurements were made several times and the order of stimulation was varied. Waveforms were generally repeatable with errors <0.5% of maximum tetanic tension.

Some of the data were analyzed by constructing length-tension curves during whole or partial muscle stimulation. The data were fit to Otten's equation (3, 23)

f(x)=PoexpxLoW+12.152711.28572.4649Equation 2
where f is force (in N), x is muscle length (in mm),Lo is the peak of the length-tension curve (in mm), and W is the width of the length tension curve (in mm). This empirical equation fit the soleus length tension data quite well by optimizing only three parameters: Lo,W, and Po. This procedure could detect small shifts in Lo.

To determine whether the common-elastic lumped parameter model of Fig.1A provided a reasonable account of the data, the puller was used to mimic stretch of the common elasticity. The common elasticity was assumed to be a linear Hookian element (K). When both halves of the muscle are active, the common elasticity will stretch by a distance (LAB) proportional to FAB

LAB(t)=FAB(t)/KEquation 3
When part A is stimulated alone, the muscle fibers will not shorten as much, since now the common elasticity will only be stretched by the force from part A
LA(t)=FA(t)/KEquation 4
Thus, for the muscle fibers of part A to shorten by the same amount when part A is stimulated alone, compared with when both parts of the muscle are active, the servomechanism needs to move by the difference betweenLAB(t) andLA(t). This movement can be approximated by
LB(t)=FB(t)/KEquation 5
The force waveform measured during activation of part B was divided by the estimated stiffness of the common elasticity. The resulting waveform is an approximation to the elastic stretch attributed to part B. It was used to drive the puller during activation of part A. Conversely, the force from part A was used to estimate the tendon stretch attributed to part A and used to drive the puller during activation of part B (Eq. 4). The summation of these two new waveforms was compared with the force from whole muscle stimulation to determine whether this could account for the nonlinearity. This procedure was repeated at both one-half and two times the initial estimate ofK.

Quick stretches (0.5 mm in 5 ms) were used to measure muscle stiffness. Within this distance and time, the muscle acts as a linear spring, so the change in the force waveform is equal to the change in the length waveform multiplied by the stiffness (s) (22). A computer program calculated the s that minimized the difference between the force and scaled length waveforms for the 5-ms period after the initiation of the stretch. This procedure provided a more accurate estimate of s than that obtained using a single-point measure 5 ms after stretch initiation.

The three-element model in Fig. 1A, coupled with the experimental stiffness measurements, can be used to determine the stiffness of the common elasticity. An algebraic solution was obtained by assuming that the common elasticity was to be linear over the range of forces measured (FA to FAB)

K=sAsB(sAsB)2+(sABsAsB)(sABsAsB)sABsAsBEquation 6
where K is the stiffness of the common-elasticity andsAB, sA, andsB are the experimentally measured stiffness of the muscle when parts AB, A, andB, respectively, are stimulated with the same stimulus train.

RESULTS

Isometric contractions with equal-duration tetani.

Figure 2 shows nonlinear summation during an isometric contraction near Lo. Tetani of 100 Hz, from 200 to 600 ms, were delivered to parts A,B, and AB together. The resulting force waveforms are shown in Fig. 2, top. The mathematical subtraction of these waveforms (Fnl) is shown in Fig. 2,middle. Overall Fnl(t) is quite small. The largest Fnl (0.65 N) is ∼3% of Po(25 N) and occurs during tension development (200–350 ms) and tension relaxation (650–850 ms). Note that Fnl during the force plateau (500–600 ms) was near zero.

Fig. 2.

Fig. 2.Top: typical example of nonlinear summation during an isometric contraction. The heavy line AB, thinnerline A, and dotted line B denote the forces from the whole muscle when ventral roots bundle A and B together, A alone, and B alone were stimulated, respectively (100 Hz, 0.2–0.6 s).Middle depicts nonlinear summation of force (Fnl) as defined by Eq. 1 in the text.Bottom shows the length of the muscle-tendon complex. In this muscle, length where peak isometric tetanic tension occurs (Lo) was −6.5 mm.


These results are consistent with the common elasticity model of Fig.1A. During the rise of force, the common elasticity stretches more when both parts of the muscle are active, resulting in a higher velocity of shortening and hence a lower force. During the plateau, fiber velocity is near zero; therefore, Fnl will depend primarily on altered fiber lengths and the resulting shift along the length-tension curve. Because this contraction occurred nearLo, Fnl is small. As shown below, steady-state Fnl increases at shorter or longer muscle lengths. During the fall in force, Fnl again becomes negative. This is probably due to a greater stretch velocity when both parts of the muscle relax together. Stretch during relaxation has been shown to cause muscle force to fall more rapidly (6).

To determine whether these results are quantitatively consistent with the common-elasticity model, the servomechanism was used to mimic stretch of the common elasticity. The common elasticity was estimated and used to compute movements for the servomechanism as described inEqs. 4 and 5. Preliminary experiments showed thatK was approximately equal to 20 N/mm. Values of Kat 10 and 40 N/mm were also tried. Figure3 shows the result using the same muscle as in Fig. 2. The solid light line shows Fnl when the estimated tendon stretch was based on a K of 20 N/mm. Note that the compensating movements of the servomechanism almost eliminated Fnl during the force-development and plateau phases. During the relaxation phase, substantial Fnl remained. The stiffK (40 N/mm) undercompensated. The compliant K (5 N/mm) overcompensated, except during the relaxation phase, where it minimized Fnl. Similar results were obtained on the one other cat subjected to this procedure.

Fig. 3.

Fig. 3.Servo movement to compensate for stretch of common elasticity during equal duration tetani. Heavy line is Fnlredrawn from Fig. 2. Light lines show puller compensation for 3 different values of K (10, 20, and 40 N/mm).


The effect of muscle length on Fnl was studied using the same equal-duration tetani. At different muscle lengths, the shape of Fnl changed. Figure 4 shows typical Fnl waveforms during isometric contractions with the same 100-Hz trains described above. During force-development and plateau region phases, Fnl(t) became increasingly negative at shorter muscle lengths. This was observed in all muscles. The response was more variable during the relaxation phase. During relaxation, most muscles showed the greatest negative Fnl at long lengths. In general, Fnl was less negative and, in some muscles, even positive at short muscle lengths.

Fig. 4.

Fig. 4.Example of nonlinear summation during isometric contractions with equal duration tetani at different muscle lengths. Same duration tetani as in Fig. 2 were used.


Partial activation of a muscle may not stretch the common elasticity as much as full activation. Therefore, when the muscle-tendon complex is held at a specific length, the muscle fibers will be longer during partial activation compared with during full activation; hence, a shift in the length-tension relationship to shorter lengths will occur during partial activation. To determine whether stretch of the common elasticity can be detected as a shift in the length-tension relationship, length-tension curves were constructed from data measured during whole and partial muscle stimulation. Force was measured every 2 mm, over a range from −16 to 0 mm, using 100-Hz trains 400 ms in duration. Force was measured on the plateau to minimize the dynamic effects of tendon stretch. The data were fit with Eq. 2. During partial muscle stimulation, Lo occurred at lengths 0.76 ± 0.38 mm shorter than that for whole muscle. The shift in Lo was statistically significant (P < 0.001, t-test, 1 tailed). When it is assumed that the tendon is linear in the region between FAand FAB, the shift can be used to estimate the stiffness of the common elasticity

K=(FABFA)/(LoABLoA)Equation 7
The data for all six cats provided an estimate of common elasticity of 11.7 ± 4.94 N/mm. Note that K estimated by using a shift in the length-tension curve is less than the Kthat provided the best fit when the puller was used to mimic tendon stretch. This is probably due to the nonlinear and history-dependent properties of the common elasticity (seediscussion). In addition, the very small shifts in the length-tension curves make this estimate susceptible to noise.

Isometric contractions with unequal-duration tetani.

The data above show that the maximum nonlinearities occurred during tension development and relaxation. Similar measurements were made during isometric contractions, except with unequal-duration tetani to each part of the muscle. The idea behind these experiments is that, if one part of a muscle contracts isometrically, stimulation of the second part will further stretch the common elasticity, transiently shortening the already active fibers. The force-velocity properties of muscle suggest this should decrease force from the active part. In contrast, during relaxation, stretch of the active fibers will lead to a positive velocity and increased force. Part A was stimulated at 100 Hz from 0.05 to 1.0 s, and part B was stimulated at 100 Hz from 0.20 to 0.50 s. Parts A and B were again stimulated together and then by themselves. Figure5 shows a typical result. Fnl(t) is quite small, with the peak <5% of Po. The peak Fnl occurred during the rise and fall of tension in part B. Note that, unlike the case with equal-duration tetani, here Fnl is positive during relaxation. The average peak Fnl(t) in six cats was 3.6 ± 1.8% of Po.

Fig. 5.

Fig. 5.Typical example of nonlinear summation during isometric contractions with different duration tetani. Part A was stimulated at 100 Hz from 0.05 to 1.0 s, and part B was stimulated at 100-Hz train from 0.2 to 0.5 s.


As before, the common-elasticity model was tested by using the servomechanism to mimic tendon stretch. Figure6 depicts this procedure. The puller movements based on the estimated tendon stretch (K = 20 N/mm) (bottom) and the force waveforms resulting from these movements (top) are both shown. The error was recalculated using these two waveforms and appears with the original Fnl(middle). Note that the error was substantially reduced. Three different values for K were tried on the same muscle, and the results are shown in Fig. 7. Compensatory movements based on the stiffest tendon (K =40 N/mm) were not large enough to minimize the error. Compensation based on K = 20 N/mm did a good job reducing the error during force development and the plateau region but only partially reduced the error during relaxation. Compensation using the most compliant value (K = 10 N/mm) overcorrected during force development but provided a good match during force relaxation. Overall, for the example shown in Fig. 7, the mean square error was reduced 55% by K = 40, 81% by K = 20, and 90% by K = 10 N/mm. Similar results were observed in muscles of all six cats. Generally, the error was most erratic during force relaxation.

Fig. 6.

Fig. 6.Demonstration of servo movement to compensate for stretch of common elasticity during different duration tetani. Same muscle and stimulus patterns as in Fig. 5 were used. See text for details.


Fig. 7.

Fig. 7.Servo movement to compensate for stretch of common elasticity during different duration tetani. Same muscle as in Figs. 5and 6 was used.


Because a single value for K did not optimally correct for the complete waveform, three different regions were analyzed:1) tension development (peak error between 200 and 300 ms),2) plateau region (mean error between 500 and 550 ms), and3) tension relaxation (peak error between 600 and 700 ms). The difference between the uncompensated Fnl and the compensated Fnl were measured at each of the three points and plotted as a function of K. The optimal value of Kwas selected by linear interpolation between points. The results based on six cats are as follows: 1) tension developmentK = 24.5 ± 10.6, 2) plateau regionK = 9.7 ± 18.1, and 3) tension relaxation K = 10.4 ± 4.8. These data are consistent with Fig. 7. Therefore, although the common-elasticity model qualitatively explains some of the error, the elasticity is not fully described by a linear elastic element.

Nonlinear summation during muscle stretch or release.

Nonlinear summation was also measured during imposed increases or decreases in muscle length during stimulation at 100 Hz. Initial tests showed that moderate ramps (±40 mm/s or less) did not produce Fnl larger than isometric contractions with unequal-duration tetani. However, at speeds of 100 mm/s or greater, larger values of Fnl were produced. A typical example is shown in Fig. 8 in which the muscle was shortened by a ramp of 2 mm in 20 ms. Note the largest Fnl (−6% of Po) occurred during the ramp. This is consistent with the common-elasticity model because more of the stretch will occur in the muscle fibers when only part of the muscle is active (the common elasticity is proportionally stiffer). Figure9 shows a similar result when the muscle was stretched by 2 mm. The average peak nonlinearity in six cats, for a positive or negative 2-mm ramp, was 6.5 ± 1.2% of Po. Attempts to use the servomechanism to compensate for stretch of the common elasticity (as above in Fig. 6) were not successful due to limitations of the servomechanism.

Fig. 8.

Fig. 8.Nonlinear summation during muscle movement. Similar to Fig. 2 except the stimulus is 100 Hz, 0.2–0.6 s, and the muscle is released with a ramp covering 2 mm in 20 ms.


Fig. 9.

Fig. 9.Nonlinear summation during muscle movement. Experiment was similar to that in Fig. 8, except that the muscle was stretched 2 mm in 20 ms.


Quick stretches were used to measure muscle stiffness when the whole muscle, part A, or part B was stimulated. Applying the model of Fig. 1A, these data enabled estimation of the common-elasticity stiffness using Eq. 6. On the basis of six cats, whole muscle stiffness was measured at 10.38 ± 3.16 N/mm with common elasticity of 22.25 ± 5.08 N/mm. By these calculations, 0.46 ± 0.07% of the total muscle compliance can be attributed to the common elasticity during 100-Hz stimulation. The remainder is attributed to compliance of the cross bridges and the independent portion of the elasticity.

DISCUSSION

Cat soleus was divided into two functional parts by splitting the ventral roots into two bundles, each producing approximately the same force. Nonlinear summation between these parts was measured under three different experimental conditions. In all experiments, nonlinear summation was small, <7% of Po. The largest value of Fnl occurred during changes in muscle force, either due to activation or deactivation of part of the muscle or by imposed movements on the whole muscle.

There are surprisingly few studies that have examined the summation of force in muscle. This may result, on one hand, from viewing muscle fibers as parallel independent force generators, and, on the other hand, from knowing that the actual anatomy is so complex that a completely accurate anatomic model is out of reach. Early studies concerned with polyneuronal innervation of muscle fibers examined summation when the muscle was stimulated by the ventral roots (4,15). These investigators found twitch summation less that linear. Brown and Mathews (4) correctly attributed this to stretch of the series elasticity. Rack and Westbury (29) used distributed stimulation to study the mechanical properties of cat soleus at low stimulation rates. They divided the ventral roots into five bundles and asynchronously stimulated each bundle so total muscle force was quite smooth compared with the unfused tetanus of each bundle. They noted that the average muscle force from the asynchronous stimulation was greater than the sum of each bundle stimulated separately. They speculated that the asynchronous stimulation reduced internal shortening of the fibers, preventing the breaking and reforming of the cross bridges, thus producing more tension. These studies are consistent with the data from this study.

Recent studies have demonstrated the nonlinear summation of force from individual motor units (7, 9, 11, 27). The underlying cause of the nonlinearity is not fully understood. It is likely caused by friction or elastic links between the relatively small number of active fibers in the motor units and the surrounding mass of passive fibers or other tissue. In this study, the ventral roots were divided into two bundles to form two large pseudo-motor units. Because the territory of a single motor unit can span almost the whole soleus (5), it can be assumed that the fibers for part A were randomly distributed between the fibers for part B. Thus nonlinear summation similar to that shown in motor units might have been observed in this study. However, the nature of the nonlinearities observed was much different. In the motor unit studies, activation of a second motor unit often increased the force from the first, even after stimulation stopped in the second unit (7). Powers and Binder (27) showed that the nonlinearities diminished if the muscle was moved. Stretch and relaxation effects similar to those shown in Fig. 5 were not observed. Therefore, it seems that the cause of the nonlinearities in motor units relates more to the large mass of passive tissue and not to the stretch of the common elastic elements. Thus nonlinearities measured between single motor units probably do not have relevance to this study.

The angle of pennation of the muscle fibers may contribute to nonlinear summation, but it does not appear to be a major factor in this study. In a pennate muscle, such as the soleus, as more of the muscle contracts, the tendon stretches and the angle of pennation increases. The resulting force from the muscle is believed to decrease by the cosine of the angle of pennation. Therefore, pennation effects would be expected to produce a negative Fnl even during constant force contractions. In this study, the steady-state Fnl was small (plateau in Figs. 2, 4, 5, 8, and 9) and at least partially accounted for by the shift in the length-tension curve. Therefore, pennation effects do not account for substantial nonlinearity in cat soleus.

The mechanical links between fibers within a muscle may contribute to nonlinear summation of force. Recent research has shown that the force transmission in muscle is complex and incompletely understood (14, 26, 33, 35). It is generally agreed that cross-bridge connections between actin and myosin form the basis for active force generation (16). These cross bridges are an important part of the total muscle compliance. Titin also appears to be a significant load-bearing structure. It connects the myosin filaments to the Z disk and is probably responsible for most of the passive tension within a muscle (12, 13, 19). Other cytoskeletal proteins probably play roles in maintaining lattice structure (26). Intermediate filaments, desmin and skelemin, probably cross-link the myofibrils, keeping them in register (17). Regularly spaced costameres are probably the structures that transmit force between the myofibrils and the extracellular space (25). Thus, within the sarcolemma, the actin-myosin cross bridges and titin support longitudinal force transmission and other proteins are involved in transverse transmission of force. Outside the muscle fibers, there is a well-developed endomysial connective tissue matrix capable of transmitting muscle force (33, 34). Muscle fibers taper at their ending. Although there are specialized structures where they insert to the tendon plate, because the fibers always taper, the connective tissue matrix is likely to transmit some of the force. Furthermore, recent evidence shows that many muscles have fibers that do not run the length of the muscle (18, 20, 24). These fibers by necessity transmit their force to this endomysial connective tissue matrix. Street (32) recently demonstrated this transmission in frog fibers, which have a weaker endomysial matrix than mammalian fibers. Thus force transmission in skeletal muscle is complex, far more complex than the common view of independent parallel fibers. It appears that muscle fibers are tightly linked to each other, which may affect the way force from parts of a muscle sum together. Despite the links between muscle fibers, the simple common-elasticity model of Fig. 1A fits the data reasonably well. A portion of the common elasticity measured in this study may be attributable to the links between muscle fibers.

Aponeurosis and tendon characteristics must also be considered when examining nonlinear summation. The external tendons of muscles have been studied in vitro by attaching them to pullers and studying their stress-strain characteristics. Generally, their properties have proven to be complex. The stress-strain relationship shows a toe region, where the tendon is quite compliant, followed by either an exponential or piecewise linear relationship until the stress becomes high enough that the tendon fails (36). The common-elasticity model used in this study assumed a linear K that may have contributed to the error. Studies have shown that tendons, when stretched to a length and held, exhibit creep over time, stretching with at least two different time constants (1, 36). In cycles of stretch and release, tendons show marked hysteresis, with release exhibiting less stiffness. This may explain why different values of K were needed to match the different regions of Fnl in Fig. 7. The tendon may be stiffer during stretch than during release. The problem with studying excised tendons is the difficulty attaching a cut tendon to a puller and knowing what part of the tendon is held firm and what part is allowed to stretch. In a careful study by Bennett et al. (2), attempts were made to determine whether the stress-strain relationship is constant for tendons of different mammalian species. They concluded that the difference between species was small compared with the error of measurement. They estimated compliance measurements were in error by a factor of two due to problems attaching the tendon to a measuring device. Scott and Loeb (31) measured the mechanical properties of aponeurosis and tendon in cat soleus using sonomicrometers and concluded that the mechanical properties were the same. This fits with the conclusion of Bennet et al., although there is probably at least as much problem attaching sonomicrometers to the tendon as Bennet et al. encountered in securing the cut end of a tendon. Thus uncertainty remains as to tendon and aponeurosis properties.

In this study, when the common-elasticity model was fit to the data, different values for K were derived, which were dependent on whether the movements were dynamic or static. The shift in the length-tension curve and compensation for tendon stretch during the plateau regions yielded an estimate of K of ∼10 N/mm. Estimates of K during tension development and quick stretches gave a K value >20 N/mm. These methods of estimating K are advantageous because they are noninvasive and because the parameter can be determined in a way muscle is normally used. The disadvantage is that these methods are not direct measures of tendon strain; rather, they use the measured nonlinearity of summation to estimate common-elastic stretch. Small errors in measurement lead to large errors in stiffness estimates. This may explain why the common-elastic stiffness estimated via a shift in the length-tension curve was one-half that of the other two methods. Alternately, as mentioned above, tendon stiffness may be greater for dynamic stretch and thus a lower value during static length-tension measurements.

Estimates of common-elastic stiffness in this study agree closely with tendon stiffness measurements made using other techniques. Morgan (22) developed the alpha method to measure all the elasticity in a muscle other than that attributed to the cross bridges. In cat soleus, the total tendon stiffness was estimated to be 16.9 N/mm (22). This method is likely to measure some elastic elements that are not shared between the fibers and, therefore, not measured in this study. Rack and Westbury (28) used muscle spindles to estimate what they called the “entire tendon” in cat soleus. At 11 N, they estimated the tendon compliance to be 25 N/mm. This is close to the value of 24.5 N/mm obtained by mimicking tendon stretch during tension development and 22.5 N/mm estimated using quick stretches and algebraic calculations.

In this study, the common elasticity was measured during tetanic stimulation. Thus it was measured at relatively high force levels, ranging from about 50% of maximal tetanic force when half the muscle was stimulated to 100% of maximal tetanic force when both halves were active. As mentioned above, tendon properties are known to be nonlinear with greater compliance at low force levels (36). The elastic elements comprising the common elasticity may have similar properties. Thus, during natural activation at low firing rates, the common elasticity will probably be more compliant. However, the force during partial activation will also be less; thus fiber shortening, and therefore nonlinear summation, will not necessarily be greater than that measured in this study. Further work examining partial activation at subtetanic rates is necessary.

In conclusion, 1) the interaction (nonlinear summation of force) between two parts of cat soleus muscle was small. During isometric contractions the Fnl was <3% of Po.2) The largest nonlinearities occurred during dynamic changes in muscle force, resulting from the onset or offset of a stimulus train or movement of the muscle. Still, the Fnlwas <7% of Po. 3) Most of the nonlinear interaction can be explained by a simple lumped parameter model in which the interaction occurs through the stretch of a shared elastic elements (common elasticity). The stiffness of the common elasticity at maximum tetanic force was estimated at two times the stiffness of the whole muscle.

This work was supported by National Institute of Arthritis and Musculoskeletal and Skin Diseases Grant AR-34382

FOOTNOTES

  • Address for reprint requests and other correspondence: T. G. Sandercock, Dept. of Physiology, M211, Ward 5–295, Northwestern Univ. School of Medicine, 303 E. Chicago Ave., Chicago, IL 60611 (E-mail: ).

  • The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked “advertisement” in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

REFERENCES

  • 1 Alexander RM.Mechanics of skeleton and tendon.Handbook of Physiology. The Nervous System. Motor Control.1981Am. Physiol. SocBethesda, MD, sect. 1, pt. 1, chapt. 2, p. 17–42.
    Google Scholar
  • 2 Bennett MB, Ker RF, Dimery NJ, Alexander RM.Mechanical properties of various mammalian tendons.J Zool Lond2091986537548
    Crossref | ISI | Google Scholar
  • 3 Brown IE, Scott SH, Loeb GE.Mechanics of feline soleus. II. Design and validation of a mathematical model.J Muscle Res Cell Motil171996221233
    Crossref | ISI | Google Scholar
  • 4 Brown MC, Mathews PBCAn investigation into the possible existence of polyneuronal innervation of individual skeletal muscle fibers in certain hind-limb muscles of the cat.J Physiol (Lond)1511960436457
    Crossref | Google Scholar
  • 5 Burke RE.Motor units: anatomy, physiology, and functional organization.Handbook of Physiology. The Nervous System. Motor Control.1981Am. Physiol. SocBethesda, MD, sect. 1, pt. 1, chapt. 10, p. 345–422.
    Google Scholar
  • 6 Caputo C, Edman KA, Lou F, Sun YB.Variation in myoplasmic Ca2+ concentration during contraction and relaxation studied by the indicator fluo-3 in frog muscle fibres.J Physiol (Lond)4781994137148
    Crossref | Google Scholar
  • 7 Clamann HP, Schelhorn TB.Nonlinear force addition of newly recruited motor units in the cat hindlimb.Muscle Nerve11198810791089
    Crossref | PubMed | ISI | Google Scholar
  • 8 Edman KA, Reggiani C.Redistribution of sarcomere length during isometric contraction of frog muscle fibres and its relation to tension creep.J Physiol (Lond)3511984169198
    Crossref | Google Scholar
  • 9 Emonet-Denand F, Laporte Y, Proske U.Summation of tension in motor units of the soleus muscle of the cat.Neurosci Lett1161990112117
    Crossref | ISI | Google Scholar
  • 10 Gans C, Gaunt AS.Muscle architecture in relation to function.J Biomech2419915365
    Crossref | ISI | Google Scholar
  • 11 Heckman CJ, Weytjens JL, Loeb GE.Effect of velocity and mechanical history on the forces of motor units in the cat medial gastrocnemius muscle.J Neurophysiol68199215031515
    Link | ISI | Google Scholar
  • 12 Horowits R, Kempner ES, Bisher ME, Podolsky RJ.A physiological role for titin and nebulin in skeletal muscle.Nature3231986160164
    Crossref | PubMed | ISI | Google Scholar
  • 13 Horowits R, Podolsky RJ.The positional stability of thick filaments in activated skeletal muscle depends on sarcomere length: evidence for the role of titin filaments.J Cell Biol105198722172223
    Crossref | PubMed | ISI | Google Scholar
  • 14 Huijing PA.Muscle as a collagen fiber reinforced composite: a review of force transmission in muscle and whole limb.J Biomech321999329345
    Crossref | PubMed | ISI | Google Scholar
  • 15 Hunt CC, Kuffler SW.Motor innervation of skeletal muscle: multiple innervation of individual muscle fibres and motor unit function.J Physiol (Lond)1261954293303
    Crossref | Google Scholar
  • 16 Huxley AF.Muscle structure and theories of contraction.Prog Biophys Biophys Chem71957257318
    Crossref | Google Scholar
  • 17 Lazarides E.Intermediate filaments as mechanical integrators of cellular space.Nature2831980249256
    Crossref | PubMed | ISI | Google Scholar
  • 18 Loeb GE, Pratt CA, Chanaud CM, Richmond FJ.Distribution and innervation of short, interdigitated muscle fibers in parallel-fibered muscles of the cat hindlimb.J Morphol1911987115
    Crossref | PubMed | ISI | Google Scholar
  • 19 Magid A, Law DJ.Myofibrils bear most of the resting tension in from skeletal muscle.Science230198512801282
    Crossref | PubMed | ISI | Google Scholar
  • 20 Monti RJ, Roy RR, Hodgson JA, Edgerton VR.Transmission of forces within mammalian skeletal muscles.J Biomech321999371380
    Crossref | PubMed | ISI | Google Scholar
  • 21 Morgan DL.New insights into the behavior of muscle during active lengthening.Biophys J571990209221
    Crossref | PubMed | ISI | Google Scholar
  • 22 Morgan DL.Separation of active and passive components of short-range stiffness of muscle.Am J Physiol Cell Physiol2321977C45C49
    Link | ISI | Google Scholar
  • 23 Otten E.Concepts and models of functional architecture in skeletal muscle.Exerc Sport Sci Rev16198889137
    Crossref | Google Scholar
  • 24 Ounjian M, Roy RR, Eldred E, Garfinkel A, Payne JR, Armstrong A, Toga AW, Edgerton VR.Physiological and developmental implications of motor unit anatomy.J Neurobiol221991547559
    Crossref | PubMed | Google Scholar
  • 25 Pardo JV, Siliciano JD, Craig SW.A vinculin-containing cortical lattice in skeletal muscle: transverse lattice elements (“costameres”) mark sites of attachments between myofibrils and sarcolemma.Proc Natl Acad Sci USA80198310081012
    Crossref | PubMed | ISI | Google Scholar
  • 26 Patel TJ, Lieber RL.Force transmission in skeletal muscle: from actomyosin to external tendons.Exerc Sport Sci Rev251997321363
    Crossref | PubMed | Google Scholar
  • 27 Powers RK, Binder MD.Summation of motor unit tensions in the tibialis posterior muscle of the cat under isometric and nonisometric conditions.J Neurophysiol66199118381846
    Link | ISI | Google Scholar
  • 28 Rack PM, Westbury DR.Elastic properties of the cat soleus tendon and their functional importance.J Physiol (Lond)3471984479495
    Crossref | Google Scholar
  • 29 Rack PMH, Westbury DR.The effects of length and stimulus rate on tension in the isometric cat soleus muscle.J Physiol (Lond)2041969443460
    Crossref | Google Scholar
  • 30 Sandercock TG, Heckman CJ.Force from cat soleus muscle during locomotor-like movements: experimental data versus Hill-type model predictions.J Neurophysiol77199715381552
    Link | ISI | Google Scholar
  • 31 Scott SH, Loeb GE.The mechanical properties of the aponeurosis and tendon of the cat soleus muscle during whole-muscle isometric contractions.J Morphol22419957386
    Crossref | PubMed | ISI | Google Scholar
  • 32 Street SF.Lateral transmission of tension in frog myofibers: a myofibrillar network and transverse cytoskeletal connections are possible transmitters.J Cell Physiol1141983346364
    Crossref | PubMed | ISI | Google Scholar
  • 33 Trotter JA.Functional morphology of force transmission in skeletal muscle.Acta Anat (Basel)1461993205222
    Crossref | PubMed | Google Scholar
  • 34 Trotter JA, Purslow PP.Functional morphology of the endomysium in series fibered muscles.J Morphol2121992109122
    Crossref | PubMed | ISI | Google Scholar
  • 35 Trotter JA, Richmond FJR, Purslow PP.Functional morphology and motor control of series-fibered muscle.Exerc Sport Sci Rev231995167213
    Crossref | Google Scholar
  • 36 Zajac FE.Muscle and tendon: properties, models, scaling, and application to biomechanics and motor control.Crit Rev Biomed Eng171989359411
    PubMed | ISI | Google Scholar