Oxidative ATP synthesis in skeletal muscle is controlled by substrate feedback

AW, Dash RK. Computational Modeling of Substrate-Dependent Mitochondrial Respiration and Bioenergetics in the Heart and Kidney Cortex and Outer Medulla

Integrated computational modeling provides a mechanistic and quantitative framework to characterize alterations in mitochondrial respiration and bioenergetics in response to different metabolic substrates in-silico. These alterations play critical roles in the pathogenesis of diseases affecting metabolically active organs such as heart and kidney. Therefore, the present study aimed to develop and validate thermodynamically constrained integrated computational models of mitochondrial respiration and bioenergetics in the heart and kidney cortex and outer medulla (OM). The models incorporated the kinetics of major biochemical reactions and transport processes as well as regulatory mechanisms in the mitochondria of these tissues. Intrinsic model parameters such as Michaelis-Menten constants were fixed at previously estimated values, while extrinsic model parameters such as maximal reaction and transport velocities were estimated separately for each tissue. This was achieved by fitting the model solutions to our recently published respirometry data measured in isolated rat heart and kidney cortex and OM mitochondria utilizing various NADH- and FADH2-linked metabolic substrates. The models were validated by predicting additional respirometry and bioenergetics data, which were not used for estimating the extrinsic model parameters. The models were able to predict tissue-specific and substrate-dependent mitochondrial emergent metabolic system properties such as redox states, enzyme and transporter fluxes, metabolite concentrations, membrane potential, and respiratory control index under diverse physiological and pathological conditions. The models were also able to quantitatively characterize differential regulations of NADH- and FADH2-linked metabolic pathways, which contribute differently toward regulations of oxidative phosphorylation and ATP synthesis in the heart and kidney cortex and OM mitochondria.

Transient changes to metabolic homeostasis initiate mitochondrial adaptation to endurance exercise

Endurance exercise is well established to increase mitochondrial content and function in skeletal muscle, a process termed mitochondrial biogenesis. Current understanding is that exercise initiates skeletal muscle mitochondrial remodeling via modulation of cellular nutrient, energetic and contractile stress pathways. These subtle changes in the cellular milieu are sensed by numerous transduction pathways that serve to initiate and coordinate an increase in mitochondrial gene transcription and translation. The result of these acute signaling events is the promotion of growth and assembly of mitochondria, coupled to a greater capacity for aerobic ATP provision in skeletal muscle. The aim of this review is to highlight the acute metabolic events induced by endurance exercise and the subsequent molecular pathways that sense this transient change in cellular homeostasis to drive mitochondrial adaptation and remodeling.

How does exercise affect energy metabolism? An in silico approach for cardiac muscle

We explored an in silico model of muscle energy metabolism and demonstrated its theoretical plausibility. Results indicate that energy metabolism triggered by activation can capture the muscle condition, rest, or exercise, and can respond accordingly adjusting the rates of their respiration and energy utilization for efficient use of the nutrients. Our study demonstrated during exercise higher respiratory activity causes a substantial increase in exergy release with an increase in exergy destruction, and entropy generation rate. The thermodynamic analysis showed that at the resting state when the exergy destruction rate was 0.66 W/kg and the respiratory metabolism energetic efficiency was 36% and exergetic efficiency was 32%; whereas, when the exergy destroyed was 1.24 W/kg, the energetic efficiency was 58% and exergetic efficiency was 50% during exercise. The efficiency results suggest the ability of the system to regulate itself in response to higher work demand and become more efficient in terms of converting energy coming from nutrients to useable energy when the circulating medium has sufficient energy precursor.

Computational Model of the Effect of Mitochondrial Dysfunction on Excitation–Contraction Coupling in Skeletal Muscle

It has become well established that mitochondria not only regulate myoplasmic calcium in skeletal muscle, but also use that calcium to stimulate oxidative phosphorylation (OXPHOS). While experimental approaches have allowed for imaging of mitochondrial calcium and membrane potentials in isolated fibers, capturing the role of mitochondria and the impact of mitochondrial impairments on excitation–contraction coupling (ECC) remains difficult to explore in intact muscle. Computational models have been widely used to examine the structure and function of skeletal muscle contraction; however, models of ECC to date lack communication between the myoplasm and mitochondria for regulating calcium and ATP during sustained contractions. To address this, a mathematical model of mitochondrial calcium handling and OXPHOS was integrated into a physiological model of ECC incorporating action potential propagation, calcium handling between the sarcoplasmic reticulum (SR) and the myoplasm, and crossbridge cycling. The model was used to examine the protective role of mitochondria during repeated stimulation and the impact of mitochondrial dysfunction on ECC resulting from progressive OXPHOS inhibition. Pathological myoplasmic calcium accumulation occurred through distinct mechanisms in the model in the case of either electron transport chain, F1F0 ATP synthase, or adenine nucleotide transporter impairments. To investigate the effect of each impairment on force, a model of calcium-stimulated apoptosis was utilized to capture dysfunction-induced reductions in muscle mass, driving whole muscle force loss. The model presented in this study can be used to examine the role of mitochondria in the regulation of calcium, ATP, and force generation during voluntary contraction.

Substrate- and Calcium-Dependent Differential Regulation of Mitochondrial Oxidative Phosphorylation and Energy Production in the Heart and Kidney

Mitochondrial dehydrogenases are differentially stimulated by Ca2+. Ca2+ has also diverse regulatory effects on mitochondrial transporters and other enzymes. However, the consequences of these regulatory effects on mitochondrial oxidative phosphorylation (OxPhos) and ATP production, and the dependencies of these consequences on respiratory substrates, have not been investigated between the kidney and heart despite the fact that kidney energy requirements are second only to those of the heart. Our objective was, therefore, to elucidate these relationships in isolated mitochondria from the kidney outer medulla (OM) and heart. ADP-induced mitochondrial respiration was measured at different CaCl2 concentrations in the presence of various respiratory substrates, including pyruvate + malate (PM), glutamate + malate (GM), alpha-ketoglutarate + malate (AM), palmitoyl-carnitine + malate (PCM), and succinate + rotenone (SUC + ROT). The results showed that, in both heart and OM mitochondria, and for most complex I substrates, Ca2+ effects are biphasic: small increases in Ca2+ concentration stimulated, while large increases inhibited mitochondrial respiration. Furthermore, significant differences in substrate- and Ca2+-dependent O2 utilization towards ATP production between heart and OM mitochondria were observed. With PM and PCM substrates, Ca2+ showed more prominent stimulatory effects in OM than in heart mitochondria, while with GM and AM substrates, Ca2+ had similar biphasic regulatory effects in both OM and heart mitochondria. In contrast, with complex II substrate SUC + ROT, only inhibitory effects on mitochondrial respiration was observed in both the heart and the OM. We conclude that the regulatory effects of Ca2+ on mitochondrial OxPhos and ATP synthesis are biphasic, substrate-dependent, and tissue-specific.

Quantitative analysis of mitochondrial ATP synthesis

We present a computational framework for analyzing and simulating mitochondrial ATP synthesis using basic thermodynamic and kinetic principles. The framework invokes detailed descriptions of the thermodynamic driving forces associated with the processes of the electron transport chain, mitochondrial ATP synthetase, and phosphate and adenine nucleotide transporters. Assembling models of these discrete processes into an integrated model of mitochondrial ATP synthesis, we illustrate how to analyze and simulate in vitro respirometry experiments and how models identified from in vitro experimental data effectively explain cardiac respiratory control in vivo. Computer codes for these analyses are embedded as Python scripts in a Jupyter Book to facilitate easy adoption and modification of the concepts developed here. This accessible framework may also prove useful in supporting educational applications. All source codes are available on at https://beards-lab.github.io/QAMAS_book/.

Mitochondria-localized AMPK responds to local energetics and contributes to exercise and energetic stress-induced mitophagy

Mitochondria form a complex, interconnected reticulum that is maintained through coordination among biogenesis, dynamic fission, and fusion and mitophagy, which are initiated in response to various cues to maintain energetic homeostasis. These cellular events, which make up mitochondrial quality control, act with remarkable spatial precision, but what governs such spatial specificity is poorly understood. Herein, we demonstrate that specific isoforms of the cellular bioenergetic sensor, 5' AMP-activated protein kinase (AMPKα1/α2/β2/γ1), are localized on the outer mitochondrial membrane, referred to as mitoAMPK, in various tissues in mice and humans. Activation of mitoAMPK varies across the reticulum in response to energetic stress, and inhibition of mitoAMPK activity attenuates exercise-induced mitophagy in skeletal muscle in vivo. Discovery of a mitochondrial pool of AMPK and its local importance for mitochondrial quality control underscores the complexity of sensing cellular energetics in vivo that has implications for targeting mitochondrial energetics for disease treatment.

Impaired aerobic capacity and premature fatigue preceding muscle weakness in the skeletal muscle Tfam-knockout mouse model

Mitochondrial diseases are genetic disorders that lead to impaired mitochondrial function, resulting in exercise intolerance and muscle weakness. In patients, muscle fatigue due to defects in mitochondrial oxidative capacities commonly precedes muscle weakness. In mice, deletion of the fast-twitch skeletal muscle-specific Tfam gene (Tfam KO) leads to a deficit in respiratory chain activity, severe muscle weakness and early death. Here, we performed a time-course study of mitochondrial and muscular dysfunctions in 11- and 14-week-old Tfam KO mice, i.e. before and when mice are about to enter the terminal stage, respectively. Although force in the unfatigued state was reduced in Tfam KO mice compared to control littermates (wild type) only at 14 weeks, during repeated submaximal contractions fatigue was faster at both ages. During fatiguing stimulation, total phosphocreatine breakdown was larger in Tfam KO muscle than in wild-type muscle at both ages, whereas phosphocreatine consumption was faster only at 14 weeks. In conclusion, the Tfam KO mouse model represents a reliable model of lethal mitochondrial myopathy in which impaired mitochondrial energy production and premature fatigue occur before muscle weakness and early death.

Summary: A time-course study of mitochondrial and muscular dysfunctions in a mouse model of mitochondrial myopathy reveals that decreased resistance to fatigue together with decreased oxidative capacities arise ahead of muscle weakness.

Effects of extracellular potassium on calcium handling and force generation in a model of excitation-contraction coupling in skeletal muscle

It is well-established that extracellular potassium (K_o⁺) accumulation reduces muscle fiber excitability, however the effects of K_o⁺ on the excitation-contraction coupling (ECC) pathway are less understood. In vivo and in vitro studies following fatiguing stimulation protocols are limited in their ability to capture the effects of K_o⁺ on force production in combination with other simultaneously changing factors. To address this, a computational model of ECC for slow and fast twitch muscle is presented to explore the relative contributions of excitability-induced and metabolic-induced changes in force generation in response to increasing [Formula: see text] . The model incorporates mechanisms previously unexplored in modelling studies, including the effects of extracellular calcium on excitability, calcium-dependent inhibition of calcium release, ATP-dependent ionic pumping, and the contribution of ATP hydrolysis to intracellular phosphate accumulation rate. The model was able to capture the frequency-dependent biphasic Force- [Formula: see text] response observed experimentally. Force potentiation for moderately elevated [Formula: see text] was driven by increased action potential duration, myoplasmic calcium potentiation, and phosphate accumulation rate, while attenuation of force at higher [Formula: see text] was due to action potential failure resulting in reduced calcium release. These results suggest that altered calcium release and phosphate accumulation work together with elevated K_o⁺ to affect force during sustained contractions.

Magnetic resonance spectroscopy reveals mitochondrial dysfunction in amyotrophic lateral sclerosis

Mitochondrial dysfunction is postulated to be central to amyotrophic lateral sclerosis (ALS) pathophysiology. Evidence comes primarily from disease models and conclusive data to support bioenergetic dysfunction in vivo in patients is currently lacking. This study is the first to assess mitochondrial dysfunction in brain and muscle in individuals living with ALS using 31P-magnetic resonance spectroscopy (MRS), the modality of choice to assess energy metabolism in vivo. We recruited 20 patients and 10 healthy age and gender-matched control subjects in this cross-sectional clinico-radiological study. 31P-MRS was acquired from cerebral motor regions and from tibialis anterior during rest and exercise. Bioenergetic parameter estimates were derived including: ATP, phosphocreatine, inorganic phosphate, adenosine diphosphate, Gibbs free energy of ATP hydrolysis (ΔGATP), phosphomonoesters, phosphodiesters, pH, free magnesium concentration, and muscle dynamic recovery constants. Linear regression was used to test for associations between brain data and clinical parameters (revised amyotrophic functional rating scale, slow vital capacity, and upper motor neuron score) and between muscle data and clinico-neurophysiological measures (motor unit number and size indices, force of contraction, and speed of walking). Evidence for primary dysfunction of mitochondrial oxidative phosphorylation was detected in the brainstem where ΔGATP and phosphocreatine were reduced. Alterations were also detected in skeletal muscle in patients where resting inorganic phosphate, pH, and phosphomonoesters were increased, whereas resting ΔGATP, magnesium, and dynamic phosphocreatine to inorganic phosphate recovery were decreased. Phosphocreatine in brainstem correlated with respiratory dysfunction and disability; in muscle, energy metabolites correlated with motor unit number index, muscle power, and speed of walking. This study provides in vivo evidence for bioenergetic dysfunction in ALS in brain and skeletal muscle, which appears clinically and electrophysiologically relevant. 31P-MRS represents a promising technique to assess the pathophysiology of mitochondrial function in vivo in ALS and a potential tool for future clinical trials targeting bioenergetic dysfunction.

Complex I is bypassed during high intensity exercise

Human muscles are tailored towards ATP synthesis. When exercising at high work rates muscles convert glucose to lactate, which is less nutrient efficient than respiration. There is hence a trade-off between endurance and power. Metabolic models have been developed to study how limited catalytic capacity of enzymes affects ATP synthesis. Here we integrate an enzyme-constrained metabolic model with proteomics data from muscle fibers. We find that ATP synthesis is constrained by several enzymes. A metabolic bypass of mitochondrial complex I is found to increase the ATP synthesis rate per gram of protein compared to full respiration. To test if this metabolic mode occurs in vivo, we conduct a high resolved incremental exercise tests for five subjects. Their gas exchange at different work rates is accurately reproduced by a whole-body metabolic model incorporating complex I bypass. The study therefore shows how proteome allocation influences metabolism during high intensity exercise.

Effects of altered pyruvate dehydrogenase activity on contracting skeletal muscle bioenergetics

During aerobic exercise (>65% of maximum oxygen consumption), the primary source of acetyl-CoA to fuel oxidative ATP synthesis in muscle is the pyruvate dehydrogenase (PDH) reaction. This study investigated how regulation of PDH activity affects muscle energetics by determining whether activation of PDH with dichloroacetate (DCA) alters the dynamics of the phosphate potential of rat gastrocnemius muscle during contraction. Twitch contractions were induced in vivo over a broad range of intensities to sample submaximal and maximal aerobic workloads. Muscle phosphorus metabolites were measured in vivo before and after DCA treatment by phosphorus nuclear magnetic resonance spectroscopy. At rest, DCA increased PDH activation compared with control (90 ± 12% vs. 23 ± 3%, P < 0.05), with parallel decreases in inorganic phosphate (P_i) of 17% (1.4 ± 0.2 vs. 1.7 ± 0.1 mM, P < 0.05) and an increase in the free energy of ATP hydrolysis (ΔG_ATP) (-66.2 ± 0.3 vs. -65.6 ± 0.2 kJ/mol, P < 0.05). During stimulation DCA increased steady-state phosphocreatine (PCr) and the magnitude of ΔG_ATP, with concomitant reduction in P_i and ADP concentrations. These effects were not due to kinetic alterations in PCr hydrolysis, resynthesis, or glycolytic ATP production and altered the flow-force relationship between mitochondrial ATP synthesis rate and ΔG_ATP. DCA had no significant effect at 1.0- to 2.0-Hz stimulation because physiological mechanisms at these high stimulation levels cause maximal activation of PDH. These data support a role of PDH activation in the regulation of the energetic steady state by altering the phosphate potential (ΔG_ATP) at rest and during contraction.

Open-Loop Control of Oxidative Phosphorylation in Skeletal and Cardiac Muscle Mitochondria by Ca(2.)

In cardiac muscle, mitochondrial ATP synthesis is driven by demand for ATP through feedback from the products of ATP hydrolysis. However, in skeletal muscle at higher workloads there is an apparent contribution of open-loop stimulation of ATP synthesis. Open-loop control is defined as modulation of flux through a biochemical pathway by a moiety, which is not a reactant or a product of the biochemical reactions in the pathway. The role of calcium, which is known to stimulate the activity of mitochondrial dehydrogenases, as an open-loop controller, was investigated in isolated cardiac and skeletal muscle mitochondria. The kinetics of NADH synthesis and respiration, feedback from ATP hydrolysis products, and stimulation by calcium were characterized in isolated mitochondria to test the hypothesis that calcium has a stimulatory role in skeletal muscle mitochondria not apparent in cardiac mitochondria. A range of respiratory states were obtained in cardiac and skeletal muscle mitochondria utilizing physiologically relevant concentrations of pyruvate and malate, and flux of respiration, NAD(P)H fluorescence, and rhodamine 123 fluorescence were measured over a range of extra mitochondrial calcium concentrations. We found that under these conditions calcium stimulates NADH synthesis in skeletal muscle mitochondria but not in cardiac mitochondria.

Living on the edge: substrate competition explains loss of robustness in mitochondrial fatty-acid oxidation disorders

Background

Defects in genes involved in mitochondrial fatty-acid oxidation (mFAO) reduce the ability of patients to cope with metabolic challenges. mFAO enzymes accept multiple substrates of different chain length, leading to molecular competition among the substrates. Here, we combined computational modeling with quantitative mouse and patient data to investigate whether substrate competition affects pathway robustness in mFAO disorders.

Results

First, we used comprehensive biochemical analyses of wild-type mice and mice deficient for medium-chain acyl-CoA dehydrogenase (MCAD) to parameterize a detailed computational model of mFAO. Model simulations predicted that MCAD deficiency would have no effect on the pathway flux at low concentrations of the mFAO substrate palmitoyl-CoA. However, high concentrations of palmitoyl-CoA would induce a decline in flux and an accumulation of intermediate metabolites. We proved computationally that the predicted overload behavior was due to substrate competition in the pathway. Second, to study the clinical relevance of this mechanism, we used patients’ metabolite profiles and generated a humanized version of the computational model. While molecular competition did not affect the plasma metabolite profiles during MCAD deficiency, it was a key factor in explaining the characteristic acylcarnitine profiles of multiple acyl-CoA dehydrogenase deficient patients. The patient-specific computational models allowed us to predict the severity of the disease phenotype, providing a proof of principle for the systems medicine approach.

Conclusion

We conclude that substrate competition is at the basis of the physiology seen in patients with mFAO disorders, a finding that may explain why these patients run a risk of a life-threatening metabolic catastrophe.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-016-0327-5) contains supplementary material, which is available to authorized users.

Synthesis, characterization, and in-vitro cytocompatibility of amorphous β-tri-calcium magnesium phosphate ceramics

Biphasic mixtures of crystalline β-tricalcium magnesium phosphate (β-TCMP) and an amorphous calcium magnesium phosphate have been synthesized and reported to support enhanced hMSC differentiation in comparison to β-tricalcium phosphate (β-TCP) due to the release of increased amounts of bioactive ions. In the current study, completely amorphous β-TCMP has been synthesized which is capable of releasing increased amounts of Mg(2+) and PO4(3-) ions, rather than a biphasic mixture as earlier reported. The amorphous phase formed was observed to crystallize between temperatures of 400-600°C. The scaffolds prepared with amorphous β-TCMP were capable of supporting enhanced hMSC proliferation and differentiation in comparison to commercially available β-TCP. However, a similar gene expression of mature osteoblast markers, OCN and COL-1, in comparison to biphasic β-TCMP was observed. To further study the role of Mg(2+) and PO4(3-) ions in regulating hMSC osteogenic differentiation, the capability of hMSCs to mineralize in growth media supplemented with Mg(2+) and PO4(3-) ions was studied. Interestingly, 5mM PO4(3-) supported mineralization while the addition of 5mM Mg(2+) to 5mM PO4(3-) inhibited mineralization. It was therefore concluded that the release of Ca(2+) ions from β-TCMP scaffolds also plays a role in regulating osteogenic differentiation on these scaffolds and it is noted that further work is required to more accurately determine the exact role of Mg(2+) in regulating hMSC osteogenic differentiation.

Improving the physiological realism of experimental models

The Virtual Physiological Human (VPH) project aims to develop integrative, explanatory and predictive computational models (C-Models) as numerical investigational tools to study disease, identify and design effective therapies and provide an in silico platform for drug screening. Ultimately, these models rely on the analysis and integration of experimental data. As such, the success of VPH depends on the availability of physiologically realistic experimental models (E-Models) of human organ function that can be parametrized to test the numerical models. Here, the current state of suitable E-models, ranging from in vitro non-human cell organelles to in vivo human organ systems, is discussed. Specifically, challenges and recent progress in improving the physiological realism of E-models that may benefit the VPH project are highlighted and discussed using examples from the field of research on cardiovascular disease, musculoskeletal disorders, diabetes and Parkinson's disease.

Exercise oxidative skeletal muscle metabolism in adolescents with cystic fibrosis

NEW FINDINGS

What is the central question of this study? Do intrinsic abnormalities in oxygenation and/or muscle oxidative metabolism contribute to exercise intolerance in adolescents with mild cystic fibrosis? What is the main finding and its importance? This study found no evidence that in adolescents with mild cystic fibrosis in a stable clinical state intrinsic abnormalities in skeletal muscle oxidative metabolism seem to play a clinical significant role. Based on these results, we concluded that there is no metabolic constraint to benefit from exercise training. Patients with cystic fibrosis (CF) are reported to have limited exercise capacity. There is no consensus about a possible abnormality in skeletal muscle oxidative metabolism in CF. Our aim was to test the hypothesis that abnormalities in oxygenation and/or muscle oxidative metabolism contribute to exercise intolerance in adolescents with mild CF. Ten adolescents with CF (12-18 years of age; forced expiratory volume in 1 s >80% of predicted; and resting oxygen saturation >94%) and 10 healthy age-matched control (HC) subjects were tested with supine cycle ergometry using near-infrared spectroscopy and (31)P magnetic resonance spectroscopy to study skeletal muscle oxygenation and oxidative metabolism during rest, exercise and recovery. No statistically significant (P > 0.1) differences in peak workload and peak oxygen uptake per kilogram lean body mass were found between CF and HC subjects. No differences were found between CF and HC subjects in bulk changes of quadriceps phosphocreatine (P = 0.550) and inorganic phosphate (P = 0.896) content and pH (P = 0.512) during symptom-limited exercise. Furthermore, we found statistically identical kinetics for phosphocreatine resynthesis during recovery for CF and HC subjects (P = 0.53). No statistically significant difference in peak exercise arbitrary units for total haemoglobin content was found between CF and HC subjects (P = 0.66). The results of this study provide evidence that in patients with mild CF and a stable clinical status (without signs of systemic inflammation and/or chronic Pseudomonas aeruginosa colonization), no intrinsic metabolic constraints and/or abnormalities in oxygenation and/or muscle oxidative metabolism contribute to exercise intolerance.

Altered Energetics of Exercise Explain Risk of Rhabdomyolysis in Very Long-Chain Acyl-CoA Dehydrogenase Deficiency

Rhabdomyolysis is common in very long-chain acyl-CoA dehydrogenase deficiency (VLCADD) and other metabolic myopathies, but its pathogenic basis is poorly understood. Here, we show that prolonged bicycling exercise against a standardized moderate workload in VLCADD patients is associated with threefold bigger changes in phosphocreatine (PCr) and inorganic phosphate (Pi) concentrations in quadriceps muscle and twofold lower changes in plasma acetyl-carnitine levels than in healthy subjects. This result is consistent with the hypothesis that muscle ATP homeostasis during exercise is compromised in VLCADD. However, the measured rates of PCr and Pi recovery post-exercise showed that the mitochondrial capacity for ATP synthesis in VLCADD muscle was normal. Mathematical modeling of oxidative ATP metabolism in muscle composed of three different fiber types indicated that the observed altered energy balance during submaximal exercise in VLCADD patients may be explained by a slow-to-fast shift in quadriceps fiber-type composition corresponding to 30% of the slow-twitch fiber-type pool in healthy quadriceps muscle. This study demonstrates for the first time that quadriceps energy balance during exercise in VLCADD patients is altered but not because of failing mitochondrial function. Our findings provide new clues to understanding the risk of rhabdomyolysis following exercise in human VLCADD.

Each-step activation of oxidative phosphorylation is necessary to explain muscle metabolic kinetic responses to exercise and recovery in humans

KEY POINTS

The basic control mechanisms of oxidative phosphorylation (OXPHOS) and glycolysis during work transitions in human skeletal muscle are still a matter of debate. We used simulations of skeletal muscle bioenergetics to identify key system features that contribute to this debate, by comparing kinetic model outputs with experimental human data, including phosphocreatine, pH, pulmonary oxygen uptake and fluxes of ATP production by OXPHOS (vOX), anaerobic glycolysis and creatine kinase in moderate and severe intensity exercise transitions. We found that each-step activation of particular OXPHOS complexes, NADH supply and glycolysis, and strong (third-order) glycolytic inhibition by protons was required to reproduce observed phosphocreatine, pH and vOX kinetics during exercise. A slow decay of each-step activation during recovery, which was slowed further following severe exercise, was necessary to reproduce the experimental findings. Well-tested computer models offer new insight in the control of the human skeletal muscle bioenergetic system during physical exercise.

ABSTRACT

To better understand muscle bioenergetic regulation, a previously-developed model of the skeletal muscle cell bioenergetic system was used to simulate the influence of: (1) each-step activation (ESA) of NADH supply (including glycolysis) and oxidative phosphorylation (OXPHOS) complexes and (2) glycolytic inhibition by protons on the kinetics of ATP synthesis from OXPHOS, anaerobic glycolysis and creatine kinase. Simulations were fitted to previously published experimental data of ATP production fluxes and metabolite concentrations during moderate and severe intensity exercise transitions in bilateral knee extension in humans. Overall, the computer simulations agreed well with experimental results. Specifically, a large (>5-fold) direct activation of all OXPHOS complexes was required to simulate measured phosphocreatine and OXPHOS responses to both moderate and severe intensity exercise. In addition, slow decay of ESA was required to fit phosphocreatine recovery kinetics, and the time constant of ESA decay was slower following severe (180 s) than moderate (90 s) exercise. Additionally, a strong inhibition of (anaerobic) glycolysis by protons (glycolytic rate inversely proportional to the cube of proton concentration) provided the best fit to the experimental pH kinetics, and may contribute to the progressive increase in oxidative ATP supply during acidifying contractions. During severe-intensity exercise, an 'additional' ATP usage (a 27% increase at 8 min, above the initial ATP supply) was necessary to explain the observed V̇O2 slow component. Thus, parallel activation of ATP usage and ATP supply (ESA), and a strong inhibition of ATP supply by anaerobic glycolysis, were necessary to simulate the kinetics of muscle bioenergetics observed in humans.

ATGL-mediated triglyceride turnover and the regulation of mitochondrial capacity in skeletal muscle

Emerging evidence indicates that skeletal muscle lipid droplets are an important control point for intracellular lipid homeostasis and that regulating fatty acid fluxes from lipid droplets might influence mitochondrial capacity. We used pharmacological blockers of the major triglyceride lipases, adipose triglyceride lipase (ATGL) and hormone-sensitive lipase, to show that a large proportion of the fatty acids that are transported into myotubes are trafficked through the intramyocellular triglyceride pool. We next tested whether increasing lipolysis from intramyocellular lipid droplets could activate transcriptional responses to enhance mitochondrial and fatty acid oxidative capacity. ATGL was overexpressed by adenoviral and adenoassociated viral infection in C2C12 myotubes and the tibialis anterior muscle of C57Bl/6 mice, respectively. ATGL overexpression in C2C12 myotubes increased lipolysis, which was associated with increased peroxisome proliferator-activated receptor (PPAR)-∂ activity, transcriptional upregulation of some PPAR∂ target genes, and enhanced mitochondrial capacity. The transcriptional responses were specific to ATGL actions and not a generalized increase in fatty acid flux in the myotubes. Marked ATGL overexpression (20-fold) induced modest molecular changes in the skeletal muscle of mice, but these effects were not sufficient to alter fatty acid oxidation. Together, these data demonstrate the importance of lipid droplets for myocellular fatty acid trafficking and the capacity to modulate mitochondrial capacity by enhancing lipid droplet lipolysis in vitro; however, this adaptive program is of minor importance when superimposing the normal metabolic stresses encountered in free-moving animals.

'Idealized' state 4 and state 3 in mitochondria vs. rest and work in skeletal muscle

A computer model of oxidative phosphorylation (OXPHOS) in skeletal muscle is used to compare state 3, intermediate state and state 4 in mitochondria with rest and work in skeletal muscle. 'Idealized' state 4 and 3 in relation to various 'experimental' states 4 and 3 are defined. Theoretical simulations show, in accordance with experimental data, that oxygen consumption (V'O2), ADP and Pi are higher, while ATP/ADP and Δp are lower in rest than in state 4, because of the presence of basal ATP consuming reactions in the former. It is postulated that moderate and intensive work in skeletal muscle is very different from state 3 in isolated mitochondria. V'O2, ATP/ADP, Δp and the control of ATP usage over V'O2 are much higher, while ADP and Pi are much lower in the former. The slope of the phenomenological V'O2-ADP relationship is much steeper during the rest-work transition than during the state 4-state 3 transition. The work state in intact muscle is much more similar to intermediate state than to state 3 in isolated mitochondria in terms of ADP, ATP/ADP, Δp and metabolic control pattern, but not in terms of V'O2. The huge differences between intact muscle and isolated mitochondria are proposed to be caused by the presence of the each-step activation (ESA) mechanism of the regulation of OXPHOS in intact skeletal muscle. Generally, the present study suggests that isolated mitochondria (at least in the absence of Ca2+) cannot serve as a good model of OXPHOS regulation in intact skeletal muscle.

Quantification of skeletal muscle mitochondrial function by 31P magnetic resonance spectroscopy techniques: a quantitative review

Magnetic resonance spectroscopy (MRS) can give information about cellular metabolism in vivo which is difficult to obtain in other ways. In skeletal muscle, non-invasive (31) P MRS measurements of the post-exercise recovery kinetics of pH, [PCr], [Pi] and [ADP] contain valuable information about muscle mitochondrial function and cellular pH homeostasis in vivo, but quantitative interpretation depends on understanding the underlying physiology. Here, by giving examples of the analysis of (31) P MRS recovery data, by some simple computational simulation, and by extensively comparing data from published studies using both (31) P MRS and invasive direct measurements of muscle O2 consumption in a common analytical framework, we consider what can be learnt quantitatively about mitochondrial metabolism in skeletal muscle using MRS-based methodology. We explore some technical and conceptual limitations of current methods, and point out some aspects of the physiology which are still incompletely understood.

MITOCHONDRIA: investigation of in vivo muscle mitochondrial function by 31P magnetic resonance spectroscopy

The most important function of mitochondria is the production of energy in the form of ATP. The socio-economic impact of human diseases that affect skeletal muscle mitochondrial function is growing, and improving their clinical management critically depends on the development of non-invasive assays to assess mitochondrial function and monitor the effects of interventions. 31P magnetic resonance spectroscopy provides two approaches that have been used to assess in vivo ATP synthesis in skeletal muscle: measuring Pi→ATP exchange flux using saturation transfer in resting muscle, and measuring phosphocreatine recovery kinetics after exercise. However, Pi→ATP exchange does not represent net mitochondrial ATP synthesis flux and has no simple relationship with mitochondrial function. Post-exercise phosphocreatine recovery kinetics, on the other hand, yield reliable measures of muscle mitochondrial capacity in vivo, whose ability to define the site of functional defects is enhanced by combination with other non-invasive techniques.

Optimal experiment design for model selection in biochemical networks

Background

Mathematical modeling is often used to formalize hypotheses on how a biochemical network operates by discriminating between competing models. Bayesian model selection offers a way to determine the amount of evidence that data provides to support one model over the other while favoring simple models. In practice, the amount of experimental data is often insufficient to make a clear distinction between competing models. Often one would like to perform a new experiment which would discriminate between competing hypotheses.

Results

We developed a novel method to perform Optimal Experiment Design to predict which experiments would most effectively allow model selection. A Bayesian approach is applied to infer model parameter distributions. These distributions are sampled and used to simulate from multivariate predictive densities. The method is based on a k-Nearest Neighbor estimate of the Jensen Shannon divergence between the multivariate predictive densities of competing models.

Conclusions

We show that the method successfully uses predictive differences to enable model selection by applying it to several test cases. Because the design criterion is based on predictive distributions, which can be computed for a wide range of model quantities, the approach is very flexible. The method reveals specific combinations of experiments which improve discriminability even in cases where data is scarce. The proposed approach can be used in conjunction with existing Bayesian methodologies where (approximate) posteriors have been determined, making use of relations that exist within the inferred posteriors.

Exercise: Kinetic considerations for gas exchange

The activities of daily living typically occur at metabolic rates below the maximum rate of aerobic energy production. Such activity is characteristic of the nonsteady state, where energy demands, and consequential physiological responses, are in constant flux. The dynamics of the integrated physiological processes during these activities determine the degree to which exercise can be supported through rates of O₂ utilization and CO₂ clearance appropriate for their demands and, as such, provide a physiological framework for the notion of exercise intensity. The rate at which O₂ exchange responds to meet the changing energy demands of exercise--its kinetics--is dependent on the ability of the pulmonary, circulatory, and muscle bioenergetic systems to respond appropriately. Slow response kinetics in pulmonary O₂ uptake predispose toward a greater necessity for substrate-level energy supply, processes that are limited in their capacity, challenge system homeostasis and hence contribute to exercise intolerance. This review provides a physiological systems perspective of pulmonary gas exchange kinetics: from an integrative view on the control of muscle oxygen consumption kinetics to the dissociation of cellular respiration from its pulmonary expression by the circulatory dynamics and the gas capacitance of the lungs, blood, and tissues. The intensity dependence of gas exchange kinetics is discussed in relation to constant, intermittent, and ramped work rate changes. The influence of heterogeneity in the kinetic matching of O₂ delivery to utilization is presented in reference to exercise tolerance in endurance-trained athletes, the elderly, and patients with chronic heart or lung disease.

Regulation of oxidative phosphorylation during work transitions results from its kinetic properties

The regulation of oxidative phosphorylation (OXPHOS) during work transitions in skeletal muscle and heart is still not well understood. Different computer models of this process have been developed that are characterized by various kinetic properties. In the present research-polemic theoretical study it is argued that models belonging to one group (Model A), which predict that among OXPHOS complexes complex III keeps almost all of the metabolic control over oxygen consumption (Vo2) and involve a strong complex III activation by inorganic phosphate (Pi), lead to the conclusion that an increase in Pi is the main mechanism responsible for OXPHOS activation (feedback-activation mechanism). Models belonging to another group (Model B), which were developed to take into account an approximately uniform distribution of metabolic control over Vo2 among particular OXPHOS complexes (complex I, complex III, complex IV, ATP synthase, ATP/ADP carrier, phosphate carrier) encountered in experimental studies in isolated mitochondria, predict that all OXPHOS complexes are directly activated in parallel with ATP usage and NADH supply by some external cytosolic factor/mechanism during rest-to-work or low-to-high work transitions in skeletal muscle and heart ("each-step-activation" mechanism). Model B demonstrates that different intensities of each-step activation can account for the very different (slopes of) phenomenological Vo2-ADP relationships observed in various skeletal muscles and heart. Thus they are able to explain the differences in the regulation of OXPHOS during work transitions between skeletal muscle (where moderate changes in ADP take place) and intact heart in vivo (where ADP is essentially constant).

31P MR spectroscopy and computational modeling identify a direct relation between Pi content of an alkaline compartment in resting muscle and phosphocreatine resynthesis kinetics in active muscle in humans

The assessment of mitochondrial properties in skeletal muscle is important in clinical research, for instance in the study of diabetes. The gold standard to measure mitochondrial capacity non-invasively is the phosphocreatine (PCr) recovery rate after exercise, measured by (31)P Magnetic Resonance spectroscopy ((31)P MRS). Here, we sought to expand the evidence base for an alternative method to assess mitochondrial properties which uses (31)P MRS measurement of the Pi content of an alkaline compartment attributed to mitochondria (Pi2; as opposed to cytosolic Pi (Pi1)) in resting muscle at high magnetic field. Specifically, the PCr recovery rate in human quadriceps muscle was compared with the signal intensity of the Pi2 peak in subjects with varying mitochondrial content of the quadriceps muscle as a result of athletic training, and the results were entered into a mechanistic computational model of mitochondrial metabolism in muscle to test if the empirical relation between Pi2/Pi1 ratio and the PCr recovery was consistent with theory. Localized (31)P spectra were obtained at 7T from resting vastus lateralis muscle to measure the intensity of the Pi2 peak. In the endurance trained athletes a Pi2/Pi1 ratio of 0.07 ± 0.01 was found, compared to a significantly lower (p<0.05) Pi2/Pi1 ratio of 0.03 ± 0.01 in the normally active group. Next, PCr recovery kinetics after in magnet bicycle exercise were measured at 1.5T. For the endurance trained athletes, a time constant τPCr 12 ± 3 s was found, compared to 24 ± 5s in normally active subjects. Without any parameter optimization the computational model prediction matched the experimental data well (r(2) of 0.75). Taken together, these results suggest that the Pi2 resonance in resting human skeletal muscle observed at 7T provides a quantitative MR-based functional measure of mitochondrial density.

Slow V̇o2 off-kinetics in skeletal muscle is associated with fast PCr off-kinetics—and inversely

The computer model of the bioenergetic system in skeletal muscle, developed previously, was used to study the effect of the characteristic decay time of the parallel activation of oxidative phosphorylation [τ(OFF)] during muscle recovery on the muscle oxygen consumption rate (V̇o2) and phosphocreatine (PCr) work-to-rest transition (off)-kinetics and on the relationship between the V̇o2 and PCr rest-to-work transition (on)- and off-kinetics in moderate and heavy exercise. An increase in τ(OFF) slows down the initial phase of the muscle V̇o2 off-kinetics and accelerates the PCr off-kinetics. As a result, the relationship between the initial phase of the V̇o2 off-kinetics (lasting approximately 3–60 s in computer simulations) and the PCr off-kinetics is inverse: the slower the former, the faster the latter. A faster initial phase of the V̇o2 off-kinetics is associated with a slower late phase of the V̇o2 off-kinetics, and as a result, the integral of V̇o2 above baseline during recovery, representing the oxygen debt, is identical in all cases [values of τ(OFF)] for a given PCr decrease. Depending on τ(OFF), the muscle V̇o2 on-kinetics was either equally fast or slower than the V̇o2 off-kinetics in moderate exercise and always slower in heavy exercise. PCr on-kinetics was always faster than PCr off-kinetics. This study clearly demonstrates that τ(OFF) has a pronounced impact on the mutual relations between the muscle V̇o2 and PCr on- and off-kinetics.

Parameter uncertainty in biochemical models described by ordinary differential equations

Improved mechanistic understanding of biochemical networks is one of the driving ambitions of Systems Biology. Computational modeling allows the integration of various sources of experimental data in order to put this conceptual understanding to the test in a quantitative manner. The aim of computational modeling is to obtain both predictive as well as explanatory models for complex phenomena, hereby providing useful approximations of reality with varying levels of detail. As the complexity required to describe different system increases, so does the need for determining how well such predictions can be made. Despite efforts to make tools for uncertainty analysis available to the field, these methods have not yet found widespread use in the field of Systems Biology. Additionally, the suitability of the different methods strongly depends on the problem and system under investigation. This review provides an introduction to some of the techniques available as well as gives an overview of the state-of-the-art methods for parameter uncertainty analysis.

Modeling to link regional myocardial work, metabolism and blood flows

Given the mono-functional, highly coordinated processes of cardiac excitation and contraction, the observations that regional myocardial blood flows, rMBF, are broadly heterogeneous has provoked much attention, but a clear explanation has not emerged. In isolated and in vivo heart studies the total coronary flow is found to be proportional to the rate-pressure product (systolic mean blood pressure times heart rate), a measure of external cardiac work. The same relationship might be expected on a local basis: more work requires more flow. The validity of this expectation has never been demonstrated experimentally. In this article we review the concepts linking cellular excitation and contractile work to cellular energetics and ATP demand, substrate utilization, oxygen demand, vasoregulation, and local blood flow. Mathematical models of these processes are now rather well developed. We propose that the construction of an integrated model encompassing the biophysics, biochemistry and physiology of cardiomyocyte contraction, then combined with a detailed three-dimensional structuring of the fiber bundle and sheet arrangements of the heart as a whole will frame an hypothesis that can be quantitatively evaluated to settle the prime issue: Does local work drive local flow in a predictable fashion that explains the heterogeneity? While in one sense one can feel content that work drives flow is irrefutable, the are no cardiac contractile models that demonstrate the required heterogeneity in local strain-stress-work; quite the contrary, cardiac contraction models have tended toward trying to show that work should be uniform. The object of this review is to argue that uniformity of work does not occur, and is impossible in any case, and that further experimentation and analysis are necessary to test the hypothesis.

Computational modeling of mitochondrial energy transduction

Mitochondria are the power plant of the heart, burning fat and sugars to supply the muscle with the adenosine triphosphate (ATP) free energy that drives contraction and relaxation during each heart beat. This function was first captured in a mathematical model in 1967. Today, interest in such a model has been rekindled by ongoing in silico integrative physiology efforts such as the Cardiac Physiome project. Here, the status of the field of computational modeling of mitochondrial ATP synthetic function is reviewed.

Regulation of mitochondrial biogenesis by lipoprotein lipase in muscle of insulin-resistant offspring of parents with type 2 diabetes

Recent studies reveal a strong relationship between reduced mitochondrial content and insulin resistance in human skeletal muscle, although the underlying factors responsible for this association remain unknown. To address this question, we analyzed muscle biopsy samples from young, lean, insulin resistant (IR) offspring of parents with type 2 diabetes and control subjects by microarray analyses and found significant differences in expression of ~512 probe pairs. We then screened these genes for their potential involvement in the regulation of mitochondrial biogenesis using RNA interference and found that mRNA and protein expression of lipoprotein lipase (LPL) in skeletal muscle was significantly decreased in the IR offspring and was associated with decreased mitochondrial density. Furthermore, we show that LPL knockdown in muscle cells decreased mitochondrial content by effectively decreasing fatty acid delivery and subsequent activation of peroxisome proliferator-activated receptor (PPAR)-δ. Taken together, these data suggest that decreased mitochondrial content in muscle of IR offspring may be due in part to reductions in LPL expression in skeletal muscle resulting in decreased PPAR-δ activation.

Prediction of muscle energy states at low metabolic rates requires feedback control of mitochondrial respiratory chain activity by inorganic phosphate

The regulation of the 100-fold dynamic range of mitochondrial ATP synthesis flux in skeletal muscle was investigated. Hypotheses of key control mechanisms were included in a biophysical model of oxidative phosphorylation and tested against metabolite dynamics recorded by (31)P nuclear magnetic resonance spectroscopy ((31)P MRS). Simulations of the initial model featuring only ADP and Pi feedback control of flux failed in reproducing the experimentally sampled relation between myoplasmic free energy of ATP hydrolysis (ΔG(p) = ΔG(p)(o')+RT ln ([ADP][Pi]/[ATP]) and the rate of mitochondrial ATP synthesis at low fluxes (<0.2 mM/s). Model analyses including Monte Carlo simulation approaches and metabolic control analysis (MCA) showed that this problem could not be amended by model re-parameterization, but instead required reformulation of ADP and Pi feedback control or introduction of additional control mechanisms (feed forward activation), specifically at respiratory Complex III. Both hypotheses were implemented and tested against time course data of phosphocreatine (PCr), Pi and ATP dynamics during post-exercise recovery and validation data obtained by (31)P MRS of sedentary subjects and track athletes. The results rejected the hypothesis of regulation by feed forward activation. Instead, it was concluded that feedback control of respiratory chain complexes by inorganic phosphate is essential to explain the regulation of mitochondrial ATP synthesis flux in skeletal muscle throughout its full dynamic range.

Molecular dynamics simulations of creatine kinase and adenine nucleotide translocase in mitochondrial membrane patch

Interaction between mitochondrial creatine kinase (MtCK) and adenine nucleotide translocase (ANT) can play an important role in determining energy transfer pathways in the cell. Although the functional coupling between MtCK and ANT has been demonstrated, the precise mechanism of the coupling is not clear. To study the details of the coupling, we turned to molecular dynamics simulations. We introduce a new coarse-grained molecular dynamics model of a patch of the mitochondrial inner membrane containing a transmembrane ANT and an MtCK above the membrane. The membrane model consists of three major types of lipids (phosphatidylcholine, phosphatidylethanolamine, and cardiolipin) in a roughly 2:1:1 molar ratio. A thermodynamics-based coarse-grained force field, termed MARTINI, has been used together with the GROMACS molecular dynamics package for all simulated systems in this work. Several physical properties of the system are reproduced by the model and are in agreement with known data. This includes membrane thickness, dimension of the proteins, and diffusion constants. We have studied the binding of MtCK to the membrane and demonstrated the effect of cardiolipin on the stabilization of the binding. In addition, our simulations predict which part of the MtCK protein sequence interacts with the membrane. Taken together, the model has been verified by dynamical and structural data and can be used as the basis for further studies.

Modeling of mitochondria bioenergetics using a composable chemiosmotic energy transduction rate law: theory and experimental validation

Mitochondrial bioenergetic processes are central to the production of cellular energy, and a decrease in the expression or activity of enzyme complexes responsible for these processes can result in energetic deficit that correlates with many metabolic diseases and aging. Unfortunately, existing computational models of mitochondrial bioenergetics either lack relevant kinetic descriptions of the enzyme complexes, or incorporate mechanisms too specific to a particular mitochondrial system and are thus incapable of capturing the heterogeneity associated with these complexes across different systems and system states. Here we introduce a new composable rate equation, the chemiosmotic rate law, that expresses the flux of a prototypical energy transduction complex as a function of: the saturation kinetics of the electron donor and acceptor substrates; the redox transfer potential between the complex and the substrates; and the steady-state thermodynamic force-to-flux relationship of the overall electro-chemical reaction. Modeling of bioenergetics with this rate law has several advantages: (1) it minimizes the use of arbitrary free parameters while featuring biochemically relevant parameters that can be obtained through progress curves of common enzyme kinetics protocols; (2) it is modular and can adapt to various enzyme complex arrangements for both in vivo and in vitro systems via transformation of its rate and equilibrium constants; (3) it provides a clear association between the sensitivity of the parameters of the individual complexes and the sensitivity of the system's steady-state. To validate our approach, we conduct in vitro measurements of ETC complex I, III, and IV activities using rat heart homogenates, and construct an estimation procedure for the parameter values directly from these measurements. In addition, we show the theoretical connections of our approach to the existing models, and compare the predictive accuracy of the rate law with our experimentally fitted parameters to those of existing models. Finally, we present a complete perturbation study of these parameters to reveal how they can significantly and differentially influence global flux and operational thresholds, suggesting that this modeling approach could help enable the comparative analysis of mitochondria from different systems and pathological states. The procedures and results are available in Mathematica notebooks at http://www.igb.uci.edu/tools/sb/mitochondria-modeling.html.

Fuel oxidation at the walk-to-run-transition in humans

Multiple factors (including anthropometric, kinetic, mechanical, kinematic, perceptual, and energetic factors) are likely to play a role in the walk-to-run transition in humans. The primary purpose of the present study was to consider an additional factor, the metabolic fuel source. Indirect calorimetry was used to measure fuel oxidation, and perception of effort was recorded as 10 overnight-fasted adults locomoted on a level treadmill at speeds progressing from 1.56 to 2.46 m s(-1) in increments of 0.11 m s(-1) and 10.0 minutes under 3 conditions: (1) unconstrained choice of gait, (2) walking at all speeds, and (3) running at all speeds. The preferred transition speed was 2.08 ± 0.03 m s(-1). Gait transition from walking to running increased oxygen consumption rate, decreased the perception of effort, and decreased the rate of carbohydrate oxidation. We propose that, in an evolutionary context, gait transition, guided by the perception of effort, can be viewed as a carbohydrate-sparing strategy.

Similar mitochondrial activation kinetics in wild-type and creatine kinase-deficient fast-twitch muscle indicate significant Pi control of respiration

Past simulations of oxidative ATP metabolism in skeletal muscle have predicted that elimination of the creatine kinase (CK) reaction should result in dramatically faster oxygen consumption dynamics during transitions in ATP turnover rate. This hypothesis was investigated. Oxygen consumption of fast-twitch (FT) muscle isolated from wild-type (WT) and transgenic mice deficient in the myoplasmic (M) and mitochondrial (Mi) CK isoforms (MiM CK(-/-)) were measured at 20°C at rest and during electrical stimulation. MiM CK(-/-) muscle oxygen consumption activation kinetics during a step change in contraction rate were 30% faster than WT (time constant 53 ± 3 vs. 69 ± 4 s, respectively; mean ± SE, n = 8 and 6, respectively). MiM CK(-/-) muscle oxygen consumption deactivation kinetics were 380% faster than WT (time constant 74 ± 4 s vs. 264 ± 4 s, respectively). Next, the experiments were simulated using a computational model of the oxidative ATP metabolic network in FT muscle featuring ADP and Pi feedback control of mitochondrial respiration (J. A. L. Jeneson, J. P. Schmitz, N. A. van den Broek, N. A. van Riel, P. A. Hilbers, K. Nicolay, J. J. Prompers. Am J Physiol Endocrinol Metab 297: E774-E784, 2009) that was reparameterized for 20°C. Elimination of Pi control via clamping of the mitochondrial Pi concentration at 10 mM reproduced past simulation results of dramatically faster kinetics in CK(-/-) muscle, while inclusion of Pi control qualitatively explained the experimental observations. On this basis, it was concluded that previous studies of the CK-deficient FT muscle phenotype underestimated the contribution of Pi to mitochondrial respiratory control.

Thermodynamic calculations for biochemical transport and reaction processes in metabolic networks

Thermodynamic analysis of metabolic networks has recently generated increasing interest for its ability to add constraints on metabolic network operation, and to combine metabolic fluxes and metabolite measurements in a mechanistic manner. Concepts for the calculation of the change in Gibbs energy of biochemical reactions have long been established. However, a concept for incorporation of cross-membrane transport in these calculations is still missing, although the theory for calculating thermodynamic properties of transport processes is long known. Here, we have developed two equivalent equations to calculate the change in Gibbs energy of combined transport and reaction processes based on two different ways of treating biochemical thermodynamics. We illustrate the need for these equations by showing that in some cases there is a significant difference between the proposed correct calculation and using an approximative method. With the developed equations, thermodynamic analysis of metabolic networks spanning over multiple physical compartments can now be correctly described.

Elevated baseline VO2 per se does not slow O2 uptake kinetics during work-to-work exercise transitions

We investigated whether the characteristic slowing of pulmonary oxygen uptake (VO2) kinetics during "work-to-work" exercise is attributable to elevations in baseline metabolic rate (VO2) as opposed to the elevated baseline work rate, per se. We hypothesized that a step transition to a higher work rate from "unloaded" cycling, but with elevations in VO2 [and heart rate (HR)] reflective of a work-to-work transition, would result in a lengthened phase II time constant (τ(p)). Seven male subjects (mean ± SD age 27 ± 10 yr) completed 1) transitions to a high-intensity work rate from a moderate-intensity work rate (M→H) and 2) two consecutive bouts of high-intensity exercise (U→H and E→H, respectively) initiated from unloaded cycling, with the time separating the exercise bouts chosen such that the baseline VO2 for the second transition was similar to the baseline VO2 for the M→H transition. The τ(p) for M→H (48 ± 16 s) was significantly greater (P < 0.05) than the τ(p) for U→H (28 ± 8 s) and E→H (27 ± 6 s), which did not differ significantly. These findings suggest that the altered VO2 dynamics that are observed during work-to-work exercise are not related to the elevated baseline VO2 (or HR) per se; rather, these effects appear to be linked to the elevated baseline work rate, which would be expected to dictate the subsequent muscle fiber recruitment profile.

Comparison of in vivo postexercise phosphocreatine recovery and resting ATP synthesis flux for the assessment of skeletal muscle mitochondrial function

(31)P magnetic resonance spectroscopy (MRS) has been used to assess skeletal muscle mitochondrial function in vivo by measuring 1) phosphocreatine (PCr) recovery after exercise or 2) resting ATP synthesis flux with saturation transfer (ST). In this study, we compared both parameters in a rat model of mitochondrial dysfunction with the aim of establishing the most appropriate method for the assessment of in vivo muscle mitochondrial function. Mitochondrial dysfunction was induced in adult Wistar rats by daily subcutaneous injections with the complex I inhibitor diphenyleneiodonium (DPI) for 2 wk. In vivo (31)P MRS measurements were supplemented by in vitro measurements of oxygen consumption in isolated mitochondria. Two weeks of DPI treatment induced mitochondrial dysfunction, as evidenced by a 20% lower maximal ADP-stimulated oxygen consumption rate in isolated mitochondria from DPI-treated rats oxidizing pyruvate plus malate. This was paralleled by a 46% decrease in in vivo oxidative capacity, determined from postexercise PCr recovery. Interestingly, no significant difference in resting, ST-based ATP synthesis flux was observed between DPI-treated rats and controls. These results show that PCr recovery after exercise has a more direct relationship with skeletal muscle mitochondrial function than the ATP synthesis flux measured with (31)P ST MRS in the resting state.

Magnitude and control of mitochondrial sensitivity to ADP

The transduction function for ADP stimulation of mitochondrial ATP synthesis in skeletal muscle was reconstructed in vivo and in silico to investigate the magnitude and origin of mitochondrial sensitivity to cytoplasmic ADP concentration changes. Dynamic in vivo measurements of human leg muscle phosphocreatine (PCr) content during metabolic recovery from contractions were performed by (31)P-NMR spectroscopy. The cytoplasmic ADP concentration ([ADP]) and rate of oxidative ATP synthesis (Jp) at each time point were calculated from creatine kinase equilibrium and the derivative of a monoexponential fit to the PCr recovery data, respectively. Reconstructed [ADP]-Jp relations for individual muscles containing more than 100 data points were kinetically characterized by nonlinear curve fitting yielding an apparent kinetic order and ADP affinity of 1.9 +/- 0.2 and 0.022 +/- 0.003 mM, respectively (means +/- SD; n = 6). Next, in silico [ADP]-Jp relations for skeletal muscle were generated using a computational model of muscle oxidative ATP metabolism whereby model parameters corresponding to mitochondrial enzymes were randomly changed by 50-150% to determine control of mitochondrial ADP sensitivity. The multiparametric sensitivity analysis showed that mitochondrial ADP ultrasensitivity is an emergent property of the integrated mitochondrial enzyme network controlled primarily by kinetic properties of the adenine nucleotide translocator.

Mitochondrial reactive oxygen species production in excitable cells: modulators of mitochondrial and cell function

The mitochondrion is a major source of reactive oxygen species (ROS). Superoxide (O(2)(*-)) is generated under specific bioenergetic conditions at several sites within the electron-transport system; most is converted to H(2)O(2) inside and outside the mitochondrial matrix by superoxide dismutases. H(2)O(2) is a major chemical messenger that, in low amounts and with its products, physiologically modulates cell function. The redox state and ROS scavengers largely control the emission (generation scavenging) of O(2)(*-). Cell ischemia, hypoxia, or toxins can result in excess O(2)(*-) production when the redox state is altered and the ROS scavenger systems are overwhelmed. Too much H(2)O(2) can combine with Fe(2+) complexes to form reactive ferryl species (e.g., Fe(IV) = O(*)). In the presence of nitric oxide (NO(*)), O(2)(*-) forms the reactant peroxynitrite (ONOO(-)), and ONOOH-induced nitrosylation of proteins, DNA, and lipids can modify their structure and function. An initial increase in ROS can cause an even greater increase in ROS and allow excess mitochondrial Ca(2+) entry, both of which are factors that induce cell apoptosis and necrosis. Approaches to reduce excess O(2)(*-) emission include selectively boosting the antioxidant capacity, uncoupling of oxidative phosphorylation to reduce generation of O(2)(*-) by inducing proton leak, and reversibly inhibiting electron transport. Mitochondrial cation channels and exchangers function to maintain matrix homeostasis and likely play a role in modulating mitochondrial function, in part by regulating O(2)(*-) generation. Cell-signaling pathways induced physiologically by ROS include effects on thiol groups and disulfide linkages to modify posttranslationally protein structure to activate/inactivate specific kinase/phosphatase pathways. Hypoxia-inducible factors that stimulate a cascade of gene transcription may be mediated physiologically by ROS. Our knowledge of the role played by ROS and their scavenging systems in modulation of cell function and cell death has grown exponentially over the past few years, but we are still limited in how to apply this knowledge to develop its full therapeutic potential.

Modular regulation analysis of integrative effects of hypoxia on the energetics of contracting skeletal muscle in vivo

In the exercising muscle, acute reduction in ambient oxygen impairs muscle contraction because of the effects of hypoxia on mitochondrial ATP supply. The less marked impairment reported after long-term exposure to hypoxia points to changes in the regulation of the energetic system of contraction in HC (hypoxic conditioned) animals. This energetic system is conceptually defined here as two modules: the ATP/PCr (phosphocreatine)-producer and the ATP/PCr-consumer connected by energetic intermediates. Modular control analysis that combines top-down control analysis with non-invasive 31P-NMR spectroscopy was used to describe the effects of hypoxia on each module and their adaptation. Modulations of steady levels of ATP turnover (indirectly assessed as force output) and muscle PCr were obtained in HC rats (6 weeks at 10.5% O2) compared with N (normoxic) rats. Modular control and regulation analyses quantified the elasticity to PCr of each module in N and HC rats as well as the direct effect of acute hypoxia on the ATP/PCr-producer module. Similar elasticities in N and HC rats indicate the absence of response to long-term hypoxia in internal regulations of the ATP supply and demand pathways. The less marked impairment of contraction by acute hypoxia in HC rats (-9+/-6% versus -17+/-14% in N rats, P<0.05) was therefore fully explained by a lower direct effect (HC -31+/-13% versus N -44+/-23%, P<0.05) of acute hypoxia on mitochondrial ATP supply. This points to a positive adaptation to chronic hypoxia. Modular control analysis in vivo may provide powerful tools to find out improved function (alternatively dysfunction) at the system level in conditioned animals.

Roles of the creatine kinase system and myoglobin in maintaining energetic state in the working heart

BACKGROUND

The heart is capable of maintaining contractile function despite a transient decrease in blood flow and increase in cardiac ATP demand during systole. This study analyzes a previously developed model of cardiac energetics and oxygen transport to understand the roles of the creatine kinase system and myoglobin in maintaining the ATP hydrolysis potential during beat-to-beat transient changes in blood flow and ATP hydrolysis rate.

RESULTS

The theoretical investigation demonstrates that elimination of myoglobin only slightly increases the predicted range of oscillation of cardiac oxygenation level during beat-to-beat transients in blood flow and ATP utilization. In silico elimination of myoglobin has almost no impact on the cytoplasmic ATP hydrolysis potential (DeltaGATPase). In contrast, disabling the creatine kinase system results in considerable oscillations of cytoplasmic ADP and ATP levels and seriously deteriorates the stability of DeltaGATPase in the beating heart.

CONCLUSION

The CK system stabilizes DeltaGATPase by both buffering ATP and ADP concentrations and enhancing the feedback signal of inorganic phosphate in regulating mitochondrial oxidative phosphorylation.

Mitochondrial function in physically active elders with sarcopenia

Physical activity is reported to protect against sarcopenia and preserve mitochondrial function. Healthy normal lean (NL: n=15) and sarcopenic (SS: n=9) participants were recruited based on body composition (DXA, Lunar DPX), age, and physical activity. Gastrocnemius mitochondrial function was assessed by (31)P MRS using steady-state exercise in a 4T Bruker Biospin. Total work (429.3+/-160.2J vs. 851.0+/-211.7J, p<0.001) and muscle volume (p=0.006) were lower in SS, although these variables were not correlated (NL r=-0.31, p=0.33, SS r=(0.03, p=0.93). In the SS resting ATP/ADP was lower (p=0.03) and ATP hydrolysis higher (p=0.02) at rest. Free energy ATP hydrolysis was greater at the end of exercise (p=0.02) and [ADP] relative to total work output was higher in SS (ANCOVA, p=0.005). [PCr] recovery kinetics were not different between the groups. Adjusting these parameters for differences in total work output and muscle volume did not explain these findings. These data suggest that aerobic metabolism in physically active older adults with sarcopenia is mildly impaired at rest and during modest levels of exercise where acidosis was avoided. Muscle energetics is coordinated at multiple cellular levels and further studies are needed to determine the loci/locus of energy instability in sarcopenia.