A mathematical model of iron import and trafficking in wild-type and Mrs3/4ΔΔ yeast cells

• Published in: BMC systems biology 2019-02, Vol.13 (1), p.23-23, Article 23
• Main Authors: Wofford, Joshua D, Lindahl, Paul A
• Format: Article
• Language: English
• Subjects: Baking yeast; Biochemistry; Biological Transport - genetics; Biosynthesis; Cation Transport Proteins - genetics; Cation Transport Proteins - metabolism; Cell growth; Clusters; Complexity; Computer simulation; Cytosol; Differential equations; Dilution; Heme; Imports; Iron; Iron - metabolism; Iron deficiency; Kinetics; Mathematical models; Matrix methods; Metabolism; Microscopy; Mitochondria; Mitochondria - metabolism; Mitochondrial Proteins - genetics; Mitochondrial Proteins - metabolism; Models, Biological; Mutation; Ordinary differential equations; Organelles; Oxygen - metabolism; Parameters; Phenotype; Phenotypes; Proteins; Saccharomyces cerevisiae; Saccharomyces cerevisiae - cytology; Saccharomyces cerevisiae - genetics; Saccharomyces cerevisiae - metabolism; Saccharomyces cerevisiae Proteins - genetics; Saccharomyces cerevisiae Proteins - metabolism; Sulfur; Vacuoles; Yeast
• ISSN: 1752-0509; 1752-0509
• DOI: 10.1186/s12918-019-0702-2

Iron plays crucial roles in the metabolism of eukaryotic cells. Much iron is trafficked into mitochondria where it is used for iron-sulfur cluster assembly and heme biosynthesis. A yeast strain in which Mrs3/4, the high-affinity iron importers on the mitochondrial inner membrane, are deleted exhibits a slow-growth phenotype when grown under iron-deficient conditions. However, these cells grow at WT rates under iron-sufficient conditions. The object of this study was to develop a mathematical model that could explain this recovery on the molecular level. A multi-tiered strategy was used to solve an ordinary-differential-equations-based mathematical model of iron import, trafficking, and regulation in growing Saccharomyces cerevisiae cells. At the simplest level of modeling, all iron in the cell was presumed to be a single species and the cell was considered to be a single homogeneous volume. Optimized parameters associated with the rate of iron import and the rate of dilution due to cell growth were determined. At the next level of complexity, the cell was divided into three regions, including cytosol, mitochondria, and vacuoles, each of which was presumed to contain a single form of iron. Optimized parameters associated with import into these regions were determined. At the final level of complexity, nine components were assumed within the same three cellular regions. Parameters obtained at simpler levels of complexity were used to help solve the more complex versions of the model; this was advantageous because the data used for solving the simpler model variants were more reliable and complete relative to those required for the more complex variants. The optimized full-complexity model simulated the observed phenotype of WT and Mrs3/4ΔΔ cells with acceptable fidelity, and the model exhibited some predictive power. The developed model highlights the importance of an Fe mitochondrial pool and the necessary exclusion of O in the mitochondrial matrix for eukaryotic iron-sulfur cluster metabolism. Similar multi-tiered strategies could be used for any micronutrient in which concentrations and metabolic forms have been determined in different organelles within a growing eukaryotic cell.