Nonequilibrium thermodynamics and mitochondrial protein content predict insulin sensitivity and fuel selection during exercise in human skeletal muscle

Abstract

Introduction: Many investigators have attempted to define the molecular nature of changes responsible for insulin resistance in muscle, but a molecular approach may not consider the overall physiological context of muscle. Because the energetic state of ATP (ΔG_ATP) could affect the rate of insulin-stimulated, energy-consuming processes, the present study was undertaken to determine whether the thermodynamic state of skeletal muscle can partially explain insulin sensitivity and fuel selection independently of molecular changes.

Methods:³¹P-MRS was used with glucose clamps, exercise studies, muscle biopsies and proteomics to measure insulin sensitivity, thermodynamic variables, mitochondrial protein content, and aerobic capacity in 16 volunteers.

Results: After showing calibrated ³¹P-MRS measurements conformed to a linear electrical circuit model of muscle nonequilibrium thermodynamics, we used these measurements in multiple stepwise regression against rates of insulin-stimulated glucose disposal and fuel oxidation. Multiple linear regression analyses showed 53% of the variance in insulin sensitivity was explained by 1) VO_2max (p = 0.001) and the 2) slope of the relationship of ΔG_ATP with the rate of oxidative phosphorylation (p = 0.007). This slope represents conductance in the linear model (functional content of mitochondria). Mitochondrial protein content from proteomics was an independent predictor of fractional fat oxidation during mild exercise (R² = 0.55, p = 0.001).

Conclusion: Higher mitochondrial functional content is related to the ability of skeletal muscle to maintain a greater ΔG_ATP, which may lead to faster rates of insulin-stimulated processes. Mitochondrial protein content per se can explain fractional fat oxidation during mild exercise.

Introduction

Insulin resistance in skeletal muscle is a long-recognized hallmark feature of type 2 diabetes mellitus and obesity. Although type 2 diabetes generally only occurs when pancreatic beta cells are incapable of secreting enough insulin to maintain normoglycemia, this most often occurs in the setting of insulin resistance (DeFronzo et al., 2015). Numerous investigators using in vivo approaches in humans, rodent models, or in vitro experiments have attempted to define the molecular nature of the changes responsible for insulin resistance in skeletal muscle. Studies too numerous to cite completely have described molecular defects in insulin as action and fuel selection (Kahn and Cushman, 1985; Mandarino et al., 1986; Kelley and Mandarino, 1990; Cusi et al., 2000; Kelley and Mandarino, 2000; Patti et al., 2003). Despite learning a great deal of molecular physiology from these studies, this approach usually does not consider the overall physiological context of the skeletal muscle in which such defects may be found.

A different approach to understanding insulin action in skeletal muscle is to start from first principles, in this case from the observation that the energy state of resting and mildly exercising skeletal muscle is held far from equilibrium by the action of creatine kinase. Nonequilibrium thermodynamic theory arises from properties of thermodynamic equilibrium and can explain the mechanism and consequences of the extent to which skeletal muscle is able to maintain a high energy state when energy demand increases. In skeletal muscle, this energy state can be estimated as the energy made available by ATP hydrolysis, often termed ΔG_ATP. Mitochondria maintain this energy at rest and defend it from falling too steeply under conditions of energy demand, such as muscle contraction or other energy consuming processes, which might include insulin action. Insulin at physiologic concentrations measurably increases oxygen consumption in leg muscle, demonstrating an increase in energy demand (Kelley and Mandarino, 1990; Kelley et al., 1992; Mandarino et al., 1993; Mandarino et al., 1996).

It is well known that a higher ATP energy state is linearly related to power output of the contractile apparatus (Jeneson et al., 1995; Westerhoff et al., 1995). This relationship would apply not only to muscle contraction but also to the many processes that require movement within a cell, and we recently proposed that this principle could apply to transport processes like insulin-stimulated, kinesin-powered movement of GLUT4 vesicles along microtubules (Mandarino and Willis, 2023) due to similarities between kinesins and myosins (Vale and Milligan, 2000). Insulin sensitivity, therefore, could be a function of the ability of the cell to maintain the energy state of ATP, which depends mitochondrial functional capacity (Paganini et al., 1997; Glancy et al., 2008).

Thermodynamic variables can be measured in vivo in human skeletal muscle using calibrated ³¹P-MRS at rest or during recovery from muscle contraction (Kemp et al., 2007; Kemp et al., 2015; Sleigh et al., 2016). ³¹P-MRS allows estimation of the free energy of ATP noninvasively and continuously at rest, during exercise, and during recovery from exercise. ³¹P-MRS also can be used to estimate the rate of ATP synthesis (J_ATP) during recovery from exercise by monitoring the recovery of the concentration of phosphocreatine (Kemp et al., 2007; Kemp et al., 2015; Sleigh et al., 2016). PCr recovery from exercise follows a monoexponential time course, allowing estimation of the time constant parameter τ (Kemp et al., 2007; Kemp et al., 2015; Sleigh et al., 2016). However, ³¹P-MRS does not provide concentrations of energy phosphates unless spectra are calibrated to a standard, and unphosphorylated creatine is not visible. Because absolute concentrations of energy phosphate compounds such as PCr, ATP, ADP, as well as phosphate (Pi) are required to estimate thermodynamic variables, muscle biopsies at rest can be performed to assay [ATP], [PCr], and [Pi], allowing calculation of [ADP], the concentration of which is too low to be observed with MRS. Total creatine also can be assayed in muscle biopsies, providing all the absolute concentrations needed to calculate the energy of ATP (ΔG_ATP) and the rate of ATP synthesis (J_ATP). In skeletal muscle, mitochondria are responsible for most ATP production and maintain a high energy state of ATP. Increased ATP demand and the increased rate of ATP hydrolysis cause ΔG_ATP to fall in a nearly linear manner (Jeneson et al., 1995; Westerhoff et al., 1995). In response, J_ATP rises. The slope of this relationship depends on the cellular content of functional mitochondria (Dudley et al., 1987). Higher mitochondrial functional content enables better maintenance of ΔG_ATP (Paganini et al., 1997; Glancy et al., 2008).

A linear electrical circuit analog of the nearly linear kinetics of the relationship between the rate of oxidative phosphorylation and the energy from ATP has been used for over 30 years to simplify nonequilibrium thermodynamic considerations in skeletal muscle during non-steady state conditions when ΔG_ATP is changing (Meyer, 1989). There is theoretical and experimental evidence indicating that this model predicts linear relationships between 1) ΔG_ATP and J_ATP (the “force-flow” relationship), and 2) between the slope of the force-flow relationship and the functional mitochondrial content or “conductance” of the system. Establishment of these linear relationships is a key to validating this model in the context of the questions posed here. Conductance closely reflects mitochondrial functional content (Paganini et al., 1997; Glancy et al., 2008), while capacitance is proportional to the total creatine concentration, and affects the time course of changes during non-steady state conditions (Meyer, 1989). Finally, τ, the time constant of recovery of phosphocreatine after exercise, is the third variable needed to completely determine this model linear system. All of these parameters were estimated in the present study using a combination of ³¹P-MRS and biochemical assays in muscle biopsies. So, these studies can provide complete and detailed information about skeletal muscle nonequilibrium thermodynamics and the functional content of mitochondria needed to maintain a high energy state of the cell, which will drive energy-consuming processes more robustly (Jeneson et al., 1995). The reader interested in a more detailed theoretical treatment of calculations of thermodynamic parameters and derivation of the non-steady state treatment of nonequilibrium thermodynamics, including the meaning of the theoretical constructs of conductance and capacitance can refer to the Supplementary Material (Davis and Davis-Van Thienen, 1989; De Saedeleer and Marechal, 1984; Glancy et al., 2013; Golding et al., 1995; Harris et al., 1992; Hasselbach and Oetliker, 1983; Jeneson et al., 2009; Kushmerick, 1998; Masuda et al., 1990; Nicholls and Ferguson, 2013; Pate et al., 1998; Rottenberg, 1973; Walter et al., 1997; Westerhoff et al., 1981; Westerhoff and Van Damm, 1987; Willis et al., 2016).

We refer to mitochondrial protein abundance, determined using proteomics, as mitochondrial protein content, a biochemical measure of the abundance of all mitochondrial proteins, to distinguish it from mitochondrial functional content, which can differ from mitochondrial protein content to the extent that molecular changes in mitochondria may either increase or decrease the ability of a given mass of mitochondria to influence the rate of oxidative phosphorylation. The functional advantages of both higher mitochondrial protein and mitochondrial functional content can be seen in muscle adapted to endurance exercise, where there is a better ability to maintain the energy of ATP (Paganini et al., 1997), fuel selection that favors lipid over carbohydrate oxidation (Holloszy and Coyle, 1984; Helge et al., 2007), and better fatigue resistance (Fitts et al., 1975).

Many processes activated by insulin, such as GLUT4 translocation, require energy from ATP hydrolysis and therefore would have rates that are dependent on the energy of ATP. Compared to our knowledge of muscle contraction, we know much less about the impact of skeletal muscle nonequilibrium thermodynamics and mitochondrial functional content on skeletal muscle metabolism at rest or during insulin stimulation, especially in the setting of insulin resistance. Therefore, it is important to consider whether functional mitochondrial content and the capacity to maintain a higher ATP energy influences the rates of insulin-dependent processes, as recently conjectured (Mandarino and Willis, 2023). The primary aim of this investigation therefore was to use ³¹P-MRS, muscle biopsies, euglycemic clamps, and exercise tests in the context of the linear model to explore whether and to what extent nonequilibrium thermodynamic parameters and skeletal muscle mitochondrial functional content are related to insulin sensitivity and fuel selection (Mandarino and Willis, 2023). Finally, we measured VO_2max because of its well-described positive association with insulin sensitivity that likely is independent of mitochondrial content and more dependent on the cardiorespiratory system and delivery of oxygen to muscle and the capacity for muscle blood flow (Lundby et al., 2017). Linear models were used to determine which of these variables independently predict insulin sensitivity determined using a euglycemic clamp or fuel selection during mild exercise (Barakati et al., 2022).

Materials and methods

Participants and study design

Sixteen volunteers were screened with a history and physical examination, laboratory measurements, body composition (bioimpedance), an ECG, and a 75 g oral glucose tolerance test. Participants were not taking medications or supplements that affect glucose metabolism and were instructed to maintain their usual diet and not to engage in exercise 48 h before any testing. Volunteers were selected to have a wide range of body composition, aerobic capacity, and insulin sensitivity to reveal relationships among variables and detect the predictive value of thermodynamic, kinetic, and other variables. The design of the studies is given in Figure 1. All participants completed four study visits; a consent and screening examination, a euglycemic clamp experiment with muscle biopsy taken under resting, postabsorptive conditions, a graded cycle ergometry exercise test with indirect calorimetry and subsequent VO_2max determination, and a³¹P-MRS study with two periods of rest, leg extension exercise, and recovery from exercise, all conducted in the magnet. All studies were approved by the University of Arizona Institutional Review Board.

FIGURE 1

Euglycemic, hyperinsulinemic clamps, and muscle biopsies

All volunteers underwent a euglycemic, hyperinsulinemic clamp with a resting, basal muscle biopsy, starting at 7–8 a.m. after an overnight fast (Cusi et al., 2000). Isotopically labeled glucose (6,6-dideuteroglucose, Cambridge Laboratories) was used to trace glucose metabolism (Lefort et al., 2010), and steady state conditions were assumed for calculating the rates of glucose metabolism. Participants had a percutaneous needle biopsy of the vastus lateralis muscle taken under basal conditions 1 h before starting an insulin infusion at a rate of 80 mU·m⁻²·min⁻¹. Biopsies were frozen in liquid nitrogen before analysis (Cusi et al., 2000; Lefort et al., 2010). These biopsies were used to calibrate energy phosphates and total creatine.

Exercise testing

On another day, at least 1 week separated from the glucose clamp, gas exchange measurements were made at rest and during progressive cycle ergometer exercise using a Parvo Medics metabolic measurement system (TrueOne 2400, Parvo Medics, Salt Lake City, Utah) (Barakati et al., 2022). This system uses gas measurements from a mixing chamber with sampling four times per minute. After catheter placement for blood sampling, subjects rested on the cycle for 6 min and then exercised at power outputs of 15, 30, and 45 W for 6 minutes each. VO₂ and VCO₂ values at rest and during exercise were used to calculate ΔVO₂ and ΔVCO₂ during exercise. Delta values were used to calculate a Respiratory Exchange Ratio due to working muscle (ΔRER = ΔVCO₂/ΔVO₂). ΔRER estimates the oxidative metabolism of muscle performing mild exercise and we have validated this method previously for walking (Willis et al., 2005) and mild cycle ergometer exercise (Barakati et al., 2022). This was followed continuously by a ramp protocol to determine VO_2max.

Calibrated ³¹P magnetic resonance spectroscopy (MRS)

On a separate day, ³¹P- MRS was performed on a research-dedicated 3.0 T Skyra scanner (Siemens, Erlangen, Germany). For ³¹P data, a 60-mm-diameter coil (double tuned for ³¹P and ¹H) was positioned over the vastus lateralis. From the dimension of the coil and the size/geometry of a typical upper leg most of the signal in the unlocalized ³¹P-MRS measurements originates in the vastus (43). A ¹H-MR image was obtained (using the body coil of the 3.0 T system for transmission and the ³¹P/1H surface coil for detection) to assure correct coil position before acquiring ³¹P spectra. ³¹P data were acquired using free induction decay spectroscopy (spectral width = 2,000 Hz, 1,024 data points) (fully relaxed conditions, flip angle = 90^o, 20 s interpulse delay, 8 averages). Knee extension exercise was performed in the magnet using a weighted bag system (Sleigh et al., 2016). Bag weight was calibrated to individual exercise characteristics. Spectra were fitted in the time domain using a nonlinear least squares algorithm, jMRUI (Vanhamme et al., 1997) and fully relaxed spectra were used to quantify ATP, PCr, and Pi peak areas. Energy phosphate peak areas were calibrated to total adenylate content of muscle biopsies, assuming >99% of adenylates are present as ATP at the equilibrium established by creatine kinase. We ensured this by adding a phosphorylating system prior to luciferase detection of ATP (Lust et al., 1981; Passonneau and Lowry, 1993). Total creatine (TCr) concentration was assayed in muscle biopsies (Passonneau and Lowry, 1993). [ADP] was calculated using the creatine kinase equilibrium equation. ΔG_ATP values and J_ATP were calculated as described (Meyerspeer et al., 2020); individual time constants of recovery of PCr after exercise, τ, were derived by fitting recovery kinetics to a monoexponential function. This monoexponential function also yields the rate constant k (units of min⁻¹), which is used to calculate the instantaneous rate of oxidative phosphorylation (J_ATP), assuming J_ATP = the rate of resynthesis of PCr after exercise ceases. Thus, J_ATP (mM/min) = k ([PCR]_ss—[PCR]_t), where [PCr]_ss is the concentration of PCr at rest and [PCr]_t is the concentration of PCr at each timepoint t during recovery from exercise.

Proteomic assessment of mitochondrial protein content

Whole muscle (75 µg) lysates were resolved on 4%–20% gradient polyacrylamide gels, which were stained with Bio-Safe Coomassie G-250 Stain (#1610786; Biorad, Hercules, CA). Each lane was cut into 6 bands, which were excised and digested with trypsin. Samples were desalted and purified with C¹⁸ columns and evaporated and reconstituted in 30% acetonitrile/0.1% trifluoroacetic acid (TFA) immediately before mass spectrometry analysis. HPLC-ESI-MS/MS was performed in positive ion mode on a Thermo Scientific Orbitrap Fusion Lumos tribrid mass spectrometer fitted with an EASY-Spray Source (Thermo Scientific, San Jose, CA) (Parker et al., 2019). NanoLC was performed without a trap column using a Thermo Scientific UltiMate 3000 RSLCnano System with an EASY Spray C18 LC column (Thermo Scientific, 50 cm × 75 μm inner diameter, packed with PepMap RSLC C18 material, 2 μm, cat. #ES903); loading phase for 15 min at 0.300 μL/min; mobile phase, linear gradient of 1%–34% Buffer B in 119 min at 0.220 μL/min, followed by a step to 95% Buffer B over 4 min at 0.220 μL/min, hold 5 min at 0.250 μL/min, and then a step to 1% Buffer B over 5 min at 0.250 μL/min and a final hold for 10 min (total run 159 min); Buffer A = 0.1% FA/H₂O; Buffer B = 0.1% FA in 80% ACN. Solvents were liquid chromatography mass spectrometry grade. Spectra were acquired using XCalibur, version 2.3 (Thermo Scientific). A “top speed” data-dependent MS/MS analysis was performed. Dynamic exclusion was enabled with a repeat count of 1, a repeat duration of 30 s, and an exclusion duration of 60 s. A whole muscle proteome was quantified by label-free analysis using Progenesis QI software (Nonlinear Dynamics/Waters, Milford, MA). Normalized peak areas of all proteins annotated to mitochondria were summed into an index of total mitochondrial protein abundance. Individual normalized peak areas of mitochondrial and total proteins are given in Supplementary Table S1. Citrate synthase activity was also assayed in homogenates of muscle biopsies (Srere, 1969). The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD043032 and 10.6019/PXD043032.

Theoretical considerations

Please see Supplementary Material for theoretical background, assumptions, and calculations for thermodynamic parameters such as ΔG_ATP as well as the linear electrical circuit model of nonequilibrium thermodynamics proposed by Meyer to explain the linear relationship of τ with mitochondrial functional content (“conductance”) and capacitance represented by TCr, or total creatine (Meyer, 1989).

Analytical assays

Enrichment of deuterated glucose was determined by LC-MS (Lalia et al., 2016). Insulin was assayed using ELISA (Alpco, Salem, NH). Concentrations of total adenylates and creatine were assayed in muscle biopsies using modifications of published methods (Passonneau and Lowry, 1993) and expressed as a concentrations per volume muscle cell water.

Calculations and statistics

Rates of glucose turnover were calculated using steady state equations (Debodo et al., 1963). Statistical comparisons were performed using t-tests and stepwise multiple linear regression removing the least significant variable. Pearson’s correlation coefficient was used to assess relationships between two variables, using one or two-tailed tests as appropriate. Monoexponential curve fitting for estimating τ from PCr recovery kinetics was performed using the Solver add-in for Excel.

Results

Characteristics

Sixteen people participated in this study (4 men, 12 women). Participant characteristics are given in Table 1. Four participants had untreated type 2 diabetes as determined using HbA1c (8.1% ± 1.0%) or fasting plasma glucose (157 ± 37 mg/dL), four others had prediabetes (HbA1c 6.0% ± 1.0%), and the remainder had normal glucose tolerance (HbA1c 5.4% ± 0.2%). Endogenous glucose production in the postabsorptive state was 3.46 ± 0.13 mg/(kg-FFM^·min) and was completely suppressed during the insulin-infusion. Rates of basal and insulin-stimulated glucose metabolism are given in Table 1.

Calibrated 31P-MRS

Calibrated ³¹P-MRS was used to measure [ATP], [PCr], and [Pi] at rest, during 1.0 min of knee extension exercise, and then 4.0 min of recovery Table 2; Figure 2. This rest-exercise-recovery protocol was performed twice while the participant was positioned in the magnet. Results of biochemical analyses of total adenylates (ATP) and TCr in muscle biopsies were used to calibrate ³¹P-MRS energy phosphate concentrations on an individual level. These concentrations are given at rest and at end of exercise in Table 2. All ³¹P-MRS-determined peak areas of PCr, ATP, and Pi were calibrated to [ATP] determined using biopsies taken under basal, resting conditions. Total adenylates, a good estimate of resting [ATP] in muscle (Connett, 1989), averaged 6.7 ± 0.2 mM in the biopsies while total creatine (TCr) was 38.5 ± 7.3 mM. Free creatine was calculated as the difference between TCr and [PCr] calibrated to [ATP] in each biopsy. [ADP] was calculated from calibrated values of energy phosphates continuously using the creatine kinase equilibrium equation (Supplementary Material) (Veech et al., 1979), rising from about 18.8 μM at rest to nearly 38 µM after exercise. As expected, [PCr] fell significantly by the end of each exercise period with a nearly equimolar concomitant rise in [Pi]; [Pi] rose from 2.79 at rest to 7.05 mM at end of exercise. Resting ΔG_ATP fell about 1 kcal/mol on average, while pH rose minimally with exercise.

FIGURE 2

A typical example, in duplicate, of ³¹P-MRS-determined changes in [PCr] and [Pi] from rest to end-exercise and then recovery is shown in Figure 2. The cutout in Figure 2 also shows the monoexponential recovery from end-exercise to rest that was used to estimate τ, the time constant of recovery. Observed values were fit to a monoexponential equation to derive values of τ in duplicate for each participant as exemplified in Figure 2. The average value for τ was 30.9 ± 6 s, with a coefficient of variation of about 13.6%.

Conformation of data to the linear electrical circuit analog model of thermodynamics in healthy and insulin resistant skeletal muscle

During recovery from exercise estimates of the instantaneous values of both [PCr] andΔG_ATP were used to construct the expected near-linear relationship between ΔG_ATP and oxidative phosphorylation (J_ATP), (Meyer, 1988; Meyer, 1989). As shown in Figure 3A, averaged values were highly linear (r = 0.99). In addition, the individual force-flow slopes averaged 0.98 ± 0.01, with the lowest correlation coefficient being 0.95, emphasizing the linearity of these relationships on an average and individual level. Likewise (Figure 3B), the individual slopes of the ΔG_ATP vs. J_ATP force-flow relationship were a significant (p < 0.01) linear predictor of TCr/τ, the “conductance” term of the Meyer model (Eq. 16, Supplementary Material), that is a measure of mitochondrial functional content (Paganini et al., 1997; Glancy et al., 2008). Finally, the apparent capacitance derived as a theoretical construct of the linear model (calculated as the product of τ and the force-flow slope, or conductance, see Supplementary Material) was linearly related to biopsy-measured TCr (Figure 3C), also as predicted by the linear model of nonequilibrium thermodynamics (Eq. 15, Supplementary Material). Taken together, these results indicate that these measurements can be used validly to assess effects on insulin sensitivity or fuel oxidation in the context of the linear model of nonequilibrium thermodynamics in muscle.

FIGURE 3

Graded exercise, fuel selection, and VO_2max

Cycle ergometer exercise with indirect calorimetry consisted of a period of rest, transition to steady state exercise at 15, 30, and 45 W for 6 min each, followed continuously by a ramp protocol to determine VO_2maxFigure 4. The respiratory exchange ratio (RER) due to working muscle, or ΔRER (Ganley et al., 2011; Barakati et al., 2022), was calculated using changes in VO₂ and VCO₂ from rest to exercise (15 W) and at sequential power outputs of 30 and 45 W. The rise in whole-body RER with increasing power outputs (Figure 4A) was mirrored by the increase in muscle RER (Figure 4B) and a fall in the fraction of total muscle energy expenditure from fat (p < 0.01, Figure 4C). The average fraction of energy expenditure derived from fat (fractional fat oxidation during mild exercise) over all three power outputs was 49% ± 2%.

FIGURE 4

Multiple linear regression analysis reveals that mitochondrial functional content estimated by the slope of the force-flow relationship and VO_2max both are significant, independent predictors of insulin sensitivity

We used a stepwise multiple linear regression approach to assess independent predictors of insulin sensitivity, where the least significant independent variable was removed step-by-step for insulin-stimulated glucose disposal as the dependent variable. For the model in which insulin-stimulated Rd was the dependent variable, we used VO_2max, total mitochondrial protein abundance, fractional fat oxidation, τ, and the force-flow slope as independent variables in the stepwise analysis. For fractional fat oxidation as the dependent variable, we used VO_2max, total mitochondrial protein abundance, insulin-stimulated Rd, τ, and the force-flow slope as independent variables. The rationale for including these variables was that VO_2max principally reflects systemic processes including cardiorespiratory variables and blood hemoglobin concentrations (Lundby et al., 2017), mitochondrial protein abundance reflects the protein content of mitochondria, the force-flow slope of ΔG_ATP vs. J_ATP reflects the functional content of muscle mitochondria estimated from ³¹P-MRS, and τ is a simply-measured variable related to the rate of oxidative phosphorylation that does not require calibration of energy phosphate or creatine concentrations. Table 3 shows the correlation matrix for all these variables. We report results for the full regression model as well as the final model with stepwise removal of the least significant independent variable.

With regard to insulin-stimulated Rd as a measure of insulin action (Supplementary Table S2), the full model has an R² = 0.69 (p = 0.02), with only VO_2max (p = 0.0064) and the force-flow slope (p = 0.023) being statistically significant independent variables. Stepwise regression removing the least significant independent variable resulted in a final model with an R² = 0.54 (p = 0.0073) and again VO_2max (p = 0.012) and the force-flow slope (p = 0.048) being the only remaining statistically significant independent variables. The simple relationships of insulin-stimulated Rd with VO_2max and the force-flow slope are illustrated in Figures 5A, B, respectively. Best subsets multiple regression analysis confirmed this final stepwise model was the best fit.

FIGURE 5

Mitochondrial protein abundance is a significant independent predictor of fractional fat oxidation during mild exercise

As for fractional fat oxidation, the R² for the full model was 0.61 (p = 0.056), with no independent variable being statistically significant. However, stepwise multiple regression removing the least significant independent variable resulted in a final model with an R² = 0.55 (p = 0.001, Supplementary Table S2) and only mitochondrial protein abundance remaining as a statistically significant independent variable (p = 0.001). The simple Pearson’s correlation between these variables is shown in Figure 5C.

Relationship between mitochondrial protein abundance measured by proteomics and citrate synthase (CS) activity

Because CS activity is widely used as a measure of mitochondrial protein content, we compared this measure with the proteomics-determined global mitochondrial protein index. CS activity was 4.45 ± 1.59 μmol·min·g wet weight. CS activity was not correlated with the proteomics-based mitochondrial protein index (r = 0.12) but was highly correlated with CS protein abundance (normalized peak area) determined by proteomics (r = 0.84, p < 0.001). CS activity was not correlated with the fraction of lipid oxidized during mild exercise (r = −0.01). Interestingly, CS activity was correlated with the slope of the ΔG_ATP:J_ATP relationship.

Discussion

Because skeletal muscle insulin resistance often precedes the development of beta cell failure and type 2 diabetes mellitus, the underlying causes of insulin resistance have been the subject of intense interest. Many investigators have used a molecular approach to try to discover the biochemical and molecular factors that lead to insulin resistance. Such studies, many of which have been performed in humans, have provided valuable information about differences in insulin signal transduction, regulation of glucose transport, hexokinase and glucose uptake, glycogen synthesis, mitochondrial function, and gene expression and protein abundance measures (Kahn and Cushman, 1985; Mandarino et al., 1990; Mandarino et al., 1995; Shulman, 1999; Cusi et al., 2000; Kelley and Mandarino, 2000; Petersen and Shulman, 2002; Patti et al., 2003). Many other studies have used genetic manipulation of rodents, usually at a single gene level, to create models of altered insulin action (Wasserman and Ayala, 2005). All of these studies are predicated on the idea that by understanding the molecular details of a system, properties of the whole system can be deduced. However, there are alternate ways of understanding a system that start with fundamental principles of physiology, such as the laws of thermodynamics, and attempt to describe the system based on these principles rather than, or in accompaniment with, a description of biochemical or molecular processes. Both molecular and physiological approaches have value and should be interpreted together, yet to our knowledge, nearly all studies of the origins of skeletal muscle insulin resistance mainly have used a molecular or biochemical approach. In the present study, we asked the question of whether basic thermodynamic principles (nonequilibrium thermodynamics) can be used to explain variance in insulin sensitivity and changes in skeletal muscle fuel preference. To accomplish this, we used calibrated ³¹P-MRS to assess thermodynamic variables in skeletal muscle at rest and during recovery from mild exercise coupled with euglycemic, hyperinsulinemic clamps to assess insulin sensitivity and graded exercise tests to determine fuel selection in working muscle. The participants in the present study purposely were selected to provide a wide range of body composition, insulin action and exercise characteristics. We chose this approach in order to provide the best opportunity to determine relationships between these characteristics and the thermodynamic and fuel selection variables that we measured. Although selection of discrete groups of participants is commonly used in such studies, often the selection of discrete groups can mask the fact that these are continuous variables and participants fall along a continuous spectrum. Our approach in fact resulted in the revelation of novel relationships.

We used “calibrated” ³¹P-MRS to accomplish the purposes of this study (Kemp et al., 2007). Concentrations of energy phosphates that we calibrated to resting [ATP] assayed biochemically in skeletal muscle biopsies for each participant corresponded well to values previously reported using biopsies or ³¹P-MRS (Kemp et al., 2007; Lanza et al., 2011; Ripley et al., 2018). This also was the case for total creatine assayed in biopsies. Both resting [ATP] and total creatine were slightly lower than the often used “standard” assumptions of 8.2 mM [ATP] and 42.5 mM TCr in resting human muscle (Kemp et al., 2007; Kemp et al., 2015; Meyerspeer et al., 2020). What was most important regarding these findings, though, was that resting [ATP] varied from 5.33–8.2 mM and TCr ranged from 24.1–49.7 mM. Without calibration, this high level of variability would have had considerable impact on the calculation of individual values for energy phosphates, total creatine, J_ATP, and ΔG_ATP. Errors produced in these variables by using standard assumptions or other constant values likely would have led to an inability to determine whether thermodynamic considerations were related to insulin sensitivity or fuel selection. The importance of using calibrated values has been emphasized by Kemp and coworkers (2007, 2015) and more recently has been explored by Ripley and colleagues, who used a phosphate phantom and the water signal to calibrate energy phosphates and total creatine, respectively (Ripley et al., 2018).

We used a weighted bag system (Sleigh et al., 2016) to drive PCr down by 15%–25% and followed recovery of [PCr] to determine tau and estimate the rate of oxidative phosphorylation (J_ATP). Along with the corresponding energy state of the cell (ΔG_ATP). During the two exercise periods [ATP] remained nearly constant while [PCR] fell and [Pi] rose in a nearly equimolar fashion. The time course of [PCr] recovery, which is a function of the rate of oxidative phosphorylation (J_ATP), is characterized by the time constant τ of a monoexponential function (Meyerspeer et al., 2020). The value of τ was around 30 s in the current study, similar to what has been reported (Conley et al., 2000; Schrauwen-Hinderling et al., 2007a; Lanza et al., 2011; Cuthbertson et al., 2014), and was reproducible, with a coefficient of variation of around 13%. Since τ is estimated independently of calibrated values for energy phosphates or creatine, it is a commonly measured variable and provides some insight into muscle oxidative characteristics. However, in the absence of independent estimates of mitochondrial conductance or muscle capacitance (TCr), the value of τ potentially is misleading (see Supplementary Material for a more detailed discussion of this point). Thus, calibration of MRS signals is essential to assess the role nonequilibrium thermodynamics might play in insulin sensitivity or fuel selection.

To test the hypotheses regarding the relationship between the thermodynamic properties of muscle and insulin sensitivity or fuel selection, it is critical that our data conform to the linear model of nonequilibrium thermodynamics (Meyer, 1989). We showed that both the ΔG_ATP:J_ATP relationship and the relationship between the slope of that relationship and apparent conductance of the linear model (TCr/τ) were highly linear. In addition, the relationship between TCr measured in biopsies and apparent capacitance of the model also was linear and supports the concept that TCr is proportional to the capacitance of the electrical circuit analog model under non-steady state conditions of ΔG_ATP, such as during recovery from exercise.

We have conjectured that, on a purely theoretical basis, the ability of muscle to maintain a high ΔG_ATP in the face of a rise in energy demand could positively influence the rates of insulin-stimulated energy-consuming reactions (Mandarino and Willis, 2023). Supporting this notion, a rise in plasma insulin concentration increases oxygen consumption across the leg under resting conditions (Kelley et al., 1990; Kelley and Mandarino, 1990; Kelley et al., 1992; Kelley et al., 1993; Mandarino et al., 1996), indicating that insulin action induces a mild increase in energy demand. If lower mitochondrial content requires a lower ΔG_ATP to meet energy demand, energy demand might be met at the expense of slower rates of insulin-stimulated reactions (Mandarino and Willis, 2023). For example, force production and power output by contractile proteins decline with ΔG_ATP (Jeneson et al., 1995; Westerhoff et al., 1995). To the extent that motor proteins may share similar characteristics with proteins like kinesins (Vale and Milligan, 2000), we speculated that the rate of transit of GLUT4 vesicles along cytoskeletal tracks could be sensitive to ΔG_ATP. We used a simple linear nonequilibrium thermodynamic model to confirm that high conductance in the mitochondrial oxidative pathway predicts higher insulin sensitivity. This could indicate that higher mitochondrial functional content predicts insulin sensitivity. This finding is consistent with the concept that a higher energetic state of muscle may drive one or more insulin-stimulated energy-requiring processes at a faster rate, leading to higher insulin-stimulated glucose disposal during a euglycemic clamp. The current results do not provide direct proof of this idea, but the results provide a compelling argument to pursue additional exploration of this hypothesis. These findings appear to be related to the results of Koch and Britton (2022) that rats selected for running capability are healthier, leaner, more insulin sensitive, and long-lived than rats selected for poor running ability.

The second purpose of this study was to determine if the ability of mitochondria to maintain ΔG_ATP in skeletal muscle dictates fuel selection during mild exercise. A lower ΔG_ATP would be predicted to lead to higher carbohydrate oxidation at the expense of fat oxidation through elevations in [ADP], [Pi], and potentially [AMP] (Connett, 1989; Stanley and Connett, 1991). Lower ΔG_ATP during mild exercise would follow from a lower content of functional mitochondria that lessens the ability of muscle to defend ΔG_ATP during exercise or potentially any increase in energy demand (see Supplementary Material). The present results show that mitochondrial protein content is the best linear predictor of the fraction of fuel accounted for by lipid (r = 0.74). This is not direct evidence in favor of this mechanism but provides an argument for performing more mechanistic studies in the future to test this hypothesis.

A number of previous studies used ³¹P-MRS in vivo in humans to determine the relationship between oxidative phosphorylation defects and insulin sensitivity (Conley et al., 2000; Newcomer et al., 2001; Petersen et al., 2003; Petersen et al., 2004; Schrauwen-Hinderling et al., 2007a; Schrauwen-Hinderling et al., 2007b; Szendroedi et al., 2007; Lanza et al., 2011; Ripley et al., 2018; Toledo et al., 2018). Several studies used the saturation transfer technique (Pi →ATP) to estimate rates of oxidative phosphorylation (Petersen et al., 2003; Petersen et al., 2004; Szendroedi et al., 2007; Kacerovsky-Bielesz et al., 2009). However, this technique does not accurately measure oxidative phosphorylation (Balaban and Koretsky, 2011; From and Ugurbil, 2011; Kemp and Brindle, 2012), leaving those findings open to reinterpretation. Other studies of this topic employing ³¹P-MRS techniques used assumptions, rather than measurements, of the concentrations of ATP and total creatine (TCr) in muscle (Schrauwen-Hinderling et al., 2007a; Schrauwen-Hinderling et al., 2007b). Although these studies, mainly examining recovery of phosphocreatine after exercise, have revealed interesting information regarding the relationships between oxidative phosphorylation and insulin sensitivity, the lack of calibration of energy phosphates or creatine concentrations prohibits their use in determining a role for nonequilibrium thermodynamics in insulin action or fuel selection. Two other studies have used calibrated MRS to understand oxidative phosphorylation changes in aging (Conley et al., 2000) and type 2 diabetes (Ripley et al., 2018). Ripley et al. (2018) used ³¹P-MRS calibrated with a phosphate phantom for energy phosphate concentrations, with quantification of creatine using the water signal, to calculate [ADP] in resting skeletal muscle from insulin resistant obese and type 2 diabetic patients finding no difference in [ADP] between obese and type 2 diabetic volunteers. Rather, these investigators found lower [PCr] to be correlated with several factors including mitochondrial density determined from biopsies. We confirmed this inverse correlation between [PCr] mitochondrial protein content determined using proteomics (r = −0.42, p < 0.10). The calibrated [ATP] in that study (6–7 mM) was similar to our biopsy-determined value of about 6.8 mM, somewhat below a commonly used “standard” assumption of 8.2 mM. The reasons for these differences in [ATP] are unclear. Unfortunately, neither Conley et al. (2000) nor Ripley et al. (2018) put their findings into the context of nonequilibrium thermodynamics. Nevertheless, the results of Ripley et al. (2018) provide the valuable demonstration that to study mitochondrial function and content in vivo in human muscle requires validly calibrated values for concentrations of energy phosphates and total creatine. A unique aspect of the present study is that we used calibrated ³¹P-MRS to study the relationship of metabolic flows to thermodynamic forces and found that functional mitochondrial content and efficient oxidative energy transfer predict insulin sensitivity.

We also used a new label-free proteomics approach to calculate an index of mitochondrial protein content using lysates of whole muscle protein. This simple index was taken to be the sum of all individual abundance levels associated with proteins assigned to mitochondria. Nearly 500 mitochondrial proteins out of a total of over 4000 total proteins were quantified in this analysis. The use of hundreds rather than a single mitochondrial protein as an index of mitochondrial protein content should provide a more accurate measurement of mitochondrial content. However, mitochondrial functional content, measured as conductance in the electrical circuit linear model of nonequilibrium thermodynamics may be a better indicator of mitochondrial function than mere mitochondrial protein content with respect to insulin sensitivity. Even so, mitochondrial protein content was a good predictor of fat oxidation during mild exercise and citrate synthase activity was related to the apparent conductance of the model we used.

In conclusion, the results of the present study show that insulin resistance and a preference for carbohydrate oxidation in skeletal muscle can be explained by a lower free energy state that likely results from lower mitochondrial functional and protein content, as recently conjectured (Mandarino and Willis, 2023). The latter, in turn, may be a consequence of lower physical activity, leading to reduced mitochondrial functional content, although this study does not directly address the question of chronic physical activity. A number of studies have used various aspects of ³¹P-MRS to assess mitochondrial function in healthy and insulin resistant humans, but to our knowledge, this is the first study to use calibrated ³¹P-MRS to assess the energy state of muscle at rest and during mild exercise in the context of insulin resistance in the context of a model of skeletal muscle thermodynamics. These results show that strategies to raise skeletal muscle mitochondrial functional content, whether by increased physical activity, new pharmacological approaches, or even mitochondrial transplant (Guariento et al., 2021) likely would improve both insulin sensitivity and lipid oxidation and potentially reduce the risk of type 2 diabetes and cardiovascular disease.

References

• 1

  BalabanR. S.KoretskyA. P. (2011). Interpretation of ³¹P NMR saturation transfer experiments: What you can't see might confuse you. Focus on "standard magnetic resonance-based measurements of the Pi→ATP rate do not index the rate of oxidative phosphorylation in cardiac and skeletal muscles". Am. J. Physiol. Cell Physiol.301, C12–C15. 10.1152/ajpcell.00100.2011

  • CrossRef
  • Google Scholar
  • View reference in article

• 2

  BarakatiN.Zapata BustosR.ColettaD. K.LanglaisP. R.KohlerL. N.LuoM.et al (2022). Fuel selection in skeletal muscle exercising at low intensity; reliance on carbohydrate in very sedentary individuals. Metab. Syndr. Relat. Disord.21, 16–24. 10.1089/met.2022.0062

  • CrossRef
  • Google Scholar
  • View reference in article

• 3

  ConleyK. E.JubriasS. A.EsselmanP. C. (2000). Oxidative capacity and ageing in human muscle. J. Physiol.526 Pt 1, 203–210. 10.1111/j.1469-7793.2000.t01-1-00203.x

  • CrossRef
  • Google Scholar
  • View reference in article

• 4

  ConnettR. J. (1989). In vivo control of phosphofructokinase: System models suggest new experimental protocols. Am. J. Physiol.257, R878–R888. 10.1152/ajpregu.1989.257.4.R878

  • CrossRef
  • Google Scholar
  • View reference in article

• 5

  CusiK.MaezonoK.OsmanA.PendergrassM.PattiM. E.PratipanawatrT.et al (2000). Insulin resistance differentially affects the PI 3-kinase- and MAP kinase-mediated signaling in human muscle. J. Clin. Invest.105, 311–320. 10.1172/JCI7535

  • CrossRef
  • Google Scholar
  • View reference in article

• 6

  CuthbertsonD. J.IrwinA.SprungV. S.JonesH.PughC. J.DaousiC.et al (2014). Ectopic lipid storage in non-alcoholic fatty liver disease is not mediated by impaired mitochondrial oxidative capacity in skeletal muscle. Clin. Sci. (Lond)127, 655–663. 10.1042/CS20130404

  • CrossRef
  • Google Scholar
  • View reference in article

• 7

  DavisE. J.Davis-Van ThienenW. I. (1989). Force-flow and back-pressure relationships in mitochondrial energy transduction: An examination of extended state 3-state 4 transitions. Arch. Biochem. Biophys.275, 449–458. 10.1016/0003-9861(89)90391-3

  • CrossRef
  • Google Scholar
  • View reference in article

• 8

  De SaedeleerM.MarechalG. (1984). Chemical energy usage during isometric twitches of frog sartorius muscle intoxicated with an isomer of creatine, beta-guanidinopropionate. Pflugers Arch.402, 185–189. 10.1007/BF00583333

  • CrossRef
  • Google Scholar
  • View reference in article

• 9

  DebodoR. C.SteeleR.AltszulerN.DunnA.BishopJ. S. (1963). On the hormonal regulation of carbohydrate metabolism; studies with C14 glucose. Recent Prog. Horm. Res.19, 445–488.

  • Google Scholar
  • View reference in article

• 10

  DeFronzoR.FerraniniE.ZimmetP.AlbertiG. M. M. (2015). “Pathogenesis of type 2 diabetes mellitus,” in International textbook of diabetes mellitus. Editor DeFronzoR. (Chichester, UK: John Wiley and Sons).

  • Google Scholar
  • View reference in article

• 11

  DudleyG. A.TullsonP. C.TerjungR. L. (1987). Influence of mitochondrial content on the sensitivity of respiratory control. J. Biol. Chem.262, 9109–9114. 10.1016/s0021-9258(18)48054-4

  • CrossRef
  • Google Scholar
  • View reference in article

• 12

  FittsR. H.BoothF. W.WinderW. W.HolloszyJ. O. (1975). Skeletal muscle respiratory capacity, endurance, and glycogen utilization. Am. J. Physiol.228, 1029–1033. 10.1152/ajplegacy.1975.228.4.1029

  • CrossRef
  • Google Scholar
  • View reference in article

• 13

  FromA. H.UgurbilK. (2011). Standard magnetic resonance-based measurements of the Pi→ATP rate do not index the rate of oxidative phosphorylation in cardiac and skeletal muscles. Am. J. Physiol. Cell Physiol.301, C1–C11. 10.1152/ajpcell.00345.2010

  • CrossRef
  • Google Scholar
  • View reference in article

• 14

  GanleyK. J.StockA.HermanR. M.SantelloM.WillisW. T. (2011). Fuel oxidation at the walk-to-run-transition in humans. Metabolism60, 609–616. 10.1016/j.metabol.2010.06.007

  • CrossRef
  • Google Scholar
  • View reference in article

• 15

  GlancyB.BarstowT.WillisW. T. (2008). Linear relation between time constant of oxygen uptake kinetics, total creatine, and mitochondrial content in vitro. Am. J. Physiol. Cell Physiol.294, C79–C87. 10.1152/ajpcell.00138.2007

  • CrossRef
  • Google Scholar
  • View reference in article

• 16

  GlancyB.WillisW. T.ChessD. J.BalabanR. S. (2013). Effect of calcium on the oxidative phosphorylation cascade in skeletal muscle mitochondria. Biochemistry52, 2793–2809. 10.1021/bi3015983

  • CrossRef
  • Google Scholar
  • View reference in article

• 17

  GoldingE. M.TeagueW. E.Jr.DobsonG. P. (1995). Adjustment of K' to varying pH and pMg for the creatine kinase, adenylate kinase and ATP hydrolysis equilibria permitting quantitative bioenergetic assessment. J. Exp. Biol.198, 1775–1782. 10.1242/jeb.198.8.1775

  • CrossRef
  • Google Scholar
  • View reference in article

• 18

  GuarientoA.PiekarskiB. L.DoulamisI. P.BlitzerD.FerraroA. M.HarrildD. M.et al (2021). Autologous mitochondrial transplantation for cardiogenic shock in pediatric patients following ischemia-reperfusion injury. J. Thorac. Cardiovasc Surg.162, 992–1001. 10.1016/j.jtcvs.2020.10.151

  • CrossRef
  • Google Scholar
  • View reference in article

• 19

  HarrisR. C.SoderlundK.HultmanE. (1992). Elevation of creatine in resting and exercised muscle of normal subjects by creatine supplementation. Clin. Sci. (Lond)83, 367–374. 10.1042/cs0830367

  • CrossRef
  • Google Scholar
  • View reference in article

• 20

  HasselbachW.OetlikerH. (1983). Energetics and electrogenicity of the sarcoplasmic reticulum calcium pump. Annu. Rev. Physiol.45, 325–339. 10.1146/annurev.ph.45.030183.001545

  • CrossRef
  • Google Scholar
  • View reference in article

• 21

  HelgeJ. W.StallknechtB.RichterE. A.GalboH.KiensB. (2007). Muscle metabolism during graded quadriceps exercise in man. J. Physiol.581, 1247–1258. 10.1113/jphysiol.2007.128348

  • CrossRef
  • Google Scholar
  • View reference in article

• 22

  HolloszyJ. O.CoyleE. F. (1984). Adaptations of skeletal muscle to endurance exercise and their metabolic consequences. J. Appl. Physiol. Respir. Environ. Exerc Physiol.56, 831–838. 10.1152/jappl.1984.56.4.831

  • CrossRef
  • Google Scholar

• 23

  JenesonJ. A.WesterhoffH. V.BrownT. R.Van EchteldC. J.BergerR. (1995). Quasi-linear relationship between Gibbs free energy of ATP hydrolysis and power output in human forearm muscle. Am. J. Physiol.268, C1474–C1484. 10.1152/ajpcell.1995.268.6.C1474

  • CrossRef
  • Google Scholar

• 24

  JenesonJ. A.SchmitzJ. P.Van Den BroekN. M.Van RielN. A.HilbersP. A.NicolayK.et al (2009). Magnitude and control of mitochondrial sensitivity to ADP. Am. J. Physiol. Endocrinol. Metab.297, E774–E784. 10.1152/ajpendo.00370.2009

  • CrossRef
  • Google Scholar

• 25

  Kacerovsky-BieleszG.ChmelikM.LingC.PokanR.SzendroediJ.FarukuoyeM.et al (2009). Short-term exercise training does not stimulate skeletal muscle ATP synthesis in relatives of humans with type 2 diabetes. Diabetes58, 1333–1341. 10.2337/db08-1240

  • CrossRef
  • Google Scholar
  • View reference in article

• 26

  KahnB. B.CushmanS. W. (1985). Subcellular translocation of glucose transporters: Role in insulin action and its perturbation in altered metabolic states. Diabetes Metab. Rev.1, 203–227. 10.1002/dmr.5610010301

  • CrossRef
  • Google Scholar
  • View reference in article

• 27

  KelleyD. E.MandarinoL. J. (1990). Hyperglycemia normalizes insulin-stimulated skeletal muscle glucose oxidation and storage in noninsulin-dependent diabetes mellitus. J. Clin. Invest.86, 1999–2007. 10.1172/JCI114935

  • CrossRef
  • Google Scholar

• 28

  KelleyD. E.MandarinoL. J. (2000). Fuel selection in human skeletal muscle in insulin resistance: A reexamination. Diabetes49, 677–683. 10.2337/diabetes.49.5.677

  • CrossRef
  • Google Scholar

• 29

  KelleyD. E.ReillyJ. P.VenemanT.MandarinoL. J. (1990). Effects of insulin on skeletal muscle glucose storage, oxidation, and glycolysis in humans. Am. J. Physiol.258, E923–E929. 10.1152/ajpendo.1990.258.6.E923

  • CrossRef
  • Google Scholar
  • View reference in article

• 30

  KelleyD. E.MokanM.MandarinoL. J. (1992). Intracellular defects in glucose metabolism in obese patients with NIDDM. Diabetes41, 698–706. 10.2337/diab.41.6.698

  • CrossRef
  • Google Scholar
  • View reference in article

• 31

  KelleyD. E.MokanM.SimoneauJ. A.MandarinoL. J. (1993). Interaction between glucose and free fatty acid metabolism in human skeletal muscle. J. Clin. Invest.92, 91–98. 10.1172/JCI116603

  • CrossRef
  • Google Scholar

• 32

  KempG. J.BrindleK. M. (2012). What do magnetic resonance-based measurements of Pi-->ATP flux tell us about skeletal muscle metabolism?Diabetes61, 1927–1934. 10.2337/db11-1725

  • CrossRef
  • Google Scholar
  • View reference in article

• 33

  KempG. J.MeyerspeerM.MoserE. (2007). Absolute quantification of phosphorus metabolite concentrations in human muscle in vivo by 31P MRS: A quantitative review. NMR Biomed.20, 555–565. 10.1002/nbm.1192

  • CrossRef
  • Google Scholar

• 34

  KempG. J.AhmadR. E.NicolayK.PrompersJ. J. (2015). Quantification of skeletal muscle mitochondrial function by 31P magnetic resonance spectroscopy techniques: A quantitative review. Acta Physiol. (Oxf)213, 107–144. 10.1111/apha.12307

  • CrossRef
  • Google Scholar

• 35

  KochL. G.BrittonS. L. (2022). Biology: Motion is function. Funct. (Oxf)3, zqac030. 10.1093/function/zqac030

  • CrossRef
  • Google Scholar

• 36

  KushmerickM. J. (1998). Energy balance in muscle activity: Simulations of ATPase coupled to oxidative phosphorylation and to creatine kinase. Comp. Biochem. Physiol. B Biochem. Mol. Biol.120, 109–123. 10.1016/s0305-0491(98)00026-1

  • CrossRef
  • Google Scholar

• 37

  LaliaA. Z.DasariS.JohnsonM. L.RobinsonM. M.KonopkaA. R.DistelmaierK.et al (2016). Predictors of whole-body insulin sensitivity across ages and adiposity in adult humans. J. Clin. Endocrinol. Metab.101, 626–634. 10.1210/jc.2015-2892

  • CrossRef
  • Google Scholar

• 38

  LanzaI. R.BhagraS.NairK. S.PortJ. D. (2011). Measurement of human skeletal muscle oxidative capacity by 31P-MR spectroscopy: A cross-validation with in vitro measurements. J. Magn. Reson Imaging34, 1143–1150. 10.1002/jmri.22733

  • CrossRef
  • Google Scholar

• 39

  LawsonJ. W.VeechR. L. (1979). Effects of pH and free Mg2+ on the Keq of the creatine kinase reaction and other phosphate hydrolyses and phosphate transfer reactions. J. Biol. Chem.254, 6528–6537. 10.1016/s0021-9258(18)50400-2

  • CrossRef
  • Google Scholar

• 40

  LefortN.GlancyB.BowenB.WillisW. T.BailowitzZ.De FilippisE. A.et al (2010). Increased reactive oxygen species production and lower abundance of complex I subunits and carnitine palmitoyltransferase 1B protein despite normal mitochondrial respiration in insulin-resistant human skeletal muscle. Diabetes59, 2444–2452. 10.2337/db10-0174

  • CrossRef
  • Google Scholar

• 41

  LundbyC.MonteroD.JoynerM. (2017). Biology of VO(2) max: Looking under the physiology lamp. Acta Physiol. (Oxf)220, 218–228. 10.1111/apha.12827

  • CrossRef
  • Google Scholar

• 42

  LustW. D.FeussnerG. K.BarbehennE. K.PassonneauJ. V. (1981). The enzymatic measurement of adenine nucleotides and P-creatine in picomole amounts. Anal. Biochem.110, 258–266. 10.1016/0003-2697(81)90144-5

  • CrossRef
  • Google Scholar

• 43

  MandarinoL. J.WillisW. T. (2023). Can non-equilibrium thermodynamics explain skeletal muscle insulin resistance due to low mitochondrial content?Lancet Diabetes Endocrinol.11, 149–151. 10.1016/S2213-8587(23)00032-3

  • CrossRef
  • Google Scholar

• 44

  MandarinoL. J.MadarZ.KoltermanO. G.BellJ. M.OlefskyJ. M. (1986). Adipocyte glycogen synthase and pyruvate dehydrogenase in obese and type II diabetic subjects. Am. J. Physiol.251, E489–E496. 10.1152/ajpendo.1986.251.4.E489

  • CrossRef
  • Google Scholar

• 45

  MandarinoL. J.ConsoliA.KelleyD. E.ReillyJ. J.NurjhanN. (1990). Fasting hyperglycemia normalizes oxidative and nonoxidative pathways of insulin-stimulated glucose metabolism in noninsulin-dependent diabetes mellitus. J. Clin. Endocrinol. Metab.71, 1544–1551. 10.1210/jcem-71-6-1544

  • CrossRef
  • Google Scholar

• 46

  MandarinoL. J.ConsoliA.JainA.KelleyD. E. (1993). Differential regulation of intracellular glucose metabolism by glucose and insulin in human muscle. Am. J. Physiol.265, E898–E905. 10.1152/ajpendo.1993.265.6.E898

  • CrossRef
  • Google Scholar

• 47

  MandarinoL. J.PrintzR. L.CusiK. A.KinchingtonP.O'dohertyR. M.OsawaH.et al (1995). Regulation of hexokinase II and glycogen synthase mRNA, protein, and activity in human muscle. Am. J. Physiol.269, E701–E708. 10.1152/ajpendo.1995.269.4.E701

  • CrossRef
  • Google Scholar

• 48

  MandarinoL. J.ConsoliA.JainA.KelleyD. E. (1996). Interaction of carbohydrate and fat fuels in human skeletal muscle: Impact of obesity and NIDDM. Am. J. Physiol.270, E463–E470. 10.1152/ajpendo.1996.270.3.E463

  • CrossRef
  • Google Scholar

• 49

  MasudaT.DobsonG. P.VeechR. L. (1990). The Gibbs-Donnan near-equilibrium system of heart. J. Biol. Chem.265, 20321–20334. 10.1016/s0021-9258(17)30507-0

  • CrossRef
  • Google Scholar

• 50

  MeyerR. A. (1988). A linear model of muscle respiration explains monoexponential phosphocreatine changes. Am. J. Physiol.254, C548–C553. 10.1152/ajpcell.1988.254.4.C548

  • CrossRef
  • Google Scholar

• 51

  MeyerR. A. (1989). Linear dependence of muscle phosphocreatine kinetics on total creatine content. Am. J. Physiol.257, C1149–C1157. 10.1152/ajpcell.1989.257.6.C1149

  • CrossRef
  • Google Scholar

• 52

  MeyerspeerM.BoeschC.CameronD.DezortovaM.ForbesS. C.HeerschapA.et al (2020). (31) P magnetic resonance spectroscopy in skeletal muscle: Experts' consensus recommendations. NMR Biomed.34, e4246. 10.1002/nbm.4246

  • CrossRef
  • Google Scholar

• 53

  NewcomerB. R.Larson-MeyerD. E.HunterG. R.WeinsierR. L. (2001). Skeletal muscle metabolism in overweight and post-overweight women: An isometric exercise study using (31)P magnetic resonance spectroscopy. Int. J. Obes. Relat. Metab. Disord.25, 1309–1315. 10.1038/sj.ijo.0801673

  • CrossRef
  • Google Scholar

• 54

  NichollsD. G.FergusonS. J. (2013). Bioenergetics. Boston: Elsevier.

  • Google Scholar

• 55

  PaganiniA. T.FoleyJ. M.MeyerR. A. (1997). Linear dependence of muscle phosphocreatine kinetics on oxidative capacity. Am. J. Physiol.272, C501–C510. 10.1152/ajpcell.1997.272.2.C501

  • CrossRef
  • Google Scholar

• 56

  ParkerS. S.KrantzJ.KwakE. A.BarkerN. K.DeerC. G.LeeN. Y.et al (2019). Insulin induces microtubule stabilization and regulates the microtubule plus-end tracking protein network in adipocytes. Mol. Cell Proteomics18, 1363–1381. 10.1074/mcp.RA119.001450

  • CrossRef
  • Google Scholar

• 57

  PassonneauJ. J.LowryO. H. (1993). Enzymatic analysis. A practical guide. Totowa, New Jersey: The Humana Press, Inc.

  • Google Scholar

• 58

  PateE.Franks-SkibaK.CookeR. (1998). Depletion of phosphate in active muscle fibers probes actomyosin states within the powerstroke. Biophys. J.74, 369–380. 10.1016/S0006-3495(98)77794-X

  • CrossRef
  • Google Scholar

• 59

  PattiM. E.ButteA. J.CrunkhornS.CusiK.BerriaR.KashyapS.et al (2003). Coordinated reduction of genes of oxidative metabolism in humans with insulin resistance and diabetes: Potential role of PGC1 and NRF1. Proc. Natl. Acad. Sci. U. S. A.100, 8466–8471. 10.1073/pnas.1032913100

  • CrossRef
  • Google Scholar

• 60

  PetersenK. F.ShulmanG. I. (2002). Pathogenesis of skeletal muscle insulin resistance in type 2 diabetes mellitus. Am. J. Cardiol.90, 11G–18G. 10.1016/s0002-9149(02)02554-7

  • CrossRef
  • Google Scholar

• 61

  PetersenK. F.BefroyD.DufourS.DziuraJ.AriyanC.RothmanD. L.et al (2003). Mitochondrial dysfunction in the elderly: Possible role in insulin resistance. Science300, 1140–1142. 10.1126/science.1082889

  • CrossRef
  • Google Scholar

• 62

  PetersenK. F.DufourS.BefroyD.GarciaR.ShulmanG. I. (2004). Impaired mitochondrial activity in the insulin-resistant offspring of patients with type 2 diabetes. N. Engl. J. Med.350, 664–671. 10.1056/NEJMoa031314

  • CrossRef
  • Google Scholar

• 63

  RipleyE. M.ClarkeG. D.HamidiV.MartinezR. A.SettlesF. D.SolisC.et al (2018). Reduced skeletal muscle phosphocreatine concentration in type 2 diabetic patients: A quantitative image-based phosphorus-31 MR spectroscopy study. Am. J. Physiol. Endocrinol. Metab.315, E229–E239. 10.1152/ajpendo.00426.2017

  • CrossRef
  • Google Scholar

• 64

  RottenbergH. (1973). The thermodynamic description of enzyme-catalyzed reactions. The linear relation between the reaction rate and the affinity. Biophys. J.13, 503–511. 10.1016/S0006-3495(73)86004-7

  • CrossRef
  • Google Scholar

• 65

  Schrauwen-HinderlingV. B.KooiM. E.HesselinkM. K.JenesonJ. A.BackesW. H.Van EchteldC. J.et al (2007a). Impaired in vivo mitochondrial function but similar intramyocellular lipid content in patients with type 2 diabetes mellitus and BMI-matched control subjects. Diabetologia50, 113–120. 10.1007/s00125-006-0475-1

  • CrossRef
  • Google Scholar

• 66

  Schrauwen-HinderlingV. B.RodenM.KooiM. E.HesselinkM. K.SchrauwenP. (2007b). Muscular mitochondrial dysfunction and type 2 diabetes mellitus. Curr. Opin. Clin. Nutr. Metab. Care10, 698–703. 10.1097/MCO.0b013e3282f0eca9

  • CrossRef
  • Google Scholar

• 67

  ShulmanG. I. (1999). Cellular mechanisms of insulin resistance in humans. Am. J. Cardiol.84, 3J–10J. 10.1016/s0002-9149(99)00350-1

  • CrossRef
  • Google Scholar

• 68

  SleighA.LupsonV.ThankamonyA.DungerD. B.SavageD. B.CarpenterT. A.et al (2016). Simple and effective exercise design for assessing in vivo mitochondrial function in clinical applications using (31)P magnetic resonance spectroscopy. Sci. Rep.6, 19057. 10.1038/srep19057

  • CrossRef
  • Google Scholar

• 69

  SrereP. (1969). Citrate synthase.

  • Google Scholar

• 70

  StanleyW. C.ConnettR. J. (1991). Regulation of muscle carbohydrate metabolism during exercise. FASEB J.5, 2155–2159. 10.1096/fasebj.5.8.1827082

  • CrossRef
  • Google Scholar

• 71

  SzendroediJ.SchmidA. I.ChmelikM.TothC.BrehmA.KrssakM.et al (2007). Muscle mitochondrial ATP synthesis and glucose transport/phosphorylation in type 2 diabetes. PLoS Med.4, e154. 10.1371/journal.pmed.0040154

  • CrossRef
  • Google Scholar

• 72

  ToledoF. G. S.JohannsenD. L.CovingtonJ. D.BajpeyiS.GoodpasterB.ConleyK. E.et al (2018). Impact of prolonged overfeeding on skeletal muscle mitochondria in healthy individuals. Diabetologia61, 466–475. 10.1007/s00125-017-4496-8

  • CrossRef
  • Google Scholar

• 73

  ValeR. D.MilliganR. A. (2000). The way things move: Looking under the hood of molecular motor proteins. Science288, 88–95. 10.1126/science.288.5463.88

  • CrossRef
  • Google Scholar

• 74

  VanhammeL.Van Den BoogaartA.Van HuffelS. (1997). Improved method for accurate and efficient quantification of MRS data with use of prior knowledge. J. Magn. Reson129, 35–43. 10.1006/jmre.1997.1244

  • CrossRef
  • Google Scholar

• 75

  VeechR. L.LawsonJ. W.CornellN. W.KrebsH. A. (1979). Cytosolic phosphorylation potential. J. Biol. Chem.254, 6538–6547. 10.1016/s0021-9258(18)50401-4

  • CrossRef
  • Google Scholar

• 76

  WalterG.VandenborneK.MccullyK. K.LeighJ. S. (1997). Noninvasive measurement of phosphocreatine recovery kinetics in single human muscles. Am. J. Physiol.272, C525–C534. 10.1152/ajpcell.1997.272.2.C525

  • CrossRef
  • Google Scholar

• 77

  WassermanD. H.AyalaJ. E. (2005). Interaction of physiological mechanisms in control of muscle glucose uptake. Clin. Exp. Pharmacol. Physiol.32, 319–323. 10.1111/j.1440-1681.2005.04191.x

  • CrossRef
  • Google Scholar

• 78

  WesterhoffH. V.Van DammK. (1987). Thermodynamics and control of biological free energy transduction. Elsevier.

  • Google Scholar

• 79

  WesterhoffH. V.HellingwerfK. J.ArentsJ. C.ScholteB. J.Van DamK. (1981). Mosaic nonequilibrium thermodynamics describes biological energy transduction. Proc. Natl. Acad. Sci. U. S. A.78, 3554–3558. 10.1073/pnas.78.6.3554

  • CrossRef
  • Google Scholar

• 80

  WesterhoffH. V.Van EchteldC. J.JenesonJ. A. (1995). On the expected relationship between Gibbs energy of ATP hydrolysis and muscle performance. Biophys. Chem.54, 137–142. 10.1016/0301-4622(94)00129-8

  • CrossRef
  • Google Scholar

• 81

  WillisW. T.GanleyK. J.HermanR. M. (2005). Fuel oxidation during human walking. Metabolism54, 793–799. 10.1016/j.metabol.2005.01.024

  • CrossRef
  • Google Scholar

• 82

  WillisW. T.JackmanM. R.MesserJ. I.Kuzmiak-GlancyS.GlancyB. (2016). A simple hydraulic analog model of oxidative phosphorylation. Med. Sci. Sports Exerc48, 990–1000. 10.1249/MSS.0000000000000884

  • CrossRef
  • Google Scholar