QUICK SEARCH:   [advanced] Author: Keyword(s):
Home     Help     Feedback     Subscriptions     Archive     Search     Table of Contents

First published online March 2, 2007
Journal of Experimental Biology 210, 1092-1108 (2007)
Published by The Company of Biologists 2007
doi: 10.1242/jeb.02729

This Article Summary Figures Only Full Text (PDF) Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via HighWire Citing Articles via Google Scholar Google Scholar Articles by Guschlbauer, C. Articles by Büschges, A. Search for Related Content PubMed PubMed Citation Articles by Guschlbauer, C. Articles by Büschges, A.

The extensor tibiae muscle of the stick insect: biomechanical properties of an insect walking leg muscle

Christoph Guschlbauer, Hans Scharstein and Ansgar Büschges^*

Zoological Institute, University of Cologne, Weyertal 119, 50923 Cologne, Germany

^* Author for correspondence (e-mail: Ansgar.Bueschges{at}uni-koeln.de)

Accepted 23 January 2007

	Summary

TOP Summary Introduction Materials and methods Results Discussion References

 
We investigated the properties of the extensor tibiae muscle of the stick insect (Carausius morosus) middle leg. Muscle geometry of the middle leg was compared to that of the front and hind legs and to the flexor tibiae, respectively. The mean length of the extensor tibiae fibres is 1.41±0.23 mm and flexor fibres are 2.11±0.30 mm long. The change of fibre length with joint angle was measured and closely follows a cosine function. Its amplitude gives effective moment arm lengths of 0.28±0.02 mm for the extensor and 0.56±0.04 mm for the flexor. Resting extensor tibiae muscle passive tonic force increased from 2 to 5 mN in the maximum femur–tibia (FT)-joint working range when stretched by ramps.

Active muscle properties were measured with simultaneous activation (up to 200 pulses s^–1) of all three motoneurons innervating the extensor tibiae, because this reflects most closely physiological muscle activation during leg swing. The force–length relationship corresponds closely to the typical characteristic according to the sliding filament hypothesis: it has a plateau at medium fibre lengths, declines nearly linearly in force at both longer and shorter fibre lengths, and the muscle's working range lies in the short to medium fibre length range. Maximum contraction velocity showed a similar relationship. The force–velocity relationship was the traditional Hill curve hyperbola, but deviated from the hyperbolic shape in the region of maximum contraction force close to the isometric contraction.

Step-like changes in muscle length induced by loaded release experiments characterised the non-linear series elasticity as a quadratic spring.

Key words: pinnate insect muscle, muscle properties, contraction dynamics

	Introduction

TOP Summary Introduction Materials and methods Results Discussion References

 
When animals move, a neurally generated motor output activates the musculo-skeletal system. This precise motor action is determined by the active and passive biomechanical properties of the muscles. Depending on their contraction dynamics, muscles can respond to different temporal components of their neural inputs (e.g. Brezina et al., 2000; Ballantyne and Rathmayer, 1981; Bässler and Stein, 1996; Full, 1997; Josephson, 1993; Meyrand and Marder, 1991; Morris and Hooper, 1997) (reviewed in Hooper and Weaver, 2000; Morris et al., 2000). Understanding of how nervous systems generate motor behaviours requires investigation of muscle properties as well as neural activity. This is particularly important in complex motor systems that function via the concerted action of multiple muscle groups, e.g. for terrestrial locomotion. In these systems the mechanical arrangement and activation of muscles can make synergistic muscles perform different roles. In the cockroach, for example, one of the two leg extensor muscles acts as a motor and the other as a brake (Ahn and Full, 2002). For organisms in which the neural component is already quite well understood, detailed investigation of muscles will allow a better characterisation of the roles that neural and muscular properties play in movement generation. One such organism is the stick insect Carausius morosus (Büschges, 2005; Orlovsky et al., 1999).

Motor output of the stick insect leg muscle control system is the result of a complex interaction between local sensory feedback, central neural networks governing the individual leg joints, and coordinating signals between the legs (e.g. Bässler and Büschges, 1998; Büschges, 2005; Dürr et al., 2004). The femur–tibia (FT)-joint is the functional knee-joint of the insect leg, and detailed information has been gathered with respect to morphological organization (Bässler, 1967) and the motoneuronal innervation pattern of the muscles involved, the flexor and the extensor tibiae (Bässler et al., 1996; Bässler and Storrer, 1980; Debrodt and Bässler, 1989; Debrodt and Bässler, 1990). In addition, motor output during walking movements (Bässler, 1993; Büschges et al., 1994; Fischer et al., 2001) and some aspects of the neural control, including the activity of the central premotor networks (Bässler, 1993; Büschges, 1995; Büschges et al., 2004; Driesang and Büschges, 1996), are known. This information has been sufficient to successfully construct a neuro-mechanical simulation of the stepping stick insect (Ekeberg et al., 2004). However, the control of movement amplitude, i.e. the activation level of leg motoneuron pools and muscles, is only poorly understood, partly because neural and muscular properties interact at this level in movement generation (Blickhan et al., 2003; Brezina and Weiss, 2000; Chiel and Beer, 1997). In order to address this issue, it is necessary to investigate how the detailed aspects of neural control, e.g. changes in motoneuron activity, affect muscle activation and movements generated thereby, and the general working range of the neuron-to-movement transformation for stick insect leg muscles. Given the well-known neural aspects of stick insect leg motor control, understanding of the FT-joint control would be greatly increased by a better understanding of the biomechanics of the joint's muscular system. It would be particularly interesting to understand the control of activation of this muscle, i.e. the properties of force production and contraction related to frequency in activity of the innervating tibial extensor motoneurons, including the single twitch to maximum force ratio.

The flexor and extensor tibiae muscles are pinnate, as is typical for arthropods, and control tibia movement for posture and locomotion (Bässler, 1983; Bässler, 1993). The fibres of the two leg muscles are innervated by multiple excitatory motoneurons; fast, semifast and slow motoneurons for the flexor tibiae and one fast (the FETi) and one slow motoneuron (the SETi) for the extensor tibiae (Bässler and Storrer, 1980). In addition, the extensor muscle fibres receive innervation from the common inhibitor (CI₁) motoneuron 1 (Bässler et al., 1996; Bässler and Stein, 1996; Bässler and Storrer, 1980) and muscle fibres of the flexor tibiae receive innervation from CI₂ and CI₃ (Debrodt and Bässler, 1990). In the extensor tibiae muscle there is a systematic proximal-to-distal shift between muscle fibres innervated by the FETi alone and those that are also innervated by SETi and CI₁. During walking in the stick insect middle leg, all three extensor motoneurons are activated maximally during leg swing (Büschges et al., 1994; Schmitz and Hassfeld, 1989). CI₁ activity switches off force production of dually innervated fibres (Bässler and Stein, 1996). Depending on the walking situation, e.g. when walking on a double treadwheel, extensor motoneurons, i.e. also FETi, can also be active during stance before the initiation of the swing phase, albeit at a reduced level (Büschges et al., 1994; Graham, 1985; Schmitz and Hassfeld, 1989).

In the present study we investigated the geometrical characteristics of the stick insect FT-joint and the static and dynamic properties of the two antagonistic tibial muscles, the extensor and flexor tibae, in the middle leg. Due to the small number of excitatory motoneurons innervating it, we decided to investigate the static and dynamic properties of force generation for the extensor muscle in more detail. We first focused on static properties, i.e. the length dependence of isometric force in the muscle. We then present data on the dynamic properties of the extensor muscle, including the relationship between contraction force and velocity (Hill curve), muscle series elasticity, and the dependence of the muscle parameters on stimulation frequency.

	Materials and methods

TOP Summary Introduction Materials and methods Results Discussion References

 
Experiments were carried out on 102 adult female stick insects of the species Carausius morosus Br. from a colony maintained at the University of Cologne. 10 animals of the same size as used in the experiment had an average body length of 77.1±2.28 mm and average mass of 940±70 mg. All experiments were performed under daylight conditions and at room temperature (20–22°C).

Muscle and fibre length measurements
The extensor and flexor muscles were exposed for length measurements by cutting a small window into the proximal and distal part of the femur. Muscle length was calculated under the microscope by determining the distance from the insertion of the most proximal fibre to an orientation mark (Fig. 1, shown for the extensor tibiae) and adding the distance from this mark to the insertion of the most distal fibre into the tendon. The tibia was moved on a plastic goniometer from 30° to 180° and length measurements were taken in 10° intervals. We considered this range (150°) as the maximum working range of both tibial muscles (Storrer, 1976; Cruse and Bartling, 1995). 90° was defined as the FT-joint angle at which both muscles are at their resting length for the stick insect (Friedrich, 1932; Storrer, 1976) and for Blaberus discoidalis (Full et al., 1998; Ahn and Full, 2002). Fibre length measurements were performed on muscles fixed in situ, with the joint at the 90° position, using 2.5% glutaraldehyde in phosphate buffer, pH 7.4 (Watson and Pflüger, 1994). Fibres were pulled from the proximal and medial parts of the femur because in these locations both muscles are primarily innervated by fast motoneurons (Bässler et al., 1996; Debrodt, 1980).

View larger version (10K): [in this window] [in a new window]   Fig. 1. Schematic representation of the geometrical arrangement of the femur–tibia joint in the stick insect middle leg. For details see text.

Experiments to determine the active and passive force–length characteristics were carried out as follows: muscle length was manually set to a FT-joint angle of 90° at the beginning, our definition of the muscle's resting length (l₀). It was then released 0.75 mm (0.5l₀) and subsequently stretched in coarse steps (mostly 0.15 mm) with sequencer-generated ramps of 0.05 mm s^–1 up to 0.75 mm beyond the muscle's resting length (1.5l₀). In some experiments we examined force development within the muscle's working range using small ramps (0.05 mm) within the range of 0.8l₀ to 1.2l₀ to obtain a more accurate screening. To investigate actively generated forces, the muscle was electrically stimulated with a paradigm of different stimulation frequencies at each length position after relaxation (for details, see Results, `Passive muscle forces').

Force–velocity curves (Hill curves) were obtained according to established procedures. Extensor muscles were stimulated to reach `steady-state' contraction under isometric conditions and then allowed to shorten under isotonic conditions against a variety of sequencer-generated counterforce levels. Muscle lengthening was accomplished by stretches while applying sequencer-generated load levels larger than the tetanical steady-state contraction force (for details, see Results, `Dynamics of the muscle contraction').

Electrical stimulation of motor axons
The thorax was opened dorsally and the gut, fat and connective tissue removed to expose the mesothoracic ganglion. A bipolar hook electrode was then placed under nerve nl3 (Marquardt, 1940). The nerve was crushed proximally with a forceps and isolated with white vaseline (Engelhard Arzneimittel GmbH & CoKG, Niederdorfelden, Germany). To measure the tension generated by the middle leg extensor tibiae muscle, we electrically stimulated the axons of its innervating motoneurons, the FETi, SETi and CI₁. These axons have different diameters (Bässler and Storrer, 1980), so the three motoneurons can be sequentially activated by increasing stimulation strength; FETi has the lowest threshold (Fig. 2Ai,Aiv), SETi the next highest (Fig. 2Aii,Aiv), and CI₁ the highest (Fig. 2Aiii,Aiv). The determination of the appropriate current pulse amplitude to use in the nerve stimulation was complicated by a conflict between (1) the desire to routinely stimulate all three motoneurons (FETi, SETi and CI₁) and (2) the desire to keep the current amplitude low enough that the nerve could be repeatedly stimulated over long periods without damage. This issue was even more difficult to resolve because the dissection required to enable extracellular recordings from the extensor nerve F2 in the distal femur close to the muscle (necessary to test whether all three axons are being stimulated) inevitably damaged some more distal muscle fibres, which are mostly dually innervated (Bässler et al., 1996), and considerably lengthened the dissection procedure. It was consequently not possible both to perform the long experiments on undamaged muscles whose data are reported here in the Results, and to check if all three motor axons were being stimulated in the same preparation.

View larger version (26K): [in this window] [in a new window]   Fig. 2. Isometric forces induced in the middle leg extensor tibiae muscle by electrical stimulation of nerve nl3 with different current amplitudes. In all panels the top trace is a stimulus monitor (note pulse height changes as stimulus amplitude was increased), the second trace is an extracellular recording of nerve nl3, and the third trace is muscle force. (Ai–Aiii) Sequential recruitment of FETi (Ai), FETi and SETi (Aii) and FETi, SETi and CI1 (Aiii) recorded in extensor leg nerve F2 in response to single stimuli. (Aiv) An enlarged version of the recordings, showing the sequential addition of new units (asterisks). 1 T=0.0023 mA. (Bi–Biv) F2-recordings and forces in response to a 50 Hz pulse train. (Bi) 75% of the pulses excited FETi and 25% FETi and SETi. (Bii) 50% of the stimuli elicited FETi and 50% FETi, SETi and CI1. (Biii) Recruitment of all three motor units with every pulse. Doubling the current amplitude (Biv) induced no further increase in force. In this experiment the SETi spikes were of larger amplitude than FETi spikes. This is uncommon and likely because nerve F2 was recorded very distally in the femur. In all panels the electrical disturbance in the nerve recording that coincides with the stimulus is a stimulus artifact, not an action potential (arrow in Ai). T, threshold.

These control experiments thus showed that in 90% of preparations, stimulating nl3 with a current amplitude 1.5 times greater than the FETi threshold resulted in either activation of all three motor axons or, if not, that the failure to activate SETi and CI₁ completely had a negligible effect on muscle force production. In the remaining 10% of preparations, the lack of the SETi and CI₁ only induced a modest decrease in measured muscle force. We therefore chose to set the amplitude of the current pulses to approximately 50% above threshold generating a visible twitch (Malamud and Josephson, 1991). Pulse trains of different frequencies were applied using a SPIKE2 sequencer program at intervals of at least 30 s to allow return to the rest state. Pulse duration was 0.5 ms in all experiments (Josephson, 1985; Stevenson and Josephson, 1990; Malamud and Josephson, 1991; Full et al., 1998; Ahn and Full, 2002).

Data achievement and evaluation
Data were recorded on a PC using a MICRO1401 A/D converter with SPIKE2-software (both Cambridge Electronic Design Limited, Cambridge, UK). For isometric force experiments, the influence of filament overlap on force generation was examined by stretching the muscle in ramps of different size over a range of 1.5 mm using a SPIKE2 sequencer program. The stimulation protocol was carried out at each muscle length. In these experiments the dual-mode lever system was used exclusively as a force transducer. The stiffness of the measuring system was measured by connecting the insect pin directly to the base of the platform, and it was greater than 7500 mN mm^–1; it can therefore be neglected in our measurements [for a thorough discussion of the influence of the compliance of the measuring device, see Jewell and Wilkie (Jewell and Wilkie, 1958)]. Twitch kinetic measurements included the time to peak force (T_max), time to 50% relaxation (T_50off) and time to 90% relaxation (T_90off), all calculated relative to the force onset. For isotonic force experiments, the influence of load on contraction velocity and muscle series elasticity was determined by application of different force levels on the lever arm during tetanus using a SPIKE2 sequencer script. Custom SPIKE2 script programs were written for most of the data analysis. Plotting, curve fitting, and error evaluation were performed in ORIGIN (Microcal. Software Inc., Northampton, MA, USA).

Statistics
Mean values were compared using a modified t-test (Dixon and Massey, 1969). Means and samples were regarded as significantly different at P<0.05. The following symbols show the level of statistical significance: (–) not significant; *0.01<P<0.05; **0.001<P0.01; ***P0.001. In the text N gives the number of experiments or animals while n gives the sample size. All data were calculated as means ± s.d.

View larger version (12K): [in this window] [in a new window]   Fig. 3. (A) Relationship between muscle resting length and femur length in front leg (F; squares), middle leg (M; circles) and hind leg (H; triangles). Filled symbols, extensor tibiae; open symbols, flexor tibiae. The dotted lines give the linear fit under the assumption of pure proportionality between muscle length and femur length. The solid line indicates the 1:1 proportion, for comparison. Values are shown in Table 1. (B) Relation between middle leg muscle resting length and fibre length. Extensor muscle fibre length depends on resting muscle length (P<0.03), whereas no such dependence is present in the flexor muscle (dotted regression lines). Open, fibres in muscle medial regions; filled circles, fibres in muscle proximal regions. Values are means ± s.d. from nine experiments.

	Results

TOP Summary Introduction Materials and methods Results Discussion References

As for all orthopteran legs, both the extensor and flexor tibae of the stick insect are pinnate muscles (Storrer, 1976). This fibre arrangement markedly enhances effective muscle cross-sectional area and maximum muscle force (Hildebrand, 1988) but decreases effective muscle length and maximum contraction velocity. Given a mean fibre diameter of 0.125 mm (Bässler and Storrer, 1980), and a mean of 156 fibres per muscle (N=4; min. n=146, max. n=172), an estimated mean cross-sectional area of the extensor tibiae muscle is 1.9 mm². For pinnate muscles muscle length changes do not lead to the same changes in muscle fibre length; instead, muscle fibre length varies with muscle length times the cosine of the pinnation angle, which in turn varies as muscle length changes. For the extensor muscle, we calculated a range of angles from 8.2–12° for the proximal and 10.2–15.6° for the medial fibres within physiological muscle lengths. We ignored the cosine correction since it was only 3.7% for the largest angle (i.e. we treated the muscle fibres as though they were arranged parallel to the muscle longitudinal axis).

View larger version (14K): [in this window] [in a new window]   Fig. 4. Normalised muscle length change as a function of femur–tibia (FT)-joint angle. Extensor data are shown compared with +cos(), flexor data are compared with –cos() (solid lines).

View larger version (6K): [in this window] [in a new window]   Fig. 5. Relation between femur–tibia joint moment arm and femur length. Note that moment arm does not depend on femur length. Closed symbols, data from extensor muscles from two front legs (F; squares), two hind legs (H; triangles), and three middle legs (M; circles); open symbols, data for flexor tibiae muscles from five middle legs (circles) and two front legs (triangles).

Isometric muscle forces
Passive muscle forces
The data presented up to now show how muscle length and muscle fibre length depend on FT-joint angle. We next turned to measuring passive and active muscle tension.

The resting (passive) tension of the middle leg extensor tibiae muscle was measured by lengthening the muscle with linear ramps (described in the Materials and methods), beginning at its most relaxed state. The resting tension of a stretched muscle exhibits a phase–tonic time course due to muscle visco-elastic properties (Fig. 6A) (Bässler, 1983; Malamud, 1988). We defined the amplitude of the phasic component of the response as the peak force induced by the stretch and, somewhat simplified (see Discussion), the tonic tension as the force at which the rate of change of force was smaller than 0.3 mN min^–1. Both phasic- and tonic-component amplitude increased with increasing muscle stretch. Fig. 6B shows tonic force versus muscle length for nine extensor tibiae muscles and shows that passive tension increases from 0 mN at –0.75 mm to about 30 mN at +0.75 mm (0 mm=resting length). When leaving the innervation of the extensor tibiae muscle via nerve nl3 intact, passive tension was 2 mN higher as compared to the denervated situation (see Fig. 6B) throughout the muscle's working range (FT-joint angle between 30° and 180°, data not shown). An analysis of extensor muscle relaxation dynamics will be published elsewhere.

View larger version (13K): [in this window] [in a new window]   Fig. 6. Passive forces of the unstimulated, denervated middle leg extensor tibiae. (A) Time course of passive force induced by a ramp stretch of 0.5 mm amplitude in 1 s. Note the dynamic nature of the force change and its subsequent relaxation to a nearly constant value with maintained stretch. (B) Tonic resting force versus muscle length. Data from N=10 experiments (see text for details). Please note that the number for individual means can differ (6N10).

View larger version (10K): [in this window] [in a new window]   Fig. 7. (A) Average time course of a single twitch of the extensor tibiae (n=43). (B) Force development generated upon repetitive stimulation with differing stimulation frequencies between 30 Hz and 200 Hz (see values in Table 2). Broken line, resting force level.

At 30–50 Hz stimulation frequency, twitch contractions merged into an incomplete tetanus in all preparations (N=38; Table 2). The summated `steady-state' force amplitude of the resulting contractions varied between preparations (see Table 3 and Fig. 8). At 80 Hz and higher stimulation frequencies, tetanus was complete in all preparations. At higher stimulation frequencies the force level only increased slightly with increasing stimulation frequency (Fig. 7B, Table 2), and was maximum with no further increase at 200 Hz. Maximum force values from nine preparations are given in Table 3.

View larger version (35K): [in this window] [in a new window]   Fig. 8. Isometric force (mN, left ordinate) as a function of muscle stretch (mm, lower abscissa) and relative fibre length (top abscissa) in the middle leg extensor tibiae of nine individuals. Values for (A) single twitch (N=8), (B) 50 Hz (N=9) and (C) 200 Hz (N=9), see maximum values in Table 3. For comparison with other muscles, the right ordinate shows stress (N cm–2). The working range from 180° (fully stretched joint angle) to 30° (maximally flexed joint angle) is marked. Please note that the number for individual means can differ; for single twitch measurements 3N9, for 50 Hz and 200 Hz measurements 5N9.

Length dependence of contraction force
The isometric force generated by a muscle depends on the motoneuron firing rate and muscle length, because overlap of the filaments changes as muscle length changes (Rack and Westbury, 1969; Brown et al., 1999). We analysed the dependence of extensor tibiae muscle force upon these parameters by stimulating the muscle with single stimuli, intermediate frequency (10–80 Hz) and high frequency (200 Hz) tonic stimulations.

Fig. 8 shows (middle leg) extensor tibiae single twitch force, isometric force at 50 Hz stimulation frequency and at 200 Hz versus muscle fibre length. Besides the expected general length dependence on filament overlap (Gordon et al., 1966), the dependence of force generation on fibre length and stimulation frequency varied between experiments. In the single twitch measurements (Fig. 8A), maximum force varied eightfold (1.7–13 mN/0.09–0.68 N cm^–2); for the 50 Hz stimulation (Fig. 8B), fivefold (23.4–114.9 mN/1.22–6.00 N cm^–2); and for the 200 Hz stimulation (Fig. 8C), 2.7-fold (61.8–165.6 mN/3.23–8.65 N cm^–2) (see also Table 3).

Length dependencies at 200 Hz show maximum muscle force in the upper third of the working range, but curves with lower stimulation frequencies show a tendency to shift maximum force up to more elongated fibres (Fig. 8). This shift occurs in other muscles (Rack and Westbury, 1969; Brown et al., 1999), but it could vary between experiments; two examples with maximum and minimal frequency dependence of the length at which the muscle develops maximum force are shown in Fig. 9A,B.

View larger version (23K): [in this window] [in a new window]   Fig. 9. Two extreme examples of extensor tibiae muscle single twitch and tetanical forces as a function of muscle stretch. With increasing stimulation frequency, the maximum force in A moves markedly towards lower values of fibre length while the shift in B is much less prominent. Axes as in Fig. 8.

View larger version (9K): [in this window] [in a new window]   Fig. 10. The loaded release experiment. (A) Response of a tetanically activated muscle to a switch from isometric to isotonic conditions with counterforces less than the force the muscle had developed during its isometric contraction (see different `Force' traces). (B) The muscle shows an abrupt length change (bracket `a') followed by a smooth, initially linear contraction (b; see text for details).

(1)

where V is velocity, P is force, P₀ is the isometric force just prior to the switch to isotonic conditions, and `b' and `a' the asymptotes (lying at negative values of P and V).

View larger version (15K): [in this window] [in a new window]   Fig. 11. (A) Variability of the force–velocity characteristic (Hill hyperbola, see Materials and methods for details) of the extensor muscle (see values in Table 4). Different symbols represent different animals. (B) Deviation of the force–velocity curve from hyperbolic shape in the region of P/P0 between 0.6 and 1 for one typical experiment (open circles). In the region P/P0>1, force–velocity measurements are shown under stretch (filled circles, data from six experiments) to demonstrate the sigmoid zero crossing of the Hill hyperbola. The force axis was normalised to the maximum isometric contraction force P0. See text for details.

(2)

This normalised version of the Hill-hyperbola gives a common value for the asymptotes c=b/V₀=a/P₀; the reciprocal value 1/c can serve as a measure for the curvature of the Hill hyperbola (Josephson, 1993). Table 4 gives the values V₀, P₀, c and R² for the measurements of Fig. 11A.

Deviations of the Hill curve from the hyperbolic shape
As reported by Edman et al. (Edman et al., 1976; Edman, 1988), we also find deviations of the Hill curve from the hyperbolic shape in the region of high forces (P near P₀) at low contraction velocities. The measured contraction velocities in this region lie systematically below the hyperbolic fit function. To demonstrate this, we took one example of the experiments in Fig. 11A and normalised the force values to P₀ (open circles in Fig. 11B). If we fit a hyperbola as before, but take only data below 0.65P₀, the resulting curve (solid line) shows the deviation of the measured velocities in the range >0.65P₀ from the hyperbolic fit.

This deviation persisted for values of P/P₀>1 (negative contraction velocities, i.e. stretch) (Edman et al., 1976; Edman, 1988). We therefore conducted a set of six similar experiments, changing from tetanical isometric contraction to loads higher than P₀, resulting in a `loaded stretch', where the muscle is stretched with a certain force and stretch velocity is measured (Fig. 11B, filled circles to the right of P/P₀=1).

As maximum contraction velocity V₀ in all experiments did not scatter very much we found it reasonable to put the collected original velocity data as a function of P/P₀ in one diagram. The composite picture demonstrates that the Hill curve does not cross the abscissa with a non-zero slope (as the hyperbola would), but that the zero crossing is instead sigmoidal with an almost horizontal tangent.

Length dependence of V₀
Measuring the Hill curves at different muscle lengths and under minimal load reveals how maximum contraction velocity V₀ depends on muscle length. To investigate this issue we performed two experiments (Fig. 12), one experiment was done at 200 Hz and at 50 Hz stimulation frequency and with a rest load of 0.7 mN compared to a maximum isometric force of 69(35) mN at 200(50) Hz (filled symbols). The other experiment was done only at 200 Hz (open symbols) at a rest load of 4.5 mN compared with the maximum isometric force of 140 mN. Similar to the force–length characteristics (Fig. 8), the velocity–length curves showed a monotonic and nearly linear increase within the muscle's working range, followed by a plateau at longer fibre lengths.

View larger version (10K): [in this window] [in a new window]   Fig. 12. Length dependence of maximum contraction velocity V0 in two experiments. One experiment (filled symbols) was done at 200 Hz and at 50 Hz stimulation frequency. The other experiment (open symbols) was done only at 200 Hz (see text for details).

View larger version (8K): [in this window] [in a new window]   Fig. 13. Hill curves at different stimulation frequencies from one animal (see values in Table 5).

View larger version (19K): [in this window] [in a new window]   Fig. 14. Dependence of the parameters of the Hill curve on stimulation frequency. (A) Dependence of the isometric force P0 on stimulation frequency (N=11). (B) Dependence of maximum contraction velocity V0 on stimulation frequency (N=5). In both A and B the data from each experiment were fitted by an exponential saturation function leading to the same mean `frequency-constant' of about 60 Hz for both dependencies. The dotted lines show the shallowest and steepest fits and the bold lines show mean ± s.d. Please note that the number for individual means can differ (in A, 6N9; in B, 4N5). (C) Coupling of maximum contraction velocity V0 with isometric force P0 under variation of stimulation frequencies (N=5). The values of each set of experiments are marked by connecting straight lines. For all individual experiments the data showed a significant correlation (P<0.03). Note that the maximum values V0 and P0 at 200 Hz for each experiment roughly determine the overall slopes of the diagrams. See text for details.

Plotting V₀ versus frequency can be fitted with a similar function and shows a nearly identical dependence with an even higher mean R²=0.93±0.06 (min 0.86, max 0.98) and a similar mean f₀=57.7±17.9 Hz (Fig. 14B, N=5). Again the steepest and the shallowest fit as well as the mean ± s.d. are shown.

The paired data of V₀ and P₀ in the individual experiments in all cases showed a significant positive correlation (P<0.03, Fig. 14C, N=5). The V₀/P₀ ratios of the maximum values of V₀ and P₀ at 200 Hz correlate well with the slopes of the regression lines (P<0.05).

No correlation was found between the frequency constants f₀ of the corresponding P₀ and V₀ fits of the individual experiments, nor between the f₀ values and the P₀ and the V₀ values.

At present, a similar analysis is not possible for the curvature parameter `c', mainly because we only gathered a limited number of values for `c' in our experiments, since we looked for pairs of V₀ and P₀ and not for the complete Hill curve in most cases. However, the limited data do not show a systematic change of `c' with stimulation frequency; the mean of the 14 `c' values we evaluated was 0.5±0.22.

Series elasticity and quadratic spring property of the activated extensor tibiae muscle
To analyse the series elasticity (spring constant) of the extensor tibiae muscle, we plotted the isotonic force change in the loaded release experiment (steps in the force trace in Fig. 10A) versus the associated fast initial length change (bracket at `a' in Fig. 10B). In order to obtain a force–length characteristic comparable with similar measurements in the literature (Jewell and Wilkie, 1958) and with technical spring characteristics the data are presented in the following way (Fig. 15A): we plotted the independent, controlled force at the ordinate, the resultant length change on the abscissa. The origin length (0) for every set of data was the initial length of the isometric contracted muscle, and muscle shortening is thus plotted as negative values. The starting value of the force in each experiment is the tetanical isometric force at the given stimulation frequency. Each force step is plotted from this starting point, leading to a sequence of diminishing forces versus the ever more shortening length of the releasing muscle. Fig. 15A shows that the series elasticity of the extensor tibiae muscle is nonlinear and stiffness increases with increasing stretch. Consideration of these curves shows that the curve generated by 30 Hz stimulation, with its maximum force of 60 mN, has the same shape as the 60 mN portions of the curves generated by the 50, 80 and 200 Hz stimuli. These data can therefore be combined by translocating the 30, 50 and 80 Hz curves to the left so they overlay the 200 Hz, maximum force curve [see also Jewell and Wilkie (Jewell and Wilkie, 1958) for further discussion of this issue]. When this is performed (Fig. 15B) all the force–length curves can be fitted by a single parabola P(x)=ß*x² with ß being the quadratic spring constant (for this experiment, ß was 13900 mN mm^–2; maximum stiffness reached 3150 mN mm^–1). This procedure can be done with all series elasticity measurements at 200 Hz and led to an average ß of 15000±2600 mN mm^–2; mean R²=0.997±0.002 (N=5). It is currently unclear if the spring constant varies with muscle length.

View larger version (10K): [in this window] [in a new window]   Fig. 15. Nonlinear spring-characteristic of the series elasticity. (A) Fast length change in the loaded release experiment starting from different tetanical force levels; change in force over the resultant change in length. (B) Common parabola fit for all shifted curves in A. The respective starting values are marked by circles. See text for details.

	Discussion

TOP Summary Introduction Materials and methods Results Discussion References

 
Understanding the operation of a neuromuscular system requires investigating the system at both the neural and the muscular levels (e.g. Blickhan et al., 2003; Hooper and Weaver, 2000). We have investigated the extensor tibiae of the stick insect middle leg in great detail, because (i) the extensor tibiae muscle is a fairly simply innervated leg muscle with only two excitatory and one inhibitory motoneurons (Bässler and Storrer, 1980), (ii) its pattern of muscle fibre innervation and muscle fibre composition has been studied thoroughly (Bässler et al., 1996), and (iii) the activation of the motoneurons during active leg movements is known to some extent (Bässler and Büschges, 1998; Büschges, 2005). For example, when walking on a double treadwheel in late stance and during the initial part of the swing phase, as well as during the execution of compensatory leg reflexes, excitatory and inhibitory extensor motoneurons are activated together with instantaneous frequencies of up to 200 Hz (Büschges et al., 1994; Kittmann et al., 1996; Schmitz and Hassfeld, 1989). We therefore decided to stimulate the nerve with sufficient strength to activate all motoneurons in order to generate motoneuron activity resembling to some extent at least what is present in the swing phase of the stepping cycle. Although separate stimulation of the individual excitatory motoneurons is not possible due to FETi, SETi and CI₁ sending their axons through the same leg nerve (Bässler and Storrer, 1980), the mode of motoneuron activation, present in walking, can be mimicked by electrical stimulation of all three motoneurons at the same time. The extensor–tibiae system thus provides the opportunity to examine motoneuron output and muscular properties in a functional context. We have presented here a fairly complete picture of the geometry of the femur–tibia joint and the biomechanical properties of the extensor tibiae muscle. We will now discuss these data in terms of insect leg joint function and compare this system to previously analysed muscle–joint systems.

Geometry of the femur-tibia joint and tibial muscles
We found that the length of the flexor and extensor muscles of all three legs may be regarded as proportional to femur length. The tibial muscles of the middle legs and hind legs are 90–95% of the femur length, while for the front legs the proportion is only 73% (Fig. 3, Table 1). The latter is due to a specialized very narrow femur base in the front legs, that allows them to be held rostrally over the animal's head with virtually no angle of deviation from the animal's longitudinal axis, and thus to mimic twigs (Marquardt, 1940; Bässler, 1983). Due to this specialization there are no tibial muscle fibres present in the proximal 25% of front leg femurs.

Both femoral muscles are typical pinnate insect muscles. Due to this arrangement of the muscle fibres, their lengths are much shorter than muscle lengths (Fig. 3B). This almost parallel arrangement of many short muscle fibres, with one end attached at the inner cuticle and the other at the muscle tendon, markedly enhances the effective muscle cross-sectional area and the ability of the muscle (relative to muscle length) to generate force within a rigid, confined space (Full, 1997). For instance, the extensor tibiae muscle has an effective cross-sectional area of about 1.9 mm², while the inner lumen of the femur has a total cross-sectional area of only about 0.52 mm². This issue is even more important because the femur contains two tibial muscles, both the extensor and the flexor tibiae.

The relation between extensor and flexor tibial muscle length and the actual FT-joint angle (), follows approximately a cosine function (Fig. 4). Importantly, cosine fit of muscle length over joint angle (Fig. 4) for the flexor muscle must not be exact. If it were, the torque of the flexor would become zero at the fully extended position of =180°, and the muscle therefore would be unable to move the tibia out of this position. Therefore, we must assume that the fitted cosine function does not properly describe the range of angles close to the endpoint at 180°. Close inspection of the curves presented in Fig. 4 reveals that the muscle length does not reach the extended position at =180° with a horizontal tangent, but rather that in this range of FT-joint angle the values still show a slope different from zero. It appears that the joint mechanics thereby avoid the problems of a mechanical dead point for the tibial muscles at extended joint angles. This also means that more detailed fits must be used to model joint geometry at maximum joint angles.

The flexor's moment arm being twice as long as the extensor's seems to be an adaptation to the flexor's longer fibres (see Fig. 3B) compared to those of the extensor (Fig. 5), leading to a similar relative length change of the muscle fibres during walking. This difference also has functional consequences: since the middle leg's flexor is mainly active during stance, a longer flexor moment arm would be an advantage for torque generation, whereas the extensor benefits from its shorter moment arm to manage fast movements during swing phase (Bässler and Büschges, 1998; Graham, 1985).

Muscle forces
Passive muscle force
All muscles exhibit dynamic forces against passive lengthening. The muscle's response to a step within the stretch range has a power function characteristic (Thorson and Biedermann-Thorson, 1974). This means that stress relaxation becomes slower but will never reach a constant equilibrium (Malamud, 1989). This relaxation dynamic will be published elsewhere, so only a simplified description is provided here and we use the value of passive force after the rate of change, which was less than 0.3 mN min^–1, as the static component. We show in the present work that the static portion of resting force in the extensor tibiae muscle increases almost exponentially with muscle length (Fig. 6). This systematic change in static resting tension is in good agreement with other muscles, e.g. the tergo-coxal flight muscle in Schistocerca (Malamud, 1989), the various extensors of the cockroach Blaberus discoidalis hind leg (Full and Ahn, 1995), as well as the feline hind limb muscle (Brown et al., 1999). However, comparison of the data presented here with those of other insect muscles shows that the resting tension of the extensor tibiae muscle is much weaker than in other systems. The maximum static resting force of the extensor tibiae is only about 20% of the active muscle force (maximum static resting tensions at 140% muscle length stay below 20 mN, whereas the muscle's maxium tetanical contraction force was as large as 150–200 mN), much lower than the values published (Malamud, 1989; Brown et al., 1999). The relative weakness of extensor muscle static resting force is even more striking when considered within the muscle's working range, where it does not exceed 5 mN. Of course, during walking it is not the static but the dynamic resistance to stretch that is likely of greatest functional importance. Fig. 5A shows that the dynamic force component can be as much as twice the static component. None the less, our data suggest that passive resistance to stretch forces is likely to be much smaller than active forces (i.e. from the muscle's antagonist) at the femur–tibia joint.

Time course of single twitch and maximum tetanical forces
Single stimulation of the three motoneurons that innervate the extensor tibiae muscle generates a single twitch contraction, whose amplitude and time course can be compared with those in other insect muscles. First, the maximum force of the single twitch is only a few percent (ca 1%–10%) of the muscle's maximum force output (Table 3). This clearly differs from Schistocerca flight muscle Tcx₂ [in which single twitch amplitude is about 60% of maximum force (Malamud and Josephson, 1991)] and cockroach leg extensor muscle 177c (27%), whereas another leg extensor muscle 179 of the same animal (both 2%) lies within the range of our results (Ahn and Full, 2002).

The time course of the single twitch contraction (Fig. 7A) is much slower than that of the flight muscle from Schistocerca (Malamud and Josephson, 1991) and the flight muscles of Manduca (Stephenson and Josephson, 1990). The values for the extensors of the cockroach (Ahn and Full, 2002) and the front and middle leg extensors of the locust (Burns and Usherwood, 1978) are very similar to ours. Although the maximum tetanical force of the extensor tibiae is unexpectedly high (ranging from 80–200 mN), the normalized stress values (4.2–10.5 N cm^–2) are in the lower range of other insect muscles, much lower than e.g. the values for Tcx₂ (36 N cm^–2) of the locust or the leg muscles of the cockroach (25 and 47 N cm^–2) (Ahn and Full, 2002). Tettigoniid wing muscles [4–15 N cm^–2 (Josephson, 1993)] lie within the range of our results.

Comparing our values with the middle and front leg extensors of the locust (Burns and Usherwood, 1978) reveals, surprisingly, that the absolute maximum force values of this `normal' walking leg at 30–35 mN are much below the values for the stick insect. Unfortunately, we could not find published values for the effective cross-sectional areas of these muscles, so we cannot compare the values for maximum stress. Taken together, our results on the stick insect extensor muscle appear to be in good agreement with the expected requirements of a leg muscle in a slow walking and climbing insect: muscle forces reaching relatively high maximum values can be controlled by a broad range of motoneuron frequencies. The variability of muscle forces between different animals is discussed below.

Force–length characteristics
Due to the short length of pinnate muscle fibres, large length changes occur during contraction and relaxation. Therefore, it is particularly important to know the influence of fibre length on force in these muscles (cf. Gordon et al., 1966). Our data show that the extensor tibiae muscle shows a typical ascending-plateau-descending force versus muscle length relationship, that the muscle's working range lies on the ascending limb of this curve, and that considerable animal-to-animal variability exists with respect to the length and activation dependence of force production (Fig. 9, Table 3). A striking observation was the large change in force produced throughout the muscle working range, which sometimes showed a fourfold change in force as FT-joint angle changed from 30 to 180°, but generally resulted in a twofold change. An additional important point is that the muscle's maximum contraction force was normally generated in the upper third of its working range, thus the muscle was not generally exposed to any possible instabilities resulting from working with a non-monotonic force–length characteristic.

The dependence of maximum contraction force on fibre length at different activation levels is known from other muscles, for example feline SOL (Rack and Westbury, 1969), amphibian muscles (Stephenson and Wendt, 1984), and feline CF muscle (Brown et al., 1999). Two possible mechanisms have been proposed for this fact: length-dependent Ca²⁺-release and length-dependent change of cross bridge attachment/detachment rate (Brown and Loeb, 1999). For the stick insect extensor tibiae, this particular muscle property is variable and does not correlate with the variation in maximum contraction force. Comparable insect muscle data are not known to us.

Dynamics of muscle contraction: loaded release and the Hill hyperbola
To evaluate the dependence of isotonic contraction velocity on load (Hill curve), we conducted a series of loaded release experiments (Fig. 10). Immediately after the transition from isometric to isotonic recording in these experiments contraction oscillations occurred (Fig. 10B). Such oscillations have been reported previously (Jewell and Wilkie, 1958; Edman, 1988; Edman and Curtin, 2001; Siebert et al., 2003) and may reflect the fast dynamic properties of the contractile mechanism under fast length change. After these oscillations we can measure the initial contraction velocity, which declines rather rapidly to zero (Fig. 10A). This is the result of the short fibre length and of the force–length characteristics of this muscle: by contracting, the muscle moves to a shorter length with reduced maximum force, which results in a reduced contraction velocity. The contraction stops when the maximum isometric force is equal to the load set in the experiment (see the different length traces in Fig. 10A). Furthermore contraction velocity is influenced directly by muscle length (see below).