Skip to main content

Mathematical modeling of the effects of glutathione on arsenic methylation



Arsenic is a major environmental toxin that is detoxified in the liver by biochemical mechanisms that are still under study. In the traditional metabolic pathway, arsenic undergoes two methylation reactions, each followed by a reduction, after which it is exported and released in the urine. Recent experiments show that glutathione plays an important role in arsenic detoxification and an alternative biochemical pathway has been proposed in which arsenic is first conjugated by glutathione after which the conjugates are methylated. In addition, in rats arsenic-glutathione conjugates can be exported into the plasma and removed by the liver in the bile.


We have developed a mathematical model for arsenic biochemistry that includes three mechanisms by which glutathione affects arsenic methylation: glutathione increases the speed of the reduction steps; glutathione affects the activity of arsenic methyltranferase; glutathione sequesters inorganic arsenic and its methylated downstream products. The model is based as much as possible on the known biochemistry of arsenic methylation derived from cellular and experimental studies.


We show that the model predicts and helps explain recent experimental data on the effects of glutathione on arsenic methylation. We explain why the experimental data imply that monomethyl arsonic acid inhibits the second methylation step. The model predicts time course data from recent experimental studies. We explain why increasing glutathione when it is low increases arsenic methylation and that at very high concentrations increasing glutathione decreases methylation. We explain why the possible temporal variation of the glutathione concentration affects the interpretation of experimental studies that last hours.


The mathematical model aids in the interpretation of data from recent experimental studies and shows that the Challenger pathway of arsenic methylation, supplemented by the glutathione effects described above, is sufficient to understand and predict recent experimental data. More experimental studies are needed to explicate the detailed mechanisms of action of glutathione on arsenic methylation. Recent experimental work on the effects of glutathione on arsenic methylation and our modeling study suggest that supplements that increase hepatic glutathione production should be considered as strategies to reduce adverse health effects in affected populations.


Arsenic is a naturally occuring metalloid that finds its way into the food chain through water, plants, and animals. In many parts of the world, arsenic is a major health hazard [13]. Chronic arsenic exposure has been associated with cancer, heart disease, neuropathies, and with deficits in intelligence in children [4, 5]. Arsenic is mainly ingested as inorganic arsenic, iAs. The metabolism of arsenic in the liver has traditionally been thought to proceed via successive enzymatic methylations to methylarsonic acid, MMAs V, and dimethylarsinic acid DMAs V, with two intervening reduction steps [68]. This is known as the Challenger pathway and the methylations are catalyzed by arsenic methyltransferase, AS3MT. The Challenger pathway has been considered a detoxification pathway because reactive oxygens are replaced by methyl groups and DMAs V is readily exported from the liver and excreted in urine. However, there is considerable evidence that the intermediate trivalent MMAs is equally or more toxic than inorganic arsenic or DMAs V [911].

In recent years, evidence has been accumulating that the tripeptide glutathione, GSH, plays an important role in the Challenger pathway. Since GSH is a reductant, it increases the rates of the reduction steps [1214] and glutathione S-transferase has been shown to help convert MMAs V to MMAs III in different tissues [15, 16]. Thomas, Styblo and colleagues [1719] have studied methylation in the presence of other reductants as well as GSH. Even in the presence of other reductants, GSH increases methylation yield, and Song et al. [20] suggested that GSH increases the activity of AS3MT. In addition, in the experiments of both [20] and [19], it is shown that increasing GSH concentration when the concentration is low increases methylation rate, but increasing GSH concentration when GSH concentration is high decreases methylation rate. Finally, Hayakawa et al. [21] have proposed an alternate pathway for methylation in which only the arsenicals bound to GSH can be methylated.

In a complicated physiological and biochemical situation such as this, mathematical modeling can be a useful tool for sorting out the consequences of different hypotheses and for helping to interpret experimental data. We have made a mathematical model of the arsenic methylation pathway which incorporates three different roles for GSH. First, GSH, as well as other reductants, drives the reduction from valence 5 to valence 3 arsenicals. Second, GSH activates AS3MT. Third, GSH binds reversibly to and sequesters all three arsenic species, iAs, MMAs III, and DMAs III. In our model, the GSH-conjugated arsenicals are not further methylated. We include in the model the many known inhibitions of the methylation reactions and include a new inhibition suggested by the data in [19]. The model is depicted in Figure 1; details are given in Methods. We use the model to explain and interpret the experimental data in [20] and [19]. In particular, we show that the experimental data can be explained by the traditional Challenger pathway with the GSH effects outlined above, so it is not necessary to assume that only GSH-conjugated arsenicals can by methylated as proposed by Hayakawa [21].

Figure 1
figure 1

The reaction diagram. The diagram depicts the traditional Challenger pathway [6] augmented by three effects of glutathione. V 1and V 2are the velocities of the methylation steps. GSH and other reductants increase the velocity of the reduction reactions. GSH increases the activation of AS3MT. GSH conjugates and sequesters the arsenicals iAs, MMAs III, and DMAs III. Although not indicated, the half-life of GSH in reaction mixtures is taken into account. The various inhibitions of iAs and MMAs III are indicated. The kinetics of V 1 and V 2 and the functional form, U (GSH) by which GSH affects AS3MT are given in the Methods.

Cullen [22] discusses the current state of knowledge of methylation of arsenic and outlines four different detailed mechanisms. Considerable knowledge is now available on which cysteine residues in AS3MT are necessary for methylation and on the order of the reaction steps [23, 24]. And, it is known [19] that other thiols besides GSH affect methylation and that there may be an interaction between these thiols and GSH. On the physiological level, GSH is in high concentration in cells and can effect transport processes that control arsenic uptake and removal from cells, as well as the availability of other thiols. Furthermore, GSH is known to bind to xenobiotics, including metals, and, indeed, arsenic-glutathione conjugates appear in the bile of rats fed arsenic containing diets [25], so arsenic conjugation may be an important arsenic excretion pathway. None of these details is in our model.

There are different kinds of experiments, and corresponding models, that shed light on arsenic methylation and arsenic detoxification. There are studies in humans where arsenic metabolites are measured in urine and blood [8, 2629]. There are cell culture experiments in which arsenicals are typically measured in the external medium [10]. And, there are experiments in which reaction mixtures of arsenicals, AS3MT, and various other metabolites are prepared [1721]. A number of pharmacokinetic models have been used to interpret data in these different experimental situations. We have previously constructed a whole body model of arsenic methylation [30] and compared the results to the clinical results of Buchet et al. [26, 27] and the clinical trial of Gamble et al. in Bangladesh [28, 29]. There are other whole body models [3135]. We used a reduced version of our whole body model to study the cell culture experiments in [36]. Previous models for these cell culture experiments were created in [37, 38]. In two recent papers, Georgopoulos and coworkers create mathematical models based on the Hayakawa pathway to study hepatocyte culture experiments including GSH conjugation, reactive oxygen species, and DNA damage [39, 40].

Our model, which investigates the three “effects” of GSH described above and depicted in Figure 1, builds on our previous model of arsenic detoxicfication [30]. Although the model simplifies complicated and interesting biochemical and physiological questions that are the object of current investigations, it enables us to understand three important effects of GSH on arsenic methylation. It is vital to understand the effects of GSH on detoxification mechanisms in hepatocytes, because such understanding may give important information on whether substrates like N-acetyl-cysteine that increase liver GSH may be useful supplements in regions of the world where arsenic is endemic in the water or food supply.


A diagram of the reactions in our model is depicted in Figure 1. The variables in the model are defined in Table 1, followed by the differential equations and Table 2, which gives the values of the rate constants. After Table 2, the functions U,V 1,V 2 in the differential equations are defined and important modeling issues are discussed.

d [ iAs ] dt = U ( [GSH] ) · V 1 ( [ iAs ] , [ MMAs V ] ) - k 1 [ GSH ] 3 [ iAs ] + k - 1 [ ATG ] d [ MMAs V ] dt = k 6 U ( GSH ) · V 1 ( [ iAs ] , [ MMAs V ] ) - ( k 5 + k 6 [ GSH ] ) [ MMAs V ] d [ MMAs III ] dt = ( k 5 + k 6 [ GSH ] ) [ MMAs V ] - U ( [GSH] ) · V 2 ( [ MMAs III ] , [ iAs ] ) - k 2 [ MMAs III ] [ GSH ] 2 + k - 2 [ MMAsG ] d [ DMAs V ] dt = U ( GSH ) · V 2 ( [ MMAs III ] , [ iAs ] ) - ( k 7 + k 8 [ GSH ] ) [ DMAs V ] d [ DMAs III ] dt = ( k 7 + k 8 [ GSH ] ) [ DMAs V ] - k 3 [ DMAs III ] [ GSH ] + k - 3 [ DMAsG ] d [ ATG ] dt = k 1 [ GSH ] 3 [ iAs ] - k - 1 [ ATG ] d [ GSH ] dt = - k 4 [ GSH ] - 3 k 1 [ GSH ] 3 [ iAs ] + 3 k - 1 [ ATG ] - 2 k 2 [ GSH ] 2 [ MMAs III ] + 2 k - 2 [ MMAG ] - k 3 [ GSH ] [ DMAs III ] + k - 3 [ DMAG ] d [ MMAG ] dt = k 2 [ GSH ] 2 [ MMAs III ] - k - 2 [ MMAG ] d [ DMAG ] dt = k 3 [ GSH ] [ DMAs III ] + k - 3 [ DMAG ]
Table 1 Variables in the model ( μ M)
Table 2 Rate constants in the model ( μ M/hr)

Methylation reactions

The velocity, V 1, of the reaction in which iAsIII is methylated to become MMAsV is given by:

V 1 ( [ iAs III ] , [ MMAs V ] ) = V max [ iAs III ] K m + [ iAs III ] 1 + [ iAs III ] K i A 1 + [ MMAs V ] K i M .

We use the value K m = 4.6 μ M for the Michaelis-Menten constant for AS3MT for iAsIII as found in [41]. The reaction has substrate inhibition by iAsIII; we take the inhibition constant to be K i A =1.26μM as found in [42]. It is known that this reaction is inhibited by the product MMAsV and we take the inhibition constant, K i M =40μM from [35] and [43]. We note that it is not certain that the enzyme investigated in [43] is identical to AS3MT.

The velocity, V 2, of the reaction in which MMAsIII is methylated to become DMAsV is given by:

V 2 ( [ MMAs III ] , [ iAs III ] ) = V max [ MMAs III ] K m + [ MMAs III ] 1 + [ iAs III ] K i A 2 1 + [ MMAs III ] 2 ( K i M 2 ) 2 .

As above we take K m = 4.6 μ M and we set K i A 2 =40μM as in [35] and [43]. The inhibition of V 2 by MMAsIII is proposed in this paper; the inhibition constant K i M 2 = 6 μM was obtained by fitting the data in [19]. It is reasonable that the second methylation reaction be inhibited by MMAsIII since the first methylation reaction is inhibited by MMAsV [35, 43], though this doesn’t seem to have been remarked on before. We were driven to include this inhibition by the data in [19], their Figure six, which is discussed in detail under Results. The square gave a much better fit of the data, which suggests that the inhibition is cooperative.

Glutathione as a reductant

It has been known since [12, 44, 45] that GSH acts to reduce pentavalent to trivalent arsenicals. In cells or in vivo other thiols can also act as reductants. We take the rate of the reaction from MMAsV to MMAsIII to be k 5 + k 6[GSH], the k 5 term representing the reduction by other endogenous thiols and the second term representing the reduction by GSH. The concentration of GSH is varied in some of the experiments in [19, 20] and in some of our simulations. We take the rate of the reaction from DMAsV to DMAsIII to be k 7 + k 8[GSH] for similar reasons.

Glutathione affects arsenic methyltransferase

It has been known for a long time that the presence of GSH helps the reduction steps in the methylation chain. The importance of the Styblo data, in [19], Figure six, is that both DMAsIII and total DMAs go up by a factor of about four in the presence of GSH. This shows conclusively that GSH increases substantially the activity of AS3MT. We chose a Hill function for the effect of GSH on AS3MT, U(GSH)= k 9 [ GSH ] 5 ( k 10 ) 5 + [ GSH ] 5 , and the rate constants because they gave a good fit of the data in [19, 20].

Glutathione sequesters arsenic

Arsenic has an affinity for sulfur [12], so it is not surprising that it binds to GSH, especially since a major role of GSH in the liver is to remove xenobiotics including metals. Indeed, arsenic-glutathione compounds can be found in the bile of rats fed arsenic diets [25]. We include in our model the formation of arsenic triglutathione, AsTG, monomethylarsenic diglutathione, MMAsG, and dimethylarsenic glutathione, DMAsG, from iAs III, MMAs III, and DMAs III, respectively. We assume mass-action kinetics and that the reactions are reversible; rate constants are given in Table 2.


The Styblo experiments on MMAsIII

Styblo and colleagues conducted test tube experiments in which MMAsIII was methylated to DMAsV and DMAsIII in the presence of AS3MT and SAM, both with and without 1 mM GSH [19]. Such experiments on intermediates are particularly valuable for understanding the details of a reaction chain. Both DMAsIII and total DMAs were measured and the results are shown in [19], Figure six. Both quantities rise as MMAsIII increases when MMAsIII is low. However, DMAsIII and total DMAs start to decrease as MMAsIII gets still larger, showing clear evidence of inhibition of the second methylation step by MMAsIII. This is reasonable, of course, since MMAsV inhibits the first methylation step, but does not seem to have been commented on before. Our model gives quite good fits to these experiments (see Figure 2), both in the presence and absence of 1 mM GSH.

Figure 2
figure 2

MMAs III inhibits the second methylation step. The red dots are data regraphed from Panels C and D (WT) in Figure six of [19]. The blue curves were computed from the mathematical model. For both the data and model curves, the reaction mixture had either 0 mM GSH (left panel) or 1 mM GSH (right panel).

It has been known for a long time that the presence of GSH helps the reduction steps in the methylation chain. The importance of the data in [19], Figure six is that both DMAsIII and total DMAs go up by a factor of about four in the presence of GSH. This shows conclusively that GSH increases substantially the activity of AS3MT.

Time-course data

In one set of experiments in [19], 1 μM of iAs was introduced into a reaction mixture of volume 100 μl that contained 5 μg of recombinant AS3MT, 1 μM of SAM, and either 0 mM or 1 mM GSH. Other reductants were also in the reaction mixture. Over 40 minutes, the concentrations of MMAsIII, total MMAs, DMAsIII, and total DMAs were measured. In Figure 3 below, the red and green dots reproduce the data from those experiments originally reported in Figure three, panels A and B in [19] with wild type enzyme. In the left panels of our Figure 3, the data points for iAs were calculated by subtracting total MMAs plus total DMAs from the original amount of iAs, namely 1 μM. We note that the measurements of arsenicals in [19] do not distinguish between arsenicals and arsenicals bound to GSH. The blue and black curves in Figure 3 are model calculations of this experimental situation. The blue curves in the left panels are [iAs] + [AsTG]. In the middle panels, the black curve is [ MMAsIII] + [ MMAsV] + [MMAsG], and the blue curve is [ MMAsIII] + [MMAsG]. In the right panels, the black curve is [ DMAsIII] + [ DMAsV] + [DMAsG], and the blue curve is [ DMAsIII] + [DMAsG]. The top row has 0 mM GSH in the reaction mixture and the bottom row has 1 mM GSH in the reaction mixture.

Figure 3
figure 3

Time course experiments. 1 μM iAs was introduced into a reaction mixture containing 0 mM GSH (upper panels) or 1 mM GSH (lower panels). The red and green dots are data points taken from [19]; the curves are the result of model computations. The experimental measurements did not distinguish between arsenicals and arsenicals bound to GSH, so, likewise, the model curves represent the named arsenicals plus their GSH conjugates. For more detail, see the text.

The model curves fit the data points in each panel very well. Note that in the presence of 1 mM GSH more total DMAs is formed and also that there is almost no DMAsV present because, in this experimental context, it is immediately reduced to DMAsIII, most of which is conjugated with GSH.

The influence of GSH on methylation

Both the Styblo group ([19], Figure two) and the Wang group ([20], Figure six) conducted in vitro experiments in which different amounts of GSH were incubated in a reaction mixture for two hours and then the amounts of MMAs and DMAs and their GSH conjugates were measured. As reported above, they did not distinguish between the arsenicals and their GSH conjugates. The reaction mixtures were quite similar except that the Wang group had more AS3MT. The Wang group collected data for 1,3,5,7,10,20 mM GSH and the Styblo group for 1,5,10,20 mM GSH. Their results, which are quite similar, are shown as green dots (Styblo) and red dots (Song) in Figure 4. The connected blue dots are the predictions of our model. As one can see, the model predictions capture well the qualitative behavior of both data sets. At low GSH values and at very high GSH values methylation proceeds slowly, but at intermediate values in the range 5–10 μM methylation proceeds much more quickly. Interestingly, this intermediate level is the physiological range of GSH in human hepatocytes [46, 47].

Figure 4
figure 4

The influence of GSH on methylation. 1 μM iAs was introduced into reactions mixtures with the indicated amount of GSH. MMAs and DMAs (and their GSH conjugates) were measured after two hours. The green dots and red dots are redrawn data from the experiments in [19], Figure two and [20], Figure six, respectively. The connected blue dots, which are the predictions of our mathematical model, capture the qualitative features of the experimental data.

In silico experiments

One of the advantages of mathematical models is that one can add or remove effects of some variables on other variables to see what difference those effects have. Often it is difficult or impossible to perform the corrresponding biological experiments. That is the case here. In our model, GSH has three effects: (i) reduction of arsenicals with valence 5 to valence 3; (ii) activation of AS3MT; (iii) sequestration of arsenicals by binding to GSH. We showed in the previous section (Figure 4) that with these three effects present, the model reproduces well the experimental data of [20] and [19] on the effect of the amount of GSH in the reaction mixture on MMAs, DMAs, and MMAs + DMAs concentrations (left, middle, and right panels in Figure 4). Figure 5 reproduces the right panel of Figure 4 and shows what would happen if the activation of AS3MT by GSH or the sequestration of arsenicals by GSH is eliminated. In our model, the excitation of AS3MT is given by U(GSH) (see Methods), a Hill function. If, instead, we make U (GSH) a constant equal to 5000, a number in the midrange of the values of U, then our computed model curve for MMAs + DMAs would be the black curve in the left panel of Figure 5. This curve, which is monotone decreasing because U does not increase with GSH and more and more of the arsenicals are sequestered by GSH, clearly does not fit the data. On the other hand, if we keep our usual model function U and remove the binding of GSH to arsenicals then the model produces the black curve in the right panel of Figure 5 for MMAs + DMAs. This curve, which is monotone increasing because the activation of AS3MT increases with GSH concentration and there is no sequestration effect, clearly does not fit the data. These results show conclusively that both effects, activation of AS3MT and sequestration of arsenicals by GSH are necessary to explain the experimental data. Similar results hold for the individual curves for MMAs and DMAs (simulations not shown).

Figure 5
figure 5

Removing influences of GSH on methylation. The left and right panels both reproduce the experimental data, [20] red and [19] green, and blue model curve for MMAs + DMAs with varying amounts of GSH from the right panel of Figure 4. The black curve in the left panel shows what the model fit would be if we removed from the model the excitatory influence of GSH on AS3MT. The black curve in the right panel shows what our model fit would be if we removed the binding of arsenicals to GSH from the model. Clearly, neither black curve fits the data. Both influences of GSH (blue curve) are necessary to explain the experimental data.

Temporal variation of GSH

In hepatocytes, GSH has concentrations in the mM range but is exported rapidly and turns over with a half-life of 1.5 to 2.5 hours [48]. The solutions in which cells are maintained typically contain the amino acids (cysteine, glycine, and glutamate or glutamine) necessary for the cells to resynthesize GSH. Nevertheless, the GSH concentration may vary considerably. For example, in the human hepatocytes used in [36], the cellular GSH concentration increased by 80% from day 1 to day 7. The experiments in [20] and [19] that we discussed above were conducted with purified enzymes in solution and not with living cells in vitro. The half-life of GSH in solution was found to vary from.2 to 70 hours depending on pH and temperature [49]. This raises the question of whether GSH degradation might play a role in experiments with purified enzyme in solution.

We investigated the effect of GSH degradation by a model simulation in which the half-life of GSH was assumed to be 2.5 hours. The initial amount of GSH in the reaction mixture was 20 mM and the initial amount of iAs was 1 μM. Figure 6 shows the time courses of the rate of first methylation step (V 1), the rate of the second methylation step (V 2), and the total rate of methylation as a function of time over 8 hours. In the beginning, methylation is very slow because there is so much GSH in the mixture and the GSH sequesters the arsenicals. Later, around 4–5 hours, there is much less sequestering and the activation of AS3MT causes methylation to proceed rapidly. Finally, at 7–8 hours, much of the GSH has degraded and so the methylation reaction runs slowly because of the lack of activation of AS3MT and the fact that many of the arsenicals have already been methylated. The main point is that, in both in vitro experiments and experiments in solution, the amount of GSH in the solution, the medium, and the cells should be tracked in time course experiments, because if it varies considerably that would affect the interpretation of the results.

Figure 6
figure 6

The effect of the half-life of GSH. The reactions mixture starts with 20 mM GSH and 1 μM iAs. We assume the half-life of GSH in the mixture is 2.5 hours. The rates of the first methylation step, V 1, the second methylation step, V 2, and total methylation vary dramatically throughout an 8 hour period. For explanations, see the text.


The main point of this study was to explore the different ways that GSH could affect the Challenger pathway [6] for oxidative methylation of inorganic arsenic. Three effects were included: (i) reduction of arsenicals with valence 5 to valence 3; (ii) activation of AS3MT; (iii) sequestration of arsenicals by binding to GSH. We used the model to analyze the experimental data in [20] and [19]. First we showed that experiments with MMAs as a substrate in [19] show clearly that MMAs III is an inhibitor of the second methylation step. Next we showed that the model predictions, Figure 3, match well the experiments in [19] where the amounts of iAs, MMAs, and DMAs were followed over time. Both [20] and [19] show that methylation proceeds slowly at low GSH and high GSH, but quite quickly at intermediate GSH ranges. This important finding is reproduced by the model, Figure 4, and we show that the reason for this is the combined effect of AS3MT activation by GSH and the sequestration of arsenicals by GSH. Finally, we pointed out that temporal variation in the amount GSH in reaction mixtures or cells needs to be taken into account in interpreting experimental data.

An important consequence of these findings is that recent experimental data can be explained well by the Challenger pathway augmented with these effects of GSH. This does not prove that the Hayakawa pathway [21], in which only GSH-conjugated arsenicals are methylated, is wrong. It just shows that the methylation of GSH-conjugated arsenicals is not necessary to explain the effects of GSH seen in [20] and [19]. Indeed, it is possible that both GSH free arsenicals and GSH bound arsenicals can be methylated, perhaps at different rates. There is some evidence that for a methyl transferase that is orthologous to AS3MT that GSH-conjugated arsenicals are preferred substrates for binding to the enzyme’s active site [23].

Easterling et al. [38] created a pharmacokinetic model to study the hepatocyte data in [36]. In order to fit the data, they needed to introduce a storage compartment for arsenicals in cells. Likewise, in our whole body model and hepatocyte model [30] we needed to introduce cellular storage compartments. It is tempting to speculate that the binding of arsenicals to GSH was an important part of the“storage mechanism” in both cases.

S-adenosylmethionine (SAM) is the methyl group donor in the methylation reactions. It is not included explicitly in our model because SAM was not varied in the experiments that we were trying to explain. The SAM concentration occurs implicitly in the V max values of the first and second methylation reactions. The K m of AS3MT for SAM was measured as 11.8 μM in [41], but the data in [20] imply that the K m is 50 μM. This is an important issue for the applications of arsenic biochemistry to human toxicity studies. Gamble and coworkers [28, 29] showed that folate supplementation of folate-deficient individuals in Bangladesh lowers blood arsenic levels. Raising folate levels can raise SAM concentrations in folate deficient individuals [50], so the presumed mechanism was that SAM levels were raised, thus making more methyl groups available for the methylation reactions. However, once SAM levels are back into the normal range (50 – 100 μM for rats), raising SAM more by further folate supplementation won’t help if the K m = 11.8 because the reaction will already be saturated, whereas if the K m = 50 μM then further supplementation should help.

The binding of GSH to arsenicals may be a significant detoxification mechanism as there is evidence that arsenic binds to GSH and then is removed in the bile [25, 51] and sequestration might also reduce the toxicity of trivalent arsenicals. Thus, whole body models of arsenic detoxification need to take into account this removal mechanism as well as the removal of arsenic-GSH conjugates from the liver to the blood and urine. This will be the subject of future work.

The effects of GSH on arsenic methylation discussed in this study and the removal of arsenic-GSH complexes in the bile and urine imply that increasing GSH might be a way to reduce As toxicity. GSH levels are under strong regulatory control in the liver [52]. Nevertheless, supplementation strategies have proven useful in several circumstances where GSH liver levels are low. N-acetyl cysteine is the antidote given in emergency departments in cases of acetaminophen overdose [53, 54] and glutamine is often given after surgery or other trauma to decrease inflammation [55, 56]. In both cases the intent is to increase GSH production in the liver. Plasma GSH levels in Bangladesh are quite low, 2.6 μM [57] as compared to the normal range, 2 – 20 μM [47, 58, 59]. This suggests that supplementation by N-acetyl-cysteine may be a viable strategy for reducing arsenic toxicity.


  •  The Challenger pathway, supplemented by three effects of glutathione, is sufficient to explain recent data on arsenic methylation.

  •  Monomethylarsonous acid inhibits the second methylation step.

  •  The three different effects of glutathione on arsenic methylation make the interpretation of experimental results difficult.

  •  Mathematical modeling of arsenic methylation can aid in the interpretation of experimental data.

  •  Supplementation by N-acetyl-cysteine may be a viable strategy for reducing arsenic toxicity.


  1. Ahmed MF, Ahuja S, Alauddin M, Hug SJ, Lloyd JR, Pfaff A, Pichler T, Saitikov C, Stute M, vanGeen A: Ensuring safe drinking water in Bangladesh . Science. 2006, 314: 1687-1688. 10.1126/science.1133146.

    Article  CAS  PubMed  Google Scholar 

  2. Kinniburgh DG, Smedley PL: Arsenic contamination of groundwater in Bangladesh. Final report . Tech. Rep. WC/00/19, British Geological Survey, Keyworth, UK, 2001,

  3. Jackson BP, Karagas VFTMR, Punshon T, Cottingham KL: Arsenic, organic foods, and brown rice syrup . Environ Health Perspect. 2012, 120: 623-626. 10.1289/ehp.1104619.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Wasserman GA, Liu X, Parvez F, Ahsan H, Factor-Litvak P, vanGeen A, Slavkovic V, Lolacano NJ, Cheng Z, Iftikhar H, Momotaj H, Graziono JH: Water arsenic exposure and children’s intellectual function in Araihazar, Bangladesh . Environ Health Perspect. 2004, 112: 1329-1333. 10.1289/ehp.6964.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  5. Wasserman GA, Liu X, Parvez F, Ahsan H, Factor-Litvak P, Kline J, vanGeen A, Slavkovic V, Lolacano NJ, Levy D, Cheng Z, Graziono JH: Water arsenic exposure and intellectual function in 6-year-old children in Araihazar, Bangladesh . Environ Health Perspect. 2007, 115: 285-289.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  6. Challenger F: Biological methylation . Adv Enzymol Biochem. 1951, 12: 429-491.

    CAS  Google Scholar 

  7. Styblo M, Yamauchi H, Thomas DJ: Comparative methylation of trivalent and pentavalent arsenicals . Toxicol Appl Pharmacol. 1995, 135: 172-178. 10.1006/taap.1995.1220.

    Article  CAS  PubMed  Google Scholar 

  8. Vahter M: Methylation of inorganic arsenic in different mammalian species and population groups . Sci Prog. 1999, 82: 69-88.

    CAS  PubMed  Google Scholar 

  9. Petrick JS, Ayala-Fierro F, Cullen WR, Carter DE, Aposhian HV: Monomethylarsonous acid (MMA(III)) is more toxic than arsenite in Chang human hepatocytes . Toxicol Appl Pharmacol. 2000, 163: 203-207. 10.1006/taap.1999.8872.

    Article  CAS  PubMed  Google Scholar 

  10. Styblo M, DelRazo LM, Vega L, Germolec DR, LeCluyse EL, Hamilton GA, Reed W, Wang C, Cullen WR, Thomas DJ: Comparative toxicity of trivalent and pentavalent inorganic and methylated arsenicals in rat and human cells . Arch Toxicol. 2000, 74: 289-299. 10.1007/s002040000134.

    Article  CAS  PubMed  Google Scholar 

  11. Cohen SM, Arnold LL, Uzvolgyi E, Cano M, John MS, Yamamoto S, Lu X, Li XC: Possible role of dimethylarsinous acid in dimethylarsinic acid-induceD urothelial toxicity and regeneration in the rat . Chem Res Toxicol. 2002, 15: 1150-1157. 10.1021/tx020026z.

    Article  CAS  PubMed  Google Scholar 

  12. Cullen WR, McBride BC, Reglinski J: The reaction of methylarsenicals with thiols: some biological implications . J Inorganic Chem. 1984, 21: 179-194.

    CAS  Google Scholar 

  13. Delnomdedieu M, Basti MM, Otvos JD, Thomas DJ: Reduction and binding of arsenate and dimethylarsinate by glutathione: a magnetic resonance study . Chemico-Biological Inter. 1994, 90: 139-155. 10.1016/0009-2797(94)90099-X.

    Article  CAS  Google Scholar 

  14. Delnomdedieu M, Basti MM, Styblo M, Otvos JD, Thomas DJ: Complexation of arsenic species in rabbit erythrocytes . Chem Res Toxicol. 1994, 7: 621-627. 10.1021/tx00041a006.

    Article  CAS  PubMed  Google Scholar 

  15. Sampayo-Reyes A, Zakharyan RA, Healy SM, Aposhian HV: Monomethylarsonic acid reductase and monomethylarsonous acid in hamster tissue . Chem Res Toxicol. 2000, 13: 1181-1186. 10.1021/tx000154s.

    Article  CAS  PubMed  Google Scholar 

  16. Zakharyan RA, Sampayo-Reyes A, Healy SM, Tsaprailis G, Board PG, Liebler DC, Aposhian HV: Human mono-methyl arsonic acid ( M M A s V ) reductase is a member of the glutathione S-transferase superfamily . Chem Res Toxicol. 2001, 14: 1051-1057. 10.1021/tx010052h.

    Article  CAS  PubMed  Google Scholar 

  17. Waters SB, Devesa V, Razo LMD, Styblo M, Thomas DJ: Endogenous reductants support the catalytic function of recombinant rat cyt19, an arsenic methyltransferase . Chem Res Toxicol. 2004, 17: 404-409. 10.1021/tx0342161.

    Article  CAS  PubMed  Google Scholar 

  18. Waters SB, Devesa V, Fricke MW, Creed JT, Styblo M, Thomas DJ: Glutathione modulates recombinant rat arsenic (+3 oxidation state) methyltransferase-catalyzed formation of trimethylarsine oxide and trimethylarsine . Chem Res Toxicol. 2004, 17: 1621-1629. 10.1021/tx0497853.

    Article  CAS  PubMed  Google Scholar 

  19. Ding L, Saunders RJ, Drobna Z, Walton FS, Xun P, Thomas DJ, Styblo M: Methylation of arsenic by recombinant human wild-type arsenic (+3 oxidation state) methyltransferase and its methionine 287 threonine (M287T) polymorph: Role of glutathione . Toxicol Appl Pharmacol. 2012, 264: 121-130. 10.1016/j.taap.2012.07.024.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  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-1406. 10.1016/j.biochi.2010.07.002.

    Article  CAS  PubMed  Google Scholar 

  21. Hayakawa T, Kobayashi Y, Cui X, Hirano S: A new metabolic pathway of arsenite: arsenic-glutathione complexes are substrates for human arsenic methyltransferase Cyt19 . Arch Toxicol. 2005, 79 (4): 183-191. 10.1007/s00204-004-0620-x.

    Article  CAS  PubMed  Google Scholar 

  22. Cullen WR: Chemical mechanism of arsenic biomethylation . Chem Res Toxicol. 2014, 27 (4): 457-461. 10.1021/tx400441h.

    Article  CAS  PubMed  Google Scholar 

  23. Marakapala K, Qin J, Rosen BP: Identification of Catalytic Residues in the As(III)S-Adenosylmethionine Methyltransferase . Biochemistry. 2012, 51: 944-951. 10.1021/bi201500c.

    Article  Google Scholar 

  24. Wang S, Li X, Song X, Geng Z, Hin X, Wang Z: Rapid equilibrium kinetic analysis of arsenite methylation catalyzed by recombinant human arsenic (+3 oxidation state) methyltransferase (hAS3MT) . J Biol Chem. 2012, 287: 38790-38799. 10.1074/jbc.M112.368050.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Kala SV, Neely MW, Kala G, Prater CI, Atwood DW, Rice JS, Liebermann MW: The MRP2/cMOAT transporter and arsenic-glutathione complex formation are required for biliary excretion of arsenic . J Biol Chem. 2000, 275: 33404-33408. 10.1074/jbc.M007030200.

    Article  CAS  PubMed  Google Scholar 

  26. Buchet JP, Lauwerys R, Roels H: Comparison of the urinary excretion of arsenic metabolites after a single oral dose of sodium arsenite, monomethylarsonate, or dimethylarsinate in man . Occup Environ Health. 1981, 48: 71-79. 10.1007/BF00405933.

    Article  CAS  Google Scholar 

  27. Buchet JP, Lauwerys R, Roels H: Urinary excretion of inorganic arsenic and its metabolites after repeated ingestion of sodium metaarsenite by volunteers . Occup Environ Health. 1981, 48: 111-118. 10.1007/BF00378431.

    Article  CAS  Google Scholar 

  28. Gamble MV, Liu X, Ahsan H, Pilsner JR, Ilievski V, Slavkovich V, Parvez F, Chen Y, Levy D, Factor-Litvak P, Graziano JH: Folate and arsenic metabolism: a double-blind, placebo-controlled folic acid-supplementation trial in Bangladesh . Am J Clin Nutr. 2006, 84: 1093-1101.

    PubMed Central  CAS  PubMed  Google Scholar 

  29. Gamble MV, Liu X, Slavkovich V, Pilsner JR, Ilievski V, Factor-Litvak P, Levy D, Alam S, Islam M, Parvez F, Ahsan H, Graziano JH: Folic acid supplementation lowers blood arsenic . TAm J Clin Nutr. 2007, 86: 1202-09.

    CAS  Google Scholar 

  30. Lawley SD, Cinderella M, Hall MH, Gamble MV, Nijhout HF, Reed MC: Mathematical model insights into arsenic detoxicfication . Theor Biol Med Model. 2011, 8: 1-17. 10.1186/1742-4682-8-1.

    Article  Google Scholar 

  31. Mann S, Droz PO, Vahter M: A physiologically based pharmacokinetic model for arsenic exposure I. Development in hamsters and rabbits . Toxicol Appl Pharamacol. 1996, 137: 8-22. 10.1006/taap.1996.0052.

    Article  CAS  Google Scholar 

  32. Mann S, Droz PO, Vahter M: A physiologically based pharmacokinetic model for arsenic exposure II. Humans . Toxicol Appl Pharamacol. 1996, 140: 471-486. 10.1006/taap.1996.0244.

    Article  CAS  Google Scholar 

  33. Gentry PR, Covington TR, Mann S, Shipp AM, Yager JW, Clewell HJ III: Physiologically based pharmocokinetic modeling of arsenic in the mouse . J Toxicol Environ Health. 2004, 67: 43-71. 10.1080/15287390490253660.

    Article  CAS  Google Scholar 

  34. Yu D: A physiologically based pharmacokinetic model of inorganic arsenic . Regul Toxicol Pharmacol. 1999, 29: 128-141. 10.1006/rtph.1999.1282.

    Article  CAS  PubMed  Google Scholar 

  35. El-Masri HA, Kenyon EM: Development of a human physiologically based pharmacokinetic (PBPK) model for inorganic arsenic and its mono- and di-methylated metabolites . J Pharmacokinet Pharmacodyn. 2008, 35: 31-68. 10.1007/s10928-007-9075-z.

    Article  CAS  PubMed  Google Scholar 

  36. Styblo M, Razo LMD, LeCluyse EL, Hamilton GA, Wang C, Cullen WR, Thomas DJ: Metabolism of arsenic in primary cultures of human and rat hepatocytes . Chem Res Toxicol. 1999, 12: 560-565. 10.1021/tx990050l.

    Article  CAS  PubMed  Google Scholar 

  37. Kenyon EM, Fea M, Styblo M, Evans MV: Application of modeling techniques to planning of in vitro arsenic pharmacokinetic studies . Alter Lab Animals. 2001, 29: 15-33.

    CAS  Google Scholar 

  38. Easterling MR, Styblo M, Evans MV, Kenyon EM: Pharmacokinetic modeling of arsenite uptake and metabolism in hepatocytes- mechanistic insights and implications for further experiments . J Pharmacokin Pharmacodyn. 2002, 29 (3): 207-234. 10.1023/A:1020248922689.

    Article  CAS  Google Scholar 

  39. Stamatelos SK, Brinkerhof CJ, Georgeopoulos SSIPG: Mathematical model of uptake and metabolism of arsenic(III) in human hepatocytes - Incorporation of cellular antioxidant response and thresholddependent behavior . BMC Syst Biol. 2011, 5: 1-15. 10.1186/1752-0509-5-1.

    Article  Google Scholar 

  40. Stamatelos SK, Androulakis IP, Kong AT, Georopoulos PG: A semi-mechanistic integrated toxicokinetic–toxicodynamic (TK/TD) model for arsenic(III) in hepatocytes . J Theor Biol. 2013, 317: 244-256.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. Wood TC, Salavigionne OE, Mukherjee B, Wang L, Klumpp AF, Thomae BA, Eckloff BW, Schaid DJ, Wieben ED, Weinshilboum RM: Human arsenic methyltransferase (AS3MT) pharmacogenetics . J Biol Chem. 2006, 281 (11): 7364-7373. 10.1074/jbc.M512227200.

    Article  CAS  PubMed  Google Scholar 

  42. Kedderis GL, Elmore AR, Crecelius EA, Yager JW, Goldsworthy TL: Kinetics of arsenic methylation by freshly isolated B6C3F1 hepatocytes . Chem-Biol Interact. 2006, 161: 139-145. 10.1016/j.cbi.2006.04.001.

    Article  CAS  PubMed  Google Scholar 

  43. Zakharyan RA, Ayala-Fierro F, Cullen WR, Carter DM, Aposhian HV: Enzymatic methylation of arsenic compunds. VII. Monomethylaraonous acid (MMAIII) Is the substrate for MMA Methyltransferase of rabbit liver and human hepatocytes . Toxicol Appl Pharmacol. 1999, 158: 9-15. 10.1006/taap.1999.8687.

    Article  CAS  PubMed  Google Scholar 

  44. Buchet JP, Lauwerys R: Role of thiols in the in vitro methylation of inorganic arsenic . Biochem Pharmacol. 1988, 37: 3149-3153. 10.1016/0006-2952(88)90313-9.

    Article  CAS  PubMed  Google Scholar 

  45. Scott N, Hatlelid KM, MacKenzie NE, Carter DE: Reactions of arsenic(III) and arsenic(V) species with glutathione . Chem Res Toxicol. 1993, 6: 102-106. 10.1021/tx00031a016.

    Article  CAS  PubMed  Google Scholar 

  46. Inoue M, Kinne R, Trans T, Arias I: Glutathione transport across hepatocyte plasma membranes . Eur J Biochem. 1984, 138: 491-495. 10.1111/j.1432-1033.1984.tb07943.x.

    Article  CAS  PubMed  Google Scholar 

  47. Wu G, Fang Y, Lupton JR, Turner ND: Glutathione metabolism and its implications for health . J Nutr. 2004, 134: 489-492.

    CAS  PubMed  Google Scholar 

  48. Baudouin-Cornu P, Lagniel G, Kumar C, Huang ME, Labarre J: Glutathione degradation is a key determinant of glutathione homeostasis . J Biol Chem. 2012, 287: 4552-4561. 10.1074/jbc.M111.315705.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Stevens R, Stevens L, Price NC: The stabilities of various thiol compounds used in protein purifications . Biochem Ed. 1983, 11 (2): 70-10.1016/0307-4412(83)90048-1.

    Article  CAS  Google Scholar 

  50. Reed M, Nijhout HF, Neuhouser M, Shane B, James SJ, Boynton A, Ulrich CM, III JFG: A mathematical model gives insights into nutritional and genetic aspects of folate-mediated one-carbon metabolism . J Nutr. 2006, 136: 2653-2661.

    CAS  PubMed  Google Scholar 

  51. Gyurasics A, Varga F, Gregus Z: Glutathione-dependent biliary excretion of arsenic . Biochem Pharmacol. 1991, 42: 465-468. 10.1016/0006-2952(91)90306-P.

    Article  CAS  PubMed  Google Scholar 

  52. 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-10.1186/1742-4682-5-8.

    Article  PubMed Central  PubMed  Google Scholar 

  53. Daly FF, Fountain JS, Murray L, Graudins A, Buckley NA: Guidelines for the management of paracetamol poinsoning in Australia and New Zealand - explanation and elaboration . Med J Australia. 2008, 188: 296-301.

    PubMed  Google Scholar 

  54. Ben-shachar R, Chen Y, Luo S, Hartman C, Reed M, Nijhout HF: The biochemistry of acetaminophen hepatotoxicity and rescue . Theor Biol Med Model. 2012, 9: 55-10.1186/1742-4682-9-55.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  55. Wilmore D: The effect of glutamine supplementation in patients following elective surgery and accidental injury . J Nutr. 2001, 131: 2543S-2549S.

    CAS  PubMed  Google Scholar 

  56. Roth E, Oehler R, Manhart N, Exner R, Wessner B, Strasser E, Spittler A: Regulative potential of glutathione - relation to glutathione metabolism . Nutrition. 2002, 18: 217-221. 10.1016/S0899-9007(01)00797-3.

    Article  CAS  PubMed  Google Scholar 

  57. Hall MN, Niedzwiecki M, Liu X, Harper KN, Alam S, Slavowich V, Ilievski V, Levy D, Siddique AB, Parvez F, Mey JL, Geen AV, Graziano J, Gamble MV: Chronic arsenic exposure and blood glutathione and glutathione disulfide concentrations in Bangladeshi adults . Environ Health Per. 2013, 121: 1068-1074.

    Google Scholar 

  58. Pogribna M, Melnyk S, Pogribny I, Chango A, Yi P, James SJ: Homocysteine metabolism in children with down syndrome: In Vitro modulation . Am J Hum Genet. 2001, 69: 88-95. 10.1086/321262.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. James SJ, Cutler P, Melnyk S, Jernigan S, Janak L, Gaylor DW, Neubrander JA: Metabolic biomarkers of increased oxidative stress and impaired methylation capacity in children with autism . Am J Clin Nutr. 2004, 80: 1611-1617.

    CAS  PubMed  Google Scholar 

Download references


This research was partially supported by NSF EF-1038593 (HFN,MR), and NIH grants R01 ES019876 (D. Thomas), R01 ES011601 (MVG), RO1 CA133595, and R00ES018890 (MNH).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Michael C Reed.

Additional information

Competing interests

The authors declare that they have no competing interest.

Authors’ contributions

MG and MH initiated the project and gave advice on the epidemiology and biochemistry of arsenic. FN directed the project. JY and SL developed the mathematical model and wrote the code with the advice of MR and FN. MR wrote the manuscript with the advice of FN, MG, MH, and SL. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lawley, S.D., Yun, J., Gamble, M.V. et al. Mathematical modeling of the effects of glutathione on arsenic methylation. Theor Biol Med Model 11, 20 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: