Metabolic control over the oxygen consumption flux in intact skeletal muscle: in silico studies

Integrated Computational Model of Lung Tissue Bioenergetics

Altered lung tissue bioenergetics plays a key role in the pathogenesis of lung diseases. A wealth of information exists regarding the bioenergetic processes in mitochondria isolated from rat lungs, cultured pulmonary endothelial cells, and intact rat lungs under physiological and pathophysiological conditions. However, the interdependence of those processes makes it difficult to quantify the impact of a change in a single or multiple process(es) on overall lung tissue bioenergetics. Integrated computational modeling provides a mechanistic and quantitative framework for the bioenergetic data at different levels of biological organization. The objective of this study was to develop and validate an integrated computational model of lung bioenergetics using existing experimental data from isolated perfused rat lungs. The model expands our recently developed computational model of the bioenergetics of mitochondria isolated from rat lungs by accounting for glucose uptake and phosphorylation, glycolysis, and the pentose phosphate pathway. For the mitochondrial region of the model, values of kinetic parameters were fixed at those estimated in our recent model of the bioenergetics of mitochondria isolated from rat lungs. For the cytosolic region of the model, intrinsic parameters such as apparent Michaelis constants were determined based on previously published enzyme kinetics data, whereas extrinsic parameters such as maximal reaction and transport velocities were estimated by fitting the model solution to published data from isolated rat lungs. The model was then validated by assessing its ability to predict existing experimental data not used for parameter estimation, including relationships between lung nucleotides content, lung lactate production rate, and lung energy charge under different experimental conditions. In addition, the model was used to gain novel insights on how lung tissue glycolytic rate is regulated by exogenous substrates such as glucose and lactate, and assess differences in the bioenergetics of mitochondria isolated from lung tissue and those of mitochondria in intact lungs. To the best of our knowledge, this is the first model of lung tissue bioenergetics. The model provides a mechanistic and quantitative framework for integrating available lung tissue bioenergetics data, and for testing novel hypotheses regarding the role of different cytosolic and mitochondrial processes in lung tissue bioenergetics.

Muscle V˙O2-power output nonlinearity in constant-power, step-incremental, and ramp-incremental exercise: magnitude and underlying mechanisms

A computer model of the skeletal muscle bioenergetic system was used to simulate time courses of muscle oxygen consumption (V˙O2), cytosolic metabolite (ADP, PCr, P_i, and ATP) concentrations, and pH during whole‐body constant‐power exercise (CPE) (6 min), step‐incremental exercise (SIE) (30 W/3 min), and slow (10 W/min), medium (30 W/min), and fast (50 W/min) ramp‐incremental exercise (RIE). Different ESA (each‐step activation) of oxidative phosphorylation (OXPHOS) intensity‐ATP usage activity relationships, representing different muscle fibers recruitment patterns, gave best agreement with experimental data for CPE, and for SIE and RIE. It was assumed that the muscle V˙O2‐power output (PO) nonlinearity is related to a time‐ and PO‐dependent increase in the additional ATP usage underlying the slow component of the V˙O2 on‐kinetics minus the increase in ATP supply by anaerobic glycolysis leading to a decrease in V˙O2. The muscle V˙O2‐PO relationship deviated upward (+) or downward (−) from linearity above critical power (CP), and the nonlinearity equaled +16% (CPE),+12% (SIE), +8% (slow RIE), +1% (moderate RIE), and −2% (fast RIE) at the end of exercise, in agreement with experimental data. During SIE and RIE, changes in PCr and P_i accelerated moderately above CP, while changes in ADP and pH accelerated significantly with time and PO. It is postulated that the intensity of the additional ATP usage minus ATP supply by anaerobic glycolysis determines the size of the muscle V˙O2‐PO nonlinearity. It is proposed that the extent of the additional ATP usage is proportional to the time integral of PO ‐ CP above CP.

Regulation of oxidative phosphorylation is different in electrically- and cortically-stimulated skeletal muscle

A computer model of the skeletal muscle bioenergetic system was used to study the regulation of oxidative phosphorylation (OXPHOS) in electrically-stimulated and cortically-stimulated skeletal muscle. Two types of the dependence of the intensity of each-step activation (ESA) of OXPHOS complexes on ATP usage activity were tested: power-type dependence and saturating-type dependence. The dependence of muscle oxygen consumption ([Formula: see text]), phosphocreatine (PCr), cytosolic ADP, ATP, inorganic phosphate (Pi), pH and τp (characteristic transition time) of the principal component of the muscle [Formula: see text] on-kinetics on the ATP usage activity was simulated for both types of the ESA intensity-ATP usage activity dependence. Computer simulations involving the power-type dependence predict system properties that agree well with experimental data for electrically-stimulated muscle. On the other hand, model predictions for the saturating-type dependence in the presence of the 'additional' ATP usage (postulated previously to underlie the slow component of the VO2 on-kinetics) reproduce well system properties encountered in human skeletal muscle during voluntary exercise. It is postulated that the difference between the regulation and kinetic properties of the system in electrically- and cortically-stimulated muscle is mostly due to the different muscle fibers recruitment pattern. In the former, all fiber types are recruited in parallel already at low power output (PO) values, while in the latter type I fibers (with higher ESA intensity) are stimulated at low PO values, while type II fibers (especially type II b and IIx fibers) with low ESA intensity are recruited predominantly at high PO values.

Contribution of proton leak to oxygen consumption in skeletal muscle during intense exercise is very low despite large contribution at rest

A computer model was used to simulate the dependence of protonmotive force (Δp), proton leak and phenomenological (involving proton leak) ATP/O₂ ratio on work intensity in skeletal muscle. Δp, NADH and proton leak decreased with work intensity. The contribution of proton leak to oxygen consumption (V˙O2) decreased from about 60% at rest to about 3 and 1% at moderate and heavy/severe exercise, respectively, while the ATP/O₂ ratio increased from 2.1 to 5.5 and 5.7. A two-fold increase in proton leak activity or its decrease to zero decreased/increased the ATP/O₂ ratio by only about 3 and 1% during moderate and heavy/severe exercise, respectively. The low contribution of proton leak to V˙O2 in intensively working skeletal muscle was mostly caused by a huge increase in ATP usage intensity during rest-to-work transition, while OXPHOS, and thus oxidative ATP supply and V˙O2 related to it, was mostly stimulated by high each-step activation (ESA) of OXPHOS complexes. The contribution of proton leak to V˙O2 and ATP/O₂ ratio in isolated mitochondria should not be directly extrapolated to working muscle, as mitochondria lack ESA, at least in the absence of Ca²⁺, and therefore V˙O2 cannot be elevated as much as in intact muscle.

Faster and stronger manifestation of mitochondrial diseases in skeletal muscle than in heart related to cytosolic inorganic phosphate (Pi) accumulation

A model of the cell bioenergetic system was used to compare the effect of oxidative phosphorylation (OXPHOS) deficiencies in a broad range of moderate ATP demand in skeletal muscle and heart. Computer simulations revealed that kinetic properties of the system are similar in both cases despite the much higher mitochondria content and "basic" OXPHOS activity in heart than in skeletal muscle, because of a much higher each-step activation (ESA) of OXPHOS in skeletal muscle than in heart. Large OXPHOS deficiencies lead in both tissues to a significant decrease in oxygen consumption (V̇o2) and phosphocreatine (PCr) and increase in cytosolic ADP, Pi, and H(+) The main difference between skeletal muscle and heart is a much higher cytosolic Pi concentration in healthy tissue and much higher cytosolic Pi accumulation (level) at low OXPHOS activities in the former, caused by a higher PCr level in healthy tissue (and higher total phosphate pool) and smaller Pi redistribution between cytosol and mitochondria at OXPHOS deficiency. This difference does not depend on ATP demand in a broad range. A much greater Pi increase and PCr decrease during rest-to-moderate work transition in skeletal muscle at OXPHOS deficiencies than at normal OXPHOS activity significantly slows down the V̇o2 on-kinetics. Because high cytosolic Pi concentrations cause fatigue in skeletal muscle and can compromise force generation in skeletal muscle and heart, this system property can contribute to the faster and stronger manifestation of mitochondrial diseases in skeletal muscle than in heart. Shortly, skeletal muscle with large OXPHOS deficiencies becomes fatigued already during low/moderate exercise.

Mechanisms of Attenuation of Pulmonary V'O2 Slow Component in Humans after Prolonged Endurance Training

In this study we have examined the effect of prolonged endurance training program on the pulmonary oxygen uptake (V'O2) kinetics during heavy-intensity cycling-exercise and its impact on maximal cycling and running performance. Twelve healthy, physically active men (mean±SD: age 22.33±1.44 years, V'O2peak 3198±458 mL ∙ min-1) performed an endurance training composed mainly of moderate-intensity cycling, lasting 20 weeks. Training resulted in a decrease (by ~5%, P = 0.027) in V'O2 during prior low-intensity exercise (20 W) and in shortening of τp of the V'O2 on-kinetics (30.1±5.9 s vs. 25.4±1.5 s, P = 0.007) during subsequent heavy-intensity cycling. This was accompanied by a decrease of the slow component of V'O2 on-kinetics by 49% (P = 0.001) and a decrease in the end-exercise V'O2 by ~5% (P = 0.005). An increase (P = 0.02) in the vascular endothelial growth factor receptor 2 mRNA level and a tendency (P = 0.06) to higher capillary-to-fiber ratio in the vastus lateralis muscle were found after training (n = 11). No significant effect of training on the V'O2peak was found (P = 0.12). However, the power output reached at the lactate threshold increased by 19% (P = 0.01). The power output obtained at the V'O2peak increased by 14% (P = 0.003) and the time of 1,500-m performance decreased by 5% (P = 0.001). Computer modeling of the skeletal muscle bioenergetic system suggests that the training-induced decrease in the slow component of V'O2 on-kinetics found in the present study is mainly caused by two factors: an intensification of the each-step activation (ESA) of oxidative phosphorylation (OXPHOS) complexes after training and decrease in the ''additional" ATP usage rising gradually during heavy-intensity exercise.

Each-step activation of oxidative phosphorylation is necessary to explain muscle metabolic kinetic responses to exercise and recovery in humans

KEY POINTS

The basic control mechanisms of oxidative phosphorylation (OXPHOS) and glycolysis during work transitions in human skeletal muscle are still a matter of debate. We used simulations of skeletal muscle bioenergetics to identify key system features that contribute to this debate, by comparing kinetic model outputs with experimental human data, including phosphocreatine, pH, pulmonary oxygen uptake and fluxes of ATP production by OXPHOS (vOX), anaerobic glycolysis and creatine kinase in moderate and severe intensity exercise transitions. We found that each-step activation of particular OXPHOS complexes, NADH supply and glycolysis, and strong (third-order) glycolytic inhibition by protons was required to reproduce observed phosphocreatine, pH and vOX kinetics during exercise. A slow decay of each-step activation during recovery, which was slowed further following severe exercise, was necessary to reproduce the experimental findings. Well-tested computer models offer new insight in the control of the human skeletal muscle bioenergetic system during physical exercise.

ABSTRACT

To better understand muscle bioenergetic regulation, a previously-developed model of the skeletal muscle cell bioenergetic system was used to simulate the influence of: (1) each-step activation (ESA) of NADH supply (including glycolysis) and oxidative phosphorylation (OXPHOS) complexes and (2) glycolytic inhibition by protons on the kinetics of ATP synthesis from OXPHOS, anaerobic glycolysis and creatine kinase. Simulations were fitted to previously published experimental data of ATP production fluxes and metabolite concentrations during moderate and severe intensity exercise transitions in bilateral knee extension in humans. Overall, the computer simulations agreed well with experimental results. Specifically, a large (>5-fold) direct activation of all OXPHOS complexes was required to simulate measured phosphocreatine and OXPHOS responses to both moderate and severe intensity exercise. In addition, slow decay of ESA was required to fit phosphocreatine recovery kinetics, and the time constant of ESA decay was slower following severe (180 s) than moderate (90 s) exercise. Additionally, a strong inhibition of (anaerobic) glycolysis by protons (glycolytic rate inversely proportional to the cube of proton concentration) provided the best fit to the experimental pH kinetics, and may contribute to the progressive increase in oxidative ATP supply during acidifying contractions. During severe-intensity exercise, an 'additional' ATP usage (a 27% increase at 8 min, above the initial ATP supply) was necessary to explain the observed V̇O2 slow component. Thus, parallel activation of ATP usage and ATP supply (ESA), and a strong inhibition of ATP supply by anaerobic glycolysis, were necessary to simulate the kinetics of muscle bioenergetics observed in humans.

Exercise: Kinetic considerations for gas exchange

The activities of daily living typically occur at metabolic rates below the maximum rate of aerobic energy production. Such activity is characteristic of the nonsteady state, where energy demands, and consequential physiological responses, are in constant flux. The dynamics of the integrated physiological processes during these activities determine the degree to which exercise can be supported through rates of O₂ utilization and CO₂ clearance appropriate for their demands and, as such, provide a physiological framework for the notion of exercise intensity. The rate at which O₂ exchange responds to meet the changing energy demands of exercise--its kinetics--is dependent on the ability of the pulmonary, circulatory, and muscle bioenergetic systems to respond appropriately. Slow response kinetics in pulmonary O₂ uptake predispose toward a greater necessity for substrate-level energy supply, processes that are limited in their capacity, challenge system homeostasis and hence contribute to exercise intolerance. This review provides a physiological systems perspective of pulmonary gas exchange kinetics: from an integrative view on the control of muscle oxygen consumption kinetics to the dissociation of cellular respiration from its pulmonary expression by the circulatory dynamics and the gas capacitance of the lungs, blood, and tissues. The intensity dependence of gas exchange kinetics is discussed in relation to constant, intermittent, and ramped work rate changes. The influence of heterogeneity in the kinetic matching of O₂ delivery to utilization is presented in reference to exercise tolerance in endurance-trained athletes, the elderly, and patients with chronic heart or lung disease.

Regulation of oxidative phosphorylation during work transitions results from its kinetic properties

The regulation of oxidative phosphorylation (OXPHOS) during work transitions in skeletal muscle and heart is still not well understood. Different computer models of this process have been developed that are characterized by various kinetic properties. In the present research-polemic theoretical study it is argued that models belonging to one group (Model A), which predict that among OXPHOS complexes complex III keeps almost all of the metabolic control over oxygen consumption (Vo2) and involve a strong complex III activation by inorganic phosphate (Pi), lead to the conclusion that an increase in Pi is the main mechanism responsible for OXPHOS activation (feedback-activation mechanism). Models belonging to another group (Model B), which were developed to take into account an approximately uniform distribution of metabolic control over Vo2 among particular OXPHOS complexes (complex I, complex III, complex IV, ATP synthase, ATP/ADP carrier, phosphate carrier) encountered in experimental studies in isolated mitochondria, predict that all OXPHOS complexes are directly activated in parallel with ATP usage and NADH supply by some external cytosolic factor/mechanism during rest-to-work or low-to-high work transitions in skeletal muscle and heart ("each-step-activation" mechanism). Model B demonstrates that different intensities of each-step activation can account for the very different (slopes of) phenomenological Vo2-ADP relationships observed in various skeletal muscles and heart. Thus they are able to explain the differences in the regulation of OXPHOS during work transitions between skeletal muscle (where moderate changes in ADP take place) and intact heart in vivo (where ADP is essentially constant).

Ischemia/Reperfusion-inducible protein modulates the function of organic cation transporter 1 and multidrug and toxin extrusion 1

The recently identified ischemia/reperfusion-inducible protein (IRIP) has been reported to negatively modulate the activities of several transporters in cell culture systems. The goal of this study is to determine whether IRIP regulates the activities of OCT1 and MATE1, and hence the disposition in vivo of their substrate metformin, a therapeutic drug for diabetes and other obesity-related syndromes. In the uptake studies in the human embryonic kidney 293 cells overexpressing IRIP with and without OCT1 or MATE1, IRIP overexpression was found to significantly inhibit the uptake of 1-methyl-4-phenylpyridinium mediated by OCT1 or MATE1. In contrast, knockdown of IRIP by small hairpin RNA (shRNA) increased the transporter activities in vitro. IRIP overexpression decreased the membrane localization of transporter proteins without any changes in transcript levels in cells. By overexpressing IRIP in mouse liver via hydrodynamic tail vein injection, we demonstrated that increased IRIP expression could cause a significant reduction in hepatic accumulation of metformin (P < 0.01). In addition, we observed that the expression of IRIP was approximately half (P < 0.01) in ob/ob mice when compared to their lean littermates, with significant increases in hepatic Oct1 protein expression and metformin accumulation. In conclusion, IRIP negatively modulates the function of OCT1 and MATE1 in cells. Importantly, we provide in vivo evidence for such modulation that may cause an alteration in drug disposition. The regulation by IRIP on transporter activities likely occurs at a post-transcriptional level, and future studies are needed to characterize the exact mechanism.

Model analysis of the relationship between intracellular PO2 and energy demand in skeletal muscle

On the basis of experimental studies, the intracellular O(2) (iPo(2))-work rate (WR) relationship in skeletal muscle is not unique. One study found that iPo(2) reached a plateau at 60% of maximal WR, while another found that iPo(2) decreased linearly at higher WR, inferring capillary permeability-surface area (PS) and blood-tissue O(2) gradient, respectively, as alternative dominant factors for determining O(2) diffusion changes during exercise. This relationship is affected by several factors, including O(2) delivery and oxidative and glycolytic capacities of the muscle. In this study, these factors are examined using a mechanistic, mathematical model to analyze experimental data from contracting skeletal muscle and predict the effects of muscle contraction on O(2) transport, glycogenolysis, and iPo(2). The model describes convection, O(2) diffusion, and cellular metabolism, including anaerobic glycogenolysis. Consequently, the model simulates iPo(2) in response to muscle contraction under a variety of experimental conditions. The model was validated by comparison of simulations of O(2) uptake with corresponding experimental responses of electrically stimulated canine muscle under different O(2) content, blood flow, and contraction intensities. The model allows hypothetical variation of PS, glycogenolytic capacity, and blood flow and predictions of the distinctive effects of these factors on the iPo(2)-contraction intensity relationship in canine muscle. Although PS is the main factor regulating O(2) diffusion rate, model simulations indicate that PS and O(2) gradient have essential roles, depending on the specific conditions. Furthermore, the model predicts that different convection and diffusion patterns and metabolic factors may be responsible for different iPo(2)-WR relationships in humans.

Metabolic control and regulation of the tricarboxylic acid cycle in photosynthetic and heterotrophic plant tissues

The tricarboxylic acid (TCA) cycle is a crucial component of respiratory metabolism in both photosynthetic and heterotrophic plant organs. All of the major genes of the tomato TCA cycle have been cloned recently, allowing the generation of a suite of transgenic plants in which the majority of the enzymes in the pathway are progressively decreased. Investigations of these plants have provided an almost complete view of the distribution of control in this important pathway. Our studies suggest that citrate synthase, aconitase, isocitrate dehydrogenase, succinyl CoA ligase, succinate dehydrogenase, fumarase and malate dehydrogenase have control coefficients flux for respiration of -0.4, 0.964, -0.123, 0.0008, 0.289, 0.601 and 1.76, respectively; while 2-oxoglutarate dehydrogenase is estimated to have a control coefficient of 0.786 in potato tubers. These results thus indicate that the control of this pathway is distributed among malate dehydrogenase, aconitase, fumarase, succinate dehydrogenase and 2-oxoglutarate dehydrogenase. The unusual distribution of control estimated here is consistent with specific non-cyclic flux mode and cytosolic bypasses that operate in illuminated leaves. These observations are discussed in the context of known regulatory properties of the enzymes and some illustrative examples of how the pathway responds to environmental change are given.

The control systems structures of energy metabolism

The biochemical regulation of energy metabolism (EM) allows cells to modulate their energetic output depending on available substrates and requirements. To this end, numerous biomolecular mechanisms exist that allow the sensing of the energetic state and corresponding adjustment of enzymatic reaction rates. This regulation is known to induce dynamic systems properties such as oscillations or perfect adaptation. Although the various mechanisms of energy regulation have been studied in detail from many angles at the experimental and theoretical levels, no framework is available for the systematic analysis of EM from a control systems perspective. In this study, we have used principles well known in control to clarify the basic system features that govern EM. The major result is a subdivision of the biomolecular mechanisms of energy regulation in terms of widely used engineering control mechanisms: proportional, integral, derivative control, and structures: feedback, cascade and feed-forward control. Evidence for each mechanism and structure is demonstrated and the implications for systems properties are shown through simulations. As the equivalence between biological systems and control components presented here is generic, it is also hypothesized that our work could eventually have an applicability that is much wider than the focus of the current study.

Ischemia-reperfusion-inducible protein modulates cell sensitivity to anticancer drugs by regulating activity of efflux transporter

Human ischemia-reperfusion-inducible protein (hIRIP) or hYrdC belongs to the SUA5/YrdC/YciO protein family and affects activity of a variety of cellular transporters. We observed that overexpression of wild-type or dominant-negative mutant of hIRIP protein affects the cellular sensitivity to anticancer drugs with different mechanisms of toxicity. Here we investigated in detail the effect of hIRIP on cell sensitivity to doxorubicin and show that hIRIP inhibits the drug efflux. Multidrug-resistant P-glycoprotein was identified as one of the target transporters. IRIP does not influence P-glycoprotein biosynthesis but affects its processing and promotes degradation. We also show that P-glycoprotein is associated with COP-α, one of the proteins of the COPI complex. This interaction is sensitive to the level of hIRIP expression. These findings suggest that hIRIP expression can regulate cargo assembly and function of efflux transporters, including P-glycoprotein, which mediates one of the most common mechanisms of the multidrug resistance.

In vivo modular control analysis of energy metabolism in contracting skeletal muscle

We used (31)P MRS (magnetic resonance spectroscopy) measurements of energetic intermediates [ATP, P(i) and PCr (phosphocreatine)] in combination with the analytical tools of metabolic control analysis to study in vivo energy metabolism in the contracting skeletal muscle of anaesthetized rats over a broad range of workload. According to our recent MoCA (modular control analysis) used to describe regulatory mechanisms in beating heart, we defined the energetic system of muscle contraction as two modules (PCr-Producer and PCr-Consumer) connected by the energetic intermediates. Hypoxia and electrical stimulation were used in this in vivo study as reasonably selective modulations of Producer and Consumer respectively. As quantified by elasticity coefficients, the sensitivities of each module to PCr determine the control of steady-state contractile activity and metabolite concentrations. The magnitude of the elasticity of the producer was high (4.3+/-0.6) at low workloads and decreased 5-fold (to 0.9+/-0.2) at high workloads. By contrast, the elasticity of the consumer remained low (0.5-1.2) over the range of metabolic rates studied. The control exerted by each module over contraction was calculated from these elasticities. The control of contraction was found on the consumer at low workloads and then swung to the producer, due to the workload-dependent decrease in the elasticity of producer. The workload-dependent elasticity and control pattern of energy production in muscle is a major difference from heart. Since module rate and elasticity depend on the concentrations of substrates and products, the absence of homoeostasis of the energetic intermediates in muscle, by contrast with heart, is probably the origin of the workload-dependent elasticity of the producer module.

Role of NADH/NAD+ transport activity and glycogen store on skeletal muscle energy metabolism during exercise: in silico studies

Skeletal muscle can maintain ATP concentration constant during the transition from rest to exercise, whereas metabolic reaction rates may increase substantially. Among the key regulatory factors of skeletal muscle energy metabolism during exercise, the dynamics of cytosolic and mitochondrial NADH and NAD+ have not been characterized. To quantify these regulatory factors, we have developed a physiologically based computational model of skeletal muscle energy metabolism. This model integrates transport and reaction fluxes in distinct capillary, cytosolic, and mitochondrial domains and investigates the roles of mitochondrial NADH/NAD+ transport (shuttling) activity and muscle glycogen concentration (stores) during moderate intensity exercise (60% maximal O2 consumption). The underlying hypothesis is that the cytosolic redox state (NADH/NAD+) is much more sensitive to a metabolic disturbance in contracting skeletal muscle than the mitochondrial redox state. This hypothesis was tested by simulating the dynamic metabolic responses of skeletal muscle to exercise while altering the transport rate of reducing equivalents (NADH and NAD+) between cytosol and mitochondria and muscle glycogen stores. Simulations with optimal parameter estimates showed good agreement with the available experimental data from muscle biopsies in human subjects. Compared with these simulations, a 20% increase (or approximately 20% decrease) in mitochondrial NADH/NAD+ shuttling activity led to an approximately 70% decrease (or approximately 3-fold increase) in cytosolic redox state and an approximately 35% decrease (or approximately 25% increase) in muscle lactate level. Doubling (or halving) muscle glycogen concentration resulted in an approximately 50% increase (or approximately 35% decrease) in cytosolic redox state and an approximately 30% increase (or approximately 25% decrease) in muscle lactate concentration. In both cases, changes in mitochondrial redox state were minimal. In conclusion, the model simulations of exercise response are consistent with the hypothesis that mitochondrial NADH/NAD+ shuttling activity and muscle glycogen stores affect primarily the cytosolic redox state. Furthermore, muscle lactate production is regulated primarily by the cytosolic redox state.

Identification of cysteines in rat organic cation transporters rOCT1 (C322, C451) and rOCT2 (C451) critical for transport activity and substrate affinity

Effects of the sulfhydryl reagent methylmethanethiosulfonate (MMTS) on functions of organic cation transporters (OCTs) were investigated. Currents induced by 10 mM choline [ Imax(choline)] in Xenopus laevis oocytes expressing rat OCT1 (rOCT1) were increased four- to ninefold after 30-s incubation with 5 mM MMTS whereas Imax(choline) by rat OCT2 was 70% decreased. MMTS activated the rOCT1 transporter within the plasma membrane without changing stoichiometry between translocated charge and cation. After modification of oocytes expressing rOCT1 or rOCT2 with MMTS, I0.5(choline) values for choline-induced currents were increased. For rOCT1 it was shown that MMTS increased I0.5 values for different cations by different degrees. Mutagenesis of individual cysteine residues in rOCT1 revealed that modification of cysteine 322 in the large intracellular loop, and of cysteine 451 at the transition of the transmembrane α-helix (TMH) 10 to the short intracellular loop between the TMH 10 and 11 is responsible for the observed effects of MMTS. After replacement of cysteine 451 by methionine, the IC50(choline) for choline to inhibit MPP uptake by rOCT1 was increased whereas the I0.5(choline) value for choline-induced current remained unchanged. At variance, in double mutant Cys322Ser, Cys451Met, I0.5(choline) was increased compared with rOCT1 wild-type whereas in the single mutant Cys322Ser I0.5(choline) was not changed. The data suggest that modification of rOCT1 at cysteines 322 and 451 leads to an increase in turnover. They indicate that cysteine 451 in rOCT1 interacts with the large intracellular loop and that cysteine 451 in both rOCT1 and rOCT2 is critical for the affinity of choline.

Tripeptides of RS1 (RSC1A1) inhibit a monosaccharide-dependent exocytotic pathway of Na+-D-glucose cotransporter SGLT1 with high affinity

The human gene RSC1A1 codes for a 67-kDa protein named RS1 that mediates transcriptional and post-transcriptional regulation of Na(+)-D-glucose cotransporter SGLT1. The post-transcriptional regulation occurs at the trans-Golgi network (TGN). We identified two tripeptides in human RS1 (Gln-Cys-Pro (QCP) and Gln-Ser-Pro (QSP)) that induce posttranscriptional down-regulation of SGLT1 at the TGN leading to 40-50% reduction of SGLT1 in plasma membrane. For effective intracellular concentrations IC(50) values of 2.0 nM (QCP) and 0.16 nm (QSP) were estimated. Down-regulation of SGLT1 by tripeptides was attenuated by intracellular monosaccharides including non-metabolized methyl-alpha-D-glucopyranoside and 2-deoxyglucose. In small intestine post-transcriptional regulation of SGLT1 may contribute to glucose-dependent regulation of liver metabolism and intestinal mobility. QCP and QSP are transported by the H(+)-peptide cotransporter PepT1 that is colocated with SGLT1 in small intestinal enterocytes. Using coexpression of SGLT1 and PepT1 in Xenopus oocytes or polarized Caco-2 cells that contain both transporters we demonstrated that the tripeptides were effective when applied to the extracellular compartment. After a 1-h perfusion of intact rat small intestine with QSP, glucose absorption was reduced by 30%. The data indicate that orally applied tripeptides can be used to down-regulate small intestinal glucose absorption, e.g. in diabetes mellitus.

Interleukin-6 stimulates alpha-MG uptake in renal proximal tubule cells: involvement of STAT3, PI3K/Akt, MAPKs, and NF-kappaB

Recent studies have shown that interleukin 6 (IL-6) acts on the cellular proliferation-activating transduction signals during cellular regeneration. Therefore, this study examined the effect of IL-6 on the activation of Na(+)/glucose cotransporters (SGLTs) and its related signaling pathways in primary cultured renal proximal tubule cells (PTCs). IL-6 increased the level of alpha-methyl-d-[(14)C]glucopyranoside (alpha-MG) uptake in time- and dose-dependent manners. IL-6 also increased SGLT1 plus SGLT2 mRNA and protein expression level. The IL-6 receptors (IL-6Ralpha and gp 130) were expressed in PTCs. In addition, genistein and herbimycin A completely blocked the IL-6-induced increases in alpha-MG uptake and the protein expression level of SGLTs. On the other hand, IL-6 increased the level of 5-(and-6)-chloromethyl-2',7'-dichlorodihydrofluorescein diacetate-sensitive cellular reactive oxygen species (ROS), and IL-6-induced increases in alpha-MG uptake and the protein expression level of SGLTs were blocked by ascorbic acid or taurine (antioxidants). IL-6 also increased the phosphorylation of signal transducer and activator of transcription-3 (STAT3), phosphoinositide-3 kinase (PI3K)/Akt, and mitogen-activated protein kinases (MAPKs) in a time-dependent manner. A pretreatment with STAT3 inhibitor LY 294002, an Akt inhibitor, or MAPK inhibitors significantly blocked the IL-6-induced increase in alpha-MG uptake. In addition, IL-6 increased the level of nuclear factor-kappaB (NF-kappaB) phosphorylation. A pretreatment with SN50 or BAY 11-7082 also blocked the IL-6-induced increase in alpha-MG uptake. In conclusion, IL-6 increases the SGLT activity through ROS, and its action in renal PTCs is associated with the STAT3, PI3K/Akt, MAPKs, and NF-kappaB signaling pathways.

Intestinal d-Galactose Transport in an Endotoxemia Model in the Rabbit

Lipopolysaccharide (LPS) is an endotoxin causing sepsis. Studies from our laboratory revealed impaired intestinal absorption of L-leucine and D-fructose in LPS-treated rabbits. The aim of this study was to examine intestinal D-galactose transport following intravenous administration of LPS in the rabbit and to identify the cellular mechanisms driving this process. Endotoxin treatment diminished the buildup of D-galactose in intestinal tissue, the mucosal to serosal transepithelial flux of the sugar and its uptake by brush border membrane vesicles (BBMVs). Intracellular signaling pathways associated with protein kinase C (PKC), protein kinase A (PKA), p38 mitogen-activated protein kinase (p38MAPK), Jun N-terminal kinase (JNK), MAPK/extracellular signal-regulated kinases 1 and 2 (MEK1/2) and proteasome were found to be involved in this reduction in sugar uptake. Na(+)/glucose cotransporter 1 (SGLT1) protein levels in BBMVs were lower for LPS-treated animals than control animals. These findings indicate that LPS inhibits the intestinal absorption of D-galactose via a complex cellular mechanism that could involve posttranscriptional regulation of the SGLT1 transporter.

Regulation of oxidative phosphorylation through parallel activation

When the mechanical work intensity in muscle increases, the elevated ATP consumption rate must be matched by the rate of ATP production by oxidative phosphorylation in order to avoid a quick exhaustion of ATP. The traditional mechanism of the regulation of oxidative phosphorylation, namely the negative feedback involving [ADP] and [Pi] as regulatory signals, is not sufficient to account for various kinetic properties of the system in intact skeletal muscle and heart in vivo. Theoretical studies conducted using a dynamic computer model of oxidative phosphorylation developed previously strongly suggest the so-called each-step-activation (or parallel activation) mechanism, due to which all oxidative phosphorylation complexes are directly activated by some cytosolic factor/mechanism related to muscle contraction in parallel with the activation of ATP usage and substrate dehydrogenation by calcium ions. The present polemic article reviews and discusses the growing evidence supporting this mechanism and compares it with alternative mechanisms proposed in the literature. It is concluded that only the each-step-activation mechanism is able to explain the rich set of various experimental results used as a reference for estimating the validity and applicability of particular mechanisms.

Oxygen delivery by blood determines the maximal VO2 and work rate during whole body exercise in humans: in silico studies

It has been proposed by Saltin (J Exp Biol 115: 345-354, 1985) that oxygen delivery by blood is limiting for maximal work and oxygen consumption in humans during whole body exercise but not during single-muscle exercise. To test this prediction quantitatively, we developed a static (steady-state) computer model of oxygen transport to and within human skeletal muscle during single-muscle (quadriceps) exercise and whole body (cycling) exercise. The main system fluxes, namely cardiac output and oxygen consumption by muscle, are described as a function of the "primary" parameter: work rate. The model is broadly validated by comparison of computer simulations with various experimental data. In silico studies show that, when all other parameters and system properties are kept constant, an increase in the working muscle mass from 2.5 kg (single quadriceps) to 15 kg (two legs) causes, at some critical work intensity, a drop in oxygen concentration in muscle cells to (very near) zero, and therefore oxygen supply by blood limits maximal oxygen consumption and oxidative ATP production. Therefore, the maximal oxygen consumption per muscle mass is significantly higher during single-muscle exercise than during whole body exercise. The effect is brought about by a distribution of a limited amount of oxygen transported by blood in a greater working muscle mass during whole body exercise.

Transporter regulator RS1 (RSC1A1) coats the trans-Golgi network and migrates into the nucleus

The product of gene RSC1A1, named RS1, is involved in transcriptional and posttranscriptional regulation of sodium-d-glucose cotransporter SGLT1, and removal of RS1 in mice led to an increase of SGLT1 expression in small intestine and to obesity (Osswald C, Baumgarten K, Stümpel F, Gorboulev V, Akimjanova M, Knobeloch K-P, Horak I, Kluge R, Joost H-G, and Koepsell H. Mol Cell Biol 25: 78-87, 2005). Previous data showed that RS1 inhibits transcription of SGLT1 in LLC-PK1 cells derived from porcine kidney. A decrease of the intracellular amount of RS1 protein was observed during cell confluence, which was paralleled by transcriptional upregulation of SGLT1. In the present study, the subcellular distributions of endogenously expressed RS1 and SGLT1 were compared in LLC-PK1 cells and human embryonic kidney (HEK)-293 cells using immunofluorescence microscopy. RS1 was located at the plasma membrane, at the entire trans-Golgi network (TGN), and within the nucleus. Treatment of LLC-PK1 cells with brefeldin A induced rapid release of RS1 from the TGN, and confluence of LLC-PK1 cells was accompanied by reduction of nuclear location of RS1; 84-90% of subconfluent cells and 5-34% of confluent cells contained RS1 in the nuclei. This suggests that confluence-dependent transcriptional inhibition by RS1 is partially regulated by nuclear migration. Furthermore, we assigned SGLT1 to microtubule-associated tubulovesicular structures and dynamin-containing parts of the TGN. The data indicate that RS1 inhibits the dynamin-dependent release of SGLT1-containing vesicles from the TGN.