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

It has been postulated previously that a direct activation of all oxidative phosphorylation complexes in parallel with the activation of ATP usage and substrate dehydrogenation (the so-called each-step activation) is the main mechanism responsible for adjusting the rate of ATP production by mitochondria to the current energy demand during rest-to-work transition in intact skeletal muscle in vivo. The present in silico study, using a computer model of oxidative phosphorylation developed previously, analyzes the impact of the each-step-activation mechanism on the distribution of control (defined within Metabolic Control Analysis) over the oxygen consumption flux among the components of the bioenergetic system in intact oxidative skeletal muscle at different energy demands. It is demonstrated that in the absence of each-step activation, the oxidative phosphorylation complexes take over from ATP usage most of the control over the respiration rate and oxidative ATP production at higher (but still physiological) energy demands. This leads to a saturation of oxidative phosphorylation, impossibility of a further acceleration of oxidative ATP synthesis, and dramatic drop in the phosphorylation potential. On the other hand, the each-step-activation mechanism allows maintenance of a high degree of the control exerted by ATP usage over the ATP turnover and oxygen consumption flux even at high energy demands and thus enables a potentially very large increase in ATP turnover. It is also shown that low oxygen concentration shifts the metabolic control from ATP usage to cytochrome oxidase and thus limits the oxidative ATP production. respiration rate; parallel activation; oxidative phosphorylation; metabolic control analysis; flux control coefficient; muscle contraction     

A model of oxidative phosphorylation in mammalian skeletal muscle

A dynamic computer model of oxidative phosphorylation in oxidative mammalian skeletal muscle was developed. The previously published model of oxidative phosphorylation in isolated skeletal muscle mitochondria was extended by incorporation of the creatine kinase system (creatine kinase plus phosphocreatine/creatine pair), cytosolic proton production/consumption system (proton production/consumption by the creatine kinase-catalysed reaction, efflux/influx of protons), physiological size of the adenine nucleotide pool and some additional minor changes. Theoretical studies performed by means of the extended model demonstrated that the CK system, which allows for large changes in P(i) in relation to isolated mitochondria system, has no significant influence on the kinetic properties of oxidative phosphorylation, as inorganic phosphate only slightly modifies the relationship between the respiration rate and [ADP]. Computer simulations also suggested that the second-order dependence of oxidative phosphorylation on [ADP] proposed in the literature refers only to the ATP synthesis flux, but not to the oxygen consumption flux (the difference between these two fluxes being due to the proton leak). Next, time courses of changes in fluxes and metabolite concentrations during transition between different steady-states were simulated. The model suggests, in accordance with previous theoretical predictions, that activation of oxidative phosphorylation by an increase in [ADP] can (roughly) explain the behaviour of the system only at low work intensities, while at higher work intensities parallel activation of different steps of oxidative phosphorylation is involved.     

Influence of substrate activation (hydrolysis of ATP by first steps of glycolysis and beta-oxidation) on the effect of enzyme deficiencies,inhibitors, substrate shortage and energy demand on oxidative phosphorylation

In intact tissues respiratory substrates (glucose, fatty acids) must be activated with the use of ATP before they may be oxidised and used for energy (ATP) production. This activation by product constitutes an example of a typical positive feedback. In the present paper, the influence of substrate activation on the effect of inborn enzyme deficiencies, inhibitors, lowered oxygen tension, respiratory fuel shortage and increased energy demand on respiration and ATP synthesis is studied with the aid of the dynamic computer model of oxidative phosphorylation in isolated mitochondria developed previously. Computer simulations demonstrate that, in the case where oxidative phosphorylation in the whole organism is partially inhibited, the necessity of substrate activation can have significant impact on the relationship between the activity of (particular steps of) oxidative phosphorylation (or the value of energy demand) and the respiration rate. Depending on the sensitivity of ATP usage to ATP concentration, substrate activation may either slightly enhance the effect of the decrease in the oxidative phosphorylation activity (increase in energy demand) or may lead to a non-stability and sudden collapse of the respiration rate and phosphorylation potential below (above) a certain threshold value of oxidative phosphorylation activity (energy demand). This theoretical finding suggests a possible causal relationship between the affinity of ATP usage to [ATP] and the tissue specificity of mitochondrial diseases.     

Regulation of oxidative phosphorylation in intact mammalian heart in vivo

A dynamic computer model of oxidative phosphorylation in intact heart was developed by modifying the model of oxidative phosphorylation in intact skeletal muscle published previously. Next, this model was used for theoretical studies on the regulation of oxidative phosphorylation in intact heart in vivo during transition between different work intensities. It is shown that neither a direct activation of ATP usage alone nor a direct activation of both ATP usage and substrate dehydrogenation, including the calcium-activated tricarboxylate acid cycle dehydrogenases, can account for the constancy of [ADP], [PCr], [P(i)] and [NADH] during a significant increase in oxygen consumption and ATP turnover encountered in intact heart in vivo. Only a direct activation of all oxidative phosphorylation complexes in parallel with a stimulation of ATP usage and substrate dehydrogenation enabled to reproduce the experimental data concerning the constancy of metabolite concentrations. The molecular background of the differences between heart and skeletal muscle in the kinetic behaviour of the oxidative phosphorylation system is also discussed.     

Regulation of ATP supply in mammalian skeletal muscle during resting state-->intensive work transition

In the present debating paper, the problem how the rate of ATP supply by oxidative phosphorylation in mitochondria is adjusted to meet a greatly increased demand for ATP during intensive exercise of skeletal muscle is discussed. Different experimental results are collected from different positions of the literature and confronted with five conceptual models of the regulation of the oxidative phosphorylation system. The previously performed computer simulations using a dynamic model of oxidative phosphorylation are also discussed in this context. The possible regulatory mechanisms considered in the present article are: (A) output activation: an external effector activates directly only the output of the system (ATP turnover); (B) input/output activation: an external effector activates directly the output (ATP usage) and input (substrate dehydrogenation) of the system; (C) removal of substrate shortage: only ATP consumption and substrate supply by blood are directly activated; (D) removal of oxygen shortage: only ATP consumption and oxygen supply by blood are directly activated; (E) each step activation: an external effector activates both the ATP-consuming subsystem and all the steps in the ATP-producing subsystem (particular enzymes/carriers/blocks of oxidative phosphorylation, substrate supply, oxygen supply). The performed confrontation of the considered mechanisms with the presented results leads to the conclusion that only the each step activation model is quantitatively consistent with the whole set of experimental data discussed. It is therefore postulated that a universal effector/regulatory mechanism of a still unknown nature which activates all steps of oxidative phosphorylation should exist and be discovered. A possible nature of such an effector is shortly discussed.     

Effect of enzyme deficiencies on oxidative phosphorylation: from isolated mitochondria to intact tissues. Theoretical studies

The present article briefly summarizes the theoretical studies made by the authors and co-workers on the effect of inborn enzyme deficiencies on oxidative phosphorylation in intact tissues and on the genesis of mitochondrial diseases. The dynamic computer model of oxidative phosphorylation developed previously allowed to extrapolate experimental data (especially: threshold curves describing the dependence of oxygen consumption and ATP turnover on activities/concentrations of particular oxidative phosphorylation enzymes) obtained for isolated muscle mitochondria in state 3 at saturating oxygen concentrations to more physiological conditions prevailing in intact tissues. In particular, theoretical studies demonstrated that the threshold value of the relative activity/concentration of a given mitochondrial complex, below which a significant decrease in the respiration rate takes place, increases with an increase in energy demand. This fact was proposed as a possible explanation of the tissue specificity of mitochondrial diseases. Additionally, a decreased oxygen concentration was shown to increase the threshold value (and flux control coefficient) for cytochrome oxidase. We subsequently developed a model called binary mitochondria heteroplasmy, in which there are only two subpopulations of mitochondria: one wild-type and one containing only defected molecules of a given enzyme. In this model we show that a defect has a pronounced effect on oxidative phosphorylation, significantly increasing the threshold value. It was also proposed that a parallel activation in the ATP supply-demand system during an increased energy demand significantly lessens the effect of enzyme deficiencies on oxidative phosphorylation (decreases the threshold value). Finally, the necessity of substrate activation may lead to an instability in the system and to appearance of a second threshold, below which respiration suddenly drops to zero, which is equivalent to the energetic death of a cell.     

Theoretical studies on the regulation of oxidative phosphorylation in intact tissues

The theoretical studies on the regulation of oxidative phosphorylation that were performed with the aid of kinetic models of this process are overviewed. A definition of the regulation of the flux through a metabolic pathway is proposed and opposed to the control exerted by particular enzymes over this flux. Different kinetic models of oxidative phosphorylation proposed in the literature are presented, of which only the model proposed by myself and co-workers was extensively used in theoretical studies on the regulation and compensation in the oxidative phosphorylation system. These theoretical studies have led to the following conclusions: (1) in isolated mitochondria, an increase in the activity of an artificial ATP-using system stimulates mitochondria mainly via changes in [ADP], while changes in [ATP] and [P(i)] play only a minor role; (2) in non-excitable tissues (e.g. liver), hormones (acting via some cytosolic factor(s)) activate directly both ATP usage and at least some enzymes of the ATP-producing block; (3) in excitable tissues (e.g. skeletal muscle), neural signals stimulate (via some cytosolic factor(s)) in parallel all the steps of oxidative phosphorylation together with ATP usage and substrate dehydrogenation; (4) the decrease in the flux through cytochrome oxidase caused by a decrease in oxygen concentration is, at least partially, compensated by a decrease in Delta p and increase in the reduction level of cytochrome c. A theoretical prediction is formulated that there should exist and be observable a universal cytosolic factor/regulatory mechanism which directly activates (at least in excitable tissues) all complexes of oxidative phosphorylation during an increased energy demand.     

Oxygen consumption and metabolite concentrations during transitions between different work intensities in heart

Steady-state metabolite (ADP, ATP, P(i), PCr, and NADH) concentrations usually differ little between different workloads with significantly different oxygen consumption rates in the heart. However, during transitions between steady states, metabolite concentrations may in some cases change transiently, exhibiting a significant overshoot or undershoot, whereas in other cases they approach near-exponentially new steady-state values. Oxygen consumption rate usually reaches the new steady-state value very quickly (within a few seconds). The present in silico studies, performed using a previously developed computer model of oxidative phosphorylation in the heart, demonstrate that such a behavior of the oxidative phosphorylation system can be reproduced only under the assumption that ATP usage, substrate dehydrogenation, and (particular steps of) oxidative phosphorylation are directly activated to a similar extend by some cytosolic factor/mechanism during transition from low work to high work (the so-called parallel-activation mechanism). Computer simulations show that some differences observed between different experimental systems can be explained by a slightly different balance of the activation of particular components of the system and/or by a delay in time of the activation/inactivation of substrate dehydrogenation and oxidative phosphorylation during low-to-high and high-to-low work transitions. Thus the presented theoretical approach offers a general idea that is able to unify, at least semiquantitatively, different experimental data available in the literature.     

Computer-aided studies on the regulation of oxidative phosphorylation during work transitions

In the present polemic paper the application of computer models of oxidative phosphorylation (OXPHOS) in heart, skeletal muscle and liver to the studies on the regulation of the bioenergetic system in intact cells during work transitions is discussed. Two groups of such models are compared: group I models that involve only a direct activation of ATP usage by Ca²⁺, and group II models that assume a direct activation by some (probably) Ca²⁺-related mechanism of essentially all steps of the system. It is argued that group II models reproduce much better a broad range of variable values and system properties encountered in experimental studies. The consequences of the theoretical and experimental development of Metabolic Control Analysis, within the framework of which it has been shown that the control over the flux through the oxidative phosphorylation system is shared by essentially all components of this system, are analyzed. In particular, it is argued that in order to increase the flux very significantly, and at the same time to maintain relatively constant concentrations of such intermediate metabolites as ADP, ATP, P_i, PCr and NADH, it is necessary to activate directly many, if not all components of the system (the ‘multi-step parallel activation’ mechanism). Generally, it is suggested that this is not a particular form or complexity of computer models, but rather their agreement with a broad range of experimental data concerning ‘macroscopic’ system properties that really matters.The specificity of the regulation of the energetic system of pancreatic β-cells is discussed.     

Theoretical studies on control of oxidative phosphorylation in muscle mitochondria at different energy demands and oxygen concentrations

The mathematical dynamic model of oxidative phosphorylation in muscle mitochondria developed previously was used to calculate the flux control coefficients of particular steps of this process in isolated mitochondria at different amounts of hexokinase and oxygen concentrations. The pattern of control was completely different under different conditions. For normoxic concentration, the main controlling steps in state 4, state 3.5 and state 3 were proton leak, ATP usage (hexokinase) and complex III, respectively. The pattern of control in state 4 was not changed at hypoxic oxygen concentration, while in state 3.5 and state 3 much of the control was shifted from other steps to cytochrome oxidase. The implications of the theoretical results obtained for the regulation of oxidative phosphorylation in intact muscle are discussed.     

Energy metabolism transition in multi-cellular human tumor spheroids

It is thought that glycolysis is the predominant energy pathway in cancer, particularly in solid and poorly vascularized tumors where hypoxic regions develop. To evaluate whether glycolysis does effectively predominate for ATP supply and to identify the underlying biochemical mechanisms, the glycolytic and oxidative phosphorylation (OxPhos) fluxes, ATP/ADP ratio, phosphorylation potential, and expression and activity of relevant energy metabolism enzymes were determined in multi-cellular tumor spheroids, as a model of human solid tumors. In HeLa and Hek293 young-spheroids, the OxPhos flux and cytochrome c oxidase protein content and activity were similar to those observed in monolayer cultured cells, whereas the glycolytic flux increased two- to fourfold; the contribution of OxPhos to ATP supply was 60%. In contrast, in old-spheroids, OxPhos, ATP content, ATP/ADP ratio, and phosphorylation potential diminished 50-70%, as well as the activity (88%) and content (3 times) of cytochrome c oxidase. Glycolysis and hexokinase increased significantly (both, 4 times); consequently glycolysis was the predominant pathway for ATP supply (80%). These changes were associated with an increase (3.3 times) in the HIF-1alpha content. After chronic exposure, both oxidative and glycolytic inhibitors blocked spheroid growth, although the glycolytic inhibitors, 2-deoxyglucose and gossypol (IC(50) of 15-17 nM), were more potent than the mitochondrial inhibitors, casiopeina II-gly, laherradurin, and rhodamine 123 (IC(50) > 100 nM). These results suggest that glycolysis and OxPhos might be considered as metabolic targets to diminish cellular proliferation in poorly vascularized, hypoxic solid tumors.     

Paraquat increases cyanide-insensitive respiration in murine lung epithelial cells by activating an NAD(P)H:paraquat oxidoreductase: identification of the enzyme as thioredoxin reductase

Pulmonary fibrosis is one of the most severe consequences of exposure to paraquat, an herbicide that causes rapid alveolar inflammation and epithelial cell damage. Paraquat is known to induce toxicity in cells by stimulating oxygen utilization via redox cycling and the generation of reactive oxygen intermediates. However, the enzymatic activity mediating this reaction in lung cells is not completely understood. Using self-referencing microsensors, we measured the effects of paraquat on oxygen flux into murine lung epithelial cells. Paraquat (10-100 microm) was found to cause a 2-4-fold increase in cellular oxygen flux. The mitochondrial poisons cyanide, rotenone, and antimycin A prevented mitochondrial- but not paraquat-mediated oxygen flux into cells. In contrast, diphenyleneiodonium (10 microm), an NADPH oxidase inhibitor, blocked the effects of paraquat without altering mitochondrial respiration. NADPH oxidases, enzymes that are highly expressed in lung epithelial cells, utilize molecular oxygen to generate superoxide anion. We discovered that lung epithelial cells possess a distinct cytoplasmic diphenyleneiodonium-sensitive NAD(P)H:paraquat oxidoreductase. This enzyme utilizes oxygen, requires NADH or NADPH, and readily generates the reduced paraquat radical. Purification and sequence analysis identified this enzyme activity as thioredoxin reductase. Purified paraquat reductase from the cells contained thioredoxin reductase activity, and purified rat liver thioredoxin reductase or recombinant enzyme possessed paraquat reductase activity. Reactive oxygen intermediates and subsequent oxidative stress generated from this enzyme are likely to contribute to paraquat-induced lung toxicity.     

Glutamine‐driven oxidative phosphorylation is a major ATP source in transformed mammalian cells in both normoxia and hypoxia

Mammalian cells can generate ATP via glycolysis or mitochondrial respiration. Oncogene activation and hypoxia promote glycolysis and lactate secretion. The significance of these metabolic changes to ATP production remains however ill defined. Here, we integrate LC‐MS‐based isotope tracer studies with oxygen uptake measurements in a quantitative redox‐balanced metabolic flux model of mammalian cellular metabolism. We then apply this approach to assess the impact of Ras and Akt activation and hypoxia on energy metabolism. Both oncogene activation and hypoxia induce roughly a twofold increase in glycolytic flux. Ras activation and hypoxia also strongly decrease glucose oxidation. Oxidative phosphorylation, powered substantially by glutamine‐driven TCA turning, however, persists and accounts for the majority of ATP production. Consistent with this, in all cases, pharmacological inhibition of oxidative phosphorylation markedly reduces energy charge, and glutamine but not glucose removal markedly lowers oxygen uptake. Thus, glutamine‐driven oxidative phosphorylation is a major means of ATP production even in hypoxic cancer cells.     

Effects of sodium and potassium ions on oxidative phosphorylation in relation to respiratory control by a cell-membrane adenosine triphosphatase

1. A study has been made of the oxygen consumption of kidney homogenates in relation to the ADP concentration as regulated by the cell-membrane adenosine triphosphatase. Stimulation of this enzymic activity by Na(+) and K(+) caused parallel increases in oxygen consumption and ADP concentration. Similarly, inhibition with ouabain caused a parallel fall. The membrane adenosine triphosphatase concerned in active transport therefore appears to regulate respiration through its control of ADP concentration. 2. The respiration of homogenates and mitochondria was also stimulated by K(+) in a way independent of adenosine-triphosphatase activity. It was shown that K(+) facilitates oxidative phosphorylation and the respiratory response to ADP. A K(+) concentration of 25-50mm was needed for maximum oxidative phosphorylation in the presence of physiological concentration of Na(+). Na(+) counteracted K(+) in the effects on mitochondria. It is concluded that K(+) regulates cellular respiration at two structures, one directly in mitochondria, and the second indirectly through control of ADP production at the cell membrane.     

Rate-controlling steps of oxidative phosphorylation in rat liver mitochondria. A synoptic approach of model and experiment

The contribution of different steps to the control of oxidative phosphorylation in isolated rat liver mitochondria was investigated by a combination of experiments and computer simulations. The parameters of the mathematical model of phosphorylating mitochondria were derived from experimental data. The model correctly describes the competition between ATP utilization inside and outside mitochondria for the ATP generated in mitochondria. On the basis of the good agreement between experiments and simulations, the contribution of different steps to the control of respiration was estimated by computing their control strengths, i.e., the influence of their activities on the rate of respiration. The rate-controlling influences vary depending on the load of oxidative phosphorylation. The predominant steps are: in the fully active state (State 3) — the hydrogen supply to the respiratory chain; in the resting state (State 4) — the proton leak of the mitochondrial inner membrane; in states of non-maximum ATP export — the adenine nucleotide translocator. Titrations of respiration with phenylsuccinate, antimycin, oligomycin and carboxyatractyloside completely support these conclusions.     

Some factors determining the PCr recovery overshoot in skeletal muscle

It has been proposed recently that the phosphocreatine (PCr) overshoot (increase above the resting level) during muscle recovery after exercise is caused by a slow decay during this recovery of the direct activation of oxidative phosphorylation taking place during muscle work. In the present article the factors determining the appearance and size of the PCr overshoot are studied using the computer model of oxidative phosphorylation in intact skeletal muscle developed previously. It is demonstrated that the appearance and duration of this overshoot is positively correlated with the value of the characteristic decay time of the direct activation of oxidative phosphorylation. It is also shown that the size of PCr overshoot is increased by low resting PCr/Cr ratio (what is confirmed by our unpublished experimental data), by high intensity of the direct activation of oxidative phosphorylation, by high muscle work intensity and by low rate of the return of cytosolic pH to the resting value during muscle recovery.     

Nicotine activates the extracellular signal-regulated kinase 1/2 via the alpha7 nicotinic acetylcholine receptor and protein kinase A, in SH-SY5Y cells and hippocampal neurones

Neuronal nicotinic acetylcholine receptors (nAChR) can modulate many cellular mechanisms, such as cell survival and memory processing, which are also influenced by the serine/threonine protein kinases ERK1/2. In SH-SY5Y cells and hippocampal neurones, nicotine (100 microM) increased the activity of ERK1/2. This effect was Ca2+ dependent, and prevented by the alpha7 nAChR antagonist alpha-bungarotoxin (alpha-Bgt) and an inhibitor (PD98059) of the upstream kinase MEK. To determine the intervening steps linking Ca2+ entry to MEK-ERK1/2 activation, inhibitors of Ca2+-dependent kinases were deployed. In SH-SY5Y cells, selective blockers for PKC (Ro 31-8220), CaM kinase II (KN-62) or PI3 kinase (LY 294002) failed to inhibit the nicotine-evoked increase in ERK1/2 activity. In contrast, two structurally different inhibitors of PKA (KT 5720 and H-89) completely prevented the nicotine-dependent increase in ERK1/2 activity. Inhibition of the nicotine-evoked increase in ERK1/2 activity by H-89 was also observed in hippocampal cultures. Down stream of PKA, the activity of B-Raf was significantly decreased by nicotine in SH-SY5Y cells, as determined by direct measurement of MEK1 phosphorylation or in vitro kinase assays, whereas the modulation of MEK1 phosphorylation by Raf-1 tended to increase. Thus, this study provides evidence for a novel signalling route coupling the stimulation of alpha7 nAChR to the activation of ERK1/2, in a Ca2+ and PKA dependent manner.     

The nature of methanol-induced acidification of the medium by products of metabolism of methylotrophic yeasts

The washed cells of Hansenula polymorpha and Pichia pinus grown on the medium with methanol rapidly acidify the medium during incubation with the mentioned alcohol or formaldehyde. It is found that proton extrusion is coupled with formate anion efflux. Acidification is proved to be energy-dependent process since it is inhibited by respiration poisons, uncouplers of oxidative phosphorylation, and by ATPase inhibitors.     

Role of oxidative phosphorylation and ATP release in mediating birth-related pulmonary vasodilation in fetal lambs

We investigated the hypothesis that birth-related pulmonary vasodilation is mediated in part by an increase in oxidative phosphorylation and ATP release in response to oxygen exposure at birth. Studies were done in fetal lambs to evaluate the independent effects of oxygen, lung distension alone, or lung distension accompanied by oxygenation and shear stress on fetal pulmonary blood flow and resistance and plasma ATP levels in the pulmonary artery. The effect of each intervention was evaluated in lambs assigned to one of three groups: control or pretreatment with 2,4-dinitrophenol or antimycin-A, inhibitors of oxidative phosphorylation. Exposure to oxygen alone or with lung distension was associated with increases in plasma ATP levels and pulmonary blood flow and a decrease in pulmonary vascular resistance. Plasma ATP levels did not change during lung distension alone. 2,4-Dinitrophenol and antimycin-A attenuated the pulmonary vasodilator response to oxygen but did not attenuate the response to lung distension alone. An increase in oxidative phosphorylation and ATP release during oxygen exposure may contribute to birth-related pulmonary vasodilation in fetal lambs.     

Effects of thyroid hormone on mitochondrial oxidative phosphorylation.

In order to locate sites of action of thyroid hormone on mitochondrial oxidative phosphorylation we have used an experimental application of control analysis as previously described [Groen, Wanders, Westerhoff, Van der Meer & Tager (1982) J. Biol. Chem. 257, 2754-2757]. Rat-liver mitochondria were isolated from hypothyroid rats or from hypothyroid rats 24 h after treatment with a single dose of 3,3',5-triiodothyronine (T3). The amount of control exerted by four different steps on State-3 respiration with succinate as respiratory substrate was quantified by using specific inhibitors. The hormone treatment resulted in an increase in the flux control coefficient of the adenine nucleotide translocator, the dicarboxylate carrier and cytochrome c oxidase and a decrease in the flux control coefficient of the bc1-complex. The results of this analysis indicate that thyroid hormone treatment results in an activation of the bc1-complex and of at least one other enzyme, possibly succinate dehydrogenase. Measurement of the extramitochondrial ATP/ADP ratio at different rates of respiration (induced by addition of different amounts of hexokinase in the presence of glucose and ATP) showed that the adenine nucleotide translocator operates at a higher (ATP/ADP)out after T3 treatment, which supports previous reports on stimulation of this step by thyroid hormone.