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 modelling identifies the impact of subtle anatomical variations between amphibian and mammalian skeletal muscle on spatiotemporal calcium dynamics

The physical sites of calcium entry and exit in the skeletal muscle cell are distinct and highly organised in space. It was investigated whether the highly structured spatial organisation of sites of Ca(2+) release, uptake and action in skeletal muscle cells substantially impacts the dynamics of cytosolic Ca(2+) handling and thereby the physiology of the cell. Hereto, the spatiotemporal dynamics of the free calcium distribution in a fast-twitch (FT) muscle sarcomere was studied using a reaction-diffusion computational model for two genotypes with known anatomical differences. A computational model of a murine FT muscle sarcomere is developed, de novo including a closed calcium mass balance to simulate spatiotemporal high stimulation frequency calcium dynamics at 35 degrees C. Literature data on high-frequency calcium dye measurements were used as a first step towards model validation. The murine and amphibian sarcomere models were phenotypically distinct to capture known differences in positions of troponin C, actin-myosin overlap and calcium release within the sarcomere between frog and mouse. The models predicted large calcium gradients throughout the myoplasm as well as differences in calcium concentrations near the mitochondria of frog and mouse. Furthermore, the predicted Ca(2+) concentration was high at positions where Ca(2+) has a regulatory function, close to the mitochondria and troponin C.

Modeling glucose and water dynamics in human skin

BACKGROUND

Glucose is heterogeneously distributed in the different physiological compartments in the human skin. Therefore, for the development of a noninvasive measurement method, both a good quantification of the different compartments of human skin and an understanding of glucose transport processes are important.

METHODS

The composition of human skin was quantified by histology research. Based on this information a mathematical model was developed to simulate glucose dynamics in human skin.

RESULTS

The model predicts dynamically glucose concentrations in the different layers of the skin as a result of changes in blood glucose concentration. The model was validated with published time course data of blood and interstitial fluid glucose during a clamp study with three different set points for blood glucose, and model outcomes were compared to measurements for the lag time and gradient. According to the model, glucose in the interstitial fluid of the dermis best matches the amplitude and dynamics of blood glucose.

CONCLUSIONS

The new data obtained from quantitative histology appeared crucial for the model. The proposed model was successfully validated. This result was obtained without tuning or fitting of any parameter. It was shown how the model can be used to set standards for measurements and to define the best measurement depth for noninvasive glucose monitoring.