The erythrocyte is exposed to reactive oxygen species in the circulation and also to those produced by autoxidation of hemoglobin. Consequently, erythrocytes depend on protection by the antioxidant glutathione. Mathematical models based on realistic kinetic data have provided valuable insights into the regulation of biochemical pathways within the erythrocyte but none have satisfactorily accounted for glutathione metabolism. In the current model, rate equations were derived for the enzyme-catalyzed reactions, and for each equation the nonlinear algebraic relationship between the steady-state kinetic parameters and the unitary rate constants was derived. The model also includes the transport processes that supply the amino acid constituents of glutathione and the export of oxidized glutathione. Values of the kinetic parameters for the individual reactions were measured predominately using isolated enzymes under conditions that differed from the intracellular environment. By comparing the experimental and simulated results, the values of the enzyme-kinetic parameters of the model were refined to yield conformity between model simulations and experimental data. Model output accurately represented the steady-state concentrations of metabolites in erythrocytes suspended in plasma and the changing glutathione concentrations in whole and hemolyzed erythrocytes under specific experimental conditions. Analysis indicated that feedback inhibition of ␥-glutamate-cysteine ligase by glutathione had a limited effect on steady-state glutathione concentrations and was not sufficiently potent to return glutathione concentrations to normal levels in erythrocytes exposed to sustained increases in oxidative load.Erythrocyte glutathione plays a vital role in mitigating the damaging effects of reactive oxygen species (ROS) 4 encountered in the circulation (1) and produced by continuous oxidation of hemoglobin within the cytosol of the erythrocyte (2, 3). Reduced glutathione (GSH) reacts with superoxide reaction products, and via specific enzymes, it degrades hydrogen peroxide and lipid peroxides (glutathione peroxidase, EC 1.11.1.9), and covalently modifies toxic xenobiotics and endogenous electrophiles (glutathione S-transferase, EC 2.5.1.18) to form water-soluble conjugates that are exported from the erythrocyte for excretion (4). In diseases associated with increased production of ROS that results in GSH depletion, restoration of the normal erythrocyte GSH concentration has been shown to have positive therapeutic effects (5, 6).Experimental studies have provided a wealth of information on the individual reactions that together underlie the metabolic processes of the human erythrocyte (7)(8)(9)(10)(11)(12)(13)(14). This information has been used to develop mathematical models to analyze and then predict the way in which the kinetic characteristics and control mechanisms for each reaction combine to regulate particular metabolic processes. Comprehensive models of erythrocyte glycolysis, the pentose phosphate pathway, 2,3-bisphosphogl...