(Matlab) Computational Model of Cellular Metabolic Dynamics in Skeletal Muscle Fibers during Moderate Intensity Exercise.

Description

(Abstract) Human skeletal muscles have different fiber types with distinct metabolic functions and physiological properties. The quantitative metabolic responses of muscle fibers to exercise provide essential information for understanding and modifying the regulatory mechanisms of skeletal muscle. Since in vivo data from skeletal muscle during exercise is limited, a computational, physiologically based model has been developed to quantify the dynamic metabolic responses of many key chemical species. This model distinguishes type I and II muscle fibers, which share the same blood supply. An underlying hypothesis is that the recruitment and metabolic activation of the two main types of muscle fibers differ depending on the pre-exercise state and exercise protocols. Here, activation measured by metabolic response (or enzymatic activation) in single fibers is considered linked but distinct from fiber recruitment characterized by the number (or mass) of each fiber type involved during a specific exercise. The model incorporates species transport processes between blood and muscle fibers and most of the important reactions/pathways in cytosol and mitochondria within each fiber type. Model simulations describe the dynamics of intracellular species concentrations and fluxes in muscle fibers during moderate intensity exercise according to various experimental protocols and conditions. This model is validated by comparing model simulations with experimental data in single muscle fibers and in whole muscle. Model simulations demonstrate that muscle-fiber recruitment and metabolic activation patterns in response to exercise produce significantly distinctive effects depending on the exercise conditions.

Figure: (from fig 3, Li et al. 2012) Simulated dynamics and experimental data for intracellular species concentrations in skeletal muscle in response to an increased work rate equivalent to cycle ergometer exercise at 50% VO2max for subjects corresponding to Experiment 1: (left to right) Glycogen, Lactate, Phosphocreatine (PCr) in tissue; Plasma Lactate, Respiratory quotient (RQ); Glycogen in fibers, and ATP in type I and II fibers; The lines represent model simulations of whole muscle (—), type I fibers (---) and type II fibers (. . .). Experimental data are means ± standard error. Solid squares for whole tissue; upper triangles for type I fibers; diamonds for type II fibers.

Equations

The equations for this model may be viewed by opening up the Matlab model files downloaded from below.

Download Matlab model and associated files.

• Matlab files (zip file)3.01 MB
• Run Matlab simulation through 'SkelMusMet_Manu_V10_CON.m' file after extracting zip file.