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.

Combined in vivo and in silico investigations of activation of glycolysis in contracting skeletal muscle.

The hypothesis was tested that the variation of in vivo glycolytic flux with contraction frequency in skeletal muscle can be qualitatively and quantitatively explained by calcium-calmodulin activation of phosphofructokinase (PFK-1). Ischemic rat tibialis anterior muscle was electrically stimulated at frequencies between 0 and 80 Hz to covary the ATP turnover rate and calcium concentration in the tissue. Estimates of in vivo glycolytic rates and cellular free energetic states were derived from dynamic changes in intramuscular pH and phosphocreatine content, respectively, determined by phosphorus magnetic resonance spectroscopy ((31)P-MRS). Computational modeling was applied to relate these empirical observations to understanding of the biochemistry of muscle glycolysis. Hereto, the kinetic model of PFK activity in a previously reported mathematical model of the glycolytic pathway (Vinnakota KC, Rusk J, Palmer L, Shankland E, Kushmerick MJ. J Physiol 588: 1961-1983, 2010) was adapted to contain a calcium-calmodulin binding sensitivity. The two main results were introduction of regulation of PFK-1 activity by binding of a calcium-calmodulin complex in combination with activation by increased concentrations of AMP and ADP was essential to qualitatively and quantitatively explain the experimental observations. Secondly, the model predicted that shutdown of glycolytic ATP production flux in muscle postexercise may lag behind deactivation of PFK-1 (timescales: 5-10 s vs. 100-200 ms, respectively) as a result of accumulation of glycolytic intermediates downstream of PFK during contractions.

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.