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

Linking myosin heavy chain isoform shift to mechanical properties and fracture modes in skeletal muscle tissue

Muscle fibers play a crucial role in the mechanical action of skeletal muscle tissue. However, it is unclear how the histological variations affect the mechanical properties of tissues. In this study, the shift of myosin heavy chain (MHC) isoforms is used for the first time to establish a linkage between tissue histological variation and passive mechanical properties. The shift of MHC isoform is found not only to induce significant differences in skeletal muscle passive mechanical properties, but also to lead to differences in strain rate responses. Non-negligible rate dependence is observed even in the conventionally defined quasi-static regime. Fidelity in the estimated constitutive parameters, which can be impacted due to variation in MHC isoforms and hence in rate sensitivity, is enhanced using a Bayesian inference framework. Subsequently, scanning electron microscopy and fluorescence microscopy are used to characterize the fracture morphology of muscle tissues and fibers. The fracture mode of both MHC I and II muscle fibers exhibited shearing of endomysium. Results show that the increase in strain rate only leads to stronger rebounding of the muscle fibers during tissue rupture without changing fracture modes.

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.

Experimentally characterizing the spatially varying anisotropic mechanical property of cancellous bone via a Bayesian calibration method

Traditional experimental tests for characterizing bone's mechanical properties usually hypothesize a uniaxial stress condition without quantitatively evaluating the influence of spatially varying principal material orientations, which cannot accurately predict the mechanical properties distribution of bones in vivo environment. In this study, a Bayesian calibrating procedure was developed using quantified multiaxial stress to investigate cancellous bone's local anisotropic elastic performance around joints as the spatial variation of main bearing orientations. First, the bone cube specimens from the distal femur of sheep are prepared using traditional anatomical axes. The multiaxial stress state of each bone specimen is calibrated using the actual principal material orientations derived from fabric tensor at different anatomical locations. Based on the calibrated multiaxial stress state, the process of identifying mechanical properties is described as an inverse problem. Then, a Bayesian calibration procedure based on a surrogate constitutive model was developed via multiaxial stress correction to identify the anisotropic material parameters. Finally, a comparison between the experiment and simulation results is discussed by applying the optimal model parameters obtained from the Bayesian probability distribution. Compared to traditional uniaxial methods, our results prove that the calibration based on the spatial variation of the main bearing orientations can significantly improve the accuracy of characterizing regional anisotropic mechanical responses. Moreover, we determine that the actual mechanical property distribution is influenced by complicated mechanical stimulation. This study provides a novel method to evaluate the spatially varying mechanical properties of bone tissues enduring complex mechanical loading accurately and effectively. It is expected to provide more realistic mechanical design targets in vivo for a personalized artificial bone prosthesis in clinical treatment.

On multiscale tension-compression asymmetry in skeletal muscle

Skeletal muscle tissue shows a clear asymmetry with regard to the passive stresses under tensile and compressive deformation, referred to as tension-compression asymmetry (TCA). The present study is the first one reporting on TCA at different length scales, associated with muscle tissue and muscle fibres, respectively. This allows for the first time the comparison of TCA between the tissue and one of its individual components, and thus to identify the length scale at which this phenomenon originates. Not only the passive stress-stretch characteristics were recorded, but also the volume changes during the axial tension and compression experiments. The study reveals clear differences in the characteristics of TCA between fibres and tissue. At tissue level TCA increases non-linearly with increasing deformation and the ratio of tensile to compressive stresses at the same magnitude of strain reaches a value of approximately 130 at 13.5% deformation. At fibre level instead it initially drops to a value of 6 and then rises again to a TCA of 14. At a deformation of 13.5%, the tensile stress is about 6 times higher. Thus, TCA is about 22 times more expressed at tissue than fibre scale. Moreover, the analysis of volume changes revealed little compressibility at tissue scale whereas at fibre level, especially under compressive stress, the volume decreases significantly. The data collected in this study suggests that the extracellular matrix has a distinct role in amplifying the TCA, and leads to more incompressible tissue behaviour. STATEMENT OF SIGNIFICANCE: This article analyses and compares for the first time the tension-compression asymmetry (TCA) displayed by skeletal muscle at tissue and fibre scale. In addition, the volume changes of tissue and fibre specimens with application of passive tensile and compressive loads are studied. The study identifies a key role of the extracellular matrix in establishing the mechanical response of skeletal muscle tissue: It contributes significantly to the passive stress, it is responsible for the major part of tissue-scale TCA and, most probably, prevents/balances the volume changes of muscle fibres during deformation. These new results thus shed light on the origin of TCA and provide new information to be used in microstructure-based approaches to model and simulate skeletal muscle tissue.

Surrogate-based Bayesian calibration of biomechanical models with isotropic material behavior

This work introduces a computational methodology to calibrate material models in biomechanical applications under uncertainty. We adopt a Bayesian approach, which estimates the probability distributions of hyperelastic material parameters, based on force-strain measurements. We approximate the parametric biomechanical model by combining a reduced order representation of the force response with a Polynomial Chaos expansion. The surrogate model allows to employ sampling-intensive Markov chain Monte Carlo methods and provides an efficient way to estimate (generalized) Sobol coefficients. We use a Sobol sensitivity analysis to assess the influence of material parameters and present an iterative procedure to quantify the accuracy of the surrogate model as additional uncertainty during Bayesian updating. The methodology is illustrated with three cases, tensile experiments on heat-induced whey protein gel, indentation experiments for oocytes and a manufactured example. Real experimental data are used for the calibration.

Multiscale modeling of skeletal muscle to explore its passive mechanical properties and experiments verification

The research of the mechanical properties of skeletal muscle has never stopped, whether in experimental tests or simulations of passive mechanical properties. To investigate the effect of biomechanical properties of micro-components and geometric structure of muscle fibers on macroscopic mechanical behavior, in this manuscript, we establish a multiscale model where constitutive models are proposed for fibers and the extracellular matrix, respectively. Besides, based on the assumption that the fiber cross-section can be expressed by Voronoi polygons, we optimize the Voronoi polygons as curved-edge Voronoi polygons to compare the effects of the two cross-sections on macroscopic mechanical properties. Finally, the macroscopic stress response is obtained through the numerical homogenization method. To verify the effectiveness of the multi-scale model, we measure the mechanical response of skeletal muscles in the in-plane shear, longitudinal shear, and tensions, including along the fiber direction and perpendicular to the fiber direction. Compared with experimental data, the simulation results show that this multiscale framework predicts both the tension response and the shear response of skeletal muscle accurately. The root mean squared error (RMSE) is 0.0035 MPa in the tension along the fiber direction; The RMSE is 0.011254 MPa in the tension perpendicular to the fiber direction; The RMSE is 0.000602 MPa in the in-plane shear; The RMSE was 0.00085 MPa in the longitudinal shear. Finally, we obtained the influence of the component constitutive model and muscle fiber cross-section on the macroscopic mechanical behavior of skeletal muscle. In terms of the tension perpendicular to the fiber direction, the curved-edge Voronoi polygons achieve the result closer to the experimental data than the Voronoi polygons. Skeletal muscle mechanics experiments verify the effectiveness of our multiscale model. The comparison results of experiments and simulations prove that our model can accurately capture the tension and shear behavior of skeletal muscle.

Systematic review of skeletal muscle passive mechanics experimental methodology

Understanding passive skeletal muscle mechanics is critical in defining structure-function relationships in skeletal muscle and ultimately understanding pathologically impaired muscle. In this systematic review, we performed an exhaustive literature search using PRISMA guidelines to quantify passive muscle mechanical properties, summarized the methods used to create these data, and make recommendations to standardize future studies. We screened over 7500 papers and found 80 papers that met the inclusion criteria. These papers reported passive muscle mechanics from single muscle fiber to whole muscle across 16 species and 54 distinct muscles. We found a wide range of methodological differences in sample selection, preparation, testing, and analysis. The systematic review revealed that passive muscle mechanics is species and scale dependent-specifically within mammals, the passive mechanics increases non-linearly with scale. However, a detailed understanding of passive mechanics is still unclear because the varied methodologies impede comparisons across studies, scales, species, and muscles. Therefore, we recommend the following: smaller scales may be maintained within storage solution prior to testing, when samples are tested statically use 2-3 min of relaxation time, stress normalization at the whole muscle level be to physiologic cross-sectional area, strain normalization be to sarcomere length when possible, and an exponential equation be used to fit the data. Additional studies using these recommendations will allow exploration of the multiscale relationship of passive force within and across species to provide the fundamental knowledge needed to improve our understanding of passive muscle mechanics.

Modelling extracellular matrix and cellular contributions to whole muscle mechanics

Skeletal muscle tissue has a highly complex and heterogeneous structure comprising several physical length scales. In the simplest model of muscle tissue, it can be represented as a one dimensional nonlinear spring in the direction of muscle fibres. However, at the finest level, muscle tissue includes a complex network of collagen fibres, actin and myosin proteins, and other cellular materials. This study shall derive an intermediate physical model which encapsulates the major contributions of the muscle components to the elastic response apart from activation-related along-fibre responses. The micro-mechanical factors in skeletal muscle tissue (eg. connective tissue, fluid, and fibres) can be homogenized into one material aggregate that will capture the behaviour of the combination of material components. In order to do this, the corresponding volume fractions for each type of material need to be determined by comparing the stress-strain relationship for a volume containing each material. This results in a model that accounts for the micro-mechanical features found in muscle and can therefore be used to analyze effects of neuro-muscular diseases such as cerebral palsy or muscular dystrophies. The purpose of this study is to construct a model of muscle tissue that, through choosing the correct material parameters based on experimental data, will accurately capture the mechanical behaviour of whole muscle. This model is then used to look at the impacts of the bulk modulus and material parameters on muscle deformation and strain energy-density distributions.

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.

Predicting muscle tissue response from calibrated component models and histology-based finite element models

Skeletal muscle is an anisotropic soft biological tissue composed of muscle fibres embedded in a structurally complex, hierarchically organised extracellular matrix. In a recent work (Kuravi et al., 2021) we have developed 3D finite element models from series of histological sections. Moreover, based on decellularisation of fresh tissue samples, a novel set of experimental data on the direction dependent mechanical properties of collagenous ECM was established (Kohn et al., 2021). Together with existing information on the material properties of single muscle fibres, the combination of these techniques allows computing predictions of the composite tissue response. To this end, an inverse finite element procedure is proposed in the present work to calibrate a constitutive model of the extracellular matrix, and supplementary biaxial tensile tests on fresh and decellularised tissues are performed for model validation. The results of this rigorously predictive and thus unforgiving strategy suggest that the prediction of the tissue response from the individual characteristics of muscle cells and decellularised tissue is only possible within clear limits. While orders of magnitude are well matched, and the qualitative behaviour in a wide range of load cases is largely captured, the existing deviations point at potentially missing components of the model and highlight the incomplete experimental information in bottom-up multiscale approaches to model skeletal muscle tissue.

3D finite element models from serial section histology of skeletal muscle tissue - The role of micro-architecture on mechanical behaviour

In this contribution we create three-dimensional (3D) finite element models from a series of histological sections of porcine skeletal muscle tissue. Image registration is performed on the stained sections by affinely aligning them using auxiliary markers, followed by image segmentation to determine muscle fibres and the extracellular matrix in each section, with particular regard to the continuity of the fibres through the stack. With this information, 3D virtual tissue samples are reconstructed, discretised, and associated with appropriate non-linear elastic anisotropic material models. While the gross anatomy is directly obtained from the images, the local directions of anisotropy were determined by the use of an analogy with steady state diffusion. The influence of the number of histological sections considered for reconstruction on the numerically simulated mechanical response of the virtual tissue samples is then studied. The results show that muscle tissue is fairly heterogeneous along the fascicles, and that transverse isotropy is inadequate in describing their material symmetry at the typical length scale of a fascicle. Numerical simulations of different load cases suggest that ignoring the undulations of fibres and their non-uniform cross-sections only moderately affects the passive response of the tissue in tensile and compressive modes, but can become crucial when predicting the response to generic loads and activation.

Mechano-geometrical skeletal muscle fibre characterisation under cyclic and relaxation loading

In the present work, mechano-geometrical characterisations of skeletal muscle fibres in two different deformation states, namely, axial tension and axial compression, were realised. In both cases, cyclic and relaxation tests were performed. Additionally, the changes in the volume of the fibres during deformation were recorded to obtain more detailed information about the muscle fibre load transfer mechanisms. To the best of the authors' knowledge, the present experimental investigation of the mechanical and geometrical characteristics of muscle fibres provides a novel comprehensive data set that can be used to obtain a better understanding of muscle fibre load transfer mechanisms and to construct meaningful models. In the present study, it is shown that muscle fibres exhibit incompressibility (5% volume decrease at maximum deformation) under tension and that this feature is more pronounced under compression loading (37% volume decrease at maximum deformation). These findings are particularly interesting and lead to a further understanding of load transfer mechanisms and to the development of new modelling strategies.

An experimental study of nonlinear rate-dependent behaviour of skeletal muscle to obtain passive mechanical properties

Accurate modelling of biological tissues has been a significant need for analysis of the human body. In this article, a comprehensive in vitro experimental study has been done on the fresh bovine skeletal muscle before the onset of rigour mortis in order to provide an experimental description of passive skeletal muscle properties in three dimensions. Different situations including various deformation modes, different loading rates and loading directions are tested to consider all features of skeletal muscle behaviour. Based on the nonlinear continuum mechanics, a three-dimensional visco-hyperelastic model is introduced which considers all aspects of skeletal muscle's features such as nonlinear hyperelastic, time-dependent behaviour, anisotropy and quasi-incompressibility. Visco-hyperelastic material constants are obtained for passive behaviour of the muscle based on genetic algorithm optimization method via comparing the theoretical and experimental results. Experiments show that the rate of loading affects the configuration of experimental curves considerably. It could be also concluded that compression-tension asymmetry, as well as anisotropic behaviour, of the muscle is due to fibres orientation. Obtained experimental results help to achieve a better understanding of mechanical properties and nonlinear behaviour of the skeletal muscles.

An inverse model of the mechanical response of passive skeletal muscle: Implications for microstructure

The constitutive response of passive skeletal muscle is important for many human body modelling applications, but modelling the tension-compression asymmetry and the anisotropy observed in ex-vivo samples is challenging. Existing microstructural models do not capture the full three-dimensional response while models suitable for application in finite element environments mostly have a limited microstructural basis and cannot capture the observed Poisson's ratios. The aim of this paper is to derive an inverse model based on the microstructure of a skeletal muscle that can predict its passive mechanical response. The model parameters and predictions were derived and assessed by comparison with published experimental stress-strain response and Poisson's ratio data. Results show a close match for both predicted stress-strain response for fibre and cross-fibre direction deformations and similar Poisson's ratio values. Some microstructural observations which strengthen our understanding of the role of the collagen network and intramuscular pressure are also provided.

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.