Muscle force depends on the amount of transversal muscle loading

Note on hydrostatic skeletons: muscles operating within a pressurized environment

Muscles and muscle fibers are volume-constant constructs that deform when contracted and develop internal pressures. However, muscles embedded in hydrostatic skeletons are also exposed to external pressures generated by their activity. For two examples, the pressure generation in spiders and in annelids, we used simplified biomechanical models to demonstrate that high intracellular pressures diminishing the resulting tensile stress of the muscle fibers are avoided in the hydrostatic skeleton. The findings are relevant for a better understanding of the design and functionality of biological hydrostatic skeletons.

Physical exercise to promote regeneration after peripheral nerve injury in animal models: A systematic review

Peripheral nerve injuries are common injuries that often have a drastic effect on patient's activities of daily living and physical function. While techniques for the surgical repair of these injuries have improved over time, rehabilitation methods following these repairs have been non-standardized and under researched. Electronic searches were conducted in Ovid/Medline and SCOPUS to identify articles that discuss rehabilitation and exercise following peripheral nerve injury in animal models and its effects on peripheral nerve regeneration and recovery of function. Thirty-eight articles met inclusion criteria; all were prospective studies in animal models. This systematic review suggests that exercise is a useful tool in returning autonomy to the individual and has beneficial effects in the recovery from peripheral nerve injury. It is beneficial to use rehabilitation exercises following the repair of peripheral nerve injuries to promote regeneration, and timing of that exercise may be just as important as the exercise prescribed. However, further studies with standardized models and outcome measures need to be conducted before translation to clinical trials.

A surgical technique for individual control of the muscles of the rabbit lower hindlimb

Little is known regarding the precise muscle, bone and joint actions resulting from individual and simultaneous muscle activation(s) of the lower limb. An in situ experimental approach is described herein to control the muscles of the rabbit lower hindlimb, including the medial and lateral gastrocnemius, soleus, plantaris and tibialis anterior. The muscles were stimulated using nerve-cuff electrodes placed around the innervating nerves of each muscle. Animals were fixed in a stereotactic frame with the ankle angle set at 90 deg. To demonstrate the efficacy of the experimental technique, isometric plantarflexion torque was measured at the 90 deg ankle joint angle at a stimulation frequency of 100, 60 and 30 Hz. Individual muscle torque and the torque produced during simultaneous activation of all plantarflexor muscles are presented for four animals. These results demonstrate that the experimental approach was reliable, with insignificant variation in torque between repeated contractions. The experimental approach described herein provides the potential for measuring a diverse array of muscle properties, which is important to improve our understanding of musculoskeletal biomechanics.

Finite Element Model of the Shoulder with Active Rotator Cuff Muscles: Application to Wheelchair Propulsion

The rotator cuff is prone to injury, remarkably so for manual wheelchair users. To understand its pathomechanisms, finite element models incorporating three-dimensional activated muscles are needed to predict soft tissue strains during given tasks. This study aimed to develop such a model to understand pathomechanisms associated with wheelchair propulsion. We developed an active muscle model associating a passive fiber-reinforced isotropic matrix with an activation law linking calcium ion concentration to tissue tension. This model was first evaluated against known physiological muscle behavior; then used to activate the rotator cuff during a wheelchair propulsion cycle. Here, experimental kinematics and electromyography data was used to drive a shoulder finite element model. Finally, we evaluated the importance of muscle activation by comparing the results of activated and non-activated rotator cuff muscles during both propulsion and isometric contractions. Qualitatively, the muscle constitutive law reasonably reproduced the classical Hill model force-length curve and the behavior of a transversally loaded muscle. During wheelchair propulsion, the deformation and fiber stretch of the supraspinatus muscle-tendon unit pointed towards the possibility for this tendon to develop tendinosis due to the multiaxial loading imposed by the kinematics of propulsion. Finally, differences in local stretch and positions of the lines of action between activated and non-activated models were only observed at activation levels higher than 30%. Our novel finite element model with active muscles is a promising tool for understanding the pathomechanisms of the rotator cuff for various dynamic tasks, especially those with high muscle activation levels.

Influence of muscle length on the three-dimensional architecture and aponeurosis dimensions of rabbit calf muscles

The function of a muscle is highly dependent on its architecture, which is characterized by the length, pennation, and curvature of the fascicles, and the geometry of the aponeuroses. During in vivo function, muscles regularly undergo changes in length, thereby altering their architecture. During passive muscle lengthening, fascicle length (FL) generally increases and the angle of fascicle pennation (FP) and the fascicle curvature (FC) decrease, while the aponeuroses increase in length but decrease in width. Muscles are differently structured, making their change during muscle lengthening complex and multifaceted. To obtain comprehensive data on architectural changes in muscles during passive length, the present study determined the three-dimensional fascicle geometry of rabbit M. gastrocnemius medialis (GM), M. gastrocnemius lateralis (GL), and M. plantaris (PLA). For this purpose, the left and right legs of three rabbits were histologically fixed at targeted ankle joint angles of 95° (short muscle length [SML]) and 60° (long muscle length [LML]), respectively, and the fascicles were tracked by manual three-dimensional digitization. In a second set of experiments, the GM aponeurosis dimensions of ten legs from five rabbits were determined at varying muscle lengths via optical marker tracking. The GM consisted of a uni-pennated compartment, whereas the GL and PLA contained multiple compartments of differently pennated fascicles. In the LML compared to the SML, the GM, GL, and PLA had on average a 41%, 29%, and 41% increased fascicle length, and a 30%, 25%, and 33% decrease in fascicle pennation and a 32%, 11%, and 35% decrease in fascicle curvature, respectively. Architectural properties were also differentiated among the different compartments of the PLA and GL, allowing for a more detailed description of their fascicle structure and changes. It was shown that the compartments change differently with muscle length. It was also shown that for each degree of ankle joint angle reduction, the proximal GM aponeurosis length increased by 0.11%, the aponeurosis width decreased by 0.22%, and the area was decreased by 0.20%. The data provided improve our understanding of muscles and can be used to develop and validate muscle models.

Biomechanical and Cortical Control of Tongue Movements During Chewing and Swallowing

Tongue function is vital for chewing and swallowing and lingual dysfunction is often associated with dysphagia. Better treatment of dysphagia depends on a better understanding of hyolingual morphology, biomechanics, and neural control in humans and animal models. Recent research has revealed significant variation among animal models in morphology of the hyoid chain and suprahyoid muscles which may be associated with variation in swallowing mechanisms. The recent deployment of XROMM (X-ray Reconstruction of Moving Morphology) to quantify 3D hyolingual kinematics has revealed new details on flexion and roll of the tongue during chewing in animal models, movements similar to those used by humans. XROMM-based studies of swallowing in macaques have falsified traditional hypotheses of mechanisms of tongue base retraction during swallowing, and literature review suggests that other animal models may employ a diversity of mechanisms of tongue base retraction. There is variation among animal models in distribution of hyolingual proprioceptors but how that might be related to lingual mechanics is unknown. In macaque monkeys, tongue kinematics-shape and movement-are strongly encoded in neural activity in orofacial primary motor cortex, giving optimism for development of brain-machine interfaces for assisting recovery of lingual function after stroke. However, more research on hyolingual biomechanics and control is needed for technologies interfacing the nervous system with the hyolingual apparatus to become a reality.

Architectural gear ratio depends on actuator spacing in a physical model of pennate muscle

Pennate muscles are defined by the architectural arrangement of their muscle fibers, which run at an angle to the primary axis of muscle shortening. Pennation angles can vary dynamically over the course of individual contractions, influencing the speed and distance of muscle shortening. Despite their relevance to muscle performance, the physical mechanisms that drive dynamic changes in pennation angle remain poorly understood. Muscle fibers bulge radially as they shorten, a consequence of maintaining a constant internal fluid volume, and we hypothesized that radial interactions between tightly packed muscle fibers are essential to dynamic pennation angle changes. To explore this, we built physical models of pennate muscles in which the radial distance between fiber-like actuators could be experimentally altered. Models were built from pennate arrays of McKibben actuators, a type of pneumatic actuator that forcefully shortens and bulges radially when inflated with compressed air. Consistent with past studies of biological muscle and engineered pennate actuators, we found that the magnitude of pennation angle change during contraction varied with load. Importantly, however, we found that pennation angle changes were also strongly influenced by the radial distance between neighboring McKibben actuators. Increasing the radial distance between neighboring actuators reduced pennation angle change during contraction and effectively eliminated variable responses to load. Radial interactions between muscle fibers are rarely considered in theoretical and experimental analyses of pennate muscle; however, these findings suggest that radial interactions between fibers drive pennation angle changes and influence pennate muscle performance. Our results provide insight into the fundamental mechanism underlying dynamic pennation angle changes in biological muscle and highlight design considerations that can inform the development of engineered pennate arrays.

Impairments in muscle shape changes affect metabolic demands during in-vivo contractions

The uncoupling behaviour between muscle belly and fascicle shortening velocity (i.e. belly gearing), affects mechanical output by allowing the muscle to circumvent the limits imposed by the fascicles' force-velocity relationship. However, little is known about the 'metabolic effect' of a decrease/increase in belly gearing. In this study, we manipulated the plantar flexor muscles' capacity to change in shape (and hence belly gearing) by using compressive multidirectional loads. Metabolic, kinetic, electromyography activity and ultrasound data (in soleus and gastrocnemius medialis) were recorded during cyclic fixed-end contractions of the plantar flexor muscles in three different conditions: no load, +5 kg and +10 kg of compression. No differences were observed in mechanical power and electrophysiological variables as a function of compression intensity, whereas metabolic power increased as a function of it. At each compression intensity, differences in efficiency were observed when calculated based on fascicle or muscle behaviour and significant positive correlations (R2 range: 0.7-0.8 and p > 0.001) were observed between delta efficiency (ΔEff: Effmus-Efffas) and belly gearing (Vmus/Vfas) or ΔV (Vmus-Vfas). Thus, changes in the muscles' capacity to change in shape (e.g. in muscle stiffness or owing to compressive garments) affect the metabolic demands and the efficiency of muscle contraction.

A review of the efforts to develop muscle and musculoskeletal models for biomechanics in the last 50 years

Both the Hill and the Huxley muscle models had already been described by the time the International Society of Biomechanics was founded 50 years ago, but had seen little use before the 1970s due to the lack of computing. As computers and computational methods became available in the 1970s, the field of musculoskeletal modeling developed and Hill type muscle models were adopted by biomechanists due to their relative computational simplicity as compared to Huxley type muscle models. Muscle forces computed by Hill type muscle models provide good agreement in conditions similar to the initial studies, i.e. for small muscles contracting under steady and controlled conditions. However, more recent validation studies have identified that Hill type muscle models are least accurate for natural in vivo locomotor behaviours at submaximal activations, fast speeds and for larger muscles, and thus need to be improved for their use in understanding human movements. Developments in muscle modelling have tackled these shortcomings. However, over the last 50 years musculoskeletal simulations have been largely based on traditional Hill type muscle models or even simplifications of this model that neglected the interaction of the muscle with a compliant tendon. The introduction of direct collocation in musculoskeletal simulations about 15 years ago along with further improvements in computational power and numerical methods enabled the use of more complex muscle models in simulations of whole-body movement. Whereas Hill type models are still the norm, we may finally be ready to adopt more complex muscle models into musculoskeletal simulations of human movement.

Force transmission and interactions between synergistic muscles

The classical view of muscles as independent motors has been challenged over the past decades. An alternative view has emerged in which muscles are not isolated but embedded in a three-dimensional connective tissue network that links them to adjacent muscles and other non-muscular structures in the body. Animal studies showing that the forces measured at the distal and proximal ends of a muscle are not equal have provided undisputable evidence that these connective tissue linkages are strong enough to serve as an extra pathway for muscular force transmission. In this historical review, we first introduce the terminology and anatomy related to these pathways of muscle force transmission and provide a definition for the term epimuscular force transmission. We then focus on important experimental evidence indicating mechanical interactions between synergistic muscles that may affect force transmission and/or influence the muscles' force generating capacity. We illustrate that there may exist different expressions of the highly relevant force-length properties depending on whether the force is measured at the proximal or distal tendon and depending on the dynamics of surrounding structures. Changes in length, activation level or disruption of the connective tissue of neighboring muscles, can affect how muscles interact and produce force on the skeleton. While most direct evidence is from animal experiments, studies on humans also suggest functional implications of the connective tissues surrounding muscles. These implications may explain how distant segments, which are not part of the same joint system, affect force generation at a given joint, and, in clinical conditions, explain observations from tendon transfer surgeries, where a muscle transferred to act as an antagonist continues to produce agonistic moments.

Fractional-order viscoelastic model of musculoskeletal tissues: correlation with fractals

Self-similar fractals are widely obtained from biomaterials within the human musculoskeletal system, and their viscoelastic behaviours can be described by fractional-order derivatives. However, existing viscoelastic models neglect the internal correlation between the fractal structure of biomaterials and their fractional-order temporal responses. We further expanded the fractal hyper-cell (FHC) viscoelasticity theory to investigate this spatio-temporal correlation. The FHC element was first compared with other material elements and spring–dashpot viscoelastic models, thereby highlighting its discrete and fractal nature. To demonstrate the utility of an FHC, tree-like, ladder-like and triangle-like FHCs were abstracted from human cartilage, tendons and muscle cross-sections, respectively. The duality and symmetry of the FHC element were further discussed, where operating the duality transformation generated new types of FHC elements, and the symmetry breaking of an FHC altered its fractional-order viscoelastic responses. Thus, the correlations between the staggering patterns of FHCs and their rheological power-law orders were established, and the viscoelastic behaviour of the multi-level FHC elements fitted well in stress relaxation experiments at both the macro- and nano-levels of the tendon hierarchy. The FHC element provides a theoretical basis for understanding the connections between structural degeneration of bio-tissues during ageing or disease and their functional changes.

Passive skeletal muscle can function as an osmotic engine

Muscles are composite structures. The protein filaments responsible for force production are bundled within fluid-filled cells, and these cells are wrapped in ordered sleeves of fibrous collagen. Recent models suggest that the mechanical interaction between the intracellular fluid and extracellular collagen is essential to force production in passive skeletal muscle, allowing the material stiffness of extracellular collagen to contribute to passive muscle force at physiologically relevant muscle lengths. Such models lead to the prediction, tested here, that expansion of the fluid compartment within muscles should drive forceful muscle shortening, resulting in the production of mechanical work unassociated with contractile activity. We tested this prediction by experimentally increasing the fluid volumes of isolated bullfrog semimembranosus muscles via osmotically hypotonic bathing solutions. Over time, passive muscles bathed in hypotonic solution widened by 16.44 ± 3.66% (mean ± s.d.) as they took on fluid. Concurrently, muscles shortened by 2.13 ± 0.75% along their line of action, displacing a force-regulated servomotor and doing measurable mechanical work. This behaviour contradicts the expectation for an isotropic biological tissue that would lengthen when internally pressurized, suggesting a functional mechanism analogous to that of engineered pneumatic actuators and highlighting the significance of three-dimensional force transmission in skeletal muscle.

The Effect of Multidirectional Loading on Contractions of the M. Medial Gastrocnemius

Research has shown that compression of muscle can lead to a change in muscle force. Most studies show compression to lead to a reduction in muscle force, although recent research has shown that increases are also possible. Based on methodological differences in the loading design between studies, it seems that muscle length and the direction of transverse loading influence the effect of muscle compression on force production. Thus, in our current study we implement these two factors to influence the effects of muscle loading. In contrast to long resting length of the medial gastrocnemius (MG) in most studies, we use a shorter MG resting length by having participant seated with their knees at a 90° angle. Where previous studies have used unidirectional loads to compress the MG, in this study we applied a multidirectional load using a sling setup. Multidirectional loading using a sling setup has been shown to cause muscle force reductions in previous research. As a result of our choices in experimental design we observed changes in the effects of muscle loading compared to previous research. In the present study we observed no changes in muscle force due to muscle loading. Muscle thickness and pennation angle showed minor but significant increases during contraction. However, no significant changes occurred between unloaded and loaded trials. Fascicle thickness and length showed different patterns of change compared to previous research. We show that muscle loading does not result in force reduction in all situations and is possibly linked to differences in muscle architecture and muscle length.

Effects of compression garments on surface EMG and physiological responses during and after distance running

BACKGROUND

The few previous studies that focused on the effects of compression garments (CG) on distance running performance have simultaneously measured electromyogram, physiological, and perceptual parameters. Therefore, this study investigated the effects of CG on muscle activation and median frequency during and after distance running, as well as blood-lactate concentration and rating of perceived exertion (RPE) during distance running.

METHODS

Eight healthy male recreational runners were recruited to randomly perform two 40 min treadmill running trials, one with CG, and the other with control garment made of normal cloth. The RPE and the surface electromyography (EMG) of 5 lower extremity muscles including gluteus maximus (GM), rectus femoris (RF), semitendinosus (ST), tibialis anterior (TA), and gastrocnemius (GAS) were measured during the running trial. The blood-lactate levels before and after the running trial were measured.

RESULTS

Wearing CG led to significant lower muscle activation (p  <  0.05) in the GM (decreased 7.40%-14.31%), RF (decreased 4.39%-4.76%), and ST (decreased 3.42%-7.20%) muscles; moreover, significant higher median frequency (p< 0.05) in the GM (increased 5.57%) and ST (increased 10.58%) muscles. Wearing CG did not alter the RPE values or the blood-lactate levels (p > 0.05).

CONCLUSION

Wearing CG was associated with significantly lower muscle activation and higher median frequency in the running-related key muscles during distance running. This finding suggested that wearing CG may improve muscle function, which might enhance running performance and prevent muscle fatigue.

The Multi-Scale, Three-Dimensional Nature of Skeletal Muscle Contraction

Muscle contraction is a three-dimensional process, as anyone who has observed a bulging muscle knows. Recent studies suggest that the three-dimensional nature of muscle contraction influences its mechanical output. Shape changes and radial forces appear to be important across scales of organization. Muscle architectural gearing is an emerging example of this process.

Predicting Perturbed Human Arm Movements in a Neuro-Musculoskeletal Model to Investigate the Muscular Force Response

Human movement is generated by a dynamic interplay between the nervous system, the biomechanical structures, and the environment. To investigate this interaction, we propose a neuro-musculoskeletal model of human goal-directed arm movements. Using this model, we simulated static perturbations of the inertia and damping properties of the arm, as well as dynamic torque perturbations for one-degree-of freedom movements around the elbow joint. The controller consists of a feed-forward motor command and feedback based on muscle fiber length and contraction velocity representing short-latency (25 ms) or long-latency (50 ms) stretch reflexes as the first neuronal responses elicited by an external perturbation. To determine the open-loop control signal, we parameterized the control signal resulting in a piecewise constant stimulation over time for each muscle. Interestingly, such an intermittent open-loop signal results in a smooth movement that is close to experimental observations. So, our model can generate the unperturbed point-to-point movement solely by the feed-forward command. The feedback only contributed to the stimulation in perturbed movements. We found that the relative contribution of this feedback is small compared to the feed-forward control and that the characteristics of the musculoskeletal system create an immediate and beneficial reaction to the investigated perturbations. The novelty of these findings is (1) the reproduction of static as well as dynamic perturbation experiments in one neuro-musculoskeletal model with only one set of basic parameters. This allows to investigate the model's neuro-muscular response to the perturbations that-at least to some degree-represent stereotypical interactions with the environment; (2) the demonstration that in feed-forward driven movements the muscle characteristics generate a mechanical response with zero-time delay which helps to compensate for the perturbations; (3) that this model provides enough biomechanical detail to allow for the prediction of internal forces, including joint loads and muscle-bone contact forces which are relevant in ergonomics and for the development of assistive devices but cannot be observed in experiments.

An experimental and computational investigation of the effects of volumetric boundary conditions on the compressive mechanics of passive skeletal muscle

Computational modeling, such as finite element analysis, is employed in a range of biomechanics specialties, including impact biomechanics and surgical planning. These models rely on accurate material properties for skeletal muscle, which comprises roughly 40% of the human body. Due to surrounding tissues, compressed skeletal muscle in vivo likely experiences a semi-confined state. Nearly all previous studies investigating passively compressed muscle at the tissue level have focused on muscle in unconfined compression. The goals of this study were to (1) examine the stiffness and time-dependent material properties of skeletal muscle subjected to both confined and unconfined compression (2) develop a model that captures passive muscle mechanics under both conditions and (3) determine the extent to which different assumptions of volumetric behavior affect model results. Muscle in confined compression exhibited stiffer behavior, agreeing with previous assumptions of near-incompressibility. Stress relaxation was found to be faster under unconfined compression, suggesting there may be different mechanisms that support load these two conditions. Finite element calibration was achieved through nonlinear optimization (normalized root mean square error <6%) and model validation was strong (normalized root mean square error <17%). Comparisons to commonly employed assumptions of bulk behavior showed that a simple one parameter approach does not accurately simulate confined compression. We thus recommend the use of a properly calibrated, nonlinear bulk constitutive model for modeling of skeletal muscle in vivo. Future work to determine mechanisms of passive muscle stiffness would enhance the efforts presented here.

Internal fluid pressure influences muscle contractile force

Fluid fills intracellular, extracellular, and capillary spaces within muscle. During normal physiological activity, intramuscular fluid pressures develop as muscle exerts a portion of its developed force internally. These pressures, typically ranging between 10 and 250 mmHg, are rarely considered in mechanical models of muscle but have the potential to affect performance by influencing force and work produced during contraction. Here, we test a model of muscle structure in which intramuscular pressure directly influences contractile force. Using a pneumatic cuff, we pressurize muscle midcontraction at 260 mmHg and report the effect on isometric force. Pressurization reduced isometric force at short muscle lengths (e.g., -11.87% of P₀ at 0.9 L₀), increased force at long lengths (e.g., +3.08% of P₀ at 1.25 L₀), but had no effect at intermediate muscle lengths ∼1.1-1.15 L₀ This variable response to pressurization was qualitatively mimicked by simple physical models of muscle morphology that displayed negative, positive, or neutral responses to pressurization depending on the orientation of reinforcing fibers representing extracellular matrix collagen. These findings show that pressurization can have immediate, significant effects on muscle contractile force and suggest that forces transmitted to the extracellular matrix via pressurized fluid may be important, but largely unacknowledged, determinants of muscle performance in vivo.

Characterization of Electromechanical Delay Based on a Biophysical Multi-Scale Skeletal Muscle Model

Skeletal muscles can be voluntary controlled by the somatic nervous system yielding an active contractile stress response. Thereby, the active muscle stresses are transmitted to the skeleton by a cascade of connective tissue and thus enable motion. In the context of joint perturbations as well as the assessment of the complexity of neural control, the initial phase of the muscle-tendon system's stress response has a particular importance and is analyzed by means of electromechanical delay (EMD). EMD is defined as the time lag between the stimulation of a muscle and a measurable change in force output. While EMD is believed to depend on multiple structures / phenomena, it is hard to separate their contributions experimentally. We employ a physiologically detailed, three-dimensional, multi-scale model of an idealized muscle-tendon system to analyze the influence of (i) muscle and tendon length, (ii) the material behavior of skeletal muscle and tendon tissue, (iii) the chemo-electro-mechanical behavior of the muscle fibers and (iv) neural control on EMD. Comparisons with experimental data show that simulated EMD values are within the physiological range, i.e., between 6.1 and 68.6 ms, and that the model is able to reproduce the characteristic EMD-stretch curve, yielding the minimum EMD at optimal length. Simulating consecutive recruitment of motor units increases EMD by more than 20 ms, indicating that during voluntary contractions neural control is the dominant factor determining EMD. In contrast, the muscle fiber action potential conduction velocity is found to influence EMD even of a 27 cm long muscle by not more than 3.7 ms. We further demonstrate that in conditions where only little pre-stretch is applied to a muscle-tendon system, the mechanical behavior of both muscle and tendon tissue considerably impacts EMD. Predicting EMD for different muscle and tendon lengths indicates that the anatomy of a specific muscle-tendon system is optimized for its function, i.e., shorter tendon lengths are beneficial to minimize the neural control effort for muscles primary acting as motor in concentric contractions.

Investigating the passive mechanical behaviour of skeletal muscle fibres: Micromechanical experiments and Bayesian hierarchical modelling

Characterisation of the skeletal muscle's passive properties is a challenging task since its structure is dominated by a highly complex and hierarchical arrangement of fibrous components at different scales. The present work focuses on the micromechanical characterisation of skeletal muscle fibres, which consist of myofibrils, by realising three different deformation states, namely, axial tension, axial compression, and transversal compression. To the best of the authors' knowledge, the present study provides a novel comprehensive data set representing of different deformation states. These data allow for a better understanding of muscle fibre load transfer mechanisms and can be used for meaningful modelling approaches. As the present study shows, axial tension and compression experiments reveal a strong tension-compression asymmetry at fibre level. In comparison to the tissue level, this asymmetric behaviour is more pronounced at the fibre scale, elucidating the load transfer mechanism in muscle tissue and aiding in the development of future modelling strategies. Further, a Bayesian hierarchical modelling approach was used to consider the experimental fluctuations in a parameter identification scheme, leading to more comprehensive parameter distributions that reflect the entire observed experimental uncertainty. STATEMENT OF SIGNIFICANCE: This article examines for the first time the mechanical properties of skeletal muscle fibres under axial tension, axial compression, and transversal compression, leading to a highly comprehensive data set. Moreover, a Bayesian hierarchical modelling concept is presented to identify model parameters in a broad way. The results of the deformation states allow a new and comprehensive understanding of muscle fibres' load transfer mechanisms; one example is the effect of tension-compression asymmetry. On the one hand, the results of this study contribute to the understanding of muscle mechanics and thus to the muscle's functional understanding during daily activity. On the other hand, they are relevant in the fields of skeletal muscle multiscale, constitutive modelling.

Significance of epimuscular myofascial force transmission under passive muscle conditions

In the past 20 yr, force transmission via connective tissue linkages at the muscle belly surface, called epimuscular myofascial force transmission, has been studied extensively. In this article, the effects of epimuscular linkages under passive muscle conditions are reviewed. Several animal studies that included direct (invasive) measurements of force transmission have shown that different connective tissue structures serve as an epimuscular pathway and that these tissues have sufficient stiffness, especially at supraphysiological muscle lengths and relative positions, to transmit substantial passive forces (up to 15% of active optimal force). Exact values of lumped tissue stiffness for different connective tissue structures have not yet been estimated. Experiments using various imaging techniques (ultrasound, MRI, shear wave elastography) have yielded some, but weak, evidence of epimuscular myofascial force transmission for passive muscles in humans. At this point, the functional consequences of epimuscular pathways for muscle and joint mechanics in the intact body are still unknown. Potentially, however, these pathways may affect sensory feedback and, thereby, neuromuscular control. In addition, altered epimuscular force transmission in pathological conditions may also contribute to changes in passive range of joint motion.

In Vivo Flattening of the Central Aponeurosis of the Rectus Femoris Due to Knee Extension Torque in Healthy Young and Elderly Individuals With Knee Osteoarthritis

ABSTRACT

We aimed to elucidate the relationship between active force production and the curvature of the central aponeurosis (CA) of the rectus femoris in young healthy participants as fundamental data and compare the muscle CA curvature before and after straight leg raising (SLR) training in participants with knee osteoarthritis (OA). Central aponeurosis curvature was determined during submaximal and maximal voluntary contractions (MVCs) using ultrasonography. Twenty-five young healthy female volunteers underwent ultrasonographic measurements under conditions of isometric MVC. They were divided into a flat shaped CA group (flat) and an incompletely flat shaped CA group (remnant). Central aponeurosis curvature was calculated as the ratio of CA height and length in the axial view. Central aponeurosis shape and muscular strength before and after muscle training were measured in 11 participants with knee OA. In the young healthy individuals, maximal voluntary torque and changes in CA curvature were significantly higher in the flat group than in the remnant group (2.15 Nm/kg and - 17.7% vs 1.75 Nm/kg and -9.8%, respectively; P = 0.005). The rate of change of the CA curvature during contraction was significantly correlated with maximal voluntary torque corrected for body mass (r = 0.512). The CA curvature progressively decreased as %MVC increased. In the OA group, CA curvature during MVC after SLR training was significantly lower than that before SLR training (3.2% vs 7.2%; P = 0.031). Central aponeurosis curvature was associated with muscle strength, and the results supported our hypothesis that geometric observation of CA changes during contractions may reflect muscle fiber function. We aim to develop a new ultrasonographic skeletal muscle evaluation method based on our present findings.

Intramuscular Pressure of Human Tibialis Anterior Muscle Reflects in vivo Muscular Activity

Intramuscular pressure (IMP) is the fluid hydrostatic pressure generated within a muscle and reflects the mechanical forces produced by a muscle. By providing accurate quantification of interstitial fluid pressure, the measurement of IMP may be useful to detect changes in skeletal muscle function not identified with established techniques. However, the relationship between IMP and muscle activity has never been studied in vivo in healthy human muscles. To determine if IMP is able to evaluate electromechanical performance of muscles in vivo, we tested the following hypotheses on the human tibialis anterior (TA) muscle: (i) IMP increases in proportion to muscle activity as measured by electrical [Compound Muscle Action Potential (CMAP)] and mechanical (ankle torque) responses to activation by nerve stimulation and (ii) the onset delay of IMP (IMPD) is shorter than the ankle torque electromechanical delay (EMD). Twelve healthy adults [six females; mean (SD) = 28.1 (5.0) years old] were recruited. Ankle torque, TA IMP, and CMAP responses were collected during maximal stimulation of the fibular nerve at different intensity levels of electrical stimulation, and at different frequencies of supramaximal stimulation, i.e., at 2, 5, 10, and 20 Hz. The IMP response at different stimulation intensities was correlated with the CMAP amplitude (r² = 0.94). The area of the IMP response at different stimulation intensities was also significantly correlated with the area of the CMAP (r² = 0.93). Increasing stimulation intensity resulted in an increase of the IMP response (P < 0.001). Increasing stimulation frequency caused torque (P < 0.001) as well as the IMP (P < 0.001) to increase. The ankle torque EMD [median (interquartile range) = 41.8 (14.4) ms] was later than the IMPD [33.0 (23.6) ms]. These findings support the hypotheses and suggest that IMP captures active mechanical properties of muscle in vivo and can be used to detect muscular changes due to drugs, diseases, or aging.

Passive and dynamic muscle architecture during transverse loading for gastrocnemius medialis in man

External forces from our environment impose transverse loads on our muscles. Studies in rats have shown that transverse loads result in a decrease in the longitudinal muscle force. Changes in muscle architecture during contraction may contribute to the observed force decrease. The aim of this study was to quantify changes in pennation angle, fascicle dimensions, and muscle thickness during contraction under external transverse load. Electrical stimuli were elicited to evoke maximal force twitches in the right calf muscles of humans. Trials were conducted with transverse loads of 2, 4.5, and 10 kg. An ultrasound probe was placed on the medial gastrocnemius in line with the transverse load to quantify muscle characteristics during muscle twitches. Maximum twitch force decreased with increased transverse muscle loading. The 2, 4.5, and 10 kg of transverse load showed a 9, 13, and 16% decrease in longitudinal force, respectively. Within the field of view of the ultrasound images, and thus directly beneath the external load, loading of the muscle resulted in a decrease in the muscle thickness and pennation angle, with higher loads causing greater decreases. During twitches the muscle transiently increased in thickness and pennation angle, as did fascicle thickness. Higher transverse loads showed a reduced increase in muscle thickness. Smaller increases in pennation angle and fascicle thickness strain also occurred with higher transverse loads. This study shows that increased transverse loading caused a decrease in ankle moment, muscle thickness, and pennation angle, as well as transverse deformation of the fascicles.

Structural Determinants of Muscle Gearing During Dynamic Contractions

In skeletal muscle, interactions between contractile and connective tissue elements at multiple scales result in emergent properties that determine mechanical performance. One of these phenomena is architectural gearing, which is quantified as the ratio of muscle velocity to muscle fiber velocity. Many pennate muscles operate with a gear ratio greater than one because muscles shorten through a combination of muscle fiber shortening and fiber rotation. Within a muscle, gearing is variable across contractions. During low force contractions, muscles operate at high gear while muscles operate at low gear during high force contractions. This variable gearing has a significant impact on muscle performance as muscle architectural changes favor muscle speed during fast contractions and muscle force during slow, high force contractions. We hypothesize that gearing in any given contraction is determined by the dynamic interaction of fiber-generated forces, fluid force transmission, and the elastic behavior of intramuscular connective tissues. Because muscle is isovolumetric, muscle fibers must bulge radially when they shorten. Radial bulging and fiber-generated forces off-axis from the muscle line of action exert forces that load connective tissues that ensheath fibers, fascicles, and the whole muscle. The way in which fluid pressures and fiber forces interact to load connective tissues in three-dimensions remains poorly understood because of the complex and multiscale nature of these interactions. Here we review evidence for variable gearing in pennate muscles, present a conceptual model that describes the fundamental interactions that determine gearing, and discuss where gaps remain in our understanding of the determinants and consequences of muscle shape change and variable gearing.

Impact of transversal calf muscle loading on plantarflexion

Muscle compression commonly occurs in daily life (for instance wearing backpacks or compression garments, and during sitting). However, the effects of the compression on contraction dynamics in humans are not well examined. The aim of the study was to quantify the alterations of contraction dynamics and muscle architecture in human muscle with external transverse loads. The posterior tibialis nerve of 29 subjects was stimulated to obtain the maximal double-twitch force of the gastrocnemius muscle with and without transverse compression that was generated using an indentor. The muscle architecture was determined by a sonographic probe that was embedded within the indentor. Five stimulations each were conducted at 5 conditions: (1) pretest (unloaded), (2) indentor loading with 2 kg, (3) 4.5 kg, (4) 10 kg, and (5) posttest (unloaded). Compared to the pretest maximal force decreased by 9%, 13% and 16% for 2 kg, 4.5 kg and 10 kg, respectively. The half-relaxation time increased with increased transverse load whereas the rate of force development decreased from pretest to 2 kg and from 4.5 kg to 10 kg. The lifting height of the indentor increased with transverse load from 2 kg to 4.5 kg but decreased from 4.5 kg to 10 kg. Increases in pennation during the twitches were reduced at the highest transverse load. The results demonstrate changes of the contraction dynamics due to transversal muscle loading. Those alterations are associated with the applied pressure, changes in muscle architecture and partitioning of muscle force in transversal and longitudinal direction.

Contribution of individual quadriceps muscles to knee joint mechanics

Many attempts have been made to determine the contribution of individual muscles in an agonistic group to the mechanics of joints. However, previous approaches had the limitations that muscles often could not be controlled in a precise manner, that individual muscles in an agonistic group could not be activated individually, and that individual muscle contributions could not be measured in an actively contracting agonistic group. Here, we introduce a surgical approach that allows for controlled activation of individual muscles of an agonistic group. The approach is illustrated for the vastus lateralis (VL), vastus medialis (VM), and rectus femoris (RF) of the rabbit quadriceps femoris group. We provide exemplar results for potential applications of the approach, such as measuring the pressure distribution in the patellofemoral joint, and the torque–angle relationship of VL, VM and RF when activated individually and when the three muscles are activated simultaneously.

Why do muscles lose torque potential when activated within their agonistic group?

Agonistic muscles lose approximately 20% of their individual torque generating capacity when activated with their agonistic muscles compared to when stimulated in isolation. In this study, we (i) tested if this loss in torque was accompanied by a corresponding loss in force, thereby testing the potential role of changes in moment arms between conditions; (ii) removed all inter-muscular connections between the quadriceps muscles, thus determining the potential role of inter-muscular force transmission; and (iii) systematically changed the inter-muscular pressure by performing experiments at different activation/force levels, thereby exploring the possible role of inter-muscular pressure in the loss of torque capacity with simultaneous muscle activation. Experiments were performed in a New Zealand White rabbit quadriceps model (n=5). Torque and force were measured during activation of femoral nerve branches that supply the individual quadriceps muscles while activating these branches simultaneously and in isolation. Regardless of joint angle and inter-muscular connections between muscles, the differences in torque values between the simultaneous and the isolated activation of the quadriceps muscles were also observed for the directly measured force values. Mean differences in simultaneous and isolated muscle activation remained similar between the intact and separated conditions: torque difference (21±5% of maximum isometric torque of intact condition [MICtorque], versus 19±6% MICtorque respectively) and for force (18±3% MICforce versus 19±7% MICforce respectively). The absolute torque loss was independent of the force, and thus presumably the inter-muscular pressures. Based on these results, we conclude that neither moment arm, inter-muscular pressure nor inter-muscular force transmission seems to be the primary cause for the torque deficit observed during simultaneous compared to isolated muscle activation. The mechanisms underlying loss of force capacity during agonistic muscle contraction remain unknown.

Impact of Multidirectional Transverse Calf Muscle Loading on Calf Muscle Force in Young Adults

It has been demonstrated that unidirectional transversal muscle loading induced by a plunger influences muscle shape and reduces muscle force. The interaction between muscle and transversal forces may depend on specific neuromuscular properties that change during a lifetime. Compression garments, applying forces from all directions in the transverse plane, are widely used in sports for example to improve performance. Differences in the loading direction (unidirectional vs. multidirectional) may have an impact on force generating capacity of muscle and, thus, on muscle performance. The aim of this study was to examine the effect of multidirectional transversal loads, using a sling looped around the calf, on the isometric force during plantarflexions. Young male adults (25.7 ± 1.5 years, n = 15) were placed in a prone position in a calf press apparatus. The posterior tibial nerve was stimulated to obtain the maximal double-twitch force of the calf muscles with (59.4 and 108.4 N) and without multidirectional transverse load. Compared to the unloaded condition, the rate of force development (RFD) was reduced by 5.0 ± 8.1% (p = 0.048) and 6.9 ± 10.7% (p = 0.008) for the 59.4 and 108.4 N load, respectively. No significant reduction (3.2 ± 4.8%, p = 0.141) in maximum muscle force (Fm ) was found for the lower load (59.4 N), but application of the higher load (108.4 N) resulted in a significant reduction of Fm by 4.8 ± 7.0% (p = 0.008). Mean pressures induced in this study (14.3 and 26.3 mm Hg corresponding to the 59.4 and 108.4 N loads, respectively) are within the pressure range reported for compression garments. Taking the results of the present study into account, a reduction in maximum muscle force would be expected for compression garments with pressures ≥26.3 mm Hg. However, it should be noted that the loading condition (sling vs. compression garment) differs and that compression garments may influence other mechanisms contributing to force generation. For example, wearing compression garments may enhance sport performance by enhanced proprioception and reduced muscle oscillation. Thus, superposition of several effects should be considered when analyzing the impact of compression garments on more complex sport performance.

Transverse anisotropy in the deformation of the muscle during dynamic contractions

When pennate muscle fibres shorten, the transverse deformation of fibres results in an increase in pennation angle of fascicles (bundles of fibres) and transverse deformation of muscle belly. Transverse shape changes of a muscle can influence force generation. Recent modelling studies predicted asymmetrical transverse deformations in the muscle fascicles in the gastrocnemii. However, these predictions have not been tested experimentally. As muscle is a 3D entity, it is important to explore the structural changes in a 3D perspective to enhance our understanding of the underlying structural mechanisms that have functional implications. The medial and lateral gastrocnemius muscles from 12 subjects were imaged during plantarflexion movements on a dynamometer. The muscle belly was simultaneously scanned from two orthogonal directions using two ultrasound probes. Fascicle deformations were measured from the two orthogonal ultrasound scans to provide 3D information of muscle geometry. Whilst transverse deformations in the medial gastrocnemius were similar from the two directions, the data for the lateral gastrocnemius confirm that transverse anisotropy can occur in the muscle fascicles. As the lateral gastrocnemius fascicle length shortened, the pennation angle increased and the fascicles bulged transversally in one direction (closest to the typical 2D scanning plane) while thinning in the other orthogonal direction. We suggest that the transverse deformation of the muscle fascicles depends on the stiffness of the aponeuroses, properties of connective tissue structures surrounding muscle, and compressive forces both internal and external to the muscle. These results highlight that muscle fascicles do not bulge uniformly and the implications for this behaviour on muscle function remain largely unexplored.

Size, History-Dependent, Activation and Three-Dimensional Effects on the Work and Power Produced During Cyclic Muscle Contractions

Muscles undergo cycles of length change and force development during locomotion, and these contribute to their work and power production to drive body motion. Muscle fibers are typically considered to be linear actuators whose stress depends on their length, velocity, and activation state, and whose properties can be scaled up to explain the function of whole muscles. However, experimental and modeling studies have shown that a muscle's stress additionally depends on inactive and passive tissues within the muscle, the muscle's size, and its previous contraction history. These effects have not been tested under common sets of contraction conditions, especially the cyclic contractions that are typical of locomotion. Here we evaluate the relative effects of size, history-dependent, activation and three-dimensional effects on the work and power produced during cyclic contractions of muscle models. Simulations of muscle contraction were optimized to generate high power outputs: this resulted in the muscle models being largely active during shortening, and inactive during lengthening. As such, the history-dependent effects were dominated by force depression during simulated active shortening rather than force enhancement during active stretch. Internal work must be done to deform the muscle tissue, and to accelerate the internal muscle mass, resulting in reduced power and work that can be done on an external load. The effect of the muscle mass affects the scaling of muscle properties, with the inertial costs of contraction being relatively greater at larger sizes and lower activation levels.

A hill-type muscle model expansion accounting for effects of varying transverse muscle load

Recent studies demonstrated that uniaxial transverse loading (F_G) of a rat gastrocnemius medialis muscle resulted in a considerable reduction of maximum isometric muscle force (ΔF_im). A hill-type muscle model assuming an identical gearing G between both ΔF_im and F_G as well as lifting height of the load (Δh) and longitudinal muscle shortening (Δl_CC) reproduced experimental data for a single load. Here we tested if this model is able to reproduce experimental changes in ΔF_im and Δh for increasing transverse loads (0.64 N, 1.13 N, 1.62 N, 2.11 N, 2.60 N). Three different gearing ratios were tested: (I) constant G_c representing the idea of a muscle specific gearing parameter (e.g. predefined by the muscle geometry), (II) G_exp determined in experiments with varying transverse load, and (III) G_f that reproduced experimental ΔF_im for each transverse load. Simulations using G_c overestimated ΔF_im (up to 59%) and Δh (up to 136%) for increasing load. Although the model assumption (equal G for forces and length changes) held for the three lower loads using G_exp and G_f, simulations resulted in underestimation of ΔF_im by 38% and overestimation of Δh by 58% for the largest load, respectively. To simultaneously reproduce experimental ΔF_im and Δh for the two larger loads, it was necessary to reduce F_im by 1.9% and 4.6%, respectively. The model seems applicable to account for effects of muscle deformation within a range of transverse loading when using a linear load-dependent function for G.

Resistance to radial expansion limits muscle strain and work

The collagenous extracellular matrix (ECM) of skeletal muscle functions to transmit force, protect sensitive structures, and generate passive tension to resist stretch. The mechanical properties of the ECM change with age, atrophy, and neuromuscular pathologies, resulting in an increase in the relative amount of collagen and an increase in stiffness. Although numerous studies have focused on the effect of muscle fibrosis on passive muscle stiffness, few have examined how these structural changes may compromise contractile performance. Here we combine a mathematical model and experimental manipulations to examine how changes in the mechanical properties of the ECM constrain the ability of muscle fibers and fascicles to radially expand and how such a constraint may limit active muscle shortening. We model the mechanical interaction between a contracting muscle and the ECM using a constant volume, pressurized, fiber-wound cylinder. Our model shows that as the proportion of a muscle cross section made up of ECM increases, the muscle's ability to expand radially is compromised, which in turn restricts muscle shortening. In our experiments, we use a physical constraint placed around the muscle to restrict radial expansion during a contraction. Our experimental results are consistent with model predictions and show that muscles restricted from radial expansion undergo less shortening and generate less mechanical work under identical loads and stimulation conditions. This work highlights the intimate mechanical interaction between contractile and connective tissue structures within skeletal muscle and shows how a deviation from a healthy, well-tuned relationship can compromise performance.

Intermuscular pressure between synergistic muscles correlates with muscle force

The purpose of the study was to examine the relationship between muscle force generated during isometric contractions (i.e. at a constant muscle–tendon unit length) and the intermuscular (between adjacent muscles) pressure in synergistic muscles. Therefore, the pressure at the contact area of the gastrocnemius and plantaris muscle was measured synchronously to the force of the whole calf musculature in the rabbit species Oryctolagus cuniculus. Similar results were obtained when using a conductive pressure sensor, or a fibre-optic pressure transducer connected to a water-filled balloon. Both methods revealed a strong linear relationship between force and pressure in the ascending limb of the force-length relationship. The shape of the measured force–time and pressure–time traces was almost identical for each contraction (r=0.97). Intermuscular pressure ranged between 100 and 700 mbar (70,000 Pa) for forces up to 287 N. These pressures are similar to previous (intramuscular) recordings within skeletal muscles of different vertebrate species. Furthermore, our results suggest that the rise in intermuscular pressure during contraction may reduce the force production in muscle packages (compartments).

The pH heterogeneity in human calf muscle during neuromuscular electrical stimulation

PURPOSE

The aim of the study was to examine pH heterogeneity during fatigue induced by neuromuscular electrical stimulation (NMES) using phosphorus magnetic resonance spectroscopy (³¹ P-MRS). It is hypothesized that three pH components would occur in the ³¹ P-MRS during fatigue, representing three fiber types.

METHODS

The medial gastrocnemius of eight subjects was stimulated within a 3-Tesla whole body MRI scanner. The maximal force during stimulation (F_stim ) was examined by a pressure sensor. Phosphocreatine (PCr), adenosintriphosphate, inorganic phosphate (Pi), and the corresponding pH were estimated by a nonvolume-selective ³¹ P-MRS using a small loop coil at rest and during fatigue.

RESULTS

During fatigue, F_stim and PCr decreased to 27% and 33% of their initial levels, respectively. In all cases, the Pi peak increased when NMES was started and split into three different peaks. Based on the single Pi peaks during fatigue, an alkaline (6.76 ± 0.08), a medium (6.40 ± 0.06), and an acidic (6.09 ± 0.05) pH component were observed compared to the pH (7.02 ± 0.02) at rest.

CONCLUSION

It is suggested that NMES is able to induce pH heterogeneity in the medial gastrocnemius, and that the single Pi peaks represent the different muscle fiber types of the skeletal muscle. Magn Reson Med 77:2097-2106, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

A lumped stiffness model of intermuscular and extramuscular myofascial pathways of force transmission

Mechanical behavior of skeletal muscles is commonly modeled under the assumption of mechanical independence between individual muscles within a muscle group. Epimuscular myofascial force transmission via the connective tissue network surrounding a muscle challenges this assumption as it alters the force distributed to the tendons of individual muscles. This study aimed to derive a lumped estimate of stiffness of the intermuscular and extramuscular connective tissues and to assess changes in such stiffness in response to a manipulation of the interface between adjacent muscles. Based on in situ measurements of force transmission in the rat plantar flexors, before and after resection of their connective tissue network, a nonlinear estimate of epimuscular myofascial stiffness was quantified and included in a multi-muscle model with lumped parameters which allows for force transmission depending on the relative position between the muscles in the group. Such stiffness estimate was assessed for a group with normal intermuscular connective tissues and for a group with increased connectivity, mimicking scar tissue development. The model was able to successfully predict the amount of epimuscular force transmission for different experimental conditions than those used to obtain the model parameters. The proposed nonlinear stiffness estimates of epimuscular pathways could be integrated in larger musculoskeletal models, to provide more accurate predictions of force when effects of mechanical interaction or altered epimuscular connections, e.g. after surgery or injury, are substantial.

Force reduction induced by unidirectional transversal muscle loading is independent of local pressure

Transversal unidirectional compression applied to muscles via external loading affects muscle contraction dynamics in the longitudinal direction. A recent study reported decreasing longitudinal muscle forces with increasing transversal load applied with a constant contact area (i.e., leading to a simultaneous increase in local pressure). To shed light on these results, we examine whether the decrease in longitudinal force depends on the load, the local pressure, or both. To this end, we perform isometric experiments on rat M. gastrocnemius medialis without and with transversal loading (i) changing the local pressure from 1.1-3.2Ncm(-2) (n=9) at a constant transversal load (1.62N) and (ii) increasing the transversal load (1.15-3.45N) at a constant local pressure of 2.3Ncm(-2) (n=7). While we did not note changes in the decrease in longitudinal muscle force in the first experiment, the second experiment resulted in an almost-linear reduction of longitudinal force between 7.5±0.6% and 14.1±1.7%. We conclude that the observed longitudinal force reduction is not induced by local effects such as malfunction of single muscle compartments, but that similar internal stress conditions and myofilament configurations occur when the local pressure changes given a constant load. The decreased longitudinal force may be explained by increased internal pressure and a deformed myofilament lattice that is likely associated with the decomposition of cross-bridge forces on the one hand and the inhibition of cross-bridges on the other hand.

Three-Dimensional Muscle Architecture and Comprehensive Dynamic Properties of Rabbit Gastrocnemius, Plantaris and Soleus: Input for Simulation Studies

The vastly increasing number of neuro-muscular simulation studies (with increasing numbers of muscles used per simulation) is in sharp contrast to a narrow database of necessary muscle parameters. Simulation results depend heavily on rough parameter estimates often obtained by scaling of one muscle parameter set. However, in vivo muscles differ in their individual properties and architecture. Here we provide a comprehensive dataset of dynamic (n = 6 per muscle) and geometric (three-dimensional architecture, n = 3 per muscle) muscle properties of the rabbit calf muscles gastrocnemius, plantaris, and soleus. For completeness we provide the dynamic muscle properties for further important shank muscles (flexor digitorum longus, extensor digitorum longus, and tibialis anterior; n = 1 per muscle). Maximum shortening velocity (normalized to optimal fiber length) of the gastrocnemius is about twice that of soleus, while plantaris showed an intermediate value. The force-velocity relation is similar for gastrocnemius and plantaris but is much more bent for the soleus. Although the muscles vary greatly in their three-dimensional architecture their mean pennation angle and normalized force-length relationships are almost similar. Forces of the muscles were enhanced in the isometric phase following stretching and were depressed following shortening compared to the corresponding isometric forces. While the enhancement was independent of the ramp velocity, the depression was inversely related to the ramp velocity. The lowest effect strength for soleus supports the idea that these effects adapt to muscle function. The careful acquisition of typical dynamical parameters (e.g. force-length and force-velocity relations, force elongation relations of passive components), enhancement and depression effects, and 3D muscle architecture of calf muscles provides valuable comprehensive datasets for e.g. simulations with neuro-muscular models, development of more realistic muscle models, or simulation of muscle packages.