AJP - Regu AJP citation statistics
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Am J Physiol Regul Integr Comp Physiol 280: R695-R704, 2001;
0363-6119/01 $5.00
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (28)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ramakrishna, R.
Right arrow Articles by Palsson, B. O.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ramakrishna, R.
Right arrow Articles by Palsson, B. O.
Vol. 280, Issue 3, R695-R704, March 2001

Flux-balance analysis of mitochondrial energy metabolism: consequences of systemic stoichiometric constraints

Ramprasad Ramakrishna, Jeremy S. Edwards, Andrew McCulloch, and Bernhard O. Palsson

Department of Bioengineering, University of California San Diego, La Jolla, California 92093 - 0412


    ABSTRACT
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Mitochondrial metabolism is a critical component in the functioning and maintenance of cellular organs. The stoichiometry of biochemical reaction networks imposes constraints on mitochondrial function. A modeling framework, flux-balance analysis (FBA), was used to characterize the optimal flux distributions for maximal ATP production in the mitochondrion. The model predicted the expected ATP yields for glucose, lactate, and palmitate. Genetic defects that affect mitochondrial functions have been implicated in several human diseases. FBA can characterize the metabolic behavior due to genetic deletions at the metabolic level, and the effect of mutations in the tricarboxylic acid (TCA) cycle on mitochondrial ATP production was simulated. The mitochondrial ATP production is severely affected by TCA-cycle mutations. In addition, the model predicts the secretion of TCA-cycle intermediates, which is observed in clinical studies of mitochondriopathies such as those associated with fumarase deficiency. The model provides a systemic perspective to characterize the effect of stoichiometric constraints and specific metabolic fluxes on mitochondrial function.

mitochondria; flux analysis; adenosine 5'-phosphate production; mitochondrial disease


    INTRODUCTION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

CELLULAR RESPIRATION and, consequently, mitochondrial metabolism has a significant role in the functions of aerobic organs. In systems such as the heart, this process supports myofilament contraction, transmembrane ion, and intracellular calcium cycling. Normal well-perfused myocardium generates >90% of its ATP by oxidative metabolism and <10% by anaerobic glycolysis (11). Mitochondrial metabolism also plays a critical role in the function of other organs, such as the liver and the brain, where the impaired functioning of mitochondria has been implicated in several neurological disorders (16).

The link between the genetics and physiology of mitochondrial diseases is an area of intensive study. The genetic basis of mitochondrial disease is considered to arise from defects of nuclear DNA, including defects of protein import, defects of mitochondrial DNA (mtDNA), such as point mutations, deletions, and duplications, and defects of communication between nuclear and mitochondrial genomes (multiple deletions and mtDNA deletion) (5). The physiological manifestations of these genetic defects affect the normal functioning of substrate transport, substrate use, the tricarboxylic (TCA) cycle, the respiratory chain, and oxidation/phosphorylation coupling. mtDNA has no introns and has its own independent replication, transcription, and translation systems. Mutations in mtDNA accumulate 10-20 times faster than in comparable nuclear genes (17). The accumulation of mtDNA mutations is believed to be a key factor in aging and the progression of mitochondrial diseases (15).

Mitochondrial metabolism results from the concerted action of several biochemical reactions that are coordinately regulated. Genotypic and physiological factors interactively contribute to adversely affect mitochondrial metabolism (23). Mutations in the genetic content can lead to changes in enzyme expression and/or activity and thus alter mitochondrial function by changing fluxes of important metabolic reactions and consequently affecting the normal physiological behavior of cells, tissues, and organs. A genetic mutation that changes a flux or fluxes consequently forces the metabolic system to a new state and affects the attainment of its physiological objective. Therefore, an understanding of the systemic constraints on metabolism can provide significant insight into physiological function. The framework of flux-balance analysis (FBA) (2, 6, 20, 22) is a powerful methodology for the analysis of metabolic systems. This approach has been successfully applied to the analysis of adipocyte metabolism (6), hybidoma metabolism (1, 14), and bacterial growth (21), in which the constraints imposed by biochemical reaction stoichiometry are systemically modeled.

The application of systemic metabolic simulation approaches can significantly aid in understanding the metabolic basis of mitochondrial-linked disease. The biochemical reactions involved in mitochondrial energy metabolism function as a coordinated network subject to stoichiometric and regulatory constraints. The objective of this study was to first develop a model for mitochondrial function, using the methods of systems theory that can predict the effect of genotype changes on certain aspects of mitochondrial function. The model was to comprise of the metabolic pathways that constitute a functioning mitochondrion, including the reactions in the mitochondrial matrix, contributing reactions in the cytosol, and key shuttles. The analysis considers the constraints imposed by the stoichiometry of these biochemical reaction networks on the achievement of specific metabolic objectives. The model will be used to characterize energy metabolism on the standard substrates of mitochondrial metabolism.

As part of the model development, the genes(s) involved in the metabolic reactions considered will be detailed. This provides a framework for studying the effect of genetic changes on metabolism. The model will be used to characterize the effect of genetic mutations on mitochondrial metabolism. These simulations will be compared with literature that describes the physiology of mitochondrial disease.


    METHODS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

FBA. FBA is a modeling framework developed to characterize the capabilities and properties of a metabolic network. The network comprises the metabolites and the reactions they are involved in, including their formation and degradation, transport, and cellular utilization. For every metabolite (Xi), a material balance is derived as follows
<FR><NU>d<IT>X<SUB>i</SUB></IT></NU><DE>d<IT>t</IT></DE></FR><IT>=</IT><LIM><OP>∑</OP><LL><IT>j</IT></LL></LIM><IT> s<SUB>ij</SUB>v<SUB>j</SUB>−b<SUB>i</SUB></IT>
where sij is the stoichiometric coefficient associated with each flux, vj and bi, the net transport flux of Xi. This material balance under steady-state conditions will reduce to the algebraic relation
<LIM><OP>∑</OP><LL><IT>j</IT></LL></LIM><IT> s<SUB>ij</SUB>v<SUB>j</SUB>−b<SUB>i</SUB>=0</IT>
or, for all intermediates simultaneously at steady state, the individual balance equations can be rewritten in matrix form
S<IT>×</IT>v<IT>−</IT>b<IT>=0</IT> (1)
where S is the stoichiometric matrix (m rows and n columns), v is the vector of n metabolic fluxes, and b is the vector representing m transport fluxes (i.e., known consumption rates, by-product production rates, and uptake rates). The matrix S is not square (n m), and there exists a plurality of solutions for equation 1. That is, there can be a number of feasible flux distributions satisfying these stoichiometric constraints, each representing a particular metabolic state. Therefore, the null space, or the set of all feasible flux distributions, represents the capabilities of the metabolic genotype. The transport fluxes represent environmental conditions that, along with the genotype, define the metabolic state. However, obtaining all possible metabolic states for any genotype-environment interaction depends on how well the genotype and environmental factors are characterized.

The question, then, is which of these feasible metabolic states is manifested in the biological system under consideration. FBA postulates that the metabolic system exhibits a metabolic state that is optimal under some criteria. This objective is expressed as a linear combination of the fluxes contained in v. The model can then be formulated as a linear programming problem as follows
<AR><R><C>minimize (maximize) Z = <LIM><OP>∑</OP></LIM> c<SUB>i</SUB><IT>v<SUB>i</SUB></IT></C></R><R><C>such that</C></R><R><C>S<IT>×</IT>v<IT>−</IT>b<IT>=0, 0≤&bgr;</IT><SUB>i</SUB><IT>≤v<SUB>i</SUB>≤&agr;</IT><SUB>i</SUB></C></R></AR>
Z is the objective function, representing a phenotypic property, and c is a vector of weights that are either costs of or benefits derived from the fluxes. The limits, alpha i and beta i, represent known constraints on the maximum or minimum values that the fluxes assume. When the maximization of ATP production is considered, the net flux of ATP hydrolysis (vATP_PR) is maximized.
ATP <LIM><OP><ARROW>→</ARROW></OP><UL><IT>v<SUB>ATP_PR</SUB></IT></UL></LIM> ADP<IT>+</IT>P<SUB>i</SUB>

maximize Z<IT>=v<SUB>ATP_PR</SUB></IT>
The balance equations expressed by equation 1 form part of the constraints. A commercially available package was used to solve the LP problems (Lindo Systems, Chicago, IL). The objective Z is minimized or maximized, subject to the imposed constraints. There are two parameters associated with the linear programming (LP) problem that help characterize the optimal solution as well as help in post-optimality analysis. These are the shadow price (gamma i) and the reduced cost. The gamma i are associated with each metabolite and are defined as gamma i = dZ/dbi. For the objective of maximization of ATP production, if the gamma i of NADH is 3.00, it means that an additional molecule of NADH can generate three more molecules of ATP. The gamma i, therefore, represents the increase in the value of the objective with the addition of the associated intermediate. The reduced costs are associated with each flux (vi) and signify the amount by which the objective function is decreased if vi is brought into the basis solution. For instance, if the input flux of lactate shows a reduced cost of 17.5, it means that increasing that flux by one unit will increase ATP production by 17.5 units. Reduced costs and shadow prices are terms commonly used in translating LP solutions to real-life situations and carrying out an analysis of alternate solutions from the original solution. We introduce these quantities to analyze the physiological task of energy metabolism to gain insight into the solutions as well as provide a "microeconomic" perspective to the ATP-generation problem.

The mitochondrial model. The flux balance model for mitochondria presented here comprises the glycolytic pathways, TCA cycle, and the electron transport system (ETS). The pentose phosphate pathway has not been included, because its activity is believed to be quite low for mitochondria-related functions (19). The oxidative metabolism of substrates takes place in the mitochondria; thus the substrates, metabolites, and cofactors must cross the selectively permeable membrane that separates the mitochondrial space from the cytosolic space. The reactions comprising the TCA cycle occur in the mitochondrion and produce the reduced coenzymes NADH and FADH2 that transfer electrons to oxygen in a regulated manner. Shuttles play an important role in transporting reducing equivalents generated in the cytosol into the mitochondrion.

The reactions, which make up the model, are divided into three sets, based on whether they occur in the cytosol or in the mitochondria or in transporting an intermediate across the mitochondrial membrane. The enzymes considered in the FBA model and their E.C numbers are listed in Table 1. The glycolytic reactions take place in the cytosol. The reactions in the TCA cycle and the ETS take place in the mitochondrial matrix. This compartmentalization of reactions in an organelle leads to the need for special reaction sequences, or shuttles, to transport reducing equivalents generated in the cytosol into the mitochondria. There are two such shuttles that have been found to be active in mitochondria.

                              
View this table:
[in this window]
[in a new window]
 
Table 1.   The fluxes in the FBA model

The ubiquitous malate-aspartate shuttle transports external NADH into the mitochondria. The coordinated exchange of malate and alpha -ketoglutarate with aspartate and glutamate by respective antiporters is a key feature of the shuttle (19, 24). The malate dehydrogenase reaction functions in opposite directions in the cytosol and in the matrix. The cytosolic reaction forms malate from oxaloacetate and involves the oxidation of NADH. Malate is then transported into the matrix and through the matrix with concomitant efflux of alpha -ketoglutarate. Malate is used to form oxaloacetate in the matrix, and NADH is formed in this reaction. The cycle is completed by the transamination reactions involving oxaloacetate, alpha -ketoglutarate, glutamate, and aspartate. The normal functioning of this shuttle is outlined in Fig. 1.


View larger version (31K):
[in this window]
[in a new window]
 
Fig. 1.   The normal role of the malate-aspartate shuttle in mitochondria. OAA, oxaloacetate.

An alternate shuttle that occurs in the skeletal muscle and brain is the glycerol-3-phosphate shuttle (8). Here, the reducing equivalents from NADH are delivered as FADH2 into the mitochondria. Thus the energy yield is reduced compared with the malate-aspartate shuttle, because each FADH2 produces only 2 ATP compared with 3 ATP for NADH. There are two dehydrogenases involved in this shuttle, the cytosolic and mitochondrial glycerol-3-phosphate dehydrogenase. The cytosolic dehydrogenase reduces dihydroxyacetone phosphate to glycerol-3-phosphate using NADH, and the mitochondrial dehydrogenase mediates the reverse reaction, which generates FADH2.

Additionally, fluxes for substrate input and intermediate secretion are also included in the model. The input fluxes for the substrates can be constrained for the different situations that will be studied.

The postulated objective for mitochondrial energy metabolism. The flux-balance model as explained earlier, requires an objective, which the cell or the organelle, as is the case here, is attempting to achieve. This objective is generally a postulated theoretical assumption that completes the model formulation and enables us to simulate metabolic behavior. For this analysis, maximizing the production of ATP is chosen as the objective function. Therefore, a single flux that hydrolyses the net ATP produced was considered, and the objective function was mathematically represented as vATP_USE, with a corresponding cost of unity. Cairns et al. (4) considered the functional basis for control of mitochondrial oxidative phosphorylation in different organs for rat mitochondria. ATP production in the brain and heart mitochondrial systems was found to use more oxygen but produce ATP at a faster rate than liver systems. They attribute these qualities to the thermodynamic degree of oxidative coupling. The general conclusion is that maximizing the rate of energy production, rather than maintaining thermodynamic efficiency, is the important characteristic of mitochondria in physiological systems. This conclusion supports the postulate that the objective of mitochondrial metabolism is the maximization of ATP production. This postulate requires several detailed studies to be experimentally validated. In this modeling study, its role is to provide a complete theoretical framework to draw insights into the constraints imposed by stoichiometry on energy metabolism.


    RESULTS
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Substrate preferences. The flux distributions for optimal production of ATP, shown in Fig. 2, A-C, were determined under conditions in which the substrate uptake is restricted to 1 mol · unit time-1 · unit mass dry wt-1 or one unit. The simulations were carried out for the three substrates, glucose, lactate, and palmitic acid. Table 2 lists the constraints imposed on the different fluxes. The complete utilization of 1 mol of glucose results in the formation of 38 ATP with the concomitant utilization of 6 mol of O2. The utilization of 1 mol of lactate forms 17.5 ATP with the utilization of 3.0 mol of O2, and the utilization of palmitic acid (a 14-C fatty acid) produces 129 ATP but requires 23 mol of O2. Glucose is the preferred substrate compared with lactate and 14-carbon fatty acid when the oxygen flux is restricted and all three substrates are made available. Because glucose uptake produces maximum ATP per mole of oxygen consumed, ATP synthesis from glucose is the optimal strategy for energy metabolism.


View larger version (18K):
[in this window]
[in a new window]
 
Fig. 2.   The optimal metabolic flux distributions for the maximal production of ATP from glucose (Glu; A), lactate (B), and palmitate (C). The figures have been simplified by representing linear sequences of fluxes as a single flux and only showing the key branch points. The cytosolic and mitochondrial metabolites have been identified by the letters c and m, respectively. The malate-asparate shuttle is active for the metabolism of lactate and glucose, and there is an efflux of alpha -ketoglutarate (AKG) out of the mitochondria accompanied by an influx of malate (Mal). PEP, phosphoenolpyruvate; ETS, electron transport system.


                              
View this table:
[in this window]
[in a new window]
 
Table 2.   Constraints on key fluxes for the model simulations

Table 3 compares the three substrates in terms of the ATP yield on a per-oxygen basis and a per-carbon basis. Whereas the fatty acid provides more ATP on a per-carbon basis, the cost in terms of oxygen uptake makes the fatty acid a suboptimal alternative to the other substrates. The reduced costs contain information on the contribution of each flux to the optimal solution. The input fluxes for lactate and for fatty acid reflect reduced costs of -17.5 and -129.00, respectively, when glucose is the energy source. Because this is the number of ATP that can be synthesized from either substrate, these fluxes represent an alternate option for ATP synthesis. Therefore, if the oxygen uptake flux increases by 3 units and the glucose uptake flux remains constant, the uptake of 1 unit of lactate will yield an additional 17.5 units of ATP.

                              
View this table:
[in this window]
[in a new window]
 
Table 3.   A comparison of ATP yields on different substrates

The malate-aspartate shuttle is active when glucose or lactate is the substrates for energy metabolism. When glycolysis is active, the glycerol-3-phosphate shuttle is a secondary option for transporting reducing equivalents into the mitochondrion from the cytosol. This exchanges one molecule of cytosolic NADH for one mitochondrial FADH2 and is therefore an inferior option for maximizing ATP production. The shadow price for NADH is 3.00 and that for FADH2 is 2.00. Thus ATP synthesis from glucose yields 36 ATP if the glycerol-3-phosphate shuttle is utilized. The shuttles are not required during fatty acid oxidation, because the reducing equivalents are entirely derived from the TCA cycle.

The comparison of oxygen and substrate intake can be characterized by using a phenotype phase plane (PhPP). The axes of the phase plot (Fig. 3A) represent the uptake fluxes of oxygen and substrate. When the optimal flux distributions for different combinations of oxygen and the substrate uptake are evaluated, regions of qualitatively different metabolic flux patterns are observed. These regions are delineated on the PhPP. The PhPPs for glucose and oxygen input fluxes show the division of the region into two regions (A and B). Region A represents complete metabolism of glucose via the TCA cycle because the oxygen input flux is in stoichiometric excess. In region B, the glucose input flux is in stoichiometric excess, and it represents the channeling of excess glucose toward lactate formation. The line of separation between the two regions is the line of optimality, where the two substrates are in stoichiometric balance.


View larger version (39K):
[in this window]
[in a new window]
 
Fig. 3.   A: the phase plot showing the relationship between oxygen and glucose fluxes. The shaded region (A) shown represents the region where oxygen is in excess and thus complete aerobic oxidation of glucose. The unshaded region (B) represents the combination of fluxes of glucose and oxygen where glucose is in excess and partial or complete (0 oxygen flux) anaerobic breakdown occurs. B: the phase plot for lactate and glucose utilization is shown. The shaded region (A) represents the combination of input fluxes of glucose and lactate where the qualitative flux distribution remains the same. Both lactate and glucose are utilized completely to form ATP. Because simulations calculate the phase plot over a range of input fluxes, region B represents the combination of fluxes where lactate is in excess. C: the phenotype phase plane (PhPP) for palmitate and glucose uptakes shows the interaction between the 2 substrates for energy metabolism. The oxygen input is constrained to 6 units. In region A, both palmitate and glucose are completely oxidized to form ATP, because the oxygen input flux (constrained at 6 units) is in stoichiometric excess. The boundary region represents the stoichiometric ratio at which the carbon consumption and oxygen consumption are completely balanced. In region B, the total carbon flux is in stoichiometric excess. Glucose is completely consumed, and the palmitate consumption corresponds to the appropriate flux on the boundary region.

The PhPPs for lactate and glucose uptakes (Fig. 3B) show the interaction between two substrates for energy metabolism. The oxygen input is constrained to 6 units. The results show that there are two qualitatively different metabolic flux distributions. In region A, both lactate and glucose are completely oxidized to form ATP, because the oxygen input flux is in stoichiometric excess. In region B, the total carbon flux is in stoichiometric excess. Glucose is completely consumed, and only part of the lactate input is directed toward energy metabolism.

The PhPPs for palmitate and glucose uptakes (Fig. 3C) show the interaction between the two substrates for energy metabolism. The oxygen input is constrained to 6 units. In region A, both palmitate and glucose are completely oxidized to form ATP, because the oxygen input flux (constrained at 6 units) is in stoichiometric excess. The boundary region represents the stoichiometric ratio at which the carbon consumption and oxygen consumption are completely balanced. In region B, the total carbon flux is in stoichiometric excess. Glucose is completely consumed, and the palmitate consumption corresponds to the appropriate flux on the boundary region.

In normal functioning hearts under aerobic conditions, free fatty acids (FFA) are primarily used for energy production. FBA predicts that under these circumstances, the ideal substrate to maximize ATP production is glucose. However, it is known that FFA metabolism inhibits pyruvate dehydrogenase (PDH) and phosphofructokinase (PFK) activities and reduces the flux of glucose catabolism. This metabolic regulation can be explicitly modeled in the FBA model, in addition to the stoichiometric requirements, by constraining the fluxes through these enzymes. Therefore, when PDH and PFK are inhibited, the primary source for energy is through FFA metabolism. From the phase plane (Fig. 3C), it is evident that if the glucose flux were constrained to be lower than the stoichiometric maximum of 1 unit for complete oxidation, the consumption of palmitate would increase. If the glucose flux were at 0.2, because of enzyme inhibition, the corresponding flux of palmitate for complete oxidation would be relatively close to the stoichiometric maximum and hence the energy derived from palmitate would be significantly higher than that from glucose.

Alterations in enzyme activity. The FBA model can be used to evaluate the systemic consequences of reduced enzyme activity. Mathematically, this is equivalent to setting a constraint on a specific flux (e.g., vAKGDH, the flux through alpha -ketoglutarate dehydrogenase) to be either 0 or some specific value, i.e., alpha i = beta i = 0 or a specific value. The effect of alterations in the activity of different enzymes involved in the metabolism of ATP production was investigated for the different substrates. Table 2 details the constraints imposed for the following scenarios.

Glucose as substrate. When citrate synthase, aconitase, or isocitrate dehydrogenase are inactive, the ATP yield drops significantly and the model predicts the accumulation of oxaloacetate. The predicted flux distribution for these cases is shown in Fig. 4A. The glucose input flux was maintained at 1 unit to compare with a normally functioning network. The carbon flux is partially cycled through phosphoenolpyruvate carboxylase and PDH. Figure 4B shows the dependence of ATP production on the activity of these enzymes. This dependence is linear, and the qualitative metabolic map is constant in this range of enzyme activities. That is, the active fluxes involved in energy metabolism remain the same until the enzyme activity is zero. An interesting feature that is predicted in this situation is the secretion of oxaloacetate, resulting from the excess carbon flux through glycolysis entering the TCA cycle.


View larger version (15K):
[in this window]
[in a new window]
 
Fig. 4.   A: the predicted optimal flux distribution when citrate synthase, aconitase, or isocitrate dehydrogenase is inactive. The PEP carboxylase is active, and the carbon flux is shunted through this pathway into the incomplete tricarboxylic acid cycle (TCA) cycle. B: the dependence of ATP production on the activity of citrate synthase, aconitase, or isocitrate dehydrogenase.

When any of the enzymes from alpha -ketoglutarate dehydrogenase to fumarase is inactive, the TCA cycle is again not completely functional. Excess carbon in this situation is secreted as alpha -ketoglutarate. An excess of alpha -ketoglutarate has in fact been identified as a feature of mitochondrial diseases involving mutations in the genes that code for these enzymes (12). The ATP yield is again significantly lower than the optimal value. The flux distribution for this scenario is shown in Fig. 5A. The dependence of ATP production on the activity of these enzymes is also linear and is shown in Fig. 5B.


View larger version (17K):
[in this window]
[in a new window]
 
Fig. 5.   A: the predicted optimal flux distribution when AKG dehydrogenase, succinyl-CoA synthetase, succinate dehydrogenase, and fumarase are inactive. The carbon flux is partially shunted through this pathway and partially through pyruvate dehydrogenase into the incomplete TCA cycle. B: the dependence of ATP production on the activity of AKG dehydrogenase, succinyl-CoA synthetase, succinate dehydrogenase, or fumarase.

The activity of the malate dehydrogenase enzyme also affects the ATP production in an adverse manner and predicts the accumulation of alpha -ketoglutarate. Interestingly, the effect of malate dehydrogenase is limited until the metabolic flux activity is reduced to <50% of the normal activity. The flux distributions predicted in this situation are different depending on the activity. Figure 6A shows the flux distribution when the malate dehydrogenase flux is 75% active. Here, the glycerol-3-phosphate shuttle is active and picks up the slack in malate-aspartate shuttle activity. Thus the reducing equivalents are shuttled by the glycerol-3-phosphate shuttle that exchanges a cytosolic NADH for a mitochondrial FADH2. This results in a lower stoichiometric yield of ATP. Malate dehydrogenase is active in the cytosol as well as the mitochondria (18), and a loss in the activity of either or both isozymes leads to the same predicted metabolic state. When the malate dehydrogenase activity is zero, the flux distribution (Fig. 6B) shows that the malate-aspartate shuttle functions in the reverse direction compared with the normal case and the ATP yield is decreased. The dependence of ATP production on enzyme activity is therefore biphasic (Fig. 6C), and the slope of the linear dependence is steeper when malate dehydrogenase activity is <50%.


View larger version (24K):
[in this window]
[in a new window]
 
Fig. 6.   A: the predicted optimal flux distribution when malate dehydrogenase (cytosolic and mitochondrial) is partially active. The glycerol-3-phosphate shuttle is activated and there is a decrease in ATP yield corresponding to the shuttling of FADH2 instead of NADH. B: the predicted optimal flux distribution when malate dehydrogenase is completely inactive. The malate-AKG cotransport does not occur, and the cytosolic and mitochondrial AKG pools do not interact. C: the dependence of ATP production on the activity of malate dehydrogenase. There are 2 qualitatively different phases illustrated by the 2 different slopes.

Table 4 lists the differences in ATP yield on oxygen carbon and carbon during energy metabolism from glucose. The ATP/carbon drops significantly, because most of the carbon is excreted out as organic acids, thus being incompletely metabolized.

                              
View this table:
[in this window]
[in a new window]
 
Table 4.   A comparison of ATP yields for the enzyme defect case

Palmitate and lactate as substrates. When palmitate or lactate is the substrate, there is no ATP generation when any of the enzymes are fully inactive. During ATP generation from glucose, the anaplerotic reactions that are part of the TCA cycle are able to reroute the metabolites and provide some reduced cofactors for the ETS. This is not possible when lactate or palmitate is the substrate, because these substrates are directly converted into TCA cycle intermediates.


    DISCUSSION
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

Mitochondrial function is a critical component of cellular physiology, and its metabolic objectives are constrained by stoichiometry. In this contribution, we have 1) constructed a model for mitochondrial energy metabolism under the stoichiometric constraints inherent in the biochemical reaction network, 2) applied the model to characterize the ATP production from three common substrates, and 3) simulated the effect of genetic changes on mitochondrial energy metabolism.

The network of biochemical reactions and metabolites that constitute the critical elements of mitochondrial function were identified. The steady-state mass balances on each metabolite were mathematically represented as an underdetermined matrix equation, S × v = 0, where the vector v represents the set of all feasible fluxes satisfying the stoichiometric constraints. To find a unique solution, it was postulated that the role of mitochondrial metabolism is the maximization of ATP production. Mathematically, the model is then formulated as an LP problem with this objective and the stoichiometric constraints. By adding additional constraints to the problem, which specify input fluxes for the different substrates, the LP can be solved to identify a unique flux distribution.

The model was applied to predict the maximal yield of ATP from the three common substrates of mitochondrial metabolism, glucose, lactate, and palmitate. The model predicts the expected ATP yields for each. The activities of the shuttles are also predicted in accordance with their accepted roles. The malate-aspartate shuttle is the preferred shuttle to transfer reducing equivalents across the mitochondrial membrane. The glycerol-3-phosphate shuttle is a secondary option, and it is also a less optimal route to synthesize ATP.

The model predictions are in qualitative agreement with key features of shuttle operation. Safer et al. (13) showed that amino-oxyacetate, which inhibits tranasaminase reactions, is less effective in inhibiting the transport of reducing equivalents into mitochondria with glucose than with lactate as substrate in rat myocardium. With the inhibition of the malate-aspartate shuttle, the glycerol-3-phosphate shuttle was found to allow a nearly equivalent rate of redox transfer across the mitochondria. Additionally, the flux of alpha -ketoglutarate to malate is equal to the rate of acetyl-CoA entry into the TCA cycle, which is also predicted by FBA (Fig. 2A).

The mitochondrial model was used to simulate the effect of changes in enzyme activities. Mutations in TCA cycle enzymes, particularly fumarase, alpha -ketoglutarate dehydrogenase, and succinate dehydrogenase have been implicated in mitochondrial diseases (3, 12). These diseases have been characterized by the accumulation of intermediates of the TCA cycle, primarily the presence of alpha -ketoglutarate and lactate. The result of loss of activity was investigated for all the TCA enzymes, and it was observed that when the enzymes in the series of reactions leading up to alpha -ketoglutarate were inactive, the ATP yield was very low and excess oxaloacetate was produced. When the enzymes in the second span of the TCA cycle, from alpha -ketoglutarate dehydrogenase to malate dehydrogenase, are inactive, the model predicts a low yield of ATP and the excess production of alpha -ketoglutarate. This prediction is in qualitative agreement with the data presented by Rustin et. al (12). They have observed that anomalous urinary excretion of specific organic acids, particularly alpha -ketoglutarate, occurs in patients with enzyme defects affecting the TCA cycle.

The pathogenesis of mitochondrial diseases can be very complex, owing to the variety of inheritance patterns and the diversity of metabolic states (23). Mutations can lead to different metabolic states and different mutations can produce similar metabolic states. Models can play a role in understanding and characterizing these complex relationships. The model presented here provides a platform on which more biochemical networks can be overlaid and other processes in which mitochondrial function plays a role can be studied.

The production of reactive oxygen species (ROS) as a result of defects in electron transport can lead to processes that activate the mitochondrial permeability transition pore and initiate a sequence of events that ultimately leads to cell apoptosis (7). The incorporation of the biochemical networks involved in this process in the formulation of an extended FBA model is being considered.

Wallace (23) states that the phenotypic expression of mtDNA diseases may involve two factors: the predisposing mutation and an age-related factor that causes a decline in mitochondrial function, which exacerbates the inherited defect. In terms of the biochemistry, this can be viewed as a change in the constraints of key fluxes over time, which can lead to qualitative and quantitative changes in the metabolic state. FBA has the necessary capabilities to address these issues. Flux analysis models themselves are not dynamic representations of metabolic functions. However, they provide the means to obtain several static snapshots of metabolism based on steady-state flux activity. Therefore, if a certain allosteric effect is known to effect an enzymatic reaction, the outcome of the severity of the effect on the metabolic function can be characterized. This can help complement a dynamic representation of the same effect that factors in substrates and product concentrations and their effect on the enzyme flux.

Perspectives

This paper presents a model based on a systems-analysis approach to characterize the integrated energy metabolism of mitochondria. The model is completed by postulating that the objective of mitochondrial energy metabolism is to maximize the rate of ATP production. Such a postulate requires several experimental studies to be verified. However, the predictions from the postulated model can be used to drive such experiments. The model can incorporate genomic information to predict the combined effect of genotype changes and environmental perturbations on the metabolic state. For a well-characterized genotype, under certain environmental conditions, it is possible to determine metabolic states that represent all possible physiological outcomes. This paper explored the metabolic states related to mitochondrial energy metabolism under certain conditions. The model accurately predicts different scenarios of mitochondrial growth on the standard substrates. The model also is in qualitative agreement with studies on the physiology of mitochondrial diseases, particularly in the accumulation of TCA cycle intermediates. Flux analysis of mitochondrial metabolism promises to be a useful methodology to understand and characterize the pathophysiology of mitochondrial diseases. Mitochondrial function is a critical part of the physiology of several organs. The genetic defects in several enzymes can contribute to physiological diseases that manifest at the organ level. Flux analysis can link genetics and physiology by representing individual fluxes that are related to the action of a gene or a set of genes in a model and characterizing the overall flux distribution that is representative of the physiological response. The work is currently being extended to study further disease-related questions, particularly those involved in ROS formation.


    ACKNOWLEDGEMENTS

Support for R. Ramakrishna through a grant from Procter and Gamble is gratefully acknowledged.


    FOOTNOTES

Present addresses: R. Ramakrishna, Physiome Sciences, Inc., 307, College Road East, Princeton, NJ 08540; J. S. Edwards, Dept. of Chemical Engineering, University of Delaware, Newark, DE 19716.

Address for reprint requests and other correspondence: R. Ramakrishna, Physiome Sciences, Inc., 307 College Rd. East, Princeton, NJ 08536.

The costs of publication of this article were defrayed in part by the payment of page charges. The article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

Received 17 May 1999; accepted in final form 3 October 2000.


    REFERENCES
TOP
ABSTRACT
INTRODUCTION
METHODS
RESULTS
DISCUSSION
REFERENCES

1.   Bonarius, HPJ, Hatzimanikatis V, Meesters KPH, De Gooijer CD, Schmid G, and Tramper J. Metabolic flux analysis of hybridoma cells in different culture media using mass balances. Biotech Bioengin 50: 299-318, 1996.

2.   Bonarius, HPJ, Schmid G, and Tramper J. Flux analysis of underdetermined metabolic networks: the quest for the missing constraints. Trends Biotech 15: 308-314, 1997.

3.   Bourgeron, T, Chretien D, Poggi-Bach J, Doonan S, Rabier D, Letouze P, Munnich A, Roetig A, Landrieu P, and Rustin P. Mutation of the fumarase gene in two siblings with progressive encephalopathy and fumarase deficiency. J Clin Invest 93: 2514-2518, 1994.

4.   Cairns, CB, Walther J, Harken AH, and Banerjee A. Mitochondrial oxidative phosphorylation thermodynamic efficiencies reflect physiological organ roles. Am J Physiol Regulatory Integrative Comp Physiol 274: R1376-R1383, 1998[Abstract/Free Full Text].

5.   DiMauro, S, and Bonilla E. Mitochondrial encephalomyopathies. In: The Molecular and Genetic basis of Neurological Disease, edited by Rosenberg RN, Prusiner SB, DiMauro S, and Barchi RL.. Boston, MA: Butterworth-Heinemann, 1997, p. 201-235.

6.   Fell, DA, and Small JA. Fat synthesis in adipose tissue. An examination of stoichiometric constraints. J Biochem 238: 781-786, 1986.

7.   Green, DR, and Reed JC. Mitochondria and apoptosis. Science 281: 1309-1311, 1998[Abstract/Free Full Text].

8.   Lehninger, AL, Nelson DL, and Cox MM. Principles of Biochemistry. New York: Worth, 1993.

9.   Ogata, H, Goto S, Fujibuchi W, and Kanehisa M. Computation with the KEGG pathway database. Biosystems 47: 119-128, 1998[ISI][Medline].

10.   Ogata, H, Goto S, Sato K, Fujibuchi W, Bono H, and Kanehisa M. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res 27: 29-34, 1999[Abstract/Free Full Text].

11.   Opie, LH. The Heart: Physiology and Metabolism. New York: Raven, 1991.

12.   Rustin, P, Bourgeron T, Parfait B, Chretien D, Munnich A, and Rotig A. Inborn errors of the Krebs cycle: a group of unusual mitochondrial diseases in human. Biochim Biophys Acta 1361: 185-197, 1997[Medline].

13.   Safer, B, Smith CM, and Williamson JR. Control of the transport of reducing equivalents across the mitochondrial membrane in perfused rat heart. J Mol Cell Cardiol 2: 111-124, 1971[ISI][Medline].

14.   Savinell, JM, and Palsson BO. Network analysis of intermediary metabolism using linear optimization. II. Interpretation of hybridoma cell metabolism. J Theor Biol 154: 455-473, 1992[Medline].

15.   Schapira, AHV Mitochondrial dysfunction in neurodegenerative disorders and aging. In: Mitochondrial Disorders in Neurology, edited by Scapira AHV, and DiMauro S.. Boston, MA: Butterworth-Heinemann, 1994.

16.   Schapira, AHV, and DiMauro S. Mitochondrial Disorders in Neurology. Boston, MA: Butterworth-Heinemann, 1994.

17.   Schon, EA. The mitochondrial genome. In: The Molecular and Genetic Basis of Neurological Disease, edited by Rosenberg RN, Prusiner SB, DiMauro S, and Barchi RL.. Boston, MA: Butterworth-Heinemann, 1997, p. 201-235.

18.   Small, WC, and McAlister-Henn L. Metabolic effects of altering redundant targeting signals for yeast mitochondrial malate dehydrogenase. Arch Biochem Biophys 344: 53-60, 1997[Medline].

19.   Tzagoloff, A. Mitochondria. New York: Plenum, 1982.

20.   Varma, A, and Palsson BO. Metabolic capabilities of Escherichia coli. I. Synthesis of biosynthetic precursors and cofactors. J Theor Biol 165: 477-502, 1993.

21.   Varma, A, and Palsson BO. Metabolic capabilities of Escherichia coli. II. Optimal growth patterns. J Theor Biol 165: 503-522, 1993.

22.   Varma, A, and Palsson BO. Metabolic flux balancing: basic concepts, scientific and practical use. Biotechnology 12: 994-998, 1994.

23.   Wallace, DC. Mitochondrial diseases in man and mouse. Science 283: 1482-1488, 1999[Abstract/Free Full Text].

24.   Whittaker, PA, and Danks SM. Mitochondria: structure, function and assembly. London: Longman, 1978.


Am J Physiol Regul Integr Comp Physiol 280(3):R695-R704
0363-6119/01 $5.00 Copyright © 2001 the American Physiological Society



This article has been cited by other articles:


Home page
BioinformaticsHome page
A. L. Knorr, R. Jain, and R. Srivastava
Bayesian-based selection of metabolic objective functions
Bioinformatics, February 1, 2007; 23(3): 351 - 357.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Cell Physiol.Home page
T. D. Vo and B. O. Palsson
Building the power house: recent advances in mitochondrial studies through proteomics and systems biology
Am J Physiol Cell Physiol, January 1, 2007; 292(1): C164 - C177.
[Abstract] [Full Text] [PDF]


Home page
Brief BioinformHome page
J. M. Lee, E. P. Gianchandani, and J. A. Papin
Flux balance analysis in the era of metabolomics
Brief Bioinform, June 1, 2006; 7(2): 140 - 150.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Physiol. Heart Circ. Physiol.Home page
L. Zhou, J. E. Salem, G. M. Saidel, W. C. Stanley, and M. E. Cabrera
Mechanistic model of cardiac energy metabolism predicts localization of glycolysis to cytosolic subdomain during ischemia
Am J Physiol Heart Circ Physiol, May 1, 2005; 288(5): H2400 - H2411.
[Abstract] [Full Text] [PDF]


Home page
Hum Mol GenetHome page
V. van Heyningen and P. L. Yeyati
Mechanisms of non-Mendelian inheritance in genetic disease
Hum. Mol. Genet., October 1, 2004; 13(suppl_2): R225 - R233.
[Abstract] [Full Text] [PDF]


Home page
J. Appl. Physiol.Home page
A. M. Payne, S. L. Dodd, and C. Leeuwenburgh
Life-long calorie restriction in Fischer 344 rats attenuates age-related loss in skeletal muscle-specific force and reduces extracellular space
J Appl Physiol, December 1, 2003; 95(6): 2554 - 2562.
[Abstract] [Full Text]


Home page
Am. J. Physiol. Regul. Integr. Comp. Physiol.Home page
B. Drew, S. Phaneuf, A. Dirks, C. Selman, R. Gredilla, A. Lezza, G. Barja, and C. Leeuwenburgh
Effects of aging and caloric restriction on mitochondrial energy production in gastrocnemius muscle and heart
Am J Physiol Regulatory Integrative Comp Physiol, February 1, 2003; 284(2): R474 - R480.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (28)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Ramakrishna, R.
Right arrow Articles by Palsson, B. O.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Ramakrishna, R.
Right arrow Articles by Palsson, B. O.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
Visit Other APS Journals Online