BACKGROUND
due to the involvement of the mitochondrion in many essential metabolic processes it is implicated in a wide and growing variety of disorders and pathologies  <cit> .
however, it is often difficult to elucidate how a disease phenotype relates to the underlying genetic cause, as these disorders have complicated and non-intuitive phenotypes due to the inherent complexity of mitochondrial metabolism.
a change in a small part of a complex network can affect the overall capabilities of a system, such as metabolite production, which can be impossible to predict without a systemic model.
models of metabolism allow these disorders to be studied by simulating the effects of a change in enzyme metabolism at a systems level and modelling these complex phenotypes as the consequence of changes in metabolic fluxes.
several models have been created that have attempted to represent the mitochondrion in its entirety  <cit> , although these reconstructions are far from complete.
their construction has been hindered by the lack of a defined mitochondrial proteome, as available localisation data is difficult to query and use in conjunction with metabolic pathway data.
therefore these models have been limited to well-established pathways of the mitochondrion.
these models are unlikely to be representative of the full metabolic capabilities of the mitochondrion, which may result in difficulties in modelling metabolic disorders that often have complicated and unintuitive phenotypes caused by complex interactions between reactions that may not otherwise be obviously interlinked.
evaluating the mitochondrial localization of thousands of reactions and providing a large list of additional reactions and metabolites is enabled by mitominer  <cit> , a mitochondrial proteomic database that integrates experimental proteomic localization data with annotation from public resources and metabolic pathway data.
the data and functionality of mitominer allows the creation of a more comprehensive metabolic model and provides a new opportunity to investigate mitochondrial metabolism.
here we describe the construction of a model of mitochondrial metabolism in human heart cells called ias <dig> by using the mitominer database and other public resources.
we demonstrate that by using flux balance analysis  normal mitochondrial metabolism can be successfully simulated.
the fba method has been described in detail elsewhere  <cit> , but briefly the main data required is a stoichiometric matrix , which represents all the reactions and metabolites in the model.
the matrix is used with the assumption that the model is at a quasi steady state, represented as s • ν =  <dig>  where ν is a flux distribution vector.
capacity and directionality constraints are then included by defining upper and lower bounds  of individual reactions.
directionality is incorporated by setting either the upper or lower bound of a reaction to zero.
although these constraints limit the number of possible metabolic flux distributions that satisfy the steady state, a very large number of possible solutions still exist, so a unique solution representing a flux distribution for the model is found by using an objective function.
the objective function is a set of criteria chosen to match the biological objective of the system in vivo and an optimal flux distribution is found that maximises  this function by using linear programming.
if the objective function of the model matches the biological objective of the system, the optimal flux distribution should correspond with the flux distribution in vivo.
for example a common objective function used for mitochondrion metabolism is maximum atp production, so to reflect the primary role of the mitochondrion in cellular energy production.
the benefits of the fba approach are that detailed kinetic information such as concentrations are not required, allowing the simulation of large models that would be infeasible using alternative methods.
this makes the technique particularly attractive in biomedicine, where a phenotype can be the result of many interacting reactions and detailed kinetic information for each enzyme involved under disease conditions is unavailable and difficult to obtain experimentally.
furthermore, the models can be used to discover and assess the effect of interventions on a disease state that may restore the system to a normal state and thus be of therapeutic value.
in our new model we account for the most recent stoichiometric information of the atp synthase  <cit> , thermodynamic constraints and reaction directionality, and the change in protonation state between the mitochondrial matrix and the cytosol due to the difference in ph between these compartments.
these are factors that have not all been accounted for in previous flux balance analysis simulations of mitochondrial metabolism.
having developed a model that simulates normal mitochondrial metabolism, we demonstrate the predictive capabilities of our new model, by investigating three tca cycle disorders that affect core mitochondrial metabolism.
the first disorder, fumarase deficiency , is caused by impairment of the fumarate hydratase enzyme, which converts fumarate into malate in the tca cycle.
the condition is exceptionally rare with around  <dig> cases reported globally to date  <cit> .
the effects of the disorder are acute and include developmental delay, severe mental retardation, language impairment, seizures and dysmorphic facial features  <cit> .
no medical treatments currently exist and most patients do not survive past early childhood  <cit> .
in most of these cases there is almost no residual enzyme activity  <cit> .
the diagnostic markers for this disorder are the presence of fumarate and 2-oxoglutarate in the urine, and lactate in the cerebrospinal fluid  <cit> .
the second disorder, succinate dehydrogenase deficiency  affects mitochondrial complex ii, which links the tca cycle with the electron transport chain by coupling the conversion of succinate to fumarate to the quinone pool.
in cases of isolated succinate dehydrogenase deficiency the phenotype is highly variable and can include leigh syndrome  <cit> , leukodystrophy  <cit> , kearns-sayre syndrome  <cit> , cardiomyopathy  <cit>  and mental and motor skill deterioration  <cit> .
the onset of symptoms can be from childhood to adulthood and appears to correlate with residual enzyme activity, with those with higher complex ii activity having milder symptoms and diagnosed later in life  <cit> .
the metabolic profile of the disorder is elevated levels of lactate in blood plasma  <cit> , and in some cases oxoglutarate and succinate is present in the urine  <cit>  and succinate accumulates in brain  <cit> .
the third disorder, α-ketoglutarate dehydrogenase deficiency  affects the conversion of oxoglutarate to succinyl-coa in the tca cycle.
the condition is extremely rare with only  <dig> cases recorded in the literature  <cit> .
onset of the disorder occurs shortly after birth and is characterised by encephalopathy and hyperlactatemia resulting in death in early childhood.
the metabolic profile of the disorder is the excretion of oxoglutarate  <cit>  and in many cases lactate  <cit> .
the α-ketoglutarate dehydrogenase complex consists of three protein subunits; e <dig>  e <dig> and e <dig>  with the latter subunit also a constituent of the pyruvate dehydrogenase and ketoacid dehydrogenase complexes.
therefore a defect with the e <dig> subunit can affect all three enzymes, resulting in a more complicated clinical phenotype.
we show that simulations using our model can replicate the disease state phenotypes of these disorders by correctly predicting the effects of these perturbations on normal mitochondrial metabolism.
the simulations show how alternative metabolic pathways can maintain production of atp, but lead to the over-production of other metabolites that must be excreted, as observed in vivo in the blood, urine and cerebrospinal fluid of patients.
furthermore, the simulations suggest that treatments, such as a low glucose diet or amino acid supplements, may reduce the risk of lactic acidosis and increase atp production and thus be of therapeutic value.
RESULTS
construction of the mitochondrial metabolic model
the metabolic model of reactions in a human heart mitochondrion, called ias <dig> , was constructed using the principle of metabolite availability, where a reaction was only included in the model if the substrates or products of the reaction were used elsewhere in the mitochondrial matrix .
in total  <dig> reactions with evidence of mitochondrial localisation were evaluated resulting in the inclusion of  <dig> reactions that formed a connected network using  <dig> metabolites in the mitochondrial matrix compartment.
an additional  <dig> cytosolic reactions were added to a separate compartment to represent reactions that were intrinsic to the function of the mitochondrion, such as glycolysis.
<dig> boundary conditions were included to represent import and export of metabolites to the system.
to more accurately reflect the in vivo capabilities of the mitochondrion and to remove all thermodynamically disallowed flux loops from the model,  <dig> of the  <dig> reactions and  <dig> of the  <dig> transport steps in the model were assigned a directionality constraint.
reaction directionality was set using thermodynamic information, general rules of irreversibility, and information from several public resources and the primary literature.
the final model contained no orphan metabolites and all internal flux loops were eliminated.
<dig> reactions were capable of having a flux under the current boundary conditions between the cytosol and the matrix, with the remaining  <dig> reactions requiring the import of acetaldehyde or propanoate to have a flux, which is currently disallowed in the model.
details of the reactions, metabolites and constraints included in the model can be found in additional files  <dig>   <dig> and  <dig>  the xml model of ias <dig>  has been deposited at the embl european bioinformatics institute's biomodels database .
simulation of normal mitochondrial physiology in cardiomyocytes
six separate pseudo reactions were devised to denote the metabolite requirements of the mitochondrion.
this included the production of  <dig> amino acids required for protein synthesis, rna and dna, haem, atp and lipids in a ratio found in the inner mitochondrial membrane of rat liver mitochondria .
separate simulations were run setting each of these pseudo reactions as the objective function  and the model was capable of fulfilling all of them.
for example, with simulations using maximum atp production as the objective function, core metabolism worked correctly with both glucose and lactate converted to pyruvate in the cytosol as appropriate for heart  <cit> , which was imported into the mitochondrion to take part in the tca cycle.
nadh was produced by the tca cycle in preference to fatty acid and ketone body oxidation although these processes were still active, and was used by the reactions representing the electron transport chain to move protons into the cytosolic compartment.
reactive oxygen species  were produced from complex i at a rate corresponding to  <dig> % of the electrons passed through it, as is consistent with recent experimental findings  <cit> .
atp synthase allowed these protons to move back into the mitochondrial matrix compartment, generating atp in the process.
the maximum atp production in this simulation was  <dig>  μmol/min/gdw, which is similar to the experimentally measured figure of atp hydrolysis of resting heart, which is approximately  <dig> μmol/min/gdw  <cit> .
both oxygen and glucose uptake were at the maximum allowable rate  and the largest fluxes occurred through the protein complexes of the electron transport chain.
when the results of this simulation are compared to experimentally determined figures, the tca cycle had a flux rate of approximately  <dig>  μmol/min/gdw, compared to the experimentally measured figure of  <dig>  μmol/min/gdw in rat  <cit>  while the fatty acid oxidation rate of  <dig>  μmol/min/gdw is similar to the experimentally determined rate in rat of  <dig>  μmol/min/gdw  <cit> .
very low levels of alanine, succinate, and glutamine were effluxed from the system, corresponding with observations of mammalian heart  <cit> .
in addition 61% of acetyl-coa production was derived from fatty acid degradation, 30% from glycolysis and the remaining 9% from ketone body degradation, which compares well with estimates in the literature that 60-90% of acetyl-coa production is from fatty acid degradation with the remaining 10-40% from glycolysis  <cit> .
simulation of mitochondrial disease: 1) fumarase deficiency
one of the most noticeable effects of fumarase deficiency is the presence of fumarate in the urine.
however, in the default model there was no mechanism for the efflux of fumarate from the system.
therefore, to model fumarase deficiency it was necessary to add an unconstrained transport step and boundary condition to allow fumarate in the mitochondrial matrix to leave the system.
separate simulations were run under each of the six objective functions with this adapted model and no flux was observed through the transport step or boundary condition under any simulation under normal conditions, showing that these additions did not affect the model.
to model fumarase deficiency the reaction that represents the fumarase enzyme in the model  was constrained to different flux values representing reduced enzyme activity and an objective function of maximum atp production used.
compared to simulations under normal conditions, the most notable effect of these constraints was a dramatic drop in maximum atp production of up to 96% .
fumarate was effluxed from the system when the fumarase flux was constrained to just 97% of its value under normal conditions, while lactate was effluxed at 20% of normal flux or below.
flux through the fatty acid oxidation pathway and ketone body degradation became inactive at 28% and 95% of normal fumarase flux respectively, due to inability of the acetyl-coa produced to enter the tca cycle, resulting from an inadequate supply of oxaloacetate that would normally be derived from fumarate.
to determine if the import of other metabolites could increase maximum atp production while fumarase was impaired, the restrictions on each constrained boundary condition were removed in separate simulations.
this identified the import of several amino acids, bicarbonate, and metabolites involved in the malate-aspartate shuttle as having a substantial effect on atp production .
in the initial simulations of fumarase deficiency, many of these metabolites were imported at the maximum rate allowed under normal boundary constraints, but these uptake rates had a relatively minor effect on the overall system.
proline, arginine and glutamate were imported and degraded to enter the tca cycle as oxoglutarate while oxaloacetate, aspartate and malate were imported via the malate-asparate shuttle to circumvent the gap in the tca cycle caused by the fumarase deficiency.
bicarbonate was used for a metabolic bypass to directly convert pyruvate to oxaloacetate, with the addition of atp.
the routes of these metabolites into central metabolism are shown in figure  <dig>
one of the clinical phenotypes of fumarase deficiency is the presence of oxoglutarate in the urine.
in the model, an efflux of oxoglutarate from the system was observed when the uptake rates of oxaloacetate, malate, citrate and isocitrate were increased.
the efflux was the result of the oxoglutarate, and oxodicarboxylate transporters requiring oxoglutarate to be counter exchanged for the import of the other metabolites into the matrix.
however, the uptake of citrate and isocitrate had a small effect on maximum atp production, increasing it by only 16% when an artificially large uptake rate of  <dig> μmol/min/gdw was allowed.
simulation of mitochondrial disease: 2) succinate dehydrogenase deficiency
the results from the simulations of succinate dehydrogenase deficiency were very similar to that of fumarase deficiency.
the main difference was that when the succinate dehydrogenase reaction  was constrained and an objective function of maximum atp production used, succinate was effluxed from the system rather than fumarate.
a large drop in atp production was also observed that was slightly more pronounced than fumarase deficiency , due to the lack of complex ii contributing to the quinone pool of the electron transport chain.
fatty acid oxidation and ketone body degradation became inactive and lactate was effluxed at very similar levels of residual enzyme flux to that of fumarase deficiency.
to represent that some patients have some residual enzyme activity, the import of metabolites that had an effect on atp production was investigated at 33% and 0% of flux of normal conditions.
regardless of the flux level chosen proline, arginine, valine, bicarbonate, aspartate, malate and oxaloacetate were found to have a positive effect on maximum atp production .
the entry routes of these metabolites into the tca cycle were the same as under fumarase deficiency .
simulation of mitochondrial disease: 3) α-ketoglutarate dehydrogenase deficiency
the effect of α-ketoglutarate dehydrogenase deficiency on maximum atp production was minor.
a reduction in atp production of just 5% was observed when flux through the reaction was constrained to zero and an objective function of maximum atp production used.
this was due to the presence of a bypass that converted oxoglutarate to succinate via succinate semialdehyde and 4-aminobutanoate  .
an efflux of oxoglutarate from the system was observed when the flux through the bypass was constrained by any level.
an efflux of lactate was not observed under any of these conditions.
to simulate the effect of a defect in the e <dig> subunit of α-ketoglutarate dehydrogenase that is shared with pyruvate dehydrogenase and ketoacid dehydrogenase enzyme complexes, the dihydrolipoamide dehydrogenase reaction was also constrained .
this resulted in the production of lactate, as pyruvate from glycolysis could no longer enter the tca cycle, but maximum atp production remained stable, as fatty acid synthesis increased proportionally  to maintain acetyl-coa production levels.
discussion
here we present a new, manually curated metabolic model of the mitochondrion present in human heart tissue, called ias <dig> .
it contains  <dig> mitochondrial matrix reactions,  <dig> cytosolic reactions,  <dig> compartment transport steps and  <dig> boundary conditions.
each reaction was manually evaluated in conjunction with the principle of metabolite availability before incorporation into the model, using  <dig> mitochondrial proteomic datasets from the mitominer database, combined with annotation from public resources and the literature.
the model includes directionality constraints based upon general rules of irreversibility, thermodynamics and information from public resources and the literature.
the resultant model contains no orphan metabolites and its directionality constraints allow every reaction in model to potentially possess a flux, while all flux loops have been eliminated.
it is encoded in sbml format  <cit> , and uses kegg  <cit>  identifiers where possible.
therefore, this model should be readily adaptable for use by other modellers either in isolation or with its addition to models of cellular metabolism.
we believe it is the most refined model of the mitochondrion currently available as demonstrated by simulations of normal conditions, which closely correspond with experimentally determined flux figures and metabolic observations of heart.
for example, fluxes though the tca cycle and fatty acid oxidation pathway while using the objective function of maximum atp production corresponded with measurements in rat heart .
the similarities between the model and the results of experiments in vivo suggest that the modelled flux distribution of core metabolism is biologically relevant.
the minor discrepancies between the flux figures may be due to imprecise boundary constraints, the biological objective in vivo for core metabolism not being described accurately by the model's objective function of maximum atp production, or difficulties in measuring these values experimentally.
these results are in contrast to those of previous mitochondrial fba models  <cit>  where a much higher uptake of oxygen  was permitted, while maximum atp production was 60% lower when this was set as the objective function.
this discrepancy may be due to the use of capacity constraints on the transport steps between the cytosol and mitochondrion in these models.
these constraints were derived from uptake rates determined experimentally in vitro, on isolated carrier proteins in liposomes, so are unlikely to reflect the uptake rates in vivo for the mitochondrion.
our model does not use capacity constraints on these transport steps and instead constrains the system at the boundary conditions.
the constraints placed on the boundary conditions of the system critically affected the behaviour of the model.
these constraints reflect the imports and exports of the system and were chosen carefully from the primary literature .
there is some uncertainty in these figures as often the only data available are not from human and in one case not from heart tissue.
in addition many of these experiments were conducted on isolated hearts, which may display significant metabolic differences to hearts in situ  <cit> .
therefore these figures may not correspond to the maximum possible uptake values, which are necessary to set the upper bound of the constraint.
previous reconstructions have tried to account for this experimental uncertainty by increasing these figures by an arbitrary amount, such as 25%  <cit> .
this is problematic because no data exists on what scale of increase is appropriate whereas the uptake of some compounds, such as oxygen, are a critical limiting factor to the system under most circumstances, and will usually be at the maximum allowable rate.
therefore in the absence of more relevant data it was decided to use experimental figures without modification.
further experiments to determine and verify the uptake rates used would be beneficial for further refinement of the model and in particular to aid simulations of perturbed states where uptakes can have a large impact on model behaviour.
in addition we included directionality constraints on the transport steps between the cytosol and the mitochondrial compartments to prevent flux loops, while reflecting the biological role of the underlying carrier and allowing normal metabolism.
however, it is possible that under disease conditions the accumulation of metabolites may affect the direction of transport.
in the metabolic disorders investigated here all metabolites that were experimentally reported as accumulating could be effluxed without the need to alter these constraints.
the only exception was fumarate, which is known to accumulate in fumarase deficiency, where an extra transport step was included to allow its efflux.
a common objective function used with models of microorganisms is the production of biomass.
biomass represents all the metabolites that are required for growth of that particular organism, with the assumption that the organism is biologically optimised for maximum growth.
the relative proportions of each metabolite in the biomass are derived from experimental figures measured from cellular biomass.
however, a single biomass objective function cannot currently be devised for the mitochondrion, as no experimental figures exist for its composition.
it is also unlikely that the biological objective in vivo is maximum growth as one of the main roles of the mitochondrion is to produce atp and so this would also have to be incorporated in the objective function to make it biologically relevant.
this would require the ratios between atp production and growth and maintenance of the mitochondrion to be experimentally verified and this is not currently available, and is likely to be extremely difficult if not impossible to determine.
therefore to represent the different aspects of mitochondrial metabolism, six separate pseudo reactions were devised, encompassing growth and maintenance as well as haem and atp production.
although maximum atp production is most likely to be the closest to the biological objective , the inclusion of the other pseudo reactions ensures the model can produce metabolites known to be required by the mitochondrion.
however the flux distributions found when using these other pseudo reactions individually as objective functions are unlikely to be realistic, as no one objective function in isolation will match the biological objective of the system in vivo.
the pseudo reactions that are used for these objective functions have been included in the model for convenience and to allow other researchers to easily replicate our results.
to demonstrate the predictive abilities of the model we simulated and analysed the complicated phenotypes of three disorders of the tca cycle using the ias <dig> model.
the mechanism of fumarase deficiency
patients with the most severe symptoms of the disorder often have residual fumarase activities near zero  <cit> .
this would equate to near zero flux through the fumarase reaction in the model.
under these conditions both lactate and fumarate were effluxed from the system, matching the phenotype.
initial simulations would predict that such a reduction would result in a very large decrease in maximum atp production in heart , which would seem likely to have catastrophic effects.
however, heart defects in patients have not been reported  <cit> , which would imply that the predicted reduction in atp production must be compensated.
in a series of simulations the boundary constraints were removed one by one to determine if increased metabolite uptake could counteract the decrease.
this identified the uptake of valine, aspartate, malate, proline, bicarbonate, oxoaloacetate, arginine, glutamine, serine and glucose as having a positive effect on atp production .
the effect of aspartate on the system supports the hypothesis of bourgeron et al.
<cit> , which suggests a reaction in the malate-aspartate shuttle is used to convert aspartate and oxoglutarate into oxaloacetate and glutamate, acting as a metabolic bypass and allowing the impaired part of the tca cycle to be circumvented.
they suggest the low level of aspartate found in the blood plasma of two patients in their study supports the existence of such a bypass.
the ias <dig> model predicts that if aspartate is available, this bypass is used.
a relatively small uptake resulted in a large difference in maximum atp production, which suggests it could be biologically relevant.
oxaloacetate and malate also had a large effect by using the reactions of the malate-aspartate shuttle.
in the simulations oxaloacetate was converted to malate before being transported to the matrix using the oxoglutarate transporter.
malate was then converted back into oxaloacetate to complete the tca cycle.
as a consequence of the transport step oxoglutarate was counter-exchanged from the matrix into the cytosol and then effluxed from the system.
such a mechanism in vivo may explain the oxoglutarate observed in urine of patients with fumarase deficiency.
the import of proline, arginine, glutamate and valine caused a large increase in atp production due to the degradation of these metabolites into oxoglutarate, increasing the flux through part of the tca cycle, and leading to an increase in fumarate efflux.
the model predicts that a corresponding increase in the efflux of glutamine and alanine should also be observed.
the production of glutamine was required to remove the ammonia produced during the deamination of these metabolites, whereas alanine efflux was involved in the glucose-alanine cycle, where pyruvate is converted into oxoglutarate .
however, both of these metabolites are further metabolised in the kidney and liver respectively, so it is not clear whether they would accumulate to a sufficient degree to be elevated in a patient metabolite profile.
the uptake of bicarbonate had a large impact on atp production by allowing the direct conversion of pyruvate to oxaloacetate.
however, due to the production of lactate, and bicarbonate acting as a cellular buffer, it may not be available at sufficient levels to have a significant affect in vivo.
serine had an effect on increasing atp production, but only from its degradation to glycine, producing nadh and nadph in the process.
it did not directly contribute to the tca cycle, which explains why even at high uptakes its effect on maximum atp production is lower than the other metabolites except glucose.
although increased glucose uptake did increase atp production the effect was small, while dramatically increasing lactate efflux.
therefore, although increasing glucose levels may be easily achieved physiologically, a low glucose diet, while increasing the levels of the other metabolites may be more beneficial to maintain atp production, while minimising lactic acidosis.
it is not apparent if in vivo the uptakes of any of these metabolites are increased beyond normal limits and have a significant impact on atp production.
although it may be unlikely that any one metabolite could sustain a high uptake flux under physiological conditions, it may be feasible that all of these metabolites contribute collectively to increasing atp production.
further experimental data would be required to test these predictions and to further refine the model.
a common occurrence throughout these simulations was fatty acid metabolism was restricted, as acetyl-coa could not enter the tca cycle.
if this occurs in vivo, these excess fatty acids may be stored in adipocytes or alternatively be expressed as an increase in the production of ketone bodies , a known cellular method to use excess acetyl-coa, rather than a reduction in fatty acid metabolism.
this would be supported by the elevated levels of acetoacetate in the blood plasma of two patients with the disorder  in comparison to controls , although levels of hydroxybutanoate appeared to be within the normal range  <cit> .
additional data from fumarase deficient patients are needed to verify if these predictions are correct.
the mechanism of succinate dehydrogenase deficiency
the simulations of succinate dehydrogenase deficiency were very similar to fumarase deficiency in both the reduction in atp production  and the effect of increasing certain metabolite uptake rates .
the metabolites that had a positive effect on maximum atp production were shared between the two deficiencies and used the same degradation pathways and entry points into the tca cycle .
the main difference was the efflux of succinate rather than fumarate as found in fumarase deficiency and a slightly lower level of atp production, as complex ii could no longer contribute to the electron transport chain.
the similarities between the disorders are unsurprising as the enzymes are adjacent in the tca cycle.
the main physiological difference between them is that many patients with succinate dehydrogenase deficiency have residual enzyme activity in the range of 10-50% whereas many patients with fumarase deficiency have 0% activity  <cit> .
however, it is not possible to predict what level of succinate dehydrogenase flux this will correspond with, as residual enzyme activity and enzyme flux are not directly correlated  <cit> .
however, even low levels of enzyme flux result in a dramatic increase in maximum atp production .
this has the effect that the uptake rates of metabolites that impact atp production can be much lower while returning production to normal levels , which may make these interventions more physiologically feasible.
the mechanism of α-ketoglutarate dehydrogenase deficiency
the effect of α-ketoglutarate dehydrogenase deficiency on maximum atp production was minor even when no flux was allowed through the reaction representing the enzyme.
this was the result of a bypass known as the gaba shunt  <cit> , which converts oxoglutarate to succinate using succinate semialdehyde and gaba as intermediates, circumventing the impaired reaction.
although the gaba shunt is usually associated with brain, there is evidence in mitominer, and the literature  <cit>  that it is present in heart.
oxoglutarate was only effluxed from the system if the flux through the gaba shunt was limited so that not all the flux could proceed through it.
if the bypass in vivo is unable to sustain the large flux required for the tca cycle to function at full capacity, this could explain the phenotype of oxoglutarate excretion.
gene expression data for patients with the disorder may give some indication whether the enzymes that constitute the gaba shunt are up-regulated in heart to allow the bypass to function.
gaba also acts as an inhibitory neurotransmitter and its involvement in the bypass may begin to explain the encephalopathy reported in patients.
an efflux of lactate was never observed when only α-ketoglutarate dehydrogenase was constrained.
therefore to determine whether a defect in subunit e <dig> of the enzyme complex, which is also shared with pyruvate dehydrogenase and ketoacid dehydrogenase complexes, was responsible, the reactions representing the latter were also restricted.
only when pyruvate dehydrogenase was constrained was lactate produced, as not all pyruvate from glycolysis could enter the tca cycle.
this suggests that a lactic acidosis phenotype is the result of a defect in the e <dig> subunit, and isolated α-ketoglutarate dehydrogenase deficiency can be discounted.
CONCLUSIONS
we have successfully simulated metabolism of heart mitochondria under normal conditions as demonstrated by close similarity of flux values with experimentally measured figures.
we show that the model also has the ability to simulate the metabolic phenotype of three rare single gene disorders and illustrate that the model can provide an insight into the disease mechanism.
as the mitochondrion is involved in many more complicated disorders including parkinson's, diabetes and heart failure, the next step is to use this model with further development to simulate and understand these conditions.
we believe that our model can be used as framework for integrating genomics, proteomics, metabolomics and clinical data and so contribute to the development of therapies for mitochondrial related diseases.

 due to the involvement of the mitochondrion in many essential metabolic processes it is implicated in a wide and growing variety of disorders and pathologies. however, it is often difficult to elucidate how a disease phenotype relates to the underlying genetic cause, as these disorders have complicated and non-intuitive phenotypes due to the inherent complexity of mitochondrial metabolism. a change in a small part of a complex network can affect the overall capabilities of a system, such as metabolite production, which can be impossible to predict without a systemic model. models of metabolism allow these disorders to be studied by simulating the effects of a change in enzyme metabolism at a systems level and modelling these complex phenotypes as the consequence of changes in metabolic fluxes. several models have been created that have attempted to represent the mitochondrion in its entirety, although these reconstructions are far from complete. their construction has been hindered by the lack of a defined mitochondrial proteome, as available localisation data is difficult to query and use in conjunction with metabolic pathway data. therefore these models have been limited to well-established pathways of the mitochondrion. these models are unlikely to be representative of the full metabolic capabilities of the mitochondrion, which may result in difficulties in modelling metabolic disorders that often have complicated and unintuitive phenotypes caused by complex interactions between reactions that may not otherwise be obviously interlinked. evaluating the mitochondrial localization of thousands of reactions and providing a large list of additional reactions and metabolites is enabled by mitominer, a mitochondrial proteomic database that integrates experimental proteomic localization data with annotation from public resources and metabolic pathway data. the data and functionality of mitominer allows the creation of a more comprehensive metabolic model and provides a new opportunity to investigate mitochondrial metabolism. here we describe the construction of a model of mitochondrial metabolism in human heart cells called ias by using the mitominer database and other public resources. we demonstrate that by using flux balance analysis normal mitochondrial metabolism can be successfully simulated. the fba method has been described in detail elsewhere, but briefly the main data required is a stoichiometric matrix, which represents all the reactions and metabolites in the model. the matrix is used with the assumption that the model is at a quasi steady state, represented as s • ν = where ν is a flux distribution vector. capacity and directionality constraints are then included by defining upper and lower bounds of individual reactions. directionality is incorporated by setting either the upper or lower bound of a reaction to zero. although these constraints limit the number of possible metabolic flux distributions that satisfy the steady state, a very large number of possible solutions still exist, so a unique solution representing a flux distribution for the model is found by using an objective function. the objective function is a set of criteria chosen to match the biological objective of the system in vivo and an optimal flux distribution is found that maximises this function by using linear programming. if the objective function of the model matches the biological objective of the system, the optimal flux distribution should correspond with the flux distribution in vivo. for example a common objective function used for mitochondrion metabolism is maximum atp production, so to reflect the primary role of the mitochondrion in cellular energy production. the benefits of the fba approach are that detailed kinetic information such as concentrations are not required, allowing the simulation of large models that would be infeasible using alternative methods. this makes the technique particularly attractive in biomedicine, where a phenotype can be the result of many interacting reactions and detailed kinetic information for each enzyme involved under disease conditions is unavailable and difficult to obtain experimentally. furthermore, the models can be used to discover and assess the effect of interventions on a disease state that may restore the system to a normal state and thus be of therapeutic value. in our new model we account for the most recent stoichiometric information of the atp synthase, thermodynamic constraints and reaction directionality, and the change in protonation state between the mitochondrial matrix and the cytosol due to the difference in ph between these compartments. these are factors that have not all been accounted for in previous flux balance analysis simulations of mitochondrial metabolism. having developed a model that simulates normal mitochondrial metabolism, we demonstrate the predictive capabilities of our new model, by investigating three tca cycle disorders that affect core mitochondrial metabolism. the first disorder, fumarase deficiency, is caused by impairment of the fumarate hydratase enzyme, which converts fumarate into malate in the tca cycle. the condition is exceptionally rare with around cases reported globally to date. the effects of the disorder are acute and include developmental delay, severe mental retardation, language impairment, seizures and dysmorphic facial features. no medical treatments currently exist and most patients do not survive past early childhood. in most of these cases there is almost no residual enzyme activity. the diagnostic markers for this disorder are the presence of fumarate and 2-oxoglutarate in the urine, and lactate in the cerebrospinal fluid. the second disorder, succinate dehydrogenase deficiency affects mitochondrial complex ii, which links the tca cycle with the electron transport chain by coupling the conversion of succinate to fumarate to the quinone pool. in cases of isolated succinate dehydrogenase deficiency the phenotype is highly variable and can include leigh syndrome, leukodystrophy, kearns-sayre syndrome, cardiomyopathy and mental and motor skill deterioration. the onset of symptoms can be from childhood to adulthood and appears to correlate with residual enzyme activity, with those with higher complex ii activity having milder symptoms and diagnosed later in life. the metabolic profile of the disorder is elevated levels of lactate in blood plasma, and in some cases oxoglutarate and succinate is present in the urine and succinate accumulates in brain. the third disorder, α-ketoglutarate dehydrogenase deficiency affects the conversion of oxoglutarate to succinyl-coa in the tca cycle. the condition is extremely rare with only cases recorded in the literature. onset of the disorder occurs shortly after birth and is characterised by encephalopathy and hyperlactatemia resulting in death in early childhood. the metabolic profile of the disorder is the excretion of oxoglutarate and in many cases lactate. the α-ketoglutarate dehydrogenase complex consists of three protein subunits; e e and e with the latter subunit also a constituent of the pyruvate dehydrogenase and ketoacid dehydrogenase complexes. therefore a defect with the e subunit can affect all three enzymes, resulting in a more complicated clinical phenotype. we show that simulations using our model can replicate the disease state phenotypes of these disorders by correctly predicting the effects of these perturbations on normal mitochondrial metabolism. the simulations show how alternative metabolic pathways can maintain production of atp, but lead to the over-production of other metabolites that must be excreted, as observed in vivo in the blood, urine and cerebrospinal fluid of patients. furthermore, the simulations suggest that treatments, such as a low glucose diet or amino acid supplements, may reduce the risk of lactic acidosis and increase atp production and thus be of therapeutic value.  construction of the mitochondrial metabolic model the metabolic model of reactions in a human heart mitochondrion, called ias, was constructed using the principle of metabolite availability, where a reaction was only included in the model if the substrates or products of the reaction were used elsewhere in the mitochondrial matrix. in total reactions with evidence of mitochondrial localisation were evaluated resulting in the inclusion of reactions that formed a connected network using metabolites in the mitochondrial matrix compartment. an additional cytosolic reactions were added to a separate compartment to represent reactions that were intrinsic to the function of the mitochondrion, such as glycolysis. boundary conditions were included to represent import and export of metabolites to the system. to more accurately reflect the in vivo capabilities of the mitochondrion and to remove all thermodynamically disallowed flux loops from the model, of the reactions and of the transport steps in the model were assigned a directionality constraint. reaction directionality was set using thermodynamic information, general rules of irreversibility, and information from several public resources and the primary literature. the final model contained no orphan metabolites and all internal flux loops were eliminated. reactions were capable of having a flux under the current boundary conditions between the cytosol and the matrix, with the remaining reactions requiring the import of acetaldehyde or propanoate to have a flux, which is currently disallowed in the model. details of the reactions, metabolites and constraints included in the model can be found in additional files and the xml model of ias has been deposited at the embl european bioinformatics institute's biomodels database. simulation of normal mitochondrial physiology in cardiomyocytes six separate pseudo reactions were devised to denote the metabolite requirements of the mitochondrion. this included the production of amino acids required for protein synthesis, rna and dna, haem, atp and lipids in a ratio found in the inner mitochondrial membrane of rat liver mitochondria. separate simulations were run setting each of these pseudo reactions as the objective function and the model was capable of fulfilling all of them. for example, with simulations using maximum atp production as the objective function, core metabolism worked correctly with both glucose and lactate converted to pyruvate in the cytosol as appropriate for heart, which was imported into the mitochondrion to take part in the tca cycle. nadh was produced by the tca cycle in preference to fatty acid and ketone body oxidation although these processes were still active, and was used by the reactions representing the electron transport chain to move protons into the cytosolic compartment. reactive oxygen species were produced from complex i at a rate corresponding to % of the electrons passed through it, as is consistent with recent experimental findings. atp synthase allowed these protons to move back into the mitochondrial matrix compartment, generating atp in the process. the maximum atp production in this simulation was μmol/min/gdw, which is similar to the experimentally measured figure of atp hydrolysis of resting heart, which is approximately μmol/min/gdw. both oxygen and glucose uptake were at the maximum allowable rate and the largest fluxes occurred through the protein complexes of the electron transport chain. when the results of this simulation are compared to experimentally determined figures, the tca cycle had a flux rate of approximately μmol/min/gdw, compared to the experimentally measured figure of μmol/min/gdw in rat while the fatty acid oxidation rate of μmol/min/gdw is similar to the experimentally determined rate in rat of μmol/min/gdw. very low levels of alanine, succinate, and glutamine were effluxed from the system, corresponding with observations of mammalian heart. in addition 61% of acetyl-coa production was derived from fatty acid degradation, 30% from glycolysis and the remaining 9% from ketone body degradation, which compares well with estimates in the literature that 60-90% of acetyl-coa production is from fatty acid degradation with the remaining 10-40% from glycolysis. simulation of mitochondrial disease: 1) fumarase deficiency one of the most noticeable effects of fumarase deficiency is the presence of fumarate in the urine. however, in the default model there was no mechanism for the efflux of fumarate from the system. therefore, to model fumarase deficiency it was necessary to add an unconstrained transport step and boundary condition to allow fumarate in the mitochondrial matrix to leave the system. separate simulations were run under each of the six objective functions with this adapted model and no flux was observed through the transport step or boundary condition under any simulation under normal conditions, showing that these additions did not affect the model. to model fumarase deficiency the reaction that represents the fumarase enzyme in the model was constrained to different flux values representing reduced enzyme activity and an objective function of maximum atp production used. compared to simulations under normal conditions, the most notable effect of these constraints was a dramatic drop in maximum atp production of up to 96%. fumarate was effluxed from the system when the fumarase flux was constrained to just 97% of its value under normal conditions, while lactate was effluxed at 20% of normal flux or below. flux through the fatty acid oxidation pathway and ketone body degradation became inactive at 28% and 95% of normal fumarase flux respectively, due to inability of the acetyl-coa produced to enter the tca cycle, resulting from an inadequate supply of oxaloacetate that would normally be derived from fumarate. to determine if the import of other metabolites could increase maximum atp production while fumarase was impaired, the restrictions on each constrained boundary condition were removed in separate simulations. this identified the import of several amino acids, bicarbonate, and metabolites involved in the malate-aspartate shuttle as having a substantial effect on atp production. in the initial simulations of fumarase deficiency, many of these metabolites were imported at the maximum rate allowed under normal boundary constraints, but these uptake rates had a relatively minor effect on the overall system. proline, arginine and glutamate were imported and degraded to enter the tca cycle as oxoglutarate while oxaloacetate, aspartate and malate were imported via the malate-asparate shuttle to circumvent the gap in the tca cycle caused by the fumarase deficiency. bicarbonate was used for a metabolic bypass to directly convert pyruvate to oxaloacetate, with the addition of atp. the routes of these metabolites into central metabolism are shown in figure one of the clinical phenotypes of fumarase deficiency is the presence of oxoglutarate in the urine. in the model, an efflux of oxoglutarate from the system was observed when the uptake rates of oxaloacetate, malate, citrate and isocitrate were increased. the efflux was the result of the oxoglutarate, and oxodicarboxylate transporters requiring oxoglutarate to be counter exchanged for the import of the other metabolites into the matrix. however, the uptake of citrate and isocitrate had a small effect on maximum atp production, increasing it by only 16% when an artificially large uptake rate of μmol/min/gdw was allowed. simulation of mitochondrial disease: 2) succinate dehydrogenase deficiency the results from the simulations of succinate dehydrogenase deficiency were very similar to that of fumarase deficiency. the main difference was that when the succinate dehydrogenase reaction was constrained and an objective function of maximum atp production used, succinate was effluxed from the system rather than fumarate. a large drop in atp production was also observed that was slightly more pronounced than fumarase deficiency, due to the lack of complex ii contributing to the quinone pool of the electron transport chain. fatty acid oxidation and ketone body degradation became inactive and lactate was effluxed at very similar levels of residual enzyme flux to that of fumarase deficiency. to represent that some patients have some residual enzyme activity, the import of metabolites that had an effect on atp production was investigated at 33% and 0% of flux of normal conditions. regardless of the flux level chosen proline, arginine, valine, bicarbonate, aspartate, malate and oxaloacetate were found to have a positive effect on maximum atp production. the entry routes of these metabolites into the tca cycle were the same as under fumarase deficiency. simulation of mitochondrial disease: 3) α-ketoglutarate dehydrogenase deficiency the effect of α-ketoglutarate dehydrogenase deficiency on maximum atp production was minor. a reduction in atp production of just 5% was observed when flux through the reaction was constrained to zero and an objective function of maximum atp production used. this was due to the presence of a bypass that converted oxoglutarate to succinate via succinate semialdehyde and 4-aminobutanoate. an efflux of oxoglutarate from the system was observed when the flux through the bypass was constrained by any level. an efflux of lactate was not observed under any of these conditions. to simulate the effect of a defect in the e subunit of α-ketoglutarate dehydrogenase that is shared with pyruvate dehydrogenase and ketoacid dehydrogenase enzyme complexes, the dihydrolipoamide dehydrogenase reaction was also constrained. this resulted in the production of lactate, as pyruvate from glycolysis could no longer enter the tca cycle, but maximum atp production remained stable, as fatty acid synthesis increased proportionally to maintain acetyl-coa production levels. discussion here we present a new, manually curated metabolic model of the mitochondrion present in human heart tissue, called ias. it contains mitochondrial matrix reactions, cytosolic reactions, compartment transport steps and boundary conditions. each reaction was manually evaluated in conjunction with the principle of metabolite availability before incorporation into the model, using mitochondrial proteomic datasets from the mitominer database, combined with annotation from public resources and the literature. the model includes directionality constraints based upon general rules of irreversibility, thermodynamics and information from public resources and the literature. the resultant model contains no orphan metabolites and its directionality constraints allow every reaction in model to potentially possess a flux, while all flux loops have been eliminated. it is encoded in sbml format, and uses kegg identifiers where possible. therefore, this model should be readily adaptable for use by other modellers either in isolation or with its addition to models of cellular metabolism. we believe it is the most refined model of the mitochondrion currently available as demonstrated by simulations of normal conditions, which closely correspond with experimentally determined flux figures and metabolic observations of heart. for example, fluxes though the tca cycle and fatty acid oxidation pathway while using the objective function of maximum atp production corresponded with measurements in rat heart. the similarities between the model and the results of experiments in vivo suggest that the modelled flux distribution of core metabolism is biologically relevant. the minor discrepancies between the flux figures may be due to imprecise boundary constraints, the biological objective in vivo for core metabolism not being described accurately by the model's objective function of maximum atp production, or difficulties in measuring these values experimentally. these results are in contrast to those of previous mitochondrial fba models where a much higher uptake of oxygen was permitted, while maximum atp production was 60% lower when this was set as the objective function. this discrepancy may be due to the use of capacity constraints on the transport steps between the cytosol and mitochondrion in these models. these constraints were derived from uptake rates determined experimentally in vitro, on isolated carrier proteins in liposomes, so are unlikely to reflect the uptake rates in vivo for the mitochondrion. our model does not use capacity constraints on these transport steps and instead constrains the system at the boundary conditions. the constraints placed on the boundary conditions of the system critically affected the behaviour of the model. these constraints reflect the imports and exports of the system and were chosen carefully from the primary literature. there is some uncertainty in these figures as often the only data available are not from human and in one case not from heart tissue. in addition many of these experiments were conducted on isolated hearts, which may display significant metabolic differences to hearts in situ. therefore these figures may not correspond to the maximum possible uptake values, which are necessary to set the upper bound of the constraint. previous reconstructions have tried to account for this experimental uncertainty by increasing these figures by an arbitrary amount, such as 25%. this is problematic because no data exists on what scale of increase is appropriate whereas the uptake of some compounds, such as oxygen, are a critical limiting factor to the system under most circumstances, and will usually be at the maximum allowable rate. therefore in the absence of more relevant data it was decided to use experimental figures without modification. further experiments to determine and verify the uptake rates used would be beneficial for further refinement of the model and in particular to aid simulations of perturbed states where uptakes can have a large impact on model behaviour. in addition we included directionality constraints on the transport steps between the cytosol and the mitochondrial compartments to prevent flux loops, while reflecting the biological role of the underlying carrier and allowing normal metabolism. however, it is possible that under disease conditions the accumulation of metabolites may affect the direction of transport. in the metabolic disorders investigated here all metabolites that were experimentally reported as accumulating could be effluxed without the need to alter these constraints. the only exception was fumarate, which is known to accumulate in fumarase deficiency, where an extra transport step was included to allow its efflux. a common objective function used with models of microorganisms is the production of biomass. biomass represents all the metabolites that are required for growth of that particular organism, with the assumption that the organism is biologically optimised for maximum growth. the relative proportions of each metabolite in the biomass are derived from experimental figures measured from cellular biomass. however, a single biomass objective function cannot currently be devised for the mitochondrion, as no experimental figures exist for its composition. it is also unlikely that the biological objective in vivo is maximum growth as one of the main roles of the mitochondrion is to produce atp and so this would also have to be incorporated in the objective function to make it biologically relevant. this would require the ratios between atp production and growth and maintenance of the mitochondrion to be experimentally verified and this is not currently available, and is likely to be extremely difficult if not impossible to determine. therefore to represent the different aspects of mitochondrial metabolism, six separate pseudo reactions were devised, encompassing growth and maintenance as well as haem and atp production. although maximum atp production is most likely to be the closest to the biological objective, the inclusion of the other pseudo reactions ensures the model can produce metabolites known to be required by the mitochondrion. however the flux distributions found when using these other pseudo reactions individually as objective functions are unlikely to be realistic, as no one objective function in isolation will match the biological objective of the system in vivo. the pseudo reactions that are used for these objective functions have been included in the model for convenience and to allow other researchers to easily replicate our results. to demonstrate the predictive abilities of the model we simulated and analysed the complicated phenotypes of three disorders of the tca cycle using the ias model. the mechanism of fumarase deficiency patients with the most severe symptoms of the disorder often have residual fumarase activities near zero. this would equate to near zero flux through the fumarase reaction in the model. under these conditions both lactate and fumarate were effluxed from the system, matching the phenotype. initial simulations would predict that such a reduction would result in a very large decrease in maximum atp production in heart, which would seem likely to have catastrophic effects. however, heart defects in patients have not been reported, which would imply that the predicted reduction in atp production must be compensated. in a series of simulations the boundary constraints were removed one by one to determine if increased metabolite uptake could counteract the decrease. this identified the uptake of valine, aspartate, malate, proline, bicarbonate, oxoaloacetate, arginine, glutamine, serine and glucose as having a positive effect on atp production. the effect of aspartate on the system supports the hypothesis of bourgeron et al., which suggests a reaction in the malate-aspartate shuttle is used to convert aspartate and oxoglutarate into oxaloacetate and glutamate, acting as a metabolic bypass and allowing the impaired part of the tca cycle to be circumvented. they suggest the low level of aspartate found in the blood plasma of two patients in their study supports the existence of such a bypass. the ias model predicts that if aspartate is available, this bypass is used. a relatively small uptake resulted in a large difference in maximum atp production, which suggests it could be biologically relevant. oxaloacetate and malate also had a large effect by using the reactions of the malate-aspartate shuttle. in the simulations oxaloacetate was converted to malate before being transported to the matrix using the oxoglutarate transporter. malate was then converted back into oxaloacetate to complete the tca cycle. as a consequence of the transport step oxoglutarate was counter-exchanged from the matrix into the cytosol and then effluxed from the system. such a mechanism in vivo may explain the oxoglutarate observed in urine of patients with fumarase deficiency. the import of proline, arginine, glutamate and valine caused a large increase in atp production due to the degradation of these metabolites into oxoglutarate, increasing the flux through part of the tca cycle, and leading to an increase in fumarate efflux. the model predicts that a corresponding increase in the efflux of glutamine and alanine should also be observed. the production of glutamine was required to remove the ammonia produced during the deamination of these metabolites, whereas alanine efflux was involved in the glucose-alanine cycle, where pyruvate is converted into oxoglutarate. however, both of these metabolites are further metabolised in the kidney and liver respectively, so it is not clear whether they would accumulate to a sufficient degree to be elevated in a patient metabolite profile. the uptake of bicarbonate had a large impact on atp production by allowing the direct conversion of pyruvate to oxaloacetate. however, due to the production of lactate, and bicarbonate acting as a cellular buffer, it may not be available at sufficient levels to have a significant affect in vivo. serine had an effect on increasing atp production, but only from its degradation to glycine, producing nadh and nadph in the process. it did not directly contribute to the tca cycle, which explains why even at high uptakes its effect on maximum atp production is lower than the other metabolites except glucose. although increased glucose uptake did increase atp production the effect was small, while dramatically increasing lactate efflux. therefore, although increasing glucose levels may be easily achieved physiologically, a low glucose diet, while increasing the levels of the other metabolites may be more beneficial to maintain atp production, while minimising lactic acidosis. it is not apparent if in vivo the uptakes of any of these metabolites are increased beyond normal limits and have a significant impact on atp production. although it may be unlikely that any one metabolite could sustain a high uptake flux under physiological conditions, it may be feasible that all of these metabolites contribute collectively to increasing atp production. further experimental data would be required to test these predictions and to further refine the model. a common occurrence throughout these simulations was fatty acid metabolism was restricted, as acetyl-coa could not enter the tca cycle. if this occurs in vivo, these excess fatty acids may be stored in adipocytes or alternatively be expressed as an increase in the production of ketone bodies, a known cellular method to use excess acetyl-coa, rather than a reduction in fatty acid metabolism. this would be supported by the elevated levels of acetoacetate in the blood plasma of two patients with the disorder in comparison to controls, although levels of hydroxybutanoate appeared to be within the normal range. additional data from fumarase deficient patients are needed to verify if these predictions are correct. the mechanism of succinate dehydrogenase deficiency the simulations of succinate dehydrogenase deficiency were very similar to fumarase deficiency in both the reduction in atp production and the effect of increasing certain metabolite uptake rates. the metabolites that had a positive effect on maximum atp production were shared between the two deficiencies and used the same degradation pathways and entry points into the tca cycle. the main difference was the efflux of succinate rather than fumarate as found in fumarase deficiency and a slightly lower level of atp production, as complex ii could no longer contribute to the electron transport chain. the similarities between the disorders are unsurprising as the enzymes are adjacent in the tca cycle. the main physiological difference between them is that many patients with succinate dehydrogenase deficiency have residual enzyme activity in the range of 10-50% whereas many patients with fumarase deficiency have 0% activity. however, it is not possible to predict what level of succinate dehydrogenase flux this will correspond with, as residual enzyme activity and enzyme flux are not directly correlated. however, even low levels of enzyme flux result in a dramatic increase in maximum atp production. this has the effect that the uptake rates of metabolites that impact atp production can be much lower while returning production to normal levels, which may make these interventions more physiologically feasible. the mechanism of α-ketoglutarate dehydrogenase deficiency the effect of α-ketoglutarate dehydrogenase deficiency on maximum atp production was minor even when no flux was allowed through the reaction representing the enzyme. this was the result of a bypass known as the gaba shunt, which converts oxoglutarate to succinate using succinate semialdehyde and gaba as intermediates, circumventing the impaired reaction. although the gaba shunt is usually associated with brain, there is evidence in mitominer, and the literature that it is present in heart. oxoglutarate was only effluxed from the system if the flux through the gaba shunt was limited so that not all the flux could proceed through it. if the bypass in vivo is unable to sustain the large flux required for the tca cycle to function at full capacity, this could explain the phenotype of oxoglutarate excretion. gene expression data for patients with the disorder may give some indication whether the enzymes that constitute the gaba shunt are up-regulated in heart to allow the bypass to function. gaba also acts as an inhibitory neurotransmitter and its involvement in the bypass may begin to explain the encephalopathy reported in patients. an efflux of lactate was never observed when only α-ketoglutarate dehydrogenase was constrained. therefore to determine whether a defect in subunit e of the enzyme complex, which is also shared with pyruvate dehydrogenase and ketoacid dehydrogenase complexes, was responsible, the reactions representing the latter were also restricted. only when pyruvate dehydrogenase was constrained was lactate produced, as not all pyruvate from glycolysis could enter the tca cycle. this suggests that a lactic acidosis phenotype is the result of a defect in the e subunit, and isolated α-ketoglutarate dehydrogenase deficiency can be discounted.  we have successfully simulated metabolism of heart mitochondria under normal conditions as demonstrated by close similarity of flux values with experimentally measured figures. we show that the model also has the ability to simulate the metabolic phenotype of three rare single gene disorders and illustrate that the model can provide an insight into the disease mechanism. as the mitochondrion is involved in many more complicated disorders including parkinson's, diabetes and heart failure, the next step is to use this model with further development to simulate and understand these conditions. we believe that our model can be used as framework for integrating genomics, proteomics, metabolomics and clinical data and so contribute to the development of therapies for mitochondrial related diseases. 
the mechanism or these disorders can be difficult to elucidate due to the inherent complexity of mitochondrial metabolism. the results of simulations using the new model corresponded closely with phenotypic data under normal conditions and provided insight into the complicated and unintuitive phenotypes of the three disorders, including the effect of interventions that may be of therapeutic benefit, such as low glucose diets or amino acid supplements. to demonstrate the predictive abilities of the model, flux balance analysis was used to calculate metabolite fluxes under normal conditions and to simulate three metabolic disorders that affect the tca cycle: fumarase deficiency, succinate dehydrogenase deficiency and α-ketoglutarate dehydrogenase deficiency. extensive curation and manual refinement was used to create a model called ias containing reactions, metabolites and transport steps across the inner mitochondrial membrane. mitochondria are a vital component of eukaryotic cells and their dysfunction is implicated in a large number of metabolic, degenerative and age-related human diseases. a new model of mitochondrial metabolism was built on the principle of metabolite availability using mitominer, a mitochondrial proteomics database, to evaluate the subcellular localisation of reactions that have evidence for mitochondrial localisation. to understand how mitochondrial metabolic dysfunction contributes to these diseases, a metabolic model of a human heart mitochondrion was created. the model offers the ability to investigate other mitochondrial disorders and can provide the framework for the integration of experimental data in future studies.
 mitochondria are a vital component of eukaryotic cells and their dysfunction is implicated in a large number of metabolic, degenerative and age-related human diseases. the mechanism or these disorders can be difficult to elucidate due to the inherent complexity of mitochondrial metabolism. to understand how mitochondrial metabolic dysfunction contributes to these diseases, a metabolic model of a human heart mitochondrion was created.  a new model of mitochondrial metabolism was built on the principle of metabolite availability using mitominer, a mitochondrial proteomics database, to evaluate the subcellular localisation of reactions that have evidence for mitochondrial localisation. extensive curation and manual refinement was used to create a model called ias containing reactions, metabolites and transport steps across the inner mitochondrial membrane. to demonstrate the predictive abilities of the model, flux balance analysis was used to calculate metabolite fluxes under normal conditions and to simulate three metabolic disorders that affect the tca cycle: fumarase deficiency, succinate dehydrogenase deficiency and α-ketoglutarate dehydrogenase deficiency.  the results of simulations using the new model corresponded closely with phenotypic data under normal conditions and provided insight into the complicated and unintuitive phenotypes of the three disorders, including the effect of interventions that may be of therapeutic benefit, such as low glucose diets or amino acid supplements. the model offers the ability to investigate other mitochondrial disorders and can provide the framework for the integration of experimental data in future studies. 