article: Iron is essential for all known life due to its redox properties; however, these same properties can also lead to its toxicity in overload through the production of reactive oxygen species. Robust systemic and cellular control are required to maintain safe levels of iron, and the liver seems to be where this regulation is mainly located. Iron misregulation is implicated in many diseases, and as our understanding of iron metabolism improves, the list of iron-related disorders grows. Recent developments have resulted in greater knowledge of the fate of iron in the body and have led to a detailed map of its metabolism; however, a quantitative understanding at the systems level of how its components interact to produce tight regulation remains elusive. A mechanistic computational model of human liver iron metabolism, which includes the core regulatory components, is presented here. It was constructed based on known mechanisms of regulation and on their kinetic properties, obtained from several publications. The model was then quantitatively validated by comparing its results with previously published physiological data, and it is able to reproduce multiple experimental findings. A time course simulation following an oral dose of iron was compared to a clinical time course study and the simulation was found to recreate the dynamics and time scale of the systems response to iron challenge. A disease state simulation of haemochromatosis was created by altering a single reaction parameter that mimics a human haemochromatosis gene (HFE) mutation. The simulation provides a quantitative understanding of the liver iron overload that arises in this disease. This model supports and supplements understanding of the role of the liver as an iron sensor and provides a framework for further modelling, including simulations to identify valuable drug targets and design of experiments to improve further our knowledge of this system.
  Iron is an essential element from archaea to complex eukaryotes and man [1], and is required for many processes including oxygen transport, DNA synthesis and respiration. Iron deficiency is the most common nutritional deficiency affecting a large proportion of all humans [2]. The redox activity which provides iron' s utility also means poorly regulated iron metabolism can lead to highly toxic free radicals [3]. Maintaining the delicate balance of iron requires robust cellular and systemic regulation since both iron deficiency and overload can cause cell death [4]. Recent research has lead to a much greater understanding of the mechanisms controlling iron metabolism and a global view of the interactions between iron-related components is beginning to emerge [5], [6]. The liver has been proposed to play a central role in the regulation of iron homeostasis [7] through the action of the recently discovered hormone hepcidin [8]. Hepcidin is expressed predominantly in the liver [9] and distributed in the serum to control systemic iron metabolism. Hepcidin acts on ferroportin to induce its degradation. Ferroportin is the sole iron exporting protein in mammalian cells [10], therefore hepcidin expression reduces iron export into the serum from enterocytes, and reduces iron export from the liver. Intracellular iron metabolism is controlled by the action of iron response proteins (IRPs) [11]. IRPs post-transcriptionally regulate mRNAs encoding proteins involved in iron metabolism. IRPs combined with ferritin and the transferrin receptors (TfR) make up the center of cellular iron regulation. Ferritin is the iron-storage protein forming a hollow shell which counters the toxic effects of free iron by storing iron atoms in a chemically less reactive ferrihydrite [12]. Extracellular iron circulates bound to transferrin (Tf), and is imported into the cell through the action of membrane bound proteins transferrin receptors 1 and 2 (TfR1 and TfR2). Human haemochromatosis protein (HFE) competes with transferrin-bound iron for binding to TfR1 and TfR2 [13]. Systems Biology provides an excellent methodology for elucidating understanding, through computational modelling, of the complex iron metabolic network. A quantitative model of iron metabolism allows for a careful and principled examination of the effect of the various components of the network. Modelling allows one to do “what-if” experiments leading to new hypotheses that can later be put to test experimentally. However, no comprehensive model of liver iron metabolism exists to date. Models have been published that cover specific molecular events only, such as the loading of iron in ferritin [14]. A qualitative map of mammalian iron metabolism provides a detailed overview of the molecular interactions involved in iron metabolism, including in specific cell types [6]. Similarly, a detailed model of iron metabolism and oxidative stress was described but uses a Boolean approach and is specific for yeast [15]. Quantitative models of the iron network have been recently described [16], [17], yet these include only a few components of the iron network. The model from Chifman et al. suggests that the dynamics of this iron network is stable [16]. Large-scale models of the metabolism of the hepatocyte [18], [19] and a generic human metabolism stoichiometric model [20] have also been published, but these contain only four reactions relating to iron metabolism. While they include iron transport, the receptors are not considered, and regulatory details are absent altogether. Existing models are therefore at two extremes of detail: very specific and very generic — but to address questions about hepatic iron regulation, what is desirable is a model that balances coverage and detail. This is the aim of the present work. One of the problems of modelling iron metabolism quantitatively and in detail arises from the lack of parameter values for many interactions. Recently, several of those parameters have been described in the literature (Table 1), particularly using technologies like surface plasmon resonance. This has enabled us to construct a detailed mechanistic kinetic model of human hepatocyte iron metabolism. The model has been validated by being able to reproduce data from several disease conditions — importantly, these physiological data were not used in constructing the model. This validation provides a sense of confidence that the model is indeed appropriate for understanding liver iron regulation and for predicting the response to various environmental perturbations.
  Initial validation of the model was performed by assessing the ability to recreate experimentally-observed steady-state concentrations of metabolites and rates of reactions. Simulations were run to steady state using the parameters and initial conditions from Tables 1 and 2. Table 3 compares steady-state concentrations of metabolites and reaction rates with experimental observations. Chua et al. [23] injected radio-labeled transferrin-bound iron into the serum of mice and measured the total uptake of the liver after 120 minutes. The uptake rate, when expressed as mol/s, was close to that found at steady state by the computational model (Table 3). A technical aspect of note in this steady-state solution, is that it is very stiff. This originates because one section of the model is orders or magnitude faster than the rest: the cycle composed of iron binding to ferritin, internalization and release. Arguably this could be resolved by simplifying the model, but the model was left intact because this cycling is an important aspect of iron metabolism and allows the representation of ferritin saturation. Even though the stiffness is high, our software is able to cope by using an appropriate numerical method. An oral dose of iron creates a fluctuation in serum transferrin saturation of approximately 10% [24]. The fixed serum iron concentration in the simulation was replaced by a transient increase in concentration equivalent to a 10% increase in transferrin saturation as a simulation of oral iron dosage on hepatocytes. The simulated hepcidin response (Figure 2) is consistent with the hepcidin response measured by Girelli et al. [24]. The time scale and dynamics of the hepcidin response to iron challenge has been accurately replicated in the simulation presented here. Although the exact dynamics of the simulated response is not validated by either experimental technique (mass spectrometry or ELISA) the simulation appears to present an approximation of the two experimental techniques reaching a peak between 4 and 8 hours and returning to around basal levels within 24 hours. The computational model supports the proposed role of HFE and TFR2 as sensors of systemic iron. Figure 3B shows that as the concentration of HFE bound to TfR2 (HFE-TfR2) increases with serum transferrin-bound iron (Tf-Fe_intercell), at the same time the abundance of HFE bound to TfR1 (HFE-TfR1) decreases. The increase in HFE-TfR2 complex, even though of small magnitude, promotes increased expression of hepcidin (Figure 3A). It is through this mechanism that liver cells sense serum iron levels and control whole body iron metabolism through the action of hepcidin. Although the labile iron pool increases with serum transferrin-bound iron in this simulation, this is only because the model does not include the action of hepcidin in reducing duodenal export of iron. Expression and secretion of hepcidin will have a global effect of reducing the labile iron pool. Hereditary haemochromatosis is the most common hereditary disorder with a prevalence higher than 1 in 500 [25]. Type 1 haemochromatosis is the most common and is caused by a mutation in the HFE gene leading to a misregulation of hepcidin and consequent systemic iron overload. A virtual HFE knockdown was performed by reducing 100-fold the rate constant for HFE synthesis in the model, to create a simulation of type 1 hereditary haemochromatosis. The simulation was run to steady state and results were compared with experimental findings. Qualitative validation showed the in silico HFE knockdown could reproduce multiple experimental findings as shown in Table 4. Quantitatively the model was unable to reproduce accurately the finding that Hfe −/− mice have 3 times higher hepatic iron levels [26]. This was due to the fixed intercellular transferrin-bound iron concentration in the model, unlike in Hfe −/− mice where there is an increase in transferrin saturation as a result of increased intestinal iron absorption [26]. Despite fixed extracellular conditions the model predicted an intracellular hepatocyte iron overload which would be further compounded by the systemic effects of the misregulation of hepcidin. The simulation recreated increased ferroportin levels despite the expression of ferroportin remaining the same as wild type which was consistent with mRNA measurements from Ludwiczek et al. [27]. mRNA based experiments can be used to validate expression rates and protein assays are able to validate steady state protein concentrations as both expression rates and steady state protein concentrations are available as results from the computational model. The model of haemochromatosis was also able to replicate the dynamics of experimental responses to changing dietary iron conditions. An approximate 2-fold increase in hepatic ferroportin expression is caused by increased dietary iron in both haemochromatosis and healthy mice [27]. The model presented here recreated this increase with increasing intercellular iron as can be seen in Figure 4. HFE knockout has been shown to impair the induction of hepcidin by iron in mouse [27] and human [28] hepatocytes and this was seen in the computational model as increasing transferrin-bound iron did not induce hepcidin as strongly as in HFE knockdown. Although an increase in transferrin receptor 2 was observed in the model (healthy; type 1 haemochromatosis), the up-regulation was slightly smaller than the change observed in vivo [29]. This is due to the model having fixed extracellular transferrin-bound iron concentration, in contrast to haemochromatosis where this concentration increases due to higher absorption in the intestine. Type 3 haemochromatosis results in similar phenotype as type 1 haemochromatosis, however the mutation is found in the TfR2 gene as opposed to HFE. A virtual TfR2 knockdown mutation was performed by decreasing 100-fold the rate constant of synthesis of TfR2 from the model. Model results were then compared with the findings of Chua et al. [23]. The simulation showed a steady state decrease of liver TfR1 from to with TfR2 knockdown. This is supported by an approximate halving of TfR1 levels in TfR2 mutant mice [23]. An increase in hepcidin and consequent decrease in ferroportin as seen in mice was matched by the simulation. An iron overload phenotype with increased intracellular iron is not recreated by the model of the TfR2 mutant. This is, again, due to the fixed serum transferrin-bound iron concentration, while in the whole body there would be increased iron absorption from the diet through the effect of hepcidin. It is known that the iron sensing by the transferrin receptors is responsive over a wide range of intercellular iron concentrations [33]. The present model reproduces this well (Figure 5,). Becker et al. argued that a linear response of a receptor to its signal over a wide range could be achieved through a combination of: high receptor abundance, increased expression when required, recycling to the surface of internalised receptors and high receptor turnover [34]. This was illustrated with the behaviour of the erythropoietin (EPO) receptor [34]. Since the present model contains essentially the same type of reactions that can lead to such a behaviour, simulations were carried out to investigate to what extent this linearity of response is present here. In this case it is the response of the total amount of all forms of TfR1 and TfR2 bound to Tf-Fe against the amount of Tf-Fe_intercell that is important. A variable was created in the model to reflect the total receptor response (see Materials and Methods), and this variable was followed in a time course response to an iron pulse (Figure 6). The response to the iron pulse is remarkably similar to the response of the EPO receptor to EPO [34]. Becker et al. [34] reported that the linearity of EPO-R response, i. e. the integral of the response curve, is increased by increasing turnover rate of the receptor and this property was also observed in the simulation of TfR1 response (Figure 5). The range in which the iron response is linear is smaller than that found for EPO (Figure 5). As TfR1' s half life in the model matches the experimentally determined value [35] the non-linear receptor response seen in the simulation is expected to be accurate. This suggests that TfR1 is a poor sensor for high levels of intercellular iron. On the other hand TfR2 is more abundant than TfR1 [35] and accordingly shows an increased linearity for a greater range of intercellular iron concentrations (Figure 7). This suggests the two transferrin receptors play different roles in sensing intercellular iron levels with TfR2 providing a wide range of sensing and TfR1 sensing smaller perturbations. The activation of TfR2 directly influences the expression of hepcidin and therefore it is desirable for it to sense large systemic imbalances. TfR1 does not modulate hepcidin expression itself instead it plays a primary role as an iron transporter.
  Iron is an essential element of life, in humans it is involved in oxygen transport, respiration, biosynthesis, detoxification, and other processes. Iron regulation is essential because iron deficiency results in debilitating anemia, while iron excess leads to free radical generation and is involved in many diseases [3]. It is clear that healthy life depends on tight regulation of iron in the body. The mechanisms involved in iron absortion, transport, storage and regulation form a complex biochemical network [6]. The liver has a central role in the regulation of systemic iron metabolism through secretion of the peptide hormone hepcidin. Here we analysed the hepatic biochemical network involved in iron sensing and regulation through a mathematical model and computer simulation. The model was constructed mostly based on in vitro biochemical data, such as protein complex dissociation constants. The model was then validated by comparison with experimental data from multiple physiological studies at both steady state and during dynamic responses. Where quantitative data were available the model matched these well and also qualitatively recreated many findings from clinical and experimental investigations. The simulation accurately modelled the highly prevalent iron disorder haemochromatosis. The disease state was simulated through altering a single parameter of the model and showed quantitatively how an iron overload phenotype occurs in patients with a HFE mutation. Due to the limited availability of quantitative clinical data on human iron metabolism, various other data sources, particularly from in vitro experiments and animal models, were integrated for the parameterisation of this model. This computational modelling effort constitutes a clinical translational approach, enabling data from multiple sources to improve our understanding of human iron metabolism. Several arguments could be raised to cast doubt on this approach, such as the the failure of in vitro conditions to mimic those in vivo or the difference between animal models and humans. This means that this type of data integration must be carefully monitored in terms of establishing the validity of the resulting model. Examining the behaviour of the model by simulating it at different values of initial conditions or other parameters (parameter scans) is important to establish the limits of utility of the model. Global sensitivity analysis is another approach that determines the boundaries of parameter variation that the model tolerates before it becomes too distant from the actual system behaviour. Model validation is an essential step in modelling. Validation should be carried out by enquiring if the model is able to match experimental observations that were not used to calibrate it. Here the model was validated by the simulation of haemochromatosis disease, where the model behaviour matched the clinical data (Table 4). The precise regulatory mechanism behind transferrin receptors and HFE controlling hepcidin expression remains to be validated experimentally, however the model presented here supports current understanding that the interaction of TfR2 and HFE form the signal transduction pathway that leads to the induction of hepcidin expression [36]. The global metabolic control analysis results support the identification of the transferrin receptors, particularly TfR2, and HFE as potential therapeutic targets; a result that is robust even to inaccuracies in parameter values. Although hepcidin would be an intuitive point of high control of this system (and therefore a good therapeutic target), in the present model this is not the case. It seems that targeting the promoters of hepcidin expression may be more desirable. However this conclusion has to be expressed with some reservation that stems from the fact that the global sensitivity analysis identified the hepcidin synthesis and degradation reactions in the group of those with the largest uncertainty. By changing parameter values by no more than 10% it would be possible to have the hepcidin expression and degradation show higher control. So it seems important that the expression of hepcidin be studied in more detail. We also predict that the control of hepcidin over the system would be higher if the model had included the regulation of intestinal ferroportin by hepatic ferroportin. The global sensitivity analysis, however, allows taking strong conclusions about the reactions for which the reference steady state is not much different from the maximal and minimal values. It turns out that these are the reactions that have the largest and the smallest control over the system variables. For example, the reactions with greatest control on the labile iron pool and iron export are those of the HFE-TfR2 system. But the reactions of the HFE-TfR1 system have always low control. These two conclusions are valid under a wide range of parameter values. Construction of this model required several assumptions to be made due to lack of measured parameter values, as described in Materials and Methods. These assumptions may or may not have a large impact on the model behaviour, and it is important to identify those that have a large impact, as their measurement will improve our knowledge the most. Of all the assumptions made, the rates of expression and degradation of ferroportin are those that have a significant impact on the labile iron pool in the model (see Table 5). This means that if the values assumed for these rate parameters were to be significantly different the model prediction for labile iron pool behaviour would also be different. The model is therefore also useful by suggesting experiments that will optimally improve our knowledge about this system. Limitations on the predictive power of the model occur due to the scope of the system chosen. Fixed serum iron conditions, which were used as boundary conditions in the model, do not successfully recreate the amplifying feedbacks that occur as a result of hepcidin expression controlling enterocyte iron export. To relieve this limitation, a more advanced model should include dietary iron uptake and the action of hepcidin on that process. The model predicts a quasi-linear response to increasing pulses of serum iron, similar to what has been predicted for the erythropoietin system [34]. Our simulations display response of the transferrin receptors to pulses of extracellular transferrin-bound iron that is similar to the EPO-R response to EPO (Figure 5). The integral of this response versus the iron sensed deviates very little from linearity in the range of physiological iron (Figure 6). Computational models are research tools whose function is to allow for reasoning in a complex nonlinear system. The present model can be useful in terms of predicting properties of the liver iron system. These predictions form hypotheses that lead to new experiments. Their outcome will undoubtedly improve our knowledge and will also either confirm the accuracy of the model or refute it (in which case it then needs to be corrected). The present model and its results identified a number of predictions about liver iron regulation that should be investigated further: The present model is the most detailed quantitative mechanistic model of cellular iron metabolism to date, allowing for a comprehensive description of its regulation. It can be used to elucidate the link from genotype to phenotype, as demonstrated here with hereditary haemochromatosis. The model provides the ability to investigate scenarios for which there are currently no experimental data available — thus making predictions about the system and aiding in experimental design.
  The model is constructed using ordinary differential equations (ODEs) to represent the rate of change of each chemical species. COPASI [37] was used as the software framework for model construction, simulation and analysis. Cell Designer [38] was used for construction of a SBGN process diagram (Figure 1). The model consists of 2 compartments representing the serum and the liver. Concentrations of haeme and transferrin-bound iron in the serum were fixed to represent constant extracellular conditions. Fixed metabolites simulate a constant influx of iron through the diet as any iron absorbed by the liver is effectively replenished. A labile iron pool (LIP) consumption reaction is added to represent various uses of iron and create a flow through the system. Some of the LIP consumption reaction would be attributed to heme biosynthesis however this process was not considered explicitly in this study to avoid unnecessary complexity and because the bone marrow is the major site of heme biosynthesis [39]. Initial concentrations for metabolites were set to appropriate concentrations based on a literature survey (Table 2). All metabolites formed through complex binding were set to zero initial concentrations (Table 2). The concentration of a chemical species at a time point in the simulation is determined by integrating the system of ODEs. For some proteins a half-life was available in the literature, but sources could not be found for their synthesis rates (translation). In this occurrence, estimated steady-state concentrations were used from the literature and a synthesis rate was chosen such that at steady state the concentration of the protein would be approximately accurate, following Equation 1: (1) This is solved for where is the steady-state concentration of the protein and is the degradation rate obtained from the half-life () using: (2) Complex formation reactions, such as binding of TfR1 to Tf-Fe for iron uptake, are modelled using the on and off binding constants as a forward and reverse mass action reaction. For example: (3) is modelled using two reactions: (4) (5) There is one ODE per each chemical species. The two reactions 4 and 5 add the following terms to the set of ODEs: (6) Intracellular haeme levels are controlled by a balance between uptake, export and oxygenation. haeme import through the action of haeme carrier protein 1 (HCP1), export by ATP-binding cassette sub-family G member 2 (ABCG2) and oxygenation by haeme oxygenase-1 (HO-1) follow Michaelis-Menten kinetics. HO-1 expression is promoted by haeme through by a Hill function (Equation (7) ). (7) (8) Where is the substrate, is the modifier, is the turnover number, is the ligand concentration which produces half occupancy of the binding sites of the enzyme, and is the Hill coefficient. Values of larger than produce positive cooperativity (i. e. a sigmoidal response); when the response is the same as Michaelis-Menten kinetics. A Hill coefficient of was assumed unless there is literature evidence for a different value. Where is not known it has been estimated to be of the order of magnitude of experimentally observed concentrations for the ligand. IRP/Iron-responsive elements (IRE) regulation is represented by Hill kinetics using Equation (7) to simulate the 3′ binding of IRP promoting the translation rate, and Equation (8) to represent the 5′ binding of IRP reducing the translation rate. Ferroportin degradation is modelled using 2 reactions: one representing the standard half-life and the other representing the hepcidin-induced degradation. A Hill equation (Equation 7) is used to simulate the hepcidin-induced degradation of ferroportin. Hepcidin expression is the only reaction modelled using a Hill coefficient greater than 1. Due to the small dynamic range of HFE-TfR2 concentrations a Hill coefficient of 5 was chosen to provide the sensitivity required to produce the expected range of hepcidin concentrations. The mechanism by which HFE-TfR2 interactions induce hepcidin expression is not well understood, but is thought to involve the mitogen-activated protein kinase (MAPK) signalling pathway [40]. The stimulus/response curve of the MAPK cascade has been found to be equivalent to a cooperative enzyme with a Hill coefficient of 4–5 [41], making the steep Hill function appropriate to model hepcidin expression. Ferritin modelling follows the work of Salgado et al. [14]. Iron from the LIP binds to, and is internalised in, ferritin with mass action kinetics. Internalised iron release from ferritin occurs through 2 reactions (intact ferritin release and release due to ferritin degradation). The average amount of iron internalised per ferritin affects the iron release rate and this is modelled using Equation 9 (adapted from [14]): (9) Where is internalised iron, is the rate constant and is the ratio of iron internalised in ferritin to total ferritin available. The value 0. 0048 was obtained by dividing the value given in Salgado et al. [14] by 50 as that simulation was scaled for 50 iron atom packages. Iron is also released from ferritin when the entire ferritin cage is degraded. The kinetics of ferritin degradation are mass action, however the amount of iron released when a ferritin cage is degraded is an average based on ferritin levels and total iron internalised in ferritin. Incorporating mass action and ferritin saturation ratio gives the following rate law for: (10) Iron export rate was modelled using Equation 7 with ferroportin as the modifier and a Hill coefficient of 1. was assumed to be around the steady state concentration of IRP (). A rate () of was used from Sarkar et al. [42] and these values were substituted into the equation and solved for. Ferroportin expression rates and degradation rates are poorly understood. Ferroportin abundance data [43] lead to an estimate of ferroportin concentration around. The hepcidin induced degradation of ferroportin is represented in the model by a rate law in the form of Equation 7 with a Hill coefficient (see above) and a equal to the measured concentration of hepcidin [44] (see Table 2). We then assume a maximal rate of degradation to be, and using the steady state concentration of ferroportin, the rate constant can be estimated as. The ferroportin synthesis rate was then calculated to produce the required steady-state concentration of ferroportin at the nominal hepcidin concentration. The HFE-TfR2 binding and dissociation constants were also not available and so it was assumed that they were the same as those of TfR1-HFE. Finally, the HFE-TfR and HFE-TfR2 degradation rates are also not known; we used a value that is an order of magnitude lower than the half life for unbound TfR (i. e. we assume the complex to be more stable than the free form of TfR). Although DMT1 may contribute towards transferrin-bound iron uptake in hepatocytes this contribution has been found to be minor and DMT1 knockout has little affect on iron metabolism [45], therefore DMT1 was not included in the model. The two iron response proteins (IRP1 and IRP2), which are responsible for cellular iron regulation, were modelled as a single pool in this study as the mechanistic differences in their regulatory roles are poorly understood. Equivalent regulation by both IRPs has been found in multiple studies [46]–[48]. Global sensitivity analysis was performed using the method proposed by Sahle et al. [32], where all parameter values were allowed to vary within of their nominal value in the model and we search for the maximum and minimum value that concentration- or flux-control coefficients of interest are able to reach within that parameter space. The searches were carried out with the particle swarm optimisation algorithm [49]. In order to process these optimisations in a reasonable time a HTCondor [50] distributed computing system was used, managed through the Condor-COPASI package [51]. To perform analysis of receptor response in a similar manner to the EPO system studied by Becker et al. [34] initial conditions were adjusted to recreate a similar virtual experiment. Haeme concentration was fixed at zero to isolate transferrin-bound iron uptake. The rate constant of the labile iron pool depletion reaction was reduced to balance the reduced iron uptake (which results in iron having a similar half-life to EPO in [34]). Initial concentrations for all metabolites were set to steady-state concentrations with the exception of the labile iron pool and iron bound to all receptors which were set to zero. Extracellular transferrin-bound iron was set at increasing concentrations to determine receptor response. Time courses were calculated for Tf-Fe-TfR1,2 (Tf-Fe) -TfR1, Tf-Fe-TfR2 and 2 (Tf-Fe) -TfR2 as iron binds its two receptors in a two-staged process. Two new variables were defined in COPASI which integrated the results of the time courses corresponding to the two receptors (in their different ligand states): (11) (12)
summary: Iron is an essential nutrient required for healthy life but, in excess, is the cause of debilitating and even fatal conditions. The most common genetic disorder in humans caused by a mutation, haemochromatosis, results in an iron overload in the liver. Indeed, the liver plays a central role in the regulation of iron. Recently, an increasing amount of detail has been discovered about molecules related to iron metabolism, but an understanding of how they work together and regulate iron levels (in healthy people) or fail to do it (in disease) is still missing. We present a mathematical model of the regulation of liver iron metabolism that provides explanations of its dynamics and allows further hypotheses to be formulated and later tested in experiments. Importantly, the model reproduces accurately the healthy liver iron homeostasis and simulates haemochromatosis, showing how the causative mutation leads to iron overload. We investigate how best to control iron regulation and identified reactions that can be targets of new medicines to treat iron overload. The model provides a virtual laboratory for investigating iron metabolism and improves understanding of the method by which the liver senses and controls iron levels.
title: A Computational Model of Liver Iron Metabolism