Influence of reaction and diffusion on spatial organization of mitochondria and effectiveness factors in skeletal muscle cell design

A mandibular advancement device attenuates the abnormal morphology and function of mitochondria from the genioglossus in obstructive sleep apnea-hypopnea syndrome rabbits

BACKGROUND

Obstructive sleep apnea hypopnea syndrome (OSAHS) is a serious and potentially life-threatening disease. Mandibular advancement device (MAD) has the characteristics of non-invasive, comfortable, portable and low-cost, making it the preferred treatment for mild-to-moderate OSAHS. Our previous studies found that abnormal contractility and fibre type distribution of the genioglossus could be caused by OSAHS. However, whether the mitochondria participate in these tissue changes is unclear. The effect of MAD treatment on the mitochondria of the genioglossus in OSAHS is also uncertain.

OBJECTIVE

To examine the morphology and function of mitochondria from the genioglossus in a rabbit model of obstructive sleep apnea-hypopnea syndrome (OSAHS), as well as these factors after insertion of a mandibular advancement device (MAD).

METHODS

Thirty male New Zealand white rabbits were randomised into three groups: control, OSAHS and MAD, with 10 rabbits in each group. Animals in Group OSAHS and Group MAD were induced to develop OSAHS by injection of gel into the submucosal muscular layer of the soft palate. The rabbits in Group MAD were fitted with a MAD. The animals in the control group were not treated. Further, polysomnography (PSG) and cone-beam computed tomography (CBCT) scan were used to measure MAD effectiveness. CBCT of the upper airway and PSG suggested that MAD was effective. Rabbits in the three groups were induced to sleep for 4-6 h per day for eight consecutive weeks. The genioglossus was harvested and detected by optical microscopy and transmission electron microscopy. The mitochondrial membrane potential was determined by laser confocal microscopy and flow cytometry. Mitochondrial complex I and IV activities were detected by mitochondrial complex assay kits.

RESULTS

OSAHS-like symptoms were induced successfully in Group OSAHS and rescued by MAD treatment. The relative values of the mitochondrial membrane potential, mitochondrial complex I activity and complex IV activity were significantly lower in Group OSAHS than in the control group; however, there was no significant difference between Group MAD and the control group. The OSAHS-induced injury and the dysfunctional mitochondria of the genioglossus muscle were reduced by MAD treatment.

CONCLUSION

Damaged mitochondrial structure and function were induced by OSAHS and could be attenuated by MAD treatment.

Evolved changes in the intracellular distribution and physiology of muscle mitochondria in high-altitude native deer mice

KEY POINTS

Mitochondrial function changes over time at high altitudes, but the potential benefits of these changes for hypoxia resistance remains unclear. We used high-altitude-adapted populations of deer mice, which exhibit enhanced aerobic performance in hypoxia, to examine whether changes in mitochondrial physiology or intracellular distribution in the muscle contribute to hypoxia resistance. Permeabilized muscle fibres from the gastrocnemius muscle had higher respiratory capacities in high-altitude mice than in low-altitude mice. Highlanders also had higher mitochondrial volume densities, due entirely to an enriched abundance of subsarcolemmal mitochondria, such that more mitochondria were situated near the cell membrane and adjacent to capillaries. There were several effects of hypoxia acclimation on mitochondrial function, some of which were population specific, but they differed from the evolved changes in high-altitude natives, which probably provide a better indication of adaptive traits that improve performance and hypoxia resistance at high altitudes.

ABSTRACT

High-altitude natives that have evolved to live in hypoxic environments provide a compelling system to understand how animals can overcome impairments in oxygen availability. We examined whether these include changes in mitochondrial physiology or intracellular distribution that contribute to hypoxia resistance in high-altitude deer mice (Peromyscus maniculatus). Mice from populations native to high and low altitudes were born and raised in captivity, and as adults were acclimated to normoxia or hypobaric hypoxia (equivalent to 4300 m elevation). We found that highlanders had higher respiratory capacities in the gastrocnemius (but not soleus) muscle than lowlanders (assessed using permeabilized fibres with single or multiple inputs to the electron transport system), due in large part to higher mitochondrial volume densities in the gastrocnemius. The latter was attributed to an increased abundance of subsarcolemmal (but not intermyofibrillar) mitochondria, such that more mitochondria were situated near the cell membrane and adjacent to capillaries. Hypoxia acclimation had no significant effect on these population differences, but it did increase mitochondrial cristae surface densities of mitochondria in both populations. Hypoxia acclimation also altered the physiology of isolated mitochondria by affecting respiratory capacities and cytochrome c oxidase activities in population-specific manners. Chronic hypoxia decreased the release of reactive oxygen species by isolated mitochondria in both populations. There were subtle differences in O₂ kinetics between populations, with highlanders exhibiting increased mitochondrial O₂ affinity or catalytic efficiency in some conditions. Our results suggest that evolved changes in mitochondrial physiology in high-altitude natives are distinct from the effects of hypoxia acclimation, and probably provide a better indication of adaptive traits that improve performance and hypoxia resistance at high altitudes.

Oxygen control of intracellular distribution of mitochondria in muscle fibers

Mitochondrial density in skeletal muscle fibers is governed by the demand for aerobic ATP production, but the heterogeneous distribution of these mitochondria appears to be governed by constraints associated with oxygen diffusion. We propose that each muscle fiber has an optimal mitochondrial distribution at which it attains a near maximal rate of ATP consumption (RATPase ) while mitochondria are exposed to a minimal oxygen concentration, thus minimizing reactive oxygen species (ROS) production. We developed a coupled reaction-diffusion/cellular automata (CA) mathematical model of mitochondrial function and considered four fiber types in mouse extensor digitorum longus (EDL) and soleus (SOL) muscle. The developed mathematical model uses a reaction-diffusion analysis of metabolites including oxygen, ATP, ADP, phosphate, and phosphocreatine (PCr) involved in energy metabolism and mitochondrial function. A CA approach governing mitochondrial life cycles in response to the metabolic state of the fiber was superimposed and coupled to the reaction-diffusion approach. The model results show the sensitivity of important model outputs such as the RATPase , effectiveness factor (η) and average oxygen concentration available at each mitochondrion to local oxygen concentration in the fibers through variation in the CA model parameter θdet , which defines the sensitivity of mitochondrial death to the oxygen concentration. The predicted optimal mitochondrial distributions matched previous experimental findings. Deviations from this optimal distribution corresponding to higher CA model parameter values (a more uniform mitochondrial distribution) lead to lower aerobic rates. In contrast, distributions corresponding to lower CA model parameter values (a more asymmetric distribution) lead to an increased exposure of mitochondria to oxygen, usually without substantial increases in aerobic rates, which would presumably result in increased ROS production and thus increased risks of cytotoxicity.

The formation and functional consequences of heterogeneous mitochondrial distributions in skeletal muscle

Diffusion plays a prominent role in governing both rates of aerobic metabolic fluxes and mitochondrial organization in muscle fibers. However, there is no mechanism to explain how the non-homogeneous mitochondrial distributions that are prevalent in skeletal muscle arise. We propose that spatially variable degradation with dependence on O(2) concentration, and spatially uniform signals for biogenesis, can account for observed distributions of mitochondria in a diversity of skeletal muscle. We used light and transmission electron microscopy and stereology to examine fiber size, capillarity and mitochondrial distribution in fish red and white muscle, fish white muscle that undergoes extreme hypertrophic growth, and four fiber types in mouse muscle. The observed distributions were compared with those generated using a coupled reaction-diffusion/cellular automata (CA) mathematical model of mitochondrial function. Reaction-diffusion analysis of metabolites such as oxygen, ATP, ADP and PCr involved in energy metabolism and mitochondrial function were considered. Coupled to the reaction-diffusion approach was a CA approach governing mitochondrial life cycles in response to the metabolic state of the fiber. The model results were consistent with the experimental observations and showed higher mitochondrial densities near the capillaries because of the sometimes steep gradients in oxygen. The present study found that selective removal of mitochondria in the presence of low prevailing local oxygen concentrations is likely the primary factor dictating the spatial heterogeneity of mitochondria in a diversity of fibers. The model results also suggest decreased diffusional constraints corresponding to the heterogeneous mitochondrial distribution assessed using the effectiveness factor, defined as the ratio of the reaction rate in the system with finite rates of diffusion to that in the absence of any diffusion limitation. Thus, the non-uniform distribution benefits the muscle fiber by increasing the energy status and increasing sustainable metabolic rates.

Sensitivity analysis of reaction-diffusion constraints in muscle energetics

Theoretical and experimental studies of aerobic metabolism on a wide range of skeletal muscle fibers have shown that while all fibers normally function within the reaction control regime, some fibers operate near the transition region where reaction control switches to diffusion control. Thus, the transition region between reaction and diffusion control may define the limits of muscle function, and analysis of factors that affect this transition is therefore needed. In order to assess the role of all important model parameters, a sensitivity analysis (SA) was performed to define the parameter space where muscle fibers transition from reaction to diffusion control. SA, performed on a previously developed reaction-diffusion model, shows that the maximum rate for the ATPase reaction (V(max,ATPase)), boundary oxygen concentration in the capillary supply (O ₂⁰), the mitochondrial volume fraction (ε(mito)), and the diffusion coefficient of oxygen (DO ₂) are the most sensitive parameters affecting this transition to diffusion control. It is demonstrated that fibers are not limited by diffusion for slow reactions (V(max,ATPase)  < 25 mM/min), high oxygen supply for the capillaries (O ₂⁰ ≥ 35 µM), and large amounts of mitochondria (ε(mito) ≥ 0.1). These conditions are applicable to muscle cells spanning a very broad range of animals. Within the diffusion-controlled region, the overall metabolic rate and ATP concentrations have much higher sensitivity to the diffusion coefficient of oxygen than to the diffusion coefficients of the other metabolites (ATP, ADP, P(i)).

Facilitated diffusion of myoglobin and creatine kinase and reaction-diffusion constraints of aerobic metabolism under steady-state conditions in skeletal muscle

The roles of creatine kinase (CK) and myoglobin (Mb) on steady-state facilitated diffusion and temporal buffering of ATP and oxygen, respectively, are assessed within the context of a reaction-diffusion model of muscle energetics. Comparison of the reaction-diffusion model with experimental data from a wide range of muscle fibers shows that the experimentally observed skeletal muscle fibers are generally not limited by diffusion, and the model further indicates that while some muscle fibers operate near the edge of diffusion limitation, no detectable effects of Mb and CK on the effectiveness factor, a measure of diffusion constraints, are observed under steady-state conditions. However, CK had a significant effect on average ATP concentration over a wide range of rates and length scales within the reaction limited regime. The facilitated diffusion functions of Mb and CK become observable in the model for larger size cells with low mitochondrial volume fraction and for low boundary O(2) concentration and high ATP demand, where the fibers may be limited by diffusion. From the transient analysis it may be concluded that CK primarily functions to temporally buffer ATP as opposed to facilitating diffusion while Mb has a small temporal buffering effect on oxygen but does not play any significant role in steady-state facilitated diffusion in skeletal muscle fibers under most physiologically relevant regions.