In an integrative approach to gaining a better understanding of cell-wide molecular networks at the systems level (systems biology), many attempts at computer simulations of metabolic networks have been made. Because of the simplicity of its structure and components, and the availability of kinetic information, erythrocyte metabolism has been the focus of mathematical modeling for over three decades [1], and is a leading example of the use of mathematical modeling not only for understanding biochemical regulation, but also as the basis for the construction of other mathematical frameworks [2][3][4][5]. For example, models of erythrocyte metabolism were able to predict the importance of the de novo synthesis of glutathione and its mechanism of action in glucose-6-phosphate dehydrogenase (G6PDH)-deficient cells [6], and to analyze Methemoglobin (metHb), an oxidized form of hemoglobin, is unable to bind and carry oxygen. Erythrocytes are continuously subjected to oxidative stress and nitrite exposure, which results in the spontaneous formation of metHb. To avoid the accumulation of metHb, reductive pathways mediated by cytochrome b5 or flavin, coupled with NADH-dependent or NADPH-dependent metHb reductases, respectively, keep the level of metHb in erythrocytes at less than 1% of the total hemoglobin under normal conditions. In this work, a mathematical model has been developed to quantitatively assess the relative contributions of the two major metHbreducing pathways, taking into consideration the supply of NADH and NADPH from central energy metabolism. The results of the simulation experiments suggest that these pathways have different roles in the reduction of metHb; one has a high response rate to hemoglobin oxidation with a limited reducing flux, and the other has a low response rate with a high capacity flux. On the basis of the results of our model, under normal oxidative conditions, the NADPH-dependent system, the physiological role of which to date has been unclear, is predicted to be responsible for most of the reduction of metHb. In contrast, the cytochrome b5-NADH pathway becomes dominant under conditions of excess metHb accumulation, only after the capacity of the flavin-NADPH pathway has reached its limit. We discuss the potential implications of a system designed with two metHbreducing pathways in human erythrocytes.Abbreviations cytb5, cytochrome b5; cytb5R, (NADH-dependent) cytochrome b5 reductase; FR, (NADPH-dependent) flavin reductase; GAPDH, glyceraldehyde-3-phosphate dehydrogenase; G6PDH, glucose-6-phosphate dehydrogenase; GSSG, oxidized glutathione; GSH, glutathione (reduced form); LDH, lactate dehydrogenase; metHb, methemoglobin.