Mathematical modeling of the effects of glutathione on arsenic methylation
Theoretical Biology and Medical Modelling volume 11, Article number: 20 (2014)
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 [1–3]. 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 [6–8]. 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 [9–11].
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 [12–14] and glutathione S-transferase has been shown to help convert MMAs V to MMAs III in different tissues [15, 16]. Thomas, Styblo and colleagues [17–19] 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.  suggested that GSH increases the activity of AS3MT. In addition, in the experiments of both  and , 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.  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 . The model is depicted in Figure 1; details are given in Methods. We use the model to explain and interpret the experimental data in  and . 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 .
Cullen  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  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 , 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, 26–29]. There are cell culture experiments in which arsenicals are typically measured in the external medium . And, there are experiments in which reaction mixtures of arsenicals, AS3MT, and various other metabolites are prepared [17–21]. 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  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 [31–35]. We used a reduced version of our whole body model to study the cell culture experiments in . 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 . 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.
The velocity, V 1, of the reaction in which iAsIII is methylated to become MMAsV is given by:
We use the value K m = 4.6 μ M for the Michaelis-Menten constant for AS3MT for iAsIII as found in . The reaction has substrate inhibition by iAsIII; we take the inhibition constant to be M as found in . It is known that this reaction is inhibited by the product MMAsV and we take the inhibition constant, M from  and . We note that it is not certain that the enzyme investigated in  is identical to AS3MT.
The velocity, V 2, of the reaction in which MMAsIII is methylated to become DMAsV is given by:
As above we take K m = 4.6 μ M and we set M as in  and . The inhibition of V 2 by MMAsIII is proposed in this paper; the inhibition constant M was obtained by fitting the data in . 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 , 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 , 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, 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 , 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 . 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 . 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 , 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.
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 , 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.
In one set of experiments in , 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  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  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.
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 (, Figure two) and the Wang group (, 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].
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  and  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).
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 . 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 , the cellular GSH concentration increased by 80% from day 1 to day 7. The experiments in  and  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 . 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.
The main point of this study was to explore the different ways that GSH could affect the Challenger pathway  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  and . First we showed that experiments with MMAs as a substrate in  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  where the amounts of iAs, MMAs, and DMAs were followed over time. Both  and  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 , 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  and . 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 .
Easterling et al.  created a pharmacokinetic model to study the hepatocyte data in . 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  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 , but the data in  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 , 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 . 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  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.
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.
Kinniburgh DG, Smedley PL: Arsenic contamination of groundwater in Bangladesh. Final report . Tech. Rep. WC/00/19, British Geological Survey, Keyworth, UK, 2001,
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.
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.
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.
Challenger F: Biological methylation . Adv Enzymol Biochem. 1951, 12: 429-491.
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.
Vahter M: Methylation of inorganic arsenic in different mammalian species and population groups . Sci Prog. 1999, 82: 69-88.
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.
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.
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.
Cullen WR, McBride BC, Reglinski J: The reaction of methylarsenicals with thiols: some biological implications . J Inorganic Chem. 1984, 21: 179-194.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Cullen WR: Chemical mechanism of arsenic biomethylation . Chem Res Toxicol. 2014, 27 (4): 457-461. 10.1021/tx400441h.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Yu D: A physiologically based pharmacokinetic model of inorganic arsenic . Regul Toxicol Pharmacol. 1999, 29: 128-141. 10.1006/rtph.1999.1282.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Wu G, Fang Y, Lupton JR, Turner ND: Glutathione metabolism and its implications for health . J Nutr. 2004, 134: 489-492.
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.
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.
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.
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.
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.
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.
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.
Wilmore D: The effect of glutamine supplementation in patients following elective surgery and accidental injury . J Nutr. 2001, 131: 2543S-2549S.
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.
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.
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.
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.
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).
The authors declare that they have no competing interest.
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
About this article
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). https://doi.org/10.1186/1742-4682-11-20