Time course of hepatic gluconeogenesis during hindlimb suspension unloading.

The goal of this work was to determine the time-dependent changes in fractional hepatic gluconeogenesis (GNG) during conditions of hindlimb suspension unloading (HSU), a 'ground-based' method for inducing muscular atrophy to simulate space flight. We hypothesized that GNG would increase in HSU conditions as a result of metabolic shifts in the liver and skeletal muscle. A significant and progressive atrophy was observed in the soleus (30, 47 and 55%) and gastrocnemius muscles (0, 15 and 17%) after 3, 7 and 14 days of HSU, respectively. Fractional hepatic GNG was determined following the incorporation of deuterium from deuterated water ((2)H(2)O) into C-H bonds of newly synthesized glucose after an 8 h fast. Enrichment of plasma glucose with (2)H was measured using the classic method of Landau et al. (the 'hexamethylenetetramine (HMT) method'), based on specific (2)H labelling of glucose carbons, and the novel method of Chacko et al. ('average method'), based on the assumption of equal (2)H enrichment on all glucose carbons (except C2). After 3 days of HSU, fractional GNG was significantly elevated in the HSU group, as determined by either method (∼13%, P < 0.05). After 7 and 14 days of HSU, gluconeogenesis was not significantly different. Both analytical methods yielded similar time-dependent trends in gluconeogenic rates, but GNG values determined using the average method were consistently lower (∼30%) than those found by the HMT method. To compare and validate the average method against the HMT method further, we starved animals for 13 h to allow for hepatic GNG to contribute 100% to endogenous glucose production. The HMT method yielded 100% GNG, while the average method yielded GNG of ∼70%. As both methods used the same values of precursor enrichment, we postulated that the underestimation of gluconeogenic rate was as a result of differences in the measurements of product enrichment ((2)H labelling of plasma glucose). This could be explained by the following factors: (i) loss of deuterium via exchange between acetate and glucose; (ii) interference caused by fragment m/z 169, representing multiple isobaric species; and (iii) interference from other sugars at m/z 169. In conclusion, HSU caused a time-dependent increase in hepatic gluconeogenesis, irrespective of the analytical methods used.

Exercise intensity and oxygen uptake kinetics in African-American and Caucasian women.

The effect of exercise intensity on the on- and off-transient kinetics of oxygen uptake (VO(2)) was investigated in African American (AA) and Caucasian (C) women. African American (n = 7) and Caucasian (n = 6) women of similar age, body mass index and weight, performed an incremental test and bouts of square-wave exercise at moderate, heavy and very heavy intensities on a cycle ergometer. Gas exchange threshold (LT(GE)) was lower in AA (13.6 ± 2.3 mL kg(-1) min(-1)) than C (18.6 ± 5.6 mL kg(-1) min(-1)). The dynamic exercise and recovery VO(2) responses were characterized by mathematical models. There were no significant differences in (1) peak oxygen uptake (VO(2peak)) between AA (28.5 ± 5 mL kg(-1) min(-1)) and C (31.1 ± 6.6 mL kg(-1) min(-1)) and (2) VO(2) kinetics at any exercise intensity. At moderate exercise, the on- and off- VO(2) kinetics was described by a monoexponential function with similar time constants τ (1,on) (39.4 ± 12.5; 38.8 ± 15 s) and τ (1,off) (52.7 ± 10.1; 40.7 ± 4.4 s) for AA and C, respectively. At heavy and very heavy exercise, the VO(2) kinetics was described by a double-exponential function. The parameter values for heavy and very heavy exercise in the AA group were, respectively: τ (1,on) (47.0 ± 10.8; 44.3 ± 10 s), τ (2,on) (289 ± 63; 219 ± 90 s), τ (1,off) (45.9 ± 6.2; 50.7 ± 10 s), τ (2,off) (259 ± 120; 243 ± 93 s) while in the C group were, respectively: τ (1,on) (41 ± 12; 43.2 ± 15 s); τ (2, on) (277 ± 81; 215 ± 36 s), τ (1,off) (40.2 ± 3.4; 42.3 ± 7.2 s), τ (2,off) (215 ± 133; 228 ± 64 s). The on- and off-transients were symmetrical with respect to model order and dependent on exercise intensity regardless of race. Despite similar VO(2) kinetics, LT(GE) and gain of the VO(2) on-kinetics at moderate intensity were lower in AA than C. However, generalization to the African American and Caucasian populations is constrained by the small subject numbers.

Computational model of cellular metabolic dynamics: effect of insulin on glucose disposal in human skeletal muscle.

Identifying the mechanisms by which insulin regulates glucose metabolism in skeletal muscle is critical to understanding the etiology of insulin resistance and type 2 diabetes. Our knowledge of these mechanisms is limited by the difficulty of obtaining in vivo intracellular data. To quantitatively distinguish significant transport and metabolic mechanisms from limited experimental data, we developed a physiologically based, multiscale mathematical model of cellular metabolic dynamics in skeletal muscle. The model describes mass transport and metabolic processes including distinctive processes of the cytosol and mitochondria. The model simulated skeletal muscle metabolic responses to insulin corresponding to human hyperinsulinemic-euglycemic clamp studies. Insulin-mediated rate of glucose disposal was the primary model input. For model validation, simulations were compared with experimental data: intracellular metabolite concentrations and patterns of glucose disposal. Model variations were simulated to investigate three alternative mechanisms to explain insulin enhancements: Model 1 (M.1), simple mass action; M.2, insulin-mediated activation of key metabolic enzymes (i.e., hexokinase, glycogen synthase, pyruvate dehydrogenase); or M.3, parallel activation by a phenomenological insulin-mediated intracellular signal that modifies reaction rate coefficients. These simulations indicated that models M.1 and M.2 were not sufficient to explain the experimentally measured metabolic responses. However, by application of mechanism M.3, the model predicts metabolite concentration changes and glucose partitioning patterns consistent with experimental data. The reaction rate fluxes quantified by this detailed model of insulin/glucose metabolism provide information that can be used to evaluate the development of type 2 diabetes.

Modeling oxygenation in venous blood and skeletal muscle in response to exercise using near-infrared spectroscopy.

Noninvasive, continuous measurements in vivo are commonly used to make inferences about mechanisms controlling internal and external respiration during exercise. In particular, the dynamic response of muscle oxygenation (Sm(O(2))) measured by near-infrared spectroscopy (NIRS) is assumed to be correlated to that of venous oxygen saturation (Sv(O(2))) measured invasively. However, there are situations where the dynamics of Sm(O(2)) and Sv(O(2)) do not follow the same pattern. A quantitative analysis of venous and muscle oxygenation dynamics during exercise is necessary to explain the links between different patterns observed experimentally. For this purpose, a mathematical model of oxygen transport and utilization that accounts for the relative contribution of hemoglobin (Hb) and myoglobin (Mb) to the NIRS signal was developed. This model includes changes in microvascular composition within skeletal muscle during exercise and integrates experimental data in a consistent and mechanistic manner. Three subjects (age 25.6 +/- 0.6 yr) performed square-wave moderate exercise on a cycle ergometer under normoxic and hypoxic conditions while muscle oxygenation (C(oxy)) and deoxygenation (C(deoxy)) were measured by NIRS. Under normoxia, the oxygenated Hb/Mb concentration (C(oxy)) drops rapidly at the onset of exercise and then increases monotonically. Under hypoxia, C(oxy) decreases exponentially to a steady state within approximately 2 min. In contrast, model simulations of venous oxygen concentration show an exponential decrease under both conditions due to the imbalance between oxygen delivery and consumption at the onset of exercise. Also, model simulations that distinguish the dynamic responses of oxy-and deoxygenated Hb (HbO(2), HHb) and Mb (MbO(2), HMb) concentrations (C(oxy) = HbO(2) + MbO(2); C(deoxy) = HHb + HMb) show that Hb and Mb contributions to the NIRS signal are comparable. Analysis of NIRS signal components during exercise with a mechanistic model of oxygen transport and metabolism indicates that changes in oxygenated Hb and Mb are responsible for different patterns of Sm(O(2)) and Sv(O(2)) dynamics observed under normoxia and hypoxia.

Non-invasive estimation of metabolic flux and blood flow in working muscle: effect of blood-tissue distribution.

Muscle oxygenation measurements by near infrared spectroscopy (NIRS) are frequently obtained in humans to make inferences about mechanisms of metabolic control of respiration in working skeletal muscle. However, these measurements have technical limitations that can mislead the evaluation of tissue processes. In particular, NIRS measurements of working muscle represent oxygenation of a mix of fibers with heterogeneous activation, perfusion and architecture. Specifically, the relative volume distribution of capillaries, small arteries, and venules may affect NIRS data. To determine the effect of spatial volume distribution of components of working muscle on oxygen utilization dynamics and blood flow changes, a mathematical model of oxygen transport and utilization was developed. The model includes blood volume distribution within skeletal muscle and accounts for convective, diffusive, and reactive processes of oxygen transport and metabolism in working muscle. Inputs to the model are arterial O2 concentration, cardiac output and ATP demand. Model simulations were compared to exercise data from human subjects during a rest-to-work transition. Relationships between muscle oxygen consumption, blood flow, and the rate coefficient of capillary-tissue transport are analyzed. Blood volume distribution in muscle has noticeable effects on the optimal estimates of metabolic flux and blood flow in response to an exercise stimulus.

Modeling of cellular metabolism and microcirculatory transport.

Oxygen and other substrates, waste products, hormone messengers, and cells and other particles of the immune system are all transported in a closed-loop circulatory system in vertebrates, within which pumped blood travels to within diffusion distances of practically every cell in the body. Exchange of oxygen and carbon dioxide in the pulmonary capillaries and absorption of nutrients in the gut provide the circulating blood with biochemical reactants to sustain bioenergetic processes throughout the body. Inputs and outputs transported by the microcirculation are necessary to drive the open-system nonequilibrium chemical reactions of metabolism that are essential for cellular function. In turn, metabolically derived signals influence microcirculatory dynamics. Indeed, the microcirculation is the key system that ties processes at the whole-body level of the cardiovascular system to subcellular phenomena. This tight integration between cellular metabolism and microcirculatory transport begs for integrative simulations that span the cell, tissue, and organ scales.

Multi-scale model of O2 transport and metabolism: response to exercise.

Regulation of pulmonary oxygen uptake (VO2p) during exercise depends on cellular energy demand, blood flow, ventilation, oxygen exchange across membranes, and oxygen utilization in the contracting skeletal muscle. In human and animal studies of metabolic processes that control cellular respiration in working skeletal muscle, pulmonary VO2 dynamics is measured at the mouth using indirect calorimetry. To provide information on the dynamic balance between oxygen delivery and oxygen consumption at the microvascular level, muscle oxygenation is measured using near-infrared spectroscopy. A multi-scale computational model that links O2 transport and cellular metabolism in the skeletal muscle was developed to relate the measurements and gain quantitative understanding of the regulation of VO2 at the cellular, tissue, and whole-body level. The model incorporates mechanisms of oxygen transport from the airway openings to the cell, as well as the phosphagenic and oxidative pathways of ATP synthesis in the muscle cells.

The role of Ca2+ in coupling cardiac metabolism with regulation of contraction: in silico modeling.

The heart adapts the rate of mitochondrial ATP production to energy demand without noticeable changes in the concentration of ATP, ADP and Pi, even for large transitions between different workloads. We suggest that the changes in demand modulate the cytosolic Ca2+ concentration that changes mitochondrial Ca2+ to regulate ATP production. Thus, the rate of ATP production by the mitochondria is coupled to the rate of ATP consumption by the sarcomere cross-bridges (XBs). An integrated model was developed to couple cardiac metabolism and mitochondrial ATP production with the regulation of Ca2+ transient and ATP consumption by the sarcomere. The model includes two interrelated systems that run simultaneously utilizing two different integration steps: (1) The faster system describes the control of excitation contraction coupling with fast cytosolic Ca2+ transients, twitch mechanical contractions, and associated fluctuations in the mitochondrial Ca2+. (2) A slower system simulates the metabolic system, which consists of three different compartments: blood, cytosol, and mitochondria. The basic elements of the model are dynamic mass balances in the different compartments. Cytosolic Ca2+ handling is determined by four organelles: sarcolemmal Ca2+ influx and efflux; sarcoplasmic reticulum (SR) Ca2+ release and sequestration (SR); binding and dissociation from sarcomeric regulatory troponin complexes; and mitochondrial Ca2+ flows. Mitochondrial Ca2+ flows are determined by the Ca2+ uniporter and the mitochondrial Na+Ca2+ exchanger. The cytosolic Ca2+ determines the rate of ATP consumption by the sarcomere. Ca2+ binding to troponin regulates the rate of XBs recruitment and force development. The mitochondrial Ca2+ concentration determines the pyruvate dehydrogenase activity and the rate of ATP production by the F(1)-F(0) ATPase. The workload modulates the cytosolic Ca2+ concentration through feedback loops. The preload and afterload affect the number of strong XBs. The number of strong XBs determines the affinity of troponin for Ca2+, which alters the cytosolic Ca2+ transient. Model simulations quantify the role of Ca2+ in simultaneously controlling the power of contraction and the rate of ATP production. It explains the established empirical observation that significant changes in the metabolic fluxes can occur without significant changes in the key nucleotide (ATP and ADP) concentrations. Quantitative investigations of the mechanisms underlying the cardiac control of biochemical to mechanical energy conversion may lead to novel therapeutic modalities for the ischemic and failing myocardium.

Role of the malate-aspartate shuttle on the metabolic response to myocardial ischemia.

The malate-aspartate (M-A) shuttle provides an important mechanism to regulate glycolysis and lactate metabolism in the heart by transferring reducing equivalents from cytosol into mitochondria. However, experimental characterization of the M-A shuttle has been incomplete because of limitations in quantifying cytosolic and mitochondrial metabolites. In this study, we developed a multi-compartment model of cardiac metabolism with detailed presentation of the M-A shuttle to quantitatively predict non-observable fluxes and metabolite concentrations under normal and ischemic conditions in vivo. Model simulations predicted that the M-A shuttle is functionally localized to a subdomain that spans the mitochondrial and cytosolic spaces. With the onset of ischemia, the M-A shuttle flux rapidly decreased to a new steady state in proportion to the reduction in blood flow. Simulation results suggest that the reduced M-A shuttle flux during ischemia was not due to changes in shuttle-associated enzymes and transporters. However, there was a redistribution of shuttle-associated metabolites in both cytosol and mitochondria. Therefore, the dramatic acceleration in glycolysis and the switch to lactate production that occur immediately after the onset of ischemia is mediated by reduced M-A shuttle flux through metabolite redistribution of shuttle associated species across the mitochondrial membrane.

Modeling cellular metabolism and energetics in skeletal muscle: large-scale parameter estimation and sensitivity analysis.

Skeletal muscle plays a major role in the regulation of whole-body energy metabolism during physiological stresses such as ischemia, hypoxia, and exercise. Current experimental techniques provide relatively little in vivo data on dynamic responses of metabolite concentrations and metabolic fluxes in skeletal muscle to such physiological stimuli. As a complementary approach to experimental measurements and as a framework for quantitatively analyzing available in vivo data, a physiologically based model of skeletal muscle cellular metabolism and energetics is developed. This model, which incorporates key transport and reaction processes, is based on dynamic mass balances of 30 chemical species in capillary (blood) and tissue (cell) domains. The reaction fluxes in the cellular domain are expressed in terms of a generalized Michaelis?Menten equation involving energy controller ratios ATP/ADP and ATP/ADP and NADH/NAD+ . This formalism introduces a large number of unknown parameters ( approximately 90). Estimating these parameters from in vivo sparse data and evaluating dynamic sensitivities of the model outputs with respect to these parameters is a challenging problem. Parameter estimation is accomplished using an efficient, nonlinear, constraint-based, optimization algorithm that minimizes differences between available experimental data and corresponding model outputs by explicitly utilizing equality constraints on resting fluxes and concentrations. With the estimated parameter values, the model is able to simulate dynamic responses to reduced blood flow (ischemia) of key metabolite concentrations and metabolic fluxes, both measured and nonmeasured. A general parameter sensitivity analysis is carried out to determine and characterize the parameters having the most and least effects on the measured outputs.

Muscle oxygen uptake differs from consumption dynamics during transients in exercise.

Relating external to internal respiration during exercise requires quantitative modeling analysis for reliable inferences with respect to metabolic rate. Often, oxygen transport and metabolism based on steady-state mass balances (Fick principle) and passive diffusion between blood and tissue are applied to link pulmonary to cellular respiration. Indeed, when the work rate does not change rapidly, a quasi-steady-state analysis based on the Fick principle is sufficient to estimate the rate of O2 consumption in working muscle. During exercise when the work rate changes quickly, however, non-invasive in vivo measurements to estimate muscle O2 consumption are not sufficient to characterize cellular respiration of working muscle. To interpret transient changes of venous O2 concentration, blood flow, and O2 consumption in working muscle, a mathematical model of O2 transport and consumption based on dynamic mass balances is required. In this study, a comparison is made of the differences between simulations of O2 uptake and O2 consumption within working skeletal muscle based on a dynamic model and quasi-steady-state approximations. The conditions are specified under which the quasi-steady-state approximation becomes invalid.

Models of muscle contraction and energetics.

How does skeletal muscle manage to regulate the pathways of ATP synthesis during large-scale changes in work rate while maintaining metabolic homeostasis remains unknown. The classic model of metabolic regulation during muscle contraction states that accelerating ATP utilization leads to increasing concentrations of ADP and Pi, which serve as substrates for oxidative phosphorylation and thus accelerate ATP synthesis. An alternative model states that both the ATP demand and ATP supply pathways are simultaneously activated. Here, we review experimental and computational models of muscle contraction and energetics at various organizational levels and compare them with respect to their pros and cons in facilitating understanding of the regulation of energy metabolism during exercise in the intact organism.

Clinical stress testing in the pediatric age group: a statement from the American Heart Association Council on Cardiovascular Disease in the Young, Committee on Atherosclerosis, Hypertension, and Obesity in Youth.

This statement is an updated report of the American Heart Association's previous publications on exercise in children. In this statement, exercise laboratory requirements for environment, equipment, staffing, and procedures are presented. Indications and contraindications to stress testing are discussed, as are types of testing protocols and the use of pharmacological stress protocols. Current stress laboratory practices are reviewed on the basis of a survey of pediatric cardiology training programs.

Model of oxygen transport and metabolism predicts effect of hyperoxia on canine muscle oxygen uptake dynamics.

Previous studies have shown that increased oxygen delivery, via increased convection or arterial oxygen content, does not speed the dynamics of oxygen uptake, Vo(2m), in dog muscle electrically stimulated at a submaximal metabolic rate. However, the dynamics of transport and metabolic processes that occur within working muscle in situ is typically unavailable in this experimental setting. To investigate factors affecting Vo(2m) dynamics at contraction onset, we combined dynamic experimental data across working muscle with a mechanistic model of oxygen transport and metabolism in muscle. The model is based on dynamic mass balances for O(2), ATP, and PCr. Model equations account for changes in cellular ATPase, oxidative phosphorylation, and creatine kinase fluxes in skeletal muscle during exercise, and cellular respiration depends on [ADP] and [O(2)]. Model simulations were conducted at different levels of arterial oxygen content and blood flow to quantify the effects of convection and diffusion of oxygen on the regulation of cellular respiration during step transitions from rest to isometric contraction in dog gastrocnemius muscle. Simulations of arteriovenous O(2) differences and (.)Vo(2m) dynamics were successfully compared with experimental data (Grassi B, Gladden LB, Samaja M, Stary CM, Hogan MC. J Appl Physiol 85: 1394-1403, 1998; and Grassi B, Gladden LB, Stary CM, Wagner PD, Hogan MC. J Appl Physiol 85: 1404-1412, 1998), thus demonstrating the validity of the model, as well as its predictive capability. The main findings of this study are: 1) the estimated dynamic response of oxygen utilization at contraction onset in muscle is faster than that of oxygen uptake; and 2) hyperoxia does not accelerate the dynamics of diffusion and consequently muscle oxygen uptake at contraction onset due to the hyperoxia-induced increase in oxygen stores. These in silico derived results cannot be obtained from experimental observations alone.

Multi-organ system model of O2 and CO2 transport during isocapnic and poikilocapnic hypoxia.

A multi-organ systems model of O(2) and CO(2) transport is developed to analyze the control of ventilation and blood flow during hypoxia. Among the aspects of the control processes that this model addressed are possible mechanisms responsible for the second phase of the ventilatory hypoxic response to mild hypoxia, i.e., hypoxic ventilatory decline (HVD). Species mass transport processes are described by compartmental mass balances in brain, heart, skeletal muscle, and "other tissues" connected in parallel via the circulation. In pulmonary and systemic capillaries and in the vasculature connecting the systemic tissues, species transport processes are represented by a one-dimensional, convection-dispersion model. The effects of bicarbonate acid-base buffering, hemoglobin, and myoglobin on the transport processes are included. The model incorporates feedback control mechanisms through a cardiorespiratory control system in which peripheral and central chemoreceptors sense O(2) and CO(2) partial pressures. Model simulations of the ventilatory responses to isocapnic and poikilocapnic hypoxia show two phases with distinct dynamics. A fast phase is discernable immediately after switching from normoxic to hypoxic conditions, while a delayed slow phase (HVD) typically becomes manifested after several minutes. Model simulations allow quantitative evaluation of several proposed mechanisms to account for HVD. Under isocapnic hypoxia, simulations indicate that an increase in brain blood flow has no effect on HVD, but that HVD can be entirely described by central ventilatory depression (CVD). Under poikilocapnic hypoxia, the hypocapnia caused by hypoxic hyperventilation has no effect on HVD.

A computational model of skeletal muscle metabolism linking cellular adaptations induced by altered loading states to metabolic responses during exercise.

BACKGROUND: The alterations in skeletal muscle structure and function after prolonged periods of unloading are initiated by the chronic lack of mechanical stimulus of sufficient intensity, which is the result of a series of biochemical and metabolic interactions spanning from cellular to tissue/organ level. Reduced activation of skeletal muscle alters the gene expression of myosin heavy chain isoforms to meet the functional demands of reduced mechanical load, which results in muscle atrophy and reduced capacity to process fatty acids. In contrast, chronic loading results in the opposite pattern of adaptations. METHODS: To quantify interactions among cellular and skeletal muscle metabolic adaptations, and to predict metabolic responses to exercise after periods of altered loading states, we develop a computational model of skeletal muscle metabolism. The governing model equations - with parameters characterizing chronic loading/unloading states- were solved numerically to simulate metabolic responses to moderate intensity exercise (WR < or = 40% VO2 max). RESULTS: Model simulations showed that carbohydrate oxidation was 8.5% greater in chronically unloaded muscle compared with the loaded muscle (0.69 vs. 0.63 mmol/min), while fat oxidation was 7% higher in chronically loaded muscle (0.14 vs. 0.13 mmol/min), during exercise. Muscle oxygen uptake (VO2) and blood flow (Q) response times were 29% and 44% shorter in chronically loaded muscle (0.4 vs. 0.56 min for VO2 and 0.25 vs. 0.45 min for Q). CONCLUSION: The present model can be applied to test complex hypotheses during exercise involving the integration and control of metabolic processes at various organizational levels (cellular to tissue) in individuals who have undergone periods of chronic loading or unloading.

An efficient deconvolution algorithm for estimating oxygen consumption during muscle activities.

The reconstruction of an unknown input function from noisy measurements in a biological system is an ill-posed inverse problem. Any computational algorithm for its solution must use some kind of regularization technique to neutralize the disastrous effects of amplified noise components on the computed solution. In this paper, following a hierarchical Bayesian statistical inversion approach, we seek estimates for the input function and regularization parameter (hyperparameter) that maximize the posterior probability density function. We solve the maximization problem simultaneously for all unknowns, hyperparameter included, by a suitably chosen quasi-Newton method. The optimization approach is compared to the sampling-based Bayesian approach. We demonstrate the efficiency and robustness of the deconvolution algorithm by applying it to reconstructing the time courses of mitochondrial oxygen consumption during muscle state transitions (e.g., from resting state to contraction and recovery), from the simulated noisy output of oxygen concentration dynamics on the muscle surface. The model of oxygen transport and metabolism in skeletal muscle assumes an in vitro cylindrical structure of the muscle in which the oxygen from the surrounding oxygenated solution diffuses into the muscle and is then consumed by the muscle mitochondria. The algorithm can be applied to other deconvolution problems by suitably replacing the forward model of the system.

Role of cellular compartmentation in the metabolic response to stress: mechanistic insights from computational models.

The mechanisms controlling ATP generation in the transition from normal resting conditions to either high work states or ischemia are poorly understood. ATP generation depends upon compartmentation between the mitochondria and cytosol of metabolic pathways and key energy transfer species that cannot be easily assessed experimentally. We developed a multicompartment mathematical model of cardiac metabolism to simulate the metabolic responses to ischemia and increased workload. The model is based on mass balances, transport, and metabolic processes in cardiac tissue, and has three distinct compartments (blood, cytosol, and mitochondria). In addition to distinguishing between cytosol and mitochondria, the model includes a cytosolic subcompartment for glycolytic metabolic channeling. The model simulations predict the rapid activation of glycogenolysis and lactate production at the onset of ischemia, and support the concept of localization of glycolysis to a cytosolic subcompartment. In addition, simulations show that mitochondrial NADH/NAD(+) is primarily determined by oxygen consumption during ischemia, while cytosolic NADH/NAD(+) and lactate production are largely a function of glycolytic flux during the initial phase, and is controlled by mitochondrial NADH/NAD(+) and the malate-aspartate shuttle during the steady state. Finally, the model predicts that metabolic activation with an abrupt increase in workload requires parallel activation of ATP hydrolysis, glycolysis, mitochondrial dehydrogenases, the electron transport chain, and ADP phosphorylation. Taken together, these studies demonstrate the importance of metabolic compartmentation in the regulation of cardiac energetics in response to acute stress, and they highlight the usefulness of computational models in this line of investigation.

Gender differences in V˙O2 and HR kinetics at the onset of moderate and heavy exercise intensity in adolescents.

The majority of the studies on V˙O2 kinetics in pediatric populations investigated gender differences in prepubertal children during submaximal intensity exercise, but studies are lacking in adolescents. The purpose of this study was to test the hypothesis that gender differences exist in the V˙O2 and heart rate (HR) kinetic responses to moderate (M) and heavy (H) intensity exercise in adolescents. Twenty-one healthy African-American adolescents (9 males, 15.8 ± 1.1 year; 12 females, 15.7 ± 1 year) performed constant work load exercise on a cycle ergometer at M and H. The V˙O2 kinetics of the male group was previously analyzed (Lai et al., Appl. Physiol. Nutr. Metab. 33:107-117, 2008b). For both genders, V˙O2 and HR kinetics were described with a single exponential at M and a double exponential at H. The fundamental time constant (τ1) of V˙O2 was significantly higher in female than male at M (45 ± 7 vs. 36 ± 11 sec, P < 0.01) and H (41 ± 8 vs. 29 ± 9 sec, P < 0.01), respectively. The functional gain (G1) was not statistically different between gender at M and statistically higher in females than males at H: 9.7 ± 1.2 versus 10.9 ± 1.3 mL min-1 W-1, respectively. The amplitude of the slow component was not significantly different between genders. The HR kinetics were significantly (τ1, P < 0.01) slower in females than males at M (61 ± 16 sec vs. 45 ± 20 sec, P < 0.01) and H (42 ± 10 sec vs. 30 ± 8 sec, P = 0.03). The G1 of HR was higher in females than males at M: 0.53 ± 0.11 versus 0.98 ± 0.2 bpm W-1 and H: 0.40 ± 0.11 versus 0.73 ± 0.23 bpm W-1, respectively. Gender differences in the V˙O2 and HR kinetics suggest that oxygen delivery and utilization kinetics of female adolescents differ from those in male adolescents.

Regulation of cardiac energetics: role of redox state and cellular compartmentation during ischemia.

The heart is capable of altering its metabolic rate during exercise or ischemia. Under most state transitions, the heart maintains the concentration of adenosine triphosphate (ATP) at relatively constant values, in spite of large fluctuations in metabolic rate or in the delivery of fuels and oxygen. However, the mechanisms responsible for the regulation of cardiac energetics under conditions of increased demand or reduced supply are still under debate. To improve quantitative understanding of the regulation of glycolysis and oxidative phosphorylation under physiological and pathological conditions, it is essential to assess the dynamics of cytosolic and mitochondrial nicotinamide adenine dinucleotide (NAD(+)) and its reduced form (NADH) during stress (e.g., ischemia, exercise). However, at present there are no reliable methods to measure the dynamics of redox state in vivo in these subcellular compartments. In the present study, computer simulations with a mathematical model of myocardial energy metabolism are used to investigate the role of cytosolic and mitochondrial redox states in regulating cardiac energetics during reduced myocardial blood flow.