Application of the transition time of metabolic system as a criterion for optimization of metabolic processes

Transition time of metabolic systems in introduced as a suitable optimization criterion for biotechnological processes in which it is desirable to reduce the lag time and minimize the mass contained within the system. Lag time is the time needed for the system to attain the steady state. Results obtained from the sensitivity analysis of this steady state response are presented within the metabolic control analysis and applied to 3 case studies. In all of them the information provided by the transition time control profile allows the implementation of a strategy for biotechnological manipulations aimed at the improvement of the process. (c) 1994 John Wiley & Sons, Inc.     

Experimental determination of flux control distribution in biochemical systems: in vitro model to analyze transient metabolite concentrations

Determination of the control coefficients allows the identification of rate-controlling steps in a reaction system. However, the measurement of the flux control coefficients in a biochemical system is not a trivial task, except for some special cases. We have developed a theoretical basis for the direct determination of these coefficients from dynamic responses. In order to show the validity of this methodology experimentally, the dynamic approach is applied to an in vitro reconstituted partial glycolytic pathway to determine the flux control coefficients of hexokinase and phosphofructokinase. It is shown that the dynamic approach gives consistent results, which agree well with values obtained by the direct enzyme titration method. The detailed procedure and potential applications to other systems, such as immobilized enzyme or cell reactors, are discussed. (c) 1993 Wiley & Sons, Inc.     

Modelization and experimental studies on the control of the glycolytic-glycogenolytic pathway in rat liver

In this paper we construct a model of the glycolytic-glycogenolytic converging pathway in rat liver, by integrating experimental data obtained in anin vitro system and information available from the literature. The model takes the mathematical expression of an S-system representation within the power law formalism (Savageau, 1976. Biochemical System Analysis: A study of function and design in Molecular Biology. Addison-Wesley, Reading, Mass.). By using this theoretical framework a model analysis was carried out that allowed us a) the assessment of the quality of the model in terms of its consistency and robustness, b) the steady state analysis and control characterization of the system, and c) the study of the dynamics of the system after changes in the level of two magnitudes of biological significance: the glucose concentration and the phosphofructokinase enzyme activity. Model predictions are compared with experimental measurements referred to Logarithmic Gains through fluxes and substrates concentrations showing that there is a good correlation between the model predictions and the experimentally determined values.     

Control analysis of transition times

Summary The present theoretical basis of Control Analysis is extended with the definition of Transition Time Response Coefficients. Some new relationships between local and global coefficients defined in Control Analysis are presented. These relationships are in the form of matrix products constructed in a priori form. The use of these straightforward relationships is shown in an exemplary application corresponding to an experimental system consisting of the glycolytic degradation from glucose to glyceraldehyde-3-phosphate.     

In situ metabolic flux analysis to quantify the liver metabolic response to experimental burn injury

Trauma such as burns induces a hypermetabolic response associated with altered central carbon and nitrogen metabolism. The liver plays a key role in these metabolic changes; however, studies to date have evaluated the metabolic state of liver using ex vivo perfusions or isotope labeling techniques targeted to specific pathways. Herein, we developed a unique mass balance approach to characterize the metabolic state of the liver in situ, and used it to quantify the metabolic changes to experimental burn injury in rats. Rats received a sham (control uninjured), 20% or 40% total body surface area (TBSA) scald burn, and were allowed to develop a hypermetabolic response. One day prior to evaluation, all animals were fasted to deplete glycogen stores. Four days post-burn, blood flow rates in major vessels of the liver were measured, and blood samples harvested. We combined measurements of metabolite concentrations and flow rates in the major vessels entering and leaving the liver with a steady-state mass balance model to generate a quantitative picture of the metabolic state of liver. The main findings were: (1) Sham-burned animals exhibited a gluconeogenic pattern, consistent with the fasted state; (2) the 20% TBSA burn inhibited gluconeogenesis and exhibited glycolytic-like features with very few other significant changes; (3) the 40% TBSA burn, by contrast, further enhanced gluconeogenesis and also increased amino acid extraction, urea cycle reactions, and several reactions involved in oxidative phosphorylation. These results suggest that increasing the severity of injury does not lead to a simple dose-dependent metabolic response, but rather leads to qualitatively different responses.     

Effect of precise control of flux ratio between the glycolytic pathways on mevalonate production in Escherichia coli

Mevalonate is a useful metabolite synthesized from three molecules of acetyl-CoA, consuming two molecules of NADPH. Escherichia coli ( E. coli) catabolizes glucose to acetyl-CoA via several routes, such as the Embden–Meyerhof–Parnas (EMP) and the oxidative pentose phosphate (oxPP) pathways. Although the oxPP pathway supplies NADPH, it is disadvantageous in terms of acetyl-CoA supply, compared with the EMP pathway. In this study, the optimal flux ratio between the EMP and oxPP pathways on the mevalonate yield was investigated. Expression level of pgi was controlled by isopropyl β-D-1-thiogalactopyranoside (IPTG) inducible promoter in an engineered mevalonate-producing E. coli strain. The relationship between the flux ratio and mevalonate yield was evaluated by changing the flux ratio by varying IPTG concentration. At the stationary phase, the mevalonate yield was maximum at an EMP flux of 39.7%, and was increased by 25% compared with that with no flux control (EMP flux of 70.4%). The optimal flux ratio was consistent with the theoretical value based on the mass balance of NADPH. The flux ratio between EMP and oxPP pathways affects the synthesis fluxes of mevalonate and acetate from acetyl-CoA. Fine tuning of the flux ratio would be necessary to achieve an optimized production of metabolites that require NADPH.     

代谢流量比率分析及其在代谢工程中的应用

细胞内代谢反应流量在系统理解细胞代谢特性和指导代谢工程改造等方面都起着重要的作用。由于代谢流量难以直接测量得到,在很多情况下通过跟踪稳定同位素在代谢网络中的转移并进行相应的模型计算能有效地定量代谢流量。代谢流量比率分析法能够高度体现系统的生物化学真实性、辨别细胞代谢网络的拓扑结构,并且能够相对简单快速地定量反应速率等,因此受到代谢工程研究者越来越多的重视。以下着重介绍并讨论了利用代谢物同位体分布信息分析关键代谢节点合成途径的流量比率、基于流量比率的代谢流量解析、以及应用于代谢工程等的相关原理、实验测量、数据分析、使用条件等,以期充分发挥代谢流量比率分析法的优势,并将其拓展推广至更多细胞体系的代谢特性阐明和代谢工程改造中去。     

Metabolic flux analysis of cultured hepatocytes exposed to plasma

Hepatic metabolism can be investigated using metabolic flux analysis (MFA), which provides a comprehensive overview of the intracellular metabolic flux distribution. The characterization of intermediary metabolism in hepatocytes is important for all biotechnological applications involving liver cells, including the development of bioartificial liver (BAL) devices. During BAL operation, hepatocytes are exposed to plasma or blood from the patient, at which time they are prone to accumulate intracellular lipids and exhibit poor liver-specific functions. In a prior study, we found that preconditioning the primary rat hepatocytes in culture medium containing physiological levels of insulin, as opposed to the typical supraphysiological levels found in standard hepatocyte culture media, reduced lipid accumulation during subsequent plasma exposure. Furthermore, supplementing the plasma with amino acids restored hepatospecific functions. In the current study, we used MFA to quantify the changes in intracellular pathway fluxes of primary rat hepatocytes in response to low-insulin preconditioning and amino acid supplementation. We found that culturing hepatocytes in medium containing lower physiological levels of insulin decreased the clearance of glucose and glycerol with a concomitant decrease in glycolysis. These findings are consistent with the general notion that low insulin, especially in the presence of high glucagon levels, downregulates glycolysis in favor of gluconeogenesis in hepatocytes. The MFA model shows that, during subsequent plasma exposure, low-insulin preconditioning upregulated gluconeogenesis, with lactate as the primary precursor in unsupplemented plasma, with a greater contribution from deaminated amino acids in amino acid-supplemented plasma. Concomitantly, low-insulin preconditioning increased fatty acid oxidation, an effect that was further enhanced by amino acid supplementation to the plasma. The increase in fatty acid oxidation reduced intracellular triglyceride accumulation. Overall, these findings are consistent with the notion that the insulin level in medium culture presets the metabolic machinery of hepatocytes such that it directly impacts on their metabolic behavior during subsequent plasma culture.     

Metabolic flux balance analysis during lactate and glucose concomitant consumption in HEK293 cell cultures

At early stages of the exponential growth phase in HEK293 cell cultures, the tricarboxylic acid cycle is unable to process all the amount of NADH generated in the glycolysis pathway, being lactate the main by-product. However, HEK293 cells are also able to metabolize lactate depending on the environmental conditions. It has been recently observed that one of the most important modes of lactate metabolization is the cometabolism of lactate and glucose, observed even during the exponential growth phase. Extracellular lactate concentration and pH appear to be the key factors triggering the metabolic shift from glucose consumption and lactate production to lactate and glucose concomitant consumption. The hypothesis proposed for triggering this metabolic shift to lactate and glucose concomitant consumption is that HEK293 cells metabolize extracellular lactate as a response to both extracellular protons and lactate accumulation, by means of cotransporting them (extracellular protons and lactate) into the cytosol. At this point, there exists a considerable controversy about how lactate reaches the mitochondrial matrix: the first hypothesis proposes that lactate is converted into pyruvate in the cytosol, and afterward, pyruvate enters into the mitochondria; the second alternative considers that lactate enters first into the mitochondria, and then, is converted into pyruvate. In this study, lactate transport and metabolization into mitochondria is shown to be feasible, as evidenced by means of respirometry tests with isolated active mitochondria, including the depletion of lactate concentration of the respirometry assay. Although the capability of lactate metabolization by isolated mitochondria is demonstrated, the possibility of lactate being converted into pyruvate in the cytosol cannot be excluded from the discussion. For this reason, the calculation of the metabolic fluxes for an HEK293 cell line was performed for the different metabolic phases observed in batch cultures under pH controlled and noncontrolled conditions, considering both hypotheses. The main objective of this study is to evaluate the redistribution of cellular metabolism and compare the differences or similarities between the phases before and after the metabolic shift of HEK293 cells (shift observed when pH is not controlled). That is from a glucose consumption/lactate production phase to a glucose-lactate coconsumption phase. Interestingly, switching to a glucose and lactate cometabolization results in a better-balanced cell metabolism, with decreased glucose and amino acids uptake rates, affecting minimally cell growth. This behavior could be applied to further develop new approaches in terms of cell engineering and to develop improved cell culture strategies in the field of animal cell technology.     

13C代谢流量分析发展30年

13C代谢流量分析(13C metabolic flux analysis,13C-MFA),是通过标记实验分析蛋白氨基酸或胞内代谢物同位素标记异构体的分布情况,从而准确定量胞内反应速率.该技术在系统理解细胞代谢特性、指导代谢工程改造和揭示病理生理学等方面起着重要作用,引起研究者的广泛重视.文中重点综述了代谢流分析30...     

Optimization-based metabolic control analysis

In this work, a novel optimization-based metabolic control analysis (OMCA) method is introduced for reducing data requirement for metabolic control analysis (MCA). It is postulated that using the optimal control approach, the fluxes in a metabolic network are correlated to metabolite concentrations and enzyme activities as a state-feedback control system that is optimal with respect to a homeostasis objective. It is then shown that the optimal feedback gains are directly related to the elasticity coefficients (ECs) of MCA. This approach requires determination of the relative "importance" of metabolites and fluxes for the system, which is possible with significantly reduced experimental data, as compared with typical MCA requirements. The OMCA approach is applied to a top-down control model of glycolysis in hepatocytes. It is statistically demonstrated that the OMCA model is capable of predicting the ECs observed experimentally with few exceptions. Further, an OMCA-based model reconciliation study shows that the modification of four assumed stoichiometric coefficients in the model can explain most of the discrepancies, with the exception of elasticities with respect to the NADH/NAD ratio.     

Control of glycolysis in rat liver by glucokinase and phosphofructokinase: influence of glucose concentration

Summary Control of glucose metabolism in rat liver under different glucose concentrations was studied. Flux Control Coefficients of glucokinase, glucose 6-phosphate isomerase and phosphofructokinase were determined by the shortening and enzyme titration method. Results obtained show that glucose concentration in liver can play an important role in control of liver glycolysis by enhancing the Flux Control Coefficient of phosphofructokinase. Possible physiological significance of this fact is discussed.     

肌苷发酵过程中后期的代谢流节点分析

应用代谢流分析方法,实验测定副产物的积累速率,将数据输入计算机,应用MArllAB软件计算肌苷发酵中后期代谢流分布。通过分析代谢网络中重要的节点,提出了优化肌苷生物合成途径的建议。     

13C-metabolic flux analysis of Clostridium ljungdahlii illuminates its core metabolism under mixotrophic culture conditions

Carbon dioxide-fixing acetogenic bacteria (acetogens) utilizing the Wood-Ljungdahl Pathway (WLP) play an important role in CO₂ fixation in the biosphere and in the development of biological processes – alone or in cocultures, under both autotrophic and mixotrophic conditions – for production of chemicals and fuels. To date, limited work has been reported in experimentally validating and quantifying reaction fluxes of their core metabolic pathways. Here, the core metabolic model of the acetogen Clostridium ljungdahlii was interrogated using ¹³C-metabolic flux analysis (¹³C-MFA), which required the development of a new defined culture medium. Autotrophic, heterotrophic, and mixotrophic growth in defined medium was possible by adding 1 mM methionine to replace yeast extract. Our ¹³C-MFA found an incomplete TCA cycle and inactive core pathways/reactions, notably those of the oxidative pentose phosphate pathway, Entner-Doudoroff pathway, and malate dehydrogenase. ¹³C-MFA during mixotrophic growth using the parallel tracers [1–¹³C]fructose, [1,2–¹³C]fructose, [1,2,3–¹³C]fructose, and [U–¹³C]asparagine found that externally supplied CO₂ contributed the majority of carbon consumed. All internally-produced CO₂ from the catabolism of asparagine and fructose was consumed by the WLP. While glycolysis of fructose was active, it was not a major contributor to overall production of ATP, NADH, and acetyl-CoA. Gluconeogenic reactions were active despite the availability of organic carbon. Asparagine was catabolized equally via conversion to threonine and subsequent cleavage to produce acetaldehyde and glycine, and via deamination to fumarate and then the anaplerotic conversion of malate to pyruvate. Both pathways for asparagine catabolism produced acetyl-CoA, either directly via pyruvate or indirectly via the WLP. Cofactor stoichiometry based on our data predicted an essentially zero flux through the ferredoxin-dependent transhydrogenase (Nfn) reaction. Instead, nearly all of NADPH generated from the hydrogenase reaction was consumed by the WLP. Reduced ferredoxin produced by the hydrogenase reaction and glycolysis was mostly used for ATP generation via the RNF/ATPase system, with the remainder consumed by the WLP. NADH produced by RNF/ATPase was entirely consumed via the WLP.     

Integrating cybernetic modeling with pathway analysis provides a dynamic, systems-level description of metabolic control

Cybernetic modeling strives to uncover the inbuilt regulatory programs of biological systems and leverage them toward computational prediction of metabolic dynamics. Because of its focus on incorporating the global aims of metabolism, cybernetic modeling provides a systems-oriented approach for describing regulatory inputs and inferring the impact of regulation within biochemical networks. Combining cybernetic control laws with concepts from metabolic pathway analysis has culminated in a systematic strategy for constructing cybernetic models, which was previously lacking. The newly devised framework relies upon the simultaneous application of local controls that maximize the net flux through each elementary flux mode and global controls that modulate the activities of these modes to optimize the overall nutritional state of the cell. The modeling concepts are illustrated using a simple linear pathway and a larger network representing anaerobic E. coli central metabolism. The E. coli model successfully describes the metabolic shift that occurs upon deleting the pta-ackA operon that is responsible for fermentative acetate production. The model also furnishes predictions that are consistent with experimental results obtained from additional knockout strains as well as strains expressing heterologous genes. Because of the stabilizing influence of the included control variables, the resulting cybernetic models are more robust and reliable than their predecessors in simulating the network response to imposed genetic and environmental perturbations.     

Metabolic flux control analysis of branch points: an improved approach to obtain flux control coefficients from large perturbation data

An overview of published approaches for the metabolic flux control analysis of branch points revealed that often not all fundamental constraints on the flux control coefficients have been taken into account. This has led to contradictory statements in literature on the minimum number of large perturbation experiments required to estimate the complete set of flux control coefficients C(J) for a metabolic branch point. An improved calculation procedure, based on approximate Lin-log reaction kinetics, is proposed, providing explicit analytical solutions of steady state fluxes and metabolite concentrations as a function of large changes in enzyme levels. The obtained solutions allow direct calculation of elasticity ratios from experimental data and subsequently all C(J)-values from the unique relation between elasticity ratio's and flux control coefficients. This procedure ensures that the obtained C(J)-values satisfy all fundamental constraints. From these it follows that for a three enzyme branch point only one characterised or two uncharacterised large flux perturbations are sufficient to obtain all C(J)- values. The improved calculation procedure is illustrated with four experimental cases.     

Distributed control of the glycolytic flux in wild-type cells and catabolite repression mutants of Saccharomyces cerevisiae growing in carbon-limited chemostat cultures

The sensitivity of the control of glycolysis was studied in the wild-type (WT) strain CEN.PK122 and in isogenic catabolite-repression mutants growing in carbon-limited, aerobic chemostat cultures at different dilution rates, D. Based on a model of glycolysis in which the glucose transport step was considered reversible and inhibited by glucose 6-phosphate (G6P), the matrix method of metabolic control analysis was applied. In the present work, we report that the control of glycolysis was significantly distributed between the glucose uptake, hexokinase, and phosphofructokinase steps. The flux control properties were sensitive to the glucose gradient through the membrane and the extent of inhibition of the transport by G6P as parameters of the glucose-uptake kinetics in all strains tested. In the WT strain at low and high D, most of the control was exerted by the phosphofructokinase (PFK)-catalyzed step. In the cat1 mutant, the step catalyzed by PFK was the most rate controlling while in the cat3 strain, the control was shared between the PFK, hexokinase (HK), and glucose transport steps. On the other hand, the mig1 mutant exhibited high control by the glucose transporter depending on the glucose gradient across the membrane. The results obtained are discussed in terms of the dependence upon the type of metabolism displayed by yeast and the kinetics of the sugar transport step.     

Metabolic flux analysis with a comprehensive isotopomer model in Bacillus subtilis

Fluxes in central carbon metabolism of a genetically engineered, riboflavin-producing Bacillus subtilis strain were investigated in glucose-limited chemostat cultures at low (0.11 h(-1)) and high (0.44 h(-1)) dilution rates. Using a mixture of 10% [U-(13)C] and 90% glucose labeled at natural abundance, (13)C-labeling experiments were carried out to provide additional information for metabolic flux balancing. The resulting labeling pattern in the proteinogenic amino acids were analyzed by two-dimensional [(13)C, (1)H] nuclear magnetic resonance (NMR) spectroscopy. To account rigorously for all available data from these experiments, we developed a comprehensive isotopomer model of B. subtilis central metabolism. Using this model, intracellular carbon net and exchange fluxes were estimated on the basis of validated physiological data and biomass composition in combination with 2D NMR data from 45 individual carbon atom spectra in the amino acids. Glucose catabolism proceeded primarily via glycolysis but pentose phosphate pathway fluxes increased with increasing growth rate. Moreover, significant back fluxes from the TCA cycle to the lower part of glycolysis via the gluconeogenic PEP carboxykinase were detected. The malic enzyme reaction, in contrast, was found to be inactive. A thorough statistical analysis was performed to prove the reliability of the isotopomer balance model and the obtained results. Specifically, a chi(2) test was applied to validate the model and the chi-square criterion was used to explore the sensitivity of model predictions to the experimental data.     

A stoichiometric model of Acidithiobacillus ferrooxidans ATCC 23270 for metabolic flux analysis

A stoichiometric model of Acidithiobacillus ferrooxidans based on the sequenced genome from strain ATCC 23270 is derived and parameterized using genome/pathway databases. The model describes the main aspects of catabolism and anabolism. By the construction and utilization of the mathematical determination of the network, metabolic flux analysis is performed for such a bacterium for the first time and results are successfully verified by comparison to literature values. This first metabolic model of A. ferrooxidans is able to simulate the main aspects of metabolism and will be useful for further investigation and improvement of bioleaching procedures. Biotechnol. Bioeng. 2009;102: 1448–1459. © 2008 Wiley Periodicals, Inc.