Physiological and developmental implications of motor unit anatomy

There is increasing evidence that the architectural design and arrangement of the fibers within a motor unit have important physiological and developmental ramifications. Limited data, however, are available to directly address this issue. In the present study the physiological properties of one motor unit in each of seven cat tibialis anterior (TA) muscles were determined. Each of these units then was repetitively stimulated to deplete the glycogen in all muscle fibers within the unit. Subsequently, the length, type of ending, and spatial distribution of fibers sampled from these physiologically and histochemically typed motor units were determined. Four fast fatigable (FF), one fast, fatigue resistant (FR), and two slow (S) motor units (MU) were studied. The samples consisted of all those glycogen-depleted fibers (9-27) contained within a single fascicle or a circumscribed area of each of the motor unit territories. The mean fiber lengths for the two slow motor units were 35.9 and 45.5 mm. The mean fiber lengths for the fast motor unit samples ranged from 8.8 to 48.5 mm. Some fibers of both the fast and slow units reached lengths of 58 mm. Most of the fibers in the slow units extended the entire distance between the proximal and distal musculotendinous planes, had relatively constant cross-sectional areas, and terminated at the tendon as blunt endings. In contrast, the majority of the fibers in the fast units terminated intrafascicularly at one end, and the cross-sectional area decreased progressively along their lengths, that is, showed a tapering pattern for a significant proportion of their lengths. Therefore, the force generated by units that end midfascicularly would appear to be transmitted to connective tissue elements and/or adjacent fibers. All fibers of a fast unit within a fascicle were located at approximately the same proximo-distal location. Thus, developmentally the selection of muscle fibers by a motoneuron would seem to be influenced by their spatial distribution. The architectural complexities of motor units also have clear implications for the mechanical interactions of active and inactive motor units. For example, the tension capabilities of a motor unit may be influenced not only by the spatial arrangement of its own fibers, but also by the level of activation of neighboring motor units.     

Oxygen delivery to skeletal muscle fibers: effects of microvascular unit structure and control mechanisms

The number of perfused capillaries in skeletal muscle varies with muscle activation. With increasing activation, muscle fibers are recruited as motor units consisting of widely dispersed fibers, whereas capillaries are recruited as groups called microvascular units (MVUs) that supply several adjacent fibers. In this study, a theoretical model was used to examine the consequences of this spatial mismatch between the functional units of muscle activation and capillary perfusion. Diffusive oxygen transport was simulated in cross sections of skeletal muscle, including several MVUs and fibers from several motor units. Four alternative hypothetical mechanisms controlling capillary perfusion were considered. First, all capillaries adjacent to active fibers are perfused. Second, all MVUs containing capillaries adjacent to active fibers are perfused. Third, each MVU is perfused whenever oxygen levels at its feed arteriole fall below a threshold value. Fourth, each MVU is perfused whenever the average oxygen level at its capillaries falls below a threshold value. For each mechanism, the dependence of the fraction of perfused capillaries on the level of muscle activation was predicted. Comparison of the results led to the following conclusions. Control of perfusion by MVUs increases the fraction of perfused capillaries relative to control by individual capillaries. Control by arteriolar oxygen sensing leads to poor control of tissue oxygenation at high levels of muscle activation. Control of MVU perfusion by capillary oxygen sensing permits adequate tissue oxygenation over the full range of activation without resulting in perfusion of all MVUs containing capillaries adjacent to active fibers.     

Influence of amplitude cancellation on the simulated surface electromyogram.

The purpose of the study was to quantify the influence of selected motor unit properties and patterns of activity on amplitude cancellation in the simulated surface electromyogram (EMG). The study involved computer simulations of a motor unit population with physiologically defined recruitment and rate coding characteristics that activated muscle fibers whose potentials were recorded on the skin over the muscle. Amplitude cancellation was quantified as the percent difference in signal amplitude when motor unit potentials were summed before and after rectification. The simulations involved varying the level of activation for the motor unit population, the recording configuration, the upper limit of motor unit recruitment, peak discharge rates, the amount of motor unit synchronization, muscle fiber length, the thickness of the subcutaneous tissue, and the motor unit properties that change with advancing age. The results confirmed a previous experimental report (Day SJ and Hulliger M, J Neurophysiol 86: 2144-2158, 2001) that amplitude cancellation in the surface EMG can reach 62% at maximal activation. A decrease in the range of amplitudes of the motor unit potentials, as can occur during fatiguing contractions, increased amplitude cancellation up to approximately 85%. Differences in the amount of amplitude cancellation were observed across all simulated conditions, and resulted in substantial changes in the absolute magnitude of the EMG signal. The most profound factors influencing amplitude cancellation were the number of active motor units and the duration of the action potentials. The effects of amplitude cancellation were minimal (<5%) when the EMG amplitude was normalized to maximal values, with the exception of variations in peak discharge rate and recruitment range, which resulted in differences up to 17% in the normalized EMG signal across conditions. These results indicate the amount of amplitude cancellation that can occur in various experimental conditions and its influence on absolute and relative measures of EMG amplitude.     

Distribution and innervation of short, interdigitated muscle fibers in parallel-fibered muscles of the cat hindlimb

The cat hindlimb contains several long, biarticular strap muscles composed of parallel muscle fascicles that attach to short tendons. Three of these muscles--sartorius, tenuissimus, and semitendinosus--were studied by dissecting individual gold-stained fibers and determining the surface distribution of acetylcholinesterase-stained end-plate zones. In each muscle, fascicles were composed of muscle fibers that ran only part of the fascicle length and tapered to end as fine strands that interdigitated with other tapering fibers within the muscle mass. Most muscle fibers measured 2-3 cm in length. Fascicles of muscle fibers were crossed by short transverse bands of endplates (1 mm wide by 1-5 mm long) that were spaced at fairly regular intervals from the origin to the insertion of the muscle. The endplate pattern suggested that the fiber fascicles were organized into multiple longitudinal strips. In the sartorius, the temporospatial distribution of electromyographic (EMG) activity evoked by stimulating fine, longitudinal branches of the parent nerve confirmed that each strip was selectively innervated by a small subset of the motor axons. These axons appeared to distribute their endings throughout the entire length of the fascicles, providing for synchronous activation of their in-series fibers.     

Heterogeneous fascicle behavior within the biceps femoris long head at different muscle activation levels

Magnetic resonance and ultrasound imaging have shown hamstring strain injuries occur most often in the biceps femoris long head (BFLH), and particularly in the proximal vs. distal region of this muscle. Animal research and musculoskeletal modeling (MSK) have detected heterogeneous fascicle behavior within muscle regions, and within fascicles. Understanding architectural behavior differences during muscle contractions may help to discern possible mechanisms behind proximal BFLH injuries. The purpose of our study was to assess the magnitude of shortening of the proximal and distal fascicles of the BFLH under a range of muscle activation levels under isometric conditions using ultrasound imaging (US). Thirteen healthy adults performed targeted sustained isometric contractions while US were taken of the entire BFLH. Measurements of fascicle lengths in both muscle regions were compared at 20%, 30%, 50%, and 67% MVIC. The results showed that while both regions shortened significantly with activation, the proximal fascicles were significantly longer, regardless of activation level (~38%), and shortened significantly more than the distal fascicles overall (~40%), and cumulatively at higher activation levels (30% and above). No significant strain differences were found between the two regions. These data suggest heterogeneous fascicle behavior exists in an absolute sense; however, differences in behavior are eliminated when normalized (strain). Coupled with MSK literature, the absence of regional fascicle strain differences in this study may indicate strain heterogeneity is not detectable at the whole fascicle level. Further knowledge of this commonly strained muscle?s regional behavior during dynamic movements could provide evidence of proximal hamstring strain predisposition.     

Motor unit composition has little effect on the short-range stiffness of feline medial gastrocnemius muscle.

Studies on skinned fibers and single motor units have indicated that slow-twitch fibers are stiffer than fast-twitch fibers. This suggests that skeletal muscles with different motor unit compositions may have different short-range stiffness (SRS) properties. Furthermore, the natural recruitment of slow before fast motor units may result in an SRS-force profile that is different from electrical stimulation. However, muscle architecture and the mechanical properties of surrounding tissues also contribute to the net SRS of a muscle, and it remains unclear how these structural features each contribute to the SRS of a muscle. In this study, the SRS-force characteristics of cat medial gastrocnemius muscle were measured during natural activation using the crossed-extension reflex, which activates slow before fast motor units, and during electrical activation, in which all motor units are activated synchronously. Short, rapid, isovelocity stretches were applied using a linear puller to measure SRS across the range of muscle forces. Data were collected from eight animals. Although there was a trend toward greater stiffness during natural activation, this trend was small and not statistically significant across the population of animals tested. A simple model, in which the slow-twitch fibers were assumed to be 30% stiffer than the fast-twitch fibers, was used to simulate the experimental results. Experimental and simulated results show that motor unit composition or firing rate has little effect on the SRS property of the cat MG muscle, suggesting that architectural features may be the primary determinant of SRS.     

Nonuniform activation of the agonist muscle does not covary with index finger acceleration in old adults.

This study examined the patterns of activation in the superficial and deep parts of the first dorsal interosseus muscle and in the antagonist muscle, second palmar interosseus, during postural tasks (position holding) and slow movements (position tracking) of the index finger performed by young and old adults. The position-tracking task involved the index finger lifting light loads (2.5, 10, and 35% of maximum) with shortening and lengthening contractions as steadily as possible. Steadiness was quantified in both tasks as the standard deviation of index finger acceleration. The fluctuations in acceleration during the two tasks were greater for the old subjects (62-72 yr) compared with young subjects (19-27 yr), especially with the lightest loads. The two groups of subjects activated the superficial and deep parts of first dorsal interosseus at similar intensities during the position-holding task, whereas the deep part was more active during the shortening and lengthening contractions of the position-tracking task. The nonuniform activation of first dorsal interosseus, therefore, was not associated with the difference in the standard deviation of acceleration between the young and old subjects. Furthermore, there was no association between the average level of coactivation by the antagonist muscle and the standard deviation of acceleration for either group of subjects across these tasks. Thus the greater variability in motor output exhibited by the older adults could not be explained by either the nonuniform activation of the agonist muscle or the average level of coactivation by the antagonist muscle.     

Regenerative activity of muscle tissue and thymus status in young and old mdx mice with muscle injury and wound xenoplasty

Responses of the skeletal muscle tissue and thymus to muscle injury (complete transection) and wound xenoplasty with the minced muscle tissue of newborn rats (tissue therapy) were studied in mdx mice aged 12–16 and 40–48 weeks. The muscle tissue of mdx mice has genetic defects causing chronic dystrophic processes in it. The muscle tissue of young mdx mice proved to retain a relatively high capacity for regeneration. Under conditions of tissue therapy of the wound, the formation of muscle fibers from muscle cells of the graft and active regeneration of muscle fibers in the recipient mice were observed, and no structural defects were detected in the thymus. The capacity of posttraumatic regeneration in old mdx mice was lower. The xenogenic graft was undergoing resorption, thereby suppressing regeneration of muscle fibers and causing further tissue destruction in the injured muscle. The thymus parenchyma was subject to degenerative changes such as the formation of gaps, hemorrhagic foci, and increased numbers of macrophages and mast cells.     

Neuromuscular organization of avian flight muscle: Morphology and contractile properties of motor units in the pectoralis (pars thoracicus) of pigeon (Columba livia)

We used acid digestion and glycogen depletion to determine fascicle organization, fiber morphology, and physiological and anatomical features of individual motor units of an in-series muscle, the pectoralis (pars thoracicus) of the pigeon (Columba livia). Most fascicles are attached at one end to connective tissue. Average fiber length in the four regions examined range from 42% to 66% of average fascicle length. More than 65% of fibers are blunt at one end of a fascicle and taper intrafascicularly. Fibers with blunt–blunt endings range from 13% to 31% of the population in different regions; taper–taper fibers range from 2% to 17%. Pigeon pectoralis fibers are distinguished histochemically into fast-twitch glycolytic (FG) and fast-twitch oxidative-glycolytic (FOG) populations. Three units composed of FG fibers (FG units) contract more quickly than three units composed of FOG fibers (FOG units) (range 31–37 vs 47–62 msec), produce more tetanic force (0.11–0.32 vs 0.02–0.05 N) and are more fatigable (<18% initial force vs >50% after repeated stimulation). Most motor units are confined to one of the four muscle regions. Territory of two FOG units is <30% of parent fascicle length. Territories of other units spanned parent fascicles; most fibers in these units do not extend the full fascicle length. Compared to FG units, FOG units have lower maximum innervation ratios and density indices (ratio of depleted/total FOG fibers in territory 8–14% vs 58–76% for FG units). These differences support the hypothesis that FG units are organized to produce substantial force and power for takeoff, landing and other ballistic movements whereas FOG units are suited for sustained flight when power requirements are reduced. Implications of findings for understanding the control of in-series muscles and the use of connective tissue elastic elements during wing movements are discussed. J.Morphol. 236:179–208, 1998. © 1998 Wiley-Liss, Inc.     

Lack of fiber type selectivity during reinnervation of neonatal rabbit soleus muscle

Fast and slow contracting fibers in neonatal mammalian skeletal muscle are each innervated in a highly specific manner by motor neurons of the corresponding type, even at an age when polyinnervation is widespread. Chemospecific recognition is a possible mechanism by which this pattern of innervation could be established. We have investigated this possibility by studying the degree of specificity during reinnervation of rabbit soleus muscle following nerve crush on Postnatal Day 1 or 4. We assayed fiber type composition by measuring the twitch rise times of motor units within 2 days of the onset of functional reinnervation (5-6 days after nerve crush). In contrast to the broad, bimodal distribution of single motor unit twitch rise times seen in normal muscles, motor units in reinnervated muscles yielded a narrower, unimodal distribution of rise times. Rise times of reinnervated units were intermediate to those of normal fast and slow units, suggesting that reinnervated units were composed of a mixture of fast and slow contracting fibers. An alternative possibility, that specific reinnervation was masked by contractile dedifferentiation of muscle fibers, was examined by maintaining a transmission blockade induced by botulinum toxin poisoning for an equivalent interval. Twitch rise times of treated motor units exhibited the distinctly bimodal distribution characteristic of normal muscles, suggesting that muscle fibers can retain contractile diversity during a transient period of denervation. We carried out computer simulations to estimate the amount of rise time diversity induced by varying degrees of specificity during reinnervation. Based on this analysis, we conclude that there is little if any selective reinnervation of muscle fiber types at the ages studied.     

Fiber content and myosin heavy chain composition of muscle spindles in aged human biceps brachii.

The present study investigated potential age-related changes in human muscle spindles with respect to the intrafusal fiber-type content and myosin heavy chain (MyHC) composition in biceps brachii muscle. The total number of intrafusal fibers per spindle decreased significantly with aging, due to a significant reduction in the number of nuclear chain fibers. Nuclear chain fibers in old spindles were short and some showed novel expression of MyHC alpha-cardiac. The expression of MyHC alpha-cardiac in bag1 and bag2 fibers was greatly decreased in the A region. The expression of slow MyHC was increased in nuclear bag1 fibers and that of fetal MyHC decreased in bag2 fibers whereas the patterns of distribution of the remaining MyHC isoforms were generally not affected by aging. We conclude that aging appears to have an important impact on muscle spindle composition. These changes in muscle spindle phenotype may reflect an age-related deterioration in sensory and motor innervation and are likely to have an impact in motor control in the elderly.     

Neuromuscular organization of avian flight muscle: architecture of single muscle fibres in muscle units of the pectoralis (pars thoracicus) of pigeon (Columba livia)

The M. pectoralis (pars thoracicus) of pigeons (Columba livia) is comprised of short muscle fibres that do not extend from muscle origin to insertion but overlap ''in-series''. Individual pectoralis motor units are limited in territory to a portion of muscle length and are comprised of either fast twitch, oxidative and glycolytic fibres (FOG) or fast twitch and glycolytic fibres (FG). FOG fibres make up 88 to 90% of the total muscle population and have a mean diameter one-half of that of the relatively large FG fibres. Here we report on the organization of individual fibres identified in six muscle units depleted of glycogen, three comprised of FOG fibres and three comprised of FG fibres. For each motor unit, fibre counts revealed unequal numbers of depleted fibres in different unit cross-sections. We traced individual fibres in one unit comprised of FOG fibres and a second comprised of FG fibres. Six fibres from a FOG unit (total length 15.45 mm) ranged from 10.11 to 11.82 mm in length and averaged (± s.d.) 10.74 ± 0.79 mm. All originated bluntly (en mass) from a fascicle near the proximal end of the muscle unit and all terminated intramuscularly. Five of these ended in a taper and one ended bluntly. Fibres coursed on average for 70% of the muscle unit length. Six fibres from a FG unit (total length 34.76 mm) ranged from 8.97 to 18.38 mm in length and averaged 15.32 ± 3.75 mm. All originated bluntly and terminated intramuscularly; one of these ended in a taper and five ended bluntly. Fibres coursed on average for 44% of the muscle unit length. Because fibres of individual muscle units do not extend the whole muscle unit territory, the effective cross-sectional area changes along the motor unit length. These non-uniformities in the distribution of fibres within a muscle unit emphasize that the functional interactions within and between motor units are complex.     

Motor unit firing rates and contractile properties in tibialis anterior of young and old men.

The effects of aging on motoneuron firing rates and muscle contractile properties were studied in tibialis anterior muscle by comparing results from six young (20.8 +/- 0.8 yr) and six old men (82.0 +/- 1.7 yr). For each subject, data were collected from repeated tests over a 2-wk period. Contractile tests included maximal voluntary contraction (MVC) with twitch interpolation and stimulated twitch contractions. The old men had 26% lower MVC torque (P < 0.01) than did the young men, but percent activation was not different (99.1 and 99.3%, respectively). Twitch contraction durations were 23% longer (P < 0.01) in the old compared with the young men. During a series of repeated brief steady-state contractions at 10, 25, 50, 75, and 100% MVC, motor unit firing rates were recorded. Results from approximately 950 motor unit trains in each subject group indicated that at all relative torque levels mean firing rates were 30-35% lower (P < 0.01) in the old subjects. Comparisons between young and old subjects' mean firing rates at each of 10%, 50%, and MVC torques and their corresponding mean twitch contraction duration yielded a range of moderate-to-high correlations (r = -0.67 to -0.84). That lower firing rates were matched to longer twitch contraction durations in the muscle of old men, and relatively higher firing rates were matched with shorter contraction times from the young men, indirectly supports the neuromuscular age-related remodeling principle.     

Mechanisms that contribute to differences in motor performance between young and old adults.

This paper examines the physiological mechanisms responsible for differences in the amplitude of force fluctuations between young and old adults. Because muscle force is a consequence of motor unit activity, the potential mechanisms include both motor unit properties and the behavior of motor unit populations. The force fluctuations, however, depend not only on the age of the individual but also on the muscle group performing the task, the type and intensity of the muscle contraction, and the physical activity status of the individual. Computer simulations and experimental findings performed on tasks that involved single agonist and antagonist muscles suggest that differences in force fluctuations are not attributable to motor unit twitch force, motor unit number, or nonuniform activation of the agonist muscle, but that they are influenced by the variability and common modulation of motor unit discharge in both the agonist and antagonist muscles. Because the amplitude of the force fluctuations does not vary linearly with muscle activation, these results suggest that multiple mechanisms contribute to the differences in force fluctuations between young and old adults, although the boundary conditions for each mechanism remain to be determined.     

Training-induced alterations in young and senescent rat diaphragm muscle.

The current study sought to examine the effects of chronic endurance treadmill running on oxidative capacity and capillary density in specific diaphragm muscle fiber types in young (5 mo) and senescent (greater than or equal to 23 mo) female Fischer 344 rats. Both young and senescent animals trained at approximately 75% of maximal O2 consumption for 1 h/day 5 days/wk for 10 wk. Plantaris citrate synthase activity was significantly increased (P less than 0.01) in both young and old trained groups. Densitometric analysis of succinate dehydrogenase (SDH) activity in diaphragm type I, IIa, and IIb muscle fibers was done using a computerized image-processing system. There were no age-related differences in SDH activity between the young and old groups for any of the fiber types. In addition, SDH activity was found to be significantly increased (P less than 0.05) in all three fiber types in both the young and senescent trained animals compared with their sedentary counterparts. Fiber size and capillary density did not differ between young and senescent rats, nor did exercise affect this measure. Each fiber, irrespective of type, had an average of approximately four capillaries in contact with it. However, type IIb fibers had a significantly lower capillary density per unit area than type I or IIa muscle fibers. The results indicate that the senescent costal diaphragm maintains its ability to adapt to an increased metabolic demand brought about by locomotor exercise. Of further interest is the finding that training adaptations occurred in all three fiber types, suggesting that increased work of breathing from moderate exercise leads to recruitment of all three fiber types.(ABSTRACT TRUNCATED AT 250 WORDS)     

Cardiac,skeletal muscle and serum irisin responses to with or without water exercise in young and old male rats: Cardiac muscle produces more irisin than skeletal muscle

Irisin converts white adipose tissue (WAT) into brown adipose tissue (BAT), as regulated by energy expenditure. The relationship between irisin concentrations after exercise in rats compared humans after exercise remains controversial. We therefore: (1) measured irisin expression in cardiac and skeletal muscle, liver, kidney, peripheral nerve sheath and skin tissues, as also serum irisin level in 10 week-old rats without exercise, and (2) measured tissue supernatant irisin levels in cardiac and skeletal muscle, and in response to exercise in young and old rats to establishing which tissues produced most irisin. Young (12 months) and old rats (24 months) with or without 10 min exercise (water floating) and healthy 10 week-old Sprague-Dawley rats without exercise were used. Irisin was absent from sections of skeletal muscle of unexercised rats, the only part being stained being the perimysium. In contrast, cardiac muscle tissue, peripheral myelin sheath, liver, kidneys, and skin dermis and hypodermis were strongly immunoreactivity. No irisin was seen in skeletal muscle of unexercised young and old rats, but a slight amount was detected after exercise. Strong immunoreactivity occurred in cardiac muscle of young and old rats with or without exercise, notably in pericardial connective tissue. Serum irisin increased after exercise, being higher in younger than older rats. Irisin in tissue supernatants (cardiac and skeletal muscle) was high with or without exercise. High supernatant irisin could come from connective tissues around skeletal muscle, especially nerve sheaths located within it. Skeletal muscle is probably not a main irisin source.     

Immunohistochemical studies on regulation of alternative splicing of fast skeletal muscle troponin T: non-uniform distribution of the exon x3 epitope in a single muscle fiber

Troponin T (TnT) isoforms of chicken fast skeletal muscle are classified into two types, breast-muscle-type (B-type) and leg-muscle-type (L-type) isoforms. These isoforms are produced from a single gene by differential alternative splicing of pre-mRNA. We investigated immunohistochemically the distribution of B-type TnT isoforms in chicken leg muscle (musculus biceps femoris), using anti-exon x3 that was raised against a synthetic peptide corresponding to exon x3 and recognized B-type, but not the L-type, TnT isoforms. Mosaic patterns of immunostaining showing locally different expression of B-type TnT isoforms in a single fiber were observed among fibers, and the non-uniform distribution of the isoforms was also detected in sectioned fibers and myofibrils from the muscle. The results indicated that regulation of pre-mRNA splicing of fast skeletal muscle TnT was different not only among the muscle fibers but also within a single fiber, suggesting that heterogeneous myonuclei in regulation of alternative splicings occur in a single muscle fiber.     

Effects of load carrying on metabolic cost and hindlimb muscle dynamics in guinea fowl (Numida meleagris).

The goal of this study was to test whether the contractile patterns of two major hindlimb extensors of guinea fowl are altered by load-carrying exercise. We hypothesized that changes in contractile pattern, specifically a decrease in muscle shortening velocity or enhanced stretch activation, would result in a reduction in locomotor energy cost relative to the load carried. We also anticipated that changes in kinematics would reflect underlying changes in muscle strain. Oxygen consumption, muscle activation intensity, and fascicle strain rate were measured over a range of speeds while animals ran unloaded vs. when they carried a trunk load equal to 22% of their body mass. Our results showed that loading produced no significant (P > 0.05) changes in kinematic patterns at any speed. In vivo muscle contractile strain patterns in the iliotibialis lateralis pars postacetabularis and the medial head of the gastrocnemius showed a significant increase in active stretch early in stance (P < 0.01), but muscle fascicle shortening velocity was not significantly affected by load carrying. The rate of oxygen consumption increased by 17% (P < 0.01) during loaded conditions, equivalent to 77% of the relative increase in mass. Additionally, relative increases in EMG intensity (quantified as mean spike amplitude) indicated less than proportional recruitment, consistent with force enhancement via stretch activation, in the proximal iliotibialis lateralis pars postacetabularis; however, a greater than proportional increase in the medial gastrocnemius was observed. As a result, when averaged for the two muscles, EMG intensity increased in direct proportion to the fractional increase in load carried.     

A micromechanical model of skeletal muscle to explore the effects of fiber and fascicle geometry

Computational models of muscle generally lump the material properties of connective tissue, muscle fibers, and muscle fascicles together into one constitutive relationship that assumes a transversely isotropic microstructure. These models do not take into account how variations in the microstructure of muscle affect its macroscopic material properties. The goal of this work was to develop micromechanical models of muscle to determine the effects of variations in muscle microstructure on the macroscopic constitutive behavior. We created micromechanical models at the fiber and fascicle levels based on histological cross-sections of two rabbit muscles, the rectus femoris (RF) and the soleus, to determine the effects of microstructure geometry (fiber and fascicle shapes) on the along-fiber shear modulus of muscle. The two fiber-level models predicted similar macroscopic shear moduli (within 13.5% difference); however, the two fascicle-level models predicted very different macroscopic shear moduli (up to 161% difference). We also used the micromechanical models to test the assumption that the macroscopic properties of muscle are transversely isotropic about the fiber (or fascicle) direction. The fiber-level models exhibited behavior consistent with the transverse isotropy assumption; however, the fascicle-level models exhibited transversely anisotropic behavior. Micromechanical models, combined with fiber and fiber bundle mechanical experiments, are needed to understand how normal or pathological variations in microstructure give rise to the observed macroscopic behavior of muscle.