
==== Front
Pflugers Arch
Pflugers Arch
Pflugers Archiv
0031-6768
1432-2013
Springer Berlin Heidelberg Berlin/Heidelberg

39043889
2991
10.1007/s00424-024-02991-4
Muscle Physiology
Impact of lengthening velocity on the generation of eccentric force by slow-twitch muscle fibers in long stretches
Weidner Sven sven.weidner@inspo.uni-stuttgart.de

1
Tomalka André 1
Rode Christian 2
Siebert Tobias 13
1 https://ror.org/04vnq7t77 grid.5719.a 0000 0004 1936 9713 Department of Motion and Exercise Science, University of Stuttgart, Allmandring 28, 70569 Stuttgart, Germany
2 https://ror.org/03zdwsf69 grid.10493.3f 0000 0001 2185 8338 Institute of Sport Science, Department of Biomechanics, University of Rostock, Rostock, Germany
3 https://ror.org/04vnq7t77 grid.5719.a 0000 0004 1936 9713 Stuttgart Center of Simulation Science, University of Stuttgart, Stuttgart, Germany
24 7 2024
24 7 2024
2024
476 10 15171527
26 3 2024
1 7 2024
15 7 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
After an initial increase, isovelocity elongation of a muscle fiber can lead to diminishing (referred to as Give in the literature) and subsequently increasing force. How the stretch velocity affects this behavior in slow-twitch fibers remains largely unexplored. Here, we stretched fully activated individual rat soleus muscle fibers from 0.85 to 1.3 optimal fiber length at stretch velocities of 0.01, 0.1, and 1 maximum shortening velocity, vmax, and compared the results with those of rat EDL fast-twitch fibers obtained in similar experimental conditions. In soleus muscle fibers, Give was 7%, 18%, and 44% of maximum isometric force for 0.01, 0.1, and 1 vmax, respectively. As in EDL fibers, the force increased nearly linearly in the second half of the stretch, although the number of crossbridges decreased, and its slope increased with stretch velocity. Our findings are consistent with the concept of a forceful detachment and subsequent crossbridge reattachment in the stretch’s first phase and a strong viscoelastic titin contribution to fiber force in the second phase of the stretch. Interestingly, we found interaction effects of stretch velocity and fiber type on force parameters in both stretch phases, hinting at fiber type-specific differences in crossbridge and titin contributions to eccentric force. Whether fiber type-specific combined XB and non-XB models can explain these effects or if they hint at some not fully understood properties of muscle contraction remains to be shown. These results may stimulate new optimization perspectives in sports training and provide a better understanding of structure–function relations of muscle proteins.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00424-024-02991-4.

Keywords

Skeletal muscle
Contractile behavior
Stretch
Give
Soleus
http://dx.doi.org/10.13039/501100001659 Deutsche Forschungsgemeinschaft SI841/17-1 SI841/17-1 RO5811/1-1 (project number 405834662) SI841/17-1 Weidner Sven Tomalka André Rode Christian Siebert Tobias Universität Stuttgart (1023)Open Access funding enabled and organized by Projekt DEAL.

issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

Eccentric muscle contractions are associated with unique features compared to isometric or concentric contractions, i.e., increased force, work, and performance at decreased oxygen consumption, reduced metabolic cost (ATP), and improved energy efficiency [1, 18, 31, 43, 47, 65]. Numerous studies have shown that muscle force rises steeply during the early phase of the stretch, followed by a relatively compliant transient phase. This behavior likely depends on the stretch velocity [54, 63, 85]. More specifically, the initial steep force rise (slope1, Fig. 1) often ends in a characteristic force peak (s2), and both increase with stretch velocity [79, 85]. After this initial force increase, muscle [19, 27, 75] and fiber [10, 24, 79] forces decrease in fast stretches, called Give [19]. In long stretches, muscle [75] and fiber [76, 79] forces increase nearly linearly again after the Give phase. The corresponding force slope (slope2) is larger than the force slope of the underlying total (active + passive) isometric force–length relation [77]. Slope2 also increases with stretch velocity [79, 85]. Ample evidence supports the idea of a cumulative mechanism that combines crossbridge (XB) and non-crossbridge (non-XB) structures (e.g., titin) to the force response during active muscle lengthening [76]. While slope1, s2, and Give are primarily influenced by XB behavior [22, 54, 62, 79], slope2 depends more on non-XB structures [77, 85]. Thus, the velocity dependence of the mentioned parameters likely stems from XB and non-XB structures.Fig. 1 Exemplary force–length trace obtained in a long isovelocity stretch. Fiber forces are normalized to maximum isometric force, F0, and fiber length is normalized to optimal fiber length, lopt. Experiments start with an isometric phase until force achieves a plateau (not shown). Here, the fully activated skinned soleus fiber is then stretched from 0.85 to 1.3 lopt at 1 maximum fiber contraction velocity. s2 is the first local force maximum during the stretch. sg is the local force minimum during the stretch, and Give is the difference between s2 and sg. The force at the end of the stretch is se. slope1 and slope2 are the force slopes resulting from linear regression analysis of the force–length data between the initial isometric force and s2 and the second half of the stretch, respectively. ls2 and lsg are the lengths, where s2 and sg occurred. A raw data set for three experiments with different stretch velocities is shown in Figure S1

Specificities of these structures influence the velocity-dependent muscle force generation. Slow-twitch muscle fibers show lower myosin ATPase activity, lower XB-cycling frequencies, and lower maximal shortening velocities compared to fast-twitch muscle fibers [5, 8, 64]. Furthermore, different muscles differ in non-XB structures, e.g., collagenous structures (Endo-, Peri-, Epimysium [84]) or titin isoforms [23, 55, 83], resulting in different passive properties [70] impacting eccentric force generation. Recently, Weidner et al. [85] investigated the force produced by fast-twitch extensor digitorum longus (EDL) muscle fibers during long stretches over an extensive range of stretch velocities given as percentages of maximum contraction velocity, vmax. slope1 and slope2 increased by 195% and 377% from 0.01 to 1 vmax stretch velocity, respectively. Furthermore, Give, defined as the force decrease from the initial peak to the following local minimum in force, was absent for the slow stretch velocity (0.01 vmax), appeared at 0.1 vmax, and increased at vmax [85]. However, it is unclear if slow-twitch muscle fibers show a similar velocity dependency of stretch-related parameters within the same experimental conditions as reported for fast-twitch fibers in Weidner et al. [85].

Thus, the study aims are (1) to examine the velocity dependence of eccentric force parameters in soleus fibers, a typical slow-twitch muscle [12], and (2) to identify fiber type-specific differences by comparison with a fast-twitch muscle. The two muscles were selected for analysis because both muscles (soleus and EDL) predominately express different MHC isoforms. The soleus muscle expresses the slow type 1 MHC isoform (96.1 ± 2.9%). In contrast, the EDL expresses predominately the fast type 2A (18.8 ± 1.7%) and type 2B (75.7 ± 2.2%) MHC isoforms [4, 71]. We performed eccentric muscle fiber stretches over the same length range (0.85 to 1.3 optimal fiber length, lopt) with equal relative stretch velocities (0.01, 0.1, and 1 vmax) as Weidner et al. [85] with the fast EDL muscle fibers. While 0.01 and 0.1 vmax correspond to relative lengthening velocities respectively in slow- and fast-twitch muscles during walking [3], 1 vmax was chosen to cover sprint velocities. We normalized the stretch velocities to vmax to account for fiber type-specific differences in contraction velocity (EDL’s vmax is about fivefold higher compared to soleus). We expect fiber type-independent force characteristics in the XB-dominated short range of the stretch. Because of lower absolute velocities, we expect lower forces in slow-twitch fibers in the long range of the stretch due to the likely viscoelastic non-XB (titin) contribution to the total force [6, 11, 28, 29, 79].

Materials and methods

Animal and tissue preparation

Soleus muscles were extracted from seven male Wistar rats (age, 8 to 10 months; weight, 300–370 g; 12-h light to 12-h dark cycle; housing temperature, 22 °C) immediately after euthanization. The use of all animals for this study has been approved according to the German animal protection law (Tierschutzgesetz, §4 (3); Permit Number: T 201_21 ST).

The muscles were taken from the left hind limb. Permeabilized single muscle fiber preparation, storage, and activation followed [77, 78]. Fiber bundles were dissected, permeabilized in skinning solution at 4 °C for 30 min, and then stored in a 50% glycerol and 50% skinning solution mixture at − 20 °C for 6 to 8 weeks. On experiment day, single fibers were dissected from the muscle bundles with fine forceps under a dissecting microscope and cut to 1.2 mm in length. They were then clamped on both sides with aluminum foil T-shaped clips. The fibers were treated with a relaxing solution containing 1% (v/v) Triton-X 100 for 2 to 3 min at 4 °C to remove internal membranes [42].

Experimental setup

The first step involved transferring the muscle fibers from the skinning solution to the experimental chamber of the fiber-test apparatus (600A, Aurora Scientific, Canada). Then, the fiber-clip unit was attached to a Model 403A force transducer (Aurora Scientific, Ontario, Canada) and a Model 322 C-I length controller (Aurora Scientific, Ontario, Canada). The setup was mounted on the x–y stage of an inverted Eclipse Ti-S microscope (Nikon, Japan). The fiber ends were fixed with glutaraldehyde in rigor solution, and the T-clips were secured to the apparatus with fingernail polish to enhance stability and improve mechanical performance during the experiment [30].

The sarcomere length was measured in the middle segment of the fibers. The passive fiber length was adjusted to achieve an optimal sarcomere length (ls0) of 2.5 ± 0.05 µm (mean ± standard deviation) for maximal isometric force (F0) development [72]. The fiber width (w) and height (h) were measured at approximately 0.1-mm intervals along the entire length using a 10 × ELWD dry-objective (NA 0.60, Nikon, Japan) and a 10 × eyepiece. The fiber cross-sectional area was calculated assuming an elliptical shape (pi.h.w/4). A high-speed video system (Aurora Scientific 901B, Canada) combined with a 20 × ELWD dry-objective (NA 0.40, Nikon, Japan) and an accessory lens (2.5 × , Nikon, Japan) was used to visualize the striation pattern and track sarcomere length changes dynamically and accurately.

Solutions

The relaxing solution contained (in mM) 100 TES, 7.7 MgCl2, 5.44 Na2ATP, 25 EGTA, 19.11 Na2CP, and 10 GLH (pCa 9.0); the pre-activating solution 100 TES, 6.93 MgCl2, 5.45 Na2ATP, 0.1 EGTA, 19.49 Na2CP, 10 GLH, and 24.9 HDTA; the activating solution 100 TES, 6.76 MgCl2, 5.46 Na2ATP, 19.49 Na2CP, 10 GLH, and 25 CaEGTA (pCa 4.5); and the skinning solution 170 potassium propionate, 2.5 MgCl2, 2.5 Na2ATP, 5 EGTA, 10 IMID, and 0.2 PMSF. The storage solution is the same as the skinning solution, except for the presence of 10 mM GLH and 50% glycerol (v/v). Cysteine and cysteine/serine protease inhibitors (trans-epoxysuccinyl-L-leucylamido-(4-guanidino) butane, E-64, 10 mM; leupeptin, 20 µg ml−1) were added to all solutions to preserve lattice proteins and thus sarcomere homogeneity [42, 77]. pH (adjusted with KOH) was 7.1 at 12 °C. 450 U ml−1 of creatine kinase (CK) was added to all solutions except skinning and storage solutions. CK was obtained from Roche (Mannheim, Germany); all other chemicals were from Sigma (St. Louis, MO, USA).

Experimental protocol

All trials of the permeabilized fiber experiments were performed at a constant temperature of 12 ± 0.1 °C. At this temperature, the fibers were highly stable and endured active lengthening protocols and prolonged activations [56, 57]. The fibers were activated through calcium diffusion in the presence of ATP by immersing them subsequently in three different solutions: (1) a pre-activating solution for equilibration for 60 s; (2) an activation solution (pCa = 4.5) which resulted in a fast increase in force until a plateau was reached (increase in force less than 1% within 1.5 s) before stretching the fiber or keeping the fiber in isometric condition; (3) in a relaxing solution (pCa = 9.0) after the ramp or isometric experiment for 420 s. For the exact composition of the experimental solutions, see the “Solutions” section.

The experiments involved active eccentric ramps that stretched the fibers from a length of 0.85 to 1.3 lopt, a typical working range of muscles suitable for comparison with previous research [7, 85]. Single-skinned fibers were activated at 0.85 lopt (resulting in ∼2.0-µm sarcomere length under activation) and stretched to 1.3 lopt (resulting in ∼2.9 µm under activation) with constant stretch velocities of 0.01, 0.1, and 1 vmax (Fig. 1). vmax (0.47 ± 0.11 lopt s−1; n = 6 fibers) of skinned soleus muscle fibers was determined in a different set of muscle fibers according to [79] and agrees with literature data [16]. After each stretch, the same ramp was performed again passively. Isometric reference contractions were performed at lopt (2.5 ± 0.05 µm passive) before and after each ramp contraction to determine force degradation and F0. F0 was calculated as the average of isometric force before and after the stretches. Subsequently, the force responses during stretch were normalized to F0. In the eccentric contraction experiments, the isometric force decreased by an average of 1.55 ± 0.49% per activation (with a maximum of 10 activations per fiber). This rate of force loss is in line with other studies under similar conditions [14, 77] and demonstrates repeatable preparation routines and physiological fiber functionality. Additionally, the order of the experiments was randomized to eliminate any systematic effects of fatigue on stretch parameters. During the trials, the sarcomere length, force, temperature, and length controller position were recorded.

Data processing and statistics

An A/D interface (604A, Aurora Scientific, Canada) recorded force and length data at 1 kHz. Real-time software (600A, Aurora Scientific) was utilized for data acquisition. MATLAB R2018a (Mathworks, Natick, MA, USA) was used to analyze the collected data through a custom-written script. The sarcomere lengths were either reported in absolute values or divided by optimal sarcomere length. Forces were divided by individual F0, while fiber length l was divided by individual lopt. The statistical analysis was carried out using SPSS 29 (IBM Corp., Armonk, NY, USA) and MATLAB R2018a (Mathworks, Natick, MA, USA).

Analogous to Weidner et al. [85], this study focused on four prominent force values (s2, sg, Give, se) and two force slopes (slope1, slope2), as shown and defined in Fig. 1. We determined if these parameters, as well as the lengths ls2 and lsg (where s2 and sg occurred), varied with stretch velocity. One-way repeated measures ANOVAs explored the effect of stretch velocity on these parameters. In case of significant differences, Tukey’s HSD was used for post hoc analyses. Effect sizes for one-way ANOVA were classified as small (η2 < 0.06), medium (0.06 ≤ η2 ≤ 0.14), and large (η2 > 0.14), based on Cohen’s classification [13].

Two-way mixed ANOVAs explored the combined effects of stretch velocity and fiber type (our slow-twitch [soleus] fibers vs previously published fast-twitch [extensor digitorum longus, EDL] fibers [85]) on the force values and force slopes. In the case of a significant interaction or significant main effects, we report significant simple main effects of stretch velocity for each fiber type (one-way repeated measures ANOVA) and of fiber type for each stretch velocity (independent t-tests). Interactions are divided into the following categories according to [40]: Ordinal interactions are defined by a consistent rank order of treatment levels across all factor levels, with lines in interaction plots never crossing, and main effects are generally consistent and interpretable. Disordinal interactions exhibit a change in rank order, visualized by crossing lines on the plots, indicating non-monotonic relationships between factors, and main effects may be overshadowed by the interaction. Hybrid interactions combine elements of both, with rank order consistency in some factors and variability in others, resulting in mixed plots where lines may cross in some graphs but not in others, and main effects are partially interpretable depending on the factor investigated. To account for multiple comparisons, we used post hoc analyses. Results were expressed as mean ± standard deviation.

Results

The soleus fibers (n = 27) generated a maximum isometric tension of 98.6 ± 18.6 kN m−2 at lopt and a maximum contraction speed, vmax, of 0.47 lopt s−1 (EDL [85]: 2.42 lopt s−1). The fiber cross-sectional area was 0.006 ± 0.002 mm2. The total force during isovelocity stretches initially increased steeply and linearly; then, the force declined and increased again, showing a local minimum (Fig. 2a). All observed parameters (cf. Figure 1) increased significantly with stretch velocity (Table 1; Figure S3, S4, S5). The main effect of stretch velocity and subgroup comparisons were all significant with p < 0.001.Fig. 2 Force–length traces of fully activated slow- and fast-twitch fibers during long isovelocity stretches. Total fiber forces (active + passive force), F, are normalized to maximum isometric fiber force, F0; fiber lengths, l, are normalized to optimal fiber length, lopt. Ensemble averages (solid lines) and variances (shadowed areas) of active stretch forces and corresponding ensemble averages of passive stretch forces (dashed lines) are shown for stretch velocities of 0.01 (black), 0.1 (blue), and 1 vmax (red). Soleus data (slow-twitch, a) from this study and EDL data (fast-twitch, b) from [85] were obtained under similar experimental conditions. For orientation, figures include a schematic active isometric force–length relationship (gray dashed line). An enlarged section of the region of slope1 can be found in Figure S2

Table 1 Statistical analysis of the effect of stretch velocity on slow-twitch soleus fiber force parameters in isovelocity stretches

vmax		Mean (± standard deviation)	p-value|F-value, η2	Relative change (%)	
	0.01	0.1	1		0.01 to 0.1	0.1 to 1	0.01 to 1	
slope1	F0/ls0	9.34 (± 1.45)	16.57 (± 1.58)1	31.39 (± 3.44)1,2	 < 0.001|837, 0.97	77	89	236	
s2	F0	1.01 (± 0.03)	1.17 (± 0.04)1	1.51 (± 0.06)1,2	 < 0.001|2136, 0.988	16	29	50	
ls2	lopt	0.85 (± 0.00)	0.86 (± 0.00)1	0.87 (± 0.00)1,2	 < 0.001|2398, 0.989	1	1	2	
sg	F0	0.94 (± 0.03)	0.99 (± 0.05)1	1.07 (± 0.07)1,2	 < 0.001|89, 0.774	5	8	14	
lsg	lopt	0.87 (± 0.00)	0.9 (± 0.01)1	0.99 (± 0.02)1,2	 < 0.001|922, 0.973	3	10	14	
Give	F0	0.07 (± 0.02)	0.18 (± 0.04)1	0.44 (± 0.06)1,2	 < 0.001|1323, 0.981	157	144	529	
slope2	F0/ls0	0.34 (± 0.09)	0.93 (± 0.18)1	1.89 (± 0.36)1,2	 < 0.001|695, 0.963	174	103	456	
se	F0	1.27 (± 0.08)	1.57 (± 0.12)1	1.87 (± 0.19)1,2	 < 0.001|393, 0.938	24	19	47	
All one-way repeated measures ANOVAs and post hoc comparisons of subgroups were highly significant (p < 0.001) and effects were strong (η2 > 0.14). 1A significant difference compared to 0.01 vmax; 2a significant difference compared to 0.1 vmax. Maximum isometric force, F0; optimal sarcomere length, ls0; maximum contraction speed, vmax. The last three columns show relative changes in the dependent variable (left column) values from 0.01 to 0.1, 0.1 to 1, and 0.01 to 1 vmax. For dependent variable definitions, see Fig. 1

The slow-twitch (soleus) fiber force traces in this study showed similar qualitative behavior as the fast-twitch (EDL) fibers from a previous study [85] (Fig. 2). The force–length traces of single EDL and soleus muscle fibers during stretches with constant velocity did not reflect the changes in the slope of the underlying active isometric force–length relationship (Fig. 2b). This behavior was consistent across all stretch velocities tested (0.01, 0.1, and 1 vmax). Instead, force increased steeply and linearly, with a temporary decline observed for the two higher velocities. Beyond a local minimum (sg), force increased until the end of the stretch, exceeding initial forces by 80% (for 0.01 vmax) to 150% (for 1 vmax). No sg was observed for the lowest stretch velocity (0.01 vmax). However, unlike EDL, soleus fibers showed a force peak s2 and Give at 0.01 vmax stretch velocity. Table 2 summarizes the results of our statistical comparison of soleus and EDL parameters. Stretch velocity and fiber type interacted with four (slope1, slope2, s2, and sg) out of six variables (Fig. 3). Table 2 Statistical comparison of slow-twitch soleus fibers and fast-twitch EDL fiber data

	Fiber type	Stretch velocity	p-value|F-value	
η2	
0.01 vmax	0.1 vmax	1 vmax	Fiber type	Stretch velocity	Interaction	
slope1 (F0/ls0)	Soleus	9.34 (± 1.45)	16.57 (± 1.58)1	31.39 (± 3.44)1,2	 < 0.001|211	 < 0.001|457	 < 0.001|668	
EDL	7.62 (± 1.97)	9.69 (± 3.03)1	14.05 (± 2.62)1,2	0.847	0.923	0.762	
s2 (F0)	Soleus	1.01 (± 0.03)	1.17 (± 0.04)1	1.51 (± 0.06)1,2	0.33|0.96	 < 0.001|643	 < 0.001|16.2	
EDL	1.06 (± 0.11)	1.22 (± 0.08)1	1.41 (± 0.07)1,2	0.023	0.940	0.283	
sg (F0)	Soleus	0.94 (± 0.03)	0.99 (± 0.05)1	1.07 (± 0.07)1,2	0.74|0.12	0.008|7.76	0.003|9.74	
EDL	-	1.08 (± 0.12)	1.03 (± 0.12)	0.144	0.159	0.192	
se (F0)	Soleus	1.27 (± 0.08)	1.57 (± 0.12)1	1.87 (± 0.19)1,2	 < 0.001|89.0	 < 0.001|108	0.45|0.75	
EDL	1.87 (± 0.31)	2.22 (± 0.29)1	2.49 (± 0.47)1,2	0.685	0.724	0.018	
Give (F0)	Soleus	0.07 (± 0.02)	0.18 (± 0.04)1	0.44 (± 0.06)1,2	 < 0.001|22.8	 < 0.001|451	0.223|1.53	
EDL	-	0.14 (± 0.09)	0.38 (± 0.10)b	0.358	0.917	0.036	
slope2 (F0/ls0)	Soleus	0.34 (± 0.09)	0.93 (± 0.18)1	1.89 (± 0.36)1,2	 < 0.001|151	 < 0.001|325	 < 0.001|19.1	
EDL	0.72 (± 0.40)	2.03 (± 0.46)1	3.15 (± 0.76)1,2	0.800	0.895	0.335	
The last three columns show the p-values, F-values, and effect sizes (ES, partial eta-squared) of the two-way mixed ANOVAs pertaining to the main and interaction effects of fiber type and stretch velocity. 1A significant difference compared to 0.01 vmax; 2a significant difference compared to 0.1 vmax; values in bold, sig. different from soleus (p < 0.05) and sig. interactions. Maximum isometric force, F0. For dependent variable (left column) definitions, see Fig. 1

Fig. 3 Interaction stretch velocity × fiber type on force parameters. Subgroup means (points) and 95% confidence intervals (error bars) are shown. Stretch velocity and fiber type exhibit ordinal interaction on the initial force rise, slope1 (a); the force slope in the last half of the stretch, slope2 (d); hybrid interaction (only stretch velocity is interpretable) on the peak force s2 (b); and disordinal interaction (no main effect is interpretable) on the local force minimum sg (c). Stretch velocity has a main effect on slope1, s2, and slope2 (a, b, d right). fiber type has a main effect on slope1 and slope2 (a, d left). Maximum isometric fiber force, F0

Table 2 indicates that both the EDL and soleus exhibit similar behavior concerning the parameters slope1, s2, se, and slope2. Both muscles demonstrate a significant increase in these parameters with increasing stretch velocity. However, for sg and Give, there is no consistent behavior. The EDL only exhibits Give at 0.1 vmax and 1 vmax, while the soleus shows Give at all three tested velocities. Additionally, there are no differences in the level of sg for the EDL. In contrast, the soleus shows a significant increase in sg with increasing stretch velocity.

The significant interaction effects for slope1, s2, sg, and slope2 shown in Table 2 are visually depicted in Fig. 3. The force slopes exhibit an increase with stretch velocity. While slope1 is higher for the soleus, slope2 is higher for the EDL. For sg, the main effects cannot be interpreted due to the intersection in both graphs. Finally, s2 depends solely on the velocity of stretching and not on the fiber type.

Discussion

This research on slow-twitch muscle fibers of the rat soleus extends a recent study [85] that investigated the total force response of fast-twitch rat muscle fibers (EDL) under similar experimental conditions. Our investigations reveal three major similarities of the force produced by fully activated slow-twitch fibers compared with fast-twitch fibers [77, 85] during long isovelocity stretches (Fig. 2): (i) the forces first increase, then fall (except for EDL fibers at 0.01 vmax) and rise before or within the force–length relationship’s (FLR’s) plateau region depending on stretch velocity and increase in the range of the FLR’s descending limb; (ii) the force slopes in the range of the FLR’s descending limb increase with stretch velocity; (iii) all tested force and length parameters increase with stretch velocity (except for sg in EDL fibers).

Moreover, some differences between fiber types in short-range parameters persisted despite normalized contraction velocities that aimed at accounting for differences in absolute vmax.

Impact of stretch velocity in slow-twitch fiber experiments

Interestingly, for all stretch velocities tested, the force–length traces (colored lines) did not reflect slope changes of the underlying FLR and increased in the range of the FLR’s descending limb (Fig. 2a). In addition, force slopes in the range of the FLR’s descending limb increased with stretch velocity. These findings are in contrast with classic theories of muscle contraction [33–35] that would, e.g., predict a slope change in force when the number of XBs decreases (during the transition from the FLR’s plateau to its descending limb) and force slopes that decrease with stretch velocity in the range of the FLR’s descending limb. Commonly applied Hill-type muscle models [25, 67, 69, 81] approximating classic theories of muscle contraction represent neither the decrease of force in the FLR plateau during the stretch nor force slopes increasing with stretch velocity in the range of the FLR’s descending limb. This especially hampers simulations of movements involving fast muscle stretches induced by perturbations during locomotion [2] or large muscle stretches, e.g., accident predictions of multi-body models (e.g., OpenSim: [67], Anybody: [15, 58]).

XB contribution to eccentric force generation

We found significant increases in the initial force increase’s slope (slope1), its maximum (s2), and the length where this maximum occurred (ls2) with increased stretch velocity (Table 1, Fig. 3). Assume that the S2 region of the myosin molecule exhibits purely linear elastic behavior, and that the rate of detachment is constant for eccentric contraction as in Huxley’s classical crossbridge (XB) model (Huxley, 1957). Pulling a single XB with increasing speed, the detachment force s2 as well as ls2 would increase. However, different from our results, slope1 would remain constant. Stretching a cohort of XBs with increasing speed as we did when stretching fibers, a smaller fraction of XBs would be detached at a given length because less time passed when reaching this length, resulting in more XBs contributing in parallel to force generation. Thus, slope1 would increase in faster stretches. However, when comparing the slope within the first 20% of the slope1 period, there is no significant difference between the tested velocities (F(1.47,38.28) = 0.465, p = 0.57). This behavior suggests that a similar number of XBs are stretched for all three velocities during this period. Further experiments (with higher sample rates and more tested velocities) and modeling are required to explain the finer details of our results, like differences and similarities between fiber types.

Similar behavior has also been found in other studies [22, 54, 85] and has been explained with the viscoelasticity of XBs [9] or a non-XB component [49]. The observed length change associated with the development of s2 was less than 0.05 lopt, which is consistent with previous studies that stretched intact and skinned muscle fibers from various species [22, 45, 73, 79]. It is assumed that ls2 relates to the forcible detachment of XBs bound to actin [19, 45]. We would argue that it is within the capacity of the Huxley XB model to explain the increase of slope1, s2, and ls2 when considering the dynamic equilibration of the XB distribution from the isometric to the eccentric condition.

As in our study (Table 1, Figure S5c), long stretches of muscle fibers resulted in Give across different species and stretch velocities [19, 36, 38, 45, 85]. The detachment of XBs followed by a time-dependent restoration of a steady-state XB distribution may account for increased lsg with stretch velocity.

Non-XB contribution to eccentric force generation

Non-XB structures such as titin [37, 82] might complement the force recovery and dominate the force at longer elongation. Besides its fundamental role in organizing and maintaining sarcomeres, titin performs intricate and diverse functions in muscle contraction [21, 51, 76]. Titin interacts with many muscle proteins [44] and aids in generating force when the muscle is actively stretched. Titin’s force generation mainly relies on titin-actin binding during Ca2+ activation [6, 11, 17, 28, 46]. Several studies [75, 77, 85] report a quasilinear increase in force during the second half of a long stretch (slope2). Since XBs would contribute forces that decrease with contraction speed in the range of the FLR’s descending limb, the observed increase in slope2 (Table 1, Figure S3b) and the force se at the end of the stretch (Table 1, Figure S4c) are consistent with a strong, linear viscoelastic titin contribution to fiber force [11, 29].

Comparison between soleus (slow-twitch) and EDL (fast-twitch) muscle fiber kinetics during stretch

By normalizing to vmax (resulting in fivefold absolute stretch velocities for fast-twitch EDL fibers compared with slow-twitch soleus fibers), we aimed to balance the effects of fiber type. Though we found qualitatively similar behavior of slow- and fast-twitch fibers (Fig. 2), our experiments revealed some differences between the two fiber types (Table 2). The slow-twitch soleus force–length traces show comparably low inter-subject variance. This holds even when considering absolute stretch speeds. For example, absolute stretch speed 1 vmax in soleus is about twofold that of 0.1 vmax in EDL; regardless, the variance in the red data in Fig. 2a (soleus, 1 vmax) is lower than in the blue data in Fig. 2b (EDL, 0.1 vmax). The high variance in fast-twitch fiber forces has been documented in both individual muscle fibers and fiber bundles [55]. This points to a higher variability in passive structures and titin isoforms in fast-twitch muscles, especially because the contribution of non-XB structures to eccentric muscle force [77] and the variance (Fig. 2) increase with muscle length [59]. However, it should be noted that the EDL expresses two myosin isoforms [4, 71], which may contribute to the observed variability.

XB kinetics probably dominate the steep initial force increase and the force peak s2 and partially the local force minimum sg (and hence Give). This is supported by recent studies who investigated stretches in fiber bundles and single fibers with and without XB inhibitors [29, 63, 77]. Their results showed that, in the presence of crossbridge inhibitors, there was no steep initial increase in force, followed by a subsequent decrease. Accordingly, slope1, s2, and Give increase with stretch velocity (Table 2) for both fiber types. However, the factors fiber type and stretch velocity revealed ordinal interaction on slope1 and hybrid interaction on s2. slope1 increased stronger with stretch velocity for the soleus (Table 2, Fig. 3a) compared to EDL. Further, soleus’ s2, while similar for 0.01 vmax, was significantly lower for 0.1 vmax and higher than EDLs’ s2 for 1 vmax (Fig. 3b). Hence, normalization to vmax did not alleviate fiber type-related XB effects. Moreover, fiber type and stretch velocity revealed disordinal interaction on the local force minimum during the stretch, sg (Fig. 3c). This might be due to two nonlinear processes, the restoration of a steady-state XB distribution and viscoelastic non-XB dynamics, which contribute to force redevelopment. It remains to be shown whether such interactions can be explained by a fiber model accounting for XB and non-XB dynamics or if they hint at unknown features of muscle contraction.

Interestingly, forces do not drop significantly below the maximum isometric force (cf. sg in Table 2) at any stretching rate in both fiber types (Fig. 2). Possibly, titin-actin interactions [60] secure a certain force level to prevent damage. During activation, titin has been suggested to bind to actin [17, 39, 53, 74], thereby reducing the free titin spring length and increasing titin force in subsequent stretches [52, 60, 66]. The resulting titin force during fiber stretch could protect the muscle function when XBs tear off. With the attachment of new XB and the further increase in activation-dependent titin forces (that overcompensate the loss due to decreasing numbers of XB in the FLR’s descending limb range), the muscle can generate large forces and thus effectively avoid overstretching.

In this context, it is known that changes in titin-based stiffness are likely to play an important role in adjusting the passive and active properties of skeletal muscle in health and disease (for a detailed review on posttranslational modifications of titin, see [21]). It is also known that slow-twitch muscles (e.g., soleus [80]) usually express long titin isoforms accounting for low titin-based passive tension. In contrast, fast-twitch muscles such as the rabbit psoas muscle [20, 32, 55, 83] predominantly express shorter titin accounting for higher passive forces [50]. This is in line with higher values of slope2 in EDL than soleus fibers in quasistatic stretches at 0.01 vmax (Fig. 3d). However, further studies are needed to test the idea that the interaction of XB and non-XB structures represents a protective muscle adaptation to prevent destruction.

In the second half of the stretch, titin becomes even more important for force generation [77]. Soleus and EDL show increasing slope2 and force at the end of the stretch, se, with increasing relative stretch velocity (Table 2, Fig. 2). Hence, the viscoelastic non-XB contribution to fiber force outweighs the force slope-decreasing effect of XBs in the range where XB numbers decline [33–35] (the FLR’s descending limb). As expected due to the viscoelastic behavior of non-XBs, fivefold higher absolute stretch velocity leads to a larger slope2 and se in the EDL compared with the soleus stretch experiments (Table 2). It is also noteworthy that the interaction between crossbridges and non-crossbridge components (e.g., titin) leads to an approximately linear increase in force during the second half of the stretch (Fig. 2). However, all titin segments show a nonlinear (exponential) increase in force when stretched in isolation [41]. Further development of muscle models that take activation-dependent titin-actin interactions into account [26, 48, 52, 60, 66] is necessary. In particular, a fiber type-specific adaptation of these models could contribute to better predictions of musculoskeletal models [61, 68]. Incorporating an activation-dependent titin spring into a muscle model [60] already leads to linear titin force during stretch. How this translates to varying stretch speeds and the observed viscoelastic nature of non-XB contributions is not yet clear.

Conclusions

Both slow-twitch and fast-twitch fiber forces in long isovelocity stretches at different velocities are qualitatively similar. They increase sharply, decline (Give) except for EDL at 0.01 vmax, and recover with a positive slope in the length range where XB numbers decrease. Consistent differences in force parameters between fiber types and some interaction effects of stretch velocity and fiber type on these parameters highlight differences between slow- and fast-twitch fibers in the XB-dominated short range and the non-XB-dominated long range of the stretch. Whether a fiber type-specific combined XB and non-XB model can explain these effects or if they hint at some not fully understood properties of muscle contraction remains to be shown. Moreover, despite the well-established occurrence of Give in situ and in vitro, Give is yet to be established in vivo to understand its significance for regular muscle function.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file Plot of an exemplary raw data set. The upper graph shows the force-time data over an entire activation of all three stretch velocities (red: 1 vmax, blue: 0.1 vmax, black: 0.01 vmax). The middle graph shows the sarcomere length-time graph of the experiment. The lower graph shows the relative length-time data corresponding to the upper and middle graphs 1 (JPG 1454 KB)

Supplementary file Plot of enlarged first part of the stretch. The figure shows an enlargement of the initial stretch phase. To illustrate slope1, linear regressions between the start of the stretch and s2 were drawn (dashed lines), which are used to calculate slope1 according to [7] 2 (JPG 157 KB)

Supplementary file The force slopes increased with stretch velocity. Force slope units are maximum isometric force, F0, per optimal sarcomere length, ls0; stretch velocities are normalized to maximum contraction velocity, vmax. stretch velocity had a significant effect (p < 0.001) on the initial force slope (slope1, a) and that during the second half of the stretch (slope2, b).*** indicate significant (p < 0.001) differences between subgroups. For statistical details see Table 1 3 (JPG 1194 KB)

Supplementary file The force maxima increased with stretch velocities. Forces are normalized to maximum isometric force, F0, stretch velocities to maximum contraction velocity, vmax. stretch velocity had a significant effect (p < 0.001) on the force peak s2 (a), the corresponding sarcomere length ls2 (b), and the force at the end of the stretch, se (c). *** indicate significant (p < 0.001) differences between subgroups. For statistical details see Table 1 4 (JPG 1025 KB)

Supplementary file The local force minimum and Give increased with stretch velocity. Forces are normalized to maximum isometric force, F0, stretch velocities to maximum contraction velocity, vmax. stretch velocity had a significant effect (p < 0.001) on the local force minimum sg (a), the corresponding sarcomere length lsg (b), and the force decrease Give (c). *** indicate significant (p < 0.001) differences between subgroups. For statistical details see Table 1 5 (JPG 700 KB)

Author contribution

Conceptualization, T.S., A.T., and S.W.; methodology, S.W., A.T., T.S., and C.R.; formal analysis, S.W.; investigation, A.T. and S.W.; resources, T.S.; writing—original draft preparation, S.W. and A.T.; writing—review and editing, S.W., A.T., T.S., and C.R.; visualization, S.W.; funding acquisition, A.T., T.S., and C.R. All authors have read and agreed to the published version of the manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL. This research was funded by the Deutsche Forschungsgemeinschaft (DFG) under grants SI841/17–1 and RO5811/1–1 (project number 405834662) as well as partially funded by the DFG as part of the German Excellence Strategy—EXC 2075—390740016.

Data availability

The data presented in this study are available on request from the corresponding author.

Declarations

Ethics approval

The study was conducted according to the guidelines of ARRIVE and approved according to the German animal protection law (Tierschutzgesetz, §4 (3); Permit Number: T 201_21 ST).

Informed consent

Not applicable.

Conflict of interest

The authors declare no competing interests.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Abbott B Bigland B Ritchie J The cost of negative work J Physiol 1952 117 380 390 10.1113/jphysiol.1952.sp004755 14946742
Abbott B, Bigland B, Ritchie J (1952) The cost of negative work. J Physiol 117:380–39014946742 10.1113/jphysiol.1952.sp004755
2. Araz M Weidner S Izzi F Muscle preflex response to perturbations in locomotion: in vitro experiments and simulations with realistic boundary conditions Front Bioeng Biotechnol 2023 11 1 13 10.3389/fbioe.2023.1150170
Araz M, Weidner S, Izzi F et al (2023) Muscle preflex response to perturbations in locomotion: in vitro experiments and simulations with realistic boundary conditions. Front Bioeng Biotechnol 11:1–1310.3389/fbioe.2023.1150170
3. Arnold EM Hamner SR Seth A How muscle fiber lengths and velocities affect muscle force generation as humans walk and run at different speeds J Exp Biol 2013 216 2150 2160 23470656
Arnold EM, Hamner SR, Seth A et al (2013) How muscle fiber lengths and velocities affect muscle force generation as humans walk and run at different speeds. J Exp Biol 216:2150–216023470656
4. Ausoni S Gorza L Schiaffino S Expression of myosin heavy chain isoforms in stimulated fast and slow rat muscles J Neurosci 1990 10 153 160 10.1523/JNEUROSCI.10-01-00153.1990 2405110
Ausoni S, Gorza L, Schiaffino S et al (1990) Expression of myosin heavy chain isoforms in stimulated fast and slow rat muscles. J Neurosci 10:153–1602405110 10.1523/JNEUROSCI.10-01-00153.1990
5. Barclay CJ Mechanical efficiency and fatigue of fast and slow muscles of the mouse J Physiol 1996 497 781 794 10.1113/jphysiol.1996.sp021809 9003563
Barclay CJ (1996) Mechanical efficiency and fatigue of fast and slow muscles of the mouse. J Physiol 497:781–7949003563 10.1113/jphysiol.1996.sp021809
6. Bianco P Nagy A Kengyel A Interaction forces between F-actin and titin PEVK domain measured with optical tweezers Biophys J 2007 93 2102 2109 10.1529/biophysj.107.106153 17513381
Bianco P, Nagy A, Kengyel A et al (2007) Interaction forces between F-actin and titin PEVK domain measured with optical tweezers. Biophys J 93:2102–210917513381 10.1529/biophysj.107.106153
7. Burkholder TJ Lieber RL Review sarcomere length operating range of vertebrate muscles during movement J Exp Biol 2001 204 1529 1536 10.1242/jeb.204.9.1529 11296141
Burkholder TJ, Lieber RL (2001) Review sarcomere length operating range of vertebrate muscles during movement. J Exp Biol 204:1529–153611296141 10.1242/jeb.204.9.1529
8. Capitanio M Canepari M Cacciafesta P Two independent mechanical events in the interaction cycle of skeletal muscle myosin with actin Proc Natl Acad Sci U S A 2006 103 87 92 10.1073/pnas.0506830102 16371472
Capitanio M, Canepari M, Cacciafesta P et al (2006) Two independent mechanical events in the interaction cycle of skeletal muscle myosin with actin. Proc Natl Acad Sci U S A 103:87–9216371472 10.1073/pnas.0506830102
9. Cecchi G Bagni MA Cecchini E Crossbridge viscosity in activated frog muscle fibres Biophys Chem 1997 68 1 8 10.1016/S0301-4622(96)02263-6 9468605
Cecchi G, Bagni MA, Cecchini E et al (1997) Crossbridge viscosity in activated frog muscle fibres. Biophys Chem 68:1–89468605 10.1016/S0301-4622(96)02263-6
10. Choi SJ Widrick JJ Calcium-activated force of human muscle fibers following a standardized eccentric contraction Am J Physiol - Cell Physiol 2010 299 1409 1417 10.1152/ajpcell.00226.2010
Choi SJ, Widrick JJ (2010) Calcium-activated force of human muscle fibers following a standardized eccentric contraction. Am J Physiol - Cell Physiol 299:1409–141710.1152/ajpcell.00226.2010
11. Chung CS Methawasin M Nelson OL Titin based viscosity in ventricular physiology: an integrative investigation of PEVK–actin interactions J Mol Cell Cardiol 2011 51 428 434 10.1016/j.yjmcc.2011.06.006 21708170
Chung CS, Methawasin M, Nelson OL et al (2011) Titin based viscosity in ventricular physiology: an integrative investigation of PEVK–actin interactions. J Mol Cell Cardiol 51:428–43421708170 10.1016/j.yjmcc.2011.06.006
12. Close RI Dynamic properties of fast and slow skeletal muscles of the rat after nerve cross-union J Physiol 1969 204 331 346 10.1113/jphysiol.1969.sp008916 5824642
Close RI (1969) Dynamic properties of fast and slow skeletal muscles of the rat after nerve cross-union. J Physiol 204:331–3465824642 10.1113/jphysiol.1969.sp008916
13. Cohen J Statistical Power Analysis for the Behavioral Sciences 2013 Routledge
Cohen J (2013) Statistical Power Analysis for the Behavioral Sciences. Routledge. 10.4324/9780203771587
14. Curtin NA Diack RA West TG Skinned fibres produce the same power and force as intact fibre bundles from muscle of wild rabbits J Exp Biol 2015 218 2856 2863 26206354
Curtin NA, Diack RA, West TG et al (2015) Skinned fibres produce the same power and force as intact fibre bundles from muscle of wild rabbits. J Exp Biol 218:2856–286326206354
15. Damsgaard M Rasmussen J Christensen ST Analysis of musculoskeletal systems in the AnyBody Modeling System Simul Model Pract Theory 2006 14 1100 1111 10.1016/j.simpat.2006.09.001
Damsgaard M, Rasmussen J, Christensen ST et al (2006) Analysis of musculoskeletal systems in the AnyBody Modeling System. Simul Model Pract Theory 14:1100–111110.1016/j.simpat.2006.09.001
16. Degens H Yu F Li X Effects of age and gender on shortening velocity and myosin isoforms in single rat muscle fibres Acta Physiol Scand 1998 163 33 40 10.1046/j.1365-201x.1998.00329.x 9648621
Degens H, Yu F, Li X et al (1998) Effects of age and gender on shortening velocity and myosin isoforms in single rat muscle fibres. Acta Physiol Scand 163:33–409648621 10.1046/j.1365-201x.1998.00329.x
17. Dutta S Tsiros C Sundar SL Calcium increases titin N2A binding to F-actin and regulated thin filaments Sci Rep 2018 8 1 11 10.1038/s41598-018-32952-8 29311619
Dutta S, Tsiros C, Sundar SL et al (2018) Calcium increases titin N2A binding to F-actin and regulated thin filaments. Sci Rep 8:1–1129311619 10.1038/s41598-018-32952-8
18. Fenn WO The relation between the work performed and the energy liberated in muscular contraction J Physiol 1924 58 373 395 10.1113/jphysiol.1924.sp002141 16993634
Fenn WO (1924) The relation between the work performed and the energy liberated in muscular contraction. J Physiol 58:373–39516993634 10.1113/jphysiol.1924.sp002141
19. Flitney FW Hirst DG Cross-bridge detachment and sarcomere ‘give’ during stretch of active frog’s muscle J Physiol 1978 276 449 465 10.1113/jphysiol.1978.sp012246 306433
Flitney FW, Hirst DG (1978) Cross-bridge detachment and sarcomere ‘give’ during stretch of active frog’s muscle. J Physiol 276:449–465306433 10.1113/jphysiol.1978.sp012246
20. Freiburg A Trombitas K Hell W Series of exon-skipping events in the elastic spring region of titin as the structural basis for myofibrillar elastic diversity Circ Res 2000 86 11 1114 1121 10.1161/01.RES.86.11.1114 10850961
Freiburg A, Trombitas K, Hell W et al (2000) Series of exon-skipping events in the elastic spring region of titin as the structural basis for myofibrillar elastic diversity. Circ Res 86(11):1114–112110850961 10.1161/01.RES.86.11.1114
21. Freundt JK Linke WA Titin as a force-generating muscle protein under regulatory control J Appl Physiol 2019 126 1474 1482 10.1152/japplphysiol.00865.2018 30521425
Freundt JK, Linke WA (2019) Titin as a force-generating muscle protein under regulatory control. J Appl Physiol 126:1474–148230521425 10.1152/japplphysiol.00865.2018
22. Getz EB Cooke R Lehman SL Phase transition in force during ramp stretches of skeletal muscle Biophys J 1998 75 2971 2983 10.1016/S0006-3495(98)77738-0 9826617
Getz EB, Cooke R, Lehman SL (1998) Phase transition in force during ramp stretches of skeletal muscle. Biophys J 75:2971–29839826617 10.1016/S0006-3495(98)77738-0
23. Greaser ML Pleitner JM Titin isoform size is not correlated with thin filament length in rat skeletal muscle Front Physiol 2014 5 1 9 10.3389/fphys.2014.00035 24478714
Greaser ML, Pleitner JM (2014) Titin isoform size is not correlated with thin filament length in rat skeletal muscle. Front Physiol 5:1–924478714 10.3389/fphys.2014.00035
24. Griffiths PJ Güth K Kuhn HJ Cross bridge slippage in skinned frog muscle fibres Biophys Struct Mech 1980 7 107 124 10.1007/BF00538402 6971661
Griffiths PJ, Güth K, Kuhn HJ et al (1980) Cross bridge slippage in skinned frog muscle fibres. Biophys Struct Mech 7:107–1246971661 10.1007/BF00538402
25. Haeufle DFB Günther M Bayer A Hill-type muscle model with serial damping and eccentric force-velocity relation J Biomech 2014 47 1531 1536 10.1016/j.jbiomech.2014.02.009 24612719
Haeufle DFB, Günther M, Bayer A et al (2014) Hill-type muscle model with serial damping and eccentric force-velocity relation. J Biomech 47:1531–153624612719 10.1016/j.jbiomech.2014.02.009
26. Heidlauf T Klotz T Rode C A multi-scale continuum model of skeletal muscle mechanics predicting force enhancement based on actin – titin interaction Biomech Model Mechanobiol 2016 15 1423 1437 10.1007/s10237-016-0772-7 26935301
Heidlauf T, Klotz T, Rode C et al (2016) A multi-scale continuum model of skeletal muscle mechanics predicting force enhancement based on actin – titin interaction. Biomech Model Mechanobiol 15:1423–143726935301 10.1007/s10237-016-0772-7
27. Herzog W Leonard TR The history dependence of force production in mammalian skeletal muscle following stretch-shortening and shortening-stretch cycles J Biomech 2000 33 531 542 10.1016/S0021-9290(99)00221-3 10708773
Herzog W, Leonard TR (2000) The history dependence of force production in mammalian skeletal muscle following stretch-shortening and shortening-stretch cycles. J Biomech 33:531–54210708773 10.1016/S0021-9290(99)00221-3
28. Herzog JA Leonard TR Jinha A Titin (visco-) elasticity in skeletal muscle myofibrils Mol Cell Biomech 2014 11 1 17 25330621
Herzog JA, Leonard TR, Jinha A et al (2014) Titin (visco-) elasticity in skeletal muscle myofibrils. Mol Cell Biomech 11:1–1725330621
29. Hessel AL, Monroy JA, Nishikawa KC (2021) Non-cross bridge viscoelastic elements contribute to muscle force and work during stretch-shortening cycles: evidence from whole muscles and permeabilized fibers. Front Physiol:12. 10.3389/fphys.2021.648019
30. Hilber K Galler S Improvement of the measurements on skinned muscle fibres by fixation of the fibre ends with glutaraldehyde J Muscle Res Cell Motil 1998 19 365 372 10.1023/A:1005393519811 9635279
Hilber K, Galler S (1998) Improvement of the measurements on skinned muscle fibres by fixation of the fibre ends with glutaraldehyde. J Muscle Res Cell Motil 19:365–3729635279 10.1023/A:1005393519811
31. Hody S Croisier JL Bury T Eccentric muscle contractions: risks and benefits Front Physiol 2019 10 1 18 10.3389/fphys.2019.00536 30723415
Hody S, Croisier JL, Bury T et al (2019) Eccentric muscle contractions: risks and benefits. Front Physiol 10:1–1830723415 10.3389/fphys.2019.00536
32. Horowits R Passive force generation and titin isoforms in mammalian skeletal muscle Biophys J 1992 61 392 398 10.1016/S0006-3495(92)81845-3 1547327
Horowits R (1992) Passive force generation and titin isoforms in mammalian skeletal muscle. Biophys J 61:392–3981547327 10.1016/S0006-3495(92)81845-3
33. Huxley AF Muscle structure and theories of contraction Prog Biophys Biophys Chem 1957 7 255 318 10.1016/S0096-4174(18)30128-8 13485191
Huxley AF (1957) Muscle structure and theories of contraction. Prog Biophys Biophys Chem 7:255–31813485191 10.1016/S0096-4174(18)30128-8
34. Huxley HE Hanson J Changes in the cross-striations of muscle during contraction and stretch and their structural interpretation Nature 1954 173 973 976 10.1038/173973a0 13165698
Huxley HE, Hanson J (1954) Changes in the cross-striations of muscle during contraction and stretch and their structural interpretation. Nature 173:973–97613165698 10.1038/173973a0
35. Huxley AF Niedergerke R Structural changes in muscle during contraction; interference microscopy of living muscle fibres Nature 1954 173 971 973 10.1038/173971a0 13165697
Huxley AF, Niedergerke R (1954) Structural changes in muscle during contraction; interference microscopy of living muscle fibres. Nature 173:971–97313165697 10.1038/173971a0
36. Josephson RK Stokes DR The force-velocity properties of a crustacean muscle during lengthening J Exp Biol 1999 202 593 607 10.1242/jeb.202.5.593 9929461
Josephson RK, Stokes DR (1999) The force-velocity properties of a crustacean muscle during lengthening. J Exp Biol 202:593–6079929461 10.1242/jeb.202.5.593
37. Kellermayer MSZ Smith SB Granzier HL Folding-unfolding transitions in single titin molecules characterized with laser tweezers Science 1997 276 1112 1116 10.1126/science.276.5315.1112 9148805
Kellermayer MSZ, Smith SB, Granzier HL et al (1997) Folding-unfolding transitions in single titin molecules characterized with laser tweezers. Science 276:1112–11169148805 10.1126/science.276.5315.1112
38. Krylow AM Sandercock TG Dynamic force responses of muscle involving eccentric contraction J Biomech 1997 30 27 33 10.1016/S0021-9290(96)00097-8 8970921
Krylow AM, Sandercock TG (1997) Dynamic force responses of muscle involving eccentric contraction. J Biomech 30:27–338970921 10.1016/S0021-9290(96)00097-8
39. Labeit S Kolmerer B Titins: giant proteins in charge of muscle ultrastructure and elasticity Science 1995 270 293 296 10.1126/science.270.5234.293 7569978
Labeit S, Kolmerer B (1995) Titins: giant proteins in charge of muscle ultrastructure and elasticity. Science 270:293–2967569978 10.1126/science.270.5234.293
40. Leigh JH Kinnear TC On interaction classification Educ Psychol Meas 1980 40 841 843 10.1177/001316448004000405
Leigh JH, Kinnear TC (1980) On interaction classification. Educ Psychol Meas 40:841–84310.1177/001316448004000405
41. Li H Linke WA Oberhauser AF Reverse engineering of the giant muscle protein titin Nature 2002 418 998 1002 10.1038/nature00938 12198551
Li H, Linke WA, Oberhauser AF et al (2002) Reverse engineering of the giant muscle protein titin. Nature 418:998–100212198551 10.1038/nature00938
42. Linari M Caremani M Piperio C Stiffness and fraction of myosin motors responsible for active force in permeabilized muscle fibers from rabbit psoas Biophys J 2007 92 2476 2490 10.1529/biophysj.106.099549 17237201
Linari M, Caremani M, Piperio C et al (2007) Stiffness and fraction of myosin motors responsible for active force in permeabilized muscle fibers from rabbit psoas. Biophys J 92:2476–249017237201 10.1529/biophysj.106.099549
43. Lindstedt SL LaStayo PC Reich TE When active muscles lengthen: properties and consequences of eccentric contractions News Physiol Sci 2001 16 256 261 11719600
Lindstedt SL, LaStayo PC, Reich TE (2001) When active muscles lengthen: properties and consequences of eccentric contractions. News Physiol Sci 16:256–26111719600
44. Linke WA Titin gene and protein functions in passive and active muscle Annu Rev Physiol 2018 80 1 389 411 10.1146/annurev-physiol-021317-121234 29131758
Linke WA (2018) Titin gene and protein functions in passive and active muscle. Annu Rev Physiol 80(1):389–411. 10.1146/annurev-physiol-021317-12123429131758 10.1146/annurev-physiol-021317-121234
45. Lombardi V Piazzesi G The contractile response during steady lengthening of stimulated frog muscle fibres J Physiol 1990 431 141 171 10.1113/jphysiol.1990.sp018324 2100305
Lombardi V, Piazzesi G (1990) The contractile response during steady lengthening of stimulated frog muscle fibres. J Physiol 431:141–1712100305 10.1113/jphysiol.1990.sp018324
46. Mártonfalvi Z Bianco P Linari M Low-force transitions in single titin molecules reflect a memory of contractile history J Cell Sci 2014 127 858 870 24357719
Mártonfalvi Z, Bianco P, Linari M et al (2014) Low-force transitions in single titin molecules reflect a memory of contractile history. J Cell Sci 127:858–87024357719
47. Meyer K Steiner R Lastayo P Eccentric exercise in coronary patients: central hemodynamic and metabolic responses Med Sci Sports Exerc 2003 35 1076 1082 10.1249/01.MSS.0000074580.79648.9D 12840625
Meyer K, Steiner R, Lastayo P et al (2003) Eccentric exercise in coronary patients: central hemodynamic and metabolic responses. Med Sci Sports Exerc 35:1076–108212840625 10.1249/01.MSS.0000074580.79648.9D
48. Millard M, Franklin DW, Herzog W (2023) A three filament mechanistic model of musculotendon force and impedance. 10.1101/2023.03.27.534347
49. Mutungi G Ranatunga KW The visco-elasticity muscle fibres of resting intact, mammalian should be addressed J Muscle Res Cell Motil 1996 17 357 364 10.1007/BF00240933 8814555
Mutungi G, Ranatunga KW (1996) The visco-elasticity muscle fibres of resting intact, mammalian should be addressed. J Muscle Res Cell Motil 17:357–3648814555 10.1007/BF00240933
50. Neagoe C Opitz CA Makarenko I Gigantic variety: expression patterns of titin isoforms in striated muscles and consequences for myofibrillar passive stiffness J Muscle Res Cell Motil 2003 24 175 189 10.1023/A:1026053530766 14609029
Neagoe C, Opitz CA, Makarenko I et al (2003) Gigantic variety: expression patterns of titin isoforms in striated muscles and consequences for myofibrillar passive stiffness. J Muscle Res Cell Motil 24:175–18914609029 10.1023/A:1026053530766
51. Nishikawa K Titin: a tunable spring in active muscle Physiology 2020 35 209 217 10.1152/physiol.00036.2019 32293234
Nishikawa K (2020) Titin: a tunable spring in active muscle. Physiology 35:209–21732293234 10.1152/physiol.00036.2019
52. Nishikawa KC Monroy JA Uyeno TE Is titin a ‘ winding filament ’? A new twist on muscle contraction Proc Royal Soc B: Biol Sci 2012 279 981 990 10.1098/rspb.2011.1304
Nishikawa KC, Monroy JA, Uyeno TE et al (2012) Is titin a ‘ winding filament ’? A new twist on muscle contraction. Proc Royal Soc B: Biol Sci 279:981–99010.1098/rspb.2011.1304
53. Nishikawa K Dutta S DuVall M Calcium-dependent titin–thin filament interactions in muscle: observations and theory J Muscle Res Cell Motil 2019 41 1 125 139 10.1007/s10974-019-09540-y 31289970
Nishikawa K, Dutta S, DuVall M et al (2019) Calcium-dependent titin–thin filament interactions in muscle: observations and theory. J Muscle Res Cell Motil 41(1):125–139. 10.1007/s10974-019-09540-y31289970 10.1007/s10974-019-09540-y
54. Pinniger GJ Ranatunga KW Offer GW Crossbridge and non-crossbridge contributions to tension in lengthening rat muscle: force-induced reversal of the power stroke J Physiol 2006 573 627 643 10.1113/jphysiol.2005.095448 16627571
Pinniger GJ, Ranatunga KW, Offer GW (2006) Crossbridge and non-crossbridge contributions to tension in lengthening rat muscle: force-induced reversal of the power stroke. J Physiol 573:627–64316627571 10.1113/jphysiol.2005.095448
55. Prado LG Makarenko I Andresen C Isoform diversity of giant proteins in relation to passive and active contractile properties of rabbit skeletal muscles J Gen Physiol 2005 126 461 480 10.1085/jgp.200509364 16230467
Prado LG, Makarenko I, Andresen C et al (2005) Isoform diversity of giant proteins in relation to passive and active contractile properties of rabbit skeletal muscles. J Gen Physiol 126:461–48016230467 10.1085/jgp.200509364
56. Ranatunga KW Temperature-dependence of shortening velocity and rate of isometric tension development in rat skeletal muscle J Physiol 1982 329 465 483 10.1113/jphysiol.1982.sp014314 7143257
Ranatunga KW (1982) Temperature-dependence of shortening velocity and rate of isometric tension development in rat skeletal muscle. J Physiol 329:465–4837143257 10.1113/jphysiol.1982.sp014314
57. Ranatunga KW The force-velocity relation of rat fast- and slow-twitch muscles examined at different temperatures J Physiol 1984 351 517 529 10.1113/jphysiol.1984.sp015260 6747875
Ranatunga KW (1984) The force-velocity relation of rat fast- and slow-twitch muscles examined at different temperatures. J Physiol 351:517–5296747875 10.1113/jphysiol.1984.sp015260
58. Rasmussen J Vondrak V Damsgaard M The AnyBody project - Computer analysis of the human body Comput Anal Hum Body Biomech Man 2002 1 270 274
Rasmussen J, Vondrak V, Damsgaard M et al (2002) The AnyBody project - Computer analysis of the human body. Comput Anal Hum Body Biomech Man 1:270–274
59. Rode C Siebert T The effects of parallel and series elastic components on estimated active cat soleus muscle force J Mech Med Biol 2009 9 105 122 10.1142/S0219519409002870
Rode C, Siebert T (2009) The effects of parallel and series elastic components on estimated active cat soleus muscle force. J Mech Med Biol 9:105–12210.1142/S0219519409002870
60. Rode C Siebert T Blickhan R Titin-induced force enhancement and force depression: a ‘sticky-spring’ mechanism in muscle contractions? J Theor Biol 2009 259 350 360 10.1016/j.jtbi.2009.03.015 19306884
Rode C, Siebert T, Blickhan R (2009) Titin-induced force enhancement and force depression: a ‘sticky-spring’ mechanism in muscle contractions? J Theor Biol 259:350–36019306884 10.1016/j.jtbi.2009.03.015
61. Röhrle O Sprenger M Schmitt S A two-muscle, continuum-mechanical forward simulation of the upper limb Biomech Model Mechanobiol 2017 16 743 762 10.1007/s10237-016-0850-x 27837360
Röhrle O, Sprenger M, Schmitt S (2017) A two-muscle, continuum-mechanical forward simulation of the upper limb. Biomech Model Mechanobiol 16:743–76227837360 10.1007/s10237-016-0850-x
62. Roots H Ranatunga KW An analysis of the temperature dependence of force, during steady shortening at different velocities, in (mammalian) fast muscle fibres J Muscle Res Cell Motil 2008 29 9 24 10.1007/s10974-008-9138-9 18523851
Roots H, Ranatunga KW (2008) An analysis of the temperature dependence of force, during steady shortening at different velocities, in (mammalian) fast muscle fibres. J Muscle Res Cell Motil 29:9–2418523851 10.1007/s10974-008-9138-9
63. Roots H Offer GW Ranatunga KW Comparison of the tension responses to ramp shortening and lengthening in intact mammalian muscle fibres: crossbridge and non-crossbridge contributions J Muscle Res Cell Motil 2007 28 123 139 10.1007/s10974-007-9110-0 17610136
Roots H, Offer GW, Ranatunga KW (2007) Comparison of the tension responses to ramp shortening and lengthening in intact mammalian muscle fibres: crossbridge and non-crossbridge contributions. J Muscle Res Cell Motil 28:123–13917610136 10.1007/s10974-007-9110-0
64. Rüegg JC Calcium in Muscle Contraction 1992 2 Berlin, Heidelberg Springer
Rüegg JC (1992) Calcium in Muscle Contraction, 2nd edn. Springer, Berlin, Heidelberg. 10.1007/978-3-642-77560-4
65. Schaeffer PJ Lindstedt SL How animals move: comparative lessons on animal locomotion Compr Physiol 2013 3 289 314 10.1002/cphy.c110059 23720288
Schaeffer PJ, Lindstedt SL (2013) How animals move: comparative lessons on animal locomotion. Compr Physiol 3:289–31423720288 10.1002/cphy.c110059
66. Schappacher-Tilp G Leonard T Desch G A novel three-filament model of force generation in eccentric contraction of skeletal muscles PLoS ONE 2015 10 e0117634 10.1371/journal.pone.0117634 25816319
Schappacher-Tilp G, Leonard T, Desch G et al (2015) A novel three-filament model of force generation in eccentric contraction of skeletal muscles. PLoS ONE 10:e011763425816319 10.1371/journal.pone.0117634
67. Seth A Hicks JL Uchida TK OpenSim: simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement PLOS Comput Biol 2018 14 e1006223 10.1371/journal.pcbi.1006223 30048444
Seth A, Hicks JL, Uchida TK et al (2018) OpenSim: simulating musculoskeletal dynamics and neuromuscular control to study human and animal movement. PLOS Comput Biol 14:e100622330048444 10.1371/journal.pcbi.1006223
68. Seydewitz R Siebert T Böl M On a three-dimensional constitutive model for history effects in skeletal muscles Biomech Model Mechanobiol 2019 18 6 1665 1681 10.1007/s10237-019-01167-9 31102082
Seydewitz R, Siebert T, Böl M (2019) On a three-dimensional constitutive model for history effects in skeletal muscles. Biomech Model Mechanobiol 18(6):1665–1681. 10.1007/s10237-019-01167-931102082 10.1007/s10237-019-01167-9
69. Siebert T Rode C Herzog W Nonlinearities make a difference: comparison of two common Hill-type models with real muscle Biol Cybern 2008 98 133 143 10.1007/s00422-007-0197-6 18049823
Siebert T, Rode C, Herzog W et al (2008) Nonlinearities make a difference: comparison of two common Hill-type models with real muscle. Biol Cybern 98:133–14318049823 10.1007/s00422-007-0197-6
70. Siebert T Leichsenring K Rode C Three-dimensional muscle architecture and comprehensive dynamic properties of rabbit gastrocnemius, plantaris and soleus: input for simulation studies PLoS ONE 2015 10 e0130985 10.1371/journal.pone.0130985 26114955
Siebert T, Leichsenring K, Rode C et al (2015) Three-dimensional muscle architecture and comprehensive dynamic properties of rabbit gastrocnemius, plantaris and soleus: input for simulation studies. PLoS ONE 10:e013098526114955 10.1371/journal.pone.0130985
71. Soukup T Zachařová G Smerdu V Fibre type composition of soleus and extensor digitorum longus muscles in normal female inbred Lewis rats Acta Histochem 2002 104 399 405 10.1078/0065-1281-00660 12553710
Soukup T, Zachařová G, Smerdu V (2002) Fibre type composition of soleus and extensor digitorum longus muscles in normal female inbred Lewis rats. Acta Histochem 104:399–40512553710 10.1078/0065-1281-00660
72. Stephenson DG Williams DA Effects of sarcomere length on the force-pCa relation in fast- and slow-twitch skinned muscle fibres from the rat J Physiol 1982 333 637 653 10.1113/jphysiol.1982.sp014473 7182478
Stephenson DG, Williams DA (1982) Effects of sarcomere length on the force-pCa relation in fast- and slow-twitch skinned muscle fibres from the rat. J Physiol 333:637–6537182478 10.1113/jphysiol.1982.sp014473
73. Stienen GJM Versteeg PG Papp Z Mechanical properties of skinned rabbit psoas and soleus muscle fibres during lengthening: effects of phosphate and Ca2+ J Physiol 1992 451 503 523 10.1113/jphysiol.1992.sp019176 1403822
Stienen GJM, Versteeg PG, Papp Z et al (1992) Mechanical properties of skinned rabbit psoas and soleus muscle fibres during lengthening: effects of phosphate and Ca2+. J Physiol 451:503–5231403822 10.1113/jphysiol.1992.sp019176
74. Tahir U, Monroy JA, Rice NA et al (2020) Effects of a titin mutation on force enhancement and force depression in mouse soleus muscles. J Exp Biol:223. 10.1242/jeb.197038
75. Till O Siebert T Rode C Characterization of isovelocity extension of activated muscle: a Hill-type model for eccentric contractions and a method for parameter determination J Theor Biol 2008 255 176 187 10.1016/j.jtbi.2008.08.009 18771670
Till O, Siebert T, Rode C et al (2008) Characterization of isovelocity extension of activated muscle: a Hill-type model for eccentric contractions and a method for parameter determination. J Theor Biol 255:176–18718771670 10.1016/j.jtbi.2008.08.009
76. Tomalka A Eccentric muscle contractions: from single muscle fibre to whole muscle mechanics Pflügers Arch - Eur J Physiol 2023 475 4 421 435 10.1007/s00424-023-02794-z 36790515
Tomalka A (2023) Eccentric muscle contractions: from single muscle fibre to whole muscle mechanics. Pflügers Arch - Eur J Physiol 475(4):421–435. 10.1007/s00424-023-02794-z36790515 10.1007/s00424-023-02794-z
77. Tomalka A Rode C Schumacher J The active force-length relationship is invisible during extensive eccentric contractions in skinned skeletal muscle fibres Proc R Soc B Biol Sci 2017 284 1854 20162497 10.1098/rspb.2016.2497
Tomalka A, Rode C, Schumacher J et al (2017) The active force-length relationship is invisible during extensive eccentric contractions in skinned skeletal muscle fibres. Proc R Soc B Biol Sci 284:1854–20162497. 10.1098/rspb.2016.249710.1098/rspb.2016.2497
78. Tomalka A Weidner S Hahn D Cross-bridges and sarcomeric non-cross-bridge structures contribute to increased work in stretch-shortening cycles Front Physiol 2020 11 1 14 10.3389/fphys.2020.00921 32038307
Tomalka A, Weidner S, Hahn D et al (2020) Cross-bridges and sarcomeric non-cross-bridge structures contribute to increased work in stretch-shortening cycles. Front Physiol 11:1–1432038307 10.3389/fphys.2020.00921
79. Tomalka A Weidner S Hahn D Power amplification increases with contraction velocity during stretch-shortening cycles of skinned muscle fibers Front Physiol 2021 12 1 16 10.3389/fphys.2021.644981
Tomalka A, Weidner S, Hahn D et al (2021) Power amplification increases with contraction velocity during stretch-shortening cycles of skinned muscle fibers. Front Physiol 12:1–1610.3389/fphys.2021.644981
80. Trombitás K Greaser M Labeit S Titin extensibility in situ: entropic elasticity of permanently folded and permanently unfolded molecular segments J Cell Biol 1998 140 853 859 10.1083/jcb.140.4.853 9472037
Trombitás K, Greaser M, Labeit S et al (1998) Titin extensibility in situ: entropic elasticity of permanently folded and permanently unfolded molecular segments. J Cell Biol 140:853–8599472037 10.1083/jcb.140.4.853
81. van Soest AJ Bobbert MF The contribution of muscle properties in the control of explosive movements Biol Cybern 1993 69 195 204 10.1007/BF00198959 8373890
van Soest AJ, Bobbert MF (1993) The contribution of muscle properties in the control of explosive movements. Biol Cybern 69:195–2048373890 10.1007/BF00198959
82. Wang K McClure J Tu A Titin: major myofibrillar components of striated muscle Proc Natl Acad Sci 1979 76 3698 3702 10.1073/pnas.76.8.3698 291034
Wang K, McClure J, Tu A (1979) Titin: major myofibrillar components of striated muscle. Proc Natl Acad Sci 76:3698–3702291034 10.1073/pnas.76.8.3698
83. Wang K Mccarter R Wright J Regulation of skeletal muscle stiffness and elasticity by titin isoforms: a test of the segmental extension model of resting tension Proc Natl Acad Sci U S A 1991 88 7101 7105 10.1073/pnas.88.16.7101 1714586
Wang K, Mccarter R, Wright J et al (1991) Regulation of skeletal muscle stiffness and elasticity by titin isoforms: a test of the segmental extension model of resting tension. Proc Natl Acad Sci U S A 88:7101–71051714586 10.1073/pnas.88.16.7101
84. Ward SR Winters TM O’Connor SM Non-linear scaling of passive mechanical properties in fibers, bundles, fascicles and whole rabbit muscles Front Physiol 2020 11 1 9 10.3389/fphys.2020.00211 32038307
Ward SR, Winters TM, O’Connor SM et al (2020) Non-linear scaling of passive mechanical properties in fibers, bundles, fascicles and whole rabbit muscles. Front Physiol 11:1–932038307 10.3389/fphys.2020.00211
85. Weidner S Tomalka A Rode C How velocity impacts eccentric force generation of fully activated skinned skeletal muscle fibers in long stretches J Appl Physiol 2022 133 223 233 10.1152/japplphysiol.00735.2021 35652830
Weidner S, Tomalka A, Rode C et al (2022) How velocity impacts eccentric force generation of fully activated skinned skeletal muscle fibers in long stretches. J Appl Physiol 133:223–23335652830 10.1152/japplphysiol.00735.2021
