Acute high-caffeine exposure increases autophagic flux and reduces protein synthesis in C2C12 skeletal myotubes.

Caffeine is a highly catabolic dietary stimulant. High caffeine concentrations (1-10 mM) have previously been shown to inhibit protein synthesis and increase protein degradation in various mammalian cell lines. The purpose of this study was to examine the effect of short-term caffeine exposure on cell signaling pathways that regulate protein metabolism in mammalian skeletal muscle cells. Fully differentiated C2C12 skeletal myotubes either received vehicle (DMSO) or 5 mM caffeine for 6 h. Our analysis revealed that caffeine promoted a 40% increase in autolysosome formation and a 25% increase in autophagic flux. In contrast, caffeine treatment did not significantly increase the expression of the skeletal muscle specific ubiquitin ligases MAFbx and MuRF1 or 20S proteasome activity. Caffeine treatment significantly reduced mTORC1 signaling, total protein synthesis and myotube diameter in a CaMKKß/AMPK-dependent manner. Further, caffeine promoted a CaMKII-dependent increase in myostatin mRNA expression that did not significantly contribute to the caffeine-dependent reduction in protein synthesis. Our results indicate that short-term caffeine exposure significantly reduced skeletal myotube diameter by increasing autophagic flux and promoting a CaMKKß/AMPK-dependent reduction in protein synthesis.

Inhibiting c-Jun N-terminal kinase partially attenuates caffeine-dependent cell death without alleviating the caffeine-induced reduction in mitochondrial respiration in C2C12 skeletal myotubes.

Caffeine is a widely consumed stimulant that has previously been shown to promote cytotoxic stress and even cell death in numerous mammalian cell lines. Thus far there is little information available regarding the toxicity of caffeine in skeletal muscle cells. Our preliminary data revealed that treating C2C12 myotubes with 5 mM caffeine for 6 h increased nuclear fragmentation and reduced basal and maximal oxygen consumption rate (OCR) in skeletal myotubes. The purpose of this study was to further elucidate the pathways by which caffeine increased cell death and reduced mitochondrial respiration. We specifically examined the role of c-Jun N-terminal kinase (JNK), which has previously been shown to simultaneously increase caspase-dependent cell death and reduce mitochondrial respiration in other mammalian cell lines. We found that caffeine promoted a dose-dependent increase in cell death in multinucleated myotubes but did not in mononucleated myoblasts. The addition of 10 µM Z-DEVD-FMK, a specific inhibitor of executioner caspases, completely inhibited caffeine-dependent cell death. Further, the addition of 400 µM dantrolene, a specific ryanodine receptor (RYR) inhibitor, prevented the caffeine-dependent increase in cell death and the reduction in basal and maximal OCR. We also discovered that caffeine treatment significantly increased the phosphorylation of JNK and that the addition of 30 µM SP600125 (JNKi), a specific JNK inhibitor, partially attenuated caffeine-induced cell death without preventing the caffeine-dependent reduction in basal and maximal OCR. Our results suggest that JNK partially mediates the increase in caspase-dependent cell death but does not contribute to reduced mitochondrial respiration in caffeine-treated skeletal muscle cells. We conclude that caffeine increased cell death and reduced mitochondrial respiration in a calcium-dependent manner by activating the RYR and promoting reticular calcium release.

The sensitivity of fast muscle contractile function to the major components of the sarcomere Ca(2+)-cycling system.

A reaction-diffusion model of a muscle sarcomere was developed to evaluate the sensitivity of force characteristics to diffusion and Ca(2+)-cycling components. The model compared well to experimental force measurements. Diffusion led to Ca(2+) gradients that enhanced maximal force and accelerated relaxation compared to when diffusion was infinitely fast. However, a modest increase in sarcomere length or radius led to a decrease in maximal force. Lowering the Ca(2+) release rate caused a lower maximal force, but increasing the rate led to only modest gains in maximal force while incurring much greater ATP costs associated with reuptake. Greater parvalbumin binding rates decreased maximal force but enhanced relaxation, and this effect was magnified when Ca(2+) uptake rates were lowered as may occur during fatigue. These results show a physiological set of parameters that lead to a functional sarcomere of known dimensions and contractile function, and the effects of parameter variation on muscle function.

Caffeine promotes autophagy in skeletal muscle cells by increasing the calcium-dependent activation of AMP-activated protein kinase.

Caffeine has been shown to promote calcium-dependent activation of AMP-activated protein kinase (AMPK) and AMPK-dependent glucose and fatty acid uptake in mammalian skeletal muscle. Though caffeine has been shown to promote autophagy in various mammalian cell lines it is unclear if caffeine-induced autophagy is related to the calcium-dependent activation of AMPK. The purpose of this study was to examine the role of calcium-dependent AMPK activation in regulating caffeine-induced autophagy in mammalian skeletal muscle cells. We discovered that the addition of the AMPK inhibitor Compound C could significantly reduce the expression of the autophagy marker microtubule-associated protein 1 light chain 3b-II (LC3b-II) and autophagic vesicle accumulation in caffeine treated skeletal muscle cells. Additional experiments using pharmacological inhibitors and RNA interference (RNAi) demonstrated that the calcium/calmodulin-activated protein kinases CaMKKß and CaMKII contributed to the AMPK-dependent expression of LC3b-II and autophagic vesicle accumulation in a caffeine dose-dependent manner. Our results indicate that in skeletal muscle cells caffeine increases autophagy by promoting the calcium-dependent activation of AMPK.

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.

Novel locomotor muscle design in extreme deep-diving whales.

Most marine mammals are hypothesized to routinely dive within their aerobic dive limit (ADL). Mammals that regularly perform deep, long-duration dives have locomotor muscles with elevated myoglobin concentrations that are composed of predominantly large, slow-twitch (Type I) fibers with low mitochondrial volume densities (V(mt)). These features contribute to extending ADL by increasing oxygen stores and decreasing metabolic rate. Recent tagging studies, however, have challenged the view that two groups of extreme deep-diving cetaceans dive within their ADLs. Beaked whales (including Ziphius cavirostris and Mesoplodon densirostris) routinely perform the deepest and longest average dives of any air-breathing vertebrate, and short-finned pilot whales (Globicephala macrorhynchus) perform high-speed sprints at depth. We investigated the locomotor muscle morphology and estimated total body oxygen stores of several species within these two groups of cetaceans to determine whether they (1) shared muscle design features with other deep divers and (2) performed dives within their calculated ADLs. Muscle of both cetaceans displayed high myoglobin concentrations and large fibers, as predicted, but novel fiber profiles for diving mammals. Beaked whales possessed a sprinter's fiber-type profile, composed of ~80% fast-twitch (Type II) fibers with low V(mt). Approximately one-third of the muscle fibers of short-finned pilot whales were slow-twitch, oxidative, glycolytic fibers, a rare fiber type for any mammal. The muscle morphology of beaked whales likely decreases the energetic cost of diving, while that of short-finned pilot whales supports high activity events. Calculated ADLs indicate that, at low metabolic rates, both beaked and short-finned pilot whales carry sufficient onboard oxygen to aerobically support their dives.

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.

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.

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 2°), the mitochondrial volume fraction (ε(mito)), and the diffusion coefficient of oxygen (DO 2) 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 2° ≥ 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)).

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

A mathematical model is developed to analyze the influence of chemical reaction and diffusion processes on the intracellular organization of mitochondria in skeletal muscle cells. The mathematical modeling approach uses a reaction-diffusion analysis of oxygen, ATP, and ADP involved in energy metabolism and mitochondrial function as governed by oxygen supply, volume fraction of mitochondria, and rates of reaction. Superimposed upon and coupled to the continuum species material balances is a cellular automata (CA) approach governing mitochondrial life cycles in response to the metabolic state of the cell. The effectiveness factor (η), defined as the ratio of reaction rate in the system with finite rates of diffusion to those in the absence of any diffusion limitation is used to assess diffusional constraints in muscle cells. The model shows the dramatic effects that the governing parameters have on the mitochondrial cycle of life and death and how these effects lead to changes in the distribution patterns of mitochondria observed experimentally. The model results showed good agreement with experimental results on mitochondrial distributions in mammalian muscle fibers. The η increases as the mitochondrial population is redistributed toward the fiber periphery in response to a decreased availability of oxygen. Modification of the CA parameters so that the mitochondrial lifecycle is more sensitive to the oxygen concentration caused larger mitochondrial shifts to the edge of the cell with smaller changes in oxygen concentration, and thus also lead to increased values of η. The present study shows that variation in oxygen supply, muscle activity and mitochondrial ATP supply influence the η and are the important parameters that can cause diffusion limitations. In order to prevent diffusion constraints, the cell resorts to shifts in their mitochondrial population towards the cell periphery, thus increasing η.

Reaction-diffusion constraints in living tissue: effectiveness factors in skeletal muscle design.

A mathematical model was developed to analyze the effects of intracellular diffusion of O(2) and high-energy phosphate metabolites on aerobic energy metabolism in skeletal muscle. We tested the hypotheses that in a range of muscle fibers from different species (1) aerobic metabolism was not diffusion limited and (2) that fibers had a combination of rate and fiber size that placed them at the brink of substantial diffusion limitation. A simplified chemical reaction rate law for mitochondrial oxidative phosphorylation was developed utilizing a published detailed model of isolated mitochondrial function. This rate law was then used as a boundary condition in a reaction-diffusion model that was further simplified using the volume averaging method and solved to determine the rates of oxidative phosphorylation as functions of the volume fraction of mitochondria, the size of the muscle cell, and the amount of oxygen delivered by the capillaries. The effectiveness factor, which is the ratio of reaction rate in the system with finite rates of diffusion to those in the absence of any diffusion limitations, defined the regions where intracellular diffusion of metabolites and O(2) may limit aerobic metabolism in both very small, highly oxidative fibers as well as in larger fibers with lower aerobic capacity. Comparison of model analysis with experimental data revealed that none of the fibers was strongly limited by diffusion, as expected. However, while some fibers were near substantial diffusion limitation, most were well within the domain of reaction control of aerobic metabolic rate. This may constitute a safety factor in muscle that provides a level of protection from diffusion constraints under conditions such as hypoxia.

Diffusional constraints on energy metabolism in skeletal muscle.

Aerobic metabolic flux depends on the diffusion of high-energy phosphate molecules (e.g., ATP and phosphocreatine) from the mitochondria to cellular ATPases, as well as the diffusion of other molecules (e.g., ADP, Pi) back to the mitochondria. Here, we develop an approach for evaluating the influence of intracellular metabolite diffusion on skeletal muscle aerobic metabolism through the application of the effectiveness factor (eta). This parameter provides an intuitive and informative means of quantifying the extent to which diffusion limits metabolic flux. We start with the classical approach assuming an infinite supply of substrate at the fiber boundary, and we expand this model to ultimately include nonlinear boundary and homogeneous reactions. Comparison of the model with experimental data from a wide range of skeletal muscle types reveals that most muscle fibers are not substantially limited by diffusion (eta close to unity), but many are on the brink of rather substantial diffusion limitation. This implies that intracellular metabolite diffusion does not dramatically limit aerobic metabolic flux in most fibers, but it likely plays a role in limiting the evolution of muscle fiber design and function.

Metabolic influences of fiber size in aerobic and anaerobic locomotor muscles of the blue crab, Callinectes sapidus.

Diameters of some white locomotor muscle fibers in the adult blue crab, Callinectes sapidus, exceed 500 microm whereas juvenile white fibers are <100 microm. It was hypothesized that aerobically dependent processes, such as metabolic recovery following burst contractions, will be significantly impeded in the large white fibers. In addition, dark aerobic fibers of adults, which rely on aerobic metabolism for both contraction and recovery, grow as large as the white fibers. These large aerobic fibers are subdivided, however, thus decreasing the effective diameter of each metabolic functional unit and enabling aerobic contraction. The two goals of this study were: (1) to characterize the development of subdivisions in the dark levator muscle fibers and (2) to monitor post-contractile metabolism as a function of fiber size in aerobic and anaerobic levator muscles. Dark levator muscle fibers from crabs ranging from <0.1 g to >190 g were examined with transmission electron microscopy to determine the density of mitochondria and subdivision diameters. Across all size classes, there was a constant mitochondrial fractional area (25% of the total subdivision area) and subdivision size (mean diameter of 36.5+/-2.7 microm). Thus, blue crab dark levator fibers are unusual in having metabolic functional units (subdivisions) that do not increase in size during development while the contractile functional units (fibers) grow hypertrophically. The body mass scaling of post-contractile lactate dynamics was monitored during recovery from anaerobic, burst exercise in white and dark muscle, and in hemolymph. There were no differences among size classes in lactate accumulation during exercise in either muscle. However, in white fibers from large crabs, lactate continued to increase after exercise, and lactate removal from tissues required a much longer period of time relative to smaller crabs. Differences in lactate removal among size classes were less pronounced in dark fibers, and post-contractile lactate accumulation was significantly higher in white than in dark fibers from large animals. These data suggest that the large white fibers invoke anaerobic metabolism following contraction to accelerate certain phases of metabolic recovery that otherwise would be overly slow. This implies that, in addition to the typical mass-specific decrease in oxidative capacity that accompanies increases in animal mass, aerobic metabolic processes become increasingly limited by surface area to volume and intracellular diffusion constraints in developing white muscle fibers.

Diffusional anisotropy is induced by subcellular barriers in skeletal muscle.

The time- and orientational-dependence of phosphocreatine (PCr) diffusion was measured using pulsed-field gradient nuclear magnetic resonance (PFG-NMR) as a means of non-invasively probing the intracellular diffusive barriers of skeletal muscle. Red and white skeletal muscle from fish was used because fish muscle cells are very large, which facilitates the examination of diffusional barriers in the intracellular environment, and because they have regions of very homogeneous fiber type. Fish were cold-acclimated (5 degrees C) to amplify the contrast between red and white fibers. Apparent diffusion coefficients, D, were measured axially, D(axially) and radially, D(radially), in small muscle strips over a time course ranging from 12 to 700 ms. Radial diffusion was strongly time dependent in both fiber types, and D decreased with time until a steady-state value was reached at a diffusion time approximately 100 ms. Diffusion was also highly anisotropic, with D(axially) being higher than D(radially) for all time points. The time scale over which changes in D(radially) occurred indicated that the observed anisotropy was not a result of interactions with the thick and thin filament lattice of actin and myosin or restriction within the cylindrical sarcolemma, as has been previously suggested. Rather, the sarcoplasmic reticulum (SR) and mitochondria appear to be the principal intracellular structures that inhibit mobility in an orientation-dependent manner. This work is the first example of diffusional anisotropy induced by readily identifiable intracellular structures.

Functional and Evolutionary Implications of the Distribution of Phosphagens in Primitive-Type Spermatozoa.

External fertilization is considered to be the primitive condition in metazoans. The spermatozoa of such organisms typically display a common primitive-type morphology that is present in a range of phyla. These spermatozoa are extremely polarized cells in that the site of ATP synthesis (mitochondria in midpiece) is located at large diffusion distances from the ATP sink (dynein ATPases in the flagellum). Spermatozoa of polychaetes, sipunculids, echiuroids, echinoderms, and tunicates contain the phosphagen creatine phosphate or express the corresponding phosphagen kinase creatine kinase (or both), even when other phosphagens/phosphagen kinases are present in somatic tissues and eggs. The selective expression of the creatine kinase system in these spermatozoa may be related to potential advantages in the cellular transport of energy. To evaluate this possibility, we compared the efficacy of the major phosphagen systems for cellular transport of energy. We used a facilitated diffusion model for spatial ATP buffering, taking into account relative differences in diffusivity and thermodynamic poise. At low ratios of [total phosphagen pool]/ [total adenine nucleotide pool] (CG+P/CAd ratio), creatine phosphate carried a higher fraction of total high-energy phosphate (J) than the other phosphagens. However, J values for all phosphagens were greater than 0.9, and these differences disappeared as the CG+P/CAd ratio was increased. Thus, the functional benefit of using CP, rather than other phosphagens, in energy transport is quite limited. The creatine kinase system became associated with primitive-type spermatozoa early in metazoan evolution. This association is not necessarily related to inherent advantages of this phosphagen system for buffering of ATP, but may be linked to historical events in the evolution of the cell phenotype.

Spatial resolution of transdermal water mobility using NMR microscopy.

High resolution images were obtained using high-field nuclear magnetic resonance microscopy of in vitro preparations of hydrated hairless rat skin. The major anatomical features observed were comparable to those seen by electron microscopy and include the stratum corneum, the viable epidermis, sebaceous glands, the cell layers surrounding hair follicles (the outer and inner root sheaths), and regions of subcutaneous fatty deposits. Calculated diffusion maps demonstrated that signal intensity is sufficient to obtain quantitative water mobility data from the viable epidermis and the hair follicle/sebaceous gland regions. Images from skin immersed in D2O clearly distinguish signal contributions that arise from fat from those which arise from water, and indicate that the calculated diffusion maps include only proton mobility from water in skin. These results may lead to further applications for using quantitative nuclear magnetic resonance microscopy for examining transdermal transport processes of in vitro skin preparations.