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

Impact of lengthening velocity on the generation of eccentric force by slow-twitch muscle fibers in long stretches

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.

Unique enlargement of human soleus muscle for bipedalism at the expense of the ease of leg swing

Humans exhibit unique skeletal muscle morphologies that are known to matter in upright bipedalism. However, their relevance to the ease of leg swing, which limits locomotion performance, remains unclear. Here, we aimed to examine muscle mass distribution within the human leg and the effect of each muscle on the ease of leg swing. We calculated the mass, center of mass position, and moment of inertia around the hip extension-flexion axis for all leg muscles by using a publicly available dataset of the 3D reconstruction of the musculoskeletal components in human male and female legs. The leg muscles showed a top-heavy-bottom-light tapering trend; muscles far from the hip joint tended to have smaller masses. Interestingly, however, the soleus exhibited sizable mass for its location. Consequently, the moment of inertia of the soleus was exceptionally greatest, accounting for approximately one-quarter of that of all muscles. These results indicate that compared to the other muscles the soleus muscle has a much larger effect on the leg moment of inertia and uniquely makes humans difficult to swing the leg, although the leg muscles basically show the top-heavy bottom-light tapering trend favoring the leg swing. Our findings highlight a novel functional consequence of human body evolution, suggesting that muscular enlargement for postural stability and endurance capacity has compromised the locomotion speed during the adaptation to bipedalism.

3D ultrasound-based determination of skeletal muscle fascicle orientations

Architectural parameters of skeletal muscle such as pennation angle provide valuable information on muscle function, since they can be related to the muscle force generating capacity, fiber packing, and contraction velocity. In this paper, we introduce a 3D ultrasound-based workflow for determining 3D fascicle orientations of skeletal muscles. We used a custom-designed automated motor driven 3D ultrasound scanning system for obtaining 3D ultrasound images. From these, we applied a custom-developed multiscale-vessel enhancement filter-based fascicle detection algorithm and determined muscle volume and pennation angle. We conducted trials on a phantom and on the human tibialis anterior (TA) muscle of 10 healthy subjects in plantarflexion (157 ± 7∘ ), neutral position (109 ± 7∘ , corresponding to neutral standing), and one resting position in between (145 ± 6∘ ). The results of the phantom trials showed a high accuracy with a mean absolute error of 0.92 ± 0.59∘ . TA pennation angles were significantly different between all positions for the deep muscle compartment; for the superficial compartment, angles are significantly increased for neutral position compared to plantarflexion and resting position. Pennation angles were also significantly different between superficial and deep compartment. The results of constant muscle volumes across the 3 ankle joint angles indicate the suitability of the method for capturing 3D muscle geometry. Absolute pennation angles in our study were slightly lower than recent literature. Decreased pennation angles during plantarflexion are consistent with previous studies. The presented method demonstrates the possibility of determining 3D fascicle orientations of the TA muscle in vivo.

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.

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.

On a three-dimensional model for the description of the passive characteristics of skeletal muscle tissue

In this work, a three-dimensional model was developed to describe the passive mechanical behaviour of anisotropic skeletal muscle tissue. To validate the model, orientation-dependent axial ([Formula: see text], [Formula: see text], [Formula: see text]) and semi-confined compression experiments (mode I, II, III) were performed on soleus muscle tissue from rabbits. In the latter experiments, specimen deformation is prescribed in the loading direction and prevented in an additional spatial direction, fibre compression at [Formula: see text] (mode I), fibre elongation at [Formula: see text] (mode II) and a neutral state of the fibres at [Formula: see text] where their length is kept constant (mode III). Overall, the model can adequately describe the mechanical behaviour with a relatively small number of model parameters. The stiffest tissue response during orientation-dependent axial compression ([Formula: see text] kPa) occurs when the fibres are oriented perpendicular to the loading direction ([Formula: see text]) and are thus stretched during loading. Semi-confined compression experiments yielded the stiffest tissue ([Formula: see text] kPa) in mode II when the muscle fibres are stretched. The extensive data set collected in this study allows to study the different error measures depending on the deformation state or the combination of deformation states.

Eccentric muscle contractions: from single muscle fibre to whole muscle mechanics

Eccentric muscle loading encompasses several unique features compared to other types of contractions. These features include increased force, work, and performance at decreased oxygen consumption, reduced metabolic cost, improved energy efficiency, as well as decreased muscle activity. This review summarises explanatory approaches to long-standing questions in terms of muscular contraction dynamics and molecular and cellular mechanisms underlying eccentric muscle loading. Moreover, this article intends to underscore the functional link between sarcomeric components, emphasising the fundamental role of titin in skeletal muscle. The giant filament titin reveals versatile functions ranging from sarcomere organisation and maintenance, providing passive tension and elasticity, and operates as a mechanosensory and signalling platform. Structurally, titin consists of a viscoelastic spring segment that allows activation-dependent coupling to actin. This titin-actin interaction can explain linear force increases in active lengthening experiments in biological systems. A three-filament model of skeletal muscle force production (mediated by titin) is supposed to overcome significant deviations between experimental observations and predictions by the classic sliding-filament and cross-bridge theories. Taken together, this review intends to contribute to a more detailed understanding of overall muscle behaviour and force generation-from a microscopic sarcomere level to a macroscopic multi-joint muscle level-impacting muscle modelling, the understanding of muscle function, and disease.

Comparison of the soleus and plantaris muscles in humans and other primates: Macroscopic neuromuscular anatomy and evolutionary significance

In humans, the soleus is more developed compared to other primates and has a unique architecture composed of anterior bipennate and posterior unipennate parts, which are innervated by different nerve branches. The anterior part of the human soleus was proposed to be important for bipedalism, however, the phylogenetic process resulting in its acquisition remains unclear. Providing insights into this process, the anterior part of the soleus was suggested to be closely related to the plantaris based on the branching pattern of their nerve fascicles. To reveal the phylogeny of the soleus and plantaris in primates, the innervation patterns of the posterior crural muscles were compared among a wide range of species. From their branching pattern, posterior crural muscles could be classified into superficial and deep muscle groups. The anterior part of the soleus and plantaris both belonged to the deep muscle group. In all the examined specimens of ring-tailed lemurs and chimpanzees, as well as in one out of two specimens of siamang, the nerve branches corresponding to those innervating the anterior part of the human soleus were found. The muscular branches innervating the anterior part of the soleus and plantaris formed a common trunk or were connected in all the specimens. These results indicate that the anterior part of the soleus is closely related to the plantaris across different species of primates. In turn, this suggests that the anterior part of the soleus is maintained among primates, and especially in humans, where it develops as the characteristic bipennate structure.

Common synaptic input between motor units from the lateral and medial posterior soleus compartments does not differ from that within each compartment

The human soleus muscle is anatomically divided into four separate anatomical compartments. The functional role of this compartmentalization remains unclear. Here, we tested the hypothesis that the common synaptic input to motor units between the medial and lateral posterior compartments is less than within each compartment. Fourteen male participants performed three different heel-raise tasks that were considered to place a different mechanical demand on the medial and lateral soleus compartments. High-density electromyography (EMG) signals from the medial and lateral soleus compartments and the medial gastrocnemius of the right leg were decomposed into individual motor unit spike trains. The coherence between cumulative spike trains of the motor units was estimated. The coherence analysis was also repeated for motor units that were matched across all three tasks. Furthermore, we calculated the ratio of significant correlations between the spike trains of pairs of motor units. We observed that the coherence between motor units of the two soleus compartments was similar as the coherence between motor units within each compartment, regardless of the task. The correlation analysis performed on pairs of motor units confirmed these results. We conclude that the level of common synaptic input between the motor units innervating the medial and lateral posterior soleus compartment is not different than the common synaptic input between motor units innervating each of these compartments, which contrasts with findings from previous studies on finger muscles. This suggests that there is no independent neural control for the individual posterior soleus compartments.NEW & NOTEWORTHY The human soleus muscle is anatomically subdivided into four compartments. The functional role for this compartmentalization remains unknown. Here, we showed that, contrary to previous findings in finger muscles, the common synaptic input between motor units innervating the medial and lateral posterior soleus compartment was similar as that between motor units within the individual compartments. We suggest that the contradictory findings with other compartmentalized muscles may be explained by differences in muscle-tendon anatomy and function.

From fibre to function: are we accurately representing muscle architecture and performance?

The size and arrangement of fibres play a determinate role in the kinetic and energetic performance of muscles. Extrapolations between fibre architecture and performance underpin our understanding of how muscles function and how they are adapted to power specific motions within and across species. Here we provide a synopsis of how this 'fibre to function' paradigm has been applied to understand muscle design, performance and adaptation in animals. Our review highlights the widespread application of the fibre to function paradigm across a diverse breadth of biological disciplines but also reveals a potential and highly prevalent limitation running through past studies. Specifically, we find that quantification of muscle architectural properties is almost universally based on an extremely small number of fibre measurements. Despite the volume of research into muscle properties, across a diverse breadth of research disciplines, the fundamental assumption that a small proportion of fibre measurements can accurately represent the architectural properties of a muscle has never been quantitatively tested. Subsequently, we use a combination of medical imaging, statistical analysis, and physics-based computer simulation to address this issue for the first time. By combining diffusion tensor imaging (DTI) and deterministic fibre tractography we generated a large number of fibre measurements (>3000) rapidly for individual human lower limb muscles. Through statistical subsampling simulations of these measurements, we demonstrate that analysing a small number of fibres (n < 25) typically used in previous studies may lead to extremely large errors in the characterisation of overall muscle architectural properties such as mean fibre length and physiological cross-sectional area. Through dynamic musculoskeletal simulations of human walking and jumping, we demonstrate that recovered errors in fibre architecture characterisation have significant implications for quantitative predictions of in-vivo dynamics and muscle fibre function within a species. Furthermore, by applying data-subsampling simulations to comparisons of muscle function in humans and chimpanzees, we demonstrate that error magnitudes significantly impact both qualitative and quantitative assessment of muscle specialisation, potentially generating highly erroneous conclusions about the absolute and relative adaption of muscles across species and evolutionary transitions. Our findings have profound implications for how a broad diversity of research fields quantify muscle architecture and interpret muscle function.

How velocity impacts eccentric force generation of fully activated skinned skeletal muscle fibers in long stretches

Eccentric muscle contractions are fundamental to everyday life. They occur markedly in jumping, running, and accidents. Following an initial force rise, stretching of a fully activated muscle can result in a phase of decreasing force ('Give') followed by force redevelopment. However, how the stretch velocity affects 'Give' and force redevelopment remains largely unknown. We investigated the force produced by fully activated single skinned fibers of rat extensor digitorum longus muscles during long stretches. Fibers were pulled from length .85 to 1.3 optimal fiber length at a rate of 1, 10 and 100% of the estimated maximum shortening velocity. 'Give' was absent in slow stretches. Medium and fast stretches yielded a clear 'Give'. After the initial force peak, forces decreased by 11.2% and 27.8% relative to the initial peak force before rising again. During the last half of the stretch (from 1.07 to 1.3 optimal fiber length, which is within the range of the expected descending limb of the force-length relationship), the linear force slope tripled from slow to medium stretch and increased further by 60% from medium to fast stretch. These results are compatible with forcible cross-bridge detachment and re-development of a cross-bridge distribution, and a viscoelastic titin contribution to fiber force. Accounting for these results can improve muscle models and predictions of multi-body simulations.

Age-dependent mechanical and microstructural properties of the rabbit soleus muscle

During growth there are serious changes in the skeletal muscles to compensate for the changed requirements in terms of body weight and size. In this study, the age-dependent (between 21 and 100 days) mechanical and microstructural properties of rabbit soleus muscle tissue were investigated. For this purpose, morphological properties (animal mass, soleus muscle mass, tibial length) were measured at 5 different times during aging. On the other hand, fibre orientation-dependent axial and semi-confined compression experiments were realised. In addition, the essential components (muscle fibres, extracellular matrix, remaining components), dominating the microstructure of muscle tissue, were analysed. While the mechanical results show hardly any age-dependent differences, the morphological and microstructural results show clear age-dependent differences. All morphological parameters increase significantly (animal mass by 839.2%, muscle mass 1050.6%, tibial length 233.6%). In contrast, microstructural parameters change differently. The percentage of fibres (divided into slow-twitch (ST) and fast-twitch (FT) fibres) increases significantly (137.6%), while the proportions of the extracellular matrix and the remaining components (48.2% and 46.1%) decrease. At the same time, the cross-sectional area of the fibres increases significantly (697.9%). The totality of this age-dependent information provides a deeper understanding of age-related changes in muscle structure and function and may contribute to successful development and validation of growth models in the future. STATEMENT OF SIGNIFICANCE: This article reports the first comprehensive data set on age-dependent morphological (animal mass, soleus muscle mass, tibial length), mechanical (axial and semi-confined compression), and microstructural (muscle fibres, extracellular matrix, remaining components) properties of the rabbit soleus muscle. On the one hand, the results of this study contribute to the understanding of muscle mechanics and thus to understanding of load transfer mechanisms inside the muscle tissue during growth. On the other hand, these results are relevant to the fields of constitutive formulation of age-dependent muscle tissue.

Architectural model for muscle growth during maturation

Muscle architecture, which includes parameters like fascicle length, pennation angle, and physiological cross-sectional area, strongly influences skeletal muscles' mechanical properties. During maturation, the muscle architecture has to adapt to a growing organism. This study aimed to develop an architectural model capable of predicting the complete 3D fascicle architecture for primarily unipennate muscles of an arbitrary age, based on fascicle data for an initial age. For model development, we collected novel data on 3D muscle architecture of the rabbit (Oryctolagus cuniculus) M. plantaris of eight animals ranging in age from 29 to 106 days. Experimental results show that plantaris muscle belly length increases by 73%, whereas mean fascicle length and mean pennation angle increases by 39 and 14%, respectively. Those changes were incorporated into the model. In addition to the data collected for M. plantaris the predictions of the model were compared to existing literature data of rabbit M. soleus and M. gastrocnemius medialis. With an error of −1.0 ± 8.6% for relative differences in aponeurosis length, aponeurosis width, muscle height, and muscle mass, the model delivered good results matching interindividual differences. For future studies, the model could be utilized to generate realistic architectural data sets for simulation studies.

Influence of layer separation on the determination of stomach smooth muscle properties

Uniaxial tensile experiments are a standard method to determine the contractile properties of smooth muscles. Smooth muscle strips from organs of the urogenital and gastrointestinal tract contain multiple muscle layers with different muscle fiber orientations, which are frequently not separated for the experiments. During strip activation, these muscle fibers contract in deviant orientations from the force-measuring axis, affecting the biomechanical characteristics of the tissue strips. This study aimed to investigate the influence of muscle layer separation on the determination of smooth muscle properties. Smooth muscle strips, consisting of longitudinal and circumferential muscle layers (whole-muscle strips [WMS]), and smooth muscle strips, consisting of only the circumferential muscle layer (separated layer strips [SLS]), have been prepared from the fundus of the porcine stomach. Strips were mounted with muscle fibers of the circumferential layer inline with the force-measuring axis of the uniaxial testing setup. The force-length (FLR) and force-velocity relationships (FVR) were determined through a series of isometric and isotonic contractions, respectively. Muscle layer separation revealed no changes in the FLR. However, the SLS exhibited a higher maximal shortening velocity and a lower curvature factor than WMS. During WMS activation, the transversally oriented muscle fibers of the longitudinal layer shortened, resulting in a narrowing of this layer. Expecting volume constancy of muscle tissue, this narrowing leads to a lengthening of the longitudinal layer, which counteracted the shortening of the circumferential layer during isotonic contractions. Consequently, the shortening velocities of the WMS were decreased significantly. This effect was stronger at high shortening velocities.

Direct measurement of the direction-dependent mechanical behaviour of skeletal muscle extracellular matrix

This paper reports the first comprehensive data set on the anisotropic mechanical properties of isolated endo- and perimysial extracellular matrix of skeletal muscle, and presents the corresponding protocols for preparing and testing the samples. In particular, decellularisation of porcine skeletal muscle is achieved with caustic soda solution, and mechanical parameters are defined based on compressive and tensile testing in order to identify the optimal treatment time such that muscle fibres are dissolved whereas the extracellular matrix remains largely intact and mechanically functional. At around 18 h, a time window was found and confirmed by histology, in which axial tensile experiments were performed to characterise the direction-dependent mechanical response of the extracellular matrix samples, and the effect of lateral pre-compression was studied. The typical, large variability in the experimental stress response could be largely reduced by varying a single scalar factor, which was attributed to the variation of the fraction of extracellular matrix within the tissue. While experimental results on the mechanical properties of intact muscle tissue and single muscle fibres are increasingly available in literature, there is a lack of information on the properties of the collagenous components of skeletal muscle. The present work aims at closing this gap and thus contributes to an improved understanding of the mechanics of skeletal muscle tissue and provides a missing piece of information for the development of corresponding constitutive and computational models.

Cross-Bridges and Sarcomeric Non-cross-bridge Structures Contribute to Increased Work in Stretch-Shortening Cycles

Stretch-shortening cycles (SSCs) refer to the muscle action when an active muscle stretch is immediately followed by active muscle shortening. This combination of eccentric and concentric contractions is the most important type of daily muscle action and plays a significant role in natural locomotion such as walking, running or jumping. SSCs are used in human and animal movements especially when a high movement speed or economy is required. A key feature of SSCs is the increase in muscular force and work during the concentric phase of a SSC by more than 50% compared with concentric muscle actions without prior stretch (SSC-effect). This improved muscle capability is related to various mechanisms, including pre-activation, stretch-reflex responses and elastic recoil from serial elastic tissues. Moreover, it is assumed that a significant contribution to enhanced muscle capability lies in the sarcomeres itself. Thus, we investigated the force output and work produced by single skinned fibers of rat soleus muscles during and after ramp contractions at a constant velocity. Shortening, lengthening, and SSCs were performed under physiological boundary conditions with 85% of the maximum shortening velocity and stretch-shortening magnitudes of 18% of the optimum muscle length. The different contributions of cross-bridge (XB) and non-cross-bridge (non-XB) structures to the total muscle force were identified by using Blebbistatin. The experiments revealed three main results: (i) partial detachment of XBs during the eccentric phase of a SSC, (ii) significantly enhanced forces and mechanical work during the concentric phase of SSCs compared with shortening contractions with and without XB-inhibition, and (iii) no residual force depression after SSCs. The results obtained by administering Blebbistatin propose a titin-actin interaction that depends on XB-binding or active XB-based force production. The findings of this study further suggest that enhanced forces generated during the active lengthening phase of SSCs persist during the subsequent shortening phase, thereby contributing to enhanced work. Accordingly, our data support the hypothesis that sarcomeric mechanisms related to residual force enhancement also contribute to the SSC-effect. The preload of the titin molecule, acting as molecular spring, might be part of that mechanism by increasing the mechanical efficiency of work during physiological SSCs.

Review of the methods used for calculating physiological cross-sectional area (PCSA) for ecological questions

This review examines literature that used physiological cross-sectional area (PCSA) as a representative measure of an individual muscle's maximal isometric force production. PCSA is used to understand the muscle architecture and how a trade-off between muscle force and muscle contractile velocity reflect adaptations of the musculoskeletal system as a reflection of functional demands. Over the decades, methods have been developed to measure muscle volume, fascicle lengths, and pennation angle to calculate PCSA. The advantages and limitations of these methods (especially the inclusion/elimination of pennation angle) are discussed frequently; however, these method descriptions are scattered throughout the literature. Here, we reviewed and summarised the different approaches to collecting and recording muscle architectural properties to subsequently calculate PCSA. By critically discussing the advantages and limitations of each methodology, we aim to provide readers with an overview of repeatable methods to assess muscle architecture. This review may serve as a guide to facilitate readers searching for the appropriate techniques to calculate PCSA and measure muscle architecture to be applied in ecomorphology research. RESEARCH HIGHLIGHTS: Discuss the theories behind PCSA in a synthesised review to inform researchers about PCSA methodology.

Exhaustion of Skeletal Muscle Fibers Within Seconds: Incorporating Phosphate Kinetics Into a Hill-Type Model

Initiated by neural impulses and subsequent calcium release, skeletal muscle fibers contract (actively generate force) as a result of repetitive power strokes of acto-myosin cross-bridges. The energy required for performing these cross-bridge cycles is provided by the hydrolysis of adenosine triphosphate (ATP). The reaction products, adenosine diphosphate (ADP) and inorganic phosphate (P i ), are then used-among other reactants, such as creatine phosphate-to refuel the ATP energy storage. However, similar to yeasts that perish at the hands of their own waste, the hydrolysis reaction products diminish the chemical potential of ATP and thus inhibit the muscle's force generation as their concentration rises. We suggest to use the term "exhaustion" for force reduction (fatigue) that is caused by combined P i and ADP accumulation along with a possible reduction in ATP concentration. On the basis of bio-chemical kinetics, we present a model of muscle fiber exhaustion based on hydrolytic ATP-ADP-P i dynamics, which are assumed to be length- and calcium activity-dependent. Written in terms of differential-algebraic equations, the new sub-model allows to enhance existing Hill-type excitation-contraction models in a straightforward way. Measured time courses of force decay during isometric contractions of rabbit M. gastrocnemius and M. plantaris were employed for model verification, with the finding that our suggested model enhancement proved eminently promising. We discuss implications of our model approach for enhancing muscle models in general, as well as a few aspects regarding the significance of phosphate kinetics as one contributor to muscle fatigue.

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.

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.

Extensive eccentric contractions in intact cardiac trabeculae: revealing compelling differences in contractile behaviour compared to skeletal muscles

Force enhancement (FE) is a phenomenon that is present in skeletal muscle. It is characterized by progressive forces upon active stretching-distinguished by a linear rise in force-and enhanced isometric force following stretching (residual FE (RFE)). In skeletal muscle, non-cross-bridge (XB) structures may account for this behaviour. So far, it is unknown whether differences between non-XB structures within the heart and skeletal muscle result in deviating contractile behaviour during and after eccentric contractions. Thus, we investigated the force response of intact cardiac trabeculae during and after isokinetic eccentric muscle contractions (10% of maximum shortening velocity) with extensive magnitudes of stretch (25% of optimum muscle length). The different contributions of XB and non-XB structures to the total muscle force were revealed by using an actomyosin inhibitor. For cardiac trabeculae, we found that the force-length dynamics during long stretch were similar to the total isometric force-length relation. This indicates that no (R)FE is present in cardiac muscle while stretching the muscle from 0.75 to 1.0 optimum muscle length. This finding is in contrast with the results obtained for skeletal muscle, in which (R)FE is present. Our data support the hypothesis that titin stiffness does not increase with activation in cardiac muscle.

Locational and Directional Dependencies of Smooth Muscle Properties in Pig Urinary Bladder

The urinary bladder is a distensible hollow muscular organ, which allows huge changes in size during absorption, storage and micturition. Pathological alterations of biomechanical properties can lead to bladder dysfunction and loss in quality of life. To understand and treat bladder diseases, the mechanisms of the healthy urinary bladder need to be determined. Thus, a series of studies focused on the detrusor muscle, a layer of urinary bladder made of smooth muscle fibers arranged in longitudinal and circumferential orientation. However, little is known about whether its active muscle properties differ depending on location and direction. This study aimed to investigate the porcine bladder for heterogeneous (six different locations) and anisotropic (longitudinal vs. circumferential) contractile properties including the force-length-(FLR) and force-velocity-relationship (FVR). Therefore, smooth muscle tissue strips with longitudinal and circumferential direction have been prepared from different bladder locations (apex dorsal, apex ventral, body dorsal, body ventral, trigone dorsal, trigone ventral). FLR and FVR have been determined by a series of isometric and isotonic contractions. Additionally, histological analyses were conducted to determine smooth muscle content and fiber orientation. Mechanical and histological examinations were carried out on 94 and 36 samples, respectively. The results showed that maximum active stress (p_act ) of the bladder strips was higher in the longitudinal compared to the circumferential direction. This is in line with our histological investigation showing a higher smooth muscle content in the bladder strips in the longitudinal direction. However, normalization of maximum strip force by the cross-sectional area (CSA) of smooth muscle fibers yielded similar smooth muscle maximum stresses (165.4 ± 29.6 kPa), independent of strip direction. Active muscle properties (FLR, FVR) showed no locational differences. The trigone exhibited higher passive stress (p_pass ) than the body. Moreover, the bladder exhibited greater p_pass in the longitudinal than circumferential direction which might be attributed to its microstructure (more longitudinal arrangement of muscle fibers). This study provides a valuable dataset for the development of constitutive computational models of the healthy urinary bladder. These models are relevant from a medical standpoint, as they contribute to the basic understanding of the function of the bladder in health and disease.

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.

The active force-length relationship is invisible during extensive eccentric contractions in skinned skeletal muscle fibres

In contrast to experimentally observed progressive forces in eccentric contractions, cross-bridge and sliding-filament theories of muscle contraction predict that varying myofilament overlap will lead to increases and decreases in active force during eccentric contractions. Non-cross-bridge contributions potentially explain the progressive total forces. However, it is not clear whether underlying abrupt changes in the slope of the nonlinear force-length relationship are visible in long isokinetic stretches, and in which proportion cross-bridges and non-cross-bridges contribute to muscle force. Here, we show that maximally activated single skinned rat muscle fibres behave (almost across the entire working range) like linear springs. The force slope is about three times the maximum isometric force per optimal length. Cross-bridge and non-cross-bridge contributions to the muscle force were investigated using an actomyosin inhibitor. The experiments revealed a nonlinear progressive contribution of non-cross-bridge forces and suggest a nonlinear cross-bridge contribution similar to the active force-length relationship (though with increased optimal length and maximum isometric force). The linear muscle behaviour might significantly reduce the control effort. Moreover, the observed slight increase in slope with initial length is in accordance with current models attributing the non-cross-bridge force to titin.

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.

A novel diffusion-tensor MRI approach for skeletal muscle fascicle length measurements

Musculoskeletal (dys‐)function relies for a large part on muscle architecture which can be obtained using Diffusion‐Tensor MRI (DT‐MRI) and fiber tractography. However, reconstructed tracts often continue along the tendon or aponeurosis when using conventional methods, thus overestimating fascicle lengths. In this study, we propose a new method for semiautomatic segmentation of tendinous tissue using tract density (TD). We investigated the feasibility and repeatability of this method to quantify the mean fascicle length per muscle. Additionally, we examined whether the method facilitates measuring changes in fascicle length of lower leg muscles with different foot positions. Five healthy subjects underwent two DT‐MRI scans of the right lower leg, with the foot in 15° dorsiflexion, neutral, and 30° plantarflexion positions. Repeatability of fascicle length measurements was assessed using Bland–Altman analysis. Changes in fascicle lengths between the foot positions were tested using a repeated multivariate analysis of variance (MANOVA). Bland–Altman analysis showed good agreement between repeated measurements. The coefficients of variation in neutral position were 8.3, 16.7, 11.2, and 10.4% for soleus (SOL), fibularis longus (FL), extensor digitorum longus (EDL), and tibialis anterior (TA), respectively. The plantarflexors (SOL and FL) showed significant increase in fascicle length from plantarflexion to dorsiflexion, whereas the dorsiflexors (EDL and TA) exhibited a significant decrease. The use of a tract density for semiautomatic segmentation of tendinous structures provides more accurate estimates of the mean fascicle length than traditional fiber tractography methods. The method shows moderate to good repeatability and allows for quantification of changes in fascicle lengths due to passive stretch.

Porcine Stomach Smooth Muscle Force Depends on History-Effects

The stomach serves as food reservoir, mixing organ and absorption area for certain substances, while continually varying its position and size. Large dimensional changes during ingestion and gastric emptying of the stomach are associated with large changes in smooth muscle length. These length changes might induce history-effects, namely force depression (FD) following active muscle shortening and force enhancement (FE) following active muscle stretch. Both effects have impact on the force generating capacity of the stomach, and thus functional relevance. However, less is known about history-effects and active smooth muscle properties of stomach smooth muscle. Thus, the aim of this study was to investigate biomechanical muscle properties as force-length and force-velocity relations (FVR) of porcine stomach smooth muscle strips, extended by the analysis of history-effects on smooth muscle force. Therefore, in total n = 54 tissue strips were dissected in longitudinal direction from the ventral fundus of porcine stomachs. Different isometric, isotonic, and isokinetic contraction protocols were performed during electrical muscle stimulation. Cross-sectional areas (CSA) of smooth muscles were determined from cryo-histological sections stained with Picrosirius Red. Results revealed that maximum smooth muscle tension was 10.4 ± 2.6 N/cm². Maximum shortening velocity (V_max) and curvature factor (curv) of the FVR were 0.04 ± 0.01 [optimum muscle length/s] and 0.36 ± 0.15, respectively. The findings of the present study demonstrated significant (P < 0.05) FD [up to 32% maximum muscle force (F_im)] and FE (up to 16% F_im) of gastric muscle tissue, respectively. The FE- and FD-values increased with increasing ramp amplitude. This outstanding muscle behavior is not accounted for in existing models so far and strongly supports the idea of a holistic reflection of distinct stomach structure and function. For the first time this study provides a comprehensive set of stomach smooth muscle parameters including classic biomechanical muscle properties and history-dependent effects, offering the possibility for the development and validation of computational stomach models. Furthermore, this data set facilitates novel insights in gastric motility and contraction behavior based on the re-evaluation of existing contractile mechanisms. That will likely help to understand physiological functions or dysfunctions in terms of gastric accommodation and emptying.

A continuum-mechanical skeletal muscle model including actin-titin interaction predicts stable contractions on the descending limb of the force-length relation

Contractions on the descending limb of the total (active + passive) muscle force-length relationship (i. e. when muscle stiffness is negative) are expected to lead to vast half-sarcomere-length inhomogeneities. This is however not observed in experiments-vast half-sarcomere-length inhomogeneities can be absent in myofibrils contracting in this range, and initial inhomogeneities can even decrease. Here we show that the absence of half-sarcomere-length inhomogeneities can be predicted when considering interactions of the semi-active protein titin with the actin filaments. Including a model of actin-titin interactions within a multi-scale continuum-mechanical model, we demonstrate that stability, accurate forces and nearly homogeneous half-sarcomere lengths can be obtained on the descending limb of the static total force-length relation. This could be a key to durable functioning of the muscle because large local stretches, that might harm, for example, the transverse-tubule system, are avoided.

A diffusion-weighted imaging informed continuum model of the rabbit triceps surae complex

The NZ white rabbit is the animal of choice for much experimental work due to its muscular frame and similar response to human diseases, and is one of the few mammals that have had their genome sequenced. However, continuum-level computational models of rabbit muscle detailing fibre architecture are limited in the literature, especially the triceps surae complex (gastrocnemius, plantaris and soleus), which has similar biomechanics and translatable findings to the human. This study presents a geometrical model of the rabbit triceps surae informed with diffusion-weighted imaging (DWI)-based fibres. Passive rabbit-specific material properties are estimated using known muscle deformation inferred from magnetic resonance imaging data and dorsiflexion force measured with a custom-built rabbit rig and transducer. Muscle shape prediction is evaluated against a second rabbit. This study revealed that the triceps surae steady-state force post-rigor is close to post-mortem for small deformations but increases by a fixed ratio as the deformation increases and can be used to evaluate the passive behaviour of muscle. DWI fibre orientation significantly influences shape and mechanics during simulated computational muscle contraction. The presented triceps surae force and material properties may be used to inform the constitutive behaviour of continuum rabbit muscle models used to investigate pathology and musculotendon treatments that may be translated to the human condition.

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.

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.

A multi-scale continuum model of skeletal muscle mechanics predicting force enhancement based on actin-titin interaction

Although recent research emphasises the possible role of titin in skeletal muscle force enhancement, this property is commonly ignored in current computational models. This work presents the first biophysically based continuum-mechanical model of skeletal muscle that considers, in addition to actin-myosin interactions, force enhancement based on actin-titin interactions. During activation, titin attaches to actin filaments, which results in a significant reduction in titin's free molecular spring length and therefore results in increased titin forces during a subsequent stretch. The mechanical behaviour of titin is included on the microscopic half-sarcomere level of a multi-scale chemo-electro-mechanical muscle model, which is based on the classic sliding-filament and cross-bridge theories. In addition to titin stress contributions in the muscle fibre direction, the continuum-mechanical constitutive relation accounts for geometrically motivated, titin-induced stresses acting in the muscle's cross-fibre directions. Representative simulations of active stretches under maximal and submaximal activation levels predict realistic magnitudes of force enhancement in fibre direction. For example, stretching the model by 20 % from optimal length increased the isometric force at the target length by about 30 %. Predicted titin-induced stresses in the muscle's cross-fibre directions are rather insignificant. Including the presented development in future continuum-mechanical models of muscle function in dynamic situations will lead to more accurate model predictions during and after lengthening contractions.

Fast and Powerful: Biomechanics and Bite Forces of the Mandibles in the American Cockroach Periplaneta americana

Knowing the functionality and capabilities of masticatory apparatuses is essential for the ecological classification of jawed organisms. Nevertheless insects, especially with their outstanding high species number providing an overwhelming morphological diversity, are notoriously underexplored with respect to maximum bite forces and their dependency on the mandible opening angles. Aiming for a general understanding of insect biting, we examined the generalist feeding cockroach Periplaneta americana, characterized by its primitive chewing mouth parts. We measured active isometric bite forces and passive forces caused by joint resistance over the entire mandibular range with a custom-built 2D force transducer. The opening angle of the mandibles was quantified by using a video system. With respect to the effective mechanical advantage of the mandibles and the cross-section areas, we calculated the forces exerted by the mandible closer muscles and the corresponding muscle stress values. Comparisons with the scarce data available revealed close similarities of the cockroaches' mandible closer stress values (58 N/cm2) to that of smaller specialist carnivorous ground beetles, but strikingly higher values than in larger stag beetles. In contrast to available datasets our results imply the activity of faster and slower muscle fibres, with the latter becoming active only when the animals chew on tough material which requires repetitive, hard biting. Under such circumstances the coactivity of fast and slow fibres provides a force boost which is not available during short-term activities, since long latencies prevent a specific effective employment of the slow fibres in this case.