- Research Article
- Open Access
Mathematical analysis of the regulation of competing methyltransferases
BMC Systems Biology volume 9, Article number: 69 (2015)
Methyltransferase (MT) reactions, in which methyl groups are attached to substrates, are fundamental to many aspects of cell biology and human physiology. The universal methyl donor for these reactions is S-adenosylmethionine (SAM) and this presents the cell with an important regulatory problem. If the flux along one pathway is changed then the SAM concentration will change affecting all the other MT pathways, so it is difficult for the cell to regulate the pathways independently.
We created a mathematical model, based on the known biochemistry of the folate and methionine cycles, to study the regulatory mechanisms that enable the cell to overcome this difficulty. Some of the primary mechanisms are long-range allosteric interactions by which substrates in one part of the biochemical network affect the activity of enzymes at distant locations in the network (not distant in the cell). Because of these long-range allosteric interactions, the dynamic behavior of the network is very complicated, and so mathematical modeling is a useful tool for investigating the effects of the regulatory mechanisms and understanding the complicated underlying biochemistry and cell biology.
We study the allosteric binding of 5-methyltetrahydrofolate (5mTHF) to glycine-N-methyltransferase (GNMT) and explain why data in the literature implies that when one molecule binds, GNMT retains half its activity. Using the model, we quantify the effects of different regulatory mechanisms and show how cell processes would be different if the regulatory mechanisms were eliminated. In addition, we use the model to interpret and understand data from studies in the literature. Finally, we explain why a full understanding of how competing MTs are regulated is important for designing intervention strategies to improve human health.
We give strong computational evidence that once bound GNMT retains half its activity. The long-range allosteric interactions enable the cell to regulate the MT reactions somewhat independently. The low K m values of many MTs also play a role because the reactions then run near saturation and changes in SAM have little effect. Finally, the inhibition of the MTs by the product S-adenosylhomocysteine also stabilizes reaction rates against changes in SAM.
Methyltransferase reactions, in which methyl groups are attached to substrates, are fundamental to many aspects of cell biology and human physiology. In mammals there are at least 150 methyl transferases that methylate DNA, RNA, lipids, proteins, and small molecules [1, 2]. For example, the methylation of cytosines in DNA is the basis for epigenetic control, the methylation of guanidinoacetate to form creatine is fundamental to energy metabolism, and the methylation of arsenic is an important detoxification pathway. The universal methyl donor in cells is the molecule S-adenosylmethionine (SAM) that is synthesized from the amino acid methionine in the methionine cycle. This poses an unusual but important regulatory challenge for cells. If one methyltransferase (MT) is upregulated by the cell, then there will be a larger flux on that pathway, which should lower the SAM concentration and thereby decrease the flux through the other methylation pathways. That is, it is not clear how the cell can regulate the different methyltransferase pathways independently. Even 40–50 years ago when many methyltransferases were being discovered, researchers understood that this is an important issue and often referred to “competing methyltransferases” .
Many methyl groups used in the MT reactions come directly from dietary input of methionine, betaine, and choline (to make betaine) in the methionine cycle. See Fig. 1. However, the folate cycle can also contribute methyl groups in a process called methylneogenesis. In the folate cycle, a methyl group is removed from serine (making glycine) in the SHMT reaction and is attached to tetrahydrofolate (THF) making methylene-tetrahydrofolate (CH2-THF). See Fig. 1, and its legend for the full names of enzymes. In turn, CH2-THF is reduced by the enzyme MTHFR to 5-methyltetrahydrofolate (5mTHF). In the methionine synthase (MS) reaction, 5mTHF donates the methyl group to homocysteine (Hcy) to make methionine. As we will see, methylneogenesis is partially controlled by the concentration of SAM itself. Mudd and Poole  showed that if humans are given diets restricted in methionine, choline, and betaine, the subjects maintained the overall flux through the MT reactions and Davis et al.  determined that 100 % of these new methyl groups come from serine. Nijhout et al.  showed, using mathematical modeling, that, in addition to dietary input of serine, the mitochondria are an important source of serine in adult cells by converting glycine to serine.
We include five methyltransferases in our model. GNMT, PEMT, and GAMT are included because together they carry most of the methylation flux [7, 8]. DNMT is included because it is an example of the many methyltransferases that have a very low K m for SAM and it is a methyltransferase that is greatly up regulated when cells divide. Finally, we include AS3MT because the detoxification of arsenic is a special interest of ours. Given the biological importance of the MT reactions, it is not surprising that many regulatory mechanisms have evolved that allow the cell to regulate the methylation fluxes through the different pathways more or less independently. These are: (1) Long-range interactions The long-range interactions are depicted by red arrows in Fig. 1. “Long-range” does not refer to distance inside the cell (we assume that all species are at high enough concentration that we can consider the cell to be well-mixed). We use the term “long-range” when a substrate affects an enzyme at a distant location in the network diagram. The concentration of SAM affects the activity of enzymes at distant locations in the biochemical network. SAM activates CBS [9, 10] and inhibits BHMT . As the SAM concentration rises this has the effect of sending more Hcy down the transsulfuration pathway (CBS) and remethylating fewer Hcy (BHMT) to methionine, thus moderating the increase in SAM. In a series of papers, Finkelstein and Martin [12, 13] emphasized that these effects tend to conserve the total mass in the methionine cycle. In addition, SAM inhibits the enzyme MTHFR . Thus when the SAM concentration increases, the concentration of 5mTHF will drop and the MS reaction will go slower, so less Hcy will be remethylated to methionine.(2) The binding of 5mTHF to GNMT Zamierowski and Wagner  discovered that the GNMT is a major folate binding protein. 5mTHF binds to GNMT and tends to inactivate it; see Fig. 1. Thus, when SAM starts to rise, the inhibition of MTHFR causes the concentration of free 5mTHF to fall causing some of the boundGNMT-5mTHF complexes to dissociate. This makes more free GNMT so the GNMT pathway runs faster preventing the SAM concentration from rising too much. Conversely, if SAM starts to fall, the inhibition of MTHFR is partially relieved creating more 5mTHF. The increased 5mTHF binds to more GNMT, lowering the amount of free GNMT available for the GNMT reaction, thus moderating the decline in SAM concentration. Because of these ideas, Wagner and others refer to the GNMT reaction as a “salvage pathway”. When SAM is high many of the excess methyl groups are carried away in the GNMT reaction and when SAM is low the GNMT reaction is constricted to save the methyl groups for the other MT reactions. This view is consistent with the fact that the product of the GNMT reaction, sarcosine, has no known physiological function. It is transported into the mitochondria where the one carbon units are stripped off and (except for one C O 2) returned to the folate cycle in the cytosol.(3) SAH inhibits the MT reactions Almost all the MT enzymes are inhibited by S-adenosylhomocysteine (SAH), the universal product of the MT reactions. This also tends to maintain the SAM concentration. If one MT is highly upregulated, the increased flux will raise the concentration of SAH, which will inhibit somewhat more the other MT reactions, thus tending to maintain the concentration of SAM.(4) Many MT have very low K m values for SAM Some MT enzymes have K m values for SAM that are of the same magnitude as typical SAM concentrations. This includes GNMT, GAMT, and PEMT, three pathways that carry a major amount of the MT flux. However, many other MT enzymes, for example DNMT, have very low K m values for SAM . This ensures that even if the SAM concentration fluctuates the velocity of the particular MT reaction won’t change much because it is already on the saturated part of the Michaelis-Menten curve.
The main purpose of this study is to investigate, using mathematical modeling, these regulatory mechanisms. The first step is to use information about the binding of 5mTHF to GNMT recently discovered by Luka and Wagner [16–18] so that we can update and improve our previous models of folate and methionine metabolism in the liver [6, 19, 20]. We present, in Section The binding of 5mTHF to GNMT, evidence from the mathematical modeling that GNMT retains 50 % of its activity when it is once bound by 5mTHF but not when it is twice bound. Then we study the five competing methyltransferases indicated in Fig. 1: GNMT, the salvage pathway, PEMT and GAMT, major MT pathways with moderate K m values for SAM, DNMT, which has a low K m value for SAM, and AS3MT, the enzyme in the arsenic detoxification pathway that is a particular interest of ours. In Section Stabilization of methylation reactions by the long-range interactions we show that the long-range interactions moderate the effects on other fluxes when one flux is eliminated or highly upregulated. In Section Stabilization of SAM against changes in methionine input we show how the long-range interactions moderate the effects of changes in methionine input. And, in Section Methylation and folate deficiency we show that even though the SAM concentration depends on total liver folate, the total methylation flux is not markedly affected by folate deficiency. In Section The importance of small K m values we explain the importance of the small K m values for SAM of many MTs and in Section The effect of inhibition by SAH we study the inhibition of the MT reactions by SAH. In Section The importance of substrate availability we point out that individual methylation fluxes can be controlled by co-substrate availability as well as enzyme expression level. In Section Arsenic in Bangladesh, we explain why a full understanding of how competing MTs are regulated is important for understanding and designing intervention strategies for human health effects. We also explain that there are times when one would like to break some or all of these homeostatic mechanisms so as to achieve specific changes in the balance of the MT pathways. In addition, we use the model to interpret and understand data from studies in the literature.
A schematic diagram of the mathematical model is shown in Fig. 1. The pink boxes indicate the variable substrates in the folate cycle and the green boxes indicate the variable substrates in the methionine cycle. Each arrow represents a biochemical reaction and the blue ellipses on the arrows give the acronyms of the enzymes catalyzing the reactions. Full names of substrates and enzymes are given in the legend to Fig. 1. In addition, the concentration of the free enzyme, GNMT, is treated as a variable in the model, as well as its singly bound and doubly bound forms, 5mTHF-GNMT and 5mTHF-GNMT-5mTHF, respectively. The concentration of each variable satisfies a differential equation that represents mass-balance; that is, the rate of change of the variable (in μM/hr) is simply the sum of the rates of the arrows coming into the variable minus the sum of the rates of the arrows leaving the variable. Thus, the differential equation for [ S A M] is
Each V is a velocity and the subscript indicates which reaction. Each velocity depends on the current value of the concentrations of the indicated variables. The complicated, usually non-linear, form of this dependence depends on the detailed biochemistry of each of the enzymes. Each of the methylation reactions also depends on the concentration of the substrate being methylated, for example glycine for GNMT and guanadino-acetate for GAMT, though we do not indicate these explicitly. These methylation substrates are taken to be constant in the model because our main interest is the regulation of [ S A M] concentration.
The complete mathematical model that consists of 13 differential equations and explicit formulas for each reaction velocity is given in the Additional file 1. Large parts of the model are very similar to the liver folate and methione cycle parts of (larger) models that we have developed in the past to study different aspects of cell metabolism. The contributions of the mitochondria to one-carbon metabolism were studied in . The transulfuration pathway was studied in  and whole body folate and methione metabolism was studied in . In the rest of this section we present the new reactions and substrates that we have added in order to study the problem of competing methyltransferases. Full details of the complete model are available in the Additional file 1.
Binding of 5mTHF to GNMT
In a series of papers, Wagner, Luka, and colleagues have studied the inhibitory effect of 5mTHF on the activity of GNMT [17, 18, 23–25]. GNMT has two binding sites for 5mTHF, so we assume the simple reversible reactions:
with forward and backward rate constants, k 1 and k 2, for the first reaction and k 3 and k 4, for the second reaction. Thus the differential equations for [ 5m T H F],[ G N M T],[ 5m T H F -GMNT],[5m T H F -GNMT-5m T H F] are:
We choose the rate constants k 1=50,k 2=1,k 3=1,k 4=1.6 so that the K D values are those found in Table two of .
The kinetics of GNMT for SAM are cooperative and we take the Hill coefficient to be n=2 as suggested in  and we use K m =100 μM as indicated in . The inhibition by SAH is competitive  and has K i =35 μM . Of course, the reaction has glycine as a substrate but we take the glycine concentration to be constant so it’s effect is included in V max =4000 μ/hr. Thus,
This formula for V max resulted from our in silico experiments described in the first section of Results. The concentration of free GNMT, [GNMT], is a variable in our model. GNMT can be bound by one or two molecules of 5mTHF. Our simulations and the data in  suggest strongly that once bound GNMT, namely [ 5m T H F-GN M T], retains 50 % of it’s activity. The factor 4000 is chosen so that the GNMT has a normal reaction velocity comparable to the reaction velocities of PEMT and GAMT, the two other methyl transferases that carry much of the methylation flux.
Inorganic arsenic is metabolized in two methylation steps catalyzed by AS3MT. The first step uses utilizes a methyl group from SAM and is followed by a reduction step to produce methylarsonic acid (MMA). The second step uses utilizes a methyl group from SAM and is followed by a reduction step to produce dimethyarsinic acid (DMA), which is readily exported from the liver and cleared in the urine. We have recently studied the biochemistry of these methylation steps that are quite complicated . For, example the first step shows substrate inhibition by inorganic arsenic and product inhibition by MMA and glutathione (GSH) both sequesters the arsenic compounds and activates AS3MT. In our study here, we are mainly interested in studying the availability of methyl groups from SAM, so we take the arsenic concentrations and the GSH concentration to be constant, and model just the first methylation step. SAM shows substrate inhibition for AS3MT , but the effect is small and occurs only at very high SAM concentrations, so we ignore it. Thus, the velocity of methylation is taken to be:
We take the K m of AS3MT for SAM to be 50 μM as determined in . It is known that SAH inhibits AS3MT [29, 30], but the nature of the inhibition and the K i are not known. We’ll assume the inhibition is competitive and take K i =10 μM, which is typical of other methyltransferases. A high, but realistic arsenic load is 1 μM in liver  and a typical flux would be the order of magnitude of 1 μM/hr. So, we choose V max =5 μ/hr, which accomplishes this given that a typical SAM concentration is 24 μM.
The velocity of the PEMT reaction is given by
The inhibition by SAH is non-competitive . We choose K m =18.2 μM for SAM and K i =3.8 μM for SAH as indicated in . The reaction depends on phosphatidylethanolamine but since we take its concentration to be constant we fold that dependence into the V max . The value V max =250 μM/hr was chosen so that the flux of the PEMT reaction is comparable to the fluxes of the GNMT and GAMT reactions, the two other methyl transferases that carry much of the methylation flux.
The velocity of the GAMT reaction is given by
The inhibition by SAH is competitive [33, 34]. We choose K m =49 μM for SAM and K i =16 μM for SAH as indicated in . The reaction depends on the effective concentration of the cytosine substrates, but since we take its concentration to be constant we fold that dependence into the V max . The value V max =100 μM/hr was chosen so that the flux of the GAMT reaction is comparable to the fluxes of the GNMT and PEMT reactions, the two other methyl transferases that carry much of the methylation flux.
The velocity of the DNMT reaction is given by
The inhibition by SAH is competitive . We choose K m =1.4 μM for SAM and K i =1.4 μM for SAH as indicated in . The reaction depends on guanidinoacetate but since we take its concentration to be constant we fold that dependence into the V max . The value V max =12.5 μM/hr was chosen so that the flux of the DNMT reaction is normally (when the cell is not dividing) much less than the fluxes of GNMT, PEMT, and GAMT.
The fluxes of each of the MT reactions depends on the availability of co-substrates, for example phosphatidylethanolamine for the PEMT reaction or guanidinoacetate (gaa) for the GAMT reaction. We take these co-substrates to be constant in our simulations, so, for simplicity, the co-substrate terms are not shown in the above formulas. Full formulas with the co-substrates in the cases of AS3MT, GNMT, GAMT, and PEMT are given in the Additional file 1. The V max values given here and in the Additional file 1 differ somewhat because of the presence of these extra terms.
On the use of Michaelis-Menten kinetics
Almost all of our formulas for reactions velocities have Michaelis-Menten form, for example V=V max [ S]/(K m + [ S]), modified to include the effects of product inhibition, substrate inhibition, or long-range allosteric interactions. All of the explicit formulas are in the Additional file 1. These formulas are simplifications of underlying biochemical processes that are more complicated. For example, the actual sequence of steps of the MS reaction, worked out by Rowena Matthews and others, is very long and complicated (and involves SAM). But we are mainly interested in how overall reaction velocity depends on the concentrations of Hcy and 5mTHF, so we use the simple bi-molecular Michaelis-Menten formula until we think that the details of the mechanism are important for the questions we are trying to answer. Deciding what level of biological and biochemical detail to include is always a difficult and important issue. We use simple formulations until the data that we are trying to understand suggests to us that the details on a lower level may be relevant.
Results and discussion
The binding of 5mTHF to GNMT
It was established as early as 1985 by Wagner and co-workers  that 5mTHF inhibits the enzyme GNMT; 15 % inhibition was shown at 0.1 μM of 5mTHF pentaglutamate, 50 % inhibition at 1 μM and 90 % inhibition at 10 μM. Since then, Wagner, Luka, and co-workers have investigated the structure of GNMT and the binding of 5mTHF. In  it was shown that the binding of SAM to GNMT is cooperative and depends on the 5mTHF concentration and GNMT and shows substrate inhibition for SAM at high SAM concentrations. In [16, 37] the structure of GNMT was elucidated and it was shown that GNMT can bind two molecules of 5mTHF. In , KD values for the binding of 5mTHF to GNMT are given and inhibition curves and free versus bound 5mTHF curves are given.
Left open was the question of what the activity of the enzyme GNMT is when it is once bound or twice bound with 5mTHF. To investigate this question we used the four differential equations given in Methods for [ 5m T H F],[ G N M T],[ 5m T H F-G N M T] and [ 5m T H F-G N M T-5m T H F]. GNMT is 1–3 % of liver protein  and in  the GNMT concentration in liver was measured to be 1 μM, so we assume that GNMT concentration in our in silico experiments and take the KD values from . Our results are shown in Fig. 2. Panel a shows the concentrations of [ G N M T],[ 5m T H F-G N M T] and [ 5m T H F-G N M T-5m T H F] as a function of the total concentration of [ 5m T H F] in the medium. Panel b shows our computed model curve of bound 5mTHF versus free 5mTHF along with the data points from Figure 2 of . The match of the computed curve to the data is very good. Note that since there is 1 μM of GNMT in the mixture and the curve and the data rise to 1.7 μM, this shows conclusively that two molecules of 5mTHF can bind to each GNMT. Finally, in Panel c we show the data points from the inhibition curve given in Figure 3 of  and two different model curves. The red curve shows the % activity of GNMT assuming that free GNMT has 100 % activity and once bound and twice bound GNMT have 0 % activity. The curve is nowhere near the data. On the other hand, the green curve shows the computed % activity of GNMT assuming that free GNMT has 100 % activity, once bound GNMT has 50 % activity, and twice bound GNMT has 0 % activity. The computed curve lies very close to the data points. We conclude from these computational experiments that it is very likely that once bound GNMT retains 50 % of its activity.
Stabilization of methylation reactions by the long-range interactions
We performed various computational experiments that show the effects of the long-range interactions. Since our results are presented as percentage change from normal, we begin by briefly discussing the normal steady state concentrations and velocities in the model that are indicated in Table 1. Our mathematical model, shown schematically in Fig. 1, differs from previous models [19, 20, 22] in three ways: (1) We have added the PEMT, GAMT, and AS3MT methylation reactions; (2) We treat the binding of 5mTHF to GNMT in much more detail; (3) We have lowered the methionine input to the liver to 50 μM/hr so that the steady state concentration of SAM (23.97) is in the range indicated for humans in [38, 39]. The methionine concentration is correspondingly low. Total folate concentration in the liver is 20 μM  and the 5mTHF concentration is 5.35 μM . The fraction of the flux around the methionine cycle that enters the transsulfuration pathway (“frac”) is 0.5 .
The flux around the methionine cycle, 100.25 μM/hr, is given by the net flux of the reaction from SAH to Hcy. The three major methylation fluxes are those catalyzed by GNMT, PEMT, and GAMT [7, 8] and they carry most of the normal flux in our model. The fluxes through DNMT and AS3MT are small representing two of the many other methyltransferases. The remethylation flux from Hcy to Met in the liver is almost equally divided between MS and BHMT and the CBS flux at steady state must be 50 μM/hr since that is the methionine input.
Our first computational experiments investigate how sudden drastic changes in the amounts of PEMT, GAMT, GNMT, or DNMT affect the other fluxes. In Panel a of Fig. 3, we show the percentage change of each of the other fluxes if one flux is set to zero or if the amount of DNMT is multiplied by a factor of 20. Note that a flux could become zero if either the enzyme expression level goes to zero or the co-substrate concentration is zero. Across the bottom we indicate which of the fluxes has been altered. As expected, if one enzyme is set to zero the other fluxes go up and if DNMT is upregulated by a factor of 20, the other fluxes go down. The percentage changes when PEMT or GAMT is set to zero (appoximately 20 %) are consistent with the PEMT and GAMT knockout changes discussed in . The most dramatic changes come when GNMT is knocked out, which was also seen in . When DNMT is upregulated by a factor of 20, the three fluxes of GNMT, PEMT, and GAMT decrease modestly. One can clearly see the influence of GNMT as the “salvage” pathway first described by Wagner, Briggs, and Cook . If PEMT is knocked out GNMT flux goes up a lot so that the GAMT flux does not rise very much. Similarly, when DNMT is upregulated, the GNMT flux drops dramatically and so PEMT and GAMT do not drop as much.
Panel b of Fig. 3 shows what happens if we do the same experiments but turn off the long-range interactions. In every case the changes to the fluxes are larger when one flux is knocked out or (in the case of DNMT) upregulated. Some differences are quite large. If GNMT is knocked out, then GAMT increases by 55 % in the presence of the long-range interactions and increases by 121 % with no long-range interactions.
By examining Fig. 1 and Panel a of Fig. 4, we can see how these long-range interactions work. Consider the case where PEMT is knocked out. SAM goes up 25 % inhibiting MTHFR more and thus lowering 5mTHF by 14 %. This drives the reaction of 5mTHF with GNMT towards dissociation so there is less double bound GNMT and more single bound GNMT (by about 10 %), which increases the flux in the GNMT reaction. This is why (see Fig. 3, Panel a) the GNMT flux goes up a lot, but the GAMT flux goes up less (about 20 %) and DNMT flux goes up very little. Finally, the fraction of flux transsulfurated goes up, because increased SAM stimulates CBS and inhibits BHMT, and this moderates the increase in SAM. The analogous scenario happens when GAMT is knocked out.
The opposite scenario happens if DNMT is upregulated by a factor of 20. The steady state value of SAM goes down by 22 %, which removes inhibition from MTHFR, so 5mTHF goes up by 18 %. This drives the reaction of 5mTHF with GNMT towards binding and so the amount of once bound GNMT goes down. Thus less flux goes through the GNMT reaction allowing the the PEMT and GAMT reactions to decrease more modestly (Fig. 3, Panel a). Finally, since SAM goes down, the fraction transulfurated decreases by 10 %, which recycles more methyl groups back to methionine for the transmethylation reactions.
This beautiful mechanism is broken if GNMT is knocked out. As one can see in Panel a of Fig. 4, SAM increases dramatically and this causes very large increases in the fluxes through PEMT, GAMT, and DNMT (Panel a of Fig. 3). And if there are no long-range interactions at all (Panel b of Fig. 4), SAM goes way up or way down without changing GNMT-5mTHF, 5mTHF, or the fraction transsulfurated much. The large changes in SAM drive the increased changes in the fluxes seen in Panel b of Fig. 3.
Stabilization of SAM against changes in methionine input
Since SAM is the universal methyl group donor, its concentration determines the availability of methyl groups for the transmethylation reactions. This presents a problem for the cell since the precursor of SAM is methionine and the methionine concentration in the blood varies dramatically during the day due to meals. We have previously carried out a fluctuation analysis that showed how the long-range interactions greatly dampen the fluctuations in methionine input and cause the flux through the DNMT reaction to remain quite stable . Here we show how the steady state value of SAM depends on methionine input (Fig. 5). With the long-range interactions present (green curve) the steady state value of SAM increases fairly slowly and linearly as methionine input increases from 10 μM/hr to 100 μM/hr. When the long-range interactions are removed (red curve), the steady state value of SAM increases much more rapidly and accelerates as the methionine input increases. In both cases the SAM concentration increases dramatically at high values of methionine input. This is the behavior seen in the models of [40, 41] where they call it a “substrate switch,” and the in vitro experiments . Figure 5 shows that the switch occurs much earlier without the long-range interactions. These results are consistent with the classic experiments of Finkelstein and Martin who fed rats diets with different amounts of methionine and measured liver SAM (Table 1 in ). As the methionine content of their diets increased, rats showed relatively modest increases in the SAM concentrations in their livers until at very high methionine diets SAM increased dramatically. All homeostatic mechanisms, including the long-range interactions discussed here, can be broken if inputs become too extreme, a phenomenon that is discussed further in .
Methylation and folate deficiency
There is a positive statistical association between liver folate status and liver SAM concentration . Using a mathematical model of folate and methionine metabolism, we showed that over normal physiological ranges liver SAM concentration is a linear increasing function of total liver folate  and later more detailed models show the same behavior . Figure 6, Panel a, shows that this is true both in the presence of and without the long-range inhibitions. The reason is easy to understand. Higher total folate means higher 5mTHF, which drives the MS reaction faster and remethylates more Hcy to become Met. Much more interesting is the dependence of total methylation rate (i.e., the sum of all five fluxes in the model) on folate status as shown in Panel b of Fig. 6. In the presence of the long-range interactions, there is almost no decrease in total methylation flux (despite the decrease in SAM) as folate deficiency becomes more and more severe (green curve). In contrast, total methylation flux goes down rapidly as folate deficiency increases if the long-range interactions are absent (red curve). Thus, one evolutionary purpose of the long-range interactions may have been to protect the methylation reactions against the folate deficiencies that surely developed during Fall and Winter when green vegetables were scarce.
For all of our model simulations it is natural to ask how sensitive the results are to the exact choices of parameters. The answer is not much. There are 11 V max values for the velocities around the methionine cycle. We conducted experiments in which each these V max values was modified by multiplying it by a number chosen randomly between 0.9 and 1.1 (maximal 10 % variation). All the V max values were changed simultaneously and independently. The green and red shaded areas in Fig. 6 b are the envelopes of the values we obtained in 10 such simulations. As one can see, the fundamental conclusion of Fig. 6 b remains the same, whichever set of V max values one has chosen: the long-range interactions stabilize total methylation rate against folate deficiency. The individual variation of enzyme parameters like V max and K m values is not known. We chose 10 % not to represent real biological variation, but simply to show that the model results and conclusions do not depend sensitively on the exact choices of parameters.
The importance of small K m values
A notable feature of the simulations shown in Fig. 3 is that the knockouts of PEMT, GAMT, and GNMT induce only very small changes in DNMT. The reason is easy to understand. The K m of DNMT is very small, 1.4 μM , compared to the normal concentration of SAM, 23.94 μM, in the model. This means that moderate changes in SAM will hardly affect the velocity of the DNMT reaction, because the velocity curve is very flat for values of SAM considerably higher than the K m . One can see this effect clearly in Fig. 7, which shows hourly changes in methionine input (Panel e) due to three daily meals. GNMT (Panel c) and GAMT (Panel d) show large variations due to meals but DNMT (Panel d) shows almost no variation. This is a very simple regulatory mechanism, but it is very common. Of the 23 K m values for methyltransferases listed in , 12 have K m values ≤ 3, and 4 have K m values in the range 3–10 μM. If a methyltransferase has a low K m value, the cell can regulate its flux simply by changing the expression level of the enzyme, independent of the fluxes of the other methyltransferases or the concentration of SAM.
The effect of inhibition by SAH
Most methyltransferases are inhibited by SAH, the other product of the MT reaction. Of course, product inhibition is ubiquitous in cell metabolism and this simple mechanism is rightly regarded as a way of preventing the over-accumulation of product. However, the graph of PEMT (Panel d) in Fig. 7 suggests a different role for inhibition by SAH. Notice that GNMT and GAMT fluctuate up and down considerably with methionine input, while PEMT does not. All three MTs have high K m values so the discussion in Section The importance of small K m values is not relevant. The reason that PEMT fluctuates little is the inhibition by SAH. As indicated by Panels a and b, the SAH concentration fluctuates up and down with the SAM concentration. For PEMT, the increase in SAM (which will increase flux) is almost exactly compensated by the increase in SAH (which will decrease flux). But why doesn’t this happen for GAMT? The reason is that the K i of PEMT for SAH is quite small (3.8 μM) so the inhibition by SAH is quite powerful, while the K i of GAMT for SAH is quite large (16 μM) so the inhibition by SAH is quite weak. Thus, evolution may have tuned the K m and the K i values of certain methyltransferases for which it was important to buffer the flux against changes in methionine input due to meals.
The importance of substrate availability
In this paper we have been studying the competition between different methyltransferase reactions for the methyl groups provided by SAM. Each of the methyltransferase reactions has a co-substrate to which the methyl group is given, for example phosphotidylethanolamine for the PEMT reaction or guanadinoacetate (gaa) for the GAMT reaction. In our simulations in Section Stabilization of methylation reactions by the long-range interactions we looked at the effects of setting one of the methylation fluxes to zero. This could happen in two ways. The expression level of the enzyme could go to zero or the concentration of the co-substrate could become zero. In general, each methylation flux can be regulated by adjusting expression level and/or substrate level. Physiological systems typically have many layers of control mechanisms and the long-range interactions discussed here are just one set of such control mechanisms in liver metabolism. Each particular methyltransferase reaction will be regulated by other mechanisms also, some of which occur in other organs. This is nicely illustrated by the GAMT reaction.
Creatine is produced by a two step process. In the first step the enzyme L-arginine:glycine amidinotransferase (AGAT) makes gaa in the kidney from arginine and glycine. gaa is exported into the plasma and is taken up by the liver where it receives a methyl group from SAM and becomes creatine via the GAMT reaction. The creatine is exported from the liver and is taken up by muscle tissue. Normally mammals ingest about 50 % of the creatine they need and synthesize the other half through these reactions [46, 47]. Thus, one would expect that creatine synthesis in the liver would be sensitive to the amount of creatine in the diet, and indeed it is [46, 47]. High plasma concentrations in the plasma inhibit AGAT in the kidney so that less gaa is sent to the liver and thus less creatine is produced by the GAMT reaction. The second set of bar graphs in Fig. 3 a shows that if gaa is set to zero (for example, after a high creatine meal) and therefore the GAMT flux becomes zero, then the other MT reactions change modestly, except for GNMT that changes a lot. See the discussion in Section Stabilization of methylation reactions by the long-range interactions.
Arsenic in Bangladesh
In this section we discuss why a full understanding of the regulation of competing methyltransferases is important for human health effects and interventions. Arsenic in drinking water is a major health hazard to millions of people in South and East Asia and in other parts of the world [48, 49]. Long term arsenic exposure has been linked to cancer, heart disease, neuropathies and neurological sequelae, and to deficits in intelligence in children [50, 51]. Arsenic in water is normally ingested primarily as trivalent inorganic arsenic (iAs), which then undergoes hepatic methylation to methylarsonic acid (MMAs) and a second methylation to dimethylarsinic acid (DMAs) by the enzyme AS3MT.
This is considered a detoxification pathway because DMAs is rapidly exported from the liver to the blood and excreted in urine. Two of us (MVG and MNH) conduct studies whose purpose is to determine which nutritional supplements would increase the rate of the AS3MT reaction and thus lower the arsenic body burden in individuals in Bangladesh. It is known both from experimentation  and from modeling (Panel a of Fig. 6) that an increase in folate status increases the concentration of SAM in hepatic cells. Thus one might predict that increasing folate status would increase the rate of methylation of iAs; results from a randomized clinical trial of folate supplementation to folate-deficient adults in Bangladesh show this to be true [53, 54]. Because iAs and its methylated metabolites, MMAs and DMAs, are not measured in the livers of human subjects but in blood and in urine, it is important to have a whole body mathematical model that connects arsenic metabolism in the liver to the blood and urine concentrations of iAs, MMAs, and DMAs. We created such a model in . The model predictions matched the observation in the Bangladesh trials that folate supplementation of folate-deficient individuals reduced blood arsenic by 14 %. With the model, we were also able to predict that body stores of arsenic (which couldn’t be measured) would be reduced by roughly 26 %.
This led naturally to the question of whether supplementation with other nutrients, besides folate, would improve these results. The goal is to raise SAM concentrations and thus to increase the availability of methyl groups for AS3MT. A natural candidate is creatine supplementation, which has been shown to downregulate creatine biosynthesis through repression of arginine:glycine amidinotransferase (AGAT) in the kidney, which, in turn, lowers the concentration of guanidinoacetate (GAA), the substrate for GAMT in the liver [50, 51]. And this brings us immediately to the topic of the investigation in this paper, the competition between methyltransferases. If GAMT flux goes down, how much will SAM go up and how much will that increase the flux through AS3MT? We recently conducted some (unpublished) experiments with our mathematical models that suggested that creatine supplementation (in addition to folate supplementation) would lower blood arsenic an additional 8.5 %. The results in this paper help us to understand why the predicted change was modest. The long-range interactions buffer the flux through other methyltransferases when one methyl-transferase is downregulated. In a recent controlled trial in Bangladesh, 622 participants were randomized to receive 400 μg folic acid, 800 μg folic acid, 3 g creatine, 3 g creatine + 400 μg folic acid, or placebo daily . The 622 participants were comprised of both folate sufficient and deficient individuals. All participants received a water filter that removed arsenic at the start of the study. The results show that 800 μg folic acid was the only treatment to lower total blood arsenic to a significantly greater extent than placebo. However, 3 g creatine+400 μg folic acid did lower blood arsenic to a greater extent (14 %) than 400 μg folic acid alone (3.7 %) although this finding did not reach statistical significance (p=0.08) . Furthermore, 3 g creatine + 400 μg folic acid led to decreases in GAA that were correlated with decreases in homocysteine (manuscript in preparation). Treatment effects on As methylation patterns are presently being analyzed. In general, these findings are consistent with the model predictions. Thus, to improve the outcomes, one needs to break one or more of the regulatory mechanisms. A possible way to do this is to give very high levels of folate supplementation along with the creatine supplementation. If one can raise liver folate appreciably, one could shut down the GNMT pathway because of the binding of 5mTHF to GNMT, and then creatine supplementation should have a larger effect. The difficulty is that it is easy to raise plasma folate a lot but it is not known whether this causes liver folate and liver SAM to rise a lot in folate sufficient individuals.
Liver cells are not static but are constantly changing dynamic systems because the cells must continue to operate in the face of dramatically changing inputs due to meals. It is not surprising that liver cells have evolved many complicated and ingeneous mechanisms for accomplishing the homeostasis of important reactions like the methyltransferase reactions. The purpose of our discussion of AS3MT is to illustrate that a thorough understanding of how competing methyltransferases are regulated is necessary for designing successful interventions.
There are more than 150 methyltransferase reactions in which a methyl group is transferred from SAM to another substrate. These reactions are important steps in a wide variety of biological mechanisms such as the biosynthesis of a large variety of useful compounds, methylating the cytosines on DNA, and catabolizing neurotransmitters. In any given cell typically only a certain subset of the MTs are expressed depending on cell type, its function, and the tasks that the cell is currently doing. The primary mechanism used by the cell for regulating the MTs is to upregulate and downregulate the expression level of the genes that code for the MTs. However, since all the MTs use the same substrate, SAM, up and down regulation of a particular MT will affect the flux in all MT reactions, i.e. it is hard for the cell to regulate each of MT fluxes independently. It is not surprising that cells have developed a number of different mechanisms that help them meet this challenge. In addition, cells must cope with large changes in amino acid input (i.e. methionine) that happen on a shorter time scale than gene regulation.
In this paper, we have used a mathematical model to investigate the mechanisms that liver cells use to manage the non-independence of the MT fluxes.
We have found that long-range allosteric interactions in which substrates in one part of the network affect enzyme activity in distant parts of the network play an important role:
The long-range interactions ensure that when one MT is greatly upregulated or downregulated the GNMT pathway changes a lot so that the fluxes in the other MT pathways change only modestly (Fig. 3). SAM changes much more in the presence of MT up and down regulation without the long-range interactions (Fig. 4).
The long-range interactions act to conserve total mass in the methionine cycle. SAM changes linearly and slowly as methionine input changes (Fig. 5).
We have provided computational evidence that GNMT that is bound to one 5mTHF molecule retains substantial activity. This remains to be confirmed experimentally.
We explained why the small K m values of many MTs for SAM makes them independent of the fluxes through other MT reactions.
We explained why the small K i value of SAH for PEMT stabilizes PEMT flux.
And finally, we discussed why a full understanding of the regulation of MT fluxes is important for human health interventions to ameliorate arsenic toxicity in Bangladesh.
Aminoimidazolecarboxamide ribonucleotide (transferase)
Nonenzymatic interconversion of THF and 5,10-CH2-THF
Methionine adenosyl transferase I
Methionine adenosyl transferase III
- V max :
The maximum velocity of a reaction
- K m :
The Michaelis constant for an enzyme
- K i :
The inhibition constant for an enzyme.
Clarke S, Banfield K. In: (Carmel R, Jacobsen DW, editors.)S-Adenosylmethionine-dependent methyltransferases. Homocysteine in Health and Disease. Cambridge, UK: Cambridge University Press; 2001.
Clarke R, Lewington S, Landray M. Homocysteine, renal function, and risk of cardiovascular disease. Kidney Int Suppl. 2003; 84:131–3.
Kerr S. Competing methyltransferase systems. J Biol Chem. 1972; 247(13):4248–52.
Mudd SH, Poole JR. Labile methyl balances for normal humans on various dietary regimens. Metabolism. 1975; 24:721–35.
Davis SR, Stacpoole PW, Williamson J, Kick LS, Quinlivan EP, Coats BS, et al. Tracer-derived total and folate-dependent homocysteine remethylation and synthesis rates in humans indicate that serine is the main one-carbon donor. Am J Physiol Endocrinol Metab. 2004; 286:272–9.
Nijhout HF, Reed MC, Lam S, Shane B, III JFG, Ulrich CM. In silico experimentation with a model of hepatic mitochondrial folate metabolism. Theor Biol Med Model. 2006; 3:1–11.
Mudd SH, Brosnan JT, Brosnan ME, Jacobs RL, Stabler SP, Allen RH, et al. Methyl balance and transmethylation fluxes in humans. Am J Clin Nutr. 2007; 85:19–25.
Schalinske KTWKL. New insights into the regulation of methyl group and homocysteine metabolism. J Nutr. 2007; 137:311–4.
Janosik M, Kery V, Gaustadnes M, Maclean KN, Kraus JP. Regulation of human cystathionine beta-synthase by s-adenosyl-l-methionine: Evidence for two catalytically active conformations involving an autoinhibitory domain in the c-terminal region. Biochemistry. 2001; 40:10625–33.
Kluijtmans LAJ, Boers GHJ, Stevens EMB, Renier WO, Kraus JP, Trijbels FJM, et al. Defective cystathionine beta-synthase regulation by s-adenosylmethionine in a partially pyridoxine responsive homocystinuria patient. J Clin Invest. 1996; 98:285–9.
Finkelstein JD, Martin JJ. Inactivation of betaine-homocysteine methyltransferase by adenosylmethionine and adenosylethionine. Biochem Biophys Res Commun. 1984; 118:14–9.
Finkelstein JD, Martin JJ. Methionine metabolism is mammals. distribution of homocysteine between competing pathways. J Biol Chem. 1984; 259:9508–13.
Finkelstein JD. Regulation of homocysteine metabolism. Homocysteine Metabolism is Health and Disease. Cambridge, UK: Cambridge University Press: 2001.
Jencks DA, Matthews RG. Allosteric inhibition of methylenetetrahydrofolate reductase by adenosylmethionine effects of adenosylmethionine and nadph on the equilibrium between active and inactive forms of the enzyme and on the kinetics of approach to equilibrium. J Biol Chem. 1987; 262:2485–93.
Zamierowski MM, Wagner C. Identification of folate-binding proteins in rat liver. J Biol Chem. 1977; 252:933–8.
Luka Z, Pakhomova S, Loukachevitch L, Egli M, New Comer ME, Wagner C. 5-methyltetrahydrofolate is bound in intersubunit areas of rat liver folate-binding protein glycine n-methyltransferase. J Biol Chem. 2007; 282:4069–75.
Luka Z, Loukchevitch LV, Wagner C. Acetylation of n-terminal valine of glycine n-methyltransferase affects enzyme inhibition by folate. Biochem Biophys Acta. 2008; 1794(9):1342–6.
Luka Z, Mudd SH, Wagner C. Glycine n-methyltransferase and regulation of s-adenosylmethionine levels. J Biol Chem. 2009; 284:22507–11.
Nijhout HF, Reed MC, Budu P, Ulrich CM. A mathematical model of the folate cycle: new insights into folate homeostasis. J Biol Chem. 2004; 279:55008–16.
Nijhout HF, Reed MC, Anderson DE, Mattingly JC, James SJ, Ulrich CM. Long-range interactions between the folate and methionine cycles stabilize dna methylation. Epigenetics. 2006; 1:81–7.
Reed MC, Thomas RL, Pavisic J, James SJ, Ulrich CM, Nijhout HF. A mathematical model of glutathione metabolism. Theor Biol Med Model. 2008; 5:8.
Duncan TM, Reed MC, Nijhout HF. A population model of folate-mediated one-carbon metabolism. Nutrients. 2013; 5:2457–74.
Wagner C, Briggs WT, Cook RJ. Inhibition of glycine n-methyltransferase acitivty by folate derivatives: Implications for regulation of methyl group metabolism. Biochem Biophys Res Commun. 1985; 127:746–52.
Yeo EJ, Briggs WT, Wagner C. Inhibition of glycine n-methyltransferase by 5-methyltetrahydrofolate pentaglutamate. J Biol Chem. 1999; 274(53):37559–64.
Luka Z, Pakhomova S, loukachevitch LV, newcomer ME, Wagner C. Differences in folate–protein interactions result in differing inhibition of native rat liver and recombinant glycine n-methyltransferase by 5-methyltetrahydrofolate. Biochem Biophys Acta. 2012; 1824:286–91.
Heady JE, Kerr SJ. Purification and characterization of glycine n-methyltransferase. J Biol Chem. 1973; 248(1):69–72.
Lawley SD, Yun J, Gamble MV, Hall MN, Reed MC, Nijhout HF. Mathematical modeling of the effects of glutathione on arsenic methylation. Theor Biol Med Model. 2014; 11:20.
Song X, Geng Z, Li X, Hu X, Bian N, Zhang X, Wang Z. New insights into the mechanism of arsenite methylation with the recombinant human arsenic (?3) methyltransferase (has3mt). Biochimie. 2010; 92:1397–406.
Styblo M, Delnomdedieu M, Thomas DJ. Mono- and dimethylation of arsenic in rat liver cytosol in vitro. Chemico-biological Interact. 1996; 99:147–64.
DeKimpe J, Cornelius R, Vanderholder R. In vitro methylation of arsenite by rabbit liver cytosol: effect of metal ions, metal chelating agents, methyltransferase inhibitors and uremic toxins. Drug Chem Toxicol. 1999; 22:613–28.
Lawley SD, Cinderella M, Hall MN, Gamble MV, Nijhiout HF, Reed MC. Mathematical model insights into arsenic metabolism. Theor Biol Med Model. 2011; 8:31.
Vance DE, Ridgway ND. The methylation of phosphatidylethanolamine. Prog Lip Res. 1988; 27:61–79.
Im YS, Chiang PK, Cantoni GL. Guanidoacetate methyltransferase. purification and molecular properties. J Biol Chem. 1979; 254:11047–50.
Takata Y, Fujioka M. Identification of a tyrosine residue in rat guanidinoacetate methyltransferase that is photolabeled with s-adenosyl-l-methionine. Biochemistry. 1992; 31:4369–74.
Simon D, Grunert F, V Acken U, Doring HP, Kroger H. Dna-methylase from regenerating rat liver: purification and characterisation. Nuceic Acids Res. 1978; 5(6):2153–67.
Flynn J, Reich N. Murine dna (cytosine-5-)-methyltransferase: Steady-state and substrate trapping analyses of the kinetic mechanism. Biochemistry. 1998; 37:15162–9.
Luka Z, Wagner C. Human glycine n-methyltransferase is unfolded by urea through a compact monomer state. Arch Biochem Biophys. 2003; 420:153–60.
Cabrero C, Duce AM, Ortiz P, Alemany S, Mato JM. Specific loss of the high-molecular-weight form of sadenosyl-l-methionine synthetase in human liver cirrhosis. Hepatology. 1988; 8:1530–4.
Lee TD, Sadda MR, Mendler MH, Bottiglieri T, Kanel G, Mato JM, et al. Abnormal hepatic methionine and glutathione metabolism in patients with alcoholic hepatitis. Alcoholism. 2004; 28:173–81.
Martinov MV, Vivitsky VM, Mosharov EV, Banerjee R, Ataullakahanov FI. A substrate switch: A new mode of regulation in the methionine metabolic pathway. J Theor Biol. 2000; 204:521–32.
Korendyaseva TK, Kuvatov DN, Volkov VA, Martinov MV, Vivitsky VM, banerjee R, et al. An allosteric mechanism for switching between parallel tracks in mammalian sulfur metabolism. PLOS Comput Biol. 2008; 4:1–10.
Finkelstein JD, Martin JJ. Methionine metabolism is mammals: Adaptation to methionine excess. J Biol Chem. 1986; 261:1582–7.
Nijhout HF, Best JA, Reed MC. Escape from homeostasis. Math Biosci. 2014; 257(104-110).
Min H, Kim M. A critical evaluation of the correlation between biomarkers of folate and vitamin b12 in nutritional homocysteinemia. Korean J Nutr. 2009; 42:423–33.
Reed M, Nijhout HF, Neuhouser M, III JFG, Shane B, James SJ, et al. A mathematical model gives insights into nutritional and genetic aspects of folate-mediated one-carbon metabolism. J Nutr. 2006; 136:2653–61.
Wyss M, Kaddurah-Daouk R. Creatine and creatinine metabolism. Physiol Rev. 2000; 80:1107–213.
daSilva RP, Nissim I, Brosnan ME, Brosnan JT. Creatine synthesis: hepatic metabolism of guanidinoacetate and creatine in the rat in vitro and in vivo. Am J Physiol Endocrinol Metab. 2009; 296:256–61.
Ahmed MF, Ahuja S, Alauddin M, Hug SJ, Lloyd JR, Pfaff A, et al. Ensuring safe drinking water in bangladesh. Science. 2006; 314:1687–8.
Kinniburgh DG, Smedley PL, Arsenic contamination of groundwater in bangladesh. final report. Technical Report WC/00/19 BritishGeologicalSurvey Keyworth. 2001.
Wasserman GA, Liu X, Parvez F, Ahsan H, Factor-Litvak P, vanGeen A, et al. Water arsenic exposure and children’s intellectual function in araihazar, bangladesh. Eviron Health Pers. 2004; 112:1329–33.
Wasserman GA, Liu X, Parvez F, Ahsan H, Factor-Litvak P, Kline J, et al. Water arsenic exposure and intellectual function in 6-year-old children in araihazar, bangladesh. Environ Health Per. 2007; 115:285–9.
Miller JW, Nadeau MR, Smith J, Smith D, Selhub J. Folate-deficiency-induced homocysteinaemia in rats: disruption of s-adenosylmethionine’s co-ordinate regulation of homocysteine metabolism. Biochem J. 1994; 298:415–9.
Gamble MV, Liu X, Ahsan H, Pilsner JR, Ilievski V, Slavkovich V, et al. Folate and arsenic metabolism: a double-blind, placebo-controlled folic acid-supplementation trial in bangladesh. Am J Clin Nutr. 2006; 84:1093–101.
Gamble MV, Liu X, Slavkovich V, Pilsner JR, Ilievski V, Factor-Litvak P, et al. Folic acid supplementation lowers blood arsenic. Am J Clin Nutr. 2007; 86:1202–09.
Peters BA, Hall MN, Liu X, Factor-Litvak P, Parvez F, V Geen A, et al. Folic acid and creatine as therapeutic approaches to lower blood arsenic: A randomized-controlled trial: Environmental Health Perspectives. http://www.ncbi.nlm.nih.gov/pubmed/25978852 (in press).
This research was partially supported by NSF EF-1038593 (HFN,MR), and NIH grants R01 ES019876 (D. Thomas), RO1 CA133595, RO1 ES017875, P42 ES10349 (MGV) and K99ES018890 (MNH).
The authors declare that they have no competing interests.
The research questions to be addressed were formulated by all four authors. MCR and HFN wrote the code and conducted the simulations. All four authors helped to interpret the results. The manuscript was written by MCR and HFN except for Section The importance of substrate availability, which was written by MVG and MNH. All authors read and approved the final manuscript.
About this article
Cite this article
Reed, M.C., Gamble, M.V., Hall, M.N. et al. Mathematical analysis of the regulation of competing methyltransferases. BMC Syst Biol 9, 69 (2015). https://doi.org/10.1186/s12918-015-0215-6
- Mathematical model
- One-carbon metabolism