Skip to main content

A mathematical model of glutathione metabolism



Glutathione (GSH) plays an important role in anti-oxidant defense and detoxification reactions. It is primarily synthesized in the liver by the transsulfuration pathway and exported to provide precursors for in situ GSH synthesis by other tissues. Deficits in glutathione have been implicated in aging and a host of diseases including Alzheimer's disease, Parkinson's disease, cardiovascular disease, cancer, Down syndrome and autism.


We explore the properties of glutathione metabolism in the liver by experimenting with a mathematical model of one-carbon metabolism, the transsulfuration pathway, and glutathione synthesis, transport, and breakdown. The model is based on known properties of the enzymes and the regulation of those enzymes by oxidative stress. We explore the half-life of glutathione, the regulation of glutathione synthesis, and its sensitivity to fluctuations in amino acid input. We use the model to simulate the metabolic profiles previously observed in Down syndrome and autism and compare the model results to clinical data.


We show that the glutathione pools in hepatic cells and in the blood are quite insensitive to fluctuations in amino acid input and offer an explanation based on model predictions. In contrast, we show that hepatic glutathione pools are highly sensitive to the level of oxidative stress. The model shows that overexpression of genes on chromosome 21 and an increase in oxidative stress can explain the metabolic profile of Down syndrome. The model also correctly simulates the metabolic profile of autism when oxidative stress is substantially increased and the adenosine concentration is raised. Finally, we discuss how individual variation arises and its consequences for one-carbon and glutathione metabolism.


Glutathione is a low molecular weight tri-peptide (γ-glutamyl-cysteinyl-glycine) found at relatively high concentrations (0.5–10 mM) in all mammalian cells and relatively low concentrations (2–20 μM) in plasma [1]. Inside cells, most of the glutathione (85–90%) is in the cytosol where it primarily exists in a reduced form (GSH) and to a much lesser extent as an oxidized disulfide form (GSSG). The high GSH/GSSG ratio provides the essential reducing environment inside the cell. GSH is manufactured in the cytosol by a two-step process: the first step, which combines cysteine and glutamate, is catalyzed by γ-glutamylcysteine synthetase (GCS); the second step, which adds the glycine residue, is catalyzed by glutathione synthetase (GS). Glycine and glutamate are produced and used by many metabolic reactions and have relatively high cytosolic concentrations. Cytosolic cysteine is the limiting amino acid for GSH synthesis because it has a low concentration compared to glycine and glutamate. Cytosolic cysteine comes from only three sources: (1) from methionine via the methionine cycle and the transsulfuration pathway, (2) direct import into the cell from the plasma, and (3) from excess protein catabolism over protein synthesis. Thus the availability of cysteine and the activity of GCS are the major determinants of GSH synthesis. The enzyme cystathionine-β-synthase (CBS) that catalyzes the first step in the transsulfuration pathway is highly expressed in liver cells but not highly expressed in peripheral cells, so it is not surprising that the liver is the major producer of GSH, much of which is exported to the plasma and enzymatically broken down to cysteinylglycine and cyst(e)ine that is subsequently taken up by other cells for GSH synthesis.

Glutathione is involved in many pathways that are essential for normal intracellular homeostasis. It detoxifies xenobiotics and heavy metals through a reaction catalyzed by GSH S-transferases that bind them to the sulfhydryl group on the cysteine residue. GSH plays a role in regulating lipid, glucose, and amino acid metabolism because it is necessary for the hepatic response to insulin-sensitizing agents [2]. GSH is necessary for the interconversion of prostaglandins [3]. The removal of formaldehyde, a carcinogen and a product of one-carbon metabolism, requires glutathione, and glutathione is involved in T-lymphocyte activation and viral resistance [4]. Finally, glutathione scavenges reactive oxygen species including superoxide and hydrogen peroxide. In these reactions GSH is oxidized to GSSG and the ratio [GSH]/[GSSG], an indicator of the redox status of the cell, is known to regulate redox sensitive enzymes in the pathways for cell proliferation and cell apoptosis [5]. Thus, it is not surprising that GSH (or the [GSH]/[GSSG] ratio) plays a key role in many diseases including cancer, inflammation, Alzheimer's disease, Parkinson's disease, sickle cell anemia, liver disease, cystic fibrosis, AIDS, heart attack, stroke, and diabetes [4, 6] as well as in aging [79]. Reactive oxygen species also cause birth defects in rats, which are prevented by administration of GSH [10]. For more on glutathione chemistry and health effects, see [1, 4, 1115].

During the past several years we have created mathematical models for different parts of one-carbon metabolism [1621]. The purpose of the modeling was to answer questions posed by experiments or experimentalists and to investigate mechanisms of regulation in one-carbon metabolism. In this paper, we extend our most recent model [20] to include cysteine and glutathione metabolism (Figure 1). Since this mathematical model is quite complicated, it is useful to be clear why our model needs to include all of one carbon metabolism and not just the transsulfuration pathway. First, methionine is a major hepatic source of cysteine through the methionine cycle and the transsulfuration pathway. Secondly, the redox status of the cell affects many of the enzymes in one-carbon metabolism including MATI, MATIII, MS, BHMT, as well as CBS and GCS in the transsulfuration pathway, and therefore one cannot evaluate GSH metabolism without including the methionine and folate cycles. Thirdly, patients with Down syndrome or autism have increased oxidative stress and exhibit particular disturbed profiles of one-carbon metabolism [1315]. We would like to understand how oxidative stress (and the chromosome 21 trisomy in the case of Down syndrome) could create these disturbed profiles.

Figure 1
figure 1

One-carbon metabolism and the transsulfuration pathway. Rectangles enclose the names or acronyms of substrates that are variables in the model. Substrates not in rectangles are held constant or are the products of reactions that we do not keep track of. Arrows at the bottom of the figure represent import from the gut and other cells, and losses to other cells and to degradation. There is one differential equation for each substrate. The ellipses contain the acronyms of the enzymes that catalyze the reactions. Full names for all the enzymes and substrates, as well as a complete description of the mathematical model and the values of all parameters are given in the online Additional File 1.

Model Overview

Figure 1 shows the biochemical pathways in the hepatic cellular model used in this paper. Rectangular boxes represent the substrates that can vary in the model, and the ellipses contain the acronyms of the enzymes that catalyze particular reactions. There is one differential equation for each substrate that says that the rate of change of the concentration of the substrate is the sum of the reaction velocities (μM/hr) that produce it minus the sum of the reaction velocities that use it. Full names for all the enzymes and substrates are given in Additional File 1. Non-boxed substrates are taken to be constant or indicate products of reactions. This model extends the model for one carbon metabolism in [20] by adding the transsulfuration pathway, the synthesis of glutathione, and its transport into the blood. Here we discuss the main ideas involved in modeling the transsulfuration pathway, referring the reader to [20] (and its online supplementary material) for a discussion of the other parts of the model. A complete description of the full mathematical model and the values of all parameters are given in the online Additional File 1.

One of the most interesting features of one carbon metabolism is that reaction velocities are often affected not just by substrate and product concentrations and enzyme activities but also by the concentrations of substrates in distant parts of the reaction network that act as allosteric activators or inhibitors of the enzyme. As a result, the formulas for reaction velocities are often complicated functions of many variables. This is well illustrated by the velocity equation for the CBS reaction, the first step in the transsulfuration pathway:

V C B S = ( V m a x [ H c y ] [ S e r ] ( K m h c y + [ H c y ] ) ( K m s e r + [ S e r ] ) ) ( ( 1.2 ) C ( [ S A M ] + [ S A H ] ) 2 30 2 + ( [ S A M ] + [ S A H ] ) 2 ) ( k a + [ H 2 O 2 ] k a + [ H 2 O 2 ] n o r m ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiqaaaqaaiabdAfawnaaBaaaleaacqWGdbWqcqWGcbGqcqWGtbWuaeqaaOGaeyypa0ZaaeWaaKqbagaadaWcaaqaaiabdAfawnaaBaaabaGaemyBa0MaemyyaeMaemiEaGhabeaacqGGBbWwcqWGibascqWGJbWycqWG5bqEcqGGDbqxcqGGBbWwcqWGtbWucqWGLbqzcqWGYbGCcqGGDbqxaeaacqGGOaakcqWGlbWsdaqhaaqaaiabd2gaTbqaaiabdIgaOjabdogaJjabdMha5baacqGHRaWkcqGGBbWwcqWGibascqWGJbWycqWG5bqEcqGGDbqxcqGGPaqkcqGGOaakcqWGlbWsdaqhaaqaaiabd2gaTbqaaiabdohaZjabdwgaLjabdkhaYbaacqGHRaWkcqGGBbWwcqWGtbWucqWGLbqzcqWGYbGCcqGGDbqxcqGGPaqkaaaakiaawIcacaGLPaaaaeaacqGHflY1daqadaqcfayaamaalaaabaGaeiikaGIaeGymaeJaeiOla4IaeGOmaiJaeiykaKIaem4qamKaeiikaGIaei4waSLaem4uamLaemyqaeKaemyta0Kaeiyxa0Laey4kaSIaei4waSLaem4uamLaemyqaeKaemisaGKaeiyxa0LaeiykaKYaaWbaaeqabaGaeGOmaidaaaqaaiabiodaZiabicdaWmaaCaaabeqaaiabikdaYaaacqGHRaWkcqGGOaakcqGGBbWwcqWGtbWucqWGbbqqcqWGnbqtcqGGDbqxcqGHRaWkcqGGBbWwcqWGtbWucqWGbbqqcqWGibascqGGDbqxcqGGPaqkdaahaaqabeaacqaIYaGmaaaaaaGccaGLOaGaayzkaaWaaeWaaKqbagaadaWcaaqaaiabdUgaRnaaBaaabaGaemyyaegabeaacqGHRaWkcqGGBbWwcqWGibasdaWgaaqaaiabikdaYaqabaGaem4ta80aaSbaaeaacqaIYaGmaeqaaiabc2faDbqaaiabdUgaRnaaBaaabaGaemyyaegabeaacqGHRaWkcqGGBbWwcqWGibasdaWgaaqaaiabikdaYaqabaGaem4ta80aaSbaaeaacqaIYaGmaeqaaiabc2faDjaaxcW7daWgaaqaaiabd6gaUjabd+gaVjabdkhaYjabd2gaTbqabaaaaaGccaGLOaGaayzkaaaaaaaa@B125@

The first factor is simply Michaelis-Menten kinetics for the CBS reaction that uses homocysteine (Hcy) and serine (Ser) as substrates. We take the Michaelis constants from the literature. The second term is the activation of CBS by SAM and SAH that was discovered by Finkelstein and Martin [22]. The form of the activation was derived by nonlinear regression on the data in [23, 24]. The constant C is chosen so that the second term equals one in the normal steady state (see below). The last term is the activation of CBS by oxidative stress [25, 26], which is represented by the concentration of H2O2.

The differential equation for the concentration of cytosolic cysteine (Cys) is straightforward:

d d t [ C y s ] = V G C S ( [ C y s ] , [ G l u ] , [ G S H ] , [ H 2 O 2 ] ) + V C G T L ( [ C y s t a ] ) + V c y s i n ( [ b C y s ] ) ( ( 0.35 ) [ C y s ] 2 200 ) MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiqaaaqaaKqbaoaalaaabaGaemizaqgabaGaemizaqMaemiDaqhaaOGaei4waSLaem4qamKaemyEaKNaem4CamNaeiyxa0Laeyypa0JaeyOeI0IaemOvay1aaSbaaSqaaiabdEeahjabdoeadjabdofatbqabaGccqGGOaakcqGGBbWwcqWGdbWqcqWG5bqEcqWGZbWCcqGGDbqxcqGGSaalcqGGBbWwcqWGhbWrcqWGSbaBcqWG1bqDcqGGDbqxcqGGSaalcqGGBbWwcqWGhbWrcqWGtbWucqWGibascqGGDbqxcqGGSaalcqGGBbWwcqWGibasdaWgaaWcbaGaeGOmaidabeaakiabd+eapnaaBaaaleaacqaIYaGmaeqaaOGaeiyxa0LaeiykaKcabaGaey4kaSIaemOvay1aaSbaaSqaaiabdoeadjabdEeahjabdsfaujabdYeambqabaGccqGGOaakcqGGBbWwcqWGdbWqcqWG5bqEcqWGZbWCcqWG0baDcqWGHbqycqGGDbqxcqGGPaqkcqGHRaWkcqWGwbGvdaWgaaWcbaGaem4yamMaemyEaKNaem4CamNaemyAaKMaemOBa4gabeaakiabcIcaOiabcUfaBjabdkgaIjabdoeadjabdMha5jabdohaZjabc2faDjabcMcaPiabgkHiTmaabmaajuaGbaWaaSaaaeaacqGGOaakcqaIWaamcqGGUaGlcqaIZaWmcqaI1aqncqGGPaqkcqGGBbWwcqWGdbWqcqWG5bqEcqWGZbWCcqGGDbqxdaahaaqabeaacqaIYaGmaaaabaGaeGOmaiJaeGimaaJaeGimaadaaaGccaGLOaGaayzkaaaaaaaa@938E@

The first term on the right is the production of cysteine from cystathionine (Cysta) by CGTL and the second term is the import of cysteine into the cell, which depends on the concentration of cysteine in the blood ([bCys]). The third term is the loss of cysteine in the reaction catalyzed by GCS that makes glutamyl-cysteine and the fourth term represents the loss of cysteine to other pathways (for example to sulfate and taurine). Cysteine also used for protein synthesis and is produced by protein catabolism; in the model we assume that these two rates balance. The form for the fourth term was chosen because the data in [27] indicate that at normal cysteine concentrations (approximately 200 μM) most of the flux away from cysteine is toward GSH and only a moderate amount towards other pathways. However, as [Cys] rises an increasing fraction is sent towards the synthesis of taurine. Formulas for VCGTL and VbCYSc (the transport of cysteine into the cell from the blood) appear in Additional File 1.

The first step in the synthesis of GSH is the formation of γ-glutamyl-cysteine (GluCys) from the constituent amino acids glutamate and cysteine by the enzyme γ-glutamyl-cysteine synthetase (GCS, also called glutamate cysteine ligase (GCL)). The reaction is reversible and GSH is a competitive inhibitor of GCS against glutamate [2830].

The third factor in the following formula is the activation of GCS by H2O2 [25, 26].

V G C S = ( k a + [ H 2 O 2 ] k a + [ H 2 O 2 ] s s ) ( V m a x ( [ G l u ] [ C y s ] [ G l u C y s ] / K e ) K m c y s K m g l u + K m c y s [ G l u ] + K m g l u [ G y s ] ( 1 + [ G S H ] K i + [ G l u ] K m g l u ) + [ G l u C y s ] K p + [ G S H ] K i ) ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeaabiqaaaqaaiabdAfawnaaBaaaleaacqWGhbWrcqWGdbWqcqWGtbWuaeqaaOGaaGjbVlabg2da9maabmaabaWaaSaaaeaacqWGRbWAdaWgaaWcbaGaemyyaegabeaakiabgUcaRiabcUfaBjabdIeainaaBaaaleaacqaIYaGmaeqaaOGaem4ta80aaSbaaSqaaiabikdaYaqabaGccqGGDbqxaeaacqWGRbWAdaWgaaWcbaGaemyyaegabeaakiabgUcaRiabcUfaBjabdIeainaaBaaaleaacqaIYaGmaeqaaOGaem4ta80aaSbaaSqaaiabikdaYaqabaGccqGGDbqxdaWgaaWcbaGaem4CamNaem4CamhabeaaaaaakiaawIcacaGLPaaacqGHflY1aeaacaaMe8+aaeWaaKqbagaadaWcaaqaaiabdAfawnaaBaaabaGaemyBa0MaemyyaeMaemiEaGhabeaacqGGOaakcqGGBbWwcqWGhbWrcqWGSbaBcqWG1bqDcqGGDbqxcqGGBbWwcqWGdbWqcqWG5bqEcqWGZbWCcqGGDbqxcqGHsislcqGGBbWwcqWGhbWrcqWGSbaBcqWG1bqDcqWGdbWqcqWG5bqEcqWGZbWCcqGGDbqxcqGHVaWlcqWGlbWsdaWgaaqaaiabdwgaLbqabaGaeiykaKcabaGaem4saS0aa0baaeaacqWGTbqBaeaacqWGJbWycqWG5bqEcqWGZbWCaaGaem4saS0aa0baaeaacqWGTbqBaeaacqWGNbWzcqWGSbaBcqWG1bqDaaGaey4kaSIaem4saS0aa0baaeaacqWGTbqBaeaacqWGJbWycqWG5bqEcqWGZbWCaaGaei4waSLaem4raCKaemiBaWMaemyDauNaeiyxa0Laey4kaSIaem4saS0aa0baaeaacqWGTbqBaeaacqWGNbWzcqWGSbaBcqWG1bqDaaGaei4waSLaem4raCKaemyEaKNaem4CamNaeiyxa0LaeiikaGIaeGymaeJaey4kaSYaaSqaaeaacqGGBbWwcqWGhbWrcqWGtbWucqWGibascqGGDbqxaeaacqWGlbWsdaWgaaqaaiabdMgaPbqabaaaaiabgUcaRmaaleaabaGaei4waSLaem4raCKaemiBaWMaemyDauNaeiyxa0fabaGaem4saS0aa0baaeaacqWGTbqBaeaacqWGNbWzcqWGSbaBcqWG1bqDaaaaaiabcMcaPiabgUcaRmaaleaabaGaei4waSLaem4raCKaemiBaWMaemyDauNaem4qamKaemyEaKNaem4CamNaeiyxa0fabaGaem4saS0aaSbaaeaacqWGWbaCaeqaaaaacqGHRaWkdaWcbaqaaiabcUfaBjabdEeahjabdofatjabdIeaijabc2faDbqaaiabdUealnaaBaaabaGaemyAaKgabeaaaaGaeiykaKcaaaGccaGLOaGaayzkaaGaeyOla4caaaaa@D42F@

The second step in the synthesis of GSH is the addition of the glycine residue to GluCys by the enzyme glutathione synthase (GS). We follow [30, 31] and use a reversible bi-reactant Michaelis-Menten mechanism.

V G S = ( V m a x ( [ G l u C y s ] [ G l y ] [ G S H ] / K e ) K m g l u c y s K m g l y + K m g l y [ G l u C y s ] + K m g l u c y s [ G l y ] ( 1 + [ G l u C y s ] K m g l u c y s ) + [ G S H ] K p ) ) . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaGaemOvay1aaSbaaSqaaiabdEeahjabdofatbqabaGccaaMe8Uaeyypa0JaaGjbVpaabmaajuaGbaWaaSaaaeaacqWGwbGvdaWgaaqaaiabd2gaTjabdggaHjabdIha4bqabaGaeiikaGIaei4waSLaem4raCKaemiBaWMaemyDauNaem4qamKaemyEaKNaem4CamNaeiyxa0Laei4waSLaem4raCKaemiBaWMaemyEaKNaeiyxa0LaeyOeI0Iaei4waSLaem4raCKaem4uamLaemisaGKaeiyxa0Laey4la8Iaem4saS0aaSbaaeaacqWGLbqzaeqaaiabcMcaPaqaaiabdUealnaaDaaabaGaemyBa0gabaGaem4zaCMaemiBaWMaemyDauNaem4yamMaemyEaKNaem4CamhaaiabdUealnaaDaaabaGaemyBa0gabaGaem4zaCMaemiBaWMaemyEaKhaaiabgUcaRiabdUealnaaDaaabaGaemyBa0gabaGaem4zaCMaemiBaWMaemyEaKhaaiabcUfaBjabdEeahjabdYgaSjabdwha1jabdoeadjabdMha5jabdohaZjabc2faDjabgUcaRiabdUealnaaDaaabaGaemyBa0gabaGaem4zaCMaemiBaWMaemyDauNaem4yamMaemyEaKNaem4CamhaaiabcUfaBjabdEeahjabdYgaSjabdMha5jabc2faDjabcIcaOiabigdaXiabgUcaRmaaleaabaGaei4waSLaem4raCKaemiBaWMaemyDauNaem4qamKaemyEaKNaem4CamNaeiyxa0fabaGaem4saS0aa0baaeaacqWGTbqBaeaacqWGNbWzcqWGSbaBcqWG1bqDcqWGJbWycqWG5bqEcqWGZbWCaaaaaiabcMcaPiabgUcaRmaaleaabaGaei4waSLaem4raCKaem4uamLaemisaGKaeiyxa0fabaGaem4saS0aaSbaaeaacqWGWbaCaeqaaaaacqGGPaqkaaaakiaawIcacaGLPaaacqGHUaGlaaa@B30B@

The differential equation for cytosolic GSH is:

d d t [ G S H ] = V G S ( [ G l y ] , [ G l u C y s ] V g s h _ o u t ( [ G S H ] ) 2 V G P X ( [ G S H ] , [ H 2 O 2 ] ) + 2 V G R ( [ G S S G ] , [ N A D P H ] ) ( 0.002 ) [ G S H ] . MathType@MTEF@5@5@+=feaafiart1ev1aaatCvAUfKttLearuWrP9MDH5MBPbIqV92AaeXatLxBI9gBaebbnrfifHhDYfgasaacPC6xNi=xI8qiVKYPFjYdHaVhbbf9v8qqaqFr0xc9vqFj0dXdbba91qpepeI8k8fiI+fsY=rqGqVepae9pg0db9vqaiVgFr0xfr=xfr=xc9adbaqaaeGaciGaaiaabeqaaeqabiWaaaGcbaqbaeWabmqaaaqaaKqbaoaalaaabaGaemizaqgabaGaemizaqMaemiDaqhaaOGaei4waSLaem4raCKaem4uamLaemisaGKaeiyxa0Laeyypa0JaemOvay1aaSbaaSqaaiabdEeahjabdofatbqabaGccqGGOaakcqGGBbWwcqWGhbWrcqWGSbaBcqWG5bqEcqGGDbqxcqGGSaalcqGGBbWwcqWGhbWrcqWGSbaBcqWG1bqDcqWGdbWqcqWG5bqEcqWGZbWCcqGGDbqxcqGHsislcqWGwbGvdaWgaaWcbaGaem4zaCMaem4CamNaemiAaGMaei4xa8Laem4Ba8MaemyDauNaemiDaqhabeaakiabcIcaOiabcUfaBjabdEeahjabdofatjabdIeaijabc2faDjabcMcaPaqaaiabgkHiTiabikdaYiabdAfawnaaBaaaleaacqWGhbWrcqWGqbaucqWGybawaeqaaOGaeiikaGIaei4waSLaem4raCKaem4uamLaemisaGKaeiyxa0LaeiilaWIaei4waSLaemisaG0aaSbaaSqaaiabikdaYaqabaGccqWGpbWtdaWgaaWcbaGaeGOmaidabeaakiabc2faDjabcMcaPaqaaiabgUcaRiabikdaYiabdAfawnaaBaaaleaacqWGhbWrcqWGsbGuaeqaaOGaeiikaGIaei4waSLaem4raCKaem4uamLaem4uamLaem4raCKaeiyxa0LaeiilaWIaei4waSLaemOta4KaemyqaeKaemiraqKaemiuaaLaemisaGKaeiyxa0LaeiykaKIaeyOeI0IaeiikaGIaeGimaaJaeiOla4IaeGimaaJaeGimaaJaeGOmaiJaeiykaKIaei4waSLaem4raCKaem4uamLaemisaGKaeiyxa0LaeiOla4caaaaa@9D17@

The first term is the synthesis of GSH from glycine and glutamyl-cysteine. The second term is the transport of GSH out of the liver cell into the blood. Vgsh_out is actually the sum of two terms, one for the high affinity transporter and one for the low affinity transporter [32]. We ignore the canalicular transport into the bile because it is a relatively small percentage of total export [33]. The third term is the rate of production of oxidized GSSG from GSH via the enzyme GPX and the fourth term is the conversion of GSSG back to GSH via the enzyme GR. We ignore other reactive oxygen species besides H2O2, or, put a different way, H2O2 represents them all. The number two occurs in both these terms because two molecules of GSH combine to make one GSSG. In the fifth term we are assuming that 0.2% of the GSH is removed each hour in detoxification reactions that form conjugates [12].

The kinetics of sinusoidal efflux of GSH has been well studied in the perfused rat liver. The major part of the flux is carried by the low affinity transporter, which has sigmoidal kinetics with a Vmax in the range 900–1400 μM/hr, a Km of approximately 3000 μM, and a Hill coefficient of approximately 3 [3436]. In our model we use a Vmax of 1100 μM/hr, a Km of 3000 μM, and a Hill coefficient of 3. We use standard Michaelis-Menten kinetics for the high affinity GSH transporter and for the two GSSG transporters.

We track five variables in the blood, [GSH], [GSSG], [Gly], [Cys] and [Glut]. Glycine, glutamate, and cysteine enter blood from intestinal absorption at rates that we vary in various experiments with the model; the normal rates are 630 μM/hr, 273 μM/hr and 70 μM/hr, respectively. We assume for convenience that the volume of the blood is the same as the volume of the liver. Glycine, cysteine, and glutamate leave the blood by transport into liver cells (depending on their concentrations) and they are also formed in the blood by the breakdown of GSH and GSSG into their component amino acids [37, 38]. We also assume that normally 10% of the cysteine, glycine, and glutamate, in the blood is taken up per hour by other cells and that an additional 25% of cysteine is converted to cystine. Under normal conditions a large percentage of blood GSH and GSSG is broken down into the component amino acids and a small amount is taken up by other cells or otherwise leaves the system. As above, full details and formulas appear in Addition File 1.

For each in silico computation, the values of various constants (like H2O2) are given, as are the methionine and serine levels in the blood, and the rates of input of cysteine glutamate, and glycine into the blood. These are the "inputs" to the model. The differential equations are then solved to determine the steady-state values of the concentrations of all the variables and the steady state rates of all the reactions. Of course, if the inputs are different the steady state will be different. We experiment with the model by changing the inputs or changing parameters (for example, a parameter that gives the strength of a particular allosteric interaction) and determine what the effect is. By removing interactions we can take the model apart piece by piece so that we can understand how and why glutathione metabolism works the way it does. We also allow the inputs to vary as functions of time (for example the amino acid input will vary because of meals) and compute the time course of each concentration and reaction rate. This allows us to investigate the homeostatic mechanisms that protect the system against fluctuations in the inputs.

A number of substrate concentrations are fixed in the model and in all the simulations reported below. These include: cytosolic GAR (10), NADPH (50), betaine (50), formaldehyde (500), dUMP (20), and total cellular folate (20). All concentrations are in μM.

Limitations of the model

This model was designed to allow us to study various regulatory mechanisms in the transsulfuration pathway and the effects of oxidative stress, particularly as applied to Down syndrome and autism. No mathematical model can track all of the variables that might affect a complex biochemical system such as glutathione metabolism. This is also true, of course, in biological experimentation. This model is no exception. We ignore canalicular excretion of GSH. We use Km values in the ranges determined experimentally but there is much less information on Vmax values. Often we choose Vmax values so that the steady state concentrations of substrates and products lie within the normal published ranges. Cellular amino acid concentrations are increased by feeding and protein degradation and decreased by protein synthesis, growth and use in one-carbon metabolism. In this model we assume that protein synthesis and degradation are in balance and that no amino acids are used for growth. The consequences of this assumption are outlined in the discussion.

One-carbon metabolism and the transsulfuration pathway contain many allosteric interactions by which substrates in one part of the pathway affect the activity of distant enzymes. We use experimentally determined forms for these allosteric interactions but sometimes the details of the kinetics are not known, forcing us to make reasonable educated guesses. Similarly, many effects of oxidative stress on the enzymes of one carbon metabolism and the transsulfuration pathways are known but detailed kinetics are not available.

In this paper we are mainly interested in intracellular liver metabolism, so we take a somewhat simple view of the fates glutathione and its metabolites in the blood. Future work will include a more detailed model of the blood compartment and inter-organ regulation of glutathione and its component amino acids. Thus, we do not expect that our model will make perfect quantitative predictions. Rather, we want to use it to investigate the qualitative features of glutathione metabolism in the normal state and in various disease states.


A. Normal model steady-state concentrations and velocities

We take the normal values of inputs to be the following. Blood methionine is 30 μM and blood serine is 150 μM. The rates of cysteine, glycine, and glutamate input to the blood are 70 μM/hr, 630 μM/hr, and 273 μM/hr respectively. The normal concentration of H2O2 is 0.01 μM. With these inputs, the model computes the concentrations of the cytosolic variables given in Table 1.

Table 1 Normal model cytosolic concentrations (μM)

These model values correspond well to the values in the literature. For cytosolic folate variables and methionine cycle variables, see the discussions in our previous papers [1621]. Typical values of GSH in animal cells are in the range 500–10,000 μM or 0.5–10 mM [1, 39]. Typical values for cysteine are 150–250 μM [1] and for cystathionine 40 μM [40]. The ratio [GSH]/[GSSG] is thought to be around 100 for cells that are not under oxidative stress [41], and [GSH]/[GSSG] = 107.5 in our model cell.

The computed velocities of the cytosolic reactions are given in Table 2. There is very little information in the literature about reaction velocities because they are difficult to measure. However, the model concentration of GSH declines in the fasting state about as rapidly as observed experimentally (See Section B, below). This indicates that the overall rates of GSH production from cysteine and methionine and the transport of GSH out of the cell are in the appropriate ranges. We also note that the flux around the methionine cycle is 205 μM/hr and approximately half enters the transsulfuration pathway (VCBS = 103 μM/hr) and half is remethylated to methionine in accordance with the results of Finkelstein and Martin [22].

Table 2 Normal model cytosolic reaction velocities (μM/hr)

The computed concentrations of variables in the blood are given in Table 3. Wu et al. report that the combined cysteine and cystine concentrations are 110–325 μM [1]. In our model the computed plasma cysteine concentration is 186 μM, which is in the middle of this range. Plasma concentrations in humans are reported in the range 2–20 μM for GSH [1, 13, 14], and 0.14–0.34 μM for GSSG [13, 14, 42]. The average plasma GSH/GSSG ratio is reported to be in the range 25–28 μM with a large standard deviation [14, 15], and in the model it is 26.5. Plasma glycine levels are reported to be approximately 300 μM in [43].

Table 3 Normal model blood concentrations (μM)

The computed values of various transport rates are given in Table 4. We use the abbreviations o = outside, b = blood, c = cytosol, so, for example, VoCysb is the transport of cysteine from the outside into the blood. VoCysb, VoGlyb, and VoGlutb are inputs to the model. All other transport velocities are computed by the model. The second row shows the transport velocities of the five amino acids in the model from the blood into liver cells. The third row shows the transport velocities of GSH and GSSG from the cell into the blood. Detailed kinetic information is available on amino acid transporters [44, 45] and on the high and low affinity transporters of GSH and GSSG [32, 39, 46] and we chose our kinetics parameters from this literature.

Table 4 Normal model net transport velocities (μM/hr)

The fourth row in Table 4 requires more comment. Our main interest is to understand the synthesis and export of GSH in liver cells and how intracellular metabolite balance is affected by oxidative stress. Since GSH is exported rapidly from liver cells and much of the export is broken down into the constituent amino acids that are then reimported into liver cells, it was necessary to include the blood compartment in our model. The blood communicates with all other tissues none of which are in our model. We have therefore necessarily made a number of assumptions about the loss of GSH, GSSG, Cys, Gly, and Glu to other tissues. For example, as discussed above, we assume that normally 10% per hour of the cysteine, glycine, and glutamate in the blood is taken up by other cells and that an additional 25% of cysteine in the blood is lost by conversion to cystine. The velocities in the fourth row reflect these assumptions.

B. The Half-life of Glutathione

Ookhtens et al. [34] reported that when buthionine sulfoximine is used to inhibit the activity of GCS (which catalyzes the first step in GSH synthesis) a half life of 2–6 hours for cellular GSH is observed. This is consistent with the experiments of [47]. Moreover, the rate of sinusoidal GSH efflux in both fed and starved rats is near saturation at about 80% of Vmax, about 1000–1200 μM/h [34]. Thus, if the cytosolic GSH concentration is approximately 7000 μM, then the half life would be in the 2–3 hour range. Therefore, a variety of experimental studies and calculations consistently suggest a short half life in the 2–3 hour range.

By contrast, Aw et al. [33] report that rats fasted for 48 hours lose approximately 44% of the intracellular GSH in their hepatocytes. They also report that after 48 hours the rate of GSH transport out of the cell declined by 38%. These results are consistent with Tateishi et al. [48, 49] who reported a decline in liver GSH to a level between one half and two thirds of normal after a 48 hour fast. These experiments suggest a half-life longer than two days. One possible explanation for this long half-life under starved conditions is that the normal dietary amino acid input is partly replaced by protein catabolism. However, given the normal rate of GSH efflux, a 48 hour half-life would require that catabolism replace 94% of daily dietary input, which seems improbably high.

An alternative explanation, which could potentially explain both sets of experiments, is that exported GSH is broken down into constituent amino acids in the blood that are rapidly reimported into the liver cells. Indeed, it is known that the enzyme γ-glutamyltranspeptidase (GGT) on the external cell membrane initiates this process (called the γ-glutamyl cycle) [12, 50, 51]. In our model the computed value of GSH transport out of the cell (Table 4) is VcGSHb = 1152 and the rates of Cys, Gly, and Glut import are also high (Table 4), although we assume that 10% per hour of the amino acids in the blood are lost to non-liver cells and an additional 25% of Cys is lost by conversion to cystine. Figure 2 shows the cytosolic concentration of GSH in our model liver cells for 10 hours after the concentration of the enzyme GCS was set to zero. The computed half life of GSH is 3 hours.

Figure 2
figure 2

GSH half life after GCS is blocked. When GSH synthesis is stopped the model intracellular GSH concentration declines rapidly with a half life of approximately 3 hours.

Figure 3 shows the concentration of GSH and other metabolites in our model liver cell during a fasting experiment over a 48 hour period. We assume that during fasting, protein catabolism supplies 1/3 of the normal amino acid input. The GSH concentration declines slowly over the 48 hour period to about 50% of normal and the rate of GSH export declines to 67% of normal consistent with the experiments reported in [33]. Thus the rapid reimport hypothesis explains both sets of data. Other metabolites show interesting changes during the fast. The methionine cycle metabolites adjust very rapidly to the decreased methionine input reaching new steady states within a few hours. However, the metabolites in the GSH synthesis, export and reimport pathway decline very slowly, achieving their new steady states in 4–5 days (data not shown).

Figure 3
figure 3

GSH and GSH transport under fasting conditions. After three hours, the inputs of cysteine, methionine, glycine, glutamate, and serine are reduced to 1/3 of normal. The intracellular GSH + GSSG concentration declines slowly over the 48 hour period to about 50% of normal and the rate of GSH export declines to 67% of normal consistent with the experiments reported in [33]. The cytosolic and blood cysteine concentrations decline proportionally to GSH. The methionine cycle metabolites and fluxes equilibrate rapidly.

Mosharov et al. [26] studied the role of the transsulfuration pathway in GSH synthesis. When they blocked CTGL they observed that that the intracellular GSH concentration dropped to a new steady state of approximately 42% of control in about 24 hours. They concluded that approximately half of GSH is derived via the transsulfuration pathway. Our methionine input is computed to be 103 μM/h consistent with experimental measurements [26, 52] and our cysteine input into the system is 70 μM/h. If we remove the methionine input in the model, the GSH concentration declines to a new steady-state 47% of normal. On the other hand, if we remove cysteine input the GSH concentration declines to a new steady-state 68% of normal. These model results support the interpretation in [26] that methionine and cysteine inputs contribute equivalently to GSH synthesis. We remark that one would not expect the contributions of the two metabolites to GSH synthesis to be strictly proportional to their inputs since they are used in other reactions and the system is highly non-linear.

C. Inhibition of GCS

It is well known [27, 41] that GSH is a competitive inhibitor against glutamate of GCS, the enzyme that catalyzes the synthesis of glutamyl-cysteine. This inhibition can naturally be thought of as product inhibition, one step removed. As glutathione rises it indirectly inhibits its own synthesis and as glutathione falls the inhibition is released. Figure 4 shows that this inhibition has the effect that is expected at steady state. As sulfur amino acid input rises, so does glutathione concentration but not as fast as it would if the inhibition were not present. At low sulfur amino acid concentrations the effect is small. Thus the primary effect of the inhibition is to prevent excess accumulation of GSH. Such accumulation would sequester more amino acids, would increase transport out of the cell up to saturation, and would therefore increase the loss of cysteine to cystine in the blood.

Figure 4
figure 4

Effect of GSH feedback inhibition. The GSH concentration is plotted as a function of the intracellular concentration of sulfur amino acids. The solid curve shows the GSH concentration when the inhibition of GCS is included in the model. The dashed curve shows the GSH concentration when the inhibition by GSH is removed.

D. Stability of GSH under large input fluctuations

Hepatocytes are seldom at steady state [21] because they receive large protein inputs during and shortly after meals and relatively little protein input between meals. How do these daily fluctuations affect the intracellular and blood glutathione pools? To investigate this question, we varied the amino acid input (cysteine, glycine, glutamate, methionine, and serine) throughout the day and used the model to compute the time courses of all the concentrations in the model, and all the velocities. Let A denote the daily average input of a particular amino acid. While fasting (for example, from 12 midnight until 7 am), we assume that the input is (0.25)A. From 7 am to 10 am, we assume that the input is (1.75)A corresponding to breakfast, followed by two hours of fasting. Then, from 12 noon until 3 pm we assume that the input is (1.75)A corresponding to lunch, followed by three hours of fasting and then an input of (3.25)A for three hours corresponding to dinner. The complete daily input profile is shown in Figure 5A.

Figure 5
figure 5

The Stability of the glutathione pools in the face of input fluctuations. Panel A shows the amino acid input to the model hepatocytes throughout a 24 hour day. The input is 25% of the mean while fasting, 175% of the mean for three hours after breakfast and lunch, and 325% of the mean for three hours after dinner. For discussion, see the text. Panel B shows moderate variations in methionine concentrations and extremely large swings in SAM concentration, but GSH concentration remains stable. Panel C shows that the velocity of the CBS reaction varies dramatically, but the velocity of the GS reaction, which synthesizes glutathione, shows milder variation. As expected, the long range allosteric reactions between the folate cycle and the methionine cycle stabilize the velocity of the DNA methylation reaction (vDNMT). Panel D shows that there are large variations in 5mTHF and Hcy throughout the day, but the GSH concentration in the blood remains stable.

Panel B of Figure 5 shows that the fluctuations in methionine input cause very large fluctuations in SAM but only modest fluctuations in intracellular methionine concentration, as suggested by the early methionine loading experiments of Finkelstein and co-workers [53, 54]. Corrales et al. [55] have suggested that the relative stability of methionine is due to the complicated kinetics of MAT-I and MAT-III and we have confirmed this by model experiments [56]. Panel B also shows that the intracellular GSH concentration (divided by 100 for graphing purposes) also remains stable throughout the day.

Panel C of Figure 5 shows that the velocity of the CBS reaction tracks the methionine input as expected, but that the velocity of the GS reaction by which GSH is synthesized has milder fluctuations. These smaller fluctuations are a result of the inhibition of the GCS reaction by GSH (see Section C). In addition, the intracellular GSH pool is normally very large (6591 μM in our model, see Table 1). Both of these effects contribute to the impressive stability of the GSH concentration, seen in Panel B, in the face of large fluctuations in amino acid input. Note that the velocity of the DNA methylation reaction (DNMT in Panel C) is also quite stable despite the fact that its substrate, SAM, is undergoing very large fluctuations. It is understood that this is a result of long-range allosteric interactions between the methionine cycle and the folate cycle [18].

Panel D of Figure 5 shows that homocysteine undergoes large fluctuations in synchrony with the methionine input as expected. 5mTHF also fluctuates but in the opposite direction from homocysteine for two reasons. First, when methionine input rises dramatically, so does [SAM] and SAM inhibits MTHFR. Secondly, when homocysteine rises, it drives the MS reaction faster, which reduces 5mTHF. Finally, the blood concentration of GSH remains completely stable despite the large transient amino acid input fluctuations.

E. Oxidative Stress

In our model, oxidative stress is represented by the concentration of H2O2. When H2O2 increases there are several effects on one-carbon metabolism. First, the increased concentration of H2O2 inhibits the enzymes MS and BHMT and activates the enzymes CBS and GCS [25, 26]. Secondly, the balance of GSH and GSSG is shifted toward GSSG via the GPX and GR reactions. This affects the upstream metabolites in the methionine and transsulfuration pathways because GSSG inhibits the enzymes MAT-I and MAT-III [57, 58]. All of these influences are in the model; for details, see Additional File 1. The response to oxidative stress in the model is surprisingly complex; see Figure 6.

Figure 6
figure 6

Effect of oxidative stress. The curves show the effect on the steady-state values of cysteine and GSH in the cytosol and the blood, the rate of GSH synthesis by GS, and the cytosolic [GSH]/[GSSG] ratio as H2O2 concentration is raised from normal (0.01 μM) to 0.05 μM.

Under moderate oxidative stress there are moderate increases on blood and cytosolic GSH and blood cysteine, while cytosolic cysteine and the [GSH]/[GSSG] ratio decline. Cytosolic GSH increases because oxidative stress activates CBS and GCS increasing the flux through the GCS and GS reactions and simultaneously lowering the cytosolic cysteine concentration. Since cytosolic GSH increases, the export out of the cell will also increase thus raising the blood GSH and blood cysteine concentrations. The elevated H2O2 concentration drives the balance in the GPX and GR reactions towards GSSG thus lowering the [GSH]/[GSSG] ratio. Under high oxidative stress this balance is shifted even further towards GSSG and this has consequences for overall cysteine balance.

In the model cytosolic GSH has three fates: it is transported into the blood, there is a net flux to GSSG, and 0.2% is removed per hour corresponding to detoxification reactions and excretion into the bile. Likewise, cytosolic GSSG has two fates: it is transported into the blood, and 10% is removed per hour corresponding to excretion into the bile. Of course, removal of one GSH or GSSG results in the removal of one or two cysteines, respectively. At normal steady state concentrations the cysteine lost by these two mechanisms are about equal. However, as the oxidative stress increases and the balance between GSH and GSSG shifts toward GSSG, more cysteines are lost from the system per hour. At moderate oxidative stress this effect small. However, with high or chronic levels of oxidative stress this effect gets much larger and the loss of cysteines is quite large. This causes the rate of the GS reaction to come back down to normal despite the upregulation of CBS and GCS and cause the steady concentrations of cytosolic GSH and the blood concentrations of GSH and cysteine to decline below normal; see Figure 6.

E. The Metabolic Profile of Down Syndrome

Down syndrome is a complex metabolic and genetic disorder whose root cause, trisomy 21, is an extra copy of chromosome 21 [59]. Down syndrome is not rare; it occurs in approximately 1 out every 700–800 live births [60]. Children with Down syndrome have abnormal metabolic profiles and show increased incidence of a large number of serious diseases including leukemia and diabetes [61]. In most cases, it is not understood whether these diseases are caused by the extra chromosome, the altered metabolic profile, or both.

To investigate the metabolite profile of Down syndrome using the model, we began by increasing by 50% the Vmax of CBS, since the gene for CBS is on chromosome 21 and is expressed at 150% of normal. The first column of Table 5 shows the average percent change in the levels of six plasma metabolites in 42 Down patients compared to controls (taken from [13]). The second column shows the percentage change in these metabolites in the model when the Vmax of CBS is increased by 50%. Note that the intracellular concentrations of Hcy, SAM, SAH, and Met all change in the same direction as seen clinically. We would not expect a close match to the clinically observed percentage changes because we are comparing intracellular model changes to blood measurements. The increased dosage of CBS has almost no effect on the model plasma concentrations of bCys and bGSH. Thus these changes must come from some other effect of chromosome 21 trisomy.

Table 5 Down syndrome: Clinical and model metabolic changes (%)*

It is known that Down patients suffer from mild to moderate oxidative stress due to the overexpression of the Cu-Zn superoxide dismutase (SOD) gene that is also located on chromosome 21 [62]. Column 3 in Table 5 shows the effects on metabolite concentrations when the H2O2 concentration (normally 0.01 μM) is increased to 0.025 μM in addition to the increased dosage of CBS. The methionine cycle metabolites are further reduced compounding the effects of trisomy 21. Blood cysteine increases substantially and blood GSH declines modestly. The reason for the increase in [bCys] under mild oxidative stress is discussed in Section E.

A number of clinical observations suggest that Down patients may have a functional folate deficiency despite having normal plasma levels of folate and vitamin B12 [13]. The model results support the discussion in [13] of why this is so. The increased expression of CBS lowers the concentration of Hcy and therefore lowers the rate of the MS reaction. Thus, folate builds up in the form of 5mTHF (the "methyl trap") and there is less folate in the forms CH2-THF and 10f-THF that are the substrates for thymidylate and purine synthesis, respectively. Indeed, the last three rows of Table 5 show that the up-regulation of CBS has exactly this effect and that the addition of oxidative stress makes the effect much stronger.

F. The Metabolic Profile of Autism

Autism is a complex neurodevelopmental disorder whose cause is unknown and which is diagnosed solely on behavioral criteria [63]. The incidence of autism has been increasing and is now thought to be present in ~1 in 150 children in the United States [64]. Numerous studies suggest that both genetic and environmental influences are important in the development of autism [6567]. Some studies have identified metabolic deficiencies in autistic children [14], and it has recently been shown that there is a characteristic metabolic profile in many autistic children involving disturbances in methionine and glutathione metabolism [14, 15]. The first and second columns in Table 6 show the changes in blood levels of various metabolites in autistic children compared to normal patients in two different clinical studies [14, 15].

Table 6 Autism: Clinical and model metabolic changes (%)*

One of the working hypotheses in [14, 15] is that autistic children have a self-perpetuating cycle of inflammation and oxidative stress. Indeed, it is reported that autistic children have increased levels of brain and gut inflammation [68, 69]. The fact that blood GSH levels in autistic children were down 32–46% and the blood GSH/GSSG ratio was down 48–66% is also consistent with this hypothesis (Table 6). Oxidative stress is represented in the model by the concentration of H2O2. The H2O2 concentration drives the redox balance in the cell from GSH toward GSSG (see Figure 1) and the increased concentrations of H2O2 and GSSG can modulate many enzymes in one-carbon metabolism (see Section E). In Column 3 of Table 6 we show the result of a substantial increase of oxidative stress in the model to 5 times the normal level. The model intracellular concentrations of Hcy, Met, and SAM all fall dramatically corresponding to the clinically observed declines in their plasma concentrations. The model concentration of bGSH and the bGSH/bGSSG ratio both decline in the plasma similarly to the degree of decline seen in the autistic patients in the clinic. Finally, we note that severe oxidative stress causes a "methyl trap" such that the model concentration of 5mTHF more than doubles. There are several underlying reasons for the methyl trap: the concentration of Hcy drops, oxidative stress inhibits MS, and the drop in [SAM] releases the inhibition of MTHFR, which drives more folate towards the 5mTHF form.

There is one apparent discrepancy between the model results in column 3 (Table 6) and the clinical results in columns 1 and 2. The model intracellular concentration of SAH declines under severe oxidative stress while in patients the average concentration of SAH in the blood went up 24–49%. It is reported in [14] that the plasma adenosine concentration was significantly elevated on the average in autistic children compared to controls, though both populations showed very large standard deviations. The observed plasma increase in adenosine could be due to a decrease in adenosine deaminase activity, inhibition of adenosine kinase, or an increase in 5-nucleotidase, all of which can occur in the presence of oxidative stress [7072]. In our model, intracellular adenosine concentration is a constant. When we raised the constant by a factor of 5, we obtained the results in column 4 (Table 6). The metabolic profile is very similar to that in column 3 except that now the intracellular SAH concentration does not decline but instead increases by 24% in line with the clinical data.

Finally, we note that the model predicts modest decreases in the rate of the DNA methylation reaction (vDNMT) and this is consistent with the reported global DNA hypomethylation seen in autistic children [73]. The decrease in reaction rate is caused by the decrease in [SAM] and the increase in [SAH], which inhibits DNMT.


We have expanded our previous model of one-carbon metabolism to include the transsulfuration pathway, the synthesis and export of glutathione, and its breakdown in and export from the blood. This has enabled us to investigate a number of specific problems about glutathione metabolism from the experimental and clinical literature.

We have shown that the glutathione pool is very insensitive to the large variations inacute amino acid input caused by daily meals. In some cases we have been able to give causal explanations for observed experimental phenomena and we are able to test the quantitative effects of various allosteric inhibitions and excitations. Finally, we used to model to show how the observed metabolic profile of Down syndrome could arise from trisomy 21 and moderate oxidative stress and how the observed metabolic profile of autism patients could arise from high oxidative stress and increased adenosine concentration.

In a series of papers [3436], Ookhtens and colleagues showed that the half-life of glutathione in the liver was quite low, approximately 2–3 hours, consistent with the result of Lauterburg [47]. On the other hand, the same group showed in [33] that glutathione has a half-life of approximately two days in the livers of fasted rats. In Section B we confirmed that there is no contradiction between these results. In the first case, GSH synthesis is blocked, so the rapid export of GSH makes the GSH concentration decline rapidly. In the second case, even though the rats are fasted, the rapid reuptake of cysteine, glycine, and glutamate by the liver cells insures that the synthesis of GSH declines relatively slowly and therefore the observed half-life is long. Finally, the model results support the conclusions of Mosharov et al. [26] that both cysteine and methionine contribute approximately equally to GSH synthesis in the liver. This is true even though GSH is exported rapidly and cysteine is reimported rapidly compared to the methionine input.

Lu proposes in [12] that the high glutathione concentration in hepatocytes is a storage mechanism for cysteine. But what is the reason for the rapid cycling, i.e. fast export of GSH, breakdown by GGT, and fast reimport of cysteine? This is a futile cycle that requires a lot of energy. A reasonable hypothesis is that the rapid cycling allows the liver to respond quickly to the glutathione requirements of other tissues. This hypothesis is consistent with the idea that glutathione is a mechanism for cysteine storage [12], but also helps explain the reason for the γ-glutamyl cycle and the reason for the short half-life of hepatic GSH.

Cell metabolism is very complicated and the same substrate is often used in many different reactions. As a result the response function of a metabolite or a reaction velocity to changes in a parameter or input may be nonlinear and non-monotone. For example, in Section E we showed that moderate oxidative stress causes blood GSH and blood cysteine to rise, but severe oxidative stress causes blood GSH and blood cysteine to fall. This increase at low oxidative stress is due to the stimulation of CBS and GCS that increases GSH synthesis and concentration, and therefore the rate of export. At high or chronic oxidative stress, however, the model suggests that the balance shifts towards GSSG, and removal of cysteine in the form of GSSG dominates, resulting in a decline in cysteine.

There is increasing evidence that oxidative stress plays a role in the development of autism [74]. The metabolic profile of autistic patients has been shown to be abnormal with elevated biomarkers that indicate chronic oxidative stress and evidence that GSH synthesis may be insufficient to maintain redox homeostasis [14, 15]. Likewise, the overexpression of SOD is children with Down syndrome leads to a reduction of GSH and an increase in oxidative stress [13]. In our model oxidative stress is represented by an elevated level of H2O2 which induces many changes in one-carbon metabolism and the transsulfuration pathway. H2O2 stimulates CBS and GCS and inhibits MS and BHMT. In addition H2O2 drives the GSH/GSSG balance towards GSSG, which inhibits MAT I and MAT III. We have found that, in our model, oxidative stress alone can produce some but not all the metabolic characteristics of Down syndrome and autism. However, the addition of trisomy 21 in the first case, and raised adenosine in the second, brings the profiles much closer to those observed in patients with Down syndrome and autism, respectively.

Cellular amino acid concentrations are increased by feeding and protein degradation and decreased by protein synthesis, growth and use in one-carbon metabolism. During early growth (until about the second year after birth), about 10–20% of the amino acid pool is used in growth and is thus not available for GSH synthesis and one-carbon metabolism [7577]. This would be expected to have an effect on the rates amino-acid-requiring processes of one-carbon metabolism and glutathione synthesis. We have found, by simulation, that if we reduce the amino acid input into the system by 15%, the concentration of GSH and the synthesis rate of GSH are proportionally diminished, but there is little effect on the DNA methylation reaction, while reactions in the folate cycle are reduced by 2–9%. This reduction in GSH synthesis may contribute to excessive oxidative stress in infants.

Calculations with the model show that blood concentrations do not necessarily reflect intracellular concentrations of metabolites. For example, the increased dosage of CBS and GCS in our simulation of Down syndrome causes the intracellular concentration of cysteine to decline while the blood concentration increases. This shows that care must be taken in interpreting blood measurements, and that ideally one would like to conduct experiments in which both intracellular and extracellular concentrations are measured. By contrast, we found in the model (simulations not shown) that the blood concentrations of GSH and GSSG track the intracellular concentrations.

The purpose of this model was to study the properties of intracellular glutathione metabolism, in particular the effects of oxidative stress and trisomy 21. Of course intracellular glutathione metabolism is affected by the import of amino acids and the export and removal of GSH and GSSG. We therefore needed include a blood compartment and to keep track of bCys, bGly, bGSH, bGSSG and bGlut. The kinetics of the contents of the blood compartment are complicated because GSH and its component amino acids in the blood interact with the kidney, the brain, and other tissues. These important interactions are beyond the scope of this initial investigation, but will be incorporated in future studies.

This model provides a tool, different from laboratory experimentation, that can be used to experiment with one-carbon and glutathione metabolism. By using the tool, it is possible to quickly and easily test hypotheses, evaluate new ideas, and give causal explanations for what is seen in the clinic and the lab. In this spirit, the authors welcome questions, suggestions, and hypotheses to test from clinicians and experimentalists.


Enzymes. AICART:

aminoimidazolecarboxamide ribonucleotide transferase


betaine-homocysteine methyltransferase


cystathionine β-synthase




dihydrofolate reductase


dimethylglycine dehydrogenase




10-formyltetrahydrofolate dehydrogenase


10-formyltetrahydrofolate synthase


γ-glutamylcysteine synthetase


glycine decarboxylase (glycine cleavage system)


glycine N-methyltransferase


glutathione peroxidase


glutathione reductase


glutathione synthetase


methionine adenosyl transferase I


methionine adenosyl transferase III


methionine synthase


5,10-methenyltetrahydrofolate cyclohydrolase


5,10-methylenetetrahydrofolate dehydrogenase


5,10-methylenetetrahydrofolate reductase


non-enzymatic conversion


Phosphoribosyl glycinamidetransformalase


S-adenosylhomocysteine hydrolase


sarcosine dehydrogenase




thymidylate synthase

Metabolites. 10f-THF:





P-ribosyl-5-amino-4-imidazole carboxamide






cysteine Cysta cystathionine






deoxythymidine monophosphate


deoxyuridine monophophate


glycinamide ribonucleotide




Cys glutamyl-cysteine


glycine. GSH: glutathione


glutathione disulfide

H2C= O:

formaldehyde. H2O2 hydrogen peroxide. HCOOH formate. Hcy homocysteine. Met methionine. NADPH nicotinamide adenine dinucleotide phosphate












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

    CAS  PubMed  Google Scholar 

  2. Guarino MP, Afonso RA, Raimundo N, Raposo JF, Macedo MP: Hepatic glutathione and nitric oxide are critical for hepatic insulin-sensitizing substance action. Am J Physiol Gastrointest Liver Physiol. 2003, 284 (4): G588-G594.

    Article  CAS  PubMed  Google Scholar 

  3. Lu SC: Regulation of glutathione synthesis. Curr Top Cell Reg. 2000, 36: 95-116.

    Article  CAS  Google Scholar 

  4. Townsend DM, Tew KD, Tapiero H: The importance of glutathione in human disease. Biomed Pharmacother. 2003, 57: 145-155. 10.1016/S0753-3322(03)00043-X.

    Article  CAS  PubMed  Google Scholar 

  5. Sen CK: Cellular thiols and redox-regulated signal transduction. Curr Top Cell Regul. 36: 1-30.

  6. Klaunig JE, Kamendulis LM: The role of oxidative stress in carcinogenesis. Ann Rev Pharm Toxicol. 2004, 44: 239-267. 10.1146/annurev.pharmtox.44.101802.121851.

    Article  CAS  PubMed  Google Scholar 

  7. Hernantz A, Fernandez-Vivancos E, Montiel C, Vazquez JJ, Arnalich S: Changes in the intracellular homocysteine and glutathione content associated with aging. Life Sciences. 2000, 67: 1317-1324. 10.1016/S0024-3205(00)00722-0.

    Article  Google Scholar 

  8. Miller R, Buehner G, Chang Y, Harper JM, Sigler R, Smith-Wheelock : Methionine deficient diet extends mouse lifespan, slows immune and lens aging, alters glucose, T4, IGF-1 and insulin levels and increases hepatocyte MIF levels and stress resistence. Aging Cell. 2005, 4: 119-10.1111/j.1474-9726.2005.00152.x.

    Article  CAS  PubMed  Google Scholar 

  9. Mockett R, Sohal RS, Orr WC: Overexpression of glutathione reductase extends survival in transgenic Drosophila melanogaster under hypoxia but not normoxia. FASEB J. 1999, 13: 1733-1742.

    CAS  PubMed  Google Scholar 

  10. Ishibashi M, Akazawa S, Sakamaki H, Matsumoto K, Yamasaki H, Yamaguchi Y, Goto S, Urata Y, Kondo T, Nagataki S: Oxygen-induced embryopathy and the significance of glutathione-dependent antioxidant system in the rat embryo during early organogenesis. Free Rad Biol Med. 1997, 22: 447-454. 10.1016/S0891-5849(96)00338-3.

    Article  CAS  PubMed  Google Scholar 

  11. Meister A: Glutathione deficiency produced by inhibition of its synthesis, and its reversal; applications in research and therapy. Pharmac Ther. 1991, 51: 155-194. 10.1016/0163-7258(91)90076-X.

    Article  CAS  Google Scholar 

  12. Lu SC: Regulation of hepatic glutathione synthesis: current concepts and controversies. FASEB J. 1999, 13: 1169-1183.

    CAS  PubMed  Google Scholar 

  13. 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 

  14. James SJ, Cutler P, Melnyk S, Jernigan S, Janak L, Gaylor D, 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 

  15. James SJ, Melnyk S, Jernigan S, Cleves MA, Halsted CH, Wong DH, Cutler P, Bock K, Boris M, Bradstreet JJ, Baker SM, Gaylor DW: Metabolic endophenotype and related genotypes are associated with oxidative stress in children with autism. Am J Med Gen. 2006, 141B: 1-10. 10.1002/ajmg.b.30026.

    Article  Google Scholar 

  16. Reed MC, Nijhout HF, Sparks R, Ulrich CM: A Mathematical model of the methionine cycle. J Theor Biol. 2004, 226: 33-43. 10.1016/j.jtbi.2003.08.001.

    Article  CAS  PubMed  Google Scholar 

  17. Nijhout HF, Reed MC, Budu P, Ulrich CM: A Mathematical Model of the Folate Cycle. J Biol Chem. 2004, 279: 55008-55016. 10.1074/jbc.M410818200.

    Article  CAS  PubMed  Google Scholar 

  18. Nijhout HF, Reed MC, Anderson DA, Mattingly J, James SJ, Ulrich CM: Long-Range Allosteric Interactions between the Folate and Methionine Cycles Stabilize DNA Methylation Reaction Rate. Epigenetics. 2006, 1: 81-87.

    Article  PubMed  Google Scholar 

  19. Reed MC, Nijhout HF, Neuhouser ML, Gregory JF, Shane B, James SJ, Boynton A, Ulrich CM: A Mathematical Model Gives Insights into Nutritional and Genetic Aspects of Folate-Mediated One-Carbon Metabolism. J Nutr. 2006, 136: 2653-2661.

    CAS  PubMed  Google Scholar 

  20. Nijhout HF, Reed MC, Lam S-L, Gregory JF, Shane B, Ulrich CM: In silico experimentation with a model of hepatic mitochondrial folate metabolism. Theor Biol Med Mod. 2006, 3: 40-10.1186/1742-4682-3-40.

    Article  Google Scholar 

  21. Nijhout HF, Reed MC, Lam S-L, Gregory JF, Shane B, Ulrich CM: A day in the life of cell metabolism. J Biol Theory. 2008, 2: 2124-127.

    Google Scholar 

  22. Finkelstein J, Martin J: Methionine metabolism in mammals: distribution of homocysteine between competing pathways. J Biol Chem. 1984, 259: 9508-9513.

    CAS  PubMed  Google Scholar 

  23. Janosik M, Kery V, Gaustadnes M, Maclean KMN, 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. 10.1021/bi010711p.

    Article  CAS  PubMed  Google Scholar 

  24. Kluijtmans LAJ, Boers GHJ, Stevens EMB, Renier WO, Kraus JP, Trijbels FJM, van den Heuvel LPWJ, Blom HJ: Defective cystathionine beta-synthase regulation by S-adenosyl-methionine in a partially pyridoxine responsive homocysteinuria patient. J Clin Invest. 1996, 98: 285-289. 10.1172/JCI118791.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  25. Deplancke B, Gaskins HR: Redox control of the transsulfuration and glutathione biosynthesis pathways. Curr Opin Clin Nutr Metab Care. 2002, 5: 85-92. 10.1097/00075197-200201000-00015.

    Article  CAS  PubMed  Google Scholar 

  26. Mosharov E, Cranford MR, Banerjee R: The quantitatively important relationship between homocysteine metabolism and glutathione synthesis by the transsulfuration pathway and its regulation by redox changes. Biochem. 2000, 39: 13005-13011. 10.1021/bi001088w.

    Article  CAS  Google Scholar 

  27. Stipanuk MH, Coloso RM, Garcia RAG, Banks MF: Cysteine concentration regulates cysteine metabolism to glutathione, sulfate, and taurine in rat hepatocytes. J Nutr. 1992, 122: 420-427.

    CAS  PubMed  Google Scholar 

  28. Richman PG, Meister A: Regulation of γ-glutamyl-cysteine synthetase by non-allosteric feedback inhibition by glutathione. J Biol Chem. 1975, 250: 1422-1426.

    CAS  PubMed  Google Scholar 

  29. Seelig G, Meister A: Glutathione biosynthesis; γ-glutamylcysteine synthetase from rat kidney. Methods in Enzymology. 1985, 113: 379-399.

    Article  CAS  PubMed  Google Scholar 

  30. Mendoza-Cozatl DG, Moreno-Sanchez R: Control of glutathione and phytochelatin synthesis under cadmium stress. Pathway modeling for plants. J Theor Biol. 2006, 238: 919-936. 10.1016/j.jtbi.2005.07.003.

    Article  CAS  PubMed  Google Scholar 

  31. Hell R, Bergmann L: γ-Glutamylcysteine synthestase in higher plants: catalytic properties and subcellular localization. Planta. 1990, 180: 603-612. 10.1007/BF02411460.

    Article  CAS  PubMed  Google Scholar 

  32. Ookhtens M, Kaplowitz N: Role of the liver in interorgan homeostasis of glutathione and cyst(e)ine. Sem Liver Dis. 1998, 18: 313-329.

    Article  CAS  Google Scholar 

  33. Aw TY, Ookhtens M, Ren C, Kaplowitz N: Kinetics of glutathione efflux from isolated rat hepatocytes. Am J Physiol Gast Liver. 1986, 250: G236-G243.

    CAS  Google Scholar 

  34. Ookhtens M, Hobdy K, Corvasce MC, Aw TY, Kaplowitz M: Sinusoidal efflux of glutathione in the perfused rat liver. J Clin Invest. 1985, 75: 258-265. 10.1172/JCI111682.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Ookhtens M, Maddatu T: Mechanism of changes in hepatic sinusoidal and biliary glutathione efflux with age in rats. Am J Physiol. 1991, 261: G648-G656.

    CAS  PubMed  Google Scholar 

  36. Ookhtens M, Mittur AV, Erhart NA: Changes in plasma glutathione concentrations, turnover, and disposal in developing rats. Am J Physiol. 1994, 266: R979-R988.

    CAS  PubMed  Google Scholar 

  37. Meister A: Glutathione. The Liver: Biology and Pathobiology. Edited by: Aria IM, Jakobi WB, Poppser , Shchacter D, Shafritz D. 1983, New York: Raven Press, 401-417.

    Google Scholar 

  38. Fukagawa N, Ajami A, Young V: Plasma methionine and cysteine kinetics in response to an intravenous glutathione infusion in adult humans. Am J Physiol. 1996, 270 (2 Pt 1): E209-E214.

    CAS  PubMed  Google Scholar 

  39. Inoue M, Kinne R, Tran 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 

  40. Tallen HH, Moore S, Stein WH: L-cystathionine in human brain. J Biol Chem. 1958, 230: 707-716.

    Google Scholar 

  41. Griffith OW: Biologic and pharmacologic regulation of mammalian glutathione synthesis. Free Radic Biol Med. 1999, 27 (9-10): 922-935. 10.1016/S0891-5849(99)00176-8.

    Article  CAS  PubMed  Google Scholar 

  42. Samiec PS, Drews-Botsch C, Flagg EW, Kurtz JC, Sternberg P, Reed RL, Jones DP: Glutathione in human plasma: decline in association with aging, age-related macular degeneration, and diabetes. Free Radic Biol Med. 1998, 24 (5): 699-704. 10.1016/S0891-5849(97)00286-4.

    Article  CAS  PubMed  Google Scholar 

  43. Mackenzie B, Ahmed A, Rennie MJ: Muscle amino acid metabolism and transport. Mammalian Amino Acid Transport: Mechanisms and Control. Edited by: Kilberg MS, Haeusinger D. 1992, Plenum Press: New York, 195-231.

    Chapter  Google Scholar 

  44. Kilberg MS, Handlogten ME, Christensen HN: Characteristics of system ASC for transport of neutral amino acids in the isolated rat hepatocyte. J Biol Chem. 1980, 256: 3304-3312.

    Google Scholar 

  45. Kilberg MS, Haeussinger D, eds: Mammalian Amino Acid Transport. 1992, New York, Plenum

    Google Scholar 

  46. Kondo T, Dale GL, Beutler E: Glutathione transport by inside-out vesicles from human erythrocytes. PNAS. 1980, 77: 6359-6362. 10.1073/pnas.77.11.6359.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  47. Lauterburg BH, Adams JD, Mitchell JR: Hepatic glutahione homeostasis in the rat: efflux accounts for glutathione turnover. Hepatology. 1984, 4: 586-590.

    Article  CAS  PubMed  Google Scholar 

  48. Tateishi N, Higashi T, Shinya S, Naruse A, Sakamoto Y: Studies on the regulation of glutathionine level in rat liver. J Biochem. 1974, 75: 93-103.

    CAS  PubMed  Google Scholar 

  49. Tateishi N, Sakamoto Y: Nutritional significance of glutathione in rat liver. Glutathione: Storage, Transport, and Turnover in Mammals. Edited by: Sakamoto Y, Higashi T, Tateishi N. 1983, Tokyo: Scientific Socoetu Press, 13-38.

    Google Scholar 

  50. Meister A: Glutathione. The Liver: Biology and Pathobiology. Edited by: Aria IM, Jakobi WB, Poppser, Shchacter D, Shafritz D. 1998, New York: Raven Press, 401-417. 2

    Google Scholar 

  51. Fernandez-Checa J, Lu SC, Ookhtens M, DeLeve L, Runnegar M, Yoshida H, Saiki H, Kannan R, Garcia-Ruiz C, Kuhlenkamp JF, Kaplowitz N: The regulation of hepatic glutathione. Hepatic Anion Transport and Bile Secretion: Physiology and Pathophysiology. Edited by: Tavoloni N, Berk PD. 1992, New York: Marcel Dekker, 363-395.

    Google Scholar 

  52. Storch KJ, Wagner DA, Burke JF, Young VR: Quantitative study in vivo of methionine cycle in humans using [methyl-2H3]- and [1-13C]methionine. Am J Physiol. 1988, 255: E322-31.

    CAS  PubMed  Google Scholar 

  53. Finkelstein JD, Harris BJ, Martin JJ: Methionine metabolism in mammals: Concentrations of metabolites in rat tissues. Kinetic studies of betaine:homocysteine methyltransferase. J Nutr. 2000, 112: 1011-1018.

    Google Scholar 

  54. Finkelstein JD, Martin JJ: Methionine metabolism in mammals: Adaptation to methionine excess. J Biol Chem. 1986, 261: 1582-1587.

    CAS  PubMed  Google Scholar 

  55. Corrales FJ, Perez-Mato I, Sanchez del Pino MM, Ruiz F, Castro C, Garcia-Trevijano ER, Latasa U, Martinez-Chantar ML, Martinez-Cruz A, Avila MA, Mato JM: Regulation of mammalian liver methionine adenosyltransferase. J Nutr. 2002, 132: 2377S-2381S.

    CAS  PubMed  Google Scholar 

  56. Nijhout HF, Reed MC, Ulrich CM: Multiple and diverse mechanisms of homeostasis in a complex metabolic network. 2007,

    Google Scholar 

  57. Corrales FJ, Ruiz F, Mato J: In vivo regulation by glutathione of methionine adenosyltransferase S-nitrosylation in rat liver. J Hepatology. 1999, 31: 887-894. 10.1016/S0168-8278(99)80291-8.

    Article  CAS  Google Scholar 

  58. Pajares MA, Duran C, Corrales F, Pliego MM, Mato JM: Modulation of Rat Liver S-Adenosylmethionine Synthetase Activity by Glutathione. J Biol Chem. 1992, 267: 17598-17605.

    CAS  PubMed  Google Scholar 

  59. Epstein CJ: Down syndrome (trisomy 221). The Metabolic and Molecular Bases of Inherited Disease. Edited by: Stanbury JB, Wingarden JB, Frederickson DS. 1995, New York: McGraw-Hill, 749-795. 3

    Google Scholar 

  60. Krivchenia E, Hether CA, Edmonds LD, May DS, Guckenberger S: Comparative epidemiology of Down syndrome in two United States Populations. Amer J Epidem. 1993, 137: 815-828.

    CAS  Google Scholar 

  61. Torfs CP, Christianson RE: Anomalies in Down syndrome individuals in a large population-based registry. Amer J Med Gen. 1998, 77: 431-438. 10.1002/(SICI)1096-8628(19980605)77:5<431::AID-AJMG15>3.0.CO;2-J.

    Article  CAS  Google Scholar 

  62. Ianello RC, Crack PJ, de Hann JB, Kola I: Oxidative stress and neural dysfunction in Down syndrome. J Neural Transm. 1999, 57: 257-267.

    Google Scholar 

  63. Lord C, Risi S, Lambrecht L, Cook EH, Leventhal BL, DiLavore PC, Pickles A, Rutter M: The autism diagnostic observation schedule-generic: A standard measure of social and communication deficits associated with the spectrum of autism. J Autism Dev Disord. 2000, 30: 205-223. 10.1023/A:1005592401947.

    Article  CAS  PubMed  Google Scholar 

  64. CDC. 2005,

  65. London EA: The environment as an etiologic factor in autism: a new direction for research. Eviron Health Perspec. 2000, 108 (suppl): 401-404. 10.2307/3454527.

    Article  Google Scholar 

  66. Folstein SE, Rosen-Sheidley B: Genetics of autism: complex aetiology for a heterogeneous disorder. Nature Rev Genet. 2001, 2: 943-955. 10.1038/35103559.

    Article  CAS  PubMed  Google Scholar 

  67. Muhley R, Trentacoste SV, Rapin I: The genetics of autism. Pediatr. 2004, 113: e472-e486. 10.1542/peds.113.5.e472.

    Article  Google Scholar 

  68. Zimmerman AW, Jyonouchi H, Comi AM, Connors SL, Milstien S, Varsou A, Heyes MP: Cerebrospinal fluid and serum markers of inflammation in autism. Pediatr Neurol. 2005, 33: 195-201. 10.1016/j.pediatrneurol.2005.03.014.

    Article  PubMed  Google Scholar 

  69. Horvath K, Perman JK: Autistic disorder and gastrointestinal disease. Curr Opin Pediatr. 2002, 14: 583-587. 10.1097/00008480-200210000-00004.

    Article  PubMed  Google Scholar 

  70. Stubbs G, Litt M, Lis E: Adenosine deaminase activity decreased in autism. J Am Acad Child Psych. 1982, 21: 71-74.

    Article  CAS  Google Scholar 

  71. Decking UK, Schlieper G, Kroll K, Schrader J: Hypoxia-induced inhibition of adenosine kinase potentiates cardiac adenosine release. Circ Res. 1997, 81: 154-164.

    Article  CAS  PubMed  Google Scholar 

  72. Chen YF, Li PL, Zou AP: Oxidative stress enhances the production and actions of adenosine in the kidney. Am J Physiol Regul Integr Comp. 2001, 281: R1808-16.

    CAS  Google Scholar 

  73. Melnyk S, Jernigan S, Savenka S, James SJ: Elevation in S-adenosylhomocysteine and DNA hypomethylation in parents and children with autism. FASEB J. 2007, 21: A348-

    Google Scholar 

  74. Ming X, Cheh MA, Yochum CL, Halladay AK, Wagner GC: Evidence of oxidative stress in autism derived from animal models. Amer J Biochem Biotech. 2008, 4: 218-225.

    Article  CAS  Google Scholar 

  75. Chien PFW, Smith K, Peter WW, Scrimgeour CM, Talor DJ, Rennie MJ: Protein turnover in the human fetus studied a term using stable isotope tracer amino acids. Am J Physiol. 1993, 265: E31-E35.

    CAS  PubMed  Google Scholar 

  76. FAO/WHO/UNU: Energy and protein requirements. 1985, Geneva: WHO

    Google Scholar 

  77. Waterlow JC: Whole-body protein turnover in humans – past, present, and future. Ann Rev Nutr. 1995, 15: 57-92.

    Article  CAS  Google Scholar 

Download references


The authors benefited greatly from the advice of Jesse Gregory, III, University of Florida, Barry Shane, University of California, Berkeley, and Marian Neuhouser, Fred Hutchinson Cancer Research Center. This research was supported by NIH grant RO1 CA105437 (C.M.U.) and NSF grant DMS-0616710 (M.C.R.).

Author information

Authors and Affiliations


Corresponding author

Correspondence to Michael C Reed.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HN, MR, JJ, and CM formulated the questions and contributed substantially to the design of the project. RT and JP wrote the code and conducted the in silico experiments with HN and MR. The initial manuscript was written by RT, JP, MR, and HN. All authors read the manuscript, contributed revisions, and approved the submitted manuscript.

Electronic supplementary material


Additional file 1: Supplementary Material – Model Details for A Mathematical Model of Glutathione Metabolism. A full description of the mathematical model is given. (PDF 186 KB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Reed, M.C., Thomas, R.L., Pavisic, J. et al. A mathematical model of glutathione metabolism. Theor Biol Med Model 5, 8 (2008).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Down Syndrome
  • Autistic Child
  • Glutathione Metabolism
  • Constituent Amino Acid
  • Transsulfuration Pathway