Functional analysis of limb recovery following autograft treatment of volumetric muscle loss in the quadriceps femoris

Severe injuries to the extremities often result in muscle trauma and, in some cases, significant volumetric muscle loss (VML). These injuries continue to be challenging to treat, with few available clinical options, a high rate of complications, and often persistent loss of limb function. To facilitate the testing of regenerative strategies for skeletal muscle, we developed a novel quadriceps VML model in the rat, specifically addressing functional recovery of the limb. Our outcome measures included muscle contractility measurements to assess muscle function and gait analysis for evaluation of overall limb function. We also investigated treatment with muscle autografts, whole or minced, to promote regeneration of the defect area. Our defect model resulted in a loss of muscle function, with injured legs generating less than 55% of muscle strength from the contralateral uninjured control legs, even at 4 weeks post-injury. The autograft treatments did not result in significant recovery of muscle function. Measures of static and dynamic gait were significantly decreased in the untreated, empty defect group, indicating a decrease in limb function. Histological sections of the affected muscles showed extensive fibrosis, suggesting that this scarring of the muscle may be in part the cause of the loss of muscle function in this VML model. Taken together, these data are consistent with clinical findings of reduced muscle function in large VML injuries. This new model with quantitative functional outcome measures offers a platform on which to evaluate treatment strategies designed to regenerate muscle tissue volume and restore limb function.     

In vitro estimates of power output by epaxial muscle during feeding in largemouth bass

Recent work has employed video and sonometric analysis combined with hydrodynamic modeling to estimate power output by the feeding musculature of largemouth bass in feeding trials. The result was an estimate of approximately 69 W kg(-1) of power by the epaxial muscle during maximal feeding strikes. The present study employed in vitro measurements of force, work and power output by fast-twitch epaxial muscle bundles stimulated under activation conditions measured in vivo to evaluate the power output results of the feeding experiments. Isolated muscle bundles from the epaxial muscle, the sternohyoideus and the lateral red or slow-twitch muscle were tied into a muscle mechanics apparatus, and contractile properties during tetanic contractions and maximum shortening velocity (Vmax) were determined. For the epaxial muscles, work and power output during feeding events was determined by employing mean stimulation conditions derived from a select set of maximal feeding trials: 17% muscle shortening at 3.6 muscle lengths/s, with activation occurring 5 ms before the onset of shortening. Epaxial and sternohyoideus muscle displayed similar contractile properties, and both were considerably faster (Vmax approximately 11-13 ML s(-1)) than red muscle (Vmax approximately 5 ML s(-1)). Epaxial muscle stimulated under in vivo activation conditions generated approximately 60 W kg(-1) with a 17% strain and approximately 86 W kg(-1) with a 12% strain. These values are close to those estimated by hydrodynamic modeling. The short lag time (5 ms) between muscle activation and muscle shortening is apparently a limiting parameter during feeding strikes, with maximum power found at an offset of 15-20 ms. Further, feeding strikes employing a faster shortening velocity generated significantly higher power output. Power production during feeding strikes appears to be limited by the need for fast onset of movement and the hydrodynamic resistance to buccal expansion.     

Subject-specific knee joint geometry improves predictions of medial tibiofemoral contact forces

Estimating tibiofemoral joint contact forces is important for understanding the initiation and progression of knee osteoarthritis. However, tibiofemoral contact force predictions are influenced by many factors including muscle forces and anatomical representations of the knee joint. This study aimed to investigate the influence of subject-specific geometry and knee joint kinematics on the prediction of tibiofemoral contact forces using a calibrated EMG-driven neuromusculoskeletal model of the knee. One participant fitted with an instrumented total knee replacement walked at a self-selected speed while medial and lateral tibiofemoral contact forces, ground reaction forces, whole-body kinematics, and lower-limb muscle activity were simultaneously measured. The combination of generic and subject-specific knee joint geometry and kinematics resulted in four different OpenSim models used to estimate muscle–tendon lengths and moment arms. The subject-specific geometric model was created from CT scans and the subject-specific knee joint kinematics representing the translation of the tibia relative to the femur was obtained from fluoroscopy. The EMG-driven model was calibrated using one walking trial, but with three different cost functions that tracked the knee flexion/extension moments with and without constraint over the estimated joint contact forces. The calibrated models then predicted the medial and lateral tibiofemoral contact forces for five other different walking trials. The use of subject-specific models with minimization of the peak tibiofemoral contact forces improved the accuracy of medial contact forces by 47% and lateral contact forces by 7%, respectively compared with the use of generic musculoskeletal model.     

Effects of inter- and extramuscular myofascial force transmission on adjacent synergistic muscles: assessment by experiments and finite-element modeling

The effects of inter- and extramuscular myofascial force transmission on muscle length force characteristics were studied in rat. Connective tissues at the bellies of the experimental synergistic muscles of the anterior crural compartment were left intact. Extensor digitorium longus (EDL) muscle was lengthened distally whereas tibialis anterior (TA) and extensor hallucis longus (EHL) were kept at constant muscle–tendon complex length. Substantial differences were found in EDL force measured at the proximal and distal tendons (maximally 46% of the proximal force). EDL with intact inter- as well as extramuscular connections had an increased length range between active slack and optimum length compared to EDL with extramuscular connections exclusively: optimum muscle length was shifted by more than 2 mm. Distal EDL lengthening caused the distal force exerted by TA+EHL complex to decrease (approximately 17% of the initial force). This indicates increased intermuscular myofascial force transmission from TA+EHL muscle complex to EDL muscle.

Finite-element modeling showed that: (1) Inter- and extramuscular myofascial force transmission leads to a substantial distribution of the lengths of the sarcomeres arranged in series within muscle fibers. Distribution of stress within the muscle fibers showed that the muscle fiber cannot be considered as a unit exerting equal forces at both ends. (2) Increased heterogeneity of mean fiber sarcomere lengths (i.e., a “parallel” distribution of length of sarcomeres among different muscle fibers) is found, particularly at high muscle lengths. This also explains the shift in muscle optimum length to higher lengths.

It is concluded that inter- and extramuscular myofascial force transmission has substantial effects on muscle length–force characteristics.     

A general-purpose framework to simulate musculoskeletal system of human body: using a motion tracking approach*

Computation of muscle force patterns that produce specified movements of muscle-actuated dynamic models is an important and challenging problem. This problem is an undetermined one, and then a proper optimization is required to calculate muscle forces. The purpose of this paper is to develop a general model for calculating all muscle activation and force patterns in an arbitrary human body movement. For this aim, the equations of a multibody system forward dynamics, which is considered for skeletal system of the human body model, is derived using Lagrange–Euler formulation. Next, muscle contraction dynamics is added to this model and forward dynamics of an arbitrary musculoskeletal system is obtained. For optimization purpose, the obtained model is used in computed muscle control algorithm, and a closed-loop system for tracking desired motions is derived. Finally, a popular sport exercise, biceps curl, is simulated by using this algorithm and the validity of the obtained results is evaluated via EMG signals.     

The role of transmembrane proteins on force transmission in skeletal muscle

Lateral transmission of force from myofibers laterally to the surrounding extracellular matrix (ECM) via the transmembrane proteins between them is impaired in old muscles. Changes in geometrical and mechanical properties of ECM of skeletal muscle do not fully explain the impaired lateral transmission with aging. The objective of this study was to determine the role of transmembrane proteins on force transmission in skeletal muscle. In this study, a 2D finite element model of single muscle fiber composed of myofiber, ECM, and the transmembrane proteins between them was developed to determine how changes in spatial density and mechanical properties of transmembrane proteins affect the force transmission in skeletal muscle. We found that force transmission and stress distribution are not affected by mechanical stiffness of the transmembrane proteins due to its non-linear stress–strain relationship. Results also showed that the muscle fiber with insufficient transmembrane proteins near the end of muscle fiber transmitted less force than that with more proteins does. Higher stress was observed in myofiber, ECM, and proteins in the muscle fiber with fewer proteins.     

Qualitatively different cross-bridge attachments in fast and slow muscle fiber types

Contractile properties differ between skeletal, cardiac and smooth muscles as well as between various skeletal muscle fiber types. This functional diversity is thought to be mainly related to different speeds of myosin head pulling cycles, with the molecular mechanism of force generation being essentially the same. In this study, force-generating attachments of myosin heads were investigated by applying small perturbations of myosin head pulling cycles in stepwise stretch experiments on skeletal muscle fibers of different type. Slow fibers (frog tonic and rat slow-twitch) exhibited only a ‘slow-type’ of myosin head attachment over the entire activation range, while fast fibers (frog and rat fast-twitch) displayed a ‘slow-type’ of myosin head attachment at low levels of activation, and an up to 30-times faster type at high levels of activation. These observations indicate that there are qualitative differences between the mechanisms of myosin head attachment in slow and fast vertebrate skeletal muscle fibers.     

Predicting muscle force generation during fast-starts for the common carp Cyprinus carpio

Muscle contractile properties have been characterised for white myotomal muscle from the common carp Cyprinus carpio at 10, 15, and 20 °C. The time course of muscle force development was measured when one, two, or three stimuli were delivered at the onset of constant velocity shortening. As the shortening velocity increased several parameters decreased including the maximum force, the time course for the contraction and the relative duration of the deactivation compared to the activation. The maximum force and the relative rates of activation to deactivation for the contraction were relatively independent of temperature, whereas the duration of the contraction decreased with increasing temperature. A predictive model was developed which was based on fitting a modified Weibull distribution to these observations. The model was used to interpolate the expected contractile forces during cyclic length-changes. Measured and predicted values for force and power during such cyclic work-loop experiments showed an excellent agreement over the range of shortening regimes typically found during swimming behaviours. However, the predicted force was overestimated during the deactivation phase of the contractions when the shortening velocities exceeded those found during swimming. Accepted: 25 May 1999     

Temperature adaptation and the contractile properties of live muscle fibres from teleost fish

Summary The contractile properties of swimming muscles have been investigated in marine teleosts from Antarctic (Trematomus lepidorhinus, Pseudochaenichthys georgianus), temperate (Pollachius virens, Limanda limanda, Agonis cataphractus, Callionymus lyra), and tropical (Abudefduf abdominalis, Thalassoma duperreyi) latitudes. Small bundles of fast twitch fibres were isolated from anterior myotomes and/or the pectoral fin adductor profundis muscle (m. add. p). Live fibre preparations were viable for several days at in vivo temperatures, but became progressively inexcitable at higher or lower temperatures. The stimulation frequency required to produce fused isometric tetani increased from 50 Hz in Antarctic species at 0°C to around 400 Hz in tropical species at 25°C. Maximum isometric tension (P_o) was produced at the normal body temperature (NBT) of each species (Antarctic, 0–2°C; North Sea and Atlantic, 8–10°C; Indo-West Pacific, 23–25°C). P₀ values at physiological temperatures (200–300 kN·m^–2) were similar for Antarctic, temperate, and tropical species. A temperature induced tension hysteresis was observed in muscle fibres from some species. Exposure to <0°C in Antarctic and <2°C in temperate fish resulted in the temporary depression of tension over the whole experimental range, an effect reversed by incubation at higher temperatures. At normal body temperatures the half-times for activation and relaxation of twitch and tetanic tension increased in the order Antarctic>temperate>tropical species. Relaxation was generally much slower at temperatures <10°C in fibres from tropical than temperate fish. Q₁₀ values for these parameters at NBTs were 1.3 2.1 for tropical species, 1.7–2.6 for temperate species, and 1.6–3.5 for Antarctic species. The forcevelocity (P-V) relationship was studied in selected species using iso-velocity releases and the data below 0.8 P₀ iteratively fitted to Hill's equation. The P-V relation at NBT was found to be significantly less curved in Antarctic than temperate species. The unloaded contraction velocity (V_max) of fibres was positively correlated with NBT increasing from about 1 muscle fibre length·s^–;1 in an Antarctic fish (Trematomus lepidorhinus) at 1°C to around 16 muscle fibre lengths·s^–1 in a tropical species (Thalassoma duperreyi) at 24°C. It is concluded that although muscle contraction in Antarctic fish shows adaptations for low temperature function, the degree of compensation achieved in shortening speed and twitch kinetics is relatively modest.Abbreviations ET environmental temperature - m. add. p major adductor profundis - m. add. s. major adductor superficialis - NBT normal body temperature - P ₀ maximum isometric tension - P-V force velocity - SR sarcoplasmic reticulum - T 1/2 a half activation time - T 1/2 r half relaxation time - V _max unloaded contraction     

Strains at the myotendinous junction predicted by a micromechanical model

The goal of this work was to create a finite element micromechanical model of the myotendinous junction (MTJ) to examine how the structure and mechanics of the MTJ affect the local micro-scale strains experienced by muscle fibers. We validated the model through comparisons with histological longitudinal sections of muscles fixed in slack and stretched positions. The model predicted deformations of the A-bands within the fiber near the MTJ that were similar to those measured from the histological sections. We then used the model to predict the dependence of local fiber strains on activation and the mechanical properties of the endomysium. The model predicted that peak micro-scale strains increase with activation and as the compliance of the endomysium decreases. Analysis of the models revealed that, in passive stretch, local fiber strains are governed by the difference of the mechanical properties between the fibers and the endomysium. In active stretch, strain distributions are governed by the difference in cross-sectional area along the length of the tapered region of the fiber near the MTJ. The endomysium provides passive resistance that balances the active forces and prevents the tapered region of the fiber from undergoing excessive strain. These model predictions lead to the following hypotheses: (i) the increased likelihood of injury during active lengthening of muscle fibers may be due to the increase in peak strain with activation and (ii) endomysium may play a role in protecting fibers from injury by reducing the strains within the fiber at the MTJ.     

Combined finite-element and rigid-body analysis of human jaw joint dynamics

The jaw joint plays a crucial role in human mastication. It acts as a guidance for jaw movements and as a fulcrum for force generation. The joint is subjected to loading which causes tensions and deformations in its cartilaginous structures. These are assumed to be a major determinant for development, maintenance and also degeneration of the joint. To analyze the distribution of tensions and deformations in the cartilaginous structures of the jaw joint during jaw movement, a dynamical model of the human masticatory system has been constructed. Its movements are controlled by muscle activation. The articular cartilage layers and articular disc were included as finite-element (FE) models. As this combination of rigid-body and FE modeling had not been applied to musculoskeletal systems yet, its benefits and limitations were assessed by simulating both unloaded and loaded jaw movements. It was demonstrated that joint loads increase with muscle activation, irrespective of the external loads. With increasing joint load, the size of the stressed area of the articular surfaces was enlarged, whereas the peak stresses were much less affected. The results suggest that the articular disc enables distribution of local contact stresses over a much wider area of the very incongruent articular surfaces by transforming compressive principal stress into shear stress.     

Phosphorylation of titin modulates passive stiffness of cardiac muscle in a titin isoform-dependent manner

We investigated the effect of protein kinase A (PKA) on passive force in skinned cardiac tissues that express different isoforms of titin, i.e., stiff (N2B) and more compliant (N2BA) titins, at different levels. We used rat ventricular (RV), bovine left ventricular (BLV), and bovine left atrial (BLA) muscles (passive force: RV > BLV > BLA, with the ratio of N2B to N2BA titin, approximately 90:10, approximately 40:60, and approximately 10:90%, respectively) and found that N2B and N2BA isoforms can both be phosphorylated by PKA. Under the relaxed condition, sarcomere length was increased and then held constant for 30 min and the peak passive force, stress-relaxation, and steady-state passive force were determined. Following PKA treatment, passive force was significantly decreased in all muscle types with the effect greatest in RV, lowest in BLA, and intermediate in BLV. Fitting the stress-relaxation data to the sum of three exponential decay functions revealed that PKA blunts the magnitude of stress-relaxation and accelerates its time constants. To investigate whether or not PKA-induced decreases in passive force result from possible alteration of titin-thin filament interaction (e.g., via troponin I phosphorylation), we conducted the same experiments using RV preparations that had been treated with gelsolin to extract thin filaments. PKA decreased passive force in gelsolin-treated RV preparations with a magnitude similar to that observed in control preparations. PKA was also found to decrease restoring force in skinned ventricular myocytes of the rat that had been shortened to below the slack length. Finally, we investigated the effect of the beta-adrenergic receptor agonist isoprenaline on diastolic force in intact rat ventricular trabeculae. We found that isoprenaline phosphorylated titin and that it reduced diastolic force to a degree similar to that found in skinned RV preparations. Taken together, these results suggest that during beta-adrenergic stimulation, PKA increases ventricular compliance in a titin isoform-dependent manner.