In silico experimentation with a model of hepatic mitochondrial folate metabolism

  • H Frederik Nijhout1Email author,

    Affiliated with

    • Michael C Reed2,

      Affiliated with

      • Shi-Ling Lam1,

        Affiliated with

        • Barry Shane3,

          Affiliated with

          • Jesse F GregoryIII4 and

            Affiliated with

            • Cornelia M Ulrich5

              Affiliated with

              Theoretical Biology and Medical Modelling20063:40

              DOI: 10.1186/1742-4682-3-40

              Received: 20 October 2006

              Accepted: 06 December 2006

              Published: 06 December 2006



              In eukaryotes, folate metabolism is compartmentalized and occurs in both the cytosol and the mitochondria. The function of this compartmentalization and the great changes that occur in the mitochondrial compartment during embryonic development and in rapidly growing cancer cells are gradually becoming understood, though many aspects remain puzzling and controversial.


              We explore the properties of cytosolic and mitochondrial folate metabolism by experimenting with a mathematical model of hepatic one-carbon metabolism. The model is based on known biochemical properties of mitochondrial and cytosolic enzymes. We use the model to study questions about the relative roles of the cytosolic and mitochondrial folate cycles posed in the experimental literature. We investigate: the control of the direction of the mitochondrial and cytosolic serine hydroxymethyltransferase (SHMT) reactions, the role of the mitochondrial bifunctional enzyme, the role of the glycine cleavage system, the effects of variations in serine and glycine inputs, and the effects of methionine and protein loading.


              The model reproduces many experimental findings and gives new insights into the underlying properties of mitochondrial folate metabolism. Particularly interesting is the remarkable stability of formate production in the mitochondria in the face of large changes in serine and glycine input. The model shows that in the presence of the bifunctional enzyme (as in embryonic tissues and cancer cells), the mitochondria primarily support cytosolic purine and pyrimidine synthesis via the export of formate, while in adult tissues the mitochondria produce serine for gluconeogenesis.


              Folate and one-carbon metabolism play a central role in cellular physiology because they are intimately involved in the control of purine, pyrimidine, and glutathione synthesis, as well as the methylation of DNA, histones and a host of other key cellular components. Deficiencies in folate metabolism have been associated with a wide range of diseases and pathologies such as anemia, spina bifida, cancer, cardiovascular disease, and neuropsychiatric disorders. Aberrant folate metabolism can be caused by polymorphisms in the genes for enzymes in the folate and methionine cycles, environmental factors that increase oxidative stress, and dietary deficiencies in B vitamins. Thus, this part of cell metabolism is a locus where genetic, environmental, and behavioral variables interact to affect many aspects of health and disease [112].

              It has been known for almost 50 years that some reactions of the folate cycle in eukaryotes (Figure 1) are duplicated in the cytosol and mitochondria [13], while other reactions such as purine and pyrimidine synthesis occur only in the cytosol, and the glycine cleavage system occurs only in the mitochondria [14, 15]. Two folate substrates, dihydrofolate (DHF) and 5-methyltetrahydrofolate (5mTHF), occur only in the cytosol. Moreover, some enzymes of mitochondrial folate metabolism are highly up-regulated in embryos and cancer cells and virtually inactive in normal adult cells [16, 17].
              Figure 1

              Diagram of the reactions modeled in the present paper. Pink rectangles represent variable metabolites and blue ellipses are enzymes. Full names corresponding to acronyms are given in Additional file 1.

              Because substrates, enzymes, and function differ between the mitochondrial and cytosolic compartments, many questions arise. What specific role does the mitochondrial folate cycle play in overall cell metabolism? What is the reason for the down-regulation of mitochondrial MTD and MTCH in adult tissues? Why does SHMT occur in both compartments? What is the role of the mitochondrial GDC reaction? How does the system accommodate changes in the input of serine and glycine? What happens during protein or methionine loading? What are the roles of folate-binding proteins in the cytosol and the mitochondria? These questions have been the subject of extensive experimental investigation and theoretical discussions [1626].

              We have developed a mathematical model for mitochondrial and cytosolic one-carbon metabolism. The model extends our earlier models of cytosolic methionine and folate metabolism [2730]. We use the model to conduct in silico experiments that address many of the above questions and compare the results to experimental observations. The model gives insights into the mechanisms underlying the experimental results and allows us to test various hypotheses that have been proposed in the literature.

              In the following section we give a brief overview of our model. Full details of the model and the full names of all abbreviations used in the text and in the model are given in additional file 1.

              Model overview

              Figure 1 shows the biochemical pathways in the hepatic cellular model used in this paper. Rectangular boxes represent the substrates that can vary in the model, and the ellipses contain the acronyms of the enzymes that catalyze specific reactions. Full names of the substrates and enzymes appear in the Additional File 1. Non-boxed substrates are taken to be constant. The model consists of 23 coupled differential equations for the time courses of the boxed substrates. The formulae for the velocities of the various reactions are taken, when possible, from the experimental literature. In some cases we adjusted the rate constants within experimental ranges so the concentrations of folates in the cytosol and mitochondria would be similar to those observed experimentally. In this model, MTD, MTCH, FTD and FTS are active in the mitochondria, and TS and DHFR are up-regulated in the cytosol, so the model represents a liver cell that is actively dividing.

              This model is an extension of the model used by Reed et al. [30]. Mitochondrial folate substrates and enzymes were added, and sarcosine and dimethylglycine are new variables (each is assumed to have the same concentrations in the cytosol and mitochondria). The extracellular serine, glycine, and methionine concentrations can be specified as functions of time. The cytosolic and mitochondrial serine and glycine concentrations can vary, as can the cytosolic methionine concentration. Appropriate kinetics are used for transport between the external, cytosolic, and mitochondrial compartments. The mitochondrial and cytosolic HCOOH concentrations are allowed to vary in the present model, and formate is transported between the mitochondrial and cytosolic compartments. In addition, AICAR is a variable in the model that couples the PGT and AICART reactions. Finally, we have added a serine sink that corresponds to the use of cytosolic serine in gluconeogenesis and the tricarboxylic acid cycle. Full details of the model are given in Additional file 1.

              Table 1 shows the steady-state concentrations and velocities in the model assuming that the extracellular glycine, serine, and methionine concentrations are held constant at 300 μM, 150 μM, and 30 μM respectively. For the reversible reactions in Figure 1, the positive directions are as follows: 5,10-CH2-THF to 5,10-CH = THF to 10f-THF for MTD and MTCH; THF to 10f-THF for FTS; serine to glycine for SHMT; SAH to homocysteine for SAAH. We refer to the values shown in Table 1 as "normal" throughout this paper. The computed distributions of folates in the cytosol and mitochondria are similar to those reported by Cook [15].
              Table 1

              Cytosolic and mitochondrial concentrations and velocities at steady-state.



              Concentrations (μM)

              Velocities (μM/h)

              Concentrations (μM)

              Velocities (μM/h)









              5,10-CH = THF




              5,10-CH = THF



































































































              In the model, we express the quantities of substrates by concentrations and we assume that the mitochondria occupy 1/4 of the cell volume. So the same number of molecules will have different concentrations in the two compartments. For example, we assume that the total "normal" cellular folate concentration is 20 μM. If the folate molecules are equally divided between the cytosol and the mitochondria [15], then the total folate concentrations will be 13.3 μM and 40 μM in the cytosol and the mitochondria, respectively. We assume that sarcosine and dimethylglycine diffuse freely and thus have the same concentrations in the cytosol and mitochondria. For transport between compartments, the rates (in μM/h) are the rates of change of concentration in the cytosol. If a substrate has concentration Cc in the cytosol and concentration Cm in the mitochondria, then the observed total cellular concentration (when the compartments are combined) will be (0.75)Cc + (0.25)Cm. Table 2 shows the normal total cellular concentrations of the folate metabolites and amino acids.
              Table 2

              Total cellular concentrations of folate metabolites and amino acids at steady-state


              Concentration (μM)



              5,10-CH = THF


















              Table 3 shows the normal rates of transport of serine, glycine, and HCOOH, between the compartments. For example, at normal steady-state, the cytosol receives 799 μM/h of serine from the extracellular medium and loses 30.4 μM/h to the mitochondria (Table 3). Of course, the cytosol also loses 99.7 μM/h through the CBS reaction, 8.19 μM/h through the cytosolic SHMT reaction, and 662 μM/h to gluconeogenesis (Tables 1 and 3).
              Table 3

              Transport rates at steady-state (μM/h).





              Extracellular to cytosol




              Mitochondria to cytosol




              To gluconeogenesis



              Limitations of the model

              The model we developed is intended specifically to study the interaction among the folate cycles in the cytosol and mitochondria. No mathematical model can include the complete biological complexity of a system. Instead, a model should contain sufficient detail to allow investigators to study the phenomenon of interest without omitting features that are likely to affect the behavior of the model dramatically. The model used in this paper, and described in detail in Additional File 1, is no exception. For instance, we do not include leucovorin as a cytosolic folate substrate, we do not include the polyamine pathway, and we use only two methyltransferase reactions between SAM and SAH. We also do not include the fact that some folate substrates can regulate gene expression of some folate enzymes. We do not include the allosteric binding of folates to folate enzymes, because the main effect of these reactions is to maintain reaction velocities in the face of severe folate deficiency (studied in Nijhout et al. [28]). Our model assumes that cytosolic amino acids are used only as substrates in the folate and methionine cycles, and as a source for gluconeogenesis. We do not include protein catabolism as a source, or protein synthesis as a sink, for amino acids. Likewise, we do not include the transport of homocysteine between the blood and the cellular compartment. All these features are important and interesting aspects of folate metabolism, and several are currently under investigation, but they do not bear directly on the role of mitochondrial folate metabolism studied in this paper.


              A. Variation in serine and glycine

              There is general agreement that one of the functions of the mitochondrial folate cycle is to provide 1-C units to the cytosol for purine and pyrimidine synthesis and the methylation reactions in rapidly dividing cells [1416]. The two primary sources of 1-C units are the mSHMT reaction, which converts serine to glycine, and the GDC reaction, which breaks down glycine, producing 5,10-methyleneTHF. Additional contributions are made by the SDH and the DMGD reactions. All four of these reactions use THF to produce 5,10-methyleneTHF, which is converted to 10f-THF. The FTS reaction regenerates THF and produces free formate that is exported to the cytosol. There are a number of natural questions here. (1) How does this system operate at different levels of extracellular serine and glycine? (2) How sensitive is the production of formate to the balance between serine and glycine? (3) How sensitive is the production of purines and pyrimidines in the cytosol to the supply of formate?

              To investigate these questions, we systematically altered the external concentrations of serine and glycine from their normal values of 150 μM and 300 μM, respectively. Figure 2A shows the cytosolic and mitochondrial glycine and serine concentrations as the extracellular glycine concentration is varied from 100 μM to 1000 μM. As external glycine increases, both the cytosolic and mitochondrial glycine concentrations increase, but not as dramatically, because the reverse transport out of the mitochondria and cytosol increases as the concentrations rise. The cytosolic serine concentrations also rise because of the interconversion of glycine and serine by the SHMT reactions.
              Figure 2

              A: Response of selected model variables to variation in external glycine input. B: Response of model variables to variation in external serine input.

              As external glycine rises, both the cytosolic and the mitochondrial SHMT reactions reverse and run in the glycine to serine direction (Fig. 2A); the mitochondrial reaction is more sensitive to external glycine. This reversal was observed by Kastanos, Woldman, and Appling [21], who grew yeast in pure serine and pure glycine environments. In spite of great variations in the glycine and serine concentrations (Fig. 2A, middle panel, blue curve) and the reversals in both SHMT reactions, the rate at which formate is supplied to the cytosol remains remarkably constant (Fig. 2A). In addition, the rates of the TS and AICART reactions (for thymidylate and purine synthesis, respectively) change very little despite large changes in external glycine. The metabolites of the methionine cycle also change relatively little except for SAM, which rises at low external glycine because the GNMT reaction rate declines (Fig. 2A).

              Figure 2B shows the cytosolic and mitochondrial serine and glycine concentrations as the extracellular serine concentration is varied from 50 μM to 400 μM. As external serine increases, both the cytosolic and mitochondrial glycine concentrations increase, and there is an even greater increase in the cytosolic and mitochondrial serine concentrations. As external serine decreases from normal, both the cytosolic and the mitochondrial SHMT reactions reverse and run in the glycine to serine direction (Fig. 2B); as before, the mitochondrial reaction is more sensitive. The rates of the TS and AICART reactions and the rate of transport of formate out of the mitochondria are again remarkably stable.

              The methionine cycle metabolites are more sensitive to external serine than to external glycine, especially at low external serine. This is because serine is required for the CBS reaction, which slows down and causes methionine cycle metabolites to accumulate, particularly as SAM, owing the internal regulatory mechanisms of the methionine cycle [27, 29].

              Lewis et al. [18] radiolabeled the methyl group of SAM. Very little of the radioactivity appeared in metabolites of the folate and methionine cycles because most of these radiolabeled methyl groups were transferred to other substrates by the methylation reactions. However, when glycine was elevated, the amount of radiolabel in HCOOH, serine, and CO2 went up considerably. This can be easily explained by the model. More glycine results in acceleration of the GNMT reaction so more of the radiolabeled methyl groups are put into sarcosine. Then, in the mitochondria, the sarcosine becomes either radiolabeled 5,10-CH2-THF or CO2. The radiolabeled 5,10-CH2-THF makes radiolabeled HCOOH via the MTD, MTCH, and FTS reactions and radiolabeled CO2 via the FTD reaction. High glycine increases the glycine-to-serine rate of the SHMT reactions in both cytosol and the mitochondria, so more radiolabel appears in serine.

              B. Reduced folate status

              Figure 3 shows the percentage change in the steady-state concentrations or rates of various metabolites and reactions in the presence of a 50% reduction in folate levels. The concentration of 5mTHF drops, of course, and this releases the inhibition of GNMT so the GNMT reactions goes faster and depletes SAM. Methionine and SAM are also depleted because less homocysteine is remethylated. Because of the long-range regulations in the methionine cycle [29], the rate of the DNMT methylation reaction decreases only modestly. The transport of formate from the mitochondria decreases and thymidylate synthesis and purine synthesis rates are reduced dramatically. The rate of the GDC reaction in the mitochondria drops because mitochondrial THF is much lower. Since the GDC reaction catabolizes glycine, the concentration of glycine rises, which also drives up the concentration of serine via the SHMT reaction. It was observed by Allen et al. [31] that both sarcosine and dimethylglycine are elevated in patients with folate deficiency, clinical results that we observe also in the model (Figure 3). More sarcosine is produced in the cytosol because the GNMT rate is elevated (because there is less inhibition by 5mTHF), and more dimethylglycine is produced because the rate of the BHMT reaction is elevated. Since the concentration of THF is lower, sarcosine and dimethylglycine are used at lower rates in the mitochondria.
              Figure 3

              Change in the steady-state concentrations of selected metabolites and steady-state reaction velocities of selected reactions in the presence of a 50% reduction in folate levels.

              C. SHMT expression

              MacFarlane et al. [32] report that mice lacking liver cytosolic SHMT have a normal SAM/SAH ratio. To test this in the model, we eliminated the cytosolic SHMT reaction entirely. This resulted in a steady-state of [SAM] = 62.6 μM, [SAH] = 13.1 μM. Normal values of these metabolites are in the model [SAM] = 63.6 μM, [SAH] = 13.1 μM. This is not surprising because the normal net rate of the cSHMT reaction is very low (9.2 μM/h), as this enzyme is poised to modulate fluctuations in serine and glycine by interconverting them.

              Herbig et al. [14] found that increasing the expression of cSHMT lowers the SAM concentration in a glycine-dependent manner. They suggested two possible alternative mechanisms: (1) that cSHMT is in competition with MTHFR for 5,10-CH2-THF, and therefore higher SHMT expression should lower [5mTHF] and thus lower the rate of remethylation of homocysteine to methionine, lowering SAM; or (2) that cSHMT sequesters 5mTHF, lowering its free concentration and thereby lowering SAM as in (1). In the current model, if SHMT does not bind to 5mTHF, we found that the SAM concentration is quite insensitive to the amount of cSHMT and also insensitive to the external glycine concentration, which does not support the first hypothesis. When we add the binding of 5mTHF to SHMT to the model, and then up-regulate cSHMT, we find that the SAM concentration is substantially reduced, supporting the second hypothesis (simulations not shown).

              Several authors have observed that Chinese hamster ovary (CHO) cells that lack mSHMT are glycine auxotrophs (see, for example, [32, 33]). These findings have been interpreted as indicating that mitochondria normally supply glycine to the cytosol. The model suggests that the situation is more complicated and more interesting. In the model, under normal conditions in hepatic cells, there is a net influx of glycine into the mitochondria (Table 4). If we eliminate the mSHMT reaction, there is a modest decrease in cytosolic and mitochondrial glycine, but the net flux remains from cytosol to mitochondria. On the other, if mSHMT is normal and we set the external glycine to zero, the concentrations of glycine in the cytosol and mitochondria decline to about one quarter of their normal value but not to zero. In addition, the mitochondria become net exporters of glycine to the cytosol (Table 4). This is because the velocity of mSHMT reaction from serine to glycine increases 10-fold. However, if mSHMT is eliminated and the external glycine is set to zero, then the glycine concentration declines to almost one tenth its normal value and formate production by the mitochondria is cut in half. Thus the elimination of mSHMT and external glycine together is clearly very detrimental to the cell. Interestingly, the net flux of glycine remains from cytosol to mitochondria. If one adds glycine back into the external medium the cell becomes almost normal (Table 4).
              Table 4

              Effects of variation in mSHMT and external glycine*



              mSHMT = 0

              Extenal gly = 0

              mSHMT = 0 Extenal gly = 0

              cellular glycine





              flux of glycine from mitochondria to cytosol





              velocity of mSHMT





              velocity of cSHMT





              flux of HCOOH to cytosol





              *Concentrations are in μM and fluxes and velocities are in μM/h.

              D. The GDC reaction

              In order to study the contribution of the GDC reaction, we set its velocity to zero in the model and calculated the percentage change in concentrations and fluxes at the new steady state (Figure 4). Because the GDC reaction is turned off, the concentration of 5,10-CH2-THF in the mitochondria drops dramatically, which lowers the flux through the MTD, MTCH, and FTS reactions. Thus, much less formate is produced in the mitochondria and therefore the rate of export of formate to the cytosol declines to about 50% of its former value. Because of the reduced supply of formate, the concentration of cytosolic 10f-THF goes down. This has two effects. First, fewer purines are produced and second, the net flux from 10f-THF to 5,10-CH2-THF reverses so that the net flux is from 5,10-CH2-THF to 10f-THF. This reduces the cytosolic concentration of 5,10-CH2-THF, which causes thymidine synthesis to drop. It also makes the concentration of 5mTHF drop, which causes the homocysteine concentration to rise. The decline in [5mTHF] releases the inhibition of GNMT, which draws down [SAM].
              Figure 4

              Change in the steady-state concentrations of selected metabolites and steady-state velocities of selected reactions when the GDC reaction is eliminated.

              E. Eliminating the mitochondrial bifunctional enzyme

              The mitochondrial bifunctional enzyme is active during embryonic development and in transformed cells [17, 24], but not in the adult liver. We examined the effect of eliminating the bifunctional enzyme by setting the Vmax of the mitochondrial MTD and MTCH reactions to zero. The significant changes in mitochondrial and cytosolic one-carbon metabolism are shown in Figure 5. The mitochondrial GDC reaction slows down somewhat, and the mSHMT reverses dramatically in the direction of serine production. Formate production by the mitochondria is reduced to zero. Thus the mitochondria have switched from being formate factories to being serine factories.
              Figure 5

              Effect of eliminating the mitochondrial bifunctional enzyme. A: selected reaction velocities when bifunctional enzyme is active. B: reaction velocities when bifunctional enzyme is eliminated. Reaction velocities are indicated with red numbers, and the units are μM/h.

              In the cytosol, thymidylate and purine synthesis are markedly reduced, and the cMTD and cMTCH reactions reverse and now run strongly in the direction from 5,10-CH2-THF to 10f-THF. The export of serine to gluconeogenesis is greater than the serine import into the cell from the blood (the opposite is true when the bifunctional enzyme is present). Thus, when the bifunctional enzyme is present (the case we call normal), mitochondrial folate metabolism produces formate for the cytosol for purine and pyrimidine synthesis and methylation reactions. When the bifunctional enzyme is absent, mitochondrial folate metabolism produces serine for the cytosol and overall folate metabolism is a net producer of serine for gluconeogenesis.

              F. Eliminating the mitochondria

              We examined the significance of the mitochondrial folate cycle for overall one-carbon metabolism by eliminating the mitochondrial folate reactions entirely (setting all the velocities to zero). The significant changes in cytosolic one-carbon metabolism are indicated in Figure 6. There is now no formate at all being supplied by the mitochondria so the rates of thymidylate synthesis and purine synthesis are greatly reduced. The rate of the cytosolic SHMT increases approximately 18-fold and the cytosolic MTD and MTCH reactions reverse and now run in the direction from 5,10-CH2-THF to 10f-THF (not shown). The cytosolic concentration of 5,10-CH2-THF drops, causing the concentration of 5mTHF to drop, which in turn causes homocysteine to rise. The concentrations of methionine and SAM are not much affected. It is interesting to observe that the elimination of the mitochondrial folate metabolism does not completely disrupt cytosolic folate metabolism; in fact the only major effects are on the rates of thymidylate and purine synthesis, while the behavior of the methionine cycle (including the DNMT reaction) is largely unaffected. This is consistent with the hypothesis that the main role of mitochondrial folate metabolism is to supply extra 1-carbon units to the cytosol as formate [16, 17, 24].
              Figure 6

              Change in the steady-state concentrations of selected metabolites and steady-state reaction velocities of selected reactions when the mitochondrial reactions are eliminated form the model.

              E. Methionine and protein loading

              The mathematical model allows us not only to compute steady states but also to compute the time course of concentrations and fluxes as they respond dynamically to changing inputs. We examined how the system responded to a methionine load, which we simulated by doubling the external methionine concentration from 30 to 60 μM during hours 5–10 of a 20 hour simulation (Figure 7A). As expected [27, 29], methionine rises modestly, SAM rises substantially, but the DNA methylation rate is very stable because the extra methyl groups are carried by an increase in the rate of the GNMT reaction during loading. The rise in SAM also causes an increase in the rate of the CBS reaction during loading, so there is an increased removal of serine from the system, which causes both cytosolic and mitochondrial serine concentrations to decrease. Despite the increase in the rate of the CBS reaction, homocysteine almost doubles because SAM is inhibiting BHMT and MTHFR, which decreases 5mTHF and slows the MS reaction. Purine and thymidylate synthesis rise modestly during loading because of modest increases in cytosolic 10f-THF and 5,10-CH2-THF, and there is a decline of 5mTHF accompanied by an increase in cytosolic THF. The mitochondrial folates are virtually unchanged, as is the rate of transport of formate from the mitochondria to the cytosol.
              Figure 7

              A: Response of selected model variables to a 5-hour pulse of elevated external methionine. B: Response to a 5-hour pulse of elevated external methionine, glycine and serine. The pulses consisted of a doubling of the external amino-acid concentrations from 5 to 10 hours after initiation of the simulation.

              To simulate a protein meal, we not only doubled methionine for five hours but also doubled the external serine and glycine concentrations for the same five hours (Figure 7B). The major difference is that all the serine and glycine concentrations rise during protein loading and the mitochondrial SHMT reaction reverses direction. The increase in serine increases the CBS reaction and the fraction transsulfurated, and the increase in glycine increases the GNMT reaction. Both of these changes prevent SAM from rising as high as it did under a pure methionine load. In turn, the inhibition of MTHFR by SAM is diminished, so 5mTHF does not fall as much as during a methionine load. The joint effect is that homocysteine does not rise nearly as much as during a pure methionine load. Purine and thymidylate synthesis rise modestly during protein loading because of modest increases in cytosolic 10f-THF and 5,10-CH2-THF. Compared to methionine loading, there are smaller changes in the cytosolic folates and, as above, the rate of transport of formate from the mitochondria to the cytosol is remarkably constant.


              Understanding the function of the compartmentalization of the folate cycle between the cytosol and mitochondria poses many challenges. Measuring concentrations and especially reaction velocities is not easy in living cells, and these measurements become particularly difficult when substrates and some reactions are sequestered into different compartments. It is especially challenging, and in many cases impossible, to measure several different variables at the same time and to track changing concentrations or velocities over time. A mathematical representation of the substrates and reactions of folate metabolism allows one to conduct in silico experiments to test ideas and hypotheses about how the system as a whole operates. Of course, no such mathematical model can represent the complete physical and biological complexity of a real cell, and the usefulness of the model depends on how accurately it represents the known biology of the cell. Experiments with the mitochondrial model in this paper show that its behavior is consistent with a wide body of experimental findings.

              Among the important questions investigated or discussed by experimentalists in recent years are: the role of the mitochondrial bifunctional enzyme, the effect of cytosolic and mitochondrial SHMT expression on purine and pyrimidine synthesis, the directionality of the SHMT reactions, the relative roles of serine and glycine as one-carbon donors, the effects of protein and methionine loading, the significance of the GDC reaction and, indeed, the significance of the mitochondria themselves. Our mathematical model has allowed us to explore these questions by in silico experimentation and we hope thereby to shed light on these questions and to contribute to the ongoing discussion.

              Particularly interesting is the remarkable stability of formate production in the mitochondria in the face of large changes in serine and glycine input. As a consequence, the cytosol has an almost constant input of formate for purine and pyrimidine synthesis, despite short term and long term variations in glycine and serine availability. This effect is largely due to the efficient rebalancing of the serine and glycine concentrations by the reversible mitochondrial and cytosolic SHMT reactions.

              Another highlight that we found is that eliminating the mitochondrial bifunctional enzyme did not lead to a runaway accumulation of 5,10-CH2-THF in the mitochondria as might be expected. Instead, the mitochondrial SHMT reaction reverses direction strongly and the mitochondria become net exporters of serine. Thus, in the presence of the bifunctional enzyme (as in embryonic tissues and cancer cells), the mitochondria primarily support cytosolic purine and pyrimidine synthesis via the export of formate, while in adult tissues the mitochondria produce serine for gluconeogenesis.

              We plan to continue our investigations by conducting a complete fluctuation analysis [29] of the stability of formate production and purine and pyrimidine synthesis. In addition, we will continue to use the model to explore other important aspects of folate biochemistry and compartmentalization by studying the roles of folate-protein binding, substrate channeling, and the interactions of simultaneous genetic polymorphisms with variation in dietary input and vitamin status. We plan to investigate the dynamic properties and potential interactions among the many methyltransferases that act in parallel using SAM as a substrate. Finally, we plan to investigate the dynamic transport processes that determine the compartmentalization of the metabolites of folate and methionine cycles between the liver and the blood.


              Additional File 1: 

              The attached Additional File 1 contains a complete description of the mathematical model as well as full names of all abbreviations used in the text and in the model.



              We would like to thank Marian Neuhouser, Jill S. James, and Alanna Boynton for their many contributions to the mathematical modeling project, and Dean Appling for useful comments on a draft of this paper. This work was supported by NIH grant RO1 CA 105437(CMU) and NSF grant DMS 0109872 (MCR).

              Authors’ Affiliations

              Department of Biology, Duke University
              Department of Mathematics, Duke University
              Department of Nutrition Sciences and Toxicology, University of California
              Department of Food Science and Human Nutrition, University of Florida
              Cancer Prevention Program, Fred Hutchinson Cancer Research Center


              1. Bailey LB: Folate in Health and Disease. 1995, New York, NY: Marcel Dekker
              2. Choi SW, Mason JB: Folate and carcinogenesis: an integrated scheme. J Nutr. 2000, 130: 129-132.PubMed
              3. Lucock M: Folic acid: nutritional biochemistry, molecular biology, and role in disease processes. Molec Genet Metabol. 2000, 71: 121-138. 10.1006/mgme.2000.3027.View Article
              4. Refsum H: Folate, vitamin B12 and homocysteine in relation to birth defects and pregnancy outcome. Br J Nutr. 2001, 85 (Suppl 2): S109-S113.View ArticlePubMed
              5. Ryan BM, Weir DG: Relevance of folate metabolism in the pathogenesis of colorectal cancer. J Lab Clin Med. 2001, 138: 164-176. 10.1067/mlc.2001.117161.View ArticlePubMed
              6. Ulrich CM: Nutrigenetics in cancer research – Folate metabolism and colorectal cancer. J Nutr. 2005, 135: 2698-2702.PubMed
              7. Giovannucci E: Epidemiologic studies of folate and colorectal neoplasia: a review. J Nutr. 2002, 132 (Suppl 8): 2350S-2355S.PubMed
              8. Robien K, Ulrich CM: 5'10-methylenetetrahydrofolate reductase polymorphisms and leukemia risk. Am J Epidemiol. 2003, 157: 571-582. 10.1093/aje/kwg024.View ArticlePubMed
              9. De Bree A, Verschuren WMM, Kromhout D, Kluijtmans LAJ, Blom HJ: Homocysteine determinants and the evidence to what extent homocysteine determines the risk of coronary heart disease. Pharm Rev. 2002, 54: 599-618. 10.1124/pr.54.4.599.View ArticlePubMed
              10. Morris MS: Homocysteine and Alzheimer's disease. Lancet Neurol. 2003, 2: 425-428. 10.1016/S1474-4422(03)00438-1.View ArticlePubMed
              11. Jackman AL: Antifolate Drugs in Cancer Therapy. 1999, Totowa, NJ: Humana PressView Article
              12. McGuire JJ: Anticancer antifolates: current status and future directions. Curr Pharm Des. 2003, 9: 2593-2613. 10.2174/1381612033453712.View ArticlePubMed
              13. Wang FK, Koch J, Stokstad EL: Folate coenzyme pattern, folate linked enzymes and methionine biosynthesis in rat liver mitochondria. Biochem Z. 1967, 346: 458-466.PubMed
              14. Wagner C: Biochemical role of folate in cellular metabolism. Folate in Health and Disease. Edited by: Bailey LB. 1995, New York: Marcel Dekker, 23-42.
              15. Cook RJ: Folate Metabolism. Homocysteine in Health and Disease. Edited by: Carmel R, Jacobsen DW. 2001, Cambridge: Cambridge Univ Press, 113-134.
              16. Appling DR: Compartmentation of folate-mediated one-carbon metabolism in eukaryotes. FASEB J. 1991, 5: 2645-2651.PubMed
              17. Christensen KE, MacKenzie RE: Mitochondrial one-carbon metabolism is adapted to the specific needs of yeast, plants and mammals. BioEssays. 2006, 28: 595-605. 10.1002/bies.20420.View ArticlePubMed
              18. Lewis KF, Randolph VM, Nemeth E, Frisell WR: Oxidation of 1-carbon compounds to formate and carbon dioxide in rat liver mitochondria. Arch Biochem Biophys. 1978, 185: 443-449. 10.1016/0003-9861(78)90187-X.View ArticlePubMed
              19. Barlowe CK, Appling DR: In vitro evidence for the involvement of mitochondrial folate metabolism in the supply of cytoplasmic one carbon units. Biofactors. 1988, 1: 171-176.PubMed
              20. Garrow TA, Brenner AA, Whitehead VM, Chen XN, Duncan RG, Korenberg JR, Hane B: Cloning of human cDNAs encoding mitochondrial and cytosolic serine hydroxymethyltransferases and chromosomal localization. J Biol Chem. 1993, 268: 11910-11916.PubMed
              21. Kastanos EK, Woldman YY, Appling DR: Role of mitochondrial and cytoplasmic serine hydroxymethyltransferase isozymes in de novo purine synthesis in Saccharomyces cerevisiae. Biochemistry. 1997, 36: 14956-14964. 10.1021/bi971610n.View ArticlePubMed
              22. Stover PJ, Chen LH, Suh JR, Stover DM, Keyomarsi K, Shane B: Molecular cloning, characterization, and regulation of the human mitochondrial serine hydroxymethyltransferase gene. J Biol Chem. 1997, 272: 1842-1848. 10.1074/jbc.272.3.1842.View ArticlePubMed
              23. Herbig K, Chiang E-P, Lee L-R, Hills J, Shane B, Stover PJ: Cytoplasmic serine hydroxymethyltransferase mediates competition between folate-dependent deoxyribonucleotide and S-adenosylmethionine biosynthesis. J Biol Chem. 2002, 277: 38381-38339. 10.1074/jbc.M205000200.View ArticlePubMed
              24. Di Pietro E, Wang X-L, MacKenzie RE: The expression of mitochondrial methylenetetrahydrofolate dehydrogenase-cyclohydrolase supports a role in rapid cell growth. Biochim Biophys Acta. 2004, 1647: 78-84.View Article
              25. Christensen KE, Patel H, Kuzmanov U, Mejia NR, MacKenzie RE: Disruption of the Mthfd1 gene reveals a monofunctional 10-formyltetrahydrofolate synthetase in mammalian mitochondria. J Biol Chem. 2005, 280: 7597-7602. 10.1074/jbc.M409380200.View ArticlePubMed
              26. Walkup AS, Appling DR: Enzymatic characterization of human mitochondrial C1-tetrahydfrofolate synthase. Arch Biochem Biophys. 2005, 442: 196-205. 10.1016/ ArticlePubMed
              27. Reed MC, Nijhout HF, Sparks R, Ulrich CM: A mathematical model of the methionine cycle. J Theor Biol. 2004, 226: 33-43. 10.1016/j.jtbi.2003.08.001.View ArticlePubMed
              28. Nijhout HF, Reed MC, Budu P, Ulrich CM: A mathematical model of the folate cycle. J Biol Chem. 2004, 279: 55008-55016. 10.1074/jbc.M410818200.View ArticlePubMed
              29. Nijhout HF, Reed MC, Anderson DF, Mattingly JC, James SJ, Ulrich CM: Long-range allosteric interactions between the folate and methionine cycles stabilize DNA methylation reaction rate. Epigenetics. 2006, 1: 81-87.View ArticlePubMed
              30. Reed MC, Nijhout HF, Neuhouser ML, Gregory JF, Shane B, James SJ, Boynton A, Ulrich CM: A mathematical model gives insights into nutritional and genetic aspects of folate-mediated one-carbon metabolism. J Nutr. 2006, 136: 2653-2661.PubMed
              31. Allen R, Stabler S, Lindenbaum J: Serum betaine, N-N-dimethylglycine and N-methylglycine levels in patients with cobalamin and folate deficiency and related inborn errors of metabolism. Metabolism. 1993, 42: 1448-1460. 10.1016/0026-0495(93)90198-W.View ArticlePubMed
              32. MacFarlane AJ, Perry C, Liu X, Allen S, Stover PJ: Cytoplasmic serine hydroxymethyltransferase regulates the homocysteine remethylation pathway. Haematologica Reports. 2005, 1: 28-
              33. Titus SA, Moran RG: Retrovirally mediated complementation of the glyB phenotype. Cloning of a human gene encoding the carrier for entry of folates into mitochondria. J Biol Chem. 2000, 275: 36811-36817. 10.1074/jbc.M005163200.View ArticlePubMed


              © Nijhout et al; licensee BioMed Central Ltd. 2006

              This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.