Pedro M. B. M. Coelho (original) (raw)

Papers by Pedro M. B. M. Coelho

PLOS Computational Biology, 2009

Robustness of organisms is widely observed although difficult to precisely characterize. Performa... more Robustness of organisms is widely observed although difficult to precisely characterize. Performance can remain nearly constant within some neighborhood of the normal operating regime, leading to homeostasis, but then abruptly break down with pathological consequences beyond this neighborhood. Currently, there is no generic approach to identifying boundaries where local performance deteriorates abruptly, and this has hampered understanding of the molecular basis of biological robustness. Here we introduce a generic approach for characterizing boundaries between operational regimes based on the piecewise power-law representation of the system's components. This conceptual framework allows us to define ''global tolerance'' as the ratio between the normal value of a parameter and the value at such a boundary. We illustrate the utility of this concept for a class of moiety-transfer cycles, which is a widespread module in biology. Our results show a region of ''best'' local performance surrounded by ''poor'' regions; also, selection for improved local performance often pushes the operating values away from regime boundaries, thus increasing global tolerance. These predictions agree with experimental data from the reduced nicotinamide adenine dinucleotide phosphate (NADPH) redox cycle of human erythrocytes.

PLOS Computational Biology, 2009

Robustness of organisms is widely observed although difficult to precisely characterize. Performa... more Robustness of organisms is widely observed although difficult to precisely characterize. Performance can remain nearly constant within some neighborhood of the normal operating regime, leading to homeostasis, but then abruptly break down with pathological consequences beyond this neighborhood. Currently, there is no generic approach to identifying boundaries where local performance deteriorates abruptly, and this has hampered understanding of the molecular basis of biological robustness. Here we introduce a generic approach for characterizing boundaries between operational regimes based on the piecewise power-law representation of the system's components. This conceptual framework allows us to define ''global tolerance'' as the ratio between the normal value of a parameter and the value at such a boundary. We illustrate the utility of this concept for a class of moiety-transfer cycles, which is a widespread module in biology. Our results show a region of ''best'' local performance surrounded by ''poor'' regions; also, selection for improved local performance often pushes the operating values away from regime boundaries, thus increasing global tolerance. These predictions agree with experimental data from the reduced nicotinamide adenine dinucleotide phosphate (NADPH) redox cycle of human erythrocytes.

Proceedings of The National Academy of Sciences, 2009

One of the major unsolved problems of modern biology is deep understanding of the complex relatio... more One of the major unsolved problems of modern biology is deep understanding of the complex relationship between the information encoded in the genome of an organism and the phenotypic properties manifested by that organism. Fundamental advances must be made before we can begin to approach the goal of predicting the phenotypic consequences of a given mutation or an organism's response to a novel environmental challenge. Although this problem is often portrayed as if the task were to find a more or less direct link between genotypic and phenotypic levels, on closer examination the relationship is far more layered and complex. Although there are some intuitive notions of what is meant by phenotype at the level of the organism, it is far from clear what this term means at the biochemical level. We have described design principles that are readily revealed by representation of molecular systems in an appropriate design space. Here, we first describe a generic approach to the construction of such a design space in which qualitatively distinct phenotypes can be identified and counted. Second, we show how the boundaries between these phenotypic regions provide a method of characterizing a system's tolerance to large changes in the values of its parameters. Third, we illustrate the approach for one of the most basic modules of biochemical networks and describe an associated design principle. Finally, we discuss the scaling of this approach to large systems.

Robustness of organisms is widely observed although difficult to precisely characterize. Performa... more Robustness of organisms is widely observed although difficult to precisely characterize. Performance can remain nearly constant within some neighborhood of the normal operating regime, leading to homeostasis, but then abruptly break down with pathological consequences beyond this neighborhood. Currently, there is no generic approach to identifying boundaries where local performance deteriorates abruptly, and this has hampered understanding of the molecular basis of biological robustness. Here we introduce a generic approach for characterizing boundaries between operational regimes based on the piecewise power-law representation of the system's components. This conceptual framework allows us to define “global tolerance” as the ratio between the normal value of a parameter and the value at such a boundary. We illustrate the utility of this concept for a class of moiety-transfer cycles, which is a widespread module in biology. Our results show a region of “best” local performance surrounded by “poor” regions; also, selection for improved local performance often pushes the operating values away from regime boundaries, thus increasing global tolerance. These predictions agree with experimental data from the reduced nicotinamide adenine dinucleotide phosphate (NADPH) redox cycle of human erythrocytes.

Proceedings of the …, Jan 1, 2009

PLOS One, 2010

One of the major unsolved problems of modern biology is deep understanding of the complex relatio... more One of the major unsolved problems of modern biology is deep understanding of the complex relationship between the information encoded in the genome of an organism and the phenotypic properties manifested by that organism. Fundamental advances must be made before we can begin to approach the goal of predicting the phenotypic consequences of a given mutation or an organism's response to a novel environmental challenge. Although this problem is often portrayed as if the task were to find a more or less direct link between genotypic and phenotypic levels, on closer examination the relationship is far more layered and complex. Although there are some intuitive notions of what is meant by phenotype at the level of the organism, it is far from clear what this term means at the biochemical level. We have described design principles that are readily revealed by representation of molecular systems in an appropriate design space. Here, we first describe a generic approach to the construction of such a design space in which qualitatively distinct phenotypes can be identified and counted. Second, we show how the boundaries between these phenotypic regions provide a method of characterizing a system's tolerance to large changes in the values of its parameters. Third, we illustrate the approach for one of the most basic modules of biochemical networks and describe an associated design principle. Finally, we discuss the scaling of this approach to large systems.

Free radical biology & medicine, 2014

Hydrogen peroxide (H 2 O 2 ) metabolism in human erythrocytes has been thoroughly investigated, b... more Hydrogen peroxide (H 2 O 2 ) metabolism in human erythrocytes has been thoroughly investigated, but unclear points persist. By integrating the available data into a mathematical model that accurately represents the current understanding and comparing computational predictions to observations we sought to (a) identify inconsistencies in present knowledge, (b) propose resolutions, and (c) examine their functional implications. The systematic confrontation of computational predictions with experimental observations of the responses of intact erythrocytes highlighted the following important discrepancy. The high rate constant (10 7 -10 8 M À 1 s À 1 ) for H 2 O 2 reduction determined for purified peroxiredoxin II (Prx2) and the high abundance of this protein indicate that under physiological conditions it consumes practically all the H 2 O 2 . However, this is inconsistent with extensive evidence that Prx2's contribution to H 2 O 2 elimination is comparable to that of catalase. Models modified such that Prx2's effective peroxidase activity is just 10 5 M À 1 s À 1 agree near quantitatively with extensive experimental observations. This low effective activity is probably due to a strong but readily reversible inhibition of Prx2's peroxidatic activity in intact cells, implying that the main role of Prx2 in human erythrocytes is not to eliminate peroxide substrates. Simulations of the responses to physiological H 2 O 2 stimuli highlight that a design combining abundant Prx2 with a low effective peroxidase activity spares NADPH while improving potential signaling properties of the Prx2/thioredoxin/thioredoxin reductase system.

Free Radical Biology and Medicine, 2014

Hydrogen peroxide (H2O2) metabolism in human erythrocytes has been thoroughly investigated, but u... more Hydrogen peroxide (H2O2) metabolism in human erythrocytes has been thoroughly investigated, but unclear points persist. By integrating the available data into a mathematical model that accurately represents the current understanding and comparing computational predictions to observations we sought to (a) identify inconsistencies in present knowledge, (b) propose resolutions, and (c) examine their functional implications. The systematic confrontation of computational predictions with experimental observations of the responses of intact erythrocytes highlighted the following important discrepancy. The high rate constant (10(7)-10(8) M(-1) s(-1)) for H2O2 reduction determined for purified peroxiredoxin II (Prx2) and the high abundance of this protein indicate that under physiological conditions it consumes practically all the H2O2. However, this is inconsistent with extensive evidence that Prx2's contribution to H2O2 elimination is comparable to that of catalase. Models modified such that Prx2's effective peroxidase activity is just 10(5) M(-1) s(-1) agree near quantitatively with extensive experimental observations. This low effective activity is probably due to a strong but readily reversible inhibition of Prx2's peroxidatic activity in intact cells, implying that the main role of Prx2 in human erythrocytes is not to eliminate peroxide substrates. Simulations of the responses to physiological H2O2 stimuli highlight that a design combining abundant Prx2 with a low effective peroxidase activity spares NADPH while improving potential signaling properties of the Prx2/thioredoxin/thioredoxin reductase system.

Free Radical Biology and Medicine, 2014

H2O2 elimination in human erythrocytes is mainly carried out by catalase (Cat), glutathione perox... more H2O2 elimination in human erythrocytes is mainly carried out by catalase (Cat), glutathione peroxidase (GPx1) and the more recently discovered peroxiredoxin 2 (Prx2). However, the contribution of Prx2 to H2O2 consumption is still unclear. Prx2's high reactivity with H2O2 (kPrx2=10×10(7) M(-1)s(-1), kCat =7×10(7) M(-1)s(-1), kGPx1 =4×10(7) M(-1)s(-1)) and high abundance ([Prx2]= 570µM, [Cat]= 32µM, [GPx1]= 1µM) suggest that under low H2O2 supply rates it should consume >99% of the H2O2. However, extensive evidence indicates that in intact erythrocytes Prx2 contributes no more than Cat to H2O2 consumption. In order for this to be attained, Prx2's effective rate constant with H2O2would have to be just ~10(5) M(-1)s(-1), much lower than that determined in multiple experiments with the purified proteins. Nevertheless, nearly all Prx2 is oxidized within 1min of exposing erythrocytes to a H2O2 bolus, which is inconsistent with an irreversible inhibition. A mathematical model of the H2O2 metabolism in human erythrocytes [Benfeitas et al. (2014) Free Radic. Biol. Med.] where Prx2 either has a low kPrx2 or is subject to a strong (>99%) but readily reversible inhibition achieves quantitative agreement with detailed experimental observations of the responses of the redox status of Prx2 in human erythrocytes and suggests functional advantages of this design (see companion abstract). By contrast, a variant where Prx2 is fully active with kPrx2=10(8) M(-1)s(-1) shows important qualitative discrepancies. Altogether, these results suggest that Prx2's peroxidase activity is strongly inhibited in human erythrocytes. We acknowledge fellowship SFRH/BD/51199/2010, grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, PEst-OE/QUI/UI0313/2014, and FCOMP-01-0124-FEDER-020978 (PTDC/QUI-BIQ/119657/2010) co-financed by FEDER through the COMPETE program and by FCT.

Free Radical Biology and Medicine, 2014

Cells are occasionally exposed to high H2O2 concentrations, often preceding exposure to other ele... more Cells are occasionally exposed to high H2O2 concentrations, often preceding exposure to other electrophylic compounds. Both H2O2 and these compounds can irreversibly modify protein thiols, with deleterious consequences. Induction of enzymatic defenses against those agents is too slow to avoid significant damage. Cells may solve this conundrum by reversibly "blocking" the thiols once H2O2 concentrations begin to increase. We term this mechanism "anticipatory blocking" because it acts in anticipation of irreversible damage upon detection of early signs of stress. Here we examine the design requirements for the Peroxiredoxin/Thioredoxin/Thioredoxin-Reductase/Protein-Dithiol System (PTTRDS) to effectively integrate H2O2 signaling and anticipatory blocking of protein dithiols as disulfides, and we compared them to the designs found in cells. To that effect, we developed a minimal model of the PTTRDS, and we defined a set of quantitative performance criteria that embody the requirements for (a) efficient scavenging capacity, (b) low NADPH consumption, (c) effective signal propagation, and (d) effective anticipatory blocking. We then sought the design principles (relationships among rate constants and species concentrations) that warrant fulfillment of all these criteria. Experimental data indicates that the design of the PTTRDS in human erythrocytes fulfills these principles and thus accomplishes effective integration between anticipatory blocking, antioxidant protection and redox signaling. A more general analysis suggests that the same principles hold in a wide variety of cell types and organisms. We acknowledge grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, FCOMP-01-0124-FEDER-020978 (PTDC/QUI-BIQ/119657/2010) financed by FEDER through the "Programa Operacional Factores de Competitividade, COMPETE" and by national funds through "FCT, Fundação para a Ciência e a Tecnologia".

Free Radical Biology and Medicine, 2014

In human erythrocytes H2O2 is mainly consumed by glutathione peroxidase, catalase and peroxiredox... more In human erythrocytes H2O2 is mainly consumed by glutathione peroxidase, catalase and peroxiredoxin 2 (Prx2). Our previous analyses indicate that Prx2's peroxidase activity is subjected to a strong but quickly reversible inhibition (see companion abstract). If this activity is inhibited then the main role of Prx2 cannot be to eliminate H2O2. What functional advantages could then such an inhibition confer?We set up and validated a kinetic model of H2O2 metabolism human erythrocytes that shows quantitative agreement with extensive experimental observations. We then applied it to analyze the behavior of Prx2 and Trx under the H2O2 exposure dynamics that erythrocytes face in circulation. The significance of Prx2 inhibition was assessed by comparing the behavior of this model with that of an otherwise identical model lacking inhibition.Our analysis shows that Prx2 inhibition leads to 25-40% lower NADPH consumption under low to moderately high H2O2 supply (<0.8µM H2O2/s). Further, the inhibition extends the range where the concentrations of potential redox signaling readouts - H2O2, Prx2 sulfenic acid, Prx2 disulfide and Trx disulfide- show a proportional response to changes in H2O2 supply, covering practically the whole physiological range of the latter. This is desirable for analogic signal transduction and allows the Prx2/Trx/TrxR system to reliably transduce changes in H2O2 supply as changes in thiol oxidation. Finally, the inhibition allows other less abundant peroxiredoxins in the erythrocyte to be oxidized by H2O2 at physiological H2O2 supplies.Altogether, these results suggest that the postulated reversible inhibition of Prx2's peroxidase facilitates signal transduction by the peroxiredoxins and spares NADPH.We acknowledge: fellowship SFRH/BD/51199/2010, grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, PEst-OE/QUI/UI0313/2014, and FCOMP-01-0124-FEDER-020978 co-financed by FEDER through the COMPETE program and by FCT (project PTDC/QUI-BIQ/119657/2010).

Background: The NADPH redox cycle plays a key role in antioxidant protection of human erythrocyte... more Background: 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.

Background The NADPH redox cycle plays a key role in antioxidant protection of human erythrocytes... more Background
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.

Methodology/Principal Findings
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.

Conclusions/Significance
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.

PLOS Computational Biology, 2009

Robustness of organisms is widely observed although difficult to precisely characterize. Performa... more Robustness of organisms is widely observed although difficult to precisely characterize. Performance can remain nearly constant within some neighborhood of the normal operating regime, leading to homeostasis, but then abruptly break down with pathological consequences beyond this neighborhood. Currently, there is no generic approach to identifying boundaries where local performance deteriorates abruptly, and this has hampered understanding of the molecular basis of biological robustness. Here we introduce a generic approach for characterizing boundaries between operational regimes based on the piecewise power-law representation of the system's components. This conceptual framework allows us to define ''global tolerance'' as the ratio between the normal value of a parameter and the value at such a boundary. We illustrate the utility of this concept for a class of moiety-transfer cycles, which is a widespread module in biology. Our results show a region of ''best'' local performance surrounded by ''poor'' regions; also, selection for improved local performance often pushes the operating values away from regime boundaries, thus increasing global tolerance. These predictions agree with experimental data from the reduced nicotinamide adenine dinucleotide phosphate (NADPH) redox cycle of human erythrocytes.

PLOS Computational Biology, 2009

Robustness of organisms is widely observed although difficult to precisely characterize. Performa... more Robustness of organisms is widely observed although difficult to precisely characterize. Performance can remain nearly constant within some neighborhood of the normal operating regime, leading to homeostasis, but then abruptly break down with pathological consequences beyond this neighborhood. Currently, there is no generic approach to identifying boundaries where local performance deteriorates abruptly, and this has hampered understanding of the molecular basis of biological robustness. Here we introduce a generic approach for characterizing boundaries between operational regimes based on the piecewise power-law representation of the system's components. This conceptual framework allows us to define ''global tolerance'' as the ratio between the normal value of a parameter and the value at such a boundary. We illustrate the utility of this concept for a class of moiety-transfer cycles, which is a widespread module in biology. Our results show a region of ''best'' local performance surrounded by ''poor'' regions; also, selection for improved local performance often pushes the operating values away from regime boundaries, thus increasing global tolerance. These predictions agree with experimental data from the reduced nicotinamide adenine dinucleotide phosphate (NADPH) redox cycle of human erythrocytes.

Proceedings of The National Academy of Sciences, 2009

One of the major unsolved problems of modern biology is deep understanding of the complex relatio... more One of the major unsolved problems of modern biology is deep understanding of the complex relationship between the information encoded in the genome of an organism and the phenotypic properties manifested by that organism. Fundamental advances must be made before we can begin to approach the goal of predicting the phenotypic consequences of a given mutation or an organism's response to a novel environmental challenge. Although this problem is often portrayed as if the task were to find a more or less direct link between genotypic and phenotypic levels, on closer examination the relationship is far more layered and complex. Although there are some intuitive notions of what is meant by phenotype at the level of the organism, it is far from clear what this term means at the biochemical level. We have described design principles that are readily revealed by representation of molecular systems in an appropriate design space. Here, we first describe a generic approach to the construction of such a design space in which qualitatively distinct phenotypes can be identified and counted. Second, we show how the boundaries between these phenotypic regions provide a method of characterizing a system's tolerance to large changes in the values of its parameters. Third, we illustrate the approach for one of the most basic modules of biochemical networks and describe an associated design principle. Finally, we discuss the scaling of this approach to large systems.

Robustness of organisms is widely observed although difficult to precisely characterize. Performa... more Robustness of organisms is widely observed although difficult to precisely characterize. Performance can remain nearly constant within some neighborhood of the normal operating regime, leading to homeostasis, but then abruptly break down with pathological consequences beyond this neighborhood. Currently, there is no generic approach to identifying boundaries where local performance deteriorates abruptly, and this has hampered understanding of the molecular basis of biological robustness. Here we introduce a generic approach for characterizing boundaries between operational regimes based on the piecewise power-law representation of the system's components. This conceptual framework allows us to define “global tolerance” as the ratio between the normal value of a parameter and the value at such a boundary. We illustrate the utility of this concept for a class of moiety-transfer cycles, which is a widespread module in biology. Our results show a region of “best” local performance surrounded by “poor” regions; also, selection for improved local performance often pushes the operating values away from regime boundaries, thus increasing global tolerance. These predictions agree with experimental data from the reduced nicotinamide adenine dinucleotide phosphate (NADPH) redox cycle of human erythrocytes.

Proceedings of the …, Jan 1, 2009

PLOS One, 2010

One of the major unsolved problems of modern biology is deep understanding of the complex relatio... more One of the major unsolved problems of modern biology is deep understanding of the complex relationship between the information encoded in the genome of an organism and the phenotypic properties manifested by that organism. Fundamental advances must be made before we can begin to approach the goal of predicting the phenotypic consequences of a given mutation or an organism's response to a novel environmental challenge. Although this problem is often portrayed as if the task were to find a more or less direct link between genotypic and phenotypic levels, on closer examination the relationship is far more layered and complex. Although there are some intuitive notions of what is meant by phenotype at the level of the organism, it is far from clear what this term means at the biochemical level. We have described design principles that are readily revealed by representation of molecular systems in an appropriate design space. Here, we first describe a generic approach to the construction of such a design space in which qualitatively distinct phenotypes can be identified and counted. Second, we show how the boundaries between these phenotypic regions provide a method of characterizing a system's tolerance to large changes in the values of its parameters. Third, we illustrate the approach for one of the most basic modules of biochemical networks and describe an associated design principle. Finally, we discuss the scaling of this approach to large systems.

Free radical biology & medicine, 2014

Hydrogen peroxide (H 2 O 2 ) metabolism in human erythrocytes has been thoroughly investigated, b... more Hydrogen peroxide (H 2 O 2 ) metabolism in human erythrocytes has been thoroughly investigated, but unclear points persist. By integrating the available data into a mathematical model that accurately represents the current understanding and comparing computational predictions to observations we sought to (a) identify inconsistencies in present knowledge, (b) propose resolutions, and (c) examine their functional implications. The systematic confrontation of computational predictions with experimental observations of the responses of intact erythrocytes highlighted the following important discrepancy. The high rate constant (10 7 -10 8 M À 1 s À 1 ) for H 2 O 2 reduction determined for purified peroxiredoxin II (Prx2) and the high abundance of this protein indicate that under physiological conditions it consumes practically all the H 2 O 2 . However, this is inconsistent with extensive evidence that Prx2's contribution to H 2 O 2 elimination is comparable to that of catalase. Models modified such that Prx2's effective peroxidase activity is just 10 5 M À 1 s À 1 agree near quantitatively with extensive experimental observations. This low effective activity is probably due to a strong but readily reversible inhibition of Prx2's peroxidatic activity in intact cells, implying that the main role of Prx2 in human erythrocytes is not to eliminate peroxide substrates. Simulations of the responses to physiological H 2 O 2 stimuli highlight that a design combining abundant Prx2 with a low effective peroxidase activity spares NADPH while improving potential signaling properties of the Prx2/thioredoxin/thioredoxin reductase system.

Free Radical Biology and Medicine, 2014

Hydrogen peroxide (H2O2) metabolism in human erythrocytes has been thoroughly investigated, but u... more Hydrogen peroxide (H2O2) metabolism in human erythrocytes has been thoroughly investigated, but unclear points persist. By integrating the available data into a mathematical model that accurately represents the current understanding and comparing computational predictions to observations we sought to (a) identify inconsistencies in present knowledge, (b) propose resolutions, and (c) examine their functional implications. The systematic confrontation of computational predictions with experimental observations of the responses of intact erythrocytes highlighted the following important discrepancy. The high rate constant (10(7)-10(8) M(-1) s(-1)) for H2O2 reduction determined for purified peroxiredoxin II (Prx2) and the high abundance of this protein indicate that under physiological conditions it consumes practically all the H2O2. However, this is inconsistent with extensive evidence that Prx2's contribution to H2O2 elimination is comparable to that of catalase. Models modified such that Prx2's effective peroxidase activity is just 10(5) M(-1) s(-1) agree near quantitatively with extensive experimental observations. This low effective activity is probably due to a strong but readily reversible inhibition of Prx2's peroxidatic activity in intact cells, implying that the main role of Prx2 in human erythrocytes is not to eliminate peroxide substrates. Simulations of the responses to physiological H2O2 stimuli highlight that a design combining abundant Prx2 with a low effective peroxidase activity spares NADPH while improving potential signaling properties of the Prx2/thioredoxin/thioredoxin reductase system.

Free Radical Biology and Medicine, 2014

H2O2 elimination in human erythrocytes is mainly carried out by catalase (Cat), glutathione perox... more H2O2 elimination in human erythrocytes is mainly carried out by catalase (Cat), glutathione peroxidase (GPx1) and the more recently discovered peroxiredoxin 2 (Prx2). However, the contribution of Prx2 to H2O2 consumption is still unclear. Prx2's high reactivity with H2O2 (kPrx2=10×10(7) M(-1)s(-1), kCat =7×10(7) M(-1)s(-1), kGPx1 =4×10(7) M(-1)s(-1)) and high abundance ([Prx2]= 570µM, [Cat]= 32µM, [GPx1]= 1µM) suggest that under low H2O2 supply rates it should consume >99% of the H2O2. However, extensive evidence indicates that in intact erythrocytes Prx2 contributes no more than Cat to H2O2 consumption. In order for this to be attained, Prx2's effective rate constant with H2O2would have to be just ~10(5) M(-1)s(-1), much lower than that determined in multiple experiments with the purified proteins. Nevertheless, nearly all Prx2 is oxidized within 1min of exposing erythrocytes to a H2O2 bolus, which is inconsistent with an irreversible inhibition. A mathematical model of the H2O2 metabolism in human erythrocytes [Benfeitas et al. (2014) Free Radic. Biol. Med.] where Prx2 either has a low kPrx2 or is subject to a strong (>99%) but readily reversible inhibition achieves quantitative agreement with detailed experimental observations of the responses of the redox status of Prx2 in human erythrocytes and suggests functional advantages of this design (see companion abstract). By contrast, a variant where Prx2 is fully active with kPrx2=10(8) M(-1)s(-1) shows important qualitative discrepancies. Altogether, these results suggest that Prx2's peroxidase activity is strongly inhibited in human erythrocytes. We acknowledge fellowship SFRH/BD/51199/2010, grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, PEst-OE/QUI/UI0313/2014, and FCOMP-01-0124-FEDER-020978 (PTDC/QUI-BIQ/119657/2010) co-financed by FEDER through the COMPETE program and by FCT.

Free Radical Biology and Medicine, 2014

Cells are occasionally exposed to high H2O2 concentrations, often preceding exposure to other ele... more Cells are occasionally exposed to high H2O2 concentrations, often preceding exposure to other electrophylic compounds. Both H2O2 and these compounds can irreversibly modify protein thiols, with deleterious consequences. Induction of enzymatic defenses against those agents is too slow to avoid significant damage. Cells may solve this conundrum by reversibly "blocking" the thiols once H2O2 concentrations begin to increase. We term this mechanism "anticipatory blocking" because it acts in anticipation of irreversible damage upon detection of early signs of stress. Here we examine the design requirements for the Peroxiredoxin/Thioredoxin/Thioredoxin-Reductase/Protein-Dithiol System (PTTRDS) to effectively integrate H2O2 signaling and anticipatory blocking of protein dithiols as disulfides, and we compared them to the designs found in cells. To that effect, we developed a minimal model of the PTTRDS, and we defined a set of quantitative performance criteria that embody the requirements for (a) efficient scavenging capacity, (b) low NADPH consumption, (c) effective signal propagation, and (d) effective anticipatory blocking. We then sought the design principles (relationships among rate constants and species concentrations) that warrant fulfillment of all these criteria. Experimental data indicates that the design of the PTTRDS in human erythrocytes fulfills these principles and thus accomplishes effective integration between anticipatory blocking, antioxidant protection and redox signaling. A more general analysis suggests that the same principles hold in a wide variety of cell types and organisms. We acknowledge grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, FCOMP-01-0124-FEDER-020978 (PTDC/QUI-BIQ/119657/2010) financed by FEDER through the "Programa Operacional Factores de Competitividade, COMPETE" and by national funds through "FCT, Fundação para a Ciência e a Tecnologia".

Free Radical Biology and Medicine, 2014

In human erythrocytes H2O2 is mainly consumed by glutathione peroxidase, catalase and peroxiredox... more In human erythrocytes H2O2 is mainly consumed by glutathione peroxidase, catalase and peroxiredoxin 2 (Prx2). Our previous analyses indicate that Prx2's peroxidase activity is subjected to a strong but quickly reversible inhibition (see companion abstract). If this activity is inhibited then the main role of Prx2 cannot be to eliminate H2O2. What functional advantages could then such an inhibition confer?We set up and validated a kinetic model of H2O2 metabolism human erythrocytes that shows quantitative agreement with extensive experimental observations. We then applied it to analyze the behavior of Prx2 and Trx under the H2O2 exposure dynamics that erythrocytes face in circulation. The significance of Prx2 inhibition was assessed by comparing the behavior of this model with that of an otherwise identical model lacking inhibition.Our analysis shows that Prx2 inhibition leads to 25-40% lower NADPH consumption under low to moderately high H2O2 supply (<0.8µM H2O2/s). Further, the inhibition extends the range where the concentrations of potential redox signaling readouts - H2O2, Prx2 sulfenic acid, Prx2 disulfide and Trx disulfide- show a proportional response to changes in H2O2 supply, covering practically the whole physiological range of the latter. This is desirable for analogic signal transduction and allows the Prx2/Trx/TrxR system to reliably transduce changes in H2O2 supply as changes in thiol oxidation. Finally, the inhibition allows other less abundant peroxiredoxins in the erythrocyte to be oxidized by H2O2 at physiological H2O2 supplies.Altogether, these results suggest that the postulated reversible inhibition of Prx2's peroxidase facilitates signal transduction by the peroxiredoxins and spares NADPH.We acknowledge: fellowship SFRH/BD/51199/2010, grants PEst-C/SAU/LA0001/2013-2014, PEst-OE/QUI/UI0612/2013, PEst-OE/QUI/UI0313/2014, and FCOMP-01-0124-FEDER-020978 co-financed by FEDER through the COMPETE program and by FCT (project PTDC/QUI-BIQ/119657/2010).

Background: The NADPH redox cycle plays a key role in antioxidant protection of human erythrocyte... more Background: 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.

Background The NADPH redox cycle plays a key role in antioxidant protection of human erythrocytes... more Background
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.

Methodology/Principal Findings
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.

Conclusions/Significance
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.