Evolution of enzymes in a series is driven by dissimilar functional demands

Design Principle for Decoding Calcium Signals to Generate Specific Gene Expression Via Transcription

The second messenger calcium plays a key role in conveying specificity of signaling pathways in plant cells. Specific calcium signatures are decoded to generate correct gene expression responses and amplification of calcium signatures is vital to this process. (1) It is not known if this amplification is an intrinsic property of all calcium-regulated gene expression responses and whether all calcium signatures have the potential to be amplified, or (2) how a given calcium signature maintains specificity in cells containing a great number of transcription factors (TFs) and other proteins with the potential to be calcium-regulated. The work presented here uncovers the design principle by which it is possible to decode calcium signals into specific changes in gene transcription in plant cells. Regarding the first question, we found that the binding mechanism between protein components possesses an intrinsic property that will nonlinearly amplify any calcium signal. This nonlinear amplification allows plant cells to effectively distinguish the kinetics of different calcium signatures to produce specific and appropriate changes in gene expression. Regarding the second question, we found that the large number of calmodulin (CaM)-binding TFs or proteins in plant cells form a buffering system such that the concentration of an active CaM-binding TF is insensitive to the concentration of any other CaM-binding protein, thus maintaining specificity. The design principle revealed by this work can be used to explain how any CaM-binding TF decodes calcium signals to generate specific gene expression responses in plant cells via transcription.

Mapping the phenotypic repertoire of the cytoplasmic 2-Cys peroxiredoxin - Thioredoxin system. 1. Understanding commonalities and differences among cell types

The system (PTTRS) formed by typical 2-Cys peroxiredoxins (Prx), thioredoxin (Trx), Trx reductase (TrxR), and sulfiredoxin (Srx) is central in antioxidant protection and redox signaling in the cytoplasm of eukaryotic cells. Understanding how the PTTRS integrates these functions requires tracing phenotypes to molecular properties, which is non-trivial. Here we analyze this problem based on a model that captures the PTTRS' conserved features. We have mapped the conditions that generate each distinct response to H2O2 supply rates (vsup), and estimated the parameters for thirteen human cell types and for Saccharomyces cerevisiae. The resulting composition-to-phenotype map yielded the following experimentally testable predictions. The PTTRS permits many distinct responses including ultra-sensitivity and hysteresis. However, nearly all tumor cell lines showed a similar response characterized by limited Trx-S- depletion and a substantial but self-limited gradual accumulation of hyperoxidized Prx at high vsup. This similarity ensues from strong correlations between the TrxR, Srx and Prx activities over cell lines, which contribute to maintain the Prx-SS reduction capacity in slight excess over the maximal steady state Prx-SS production. In turn, in erythrocytes, hepatocytes and HepG2 cells high vsup depletes Trx-S- and oxidizes Prx mainly to Prx-SS. In all nucleated human cells the Prx-SS reduction capacity defined a threshold separating two different regimes. At sub-threshold vsup the cytoplasmic H2O2 concentration is determined by Prx, nM-range and spatially localized, whereas at supra-threshold vsup it is determined by much less active alternative sinks and μM-range throughout the cytoplasm. The yeast shows a distinct response where the Prx Tsa1 accumulates in sulfenate form at high vsup. This is mainly due to an exceptional stability of Tsa1's sulfenate. The implications of these findings for thiol redox regulation and cell physiology are discussed. All estimates were thoroughly documented and provided, together with analytical approximations for system properties, as a resource for quantitative redox biology.

Design principles of autocatalytic cycles constrain enzyme kinetics and force low substrate saturation at flux branch points

A set of chemical reactions that require a metabolite to synthesize more of that metabolite is an autocatalytic cycle. Here, we show that most of the reactions in the core of central carbon metabolism are part of compact autocatalytic cycles. Such metabolic designs must meet specific conditions to support stable fluxes, hence avoiding depletion of intermediate metabolites. As such, they are subjected to constraints that may seem counter-intuitive: the enzymes of branch reactions out of the cycle must be overexpressed and the affinity of these enzymes to their substrates must be relatively weak. We use recent quantitative proteomics and fluxomics measurements to show that the above conditions hold for functioning cycles in central carbon metabolism of E. coli. This work demonstrates that the topology of a metabolic network can shape kinetic parameters of enzymes and lead to seemingly wasteful enzyme usage.

DOI:http://dx.doi.org/10.7554/eLife.20667.001

Many bacteria are able to produce all the molecules they need to survive from a limited supply of nutrients. This allows the bacteria to thrive even in harsh environments where other organisms struggle to live. The bacteria act as miniature chemical factories to convert nutrients into the desired molecules via a series of chemical reactions. Some molecules are made in sets of reactions termed autocatalytic cycles. These reaction sets require a molecule to be present in the cell in order to produce more of that molecule; like how a savings account needs to contain some money before it can generate more via interest.

Bacteria have many different enzymes that each drive specific chemical reactions. In order for an autocatalytic cycle to work properly, the cell needs to maintain adequate supplies of the molecule it is trying to make and all of the “intermediate” molecules in the cycle. If less of an intermediate molecule is produced, for example, the cell needs to reduce the demand for that molecule by controlling later chemical reactions in the cycle. Bacteria control chemical reactions by regulating the activities of the enzymes involved, but it is not clear exactly how they regulate the enzymes that drive autocatalytic cycles.

Barenholz et al. combined two approaches called proteomics and fluxomics to study autocatalytic cycles in a bacterium known as E. coli. The experiments suggest several core principles allow autocatalytic cycles to work smoothly in the bacteria. The next step is to apply these principles to different kinds of molecules produced in bacterial cells. A future challenge is to search for other structures that regulate chemical reactions in E. coli and other bacteria. Extending our understanding of autocatalytic cycles and other pathways of chemical reactions is essential for designing and engineering new reactions in bacteria. Such knowledge can be used to modify bacteria to produce valuable chemicals in environmentally friendly ways.

DOI:http://dx.doi.org/10.7554/eLife.20667.002

Hydrogen peroxide metabolism and sensing in human erythrocytes: a validated kinetic model and reappraisal of the role of peroxiredoxin II

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.

Biochemical Systems Theory: A Review

Biochemical systems theory (BST) is the foundation for a set of analytical andmodeling tools that facilitate the analysis of dynamic biological systems. This paper depicts major developments in BST up to the current state of the art in 2012. It discusses its rationale, describes the typical strategies and methods of designing, diagnosing, analyzing, and utilizing BST models, and reviews areas of application. The paper is intended as a guide for investigators entering the fascinating field of biological systems analysis and as a resource for practitioners and experts.

Relating mutant genotype to phenotype via quantitative behavior of the NADPH redox cycle in human erythrocytes

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.

Enzymes involved in osmolyte synthesis: how does oxidative stress affect osmoregulation in renal cells?

Kidney medulla cells are exposed to a wide range of changes in the ionic and osmotic composition of their environment as a consequence of the urine concentrating mechanism. During antidiuresis NaCl and urea concentrations increase and an efficient urinary concentrating mechanism is accompanied by medullar hypoxia. Medullar hypotonicity increases reactive oxygen species, a byproduct of mitochondria during ATP production. High intracellular ionic strength, hypoxia and elevated ROS concentration would have deleterious effects on medulla cell function. Medulla cells respond to hypertonicity by accumulating organic osmolytes, such as glycine betaine, glycerophosphorylcholine, sorbitol, inositol, and taurine, the main functions of which are osmoregulation and osmoprotection. The accumulation of compatible osmolytes is thus crucial for the viability of renal medulla cells. Studies about the effects of reactive oxygen species (ROS) on the enzymes involved in the synthesis of osmolytes are scarce. In this review we summarize the information available on the effects of ROS on the enzymes involved in osmolyte synthesis in kidney.

Evolutionary balancing of fitness-limiting factors

Debates concerning the roles of different factors that may limit an organism's reproductive success pervade evolutionary ecology. We suggest that a broad class of limiting-factors problems involving essential resources or essential components of reproductive effort can be analyzed with an evolutionary application of Liebig's law of the minimum. We explore life-history evolution using the metaphor of an organism that must harvest two essential resources (resources 1 and 2) from a stochastically varying environment. Our models make three predictions. First, organisms should overinvest, relative to the deterministic case, in harvesting the resource whose per-offspring harvest cost is smaller. Second, at the optimum, organisms balance multiple fitness-limiting factors rather than being consistently limited by one factor. Third, the optimal investment in harvesting a resource is directly linked to the probability that the organism's fitness will be limited by that resource. Under temporal variation, the optimal proportional investment in harvesting resource 1 is equal to the probability that resource 1 will limit fitness. Our results help to explain why the responses of populations to environmental perturbations are hard to predict: as an organism transitions between different limiting factors, its responses to perturbations of those factors will likewise change.

Identifying quantitative operation principles in metabolic pathways: a systematic method for searching feasible enzyme activity patterns leading to cellular adaptive responses

BACKGROUND

Optimization methods allow designing changes in a system so that specific goals are attained. These techniques are fundamental for metabolic engineering. However, they are not directly applicable for investigating the evolution of metabolic adaptation to environmental changes. Although biological systems have evolved by natural selection and result in well-adapted systems, we can hardly expect that actual metabolic processes are at the theoretical optimum that could result from an optimization analysis. More likely, natural systems are to be found in a feasible region compatible with global physiological requirements.

RESULTS

We first present a new method for globally optimizing nonlinear models of metabolic pathways that are based on the Generalized Mass Action (GMA) representation. The optimization task is posed as a nonconvex nonlinear programming (NLP) problem that is solved by an outer-approximation algorithm. This method relies on solving iteratively reduced NLP slave subproblems and mixed-integer linear programming (MILP) master problems that provide valid upper and lower bounds, respectively, on the global solution to the original NLP. The capabilities of this method are illustrated through its application to the anaerobic fermentation pathway in Saccharomyces cerevisiae. We next introduce a method to identify the feasibility parametric regions that allow a system to meet a set of physiological constraints that can be represented in mathematical terms through algebraic equations. This technique is based on applying the outer-approximation based algorithm iteratively over a reduced search space in order to identify regions that contain feasible solutions to the problem and discard others in which no feasible solution exists. As an example, we characterize the feasible enzyme activity changes that are compatible with an appropriate adaptive response of yeast Saccharomyces cerevisiae to heat shock

CONCLUSION

Our results show the utility of the suggested approach for investigating the evolution of adaptive responses to environmental changes. The proposed method can be used in other important applications such as the evaluation of parameter changes that are compatible with health and disease states.

Quantifying global tolerance of biochemical systems: design implications for moiety-transfer cycles

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.

The ability of organisms to survive under a multitude of conditions is readily apparent. This robustness in performance is difficult to precisely characterize and quantify. At a biochemical level, it leads to physiological behavior when the parameters of the system remain within some neighborhood of their normal values. However, this behavior can change abruptly, often becoming pathological, as the boundary of the neighborhood is crossed. Currently, there is no generic approach to identifying and characterizing such boundaries. In this paper, we address the problem by introducing a method that involves quantitative concepts for boundaries between regions and “global tolerance”. To illustrate the power of these concepts, we analyzed a large class of biological modules called moiety-transfer cycles and characterized the specific case of the NADPH redox cycle in human erythrocytes, which is involved in conferring resistance to malaria. Our results show that the wild-type system operates well within a region of “best” local performance that is surrounded by “poor” regions.

Evaluating dosage compensation as a cause of duplicate gene retention in Paramecium tetraurelia

Can dosage compensation completely explain gene retention after whole-genome duplication?

The high retention of duplicate genes in the genome of Paramecium tetraurelia has led to the hypothesis that most of the retained genes have persisted because of constraints due to gene dosage. This and other possible mechanisms are discussed in the light of expectations from population genetics and systems biology.

Simulation study of methemoglobin reduction in erythrocytes. Differential contributions of two pathways to tolerance to oxidative stress

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 metHb-reducing 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 metHb-reducing pathways in human erythrocytes.

A systematic analysis of lineage-specific evolution in metabolic pathways

In a search for the lineage-specific evolution of pathways between human, chimpanzee, mouse, and rat, orthologous gene families were generated from genome sequences. For each family, a model-based ratio of nonsynonymous to synonymous nucleotide substitution rates was calculated. Where the free-ratio model of individual ratios on each branch was supported, these families were mapped to two databases of metabolic pathways (KEGG and BioCyc) and the lineage-specific evolution of pathways was evaluated. The most similar pathway evolution was seen between mouse and rat, while the evolutionary pattern between human and chimpanzee was less correlated. Individual pathways in the human lineage were observed to evolve in a faster, lineage-specific manner, including the pathway involving arachidonic acid metabolism (identified through the KEGG analysis) and pyrimidine metabolism (identified through both analyses).

Flux control and excess capacity in the enzymes of glycolysis and their relationship to flight metabolism in Drosophila melanogaster

An important question in evolutionary and physiological genetics is how the control of flux-base phenotypes is distributed across the enzymes in a pathway. This control is often related to enzyme-specific levels of activity that are reported to be in excess of that required for demand. In glycolysis, metabolic control is frequently considered vested in classical regulatory enzymes, each strongly displaced from equilibrium. Yet the contribution of individual steps to control is unclear. To assess enzyme-specific control in the glycolytic pathway, we used P-element excision-derived mutagenesis in Drosophila melanogaster to generate full and partial knockouts of seven metabolic genes and to measure tethered flight performance. For most enzymes, we find that reduction to half of the normal activity has no measurable impact on wing beat frequency. The enzymes catalyzing near-equilibrium reactions, phosphoglucose isomerase, phosphoglucomutase, and triosephosphate isomerase fail to show any decline in flight performance even when activity levels are reduced to 17% or less. At reduced activities, the classic regulatory enzymes, hexokinase and glycogen phosphorylase, show significant drops in flight performance and are nearer to saturation. Our results show that flight performance is canalized or robust to the activity variation found in natural populations. Furthermore, enzymes catalyzing near-equilibrium reactions show strong genetic dominance down to low levels of activity. This implies considerable excess enzyme capacity for these enzymes.

Use of physiological constraints to identify quantitative design principles for gene expression in yeast adaptation to heat shock

Background

Understanding the relationship between gene expression changes, enzyme activity shifts, and the corresponding physiological adaptive response of organisms to environmental cues is crucial in explaining how cells cope with stress. For example, adaptation of yeast to heat shock involves a characteristic profile of changes to the expression levels of genes coding for enzymes of the glycolytic pathway and some of its branches. The experimental determination of changes in gene expression profiles provides a descriptive picture of the adaptive response to stress. However, it does not explain why a particular profile is selected for any given response.

Results

We used mathematical models and analysis of in silico gene expression profiles (GEPs) to understand how changes in gene expression correlate to an efficient response of yeast cells to heat shock. An exhaustive set of GEPs, matched with the corresponding set of enzyme activities, was simulated and analyzed. The effectiveness of each profile in the response to heat shock was evaluated according to relevant physiological and functional criteria. The small subset of GEPs that lead to effective physiological responses after heat shock was identified as the result of the tuning of several evolutionary criteria. The experimentally observed transcriptional changes in response to heat shock belong to this set and can be explained by quantitative design principles at the physiological level that ultimately constrain changes in gene expression.

Conclusion

Our theoretical approach suggests a method for understanding the combined effect of changes in the expression of multiple genes on the activity of metabolic pathways, and consequently on the adaptation of cellular metabolism to heat shock. This method identifies quantitative design principles that facilitate understating the response of the cell to stress.