|Home | About | Journals | Submit | Contact Us | Français|
Conceived and designed the experiments: PMBMC AS MAS. Performed the experiments: PMBMC. Analyzed the data: PMBMC AS MAS. Wrote the paper: PMBMC AS MAS.
The NADPH redox cycle plays a key role in antioxidant protection of human erythrocytes. It consists of two enzymes: glucose-6-phosphate dehydrogenase (G6PD) and glutathione reductase. Over 160 G6PD variants have been characterized and associated with several distinct clinical manifestations. However, the mechanistic link between the genotype and the phenotype remains poorly understood.
We address this issue through a novel framework (design space) that integrates information at the genetic, biochemical and clinical levels. Our analysis predicts three qualitatively-distinct phenotypic regions that can be ranked according to fitness. When G6PD variants are analyzed in design space, a correlation is revealed between the phenotypic region and the clinical manifestation: the best region with normal physiology, the second best region with a pathology, and the worst region with a potential lethality. We also show that Plasmodium falciparum, by induction of its own G6PD gene in G6PD-deficient erythrocytes, moves the operation of the cycle to a region of the design space that yields robust performance.
In conclusion, the design space for the NADPH redox cycle, which includes relationships among genotype, phenotype and environment, illuminates the function, design and fitness of the cycle, and its phenotypic regions correlate with the organism's clinical status.
The NADPH redox cycle plays a key role in the oxidative stress response of human erythrocytes. It consists of two enzymes: glucose-6-phosphate dehydrogenase (G6PD, EC 220.127.116.11) and glutathione reductase (GSR, EC 18.104.22.168). Although variants of G6PD have been intensively studied and are associated with several distinct clinical manifestations, the relationship between the genotype and the phenotype is still poorly understood. To address this issue, we have constructed a “system design space” which facilitates the quantitative comparison of wild-type and variants for the redox cycle. Our results identify three different phenotypes that correlate with clinical manifestations.
G6PD catalyses the first step of the hexose-monophosphate shunt (Figure 1A), which provides pentoses for nucleic acid synthesis and regenerates NADPH. In erythrocytes, NADPH is required for various processes, but most of it is oxidized by GSR . The latter process regenerates reduced glutathione (GSH) that is oxidized in the repair of oxidative damage. In mice, and presumably in other organisms, G6PD is dispensable for pentose synthesis but essential for defense against oxidative stress . High levels of G6PD exist for this function, but under pronounced oxidative stress hexokinase (EC 22.214.171.124) becomes rate-limiting for the NADPH supply .
Previous quantitative analysis of the NADPH redox cycle ,  indicates that normal G6PD activity is sufficient but not superfluous to avoid NADPH depletion and ensure timely adaptation of the NADPH supply during pulses of oxidative load such as those that occur during adherence of erythrocytes to phagocytes. The quantitative analysis of this system has been facilitated by two recent developments: a method for constructing the “system design space”  and a related method for calculating “global tolerances”  to large variations in the values of system parameters and environmental inputs.
In this paper we utilize the design space as a framework to compare the quantitative phenotypes of wild-type and mutant variants of the NADPH redox cycle. In particular, 160 G6PD variants have been characterized , and there are several distinct phenotypes associated with G6PD deficiency . This system presents a unique opportunity to relate genotype to phenotype by focusing on the quantitative behavior of the integrated NADPH redox cycle. Our analysis of this system and its mutants requires additional background regarding its biochemistry, genetics and clinical manifestations.
G6PD in its active form is made up of two or four identical subunits, each with a molecular mass of 59 kDa. The gene for G6PD is on the X-chromosome, and the deficiency is inherited in a sex-linked fashion. Hemizygous males and homozygous females with low-activity mutant forms of the enzyme carry only G6PD-deficient erythrocytes. However, female heterozygotes carry both normal and deficient erythrocytes. The latter outcome is due to the fact that each cell inactivates one of its X-chromosomes, chosen at random. The most common G6PD form worldwide is the B form. However, in Africa up to 40% of the population carry the non-deficient A form of G6PD . The most common G6PD-deficient variant in Africa is the A- allele with a frequency between 0 and 25% .
Over 400 million people in the world suffer from G6PD deficiency; which makes it the most common known enzymopathy. The highest prevalence rates are found in tropical and subtropical regions of the world and in some areas of the Mediterranean. It would be difficult to understand the existence of such a widespread enzymatic deficiency without a counterbalancing biological advantage. The observation that the G6PD A- mutation has fitness costs such as making individuals more prone to develop hemolytic anemias  makes this issue even more intriguing. In 1960, Allison  and Motulsky  first suggested that individuals deficient in G6PD might be at a selective advantage in malaria endemic areas. However, numerous clinical studies have yielded conflicting results over the alleged protective role of G6PD deficiency against malaria. In particular, clinical studies of severe malaria have shown a protective role for G6PD deficiency in hemizygous males ,  while, in heterozygous females, there is either a protective effect  or no effect . The studies relating to uncomplicated malaria have shown more conflicting results with heterozygous females protected , , , , at increased risk , or unaffected by G6PD deficiency . In a recent paper, Johnson M.K. et al.  compared the association between uncomplicated malaria incidence and G6PD deficiency using the two different methods currently available: G6PD enzyme activity and G6PD genotype assessment. They found a 52% reduced risk of uncomplicated malaria in G6PD-deficient females when the deficiency was assessed using enzyme activity. When the assessment method was based on the genotype, the protective association was no longer significant, which is likely due to the large dispersion in enzymatic activity caused by random X-chromosome inactivation . For G6PD deficient males, there was no association with the incidence of uncomplicated malaria regardless of assessment method. Their work suggests that the conflict between prior association studies might be related to the method of G6PD assessment.
G6PD variants have been described on the basis of their biochemical properties and grouped into five classes according to the level of enzyme activity and clinical manifestations. In 1967, a committee of the World Health Organization  recommended standard techniques for the biochemical characterization of G6PD variants. Class I includes severely deficient variants that suffer from chronic non-spherocytic haemolytic anemia. Class II variants have less than 10% enzyme activity but do not suffer from chronic non-spherocytic haemolytic anemia. Class III variants are moderately deficient (10–60% activity). Class IV has normal enzyme activity (60–150%) and in Class V the enzyme activity is higher than normal (>150% activity) . Of the 400 million people worldwide that exhibit G6PD deficiency, there are, for instance, 4 million African Americans (1 in 10)  in a country where malaria has been eradicated since 1951 . Thus, the implications of G6PD deficiencies unrelated to malaria represent a significant medical problem .
In terms of clinical manifestations, G6PD deficiency can be divided into 3 groups. The first group (Class IV and V) consists of variants that do not have any clinical disorder. The second group (Class II and III) shows a hemolytic reaction that is triggered by certain drugs, by infection or by the ingestion of certain foods. The mechanism of hemolysis is not known, but it has been proposed that it results from the inability of G6PD-deficient erythrocytes to cope with the oxidative damage produced by the agents mentioned above . The third group (Class I) may be classified as having chronic non-spherocytic hemolytic anemia, which is exacerbated by oxidative stress.
In this paper, we assembled existing data into an integrated model of the NADPH redox cycle that allowed us to address the following questions: given the variability in G6PD, how do mutations affect the quantitative performance of the cycle? Is there a correlation between clinical manifestations and poor performance of the NADPH redox cycle? Finally, how does infection of erythrocytes by Plasmodium falciparum affect the redox cycle. We base our analysis of these questions on the model shown in Figure 1B.
Under physiological conditions, the steps catalyzed by G6PD and GSR are essentially irreversible . The kinetic parameters of each enzyme and the nominal concentrations of their reactants have been well characterized (Table 1).
The rate expression for G6PD follows a compulsory order mechanism in which binds first to the enzyme:
is the apparent Michaelis-Menten constant for NADP+ after taking into account the competitive inhibition by NADPH and 2,3-diphosphoglycerate (2,3-DPG). In Table 1, we present the experimental values of kinetic parameters for G6PD encoded by the B-form of the gene.
At concentrations of oxidized glutathione (GSSG) typical of low to moderate oxidative loads, GSR follows a ping-pong mechanism with rate expression:
The kinetic parameters of GSR are given in Table 1.
Using the method outlined in , we combined the fundamental enzyme kinetic rate laws and the conservation of the NADP moiety into a system of ordinary differential equations representing mass balance, which can then be solved to determine the steady-state of integrated system. With this approach, we are able to formulate the piecewise power-law representation of G6PD [Eq(3)] and GSR [Eq (4)] kinetics in normalized form:
Note that represents the apparent Michaelis constant for NADP+ scaled by NADPtotal and it takes into account product-inhibition and the inhibition by 2,3-DPG. Conversely, represents the apparent Michaelis constant for NADPH scaled by NADPtotal. The parameter ρ represents the ratio between the apparent maximal velocity of GSR and the apparent maximal velocity of G6PD.
Based on this representation, we recognize that the NADPH redox cycle in human erythrocytes can operate under three different meaningful steady-state regimes (Table 2) and one unrealistic case (with ) (see Text S1). Each systemic regime will only apply to a particular region of the design space (see Figure 2), which is obtained by using the methods previously described , .
The system representation within each regime is a simple but nonlinear S-system equation for which determination of local behavior, after appropriate transformation, reduces to conventional linear analysis . Thus, the local behavior is completely determined and readily characterized by the evaluation of the following quantitative indices.
The change of concentration (e.g., ) or flux in response to a change in value for an independent variable (e.g., ) is defined by a relative derivative of the explicit steady-state solution. For example,
The change of concentration or flux in response to a change in value for one of the parameters that define the structure of the system (e.g., Michaelis constants) is also determined by a relative derivative of the explicit steady-state solution. For example,
The eigenvalue, which in the present case is an inverse measure of the system's response time, is determined by a first-order Taylor series approximation in logarithmic space of the differential equation about the steady-state that applies for each systemic regime.
The performance of the NADPH redox cycle can be evaluated in each systemic regime according to the following quantitative criteria:
The concentration of should be well buffered against:
The supply of () should
The sensitivity of the supply of to changes in the concentration of should
The response time should
All of the statistical tests described below, except one, were performed with algorithms in Mathematica available through reference .
Single-factor Analysis of Variance (ANOVA) is a method used to decide whether differences exist among several population means. It requires all populations to be normally distributed (D'Agostino-Pearson Test for nonnormality) and exhibit the same variance (Bartlett's test for equality of variance). We tried to transform the data so as to satisfy these assumptions whenever either was violated. Whenever the ANOVA rejected equality of all means we investigated each pairwise comparison of population means for significant differences using Tukey's honestly significant difference (HSD) test.
Whenever the transformation procedure failed to satisfy the ANOVA assumptions, we used a nonparametric method: the Kruskal-Wallis One Way Analysis of Variance. If it was significant, it indicated that at least one of the populations was different from at least one of the others. We used a procedure called “kruskalmc” written for R (included in package pgirmess ) to perform a multiple comparison test between populations after the Kruskal-Wallis test (the nonparametric equivalent of a Tukey's honestly significant difference test). This test helped to determine which populations were different in pairwise comparisons adjusted appropriately .
The local performance within the three systemic regimes is determined by the above methods and evaluated according to the previously defined criteria. Our aim is to ascertain which of the systemic regimes is better suited for effective performance of the NADPH redox cycle.
In Table 3, we summarize the results from the analysis of local performance in Systemic Regimes a, b and c (for further details see Text S1). These results show that local performance in Systemic Regime a fulfills all the criteria defined above. In contrast, performance in Systemic Regimes b and c cannot fulfill criteria 4 and 5 because there is no response to changes in the concentration of GSSG. Moreover, the local robustness (criterion 1) of systems represented in Systemic Region b is significantly worse than that of systems in either Systemic Region a or c (see Text S1).
Additionally, although Systemic Regimes b and c can exhibit a fast response time (Criterion 6), it will not be with respect to changes in GSSG. Therefore, the importance of this responsiveness is questionable.
In summary, we predict that, under basal conditions, the NADPH redox cycle should operate in Systemic Regime a, which has the best overall local performance. Moreover, natural selection should maintain the operating point far from the boundary to region c, which has less desirable overall local performance, and especially from the boundary to region b, which has the worst performance.
The boundaries of Systemic Regime a are obtained by inserting the linear solution into the corresponding linear breakpoint conditions:
The result is the following set of boundaries
Equations (8) can be represented in terms of the fundamental parameters and variables by substituting the definitions given previously
Systems represented within these boundaries exhibit the best local performance and thus these boundaries provide the basis for the analysis of global tolerance.
The most common variant of G6PD worldwide is the B form. Given the detailed description of this variant form of G6PD (Table 1) and of GSR, we are able to locate the operating point for G6PD B (Figure 3, black circle) in the design space and show the normalized steady-state concentration of NADPH in the z-direction with a heat map.
The local behavior of the NADPH redox cycle can be evaluated according to the seven criteria (Table 4, upper panel) described earlier and, given the boundaries surrounding Systemic Region a (Eq. 9), we are able to determine the numerical value of global tolerance for each of the kinetic parameters and independent concentration variables (Table 4, lower panel). We will use the expression “” to describe the global tolerances, where tolerance to a fold decrease and tolerance to a fold increase (since boundaries can be crossed either by decreasing or increasing a parameter).
As can be seen from the results in Table 4 (upper panel), natural selection results in a design that has nearly optimal local performance when we consider the normal variant of G6PD (B-form). In a similar trend, the global tolerances (Table 4, lower panel) range from the smallest of 63 fold to the largest of 1.4×105. The smallest value is associated with the concentrations of and , whereas the largest is associated with the product inhibition constant of G6PD.
The biochemical properties of more than 400 putative variants have been tabulated in a review . Detailed comparisons of the tabulated values and the values in the corresponding original literature show a number of irregularities. Thus, we have systematically re-evaluated the original literature and produced a revised set of tables for the various classes that highlights the discrepancies (see Text S2).
After eliminating the anomalous cases, we focused on the remaining G6PD variants for which there are documented experimental values for: , , and . Furthermore, all variants under analysis were characterized according to the standard techniques that the World Health Organization  recommended in 1967. We excluded variants (see  and Text S2) for which one or more of the kinetic parameters is unknown, or whose kinetics were characterized according to non-standard techniques.
In order to analyze the local performance and global tolerance of the NADPH redox cycle in which variants of G6PD participate, we assumed that , , the kinetic parameters of GSR and all the independent variables retain their normal value.
In Figure 3, we display the location of the operating point of the G6PD variants in the Design Space for the NADPH redox cycle. Class III (blue circle) and Class IV (orange circle) variants have their operating points located well within Systemic Region a. Moreover, we can visually ascertain that Class IV variants are further away from Systemic Region c than Class III variants. However, Class II (green circle) and Class I (red circle) variants have operating points that span Systemic Regions a and c. This result proves important because we know that the Class I and II variants have clinical manifestations. Therefore, there is a correlation between a Class being near to or present in Systemic Region c and the existence of clinical problems. Furthermore, several considerations suggest that if the operating point of a variant of G6PD were to be located within Systemic Region b, this might be lethal to the host. First, Systemic Region b has the worst local performance. Since being present in Systemic Region c is already associated with life-threatening medical conditions, we can only expected that being present in an even worse phenotypic region could be lethal. Second, no known variant of G6PD leads to an operating point located within Systemic Region b. Moreover, the observation that all known variants of G6PD are located in the design space () such that an increase in oxidative stress would never move their operating points into Systemic Region b suggests a selective pressure to prevent even temporary excursions into this region.
The local behavior of the NADPH redox cycle can be evaluated according to the seven criteria described in the METHODS section and our aim is to distinguish among the several classes of G6PD variants (Table 5). Note that we are only considering the G6PD variants that lead to systems whose operating points lie within Systemic Region a. Moreover, in this Regime the fulfillment of Criteria 4 and 5 only depends on the concentration of GSSG and on the Michaelis constant of GSR with respect to GSSG (see Text S1), which are fixed in our analysis. Thus, our comparisons will be limited to criteria 1, 2, 3, 6 and 7.
From the results in Table 5, we conclude that the worst local performance is given by systems for which G6PD variants are in Class I or II. In detail, according to criteria 1, 2, 3 and 6, the local performance of systems in Class I or II is statistically worse than those in either Class III or IV. Furthermore, there is no criterion that allows us to distinguish between Class I and Class II, whereas the local performance of a cycle with a Class IV variant is statistically better than a cycle with a Class III variant according to criteria 2, 3 and 6.
Once again, we are only considering the G6PD variants that lead to systems operating within Systemic Region a. Therefore, we will determine how much a given kinetic parameter or independent concentration variable has to change in order to place the NADPH redox cycle in a different Systemic Region (Table 6).
From the analysis of Table 6, we find that the worst global tolerances are associated with systems of Class I or II. According to all the kinetic parameters and independent variables except , the global tolerances of Class I or II cycles are statistically worse than those of either Class III or IV. In addition, we see that there is no global tolerance that allows us to distinguish between Class I and Class II. Furthermore, the analysis of global tolerance, except for and , allows us to differentiate between Class III and Class IV and conclude that Class IV has better global tolerances than Class III.
Our analysis has shown that there is no local performance criterion or global tolerance that distinguishes a NADPH redox cycle with a Class I variant from one with a Class II variant of G6PD. However, the clinical manifestations are quite different. How might these differences be reconciled? One possibility is differences in kinetic properties of Class I and Class II variants of G6PD. In Table 7, we display the medians of the kinetic properties of the G6PD variants, when grouped into their respective classes. These results show that it is possible to make this distinction on the basis of : the median of Class I is 7.6 µM while the median of Class II is 30.5 µM. However, this statistical difference in kinetic properties does not lead to a statistically significant difference in local performance or global tolerance (Tables 5 and and6).6). Therefore, the difference in kinetic properties is unable to explain the different clinical manifestations.
Within the context of our model, there are several ways that these differences can be reconciled. A plausible example would be Class I variants having a higher value of . For instance, if the were 790 µM in Class I variants rather than the normal value of 7.9 µM, then the design space of the NADPH redox cycle would resemble Figure 4, and all Class I variants would operate in the poor region, Systemic Region c, whereas most Class II variants would operate in the good region, Systemic Region a.
It has been shown that P. falciparum can adapt to growth in G6PD-deficient erythrocytes by producing its own G6PD enzyme .
In one particular study , researchers purified G6PD from infected and uninfected erythrocytes. They found that infected erythrocytes had two different variants of G6PD: the G6PD that P. falciparum produces, responsible for 90% of the overall activity, and the host's G6PD, responsible for 10%.
In Table 8, we list the experimental values for the kinetic properties of the host G6PD and the P. falciparum G6PD. By assuming normal values for the unknown independent variables and parameters, we are able to plot (Figure 5) the location of the operating point for the NADPH redox cycle in erythrocytes that are infected (black rectangle) or uninfected (black diamond) with P. falciparum.
As can be seen in Figure 5, the uninfected erythrocytes endure a NADPH redox cycle with poor performance due to the G6PD deficiency (the operating point is located within Systemic Region c). However, once P. falciparum infects an erythrocyte and synthesizes its own G6PD, the operating point of the NADPH redox cycle will move into the region of best performance (Systemic Region a). Furthermore, we are able to assess the local performance and global tolerance of the NADPH redox cycle in human erythrocytes infected with P. falciparum (Table 9). The local performance is comparable to the performance of erythrocytes in which G6PD is of the B-form (Table 4). The major difference is in temporal responsiveness. Infected erythrocytes are more sluggish in responding to changes in GSSG. In terms of global tolerances, we observe that all tolerances for the infected erythrocyte are lower. This observation is correlated with the fact that the normal operating point of the NADPH redox cycle (G6PD B-form – black circle in Figure 5) is located further away from the boundary to Systemic Region c than the operating point of the infected erythrocytes (black rectangle). It is known that P. falciparum-infected erythrocytes are under a larger oxidative stress  and this effect may move its operating point in the design space further to the right (black rectangle).
We have to consider these results from the viewpoint of the P. falciparum. Upon infection of a G6PD-deficient erythrocyte, P. falciparum induces its own G6PD . According to our analysis, this induction tends to correct the state of the redox cycle and P. falciparum can rely on a functional NADPH moiety supply unit. This scenario raises the question: is there a reason for P. falciparum to need a reliable NADPH moiety supply unit? P. falciparum digests the erythrocyte's hemoglobin and, as a consequence, produces , and OH radicals leading to a larger oxidative stress which the otherwise deficient cycle could not cope .
The piecewise power-law representation of the rate laws for G6PD and GSR allowed us to define the design space for the NADPH redox cycle. This design space has three phenotypic regions, each defined by a single type of stable steady state. However, the redox cycle only works properly in Systemic Region a.
In comparing the local performance and global tolerance of the several possible cycles, it becomes apparent that an individual who has inherited the normal form of G6PD (B-form) displays one of the best possible NADPH redox cycles. One may argue that these functional advantages may explain the high prevalence of this form worldwide. In broader terms, we have demonstrated that there are statistically significant differences in local performance and global tolerance between cycles with a Class IV variant vs. a Class III/II/I variant. Individuals with a Class III variant of G6PD will experience a more sluggish cycle and an increased sensitivity to changes in glucose-6-phosphate and oxidized glutathione while, for example, not being able to tolerate an environment in which a large oxidative load would lead to a decrease in the concentration of glucose-6-phosphate to less than 0.3 µM. Individuals with a Class II G6PD variant will have a NADPH redox cycle that has an even slower response time as well as an even greater sensitivity to changes in the concentrations of NADPtotal, oxidized glutathione and glucose-6-phosphate. Furthermore, these patients will not tolerate an environment in which a large oxidative load would lead to a decrease in the concentration of glucose-6-phosphate to less than 2 µM.
Individuals who have inherited a Class I variant of G6PD suffer from chronic hemolytic anemia, whereas individuals who have inherited a Class II variant develop significant hemolysis only when oxidative stress is induced. An increase in oxidative stress leads to a shift in the operating point toward higher values of ρ and, concomitantly, to a possible transition into Systemic Region c. These observations have prompted us to set forth the following hypothesis:
When the NADPH redox cycle functions within Systemic Region c, the individual will suffer from hemolytic anemia. Furthermore, individuals with a Class II/III variant only transiently operate within Systemic Region c because they require an extra source of oxidative stress to develop clinical manifestations. However, individuals with a Class I variant always operate within Systemic Region c because their hemolytic anemia is chronic.
To reconcile this hypothesis with our observations, we need to understand why Class I variants aren't mostly located within Systemic Region c (Figure 3). One possibility may reside in the stability of Class I variants of G6PD. Cheng et al  suggested that a correlation exists between mutations in evolutionary conserved amino acids and clinical manifestations. They showed the following trend: “as the severity of clinical consequences reduces from Class I, II and III, the percentage of mutations in the top pentad (conservation score >0.8) decreases from 51% (21/41) for class I to 41% (12/29) for class II and to 24% (5/21) for Class III.” Furthermore, the majority of Class I mutations are clustered near the dimer interface and NADP binding site for structural stability, or are deletion mutations . In addition, of 40 Class I variants described in , 32 had highly decreased thermostability while 4 had decreased thermostability and 4 were classified as normal. Hence, one may speculate that Class I variants have a lower stability. In turn, a lower stability may lead to a sharp decrease in activity as erythrocytes age, eventually attaining a threshold (within Systemic Region c) below which erythrocytes lyse. Chronic non-spherocytic haemolytic anemia ensues from the hematopoietic system being unable to compensate for the shortened erythrocyte lifespan with sufficiently increased reticulocytosis. On the other hand, because G6PD in these patients is sampled from a younger cell population, the overall G6PD activities are similar to those of Class II variants, which have longer erythrocyte lifespans.
A systematic reexamination of our modeling assumptions by means of the design space analysis suggests other testable possibilities for reconciling the observations with our hypothesis above. The most plausible among these is the prediction that Class I mutants exhibit a higher dissociation constant of the G6PD-NADP+ complex (Figure 4).
Overall these results illustrate the importance of maintaining both good local performance and large global tolerance and suggest the following design principle: , and . In order for Systemic Regime a to exist, should be less than 1. In biological terms, this can be achieved by selecting for a GSR whose apparent Michaelis constant for NADPH is lower than the total concentration of NADP. To avoid Systemic Regime b, which characterizes the worst local performance, should be larger than 1. In this case, G6PD should display an apparent Michaelis constant for NADP+ than is higher than the total concentration of NADP. To achieve the best Systemic Region (a) should be less than the ratio . In biological terms, the apparent maximal velocity of GSR should be much lower than that for G6PD.
As one final point, it seems that natural selection has produced a parasite, P. falciparum, that is capable of adapting to G6PD-deficient erythrocytes by converting a deficient NADPH redox cycle into one that has both robust local behavior and large global tolerance. In the particular study that we have analyzed , the host inherited a G6PD variant that leads to an operating point in Systemic Regime c. In order to adapt to this hostile environment, P. falciparum has to synthesize its own G6PD and is burdened with the associated energetic costs; the host cells, however, cannot benefit from this improved cycle because their ultimate fate is destruction by the P. falciparum.
Our analysis of the NADPH redox cycle in human erythrocytes provides a clear example of how system design space concepts serve as a framework to understand quantitative physiological and pathological behavior. It allowed us to define and characterize three different phenotypic regions: physiological (Systemic regime a), pathological (Systemic regime c) and potentially lethal (Systemic regime b). This example, and others , , , illustrates the utility of creating a design space that illuminates the function, design and fitness of the system. The system design space includes relationships among genotype, phenotype and environment. It provides an important framework to characterize the quantitative behavior of a system, and to compare wild-type and mutant variants. We know of no other approach that is able efficiently to accomplish these desirable objectives.
Steady-state analysis: steady-state solutions and analysis of local performance.
(0.71 MB PDF)
Annotated database of G6PD variants.
(0.55 MB PDF)
The authors would like to acknowledge Rick Fasani, Dean Tolla, Jason Lomnitz and Phillip Seitzer (University of California, Davis) for fruitful discussions.
Competing Interests: The authors have declared that no competing interests exist.
Funding: This work was supported in part by a fellowship to PMBMC (SFRH/BD/8304/2002) from Fundacao para a Ciencia e Tecnologia, Portugal, and by grants to AS (PTDC/QUI/70523/2006) from Fundacao para a Ciencia e Tecnologia (http://alfa.fct.mctes.pt), Portugal, and to MAS (RO1-GM30054) from the US Public Health Service (http://www.nigms.nih.gov) and a Stanislaw Ulam Distinguished Scholar Award from the Center for Non-Linear Studies of the Los Alamos National Laboratory (http://cnls.lanl.gov). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.