Control Analysis of Cooperativity and Complementarity in Metabolic Regulations: The Case of NADPH Homeostasis.

Complex feedback regulation patterns shape the cellular metabolic response to external or internal perturbations. We propose here a framework consisting of a sampling-based metabolic control analysis of kinetic models to investigate the modes of regulatory interplay in metabolic functions. NADPH homeostasis, for instance in a context of oxidative stress, is an example of metabolic function that involves multiple feedback regulations which raises the issue of their concerted action. Our computational framework allows us to characterize both respective and combined effects of regulations, distinguishing between synergistic versus complementary modes of regulatory crosstalk. Synergistic regulation of G6PD enzymes and PGI enzymes is mediated by congruent effects between concentration sensitivities and reaction elasticities. Complementary regulation of pentose phosphate pathway and lower glycolysis relates to metabolic state-dependent range of regulation efficiency. These cooperative effects are shown to significantly improve metabolic flux response to support NADPH homeostasis, providing a rationale for the complex feedback regulation pattern at work.

Stochastic simulation algorithm for isotope-based dynamic flux analysis.

Carbon isotope labeling method is a standard metabolic engineering tool for flux quantification in living cells. To cope with the high dimensionality of isotope labeling systems, diverse algorithms have been developed to reduce the number of variables or operations in metabolic flux analysis (MFA), but lacks generalizability to non-stationary metabolic conditions. In this study, we present a stochastic simulation algorithm (SSA) derived from the chemical master equation of the isotope labeling system. This algorithm allows to compute the time evolution of isotopomer concentrations in non-stationary conditions, with the valuable property that computational time does not scale with the number of isotopomers. The efficiency and limitations of the algorithm is benchmarked for the forward and inverse problems of 13C-DMFA in the pentose phosphate pathways, and is compared with EMU-based methods for NMFA and MFA including the central carbon metabolism. Overall, SSA constitutes an alternative class to deterministic approaches for metabolic flux analysis that is well adapted to comprehensive dataset including parallel labeling experiments, and whose limitations associated to the sampling size can be overcome by using Monte Carlo sampling approaches.

Quantitative modeling of pentose phosphate pathway response to oxidative stress reveals a cooperative regulatory strategy.

Living cells use signaling and regulatory mechanisms to adapt to environmental stresses. Adaptation to oxidative stress involves the regulation of many enzymes in both glycolysis and pentose phosphate pathways (PPP), so as to support PPP-driven NADPH recycling for antioxidant defense. The underlying regulatory logic is investigated by developing a kinetic modeling approach fueled with metabolomics and 13C-fluxomics datasets from human fibroblast cells. Bayesian parameter estimation and phenotypic analysis of models highlight complementary roles for several metabolite-enzyme regulations. Specifically, carbon flux rerouting into PPP involves a tight coordination between the upregulation of G6PD activity concomitant to a decreased NADPH/NADP+ ratio and the differential control of downward and upward glycolytic fluxes through the joint inhibition of PGI and GAPD enzymes. Such functional interplay between distinct regulatory feedbacks promotes efficient detoxification and homeostasis response over a broad range of stress level, but can also explain paradoxical pertubation phenotypes for instance reported for 6PGD modulation in mammalian cells.

Theoretical study of the impact of adaptation on cell-fate heterogeneity and fractional killing.

Fractional killing illustrates the cell propensity to display a heterogeneous fate response over a wide range of stimuli. The interplay between the nonlinear and stochastic dynamics of biochemical networks plays a fundamental role in shaping this probabilistic response and in reconciling requirements for heterogeneity and controllability of cell-fate decisions. The stress-induced fate choice between life and death depends on an early adaptation response which may contribute to fractional killing by amplifying small differences between cells. To test this hypothesis, we consider a stochastic modeling framework suited for comprehensive sensitivity analysis of dose response curve through the computation of a fractionality index. Combining bifurcation analysis and Langevin simulation, we show that adaptation dynamics enhances noise-induced cell-fate heterogeneity by shifting from a saddle-node to a saddle-collision transition scenario. The generality of this result is further assessed by a computational analysis of a detailed regulatory network model of apoptosis initiation and by a theoretical analysis of stochastic bifurcation mechanisms. Overall, the present study identifies a cooperative interplay between stochastic, adaptation and decision intracellular processes that could promote cell-fate heterogeneity in many contexts.