<title>Introduction</title> <p>Skeletal muscle plays a major role in the regulation of whole-body substrates and energy metabolism, especially under changing physiological conditions such as ischemia (reduced blood flow), hypoxia (reduced oxygen supply), and exercise (increased energy demand). Current experimental techniques provide relatively little <italic>in vivo</italic> data on dynamic responses of muscle metabolite concentrations and metabolic fluxes to such physiological stimuli, especially in subcellular domains, such as mitochondria. To quantitatively analyze available <italic>in vivo</italic> experimental data and predict nonmeasurable dynamic responses, we developed a physiologically-based, multi-scale computational model of skeletal muscle cellular metabolism and energetics. The model is developed here from our previous model of cellular metabolism and energetics in skeletal muscle ##REF##18390321##[1]## and incorporates inter-domain transport processes and compartmentalized metabolic reactions of many key chemical species in both cytosol and mitochondria.</p> <p>Developing a mechanistic computational model of substrates and energy metabolism in complex, multi-scale metabolic systems, such as skeletal muscle, using a detailed, bottom-up systems approach with sparse <italic>in vivo</italic> experimental data—with an objective of achieving a quantitative understanding of the system to physiological perturbations—is a challenging task. Such a modeling approach requires information about the general structural features and catalytic mechanisms of the associated transporters and enzymes, subcellular metabolic pathways and fluxes and their control mechanisms, and tissue/organ specific metabolic characteristics. Such a modeling approach also requires mechanistic models for key functional components of the system (e.g., inter-domain transport processes, glycolysis, TCA cycle, oxidative phosphorylation, fatty acid β-oxidation) to be first individually developed and validated and then integrated to emulate the systems behavior at the molecular, subcellular, cellular, and tissue/organ levels. To avoid this complex approach and facilitate analysis of available sparse <italic>in vivo</italic> experimental data to understand dynamic responses of the system to physiological stresses, approximations are often made to obtain a simplified model of the system that includes key functional components regulating cellular metabolic processes at the desired level of complexity.</p> <p>A top-down systems approach is an alternative approach which has been previously applied to determine and integrate a representative set of lumped biochemical reactions <italic>in vivo</italic> metabolic systems that incorporate primary substrates and key intermediate metabolites with coupled metabolic energy controllers ATP-ADP and NADH-NAD<sup>+</sup>
##REF##18390321##[1]##–##REF##16603683##[8]##. This approach is similar to the top-down systems approach in metabolic control analysis proposed by Brand and co-workers ##REF##10583367##[9]##, ##REF##8944168##[10]## and is intended to provide an essential or minimal set of stoichiometrically balanced lumped biochemical reactions participating in ATP synthesis within mitochondria from the metabolism of nutrients (e.g., glucose, fatty acids, amino acids). Even with such simplifications, a large number of phenomenological kinetic parameters are introduced in the governing model equations, which must be estimated from available sparse <italic>in vivo</italic> experimental data. To estimate these unknown parameters, constraint-based, robust nonlinear optimization methods are needed, as established in our previous work ##REF##18390321##[1]##.</p> <p>To date, no physiologically-based, whole-organ level model of skeletal muscle cellular metabolism and energetics has been developed that can be applied to analyze available <italic>in vivo</italic> experimental data and predict dynamic metabolic responses to physiological stimuli in subcellular compartments, such as mitochondria. Previous models have incorporated some aspects of glycolysis, TCA cycle, oxidative phosphorylation, and fatty acid β-oxidation ##REF##10355547##[3]##–##REF##17448235##[5]##, ##REF##16636861##[11]##–##REF##17591785##[21]##. None of them, however, include sufficient key substrates/metabolites and/or integrate metabolic pathways/reactions that are essential in the regulation of cellular metabolic processes in skeletal muscle <italic>in vivo</italic> at whole-organ level. While our recently developed models of skeletal muscle cellular metabolism and energetics ##REF##18390321##[1]## have incorporated key metabolic pathways and reactions, the intracellular cytosolic and mitochondrial compartments were not distinguished. The energy controller pairs ATP-ADP and NADH-NAD<sup>+</sup> that modulate several key metabolic reactions in the cytosol and mitochondria have different concentrations in these two subcellular domains ##REF##15681693##[6]##–##REF##16603683##[8]##. As a consequence, the model may not accurately predict the dynamics of several metabolite concentrations and metabolic fluxes that are critical in the regulation of fuel (carbohydrate, fat, and lactate) metabolism and cellular respiration during physiological stresses such as ischemia, hypoxia, and exercise.</p> <p>In this paper, a physiologically-based, whole-organ level model of skeletal muscle cellular metabolism and energetics is developed and applied to study dynamic cellular metabolic responses to reduced blood flow and oxygen supply to mitochondria (muscle ischemia). The model, which is extended from our previous model ##REF##18390321##[1]##, is based on a multi-scale, top-down systems approach ##REF##18390321##[1]##–##REF##16603683##[8]##, accounts for subcellular compartmentalization, and includes primary substrates (carbohydrates and fats) and key intermediate metabolites and metabolic reactions specific to skeletal muscle metabolic system. The model equations are based on dynamic mass balances of chemical species in capillary blood and tissue cells (cytosol and mitochondria) domains. The model also distinguishes the free and bound forms of O<sub>2</sub> and CO<sub>2</sub> transport in the blood and cells. The inter-domain species transport processes are considered either by passive diffusion or by carrier-mediated (facilitated) transport. The metabolic reaction fluxes in the cytosolic and mitochondrial domains are represented by a general phenomenological Michaelis-Menten equation involving the compartmentalized ATP/ADP and NADH/NAD<sup>+</sup> energy controller ratios. The phenomenological kinetic parameters of the model are estimated by using our recently developed constraint-based, robust nonlinear optimization approach ##REF##18390321##[1]##. In this estimation process, we fit the model output to sparse <italic>in vivo</italic> dynamic data on glycolytic and energetic metabolite concentrations from experiments in humans with muscle ischemia previously published ##UREF##3##[22]##. With the estimated optimal parameter values, the model is able to simulate dynamic responses of key chemical species and reaction fluxes to reduced blood flow and oxygen supply to mitochondria associated with the muscle ischemia.</p>
<title>Materials and Methods</title> <title>Model Development</title> <p>The first step in our development of a multi-scale, top-down computational model of cellular metabolism and energetics in skeletal muscle was to identify the key intermediate metabolites and regulatory enzymes in the cellular metabolic pathways of skeletal muscle. We then integrated available information on cellular metabolic pathways and fluxes, cellular metabolic control mechanisms, catalytic enzyme kinetic mechanisms, subcellular compartmentation and metabolites volumes of distributions, inter-domain transport mechanisms, and skeletal muscle tissue-specific metabolic characteristics. A simplified map of the compartmentalized cellular metabolic pathways in skeletal muscle is shown in ##FIG##0##Figure 1##. The lumped biochemical reactions in the metabolic pathways are generated by stoichiometrically coupling several sequential elementary reactions. These reactions include the compartmentalized metabolic energy controller pairs ATP-ADP and NADH-NAD<sup>+</sup> whose ratios are known to modulate (fine tune) the reaction fluxes ##REF##9785955##[2]##, ##UREF##4##[23]## in the subcellular compartments. Many of these lumped reactions are considered irreversible for which the resting Gibbs free energy (Δ<italic>G</italic>) is high and negative in favor of product formation ##UREF##5##[24]##. As a part of a general formalism for modeling <italic>in vivo</italic> metabolic systems (nonequilibrium open systems), the reversible reactions like lactate dehydrogenase (LDH), creatine kinase (CK), and adenylate kinase (AK) were decomposed into two separate irreversible reactions with distinct kinetics ##REF##18390321##[1]##.</p> <title>Dynamic mass balance equations</title> <p>The dynamic mass balance equations are based on a multi-domain model structure for skeletal muscle consisting of a spatially-lumped capillary blood domain which exchanges nutrients and metabolic waste products with a spatially-lumped domain of tissue cells (##FIG##1##Figure 2##). Although these two domains are separated by the interstitial fluid (ISF) space, we assume phase-equilibrium of chemical species between the blood and ISF domains, and consider them together as the “blood” domain. Furthermore, the tissue cells domain is compartmentalized into the cytosolic and mitochondrial domains. The chemical species are assumed to be distributed in these two subcellular domains as per their mass fractions and volumes of distributions (##TAB##0##Table 1##). Here, the mass fraction of a species in a particular domain is defined as the fractional amount of the species in that domain in comparison to the total amount of the species in the whole muscle tissue cells. The volume of distribution of a species in a particular domain is defined as the anatomical volume of the domain plus the binding space of the domain for the species. In addition, the species that are common to both of these domains are assumed to have the similar dynamics in these two domains, because the species transport processes between these two domains can be sufficiently fast ##REF##17448235##[5]##. Consequently, a change in the species concentration in one domain will be proportional to the change in the species concentration in the other domain: <italic>C<sub>mit,j</sub></italic>(t) = <italic>σ<sub>j</sub>.C<sub>cyt,j</sub></italic>(t), where <italic>σ<sub>j</sub></italic> is the equilibrium concentration ratio (or partition coefficient) of species j between mitochondria and cytosol. In the present model, a total of 10 chemical species (pyruvate, fatty acyl-CoA, CoA, ATP, ADP, inorganic phosphate, NADH, NAD<sup>+</sup>, O<sub>2</sub>, and CO<sub>2</sub>) are considered to exist in both the cytosolic and mitochondrial compartments with negligible transport flux between the compartments.</p> <p>The dynamic mass balance of a chemical species <italic>j</italic> in the spatially-lumped blood domain has the following general form:where <italic>C<sub>art,j</sub></italic> is the arterial species concentration; <italic>C<sub>bl,j</sub></italic> is the capillary blood species concentration (equal to the venous species concentration <italic>C<sub>ven,j</sub></italic>); <italic>V<sub>bl,j</sub></italic> and <italic>V<sub>isf,j</sub></italic> are the volumes of distribution of species <italic>j</italic> in blood and ISF, and <italic>Q</italic> is the regional blood flow; <italic>J<sub>bl↔cyt,j</sub></italic> is the net transport flux (mass per unit time) across the blood-cytosol exchange barrier (consisting of capillary membrane, ISF, and tissue cell membrane).</p> <p>The dynamic mass balance of the chemical species j in the spatially-lumped tissue cells domain (cytosol and/or mitochondria) has the following general forms:

where
</p> <p>Here <italic>C<sub>x,j</sub></italic> is the species concentration in domain <italic>x</italic> (cytosol/mitochondria); <italic>V<sub>x,j</sub></italic> is the volume of distribution of species <italic>j</italic> in domain <italic>x</italic>; <italic>P<sub>x,j</sub></italic> and <italic>U<sub>x,j</sub></italic> are the production and utilization of species <italic>j</italic> in domain <italic>x</italic>; <italic>φ<sub>x,p</sub></italic> and <italic>φ<sub>x,u</sub></italic> are the reaction fluxes of the reactions processes that produce and utilize species <italic>j</italic> in domain <italic>x</italic>; <italic>β<sub>x,j,p</sub></italic> and <italic>β<sub>x,j,u</sub></italic> are the corresponding stoichiometric coefficients. For chemical species which are in the tissue cells (cytosol and/or mitochondria) but not in the capillary blood, the transport flux <italic>J<sub>bl↔cyt,j</sub></italic> is zero.</p> <p>The dynamic mass balance equations for all the chemical species in blood, cytosol and mitochondrial domains can be rewritten from our previous model of skeletal muscle metabolism ##REF##18390321##[1]## by accounting for the species compartmentalized volumes of distributions as laid out in Eqs. (1) and (2a–2c). The dynamic mass balance equations for O<sub>2</sub> and CO<sub>2</sub> in these domains are developed by considering their distinct transport and binding mechanisms ##REF##15682524##[25]##, ##REF##16775761##[26]##. The detailed mass balance equations of the chemical species, including O<sub>2</sub> and CO<sub>2</sub>, are provided in ##SUPPL##0##Materials S1##.</p> <title>Transport and reaction flux equations</title> <p>The reversible transport flux <italic>J<sub>bl↔cyt,j</sub></italic> (mass per unit time) of the species <italic>j</italic> across the blood-cytosol exchange barrier is related to the concentrations <italic>C<sub>bl,j</sub></italic> and <italic>C<sub>cyt,j</sub></italic> bywhere <italic>λ<sub>bl↔cyt,j</sub></italic> is the effective permeability surface area product for diffusive mass transport across the barrier (for passive transport); <italic>T<sub>bl↔cyt,j</sub></italic> is the maximal transport flux across the barrier (<italic>T</italic>
<sub>max</sub>) and <italic>M<sub>bl↔cyt,j</sub></italic> is the corresponding Michaelis-Menten (M-M) constant (<italic>M<sub>m</sub></italic>) (for facilitated transport). The flux expressions (3) satisfy the thermodynamic equilibrium conditions across the barrier.</p> <p>The species involved in the blood-cytosol exchange are glucose, lactate, pyruvate, alanine, glycerol, free fatty acid, CO<sub>2</sub>, and O<sub>2</sub> (##TAB##1##Table 2##). The transport processes may be passive or carrier-mediated (facilitated) (Eq. 3). The transport of glucose, pyruvate/lactate, and free fatty acid across the sarcolemma of skeletal muscle is facilitated via the GLUT4 ##REF##10413738##[27]##, MCT1 or MCT4 ##REF##11820324##[28]##, ##REF##10358105##[29]##, and FABP or FAT/CD36 ##REF##16198626##[30]## proteins, respectively. The sarcolemmal transport of the remaining species (alanine, glycerol, CO<sub>2</sub>, and O<sub>2</sub>) is considered to be passive. For the 10 species that exist in both cytosol and mitochondria, the inter-domain transport flux <italic>J<sub>cyt↔mit,j</sub></italic> is zero.</p> <p>The lumped metabolic reactions of skeletal muscle cellular metabolism and energetics can be considered as special cases of a general, irreversible, multi-reactant multi-product enzymatic reaction coupled with the metabolic energy controller pairs: where P<sub>1</sub> and P<sub>2</sub> are ATP and ADP or vice-versa (PS±: phosphorylation state); R<sub>1</sub> and R<sub>2</sub> are NADH and NAD<sup>+</sup> and vice-versa (RS±: redox state). The corresponding coupled phosphorylation and redox reactions are ATP→ADP (PS+) and/or NADH→NAD<sup>+</sup> (RS+) or vice-versa (PS- and RS-). Assuming a phenomenological, single-step enzyme kinetic mechanism ##UREF##6##[31]##, the flux expression for the lumped metabolic reaction can be written as (see Ref. ##REF##18390321##[1]## for detailed description):where <italic>V</italic>
<sub>max</sub> and <italic>K<sub>m</sub></italic> are the phenomenological maximal/limiting velocity and M-M parameter of the reaction;  are the phenomenological M-M parameters for the coupled phosphorylation and redox reactions; C<sub>P1</sub>/C<sub>P2</sub> is the phosphorylation ratio: C<sub>ATP</sub>/C<sub>ADP</sub> (PS+) or C<sub>ADP</sub>/C<sub>ATP</sub> (PS-), and C<sub>R1</sub>/C<sub>R2</sub> is the redox ratio: C<sub>NADH</sub>/C<sub>NAD+</sub> (RS+) or C<sub>NAD+</sub>/C<sub>NADH</sub> (RS-). So the reaction fluxes are explicitly regulated by the coupled energy controller ratios C<sub>P1</sub>/C<sub>P2</sub> and/or C<sub>R1</sub>/C<sub>R2</sub> (if coupled).</p> <p>The flux expressions for all the lumped metabolic reactions in the subcellular compartments of cytosol and mitochondria can be rewritten from our previous model of skeletal muscle metabolism ##REF##18390321##[1]## in terms of the compartmentalized metabolites concentrations and energy controller ratios ATP/ADP and NADH/NAD<sup>+</sup>. The details are provided in ##SUPPL##1##Materials S2##.</p> <p>In summary, in this mathematical model of skeletal muscle cellular metabolism and energetics, the number of chemical species (primary substrates, intermediate metabolites, and energy controllers) in the tissue cells domain (cytosol and mitochondria) is 30, which participate in 26 metabolic reactions. The number of chemical species that exist in both blood and cells is 8, and that in both cytosol and mitochondria is 10. Accordingly, the model includes 8 <italic>J<sub>bl↔cyt,j</sub></italic> transport fluxes characterized by 4 <italic>λ</italic>, 4 <italic>T</italic>
<sub>max</sub> and 4 <italic>M</italic>
<sub>m</sub> parameters, and 26 <italic>φ</italic> reaction fluxes characterized by 26 <italic>V</italic>
<sub>max</sub>, 27 <italic>K</italic>
<sub>m</sub> (1 <italic>K</italic>
<sub>m</sub> for G6P inhibition of the hexokinase reaction; see Reaction 1, ##SUPPL##1##Materials S2##),  parameters (as in our previous model ##REF##18390321##[1]##). Besides, there are 10 partition coefficients <italic>σ</italic> for the chemical species that exists in both cytosol and mitochondria (since the 10 <italic>J<sub>cyt↔mit,j</sub></italic> transport fluxes are considered negligible). Therefore, the present model is characterized by a total of 99 unknown parameters that need to be estimated to fit the model outputs to the available <italic>in vivo</italic> experimental data. For convenience, the transport and reaction fluxes are written here in vector form: , where the parameter vector for these transport and reaction fluxes is denoted by ;  is the concentration vector. Note that the 12 transport parameters (4 <italic>λ</italic>, 4 <italic>T</italic>
<sub>max</sub> and 4 <italic>M</italic>
<sub>m</sub>) can be estimated here uniquely from the resting, steady-state transport flux-concentration relationships (##TAB##2##Table 3##) ##REF##18390321##[1]##. The 10 partition coefficients (10 <italic>σ</italic>) can be estimated uniquely from the resting, steady-state species concentration ratios between the mitochondria and cytosol (<italic>σ<sub>j</sub></italic> = <italic>C<sub>mit,j</sub></italic>/<italic>C<sub>cyt,j</sub></italic>) (##TAB##2##Table 3##).</p> <title>Model simulation</title> <p>With specified parameter values, the mathematical model is solved numerically to simulate dynamic responses of the system to ischemia produced by reducing muscle blood flow. Typically, the initial conditions:  are assumed to be at a normal, resting steady-state. These are fixed based on average resting species concentrations gathered from various literature sources on skeletal muscle cellular metabolism that are consistent with the resting, steady-state flux-concentration relationships (##TAB##0##Tables 1## and ##TAB##1##2##). The species mass fractions and volumes of distributions in the subcellular cytosolic and mitochondrial domains are set to have appropriate phosphorylation and redox potentials in the cytosol and mitochondria. For numerical solution of this stiff initial-value problem, a robust implicit integrator DLSODES (<ext-link ext-link-type="uri" xlink:href="https://computation.llnl.gov/casc/odepack/odepack_home.html">https://computation.llnl.gov/casc/odepack/odepack_home.html</ext-link>; <ext-link ext-link-type="uri" xlink:href="http://www.netlib.org/odepack">http://www.netlib.org/odepack</ext-link>; ##UREF##7##[32]##) is used. Specifically, the Gear's implicit integration method based on backward difference formula (BDF) is most suitable for this problem. An absolute and relative error of tolerance of 10<sup>−10</sup> guarantees high accuracy and convergence of the iterative solutions of the ODEs. The DLSODES solver is usually very fast; a typical simulation of this problem using the DLSODES solver in a standard desktop computer (Intel Xeon or Core 2 Duo CPU 5160 @ 3 GHz) takes only about 5 seconds of the CPU time. As a check, the numerical solutions of the initial value problem were also obtained using the ODE15S solver in MATLAB (<ext-link ext-link-type="uri" xlink:href="http://www.mathworks.com">http://www.mathworks.com</ext-link>) with the similar tolerance levels. These solutions were of comparable accuracy with that obtained using the DLSODES solver.</p> <title>Parameter Estimation</title> <p>The large number of unknown parameters of this model are estimated by comparing model outputs to <italic>in vivo</italic> experimental data of Katz ##UREF##3##[22]## with the optimization procedure described previously ##REF##18390321##[1]##. The data consist of key metabolites concentration dynamics measured during circulatory occlusion and recovery (reperfusion) in human skeletal muscle at the whole tissue-organ level ##UREF##3##[22]##. Specifically, the data is based on the biopsy measurements of glucose and glycolytic/glycogenolytic intermediates (i.e., glucose 6-phosphate, pyruvate, and lactate) and creatine and high-energy phosphates (i.e., phosphocreatine, inorganic phosphate, ATP, ADP, and AMP) in quadriceps femoris muscle tissue at rest, after 30 minutes of ischemia, and after 15 minutes of reperfusion. The tissue metabolite contents or concentrations were measured in the units of mmol/kg dry weight. For our analysis, these concentrations are converted to the units of mmol/kg wet weight (mmol/L or mM) by multiplying a conversion factor of 0.25 kg dry weight/kg wet weight, corresponding to the muscle tissue ##UREF##8##[33]##. Consistent with the experimental measurements of Katz ##UREF##3##[22]##, the following ischemia protocol is used for model simulations and parameter estimation: the muscle blood flow is reduced from <italic>Q</italic> = 0.9 L/min at rest (t&lt;0 min) to <italic>Q</italic> = <italic>Q</italic>
<sub>isch</sub> (unknown) at the onset of ischemia (0≤t≤30 min), and then increased to the starting level at the onset of recovery (t&gt;30 min). The active muscle volume (<italic>V</italic>
<sub>mus</sub>) and ischemic muscle blood flow (<italic>Q</italic>
<sub>isch</sub>) responsible for the predicted metabolic dynamics in skeletal muscle during ischemia and recovery were not known from the experimental study of Katz ##UREF##3##[22]##. Therefore, <italic>V</italic>
<sub>mus</sub> and <italic>Q</italic>
<sub>isch</sub> are also included as the unknown parameters for estimation, making the total number of unknown parameters in the model to 101. This is an increase of 10 unknown parameters (10 <italic>σ</italic>) from our previous model ##REF##18390321##[1]##.</p> <p>With this relatively sparse data and large number of unknown parameters, the parameter estimation problem is ill-conditioned and under-determined. Nevertheless, an efficient estimation method is devised based on our earlier work ##REF##18390321##[1]## to obtain physiologically reasonable parameter values that minimize the sum of squared differences between the available experimental data and corresponding model outputs. The estimation procedure proceeds in two main stages. In the first stage, the published normal, resting species concentration data (##TAB##0##Tables 1## and ##TAB##1##2##) are used to evaluate the resting transport and reaction fluxes from a steady-state flux balance analysis (##TAB##1##Tables 2##, ##TAB##2##3## and ##TAB##3##4##). From the resting flux-concentration relationships: Eqs. (3) and (5), the preliminary estimates of the transport and reaction parameters are obtained. The preliminary estimates of the 10 partition coefficients are obtained from the resting, steady-state species concentration ratios between mitochondria and cytosol (<italic>σ<sub>j</sub></italic> = <italic>C<sub>mit,j</sub></italic>/<italic>C<sub>cyt,j</sub></italic>) (##TAB##2##Table 3##). In the second stage, the dynamic species concentration data ##UREF##3##[22]## together with the resting, steady-state flux balance equations as equality constraints are used to obtain the optimal parameter estimates (##TAB##2##Tables 3## and ##TAB##3##4##).</p> <p>A detailed description of this robust estimation approach for obtaining optimal parameter estimates from species concentration dynamics during muscle ischemia and recovery is presented in Ref. ##REF##18390321##[1]##. This method has been established powerful in the analysis of large-scale <italic>in vivo</italic> metabolic systems. The efficiency and robustness of this parameter estimation approach was tested with parameter sensitivity analysis as well as with repeated parameter estimation with various initial parameter estimates (initial guesses). Various estimates of the more sensitive model parameters spanned in a small neighborhood of the optimal parameter estimates (see Ref. ##REF##18390321##[1]## for details). Since the preliminary estimates of the 12 transport parameters and 10 partition coefficients based on the resting species concentrations and transport fluxes were accurate enough, these 22 parameters were not re-estimated from the dynamic species concentration data. Therefore, as in our previous paper ##REF##18390321##[1]##, a total of 101–22 = 79 parameters (77 reaction parameters+<italic>V</italic>
<sub>mus</sub>+<italic>Q</italic>
<sub>isch</sub>) are effectively estimated from the dynamic data.</p>
<title>Results</title> <title>Comparison of model simulations with experimental data</title> <p>The optimal estimates of the transport parameters (<italic>λ</italic>,<italic>T</italic>
<sub>max</sub>,<italic>M<sub>m</sub></italic>), partition coefficients (<italic>σ</italic>), active muscle volume and ischemic muscle blood flow (<italic>V</italic>
<sub>mus</sub>, <italic>Q</italic>
<sub>isch</sub>), and reaction parameters  for the model are shown in ##TAB##2##Tables 3## and ##TAB##3##4##. These parameter estimates yield the best fit of the model outputs to the published experimental data ##UREF##3##[22]## with minimal residual errors and minimal objective function. The optimal estimates of <italic>V</italic>
<sub>mus</sub> and <italic>Q</italic>
<sub>isch</sub> corresponding to the experimental data were ∼4.0 L (∼16% of normal two-legs muscle volume of ∼25 L) and ∼0.216 L/min (∼76% reduction from normal, resting two-legs muscle blood flow of <italic>Q</italic> = 0.9 L/min). These optimal parameter values were used for model simulations during the resting, ischemia and recovery periods. The blood flow reduction levels of <italic>Q</italic>
<sub>isch</sub> = 0.18, 0.27 and 0.36 L/min (80%, 70% and 60%) were used for simulating severe to moderate to mild ischemic conditions.</p> <p>The correspondence of model simulations and experimental data is demonstrated through ##FIG##2##Figures 3## and ##FIG##3##4##. Specifically, shown are the concentration dynamics of 4 glycolytic metabolites (GLC, G6P, LAC, PYR) and 6 energy metabolites (PCR, CR, PI, ATP, ADP, AMP) in the muscle tissue cells (i.e., weighted volume averages of concentrations in cytosol and mitochondria) for which experimental data were available ##UREF##3##[22]## for four different blood flow reduction levels Q<sub>isch</sub> = 0.18, 0.216, 0.27 and 0.36 L/min. ##FIG##2##Figure 3## also includes the concentration dynamics of NADH and NAD<sup>+</sup> in the muscle tissue cells. The dynamics of compartmentalized cytosolic and mitochondrial phosphorylation and redox potentials (i.e., [ATP]/[ADP] and [NAD<sup>+</sup>]/[NADH] ratios) and cytosolic [LAC]/[PYR] and [PCR]/[CR] ratios are shown in ##FIG##4##Figure 5##. The level <italic>Q</italic>
<sub>isch</sub> = 0.216 L/min corresponds to the experimental data of Katz ##UREF##3##[22]## and the corresponding model simulations match to the data reasonably well within the experimental noise (##TAB##4##Table 5##). The experimental data are normalized here with respect to the resting metabolites concentrations (control). Furthermore, the individual metabolites responses are shown in separate plots in order to distinguish metabolite responses to different levels of blood flow reductions.</p> <p>The model simulations of cellular glucose (GLC), glucose-6-phosphate (G6P), pyruvate (PYR), and lactate (LAC) concentrations at the end of 76% ischemia and reperfusion are in close agreement with the experimental data (##FIG##2##Fig. 3(A–D)##, ##TAB##4##Table 5##). Cellular [G6P] increased quickly (exponentially) by ∼60% during ischemia and returned rapidly to its resting level at the onset of reperfusion. In contrast, the cellular [GLC] increased slowly (almost linearly) by ∼35% during ischemia and remained at an elevated level even after 30 minutes of reperfusion. Blood [GLC], however, decreased rapidly, but only by ∼10%, during ischemia and returned quickly to its baseline value during reperfusion (not shown). Furthermore, the model-predicted changes in the cellular [G6P] and [GLC] during the mild 60% and 70% blood flow reduction levels were not significant when compared to the changes during the high 76% and above (i.e., 80%) blood flow reduction levels.</p> <p>Both cellular [LAC] and [PYR] increased by ∼125% during 76% ischemia and almost returned to the resting levels after 30 minutes of reperfusion. However, the dynamic response of [PYR] was different from that of [LAC]. During ischemia, [PYR] first decreased and then increased, but during reperfusion, [PYR] first sharply increased and then rapidly decreased to the resting level. On the other hand, [LAC] slowly (almost linearly) increased during ischemia and slowly (almost linearly) decreased during reperfusion. Blood [LAC] had the similar dynamics as that of cellular [LAC] (not shown). These differential dynamics of cellular [LAC] and [PYR] characterize the dynamics of cellular or cytosolic [LAC]/[PYR] ratio. The dynamics of cytosolic [LAC]/[PYR] ratio and cytosolic and mitochondrial redox states ([NADH] and [NAD<sup>+</sup>]) and redox potentials ([NAD<sup>+</sup>]/[NADH] ratios) were all similar having biphasic behaviors during the ischemia and reperfusion periods (##FIG##2##Figs. 3(E,F)## and ##FIG##4##5(A–C)##). This is in contrast to our previous model (in which the cytosolic and mitochondrial compartments were lumped into a single tissue cells compartment ##REF##18390321##[1]##) predictions that the dynamic responses of cellular [LAC]/[PYR] and [NAD<sup>+</sup>]/[NADH] ratios to ischemia and reperfusion are distinct. Thus the present compartmentalized model is able to correctly simulate the dynamic responses of cytosolic and cellular [LAC]/[PYR] ratios as well as cytosolic and mitochondrial [NAD<sup>+</sup>]/[NADH] ratios. The cytosolic [LAC]/[PYR] ratio rapidly increased from a resting level of ∼16.4 to ∼26.8 at the onset of 76% ischemia, then quickly decreased almost to the baseline value (∼17.22) as ischemia progressed, in agreement with the data (##TAB##4##Table 5##). In contrast, the cytosolic [NAD<sup>+</sup>]/[NADH] ratio decreased quickly from a resting level of ∼540 to ∼265 at the onset of 76% ischemia, then increased and reached a new steady state value of ∼332 as ischemia progressed. The LDH mass action ratio ([LAC][NAD<sup>+</sup>])/([PYR][NADH]) decreased considerably from a resting level of ∼8856 to ∼5585 during 76% ischemia. At different levels of blood flow reductions, the model-predicted changes were not proportional. The changes were negligible below 70% blood flow reduction levels and significant above 76% blood flow reduction levels.</p> <p>The model simulations of muscle phosphocreatine (PCR), creatine (CR), and inorganic phosphate (PI) concentrations at the end of 76% ischemia and recovery are in good agreement with the experimental data (##FIG##3##Fig. 4(A–C)##). After 30 minutes, 76% ischemia resulted in ∼17.5% decrease in muscle PCR content, which was fully resynthesized after 30 minutes of recovery. Muscle CR and PI contents increased by ∼35% and ∼130%, respectively, at the end of 76% ischemia and returned to their resting levels at the end of recovery. The dynamics of [PCR] drop during ischemia and rise during recovery were exponential, but faster during the recovery period. The dynamics of both [CR] and [PI] show the opposite trends (both are mirror images of [PCR] having the similar time constants). The cytosolic [PCR]/[CR] ratio decreased substantially from a resting value of ∼2 to ∼1.2 during ischemia (##FIG##4##Fig. 5D##, ##TAB##4##Table 5##). The model-predicted muscle [PCR], [CR], [PI] and [PCR]/[CR] ratio during the mild 60% and 70% blood flow reduction levels did not change appreciably from their baseline, resting levels.</p> <p>The model-simulated muscle [ATP] decreased slightly during ischemia which resulted in appropriate increases in muscle [ADP] and [AMP] (##FIG##3##Fig. 4(D–F)##). A 76% blood flow reduction resulted in ∼5% and ∼25% increases in [ADP] and [AMP], respectively, in accordance with the data. The cytosolic and mitochondrial phosphorylation potentials ([ATP]/[ADP] ratios) had the similar dynamics as of the cytosolic [PCR]/[CR] ratio (##FIG##4##Fig. 5(D–F)##). However, the magnitude of [ATP]/[ADP] decrease (∼6%) during ischemia was negligible in comparison to the magnitude of [PCR]/[CR] decrease (∼40%). The CK mass action ratio ([CR][ATP])/([PCR][ADP]) increased from a resting value of ∼166 to ∼260 during the peak ischemia. The total adenylate nucleotide pool (TAN = [ATP]+[ADP]+[AMP]) was maintained at a constant level of ∼7.0 mM in accordance with the data (##TAB##4##Table 5##). The model-predicted changes in muscle [ATP], [ADP], [AMP], and [ATP]/[ADP] ratio during the mild 60% and 70% blood flow reduction conditions were negligible. Thus, in contrast to our previous model ##REF##18390321##[1]##, the present compartmentalized model is able to simulate the dynamic responses of cytosolic and cellular [PCR]/[CR] ratios as well as cytosolic and mitochondrial [ATP]/[ADP] ratios appropriately.</p> <title>Simulated dynamics of mitochondrial metabolites</title> <p>The model-simulated dynamic responses of mitochondrial FAC, ACoA, CIT, AKG, SCoA, SUC, MAL and OXA during the ischemia and recovery periods with blood flow reduction levels of 80%, 76%, 70% and 60% (Q<sub>isch</sub> = 0.18, 0.216, 0.27 and 0.36 L/min) are shown in ##FIG##5##Figure 6##. Experimental data were not available for these key mitochondrial and TCA cycle intermediate metabolites. In fact, due to the limitations in the available experimental techniques, most of these subcellular metabolites can not be measured conveniently in skeletal muscle tissue cells <italic>in vivo</italic> at the whole tissue-organ level. However, the present compartmentalized model is able to simulate the dynamic responses of these subcellular metabolites to physiological stresses such as muscle ischemia. Model simulations show that the changes in these metabolites concentrations during muscle ischemia are not simply proportional to the extent of blood flow reduction. The changes are negligible during the mild 60% and moderate 70% blood flow reduction levels, but significant above the severe 76% blood flow reduction level.</p> <p>The model-simulated dynamic responses of mitochondrial and TCA cycle intermediate metabolites to muscle ischemia differ from each other (##FIG##5##Figure 6##). During ischemia, the mitochondrial [ACoA], [SUC] and [MAL] were increased, while [CIT], [AKG] and [SCoA] were decreased. Furthermore, the accumulation of ACoA was very significant (∼45 fold increase in [ACoA] with 76% blood flow reduction), in spite of a large decrease in mitochondrial redox ratio [NAD<sup>+</sup>]/[NADH]. The accumulation of ACoA is primarily attributed to a biphasic increase in free [CoA] and a sharp 80-fold decrease in [OXA] (##FIG##5##Fig. 6H##) resulting in a mismatch in the reaction fluxes producing and consuming ACoA in the mitochondria. The decreases in [SCoA] and increases in [SUC] are in accordance with the increases in mitochondrial phosphorylation ratio [ADP]/[ATP]. However, the decreases in [CIT] and [AKG] and increases in [MAL] were not consistent with the increases in mitochondrial redox ratio [NAD<sup>+</sup>]/[NADH]. The TCA cycle reaction fluxes decreased from a resting level of ∼0.82 to ∼0.62 mmol/min (∼25%) with a 76% blood flow reduction (not shown).</p> <title>Simulated dynamics of fluxes of exchangeable substrates</title> <p>The model-simulated dynamic responses of transport and reaction rates associated with the blood-tissue cells exchangeable substrates (i.e., glucose, lactate, pyruvate, free fatty acid, CO<sub>2</sub>, and O<sub>2</sub>) during the resting steady-state (<italic>t</italic>&lt;0 min), ischemia (0≤<italic>t</italic>≤30 min), and recovery (<italic>t</italic>&gt;30 min) periods with a blood flow reduction level of 76% (<italic>Q</italic>
<sub>isch</sub> = 0.216 L/min) are shown in ##FIG##6##Figure 7##. The corresponding capillary blood (venous blood) and tissue cells levels of O<sub>2</sub> and CO<sub>2</sub> are shown in ##FIG##7##Figure 8##. The transport rates include the species net uptake/release rates (<italic>UR<sub>j</sub></italic> = <italic>Q</italic>(<italic>C<sub>art,j</sub></italic>−<italic>C<sub>ven,j</sub></italic>)) and net blood-tissue cells exchange rates . The reaction rates include the species net production/utilization rates (<italic>R<sub>cl,j</sub></italic> = <italic>P<sub>cl,j</sub></italic>−<italic>U<sub>cl,j</sub></italic>) in the tissue cells (cytosol and/or mitochondria). Most of these model predictions are similar to those from our previous model ##REF##18390321##[1]## which do not account for the intracellular compartmentalization. The magnitudes of changes with the present compartmentalized model, however, are lower and consistent with the lower estimate of the blood flow reduction level during muscle ischemia (present <italic>Q</italic>
<sub>isch</sub> = 0.216 vs. previous <italic>Q</italic>
<sub>isch</sub> = 0.135).</p> <p>
##FIG##6##Figure 7A## shows that the net muscle glucose uptake (<italic>UR<sub>GLC</sub></italic>) changes sharply (decreases/increases) at the onset of ischemia/recovery due to a step change (decrease/increase) in muscle blood flow, while the net glucose uptake by muscle tissue cells (<italic>J<sub>bl↔cyt,GLC</sub></italic>) remains fairly constant during ischemia and recovery. This results in a fast adjustment (decrease/increase) in the venous glucose concentration and glucose A–V difference and a rapid recovery in muscle net glucose uptake. The net glucose utilization (<italic>U<sub>cl,GLC</sub></italic>−<italic>P<sub>cl,GLC</sub></italic>) through the hexokinase reaction decreases/increases quickly during ischemia/recovery due to the increase/decrease of [G6P] and G6P inhibition of hexokinase reaction. This is reflected in the predicted elevated cellular [GLC] during recovery (##FIG##2##Fig. 3A##). Several inhibition mechanisms of hexokinase reaction by G6P were tested. A general uncompetitive inhibition mechanism, which effectively modifies the <italic>K<sub>m</sub></italic> as well as <italic>V</italic>
<sub>max</sub> of the reaction ##UREF##4##[23]##, ##UREF##6##[31]##, was found to be the most suitable for the model to fit to the data.</p> <p>The net pyruvate <italic>UR</italic> and <italic>J<sub>bl↔cyt</sub></italic> rates are closely matched where both decrease/increase during ischemia/recovery, except at the onset of recovery, where <italic>UR<sub>PYR</sub></italic> and <italic>J<sub>bl↔cyt,PYR</sub></italic> rates first sharply decrease and then exponentially increase to the baseline level (##FIG##6##Fig. 7B##, dashed and dotted lines). The pyruvate metabolism switches from a net utilization to a net production and a net uptake to a net release during ischemia and vice-versa during recovery (as depicted through the solid line in ##FIG##6##Fig. 7B##). Similar to pyruvate, the net <italic>UR</italic> and <italic>J<sub>bl↔cyt</sub></italic> rates of lactate were approximately matched, through the −<italic>J<sub>bl↔cyt,LAC</sub></italic> rate was slightly higher than the −<italic>UR<sub>LAC</sub></italic> rate during ischemia and vice-versa during recovery (##FIG##6##Fig. 7C##, dashed and dotted lines). However, the net lactate release (−<italic>UR<sub>LAC</sub></italic> and −<italic>J<sub>bl↔cyt,LAC</sub></italic>) first decreases when muscle blood flow is reduced and then increases steadily (linearly) slightly above the baseline value during the remaining ischemic period. During recovery, the opposite, but more pronounced, trend occurs. The net lactate production (<italic>P<sub>cl,LAC</sub></italic>−<italic>U<sub>cl,LAC</sub></italic>) increased quickly at the onset of ischemia and reached a new steady-state level within the first 10 minutes of ischemia. At the onset of recovery, the net lactate production decreased quickly to the baseline value (solid line, ##FIG##6##Fig. 7C##).</p> <p>
##FIG##6##Figure 7D## shows that the net <italic>UR</italic> and <italic>J<sub>bl↔cyt</sub></italic> rates of free fatty acid (FFA) decreased quickly at the onset of ischemia and then increased rapidly towards the baseline value during ischemia (biphasic behavior); the opposite trend occurred during the recovery period. The net utilization of FFA (<italic>U<sub>cl,FFA</sub></italic>−<italic>P<sub>cl,FFA</sub></italic>) had exactly the opposite behavior of <italic>UR<sub>FFA</sub></italic> and <italic>J<sub>bl↔cyt,FFA</sub></italic> rates. The dynamics of CO<sub>2</sub> release and FFA uptake were similar (dashed and dotted lines, ##FIG##6##Fig. 7(D,E)##). However, the CO<sub>2</sub> production was closely correlated to the O<sub>2</sub> consumption (solid lines, ##FIG##6##Fig. 7(E,F)##). The muscle O<sub>2</sub> uptake, cellular O<sub>2</sub> uptake, and cellular O<sub>2</sub> consumption were all reduced during the ischemia period and were closely matched (all reduced by ∼22% during 76% ischemia), except at the onset of ischemia and recovery, where the muscle O<sub>2</sub> uptake decreased/increased sharply with the step change (decrease/increase) in muscle blood flow and returned to the new steady state levels as ischemia progressed (##FIG##6##Fig. 7F##). The O<sub>2</sub> partial pressure (<italic>P</italic>
<sub>O2</sub>) in the muscle tissue cells decreased from 25 mmHg to ∼0.5 mmHg (i.e., about the <italic>K<sub>m</sub></italic> value for O<sub>2</sub> consumption), while the corresponding CO<sub>2</sub> partial pressure (<italic>P</italic>
<sub>CO2</sub>) increased from 46 mmHg to ∼52 mmHg during the 76% ischemia period (##FIG##7##Fig. 8(B,D)##). Accordingly, the total O<sub>2</sub> content in the capillary blood (venous blood) decreased from ∼6.6 mM (∼36 mmHg) to ∼0.4 mM (∼8 mmHg) and the total CO<sub>2</sub> content in the blood increased from ∼25.5 mM (∼43.6 mmHg) to ∼29.6 mM (∼51 mmHg) (##FIG##7##Fig. 8(A,C)##). These model simulations indicate that the muscle tissue cells were hypoxic during the 76% ischemia period; the reduced O<sub>2</sub> consumption was maintained by the reduced O<sub>2</sub> supply during ischemia.</p> <title>Simulated dynamics of fluxes associated with ACoA and ATP</title> <p>The model-simulated dynamic responses of metabolic fluxes associated with the production and utilization of ACoA and ATP during the resting steady-state (<italic>t</italic>&lt;0 min), ischemia (0≤<italic>t</italic>≤30 min) and recovery (<italic>t</italic>&gt;30 min) periods with a blood flow reduction level of 76% (<italic>Q</italic>
<sub>isch</sub> = 0.216 L/min) are shown in ##FIG##8##Figure 9##. Specifically, ##FIG##8##Figures 9A and 9B## show the relative contributions from carbohydrates (PYR) and fats (FAC) to the production of ACoA in mitochondria, and the dynamics of ACoA utilization to CIT production in the mitochondrial TCA cycle. ##FIG##8##Figures 9C and 9D## depict the utilization rate of ATP through ATP hydrolysis in cytosol and production rate of ATP through oxidative phosphorylation (O<sub>2</sub> consumption) in mitochondria and creatine kinase buffer reaction in cytosol.</p> <p>The production of mitochondrial ACoA from PYR decreased rapidly at the onset of ischemia and then increased relatively slowly to reach a new steady-state level as ischemia progressed (biphasic behavior). During recovery, the trend was opposite (##FIG##8##Fig. 9A, dashed line##). The production of mitochondrial ACoA from FAC varied relatively little during the entire ischemic and recovery periods (##FIG##8##Fig. 9B, dotted line##). An increase in total production of ACoA from PYR and FAC occurred at the beginning of ischemia due to a transient increase in the mitochondrial fatty acid β-oxidation. However, as ischemia progressed, ACoA production decreased and reached a new steady-state level significantly below its baseline value (##FIG##8##Fig. 9B, solid line##). In contrast, the utilization of ACoA to produce CIT decreased at the onset of ischemia and then increased to match the total ACoA production from PYR and FAC (##FIG##8##Fig. 9B, dashed-dotted line##). A mismatch between the total production and utilization of ACoA at the onset of ischemia is reflected in the predicted mitochondrial accumulation of ACoA during ischemia (##FIG##5##Fig. 6B##).</p> <p>The rate of ATP utilization in cytosol through ATP hydrolysis slowly decreased (exponentially) by about 20% during 76% ischemia and returned to its baseline level during the recovery period (##FIG##8##Fig. 9C, dashed line##). However, the rate of ATP synthesis in mitochondria through oxidative phosphorylation (O<sub>2</sub> utilization) rapidly decreased (as a step) by about 22% during 76% ischemia (##FIG##8##Fig. 9C, dotted line##) due to a rapid reduction in O<sub>2</sub> delivery to mitochondria by capillary blood flow (##FIG##7##Fig. 8B##), in spite of sufficient increases in the mitochondrial redox potential [NADH]/[NAD<sup>+</sup>] (##FIG##0##Fig. 5B##) and phosphorylation potential [ADP]/[ATP] (##FIG##4##Fig. 5F##) (the activators of the lumped oxidative phosphorylation reaction). This transient deficit in ATP (##FIG##8##Fig. 9D solid line##) was mostly matched by the ATP production from the creatine kinase ATP buffer reaction (##FIG##8##Fig. 9D, dashed-dotted line##). The ATP supply from glycolysis/glycogenolysis (ATP deficit – ATP supply from creatine kinase reaction) during ischemia and recovery was almost constant, indicating that glycolysis/glycogenolysis is not a major source of ATP supply during partial ischemia. The time to reach new steady states (response time) for the ATP hydrolysis and creatine kinase reactions were about 15 minutes, while that for the ATP synthesis (or O<sub>2</sub> consumption) reaction was about 5 minutes. Thus, in contrast to our previous model ##REF##18390321##[1]##, the present compartmentalized model is able to correctly simulate the dynamic responses of cytosolic and mitochondrial metabolic reaction fluxes that are responsible in the production and utilization of ACoA and ATP. Furthermore, the model is able to link substrate metabolism (carbohydrates vs. fat) to energy metabolism in skeletal muscle during physiological perturbations such as muscle ischemia.</p>
<title>Discussion</title> <title>Modeling framework for analysis of sparse <italic>in vivo</italic> experimental data</title> <p>Analysis of dynamic changes in cellular metabolism and energetics in skeletal muscle <italic>in vivo</italic> in response to reduced blood flow and oxygen supply to mitochondria (ischemia) is limited by available <italic>in vivo</italic> experimental data. Relatively few <italic>in vivo</italic> measurements can be obtained of metabolite concentrations and metabolic fluxes in subcellular compartments, such as mitochondria. As a complementary approach to experimental measurements and as a framework for quantitatively analyzing available <italic>in vivo</italic> data, a physiologically-based, whole-organ level model of skeletal muscle cellular metabolism and energetics is developed. This model emphasizes a multi-scale, top-down systems approach ##REF##18390321##[1]##–##REF##16603683##[8]## for quantitative understanding of metabolic processes at the molecular, subcellular and cellular levels in relation to the tissue-organ responses. For this purpose, integration of information was required from cellular metabolic pathways and fluxes of substrate and energy metabolism in cytosol and mitochondria, cellular metabolic control mechanisms, catalytic enzyme kinetic mechanisms, subcellular compartmentation and metabolites volumes of distribution, inter-domain transport mechanisms, and tissue-specific skeletal muscle metabolic characteristics. Since the extent to which enzyme activities and kinetic (Michaelis-Menten) constants determined from <italic>in vitro</italic> experimental studies are applicable to <italic>in vivo</italic> conditions is uncertain, our multi-scale, top-down modeling approach provides a unified, systematic framework for modeling and analysis of <italic>in vivo</italic> complex metabolic systems. A major improvement of the present model over our previous model ##REF##18390321##[1]## is the incorporation of distinct intracellular cytosolic and mitochondrial compartments. The present model distinguishes the cytosolic and mitochondrial [ATP]/[ADP] and [NADH]/[NAD<sup>+</sup>] ratios (major modulators of several key metabolic reactions in cytosol and mitochondria) and incorporates metabolic reactions and transport processes of key chemical species, including glucose, lactate, pyruvate, free fatty acid, oxygen and carbon dioxide (##FIG##0##Fig. 1##).</p> <p>The governing model equations are based on dynamic mass balances of chemical species in spatially-lumped, capillary blood-interstitial fluid domain and two distinct intracellular cytosolic and mitochondrial domains (##FIG##1##Fig. 2##). These include the dynamic mass balances of O<sub>2</sub> and CO<sub>2</sub> which are developed by accounting for their distinct transport and binding mechanisms (hemoglobin and myoglobin-mediated) in these three domains ##REF##15682524##[25]##, ##REF##16775761##[26]##. However, since not much information is available regarding the blood flow heterogeneity and fiber type distributions in skeletal muscle, these metabolic characteristics were not considered in the present study. Instead, different fiber types were lumped together to an effective active muscle volume with an effective blood flow supplying oxygen and other nutrients to the lumped muscle volume, which is sufficient to investigate the dynamic metabolic responses to muscle ischemia (also see discussions in Section 4.4). Furthermore, since the phenomenon of muscle acidification/alkalization (i.e., proton handling) under physiological stresses, such as muscle ischemia, is highly complex ##REF##16617075##[17]##, ##REF##17591785##[21]##, this was not accounted for in the present model.</p> <p>In this model, the compartmentalized lumped biochemical reactions along the cellular metabolic pathways are the stoichiometrically coupled sequential elementary reactions, which include the compartmentalized metabolic energy controller pairs ATP-ADP and NADH-NAD<sup>+</sup> whose ratios (phosphorylation and redox potentials) modulate (or fine tune) the reaction fluxes in the subcellular domains ##REF##9785955##[2]##, ##UREF##4##[23]##. The lumped reactions for which the resting Gibbs free energy (Δ<italic>G</italic>) is large and negative in favor of product formation are considered irreversible ##UREF##5##[24]##. The reversible reactions like lactate dehydrogenase (LDH), creatine kinase (CK), and adenylate kinase (AK) are decomposed into two separate irreversible reactions with distinct kinetics. This can be justifiable based on the fact that the <italic>in vivo</italic> metabolic systems are nonequilibrium open systems and the biochemical reactions usually operate far from true equilibrium under physiological stimuli (also see discussions in Section 4.3). Besides, the consideration of the law of microscopic reversibility (thermodynamics) requires a detailed knowledge of proton handling in cytosol and mitochondria ##REF##16617075##[17]##, ##REF##17591785##[21]##, which is not considered in this phenomenological modeling study. The metabolic reaction fluxes in the subcellular cytosolic and mitochondrial domains are expressed in terms of a general phenomenological Michaelis-Menten equation with the coupled controller factors involving the compartmentalized [ATP]/[ADP] and [NADH]/[NAD<sup>+</sup>] energy controller ratios. These flux expressions characterize the behavior of saturable enzyme kinetics and <italic>in vivo</italic> control mechanisms observed experimentally ##UREF##4##[23]##, ##UREF##6##[31]##. In comparison to the previous models in the literature ##REF##10355547##[3]##–##REF##17448235##[5]##, ##REF##16636861##[11]##–##REF##17591785##[21]##, the present model is a self-consistent, physiologically-based, multi-scale computational model dealing with cellular metabolism and energetics in skeletal muscle <italic>in vivo</italic> at the whole-organ level.</p> <title>Optimal parameter estimation for the large-scale metabolic model</title> <p>A major challenge of this work was to evaluate a large number of unknown model parameters (precisely a total of 101) with relatively sparse experimental data from <italic>in vivo</italic> studies on muscle ischemia ##UREF##3##[22]##. To deal with this under-determined, ill-conditioned problem, we had to devise a special optimization strategy (see Ref. ##REF##18390321##[1]## for details). Briefly, preliminary estimates of the model parameters were obtained based on resting, steady-state flux balance analysis together with the information on resting, steady-state metabolites blood and tissue cells (cytosol and mitochondria) concentrations, metabolites uptake-release (UR) rates or arterio-venous (AV) differences, and metabolic reaction flux rates gathered from various literature sources on skeletal muscle cellular metabolism in normal human subjects (##TAB##0##Tables 1## and ##TAB##1##2##). With these preliminary parameter estimates, the model could not simulate the dynamic cellular metabolic responses to muscle ischemia. However, these estimates provided a good starting point in an iterative process of optimization to obtain the optimal parameter estimates (##TAB##2##Tables 3## and ##TAB##3##4##) that yield the best least-squares fit of the model to the data ##UREF##3##[22]## (##FIG##2##Figs. 3## and ##FIG##3##4##). A good starting point is essential for optimization of large-scale problems with gradient-types of optimization algorithms, such as the generalized reduced gradient algorithm (GRG2) ##UREF##9##[34]##, used in the present study. The first step of obtaining the preliminary parameter estimates provided accurate optimal estimates of the 12 transport parameters (4 <italic>λ</italic>, 4 <italic>T</italic>
<sub>max</sub> and 4 <italic>M</italic>
<sub>m</sub>) and 10 partition coefficients (10 <italic>σ</italic>), resulting in a reduction of 22 kinetic parameters for estimation from the dynamic experimental data from muscle ischemia.</p> <p>The key to success of our parameter estimation approach for obtaining the remaining 79 kinetic parameters (second step) was the incorporation of constraint information to the maximum extent in the GRG2 optimization algorithm ##UREF##9##[34]##. These include (1) the bound constraints on kinetic parameters from experimental studies, (2) the nonlinear equality (physiological) constraints relating fluxes and concentrations at resting, steady-state, and (3) the nonlinear inequality (thermodynamic) constraints relating fluxes and Gibbs free energy for reversible reactions: <italic>φ</italic>
<sub>net</sub>.Δ<italic>G</italic>≤0 ##REF##12080101##[35]##–##REF##17206279##[37]##. Also, the bounded estimation of active muscle volume (<italic>V</italic>
<sub>mus</sub>) and ischemic blood flow (<italic>Q</italic>
<sub>isch</sub>) (not known from the experimental protocol of Katz ##UREF##3##[22]##) provided additional degrees of freedom for fitting the model to the data. Use of such constraints reduced the number of unknown model parameters for simultaneous estimation (e.g., the 26 <italic>V</italic>
<sub>max</sub> parameters could be explicitly estimated from the flux-concentration relationships at resting, steady-state). Such constraints also reduced uncertainty and enhanced accuracy in the estimated parameter values.</p> <p>Since the available experimental data ##UREF##3##[22]## were relatively sparse in comparison to the number of unknown model parameters, the parameter estimates for this under-determined, ill-conditioned problem were expected to be non-unique. The efficiency and robustness of this parameter estimation approach was tested previously ##REF##18390321##[1]## in detailed with repeated parameter estimations with various initial parameter estimates (guesses) and with parameter sensitivity analysis. Various parameter estimates of the more sensitivity model parameters spanned in a small neighborhood of the accepted optimal parameter estimates. Re-estimation of the model parameters for the present extended compartmentalized model with the same parameter estimation approach and the quality of the model fitting to the experimental data (##FIG##2##Figs 3## and ##FIG##3##4##) further testifies the robustness of the parameter estimation approach.</p> <title>Model predictions related to muscle ischemia</title> <p>This physiologically-based, multi-scale computational model of skeletal muscle cellular metabolism and energetics can simulate dynamic metabolic changes in cellular and subcellular compartments and provide quantitative analysis of the mechanisms of metabolic regulation during physiological stresses (e.g., reduced blood flow and oxygen supply to mitochondria associated with muscle ischemia). Model simulations based on the estimated optimal parameter values correspond well to the experimental data on muscle ischemia of Katz ##UREF##3##[22]## (##FIG##2##Figs. 3## and ##FIG##3##4##). The optimal estimates of active muscle volume <italic>V</italic>
<sub>mus</sub>≈4 L and ischemic muscle blood flow <italic>Q</italic>
<sub>isch</sub>≈0.216 L/min (##TAB##2##Table 3##) indicate that only ∼16% of whole muscle volume (∼25 L) was actively participating in cellular metabolism during reduced oxygen delivery to mitochondria induced by a blood flow reduction of ∼76% from the resting level of 0.9 L/min. Though these estimates are slightly different from those estimated previously ##REF##18390321##[1]## (<italic>V</italic>
<sub>mus</sub>≈5 L and <italic>Q</italic>
<sub>isch</sub>≈0.135 L/min), when scaled to one leg with a resting blood flow of 0.45 L/min, the estimated V<sub>mus</sub> value becomes ∼2 L which lies well within the experimental range of 1.5–3 L for one leg quadriceps femoris muscle ##REF##12456845##[38]##–##REF##9729599##[41]##. The experimental data ##UREF##3##[22]## used in this computational study were also based on measurements from similar muscle types. This indicates that the whole muscle is not uniformly perfused and that only a fraction of the muscle mass, in accord with metabolic demand, participates in cellular metabolic processes.</p> <p>The model simulations of cellular [GLC], [G6P], [LAC], [PYR], and [LAC]/[PYR] with ∼76% ischemia corresponded well with the experimental data (##FIG##2##Fig. 3##). To simulate [GLC] and [G6P] dynamics, and to some extent [LAC] and [PYR] dynamics in response to muscle ischemia, the model had to account for the inhibition of hexokinase reaction by [G6P] and activation of glycogen phosphorylase reaction by [AMP]/[ATP]. [G6P] increased during ischemia due to the increased [AMP]/[ATP] ratio and increased glycogenolysis flux. The increased [G6P] decreases the flux of hexokinase and glucose utilization. As a result, the glucose taken up by the muscle tissue cells is accumulated during ischemia.</p> <p>The model-simulated muscle [PCR], [CR], [PCR]/[CR], [PI], [ATP], [ADP], [ATP]/[ADP], and [AMP] with ∼76% ischemia were in close agreement with the experimental data (##FIG##3##Fig. 4##). One of the key aspects in achieving these successful model simulations and good fittings was to model each of the reversible reactions like creatine kinase (CK) and adenylate kinase (AK) as two separate irreversible reactions with distinct kinetics. The forward and reversible reaction fluxes were considered regulated by the phosphorylation potentials [ATP]/[ADP] or [ADP]/[ATP] in the cytosol (##SUPPL##1##Materials S2##). Without considering such control mechanisms, it was not possible to simulate the experimental data on these high-energy metabolites during muscle ischemia and reperfusion.</p> <p>We tested fitting the model to the data using alternative flux expressions for the CK and AK (also LDH) reversible reactions satisfying the Haldane relationship and thermodynamic equilibrium conditions. However, we were unable to achieve good fittings for the high-energy metabolites. This may be due to the limitations in the present model, or may be the fact that the data do not support thermodynamic equilibrium of the CK reaction during ischemia/recovery. It is evident from the data of Katz ##UREF##3##[22]## that the changes in muscle [PCR]/[CR] and [ATP]/[ADP] ratios during ischemia are not proportional. The muscle [PCR]/[CR] ratio decreased substantially from a resting level of ∼2.0 to ∼1.2 (∼40%) during ischemia (##FIG##4##Fig. 5D##), while the muscle [ATP]/[ADP] ratio decreased only little from a resting level of ∼7.5 to ∼7.0 (∼7%). However, the dynamics of cytosolic [ATP]/[ADP] ratio was not available from the experimental study, which determines the thermodynamic equilibrium for the CK reaction. With our assumption of a constant cytosolic to mitochondrial ratio of [ATP] and [ADP] (and other common chemical species) during ischemia/recovery, the model predicted a comparable 7% decrease in cytosolic [ATP]/[ADP] ratio from a resting level of ∼332 to ∼310 during the ischemia period (##FIG##4##Fig. 5E##). As a consequence, the CK mass action ratio ([CR][ATP])/([PCR][ADP]) increased from a resting value of ∼166 to ∼260 during ischemia, instead of being a constant (∼166) (the CK apparent equilibrium constant at pH = 7) ##REF##16617075##[17]##.</p> <p>This issue may be resolved by including the dynamics and buffering of cytosolic and mitochondrial pH ##REF##16617075##[17]##, ##REF##17591785##[21]## and, to some extent, by distinguishing the dynamics of common metabolites within the subcellular cytosolic and mitochondrial compartments. The later will necessitate incorporation of transport mechanisms of common metabolites across the mitochondrial inner membrane, e.g., proton pumps, ATP/ADP translocase, PI/H<sup>+</sup> cotransporter, and exchange of TCA cycle substrates/intermediates ##UREF##1##[18]##, ##REF##17591785##[21]##. Since ∼98% of ADP is in mitochondria and only ∼2% of ADP is in cytosol at normal, resting steady-state conditions, a small relative change in whole cell [ADP] may correspond to a huge relative change in cytosolic [ADP], and hence a huge relative change in cytosolic [ATP]/[ADP] ratio. This can lead to the displacement of thermodynamic equilibrium of ATP/ADP translocase without displacing the thermodynamic equilibrium of CK reaction. Note that the cytosolic [ATP]/[ADP] ratio needs to decrease proportionally to the cytosolic [PCR][H<sup>+</sup>]/[CR] ratio for maintaining the thermodynamic equilibrium of the CK reaction. This information will be the mechanistic basis for the extension of the current model which can be used to analyze the data on high-energy phosphate metabolites.</p> <p>With the estimated optimal parameter values, the present model is able to correctly simulate the dynamic metabolic responses of compartmentalized cytosolic and mitochondrial phosphorylation and redox potentials (##FIG##4##Fig. 5##), key mitochondrial and TCA cycle metabolite concentrations (##FIG##5##Fig. 6##), and other key metabolite concentrations and metabolic fluxes (##FIG##6##Figs. 7##–
##FIG##8##9##) in skeletal muscle during ischemia and reperfusion for which no experimental data are available. Some of these simulations may be questionable without further validation, such as the predictions of cytosolic and mitochondrial phosphorylation and redox potentials (##FIG##4##Fig. 5##). However, these predictions are consistent with the dynamic responses of cellular or cytosolic [PCR]/[CR] and [LAC]/[PYR] ratios. Some of these simulations could not be obtained accurately with our previous model ##REF##18390321##[1]##.</p> <p>Our compartmentalized model simulations show that the metabolic changes were negligible with respect to the extent of blood flow reduction up to ∼70%, beyond which a profound derangement in substrate and energy metabolism occurred due to reduced supply of oxygen to mitochondria. The metabolic changes during severe ischemia (∼76% and beyond) include (1) a switch from a net pyruvate uptake and utilization to a net pyruvate production and release (##FIG##2##Figs. 3C## and ##FIG##6##7B##); (2) an increase in the redox potential [NADH]/[NAD<sup>+</sup>] (##FIG##2##Figs. 3(E,F)## and ##FIG##4##5(A,B)##) accompanied by a net formation and accumulation of cytosolic LAC and PYR (##FIG##2##Fig. 3(C,D)##) and mitochondrial ACoA (##FIG##5##Figs. 6B## and ##FIG##8##9(A,B)##); (3) an increase in the glycolytic and creatine kinase contributions to ATP formation (##FIG##8##Fig. 9(C,D)##) accompanied by a net depletion of cytosolic PCR (ATP stores) (##FIG##3##Fig. 4A##). In general, these predictions are in agreement with experimental results from studies under similar conditions ##REF##6420545##[42]##–##REF##3170412##[44]##. However, since our model does not differentiate between FADH<sub>2</sub> and NADH, the less efficient fuel (i.e., fats) –from the oxygen consumption point of view– is considered to have the same oxygen cost as carbohydrates for ATP production.</p> <title>Limitations of the model</title> <p>In this multi-scale, multi-compartmental, top-down integrated model of cellular metabolism and energetics in skeletal muscle, the substrate and cation transport mechanisms between cytosol and mitochondria are not incorporated. Instead, these two domains are assumed to be in fast equilibrium with each other with the common metabolites having the similar dynamics (i.e., <italic>C<sub>mit,j</sub></italic> = <italic>σ<sub>j</sub>.C<sub>cyt,j</sub></italic>, with the partition coefficients <italic>σ<sub>j</sub></italic> constant). This is projected in the model-simulated dynamic responses of cytosolic and mitochondrial phosphorylation and redox potentials (i.e., [ATP]/[ADP] and [NADH]/[NAD<sup>+</sup>] ratios) (##FIG##4##Fig. 5##). This assumption may be reasonable as a first approximation, since not much information is available regarding the dynamic responses of common metabolites to physiological stimuli (e.g., muscle ischemia) in these two subcellular domains (cytosol and mitochondria). This assumption may also be reasonable as long as the common metabolite concentrations do not vary appreciably within these two subcellular domains during physiological stresses. However, this assumption may not be valid for the [ATP]/[ADP] and [NADH]/[NAD<sup>+</sup>] ratios during ischemia/recovery, which are major modulators of several key metabolic reactions in cytosol and mitochondria. Therefore, the model may not accurately predict the dynamics of metabolite concentrations and metabolic fluxes that are critical in the regulation of fuel (carbohydrate, fat, and lactate) metabolism and cellular respiration during physiological perturbations, such as ischemia, hypoxia, and exercise. Furthermore, the inclusion of FADH<sub>2</sub> as a reducing equivalent (produced through oxidation of fatty acids) is essential when investigating the relative contribution of carbohydrate and fat oxidations in the synthesis of ATP in mitochondria.</p> <p>Another limitation in the present model is that the model does not account for blood flow heterogeneity and fiber type distributions associated with skeletal muscle cellular metabolism. Specifically, the capillary blood domain is considered spatially-lumped and the advective-diffusive transport of chemical species in capillary blood is ignored. Such compartmental modeling may not account for situations where important gradients of chemical species exist (e.g., due to distinct metabolic characteristics associated with different fiber types). In skeletal muscle, where the blood flow is highly heterogeneous due to complex capillarization and fiber type distributions, more than 70% of arterial oxygen is extracted from the blood before it leaves the microcirculation ##UREF##11##[45]##. High extraction results in large intracellular and intravascular gradients of oxygen and other nutrients. Therefore, the well-mixed assumption may not be valid and compartmental modeling may not account for the heterogeneity of oxygen transport and consumption in different fiber types. These issues can be addressed through spatially-distributed modeling of oxygen transport and consumption, coupled with other substrates and CO<sub>2</sub>
##UREF##2##[19]##, ##REF##16775761##[26]##.</p> <p>Finally, in the present model most of the lumped biochemical reactions in the cellular metabolic pathways are considered irreversible in the direction of product formation. Although the flux in the reverse direction is almost one order of magnitude smaller than the flux in the forward direction under normal, resting steady-state conditions, their magnitudes may change during physiological perturbations such as ischemia, hypoxia or exercise. Therefore, it would be advantageous to consider all the biochemical reactions to be reversible and their fluxes to satisfy the thermodynamic equilibrium constraints ##REF##17206279##[37]##, including the dynamics and buffering of pH ##REF##16617075##[17]##, ##REF##17591785##[21]##, which can affect the apparent equilibrium constants and standard Gibbs free energy of the reactions.</p> <title>Future model developments and potential applications</title> <p>Including the substrate and cation transport mechanisms between cytosol and mitochondria and distinguishing common metabolites dynamics within these subcellular domains would be important in the analysis of dynamic cellular metabolic responses to high intensity exercise. The subcellular compartmentalization could account for distinct volumes of metabolites distribution, for example, metabolic channeling for glycolysis ##REF##15681693##[6]##–##REF##16603683##[8]##. This is a way to recognize that the key glycolytic enzymes are bound together to form a multi-enzyme complex near the sarcolemma and sarcoplasmic reticulum. To investigate the effects of spatial locality due to blood flow heterogeneity and fiber type distributions on metabolic responses to muscle ischemia, hypoxia or exercise, one will require spatially-distributed models involving parallally-arranged capillaries representing different muscle fiber types. Furthermore, the use of alternative kinetic flux expressions based on a Michaelis-Menten formalism for reversible enzymatic reactions would satisfy the thermodynamic equilibrium conditions and would provide additional thermodynamic constraints for kinetic parameter estimation ##REF##17206279##[37]##. As in our previous paper ##REF##18390321##[1]##, a formal dynamic parameter sensitivity analysis can be carried out to understand the relative importance of the model parameters on the measured outputs. This information can be used to improve confidence in the estimates of key model parameters by fixing the values of the least sensitive parameters (parameter space reduction approach).</p> <p>With these enhancements and with additional metabolic pathways similar to those of the models of cellular metabolism and energetics in cardiac muscle ##REF##15681693##[6]##–##REF##16603683##[8]##, such as the inclusion of FADH<sub>2</sub> as a reducing equivalent, the advanced model of cellular metabolism and energetics in skeletal muscle can be used to analyze physiological responses to exercise (increased energy demand). Specifically, the model can assist in providing insight into long-standing or paradoxical questions, such as the extent of control that mitochondrial oxygen concentration and redox state exert over the rate of lactate production in the cytosol during heavy intensity exercise. In general, a multi-scale computational model of cellular metabolism and energetics in tissue/organ systems within the body that incorporate sufficient mechanistic processes can be used to (1) analyze and interpret <italic>in vivo</italic> experimental data; (2) provide new hypotheses as well as suggest how to test a given hypothesis; (3) design critical experiments; (4) quantify and predict dynamic responses to physiological stimuli that can not be directly measured; (5) evaluate relative importance of metabolic pathways and fluxes and their regulatory mechanisms under both normal and pathological conditions; and (6) provide the mechanistic basis for simulating integrated effects of altering enzyme activities and substrate concentrations through pharmacological agents and/or dietary inputs.</p>
<p>Conceived and designed the experiments: RKD MEC. Performed the experiments: RKD YL JK. Analyzed the data: RKD YL JK GMS MEC. Contributed reagents/materials/analysis tools: DAB GMS. Wrote the paper: RKD GMS.</p> <p>Current address: Biotechnology and Bioengineering Center, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America</p> <p>Current address: Department of Physiology, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America</p> <p>Control mechanisms of cellular metabolism and energetics in skeletal muscle that may become evident in response to physiological stresses such as reduction in blood flow and oxygen supply to mitochondria can be quantitatively understood using a multi-scale computational model. The analysis of dynamic responses from such a model can provide insights into mechanisms of metabolic regulation that may not be evident from experimental studies. For the purpose, a physiologically-based, multi-scale computational model of skeletal muscle cellular metabolism and energetics was developed to describe dynamic responses of key chemical species and reaction fluxes to muscle ischemia. The model, which incorporates key transport and metabolic processes and subcellular compartmentalization, is based on dynamic mass balances of 30 chemical species in both capillary blood and tissue cells (cytosol and mitochondria) domains. The reaction fluxes in cytosol and mitochondria are expressed in terms of a general phenomenological Michaelis-Menten equation involving the compartmentalized energy controller ratios ATP/ADP and NADH/NAD<sup>+</sup>. The unknown transport and reaction parameters in the model are estimated simultaneously by minimizing the differences between available <italic>in vivo</italic> experimental data on muscle ischemia and corresponding model outputs in coupled with the resting linear flux balance constraints using a robust, nonlinear, constrained-based, reduced gradient optimization algorithm. With the optimal parameter values, the model is able to simulate dynamic responses to reduced blood flow and oxygen supply to mitochondria associated with muscle ischemia of several key metabolite concentrations and metabolic fluxes in the subcellular cytosolic and mitochondrial compartments, some that can be measured and others that can not be measured with the current experimental techniques. The model can be applied to test complex hypotheses involving dynamic regulation of cellular metabolism and energetics in skeletal muscle during physiological stresses such as ischemia, hypoxia, and exercise.</p>
<title>Supporting Information</title>
<p>The authors are thankful to the anonymous reviewers for their constructive criticisms and suggestions which helped in improving this manuscript significantly.</p>
<fig id="pone-0003168-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g001</object-id><label>Figure 1</label><caption><title>Schematic diagram of biochemical pathways depicting various chemical reactions and species involved in the cellular metabolism of skeletal muscle.</title><p>The pathways involve 26 lumped reactions among 30 species out of which 8 species (GLC, LAC, PYR, ALA, FFA, GLR, CO<sub>2</sub> and O<sub>2</sub>) undergo blood-tissue cells exchange. The exchange arrows show the direction of net tissue cells uptake-release rates at normal, resting conditions. The lumped reactions are further compartmentalized into the cytosolic and mitochondrial reactions. These two subcellular domains are assumed to be in rapid equilibrium state (the barrier is shown schematically by double dotted lines), so that the species that are common to both these domains can have the similar dynamics, i.e., <italic>C<sub>mit,j</sub></italic>(<italic>t</italic>) = <italic>σ<sub>j</sub>.C<sub>cyt,j</sub></italic>(<italic>t</italic>), where <italic>σ<sub>j</sub></italic> is the partition coefficient of the species j between cytosol and mitochondria. In this case, the net transport flux of a species j across the cytosol-mitochondria barrier can quickly become negligible at the onset of a physiological perturbation (the transport fluxes are shown by double dotted arrows); 10 species (PYR, FAC, CoA, NADH, NAD<sup>+</sup>, ATP, ADP, PI, CO<sub>2</sub> and O<sub>2</sub>) exist in both the cytosolic and mitochondrial domains. GLC: glucose, GLY: glycogen, G6P: glucose-6-phosphate, GA3P: glyceraldehyde-3-phosphate, 13BPG: 1,3-biphosphate-glycerate, PYR: pyruvate, LAC: lactate, ALA: alanine, TGL: triglycerides, GLR: glycerol, FFA: free fatty acid, FAC: fatty acyl-CoA, ACoA: acetyl-CoA, CIT: citrate, AKG: α-ketogluterate, SCoA: succinyl-CoA, SUC: succinate, MAL: malate, OXA: oxaloacetate, CoA: coenzyme-A (free), PCR: phosphocreatine, CR: creatine, PI: inorganic phosphate, CO<sub>2</sub>: carbon dioxide, O<sub>2</sub>: oxygen, NADH: reduced nicotinamide adenine dinucleotide, NAD<sup>+</sup>: oxidized nicotinamide adenine dinucleotide, ATP: adenosine triphosphate, ADP: adenosine diphosphate, AMP: adenosine monophosphate.</p></caption></fig> <fig id="pone-0003168-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g002</object-id><label>Figure 2</label><caption><title>Schematic diagram of the structure of the model for blood-tissue cells exchange and cellular metabolism in skeletal muscle.</title><p>The compartments are assumed to be perfectly mixed, and the capillary blood and tissue ISF regions are assumed to be in phase-equilibrium with each other, so that <italic>C<sub>isf,j</sub></italic> = <italic>C<sub>bl,j</sub></italic> = <italic>C<sub>ven,j</sub></italic> for any chemical species <italic>j</italic>. The tissue cells domain is further compartmentalized into the cytosolic and mitochondrial domains with the chemical species having the similar dynamics in these two subcellular domains, so that <italic>C<sub>mit,j</sub></italic>(<italic>t</italic>) = <italic>σ<sub>j</sub>C<sub>cyt,j</sub></italic>(<italic>t</italic>), where <italic>σ<sub>j</sub></italic> is the partition coefficient of the species j between cytosol and mitochondria. The model accounts for 30 chemical species in the tissue cells. A total of 8 species (GLC, LAC, PYR, ALA, FFA, GLR, CO<sub>2</sub> and O<sub>2</sub>) undergo blood-tissue cells exchange; 10 species (PYR, FAC, CoA, NADH, NAD<sup>+</sup>, ATP, ADP, PI, CO<sub>2</sub> and O<sub>2</sub>) exist in both the cytosolic and mitochondrial domains with a negligible transport flux (<italic>J<sub>cyt↔mit,j</sub></italic>≈0). For details, see the caption of ##FIG##0##Figure 1##.</p></caption></fig> <fig id="pone-0003168-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g003</object-id><label>Figure 3</label><caption><title>Model-predicted dynamic responses of glycolytic metabolite concentrations and redox states in the muscle tissue cells during the resting, ischemia and recovery periods with varying levels of blood flow reduction and their comparison to the experimental data of Katz</title><p>
##UREF##3##[22]##. The responses were computed using the estimated optimal parameter values with the ischemia protocol of −5 to 0 min of resting, 0 to 30 min of ischemia, and 30 to 60 min of recovery. The muscle blood flow <italic>Q</italic> is reduced as a step from 0.9 L/min at rest to <italic>Q</italic>
<sub>isch</sub> = 0.36, 0.27, 0.216, 0.18 L/min at the onset of ischemia and returned to 0.9 L/min at the onset of recovery. The lines represent the model simulation results with the symbols representing the experimental data points (mean±SD) corresponding to <italic>Q</italic>
<sub>isch</sub> = 0.216 L/min (∼76% blood flow reduction). The metabolites concentrations are shown in normalized form, normalized with respect to the resting metabolites concentrations. The concentrations in the tissue cells are calculated based on the formula: <italic>C<sub>cl</sub></italic> = (<italic>V<sub>cyt</sub>C<sub>cyt</sub></italic>+<italic>V<sub>mit</sub>C<sub>mit</sub></italic>)/<italic>V<sub>cl</sub></italic>.</p></caption></fig> <fig id="pone-0003168-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g004</object-id><label>Figure 4</label><caption><title>Model-predicted dynamic responses of energy metabolite concentrations in the muscle tissue cells during the resting, ischemia and recovery periods with varying levels of blood flow reductions and their comparison to the experimental data of Katz ##UREF##3##[22]##.</title><p>The lines represent the model simulation results with the symbols representing the experimental data points (mean±SD) corresponding to <italic>Q</italic>
<sub>isch</sub> = 0.216 L/min (∼76% blood flow reduction). The simulation strategy, ischemia protocol, and blood flow reduction levels are exactly the same as those described in the caption of ##FIG##2##Figure 3##. The metabolites concentrations are shown in normalized form, normalized with respect to the resting metabolites concentrations. The concentrations in the tissue cells were calculated based on the formula: <italic>C<sub>cl</sub></italic> = (<italic>V<sub>cyt</sub>C<sub>cyt</sub></italic>+<italic>V<sub>mit</sub>C<sub>mit</sub></italic>)/<italic>V<sub>cl</sub></italic>.</p></caption></fig> <fig id="pone-0003168-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g005</object-id><label>Figure 5</label><caption><title>Model-predicted dynamic responses of cytosolic [LAC]/[PYR] and [PCR]/[CR] ratios and cytosolic and mitochondrial [ATP]/[ADP] and [NAD<sup>+</sup>]/[NADH] ratios during the resting, ischemia and recovery periods with varying levels of blood flow reductions.</title><p>The simulation strategy, ischemia protocol, and blood flow reduction levels are exactly the same as those described in the caption of ##FIG##2##Figure 3##.</p></caption></fig> <fig id="pone-0003168-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g006</object-id><label>Figure 6</label><caption><title>Model-predicted dynamic responses of FAC, ACoA, and TCA cycle intermediate concentrations in the muscle tissue cells mitochondria during the resting, ischemia and recovery periods with varying levels of blood flow reduction.</title><p>The simulation strategy, ischemia protocol, and blood flow reduction levels are exactly the same as those described in the caption of ##FIG##2##Figure 3##. The metabolites concentrations are shown in normalized form, normalized with respect to the resting metabolites concentrations. The concentrations in the tissue cells were calculated based on the formula: <italic>C<sub>cl</sub></italic> = (<italic>V<sub>cyt</sub>C<sub>cyt</sub></italic>+<italic>V<sub>mit</sub>C<sub>mit</sub></italic>)/<italic>V<sub>cl</sub></italic>.</p></caption></fig> <fig id="pone-0003168-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g007</object-id><label>Figure 7</label><caption><title>Model-predicted dynamic responses of transport and metabolic fluxes of glucose, pyruvate, lactate, free fatty acid, oxygen, and carbon dioxide during the resting, ischemia and recovery periods with a blood flow reduction of about 76%.</title><p>The responses were computed using the estimated optimal parameter values with the ischemia protocol of −5 to 0 min of resting, 0 to 30 min of ischemia, and 30 to 60 min of recovery. The muscle blood flow Q is reduced as a step from 0.9 L/min at rest to Q<sub>isch</sub> = 0.216 L/min at the onset of ischemia and returned to 0.9 L/min at the onset of recovery. Q(C<sub>a</sub>–C<sub>v</sub>): uptake-release rates, J<sub>b↔c</sub>: blood-tissue cells transport rates, P<sub>c</sub>–U<sub>c</sub>: tissue cells production-utilization rates.</p></caption></fig> <fig id="pone-0003168-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g008</object-id><label>Figure 8</label><caption><title>Model-predicted dynamic responses of O<sub>2</sub> and CO<sub>2</sub> (total contents and partial pressures) in the capillary blood (or venous blood) and tissue cells during the resting, ischemia and recovery periods with a blood flow reduction of about 76%.</title><p>The responses were computed using the estimated optimal parameter values with the ischemia protocol: if (<italic>t</italic>&lt;0 min | <italic>t</italic>&gt;30 min) <italic>Q</italic> = 0.9 L/min, else <italic>Q</italic> = <italic>Q</italic>
<sub>isch</sub> = 0.216 L/min. The total O<sub>2</sub> and CO<sub>2</sub> contents are based on various forms of O<sub>2</sub> and CO<sub>2</sub> transports in the capillary blood and tissue cells (see ##SUPPL##0##Materials S1##).</p></caption></fig> <fig id="pone-0003168-g009" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.g009</object-id><label>Figure 9</label><caption><title>Model-predicted dynamic responses of metabolic reaction fluxes linked to the production and utilization of ACoA (A,B) and ATP (C,D) during the resting, ischemia and recovery periods with a blood flow reduction of about 76%.</title><p>The responses were computed using the estimated optimal parameter values with the ischemia protocol: if (<italic>t</italic>&lt;0 min | <italic>t</italic>&gt;30 min) <italic>Q</italic> = 0.9 L/min, else <italic>Q</italic> = <italic>Q</italic>
<sub>isch</sub> = 0.216 L/min. The factors 8 and 6 in plots B and D correspond to the stoichiometries of ACoA and ATP in the production reactions of ACoA and ATP (reactions 18 and 26), respectively.</p></caption></fig>
<table-wrap id="pone-0003168-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.t001</object-id><label>Table 1</label><caption><title>Average species concentrations in the muscle tissue cells ##REF##18390321##[1]## and their compartmentalization (distribution) into cytosol and mitochondria at normal, resting steady-state conditions.<xref ref-type="table-fn" rid="nt101">*</xref>
</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">
<italic>Species</italic>
</td><td colspan="2" align="left" rowspan="1">
<italic>Tissue Cells (100%)</italic>
</td><td colspan="2" align="left" rowspan="1">
<italic>Cytosol (90%)</italic>
</td><td colspan="2" align="left" rowspan="1">
<italic>Mitochondria (10%)</italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">
<italic>Concentration (mmol/kg cell ww)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>References</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Mass</italic>
<xref ref-type="table-fn" rid="nt102">#</xref>
<italic>(%)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Concentration (mmol/kg cyto ww)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Mass</italic>
<xref ref-type="table-fn" rid="nt102">#</xref>
<italic>(%)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Concentration (mmol/kg mito ww)</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GLC</td><td align="left" rowspan="1" colspan="1">0.5</td><td align="left" rowspan="1" colspan="1">
##UREF##3##[22]##, ##UREF##8##[33]##, ##REF##12456845##[38]##, ##REF##9688624##[39]##, ##REF##7713803##[46]##–##REF##3663177##[50]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.5556</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">GLY</td><td align="left" rowspan="1" colspan="1">95.0</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##12456845##[38]##, ##REF##7713803##[46]##, ##REF##10710508##[47]##, ##UREF##13##[49]##, ##REF##3512511##[51]##–##REF##2240197##[53]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">105.56</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">G6P</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">
##UREF##3##[22]##, ##UREF##8##[33]##, ##REF##7713803##[46]##–##REF##3663177##[50]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.2778</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">GA3P</td><td align="left" rowspan="1" colspan="1">0.08</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##7713803##[46]##, ##UREF##12##[48]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.08889</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">13BPG</td><td align="left" rowspan="1" colspan="1">0.08</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##7713803##[46]##, ##UREF##12##[48]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.08889</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">PYR</td><td align="left" rowspan="1" colspan="1">0.05</td><td align="left" rowspan="1" colspan="1">
##UREF##3##[22]##, ##REF##12456845##[38]##, ##REF##9688624##[39]##, ##REF##7713803##[46]##–##UREF##12##[48]##, ##REF##3663177##[50]##
</td><td align="left" rowspan="1" colspan="1">95</td><td align="left" rowspan="1" colspan="1">0.0528</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">0.025</td></tr><tr><td align="left" rowspan="1" colspan="1">LAC</td><td align="left" rowspan="1" colspan="1">0.78</td><td align="left" rowspan="1" colspan="1">
##UREF##3##[22]##, ##UREF##8##[33]##, ##REF##12456845##[38]##, ##REF##9688624##[39]##, ##REF##7713803##[46]##, ##UREF##12##[48]##, ##REF##3663177##[50]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.8667</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">ALA</td><td align="left" rowspan="1" colspan="1">1.3</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##2240197##[53]##–##REF##15121741##[55]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">1.4444</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">TGL</td><td align="left" rowspan="1" colspan="1">15.0</td><td align="left" rowspan="1" colspan="1">
##UREF##12##[48]##, ##REF##3512511##[51]##, ##REF##8271208##[52]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">16.667</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">GLR</td><td align="left" rowspan="1" colspan="1">0.065</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##10710508##[47]##, ##UREF##13##[49]##, ##REF##8271208##[52]##, ##REF##12231658##[56]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.07222</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">FFA</td><td align="left" rowspan="1" colspan="1">0.45</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##10710508##[47]##, ##UREF##13##[49]##, ##REF##8271208##[52]##, ##REF##12231658##[56]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.50</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">FAC</td><td align="left" rowspan="1" colspan="1">0.0035</td><td align="left" rowspan="1" colspan="1">
##UREF##13##[49]##, ##REF##5129319##[57]##
</td><td align="left" rowspan="1" colspan="1">95</td><td align="left" rowspan="1" colspan="1">0.003694</td><td align="left" rowspan="1" colspan="1">5</td><td align="left" rowspan="1" colspan="1">0.00175</td></tr><tr><td align="left" rowspan="1" colspan="1">ACoA</td><td align="left" rowspan="1" colspan="1">0.002</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##10710508##[47]##, ##REF##15121741##[55]##
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.02</td></tr><tr><td align="left" rowspan="1" colspan="1">CIT</td><td align="left" rowspan="1" colspan="1">0.095</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##9688624##[39]##, ##UREF##12##[48]##, ##REF##2240197##[53]##
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.95</td></tr><tr><td align="left" rowspan="1" colspan="1">AKG</td><td align="left" rowspan="1" colspan="1">0.0125</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##9688624##[39]##
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.125</td></tr><tr><td align="left" rowspan="1" colspan="1">SCoA</td><td align="left" rowspan="1" colspan="1">0.125</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##9688624##[39]##, ##UREF##12##[48]##, ##REF##2240197##[53]##
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">1.25</td></tr><tr><td align="left" rowspan="1" colspan="1">SUC</td><td align="left" rowspan="1" colspan="1">0.095</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##9688624##[39]##, ##UREF##12##[48]##, ##REF##2240197##[53]##
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.95</td></tr><tr><td align="left" rowspan="1" colspan="1">MAL</td><td align="left" rowspan="1" colspan="1">0.095</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##9688624##[39]##, ##UREF##12##[48]##, ##REF##2240197##[53]##
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.95</td></tr><tr><td align="left" rowspan="1" colspan="1">OXA</td><td align="left" rowspan="1" colspan="1">0.003</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##9688624##[39]##, ##UREF##12##[48]##, ##REF##2240197##[53]##
</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.03</td></tr><tr><td align="left" rowspan="1" colspan="1">CoA</td><td align="left" rowspan="1" colspan="1">0.02</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##10710508##[47]##, ##REF##15121741##[55]##
</td><td align="left" rowspan="1" colspan="1">80</td><td align="left" rowspan="1" colspan="1">0.01778</td><td align="left" rowspan="1" colspan="1">20</td><td align="left" rowspan="1" colspan="1">0.04</td></tr><tr><td align="left" rowspan="1" colspan="1">CO2 (F)</td><td align="left" rowspan="1" colspan="1">1.40 (46 mmHg)</td><td align="left" rowspan="1" colspan="1">
##REF##16775761##[26]##
</td><td align="left" rowspan="1" colspan="1">90</td><td align="left" rowspan="1" colspan="1">1.40 (46 mmHg)</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">1.40 (46 mmHg)</td></tr><tr><td align="left" rowspan="1" colspan="1">CO2 (T)</td><td align="left" rowspan="1" colspan="1">15.427</td><td align="left" rowspan="1" colspan="1">
##REF##16775761##[26]##
</td><td align="left" rowspan="1" colspan="1">90</td><td align="left" rowspan="1" colspan="1">15.427</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">15.427</td></tr><tr><td align="left" rowspan="1" colspan="1">O2 (F)</td><td align="left" rowspan="1" colspan="1">0.0338 (25 mmHg)</td><td align="left" rowspan="1" colspan="1">
##REF##16775761##[26]##
</td><td align="left" rowspan="1" colspan="1">90</td><td align="left" rowspan="1" colspan="1">0.0338 (25 mmHg)</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">0.0338 (25 mmHg)</td></tr><tr><td align="left" rowspan="1" colspan="1">O2 (T)</td><td align="left" rowspan="1" colspan="1">0.49</td><td align="left" rowspan="1" colspan="1">
##REF##16775761##[26]##
</td><td align="left" rowspan="1" colspan="1">90</td><td align="left" rowspan="1" colspan="1">0.49</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">0.49</td></tr><tr><td align="left" rowspan="1" colspan="1">PCR</td><td align="left" rowspan="1" colspan="1">21.0</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##12456845##[38]##, ##REF##6420545##[42]##, ##REF##7713803##[46]##, ##UREF##12##[48]##, ##REF##3663177##[50]##, ##REF##2240197##[53]##, ##UREF##14##[58]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">23.333</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">CR</td><td align="left" rowspan="1" colspan="1">10.5</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##12456845##[38]##, ##REF##6420545##[42]##, ##REF##7713803##[46]##, ##UREF##12##[48]##, ##REF##3663177##[50]##, ##REF##2240197##[53]##, ##UREF##14##[58]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">11.667</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">PI</td><td align="left" rowspan="1" colspan="1">2.75</td><td align="left" rowspan="1" colspan="1">
##REF##6420545##[42]##, ##REF##7713803##[46]##, ##UREF##12##[48]##
</td><td align="left" rowspan="1" colspan="1">90</td><td align="left" rowspan="1" colspan="1">2.75</td><td align="left" rowspan="1" colspan="1">10</td><td align="left" rowspan="1" colspan="1">2.75</td></tr><tr><td align="left" rowspan="1" colspan="1">NADH</td><td align="left" rowspan="1" colspan="1">0.05</td><td align="left" rowspan="1" colspan="1">
##REF##7713803##[46]##, ##UREF##12##[48]##, ##REF##3663177##[50]##, ##REF##2240197##[53]##
</td><td align="left" rowspan="1" colspan="1">0.5</td><td align="left" rowspan="1" colspan="1">0.2778E-3</td><td align="left" rowspan="1" colspan="1">99.5</td><td align="left" rowspan="1" colspan="1">0.4975</td></tr><tr><td align="left" rowspan="1" colspan="1">NAD<sup>+</sup>
</td><td align="left" rowspan="1" colspan="1">0.45</td><td align="left" rowspan="1" colspan="1">
##REF##7713803##[46]##, ##UREF##12##[48]##, ##REF##3663177##[50]##, ##REF##2240197##[53]##
</td><td align="left" rowspan="1" colspan="1">30</td><td align="left" rowspan="1" colspan="1">0.15</td><td align="left" rowspan="1" colspan="1">70</td><td align="left" rowspan="1" colspan="1">3.15</td></tr><tr><td align="left" rowspan="1" colspan="1">ATP (T)</td><td align="left" rowspan="1" colspan="1">6.2</td><td align="left" rowspan="1" colspan="1">
##UREF##3##[22]##, ##UREF##8##[33]##, ##REF##7713803##[46]##–##UREF##12##[48]##, ##REF##3663177##[50]##, ##REF##2240197##[53]##, ##UREF##14##[58]##
</td><td align="left" rowspan="1" colspan="1">86</td><td align="left" rowspan="1" colspan="1">5.924</td><td align="left" rowspan="1" colspan="1">14</td><td align="left" rowspan="1" colspan="1">8.68</td></tr><tr><td align="left" rowspan="1" colspan="1">ADP (T)</td><td align="left" rowspan="1" colspan="1">0.8</td><td align="left" rowspan="1" colspan="1">
##UREF##3##[22]##, ##REF##7713803##[46]##, ##UREF##12##[48]##, ##REF##3663177##[50]##, ##REF##2240197##[53]##, ##UREF##14##[58]##
</td><td align="left" rowspan="1" colspan="1">2</td><td align="left" rowspan="1" colspan="1">0.01778</td><td align="left" rowspan="1" colspan="1">98</td><td align="left" rowspan="1" colspan="1">7.84</td></tr><tr><td align="left" rowspan="1" colspan="1">AMP (T)</td><td align="left" rowspan="1" colspan="1">0.04</td><td align="left" rowspan="1" colspan="1">
##UREF##3##[22]##, ##REF##7713803##[46]##, ##UREF##12##[48]##, ##REF##3663177##[50]##, ##REF##2240197##[53]##, ##UREF##14##[58]##
</td><td align="left" rowspan="1" colspan="1">100</td><td align="left" rowspan="1" colspan="1">0.04444</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003168-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.t002</object-id><label>Table 2</label><caption><title>Average muscle tissue cells uptake-release rates (mmol/min) and blood species concentrations (mM) at normal, resting steady-state conditions from the literature.<xref ref-type="table-fn" rid="nt104">*</xref>
</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">
<italic>Species</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Arterial Concentrations, C<sub>art</sub></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Venous Concentrations, C<sub>ven</sub></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>UR rates, Q(C<sub>art</sub>−C<sub>ven</sub>)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>References</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GLC</td><td align="left" rowspan="1" colspan="1">5.0</td><td align="left" rowspan="1" colspan="1">4.7833</td><td align="left" rowspan="1" colspan="1">0.195</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##12456845##[38]##, ##REF##8271208##[52]##, ##REF##9831734##[54]##, ##REF##5129319##[57]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">PYR</td><td align="left" rowspan="1" colspan="1">0.075</td><td align="left" rowspan="1" colspan="1">0.0617</td><td align="left" rowspan="1" colspan="1">0.012</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##12456845##[38]##, ##REF##8271208##[52]##, ##REF##9831734##[54]##, ##REF##5129319##[57]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">LAC</td><td align="left" rowspan="1" colspan="1">0.5</td><td align="left" rowspan="1" colspan="1">0.6111</td><td align="left" rowspan="1" colspan="1">−0.10</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##12456845##[38]##, ##REF##8271208##[52]##, ##REF##9831734##[54]##, ##REF##5129319##[57]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">ALA</td><td align="left" rowspan="1" colspan="1">0.25</td><td align="left" rowspan="1" colspan="1">0.3333</td><td align="left" rowspan="1" colspan="1">−0.075</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##9688624##[39]##, ##REF##9831734##[54]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">GLR</td><td align="left" rowspan="1" colspan="1">0.04</td><td align="left" rowspan="1" colspan="1">0.0483</td><td align="left" rowspan="1" colspan="1">−0.0075</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##UREF##13##[49]##, ##REF##8271208##[52]##, ##REF##12231658##[56]##, ##REF##5129319##[57]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">FFA</td><td align="left" rowspan="1" colspan="1">0.7</td><td align="left" rowspan="1" colspan="1">0.6167</td><td align="left" rowspan="1" colspan="1">0.075</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##UREF##13##[49]##, ##REF##8271208##[52]##, ##REF##12231658##[56]##, ##REF##5129319##[57]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">CO<sub>2</sub> (T)</td><td align="left" rowspan="1" colspan="1">23.405</td><td align="left" rowspan="1" colspan="1">25.47</td><td align="left" rowspan="1" colspan="1">−1.86<xref ref-type="table-fn" rid="nt103">#</xref> (−41.66 ml/min)</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##8271208##[52]##, ##REF##9831734##[54]##, ##REF##5129319##[57]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">CO<sub>2</sub> (F)</td><td align="left" rowspan="1" colspan="1">1.22 (40 mmHg)</td><td align="left" rowspan="1" colspan="1">1.33 (43.55 mmHg)</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">
##REF##15682524##[25]##, ##REF##16775761##[26]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">O<sub>2</sub> (T)</td><td align="left" rowspan="1" colspan="1">9.235</td><td align="left" rowspan="1" colspan="1">6.556</td><td align="left" rowspan="1" colspan="1">2.41<xref ref-type="table-fn" rid="nt103">#</xref> (54.0 ml/min)</td><td align="left" rowspan="1" colspan="1">
##REF##12456845##[38]##, ##REF##8271208##[52]##, ##REF##9831734##[54]##, ##REF##5129319##[57]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">O<sub>2</sub> (F)</td><td align="left" rowspan="1" colspan="1">0.135 (100 mmHg)</td><td align="left" rowspan="1" colspan="1">0.0491 (36.37 mmHg)</td><td align="left" rowspan="1" colspan="1">NA</td><td align="left" rowspan="1" colspan="1">
##REF##15682524##[25]##, ##REF##16775761##[26]##
</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003168-t003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.t003</object-id><label>Table 3</label><caption><title>Optimal model parameter values for the inter-domain transport fluxes determined from the steady-state parameter estimation process.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td colspan="5" align="left" rowspan="1">I. Blood-cytosol passive or carrier-mediated (facilitated) transport fluxes (mmol/min) and parameters:
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Species</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>λ<sub>bl↔cyt</sub> (ml/min)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>T<sub>bl↔cyt</sub> (mmol/min)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>M<sub>bl↔cyt</sub> (mmol/L)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>J<sub>bl↔cyt</sub> (mmol/min)</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GLC</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">0.4106</td><td align="left" rowspan="1" colspan="1">2.5 ##REF##10413738##[27]##
</td><td align="left" rowspan="1" colspan="1">0.195</td></tr><tr><td align="left" rowspan="1" colspan="1">PYR</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">1.5127</td><td align="left" rowspan="1" colspan="1">1.0 ##REF##11820324##[28]##, ##REF##10358105##[29]##
</td><td align="left" rowspan="1" colspan="1">0.012</td></tr><tr><td align="left" rowspan="1" colspan="1">LAC</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">2.5769</td><td align="left" rowspan="1" colspan="1">5.0 ##REF##11820324##[28]##, ##REF##10358105##[29]##
</td><td align="left" rowspan="1" colspan="1">−0.10</td></tr><tr><td align="left" rowspan="1" colspan="1">ALA</td><td align="left" rowspan="1" colspan="1">0.0675</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">−0.075</td></tr><tr><td align="left" rowspan="1" colspan="1">GLR</td><td align="left" rowspan="1" colspan="1">0.3142</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">−0.0075</td></tr><tr><td align="left" rowspan="1" colspan="1">FFA</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">1.5589</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">0.075</td></tr><tr><td align="left" rowspan="1" colspan="1">CO2</td><td align="left" rowspan="1" colspan="1">24.926</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">−1.86</td></tr><tr><td align="left" rowspan="1" colspan="1">O2</td><td align="left" rowspan="1" colspan="1">157.1</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">N/A</td><td align="left" rowspan="1" colspan="1">2.41</td></tr></tbody></table><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td colspan="11" align="left" rowspan="1">II. Cytosol-mitochondrial species partition coefficients (unitless): <italic>σ<sub>cyt↔mit,j</sub></italic> = <italic>C<sub>mit,j</sub></italic>/<italic>C<sub>cyt,j</sub></italic>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Species</italic>
</td><td align="left" rowspan="1" colspan="1">PYR</td><td align="left" rowspan="1" colspan="1">FAC</td><td align="left" rowspan="1" colspan="1">CoA</td><td align="left" rowspan="1" colspan="1">PI</td><td align="left" rowspan="1" colspan="1">ATP</td><td align="left" rowspan="1" colspan="1">ADP</td><td align="left" rowspan="1" colspan="1">NADH</td><td align="left" rowspan="1" colspan="1">NAD<sup>+</sup>
</td><td align="left" rowspan="1" colspan="1">O2</td><td align="left" rowspan="1" colspan="1">CO2</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>σ<sub>cyt↔mit,j</sub></italic>
</td><td align="left" rowspan="1" colspan="1">0.4737</td><td align="left" rowspan="1" colspan="1">0.4737</td><td align="left" rowspan="1" colspan="1">2.25</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.465</td><td align="left" rowspan="1" colspan="1">441.0</td><td align="left" rowspan="1" colspan="1">1791.0</td><td align="left" rowspan="1" colspan="1">21.0</td><td align="left" rowspan="1" colspan="1">1.0</td><td align="left" rowspan="1" colspan="1">1.0</td></tr></tbody></table><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td colspan="4" align="left" rowspan="1">III. Anatomical volumes and muscle blood flow</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Parameter</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Definition</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Value and Unit</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Reference</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>V<sub>mus</sub></italic>
</td><td align="left" rowspan="1" colspan="1">Total volume of the active muscle</td><td align="left" rowspan="1" colspan="1">4.0 kg wet weight (L)</td><td align="left" rowspan="1" colspan="1">estimated</td></tr><tr><td align="left" rowspan="1" colspan="1">(<italic>V<sub>bl</sub></italic>,<italic>V<sub>isf</sub></italic>, <italic>V<sub>cyt</sub></italic>, <italic>V<sub>mit</sub></italic>)</td><td align="left" rowspan="1" colspan="1">Volumes of blood, ISF, cytosol and mitochondria</td><td align="left" rowspan="1" colspan="1">(0.07<italic>V<sub>m</sub></italic>, 0.13<italic>V<sub>m</sub></italic>, 0.72<italic>V<sub>m</sub></italic>, 0.08<italic>V<sub>m</sub></italic>)</td><td align="left" rowspan="1" colspan="1">calculated</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Q</italic>
</td><td align="left" rowspan="1" colspan="1">Muscle blood flow at rest for two legs</td><td align="left" rowspan="1" colspan="1">0.9 L/min</td><td align="left" rowspan="1" colspan="1">
##UREF##8##[33]##, ##REF##5129319##[57]##, ##REF##4815076##[59]##
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>Q</italic>
<sub>isch</sub>
</td><td align="left" rowspan="1" colspan="1">Muscle blood flow during ischemia</td><td align="left" rowspan="1" colspan="1">0.216 L/min</td><td align="left" rowspan="1" colspan="1">estimated</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003168-t004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.t004</object-id><label>Table 4</label><caption><title>Resting, steady-state metabolic reaction flux rates and corresponding optimal estimates of the  parameters that govern the metabolic reaction flux rates.<xref ref-type="table-fn" rid="nt107">*</xref>
</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">
<italic>Metabolic Reactions</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Flux rates (mmol/min)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>V<sub>max</sub> (mmol/min)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>K<sub>m</sub> (mmol/L)<sup>n</sup></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>K<sup>PS±</sup><sub>m</sub><italic>(unitless)</italic></italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>K<sup>RS±</sup><sub>m</sub><italic>(unitless)</italic></italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">1. Glucose Utilization</td><td align="left" rowspan="1" colspan="1">0.195</td><td align="left" rowspan="1" colspan="1">1.1171</td><td align="left" rowspan="1" colspan="1">6.9E-2, 0.224<sup>##</sup>
</td><td align="left" rowspan="1" colspan="1">334.13 (+)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">2. Glycogen Synthesis</td><td align="left" rowspan="1" colspan="1">0.25<sup>#</sup>
</td><td align="left" rowspan="1" colspan="1">0.8006</td><td align="left" rowspan="1" colspan="1">0.1641</td><td align="left" rowspan="1" colspan="1">337.69 (+)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">3. Glycogen Utilization</td><td align="left" rowspan="1" colspan="1">0.25<sup>#</sup>
</td><td align="left" rowspan="1" colspan="1">50.371</td><td align="left" rowspan="1" colspan="1">307.178</td><td align="left" rowspan="1" colspan="1">0.7473 (−)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">4. G6P Breakdown</td><td align="left" rowspan="1" colspan="1">0.195</td><td align="left" rowspan="1" colspan="1">39.615</td><td align="left" rowspan="1" colspan="1">27.778</td><td align="left" rowspan="1" colspan="1">337.055 (+)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">5. GA3P Breakdown</td><td align="left" rowspan="1" colspan="1">0.3825</td><td align="left" rowspan="1" colspan="1">43.331</td><td align="left" rowspan="1" colspan="1">14.354</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">540.233(−)</td></tr><tr><td align="left" rowspan="1" colspan="1">6. Pyruvate Production</td><td align="left" rowspan="1" colspan="1">0.3825</td><td align="left" rowspan="1" colspan="1">959.85</td><td align="left" rowspan="1" colspan="1">2.1199</td><td align="left" rowspan="1" colspan="1">0.3002 (−)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">7. Pyruvate Reduction</td><td align="left" rowspan="1" colspan="1">25<sup>#</sup>
</td><td align="left" rowspan="1" colspan="1">101.707</td><td align="left" rowspan="1" colspan="1">0.1575</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">3.904E-05(+)</td></tr><tr><td align="left" rowspan="1" colspan="1">8. Lactate Oxidation</td><td align="left" rowspan="1" colspan="1">24.9<sup>#</sup>
</td><td align="left" rowspan="1" colspan="1">3494.12</td><td align="left" rowspan="1" colspan="1">59.988</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">539.170 (−)</td></tr><tr><td align="left" rowspan="1" colspan="1">9. Alanine Production</td><td align="left" rowspan="1" colspan="1">0.075</td><td align="left" rowspan="1" colspan="1">7.5752</td><td align="left" rowspan="1" colspan="1">5.2779</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">10. TGL Synthesis</td><td align="left" rowspan="1" colspan="1">0.0075</td><td align="left" rowspan="1" colspan="1">5.7349</td><td align="left" rowspan="1" colspan="1">3.214E-2</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">1.25E-2 (+)</td></tr><tr><td align="left" rowspan="1" colspan="1">11. Lipolysis</td><td align="left" rowspan="1" colspan="1">0.0075</td><td align="left" rowspan="1" colspan="1">1.516E-2</td><td align="left" rowspan="1" colspan="1">17.029</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">12. FFA Activation</td><td align="left" rowspan="1" colspan="1">0.0975</td><td align="left" rowspan="1" colspan="1">19.705</td><td align="left" rowspan="1" colspan="1">0.8896</td><td align="left" rowspan="1" colspan="1">333.08 (+)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">13. ATP Hydrolysis</td><td align="left" rowspan="1" colspan="1">15.217</td><td align="left" rowspan="1" colspan="1">35.961</td><td align="left" rowspan="1" colspan="1">15.967</td><td align="left" rowspan="1" colspan="1">333.14 (+)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">14. PCR Breakdown</td><td align="left" rowspan="1" colspan="1">100<sup>#</sup>
</td><td align="left" rowspan="1" colspan="1">19822.6</td><td align="left" rowspan="1" colspan="1">22.091</td><td align="left" rowspan="1" colspan="1">0.3025 (−)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">15. PCR Synthesis</td><td align="left" rowspan="1" colspan="1">100<sup>#</sup>
</td><td align="left" rowspan="1" colspan="1">202.01</td><td align="left" rowspan="1" colspan="1">0.1142</td><td align="left" rowspan="1" colspan="1">333.425 (+)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">16. AMP Utilization</td><td align="left" rowspan="1" colspan="1">50<sup>#</sup>
</td><td align="left" rowspan="1" colspan="1">254.64</td><td align="left" rowspan="1" colspan="1">6.872E-02</td><td align="left" rowspan="1" colspan="1">333.297 (+)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">17. AMP Production</td><td align="left" rowspan="1" colspan="1">50<sup>#</sup>
</td><td align="left" rowspan="1" colspan="1">2064.73</td><td align="left" rowspan="1" colspan="1">1.214E-02</td><td align="left" rowspan="1" colspan="1">7.063E-2 (−)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">18. Pyruvate Oxidation</td><td align="left" rowspan="1" colspan="1">0.2195</td><td align="left" rowspan="1" colspan="1">0.5043</td><td align="left" rowspan="1" colspan="1">1.441E-04</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">6.383 (−)</td></tr><tr><td align="left" rowspan="1" colspan="1">19. FAC Oxidation</td><td align="left" rowspan="1" colspan="1">0.075</td><td align="left" rowspan="1" colspan="1">25.995</td><td align="left" rowspan="1" colspan="1">1.197E-2</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">6.430 (−)</td></tr><tr><td align="left" rowspan="1" colspan="1">20. Citrate Production</td><td align="left" rowspan="1" colspan="1">0.8195</td><td align="left" rowspan="1" colspan="1">83.932</td><td align="left" rowspan="1" colspan="1">6.085E-2</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">21. AKG Production</td><td align="left" rowspan="1" colspan="1">0.8195</td><td align="left" rowspan="1" colspan="1">3.7097</td><td align="left" rowspan="1" colspan="1">1.273</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">5.917 (−)</td></tr><tr><td align="left" rowspan="1" colspan="1">22. SCoA Production</td><td align="left" rowspan="1" colspan="1">0.8195</td><td align="left" rowspan="1" colspan="1">161.242</td><td align="left" rowspan="1" colspan="1">0.5007</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">5.987 (−)</td></tr><tr><td align="left" rowspan="1" colspan="1">23. Succinate Production</td><td align="left" rowspan="1" colspan="1">0.8195</td><td align="left" rowspan="1" colspan="1">3.1152</td><td align="left" rowspan="1" colspan="1">2.5616</td><td align="left" rowspan="1" colspan="1">0.4755 (−)</td><td align="left" rowspan="1" colspan="1">0</td></tr><tr><td align="left" rowspan="1" colspan="1">24. Malate Production</td><td align="left" rowspan="1" colspan="1">0.8195</td><td align="left" rowspan="1" colspan="1">1.6616</td><td align="left" rowspan="1" colspan="1">2.024E-02</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">6.238 (−)</td></tr><tr><td align="left" rowspan="1" colspan="1">25. Oxaloacetate Production</td><td align="left" rowspan="1" colspan="1">0.8195</td><td align="left" rowspan="1" colspan="1">1.9415</td><td align="left" rowspan="1" colspan="1">7.964E-02</td><td align="left" rowspan="1" colspan="1">0</td><td align="left" rowspan="1" colspan="1">7.509 (−)</td></tr><tr><td align="left" rowspan="1" colspan="1">26. Oxygen Utilization</td><td align="left" rowspan="1" colspan="1">2.41125</td><td align="left" rowspan="1" colspan="1">258.154</td><td align="left" rowspan="1" colspan="1">7.696E-04</td><td align="left" rowspan="1" colspan="1">5.1083 (−)</td><td align="left" rowspan="1" colspan="1">2.341(+)</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0003168-t005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0003168.t005</object-id><label>Table 5</label><caption><title>Model simulation results compared with the experimental data of Katz ##UREF##3##[22]## on muscle ischemia.<xref ref-type="table-fn" rid="nt108">*</xref>
</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1">
<italic>Species or Species Ratio</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Concentration Data (mM) (t = 0 min) (mean±SD)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Concentration Model (mM) (t = 0 min)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Concentration Data (mM) (t = 30 min) (mean±SD)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Concentration Model (mM) (t = 30 min)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Concentration Data (mM) (t = 45 min) (mean±SD)</italic>
</td><td align="left" rowspan="1" colspan="1">
<italic>Concentration Model (mM) (t = 45 min)</italic>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">GLC</td><td align="left" rowspan="1" colspan="1">0.485±0.0575</td><td align="left" rowspan="1" colspan="1">0.485</td><td align="left" rowspan="1" colspan="1">0.657±0.0925</td><td align="left" rowspan="1" colspan="1">0.676</td><td align="left" rowspan="1" colspan="1">0.653±0.0925</td><td align="left" rowspan="1" colspan="1">0.630</td></tr><tr><td align="left" rowspan="1" colspan="1">G6P</td><td align="left" rowspan="1" colspan="1">0.24±0.04</td><td align="left" rowspan="1" colspan="1">0.24</td><td align="left" rowspan="1" colspan="1">0.375±0.06</td><td align="left" rowspan="1" colspan="1">0.379</td><td align="left" rowspan="1" colspan="1">0.285±0.037</td><td align="left" rowspan="1" colspan="1">0.252</td></tr><tr><td align="left" rowspan="1" colspan="1">PYR</td><td align="left" rowspan="1" colspan="1">0.0425±0.006</td><td align="left" rowspan="1" colspan="1">0.0425</td><td align="left" rowspan="1" colspan="1">0.0925±0.0075</td><td align="left" rowspan="1" colspan="1">0.0918</td><td align="left" rowspan="1" colspan="1">0.0525±0.0075</td><td align="left" rowspan="1" colspan="1">0.057</td></tr><tr><td align="left" rowspan="1" colspan="1">LAC</td><td align="left" rowspan="1" colspan="1">0.688±0.12</td><td align="left" rowspan="1" colspan="1">0.688</td><td align="left" rowspan="1" colspan="1">1.535±0.38</td><td align="left" rowspan="1" colspan="1">1.545</td><td align="left" rowspan="1" colspan="1">0.94±0.12</td><td align="left" rowspan="1" colspan="1">0.850</td></tr><tr><td align="left" rowspan="1" colspan="1">LAC/PYR<xref ref-type="table-fn" rid="nt109">#</xref>
</td><td align="left" rowspan="1" colspan="1">16.2±2.9</td><td align="left" rowspan="1" colspan="1">16.2</td><td align="left" rowspan="1" colspan="1">16.6±4.1</td><td align="left" rowspan="1" colspan="1">16.8</td><td align="left" rowspan="1" colspan="1">17.9±2.2</td><td align="left" rowspan="1" colspan="1">14.9</td></tr><tr><td align="left" rowspan="1" colspan="1">PCR</td><td align="left" rowspan="1" colspan="1">20.1±0.65</td><td align="left" rowspan="1" colspan="1">20.1</td><td align="left" rowspan="1" colspan="1">16.55±1.175</td><td align="left" rowspan="1" colspan="1">16.55</td><td align="left" rowspan="1" colspan="1">20.225±0.6</td><td align="left" rowspan="1" colspan="1">19.94</td></tr><tr><td align="left" rowspan="1" colspan="1">CR</td><td align="left" rowspan="1" colspan="1">10.45±0.7</td><td align="left" rowspan="1" colspan="1">10.45</td><td align="left" rowspan="1" colspan="1">14.0±1.1</td><td align="left" rowspan="1" colspan="1">14.12</td><td align="left" rowspan="1" colspan="1">10.325±0.7</td><td align="left" rowspan="1" colspan="1">10.62</td></tr><tr><td align="left" rowspan="1" colspan="1">PCR/CR<xref ref-type="table-fn" rid="nt109">#</xref>
</td><td align="left" rowspan="1" colspan="1">1.97±0.12</td><td align="left" rowspan="1" colspan="1">1.97</td><td align="left" rowspan="1" colspan="1">1.35±0.13</td><td align="left" rowspan="1" colspan="1">1.17</td><td align="left" rowspan="1" colspan="1">2.03±0.15</td><td align="left" rowspan="1" colspan="1">1.88</td></tr><tr><td align="left" rowspan="1" colspan="1">Pi</td><td align="left" rowspan="1" colspan="1">2.70</td><td align="left" rowspan="1" colspan="1">2.7</td><td align="left" rowspan="1" colspan="1">6.25±0.8</td><td align="left" rowspan="1" colspan="1">6.2</td><td align="left" rowspan="1" colspan="1">2.6±0.3</td><td align="left" rowspan="1" colspan="1">2.845</td></tr><tr><td align="left" rowspan="1" colspan="1">ATP</td><td align="left" rowspan="1" colspan="1">6.15±0.175</td><td align="left" rowspan="1" colspan="1">6.15</td><td align="left" rowspan="1" colspan="1">6.225±0.15</td><td align="left" rowspan="1" colspan="1">6.092</td><td align="left" rowspan="1" colspan="1">6.2±0.175</td><td align="left" rowspan="1" colspan="1">6.148</td></tr><tr><td align="left" rowspan="1" colspan="1">ADP</td><td align="left" rowspan="1" colspan="1">0.82±0.025</td><td align="left" rowspan="1" colspan="1">0.82</td><td align="left" rowspan="1" colspan="1">0.885±0.04</td><td align="left" rowspan="1" colspan="1">0.869</td><td align="left" rowspan="1" colspan="1">0.855±0.035</td><td align="left" rowspan="1" colspan="1">0.822</td></tr><tr><td align="left" rowspan="1" colspan="1">AMP</td><td align="left" rowspan="1" colspan="1">0.04±0.005</td><td align="left" rowspan="1" colspan="1">0.04</td><td align="left" rowspan="1" colspan="1">0.05±0.007</td><td align="left" rowspan="1" colspan="1">0.0494</td><td align="left" rowspan="1" colspan="1">0.042±0.006</td><td align="left" rowspan="1" colspan="1">0.0403</td></tr><tr><td align="left" rowspan="1" colspan="1">TAN</td><td align="left" rowspan="1" colspan="1">7.0±0.175</td><td align="left" rowspan="1" colspan="1">7.01</td><td align="left" rowspan="1" colspan="1">7.05±0.125</td><td align="left" rowspan="1" colspan="1">7.01</td><td align="left" rowspan="1" colspan="1">7.1±0.175</td><td align="left" rowspan="1" colspan="1">7.01</td></tr></tbody></table></alternatives></table-wrap>
<disp-formula><label>(1)</label></disp-formula> <disp-formula><label>(2a)</label></disp-formula> <disp-formula><label>(2b)</label></disp-formula> <disp-formula><label>(2c)</label></disp-formula> <disp-formula><label>(2d)</label></disp-formula> <disp-formula><label>(3)</label></disp-formula> <disp-formula></disp-formula> <disp-formula><label>(4)</label></disp-formula> <disp-formula><label>(5)</label></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <disp-formula></disp-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula> <inline-formula></inline-formula>
<supplementary-material content-type="local-data" id="pone.0003168.s001"><label>Materials S1</label><caption><p>(0.34 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0003168.s002"><label>Materials S2</label><caption><p>(0.43 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><label>*</label><p>The species concentrations (mM or mmol/kg tissue cells ww) are converted from the experimental data (mmol/kg tissue cells dw) in the literature by multiplying a factor 0.25 kg tissue cells dw/kg tissue cells ww ##UREF##8##[33]##; for unit density, kg ww = L and mmol/kg ww = mmol/L = mM. The cytosolic and mitochondrial species concentrations are calculated from the species concentrations in the muscle tissue cells according to their approximate volumes of distributions and mass fractions; species nomenclature is adopted from Ref. ##REF##18390321##[1]## (also see the caption to ##FIG##0##Figure 1##); dw denotes dry weight, ww denotes wet weight.</p></fn><fn id="nt102"><label>#</label><p>The mass fractions are set to have reasonable compartmentalized species concentrations and mass action ratios consistent with available information from the literature: C<sub>cyt,LAC</sub>/C<sub>cyt,PYR</sub> = 16.4, C<sub>cyt,NAD+</sub>/C<sub>cyt,NADH</sub> = 540, C<sub>mit,NAD+</sub>/C<sub>mit,NADH</sub> = 6.3, C<sub>cyt,PCR</sub>/C<sub>cyt,CR</sub> = 2, C<sub>cyt,ATP</sub>/C<sub>cyt,ADP</sub> = 333.2, C<sub>mit,ATP</sub>/C<sub>mit,ADP</sub> = 1.11, <italic>K</italic>
<sub>LDH</sub> = (C<sub>cyt,LAC</sub>/C<sub>cyt,PYR</sub>)* (C<sub>cyt,NAD+</sub>/C<sub>cyt,NADH</sub>) = 8856, <italic>K</italic>
<sub>CK</sub> = (C<sub>cyt,CR</sub>/C<sub>cyt,PCR</sub>)*(C<sub>cyt,ATP</sub>/C<sub>cyt,ADP</sub>) = 166.6, and <italic>K</italic>
<sub>AK</sub> = (C<sub>cyt,ADP</sub>)<sup>2</sup>/(C<sub>cyt,ATP</sub>C<sub>cyt,AMP</sub>) = 1.2E-3. The cytosolic and mitochondrial PI concentrations are set at equal value with the assumption of a negligible pH gradient across the mitochondrial membrane.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt103"><label>#</label><p>
<italic>Physiological constrains</italic>:.</p></fn><fn id="nt104"><label>*</label><p>A resting blood flow of 0.9 L/min corresponding to the two-leg quadriceps femoris muscle is used to calculate the muscle tissue cells uptake-release (UR) rates from the arterio-venous (AV) differences. The venous species concentrations and uptake-release rates are tuned further to satisfy the physiological constrains<sup>#</sup> stated above, just below the table. (These well-compiled data are adapted from ##TAB##1##Table 2##, Ref. ##REF##18390321##[1]##.)</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt105"><label/><p>The table also includes the anatomical volumes and muscle blood flow at rest (fixed) and during muscle ischemia (estimated).*</p></fn><fn id="nt106"><label>*</label><p>The blood-cytosol transport fluxes satisfies <italic>J<sub>bl↔cyt,j</sub></italic> = <italic>UR<sub>j</sub></italic> = <italic>Q(C<sub>art,j</sub>−C<sub>ven,j</sub>)</italic> at resting, steady-state (the UR rates column in ##TAB##1##Table 2##). The subcellular cytosolic and mitochondrial domains are assumed to be in rapid equilibrium state so that the cytosol-mitochondria transport fluxes are considered negligible and the metabolites common to both of these domains are considered to have the similar dynamics (i.e., <italic>C<sub>mit,j</sub></italic>(<italic>t</italic>) = <italic>σ<sub>j</sub>.C<sub>cyt,j</sub></italic>(<italic>t</italic>)). Thus the cytosol-mitochondria metabolites partition coefficients <italic>σ<sub>j</sub></italic> are calculated based on the resting, steady-state metabolite concentrations in these two subcellular domains. The estimates of the transport flux parameters differ from those obtained in our previous work ##REF##18390321##[1]## because of the re-estimation of these parameters due to subcellular compartmentalization.</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt107"><label>*</label><p>The metabolic reaction flux rates are calculated based on resting, steady-state flux balance analysis (see ##TAB##2##Table 3##, Ref. ##REF##18390321##[1]##, for details). The flux rates with “#” can not be determined uniquely, and hence are set at values consistent with the other flux values (e.g., set a large value for the fast equilibrium reversible reactions like CK, AK and LDH). For <italic>K<sup>PS±</sup><sub>m</sub><italic>and</italic><italic><italic>K<sup>RS±</sup><sub>m</sub></italic>parameters, “+” indicates that the energy controller ratio is C<sub>ATP</sub>/C<sub>ADP</sub> or C<sub>ATP</sub>/C<sub>AMP</sub> or C<sub>NADH</sub>/C<sub>NAD+</sub>, “−” indicates that the energy controller ratio is C<sub>ADP</sub>/C<sub>ATP</sub> or C<sub>AMP</sub>/C<sub>ATP</sub> or C<sub>NAD+</sub>/C<sub>NADH</sub>, and “0” indicates that no controller ratio appears in the flux expression; ## denotes the estimated K<sub>m</sub> value  for the inhibition of hexokinase reaction by G6P.</italic></italic>
</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt108"><label>*</label><p>The model simulation are based on the ischemia protocol of Katz ##UREF##3##[22]##: if (<italic>t</italic>&lt;0 min | <italic>t</italic>&gt;30 min) <italic>Q</italic> = 0.9 L/min, else <italic>Q</italic> = <italic>Q</italic>
<sub>isch</sub> = 0.216 L/min. The species concentrations (mM or mmol/kg ww) are converted from the experimental data (mmol/kg dw) by multiplying a factor 0.25 kg dw/kg ww ##UREF##8##[33]##. For unit density, kg ww = L and mmol/kg ww = mmol/L = mM. The species concentrations are based on 100% muscle tissue cells, i.e., the volume averages of species concentrations in the cytosol and mitochondria, <italic>C<sub>cl</sub></italic> = (<italic>V<sub>cyt</sub>C<sub>cyt</sub></italic>+<italic>V<sub>mit</sub>C<sub>mit</sub></italic>)/<italic>V<sub>cl</sub></italic>. The species concentrations in the model are scaled with the resting (t = 0) species concentrations in the data (i.e., ) for comparison purpose; “dw” means dry weight and “ww” means wet weight.</p></fn><fn id="nt109"><label>#</label><p>Concentration ratios [LAC]/[PYR] and [PCR]/[CR] are unitless.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>The work presented in this manuscript was supported by grants from the National Aeronautics and Space Administration (NASA Johnson Space Center-NNJ06HD81G) and the National Institute of General Medical Sciences of the National Institute of Health (GM-66309). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p></fn></fn-group>
<graphic xlink:href="pone.0003168.g001"/> <graphic xlink:href="pone.0003168.g002"/> <graphic id="pone-0003168-t001-1" xlink:href="pone.0003168.t001"/> <graphic xlink:href="pone.0003168.e001.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003168.e002.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003168.e003.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003168.e004.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003168.e005.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003168.e006.jpg" mimetype="image" position="float"/> <graphic id="pone-0003168-t002-2" xlink:href="pone.0003168.t002"/> <graphic xlink:href="pone.0003168.e007"/> <graphic xlink:href="pone.0003168.g003"/> <graphic xlink:href="pone.0003168.e021.jpg" mimetype="image" position="float"/> <graphic xlink:href="pone.0003168.e008.jpg" mimetype="image" position="float"/> <inline-graphic xlink:href="pone.0003168.e009.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003168.e010.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003168.e011.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003168.e012.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003168.e013.jpg" mimetype="image"/> <graphic id="pone-0003168-t003-3" xlink:href="pone.0003168.t003"/> <graphic xlink:href="pone.0003168.e014"/> <inline-graphic xlink:href="pone.0003168.e015.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003168.e016.jpg" mimetype="image"/> <graphic id="pone-0003168-t004-4" xlink:href="pone.0003168.t004"/> <inline-graphic xlink:href="pone.0003168.e017.jpg" mimetype="image"/> <inline-graphic xlink:href="pone.0003168.e018.jpg" mimetype="image"/> <graphic xlink:href="pone.0003168.g004"/> <graphic xlink:href="pone.0003168.g005"/> <graphic xlink:href="pone.0003168.g006"/> <graphic id="pone-0003168-t005-5" xlink:href="pone.0003168.t005"/> <inline-graphic xlink:href="pone.0003168.e019.jpg" mimetype="image"/> <graphic xlink:href="pone.0003168.g007"/> <inline-graphic xlink:href="pone.0003168.e020.jpg" mimetype="image"/> <graphic xlink:href="pone.0003168.g008"/> <graphic xlink:href="pone.0003168.g009"/>
<media xlink:href="pone.0003168.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0003168.s002.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["12"], "element-citation": ["\n"], "surname": ["Lai", "Camesasca", "Saidel", "Dash", "Cabrera"], "given-names": ["N", "M", "GM", "RK", "ME"], "year": ["2007"], "article-title": ["Linking Pulmonary Oxygen Uptake, Muscle Oxygen Utilization and Cellular Metabolism during Exercise."], "source": ["Ann Biomed Eng"]}, {"label": ["18"], "element-citation": ["\n"], "surname": ["Beard"], "given-names": ["DA"], "year": ["2005"], "article-title": ["A biophysical model of the mitochondrial respiratory system and oxidative phosphorylation."], "source": ["PLoS Comput Biol"], "volume": ["1"], "fpage": ["252"], "lpage": ["264"]}, {"label": ["19"], "element-citation": ["\n"], "surname": ["Beard"], "given-names": ["DA"], "year": ["2006"], "article-title": ["Modeling of Oxygen Transport and Cellular Energetics Explains Observations on In Vivo Cardiac Energy Metabolism."], "source": ["PLoS Comput Biol"], "volume": ["2"], "fpage": ["1093"], "lpage": ["1106"]}, {"label": ["22"], "element-citation": ["\n"], "surname": ["Katz"], "given-names": ["A"], "year": ["1988"], "article-title": ["G-1,6-P2, glycolysis, and energy metabolism during circulatory occlusion in human skeletal muscle."], "source": ["Am J Physiol Cell Physiol"], "volume": ["255"], "fpage": ["C140"], "lpage": ["C144"]}, {"label": ["23"], "element-citation": ["\n"], "surname": ["Fell"], "given-names": ["D"], "year": ["1996"], "source": ["Understanding the Control of Metabolism"], "publisher-loc": ["London"], "publisher-name": ["Portland Press"]}, {"label": ["24"], "element-citation": ["\n"], "surname": ["Nelson", "Cox"], "given-names": ["D", "M"], "year": ["2000"], "source": ["Lehninger Principles of Biochemistry (third edition)"], "publisher-loc": ["New York"], "publisher-name": ["Worth Publishers"]}, {"label": ["31"], "element-citation": ["\n"], "surname": ["Segel"], "given-names": ["I"], "year": ["1993"], "source": ["Enzyme Kinetics: Behavior and Analysis of Rapid Equilibrium and Steady-State Enzyme Systems"], "publisher-loc": ["New York"], "publisher-name": ["Wiley-Interscience"]}, {"label": ["32"], "element-citation": ["\n"], "surname": ["Radhakrishnan", "Hindmarsh"], "given-names": ["K"], "year": ["1993"], "source": ["Description and Use of LSODE, the Livermore Solver for Ordinary Differential Equations"], "publisher-loc": ["Washington DC"], "publisher-name": ["NASA RP-1327, NASA"]}, {"label": ["33"], "element-citation": ["\n"], "surname": ["Putman", "Jones", "Hultman", "Hollidge-Horvat", "Bonen"], "given-names": ["CT", "NL", "E", "MG", "A"], "year": ["1998"], "article-title": ["Effects of short-term submaximal training in humans on muscle metabolism in exercise."], "source": ["Am J Physiol Endocrinol Metab"], "volume": ["275"], "fpage": ["E132"], "lpage": ["E139"]}, {"label": ["34"], "element-citation": ["\n"], "surname": ["Lasdon", "Warren", "Jain", "Ratner"], "given-names": ["L", "A", "A", "M"], "year": ["1978"], "article-title": ["Design and testing of a generalized reduced gradient code for nonlinear programming."], "source": ["ACM Trans Math Software"], "volume": ["4"], "fpage": ["34"], "lpage": ["50"]}, {"label": ["43"], "element-citation": ["\n"], "surname": ["Timmons", "Poucher", "Constantin-Teodosiu", "Worrall", "MacDonald"], "given-names": ["JA", "SM", "D", "V", "IA"], "year": ["1996"], "article-title": ["Metabolic responses of canine gracilis muscle during contraction with partial ischemia."], "source": ["Am J Physiol Endocrinol Metab"], "volume": ["270"], "fpage": ["E400"], "lpage": ["E406"]}, {"label": ["45"], "element-citation": ["\n"], "surname": ["Beard", "Wu", "Cabrera", "Dash"], "given-names": ["DA", "F", "ME", "RK"], "year": ["2008"], "article-title": ["Modeling of Cellular Metabolism and Microcirculatory Transport."], "source": ["Microcirculation"], "volume": ["1"]}, {"label": ["48"], "element-citation": ["\n"], "surname": ["Phillips", "Green", "Tarnopolsky", "Heigenhauser", "Grant"], "given-names": ["SM", "HJ", "MA", "GJ", "SM"], "year": ["1996"], "article-title": ["Progressive effect of endurance training on metabolic adaptations in working skeletal muscle."], "source": ["Am J Physiol Endocrinol Metab"], "volume": ["270"], "fpage": ["E265"], "lpage": ["E272"]}, {"label": ["49"], "element-citation": ["\n"], "surname": ["Turcotte", "Richter", "Kiens"], "given-names": ["LP", "EA", "B"], "year": ["1992"], "article-title": ["Increased plasma FFA uptake and oxidation during prolonged exercise in trained vs. untrained humans."], "source": ["Am J Physiol Endocrinol Metab"], "volume": ["262"], "fpage": ["E791"], "lpage": ["E799"]}, {"label": ["58"], "element-citation": ["\n"], "surname": ["Sahlin", "Poortmans"], "given-names": ["K", "JR"], "year": ["2004"], "article-title": ["High-Energy Phosphates and Muscle Energetics."], "source": ["Principles of Exercise Biochemistry (Third edition)"], "publisher-loc": ["Basel, Karger"], "publisher-name": ["Med Sport Sci"], "fpage": ["87"], "lpage": ["107"]}]
PLoS One. 2008 Sep 9; 3(9):e3168