Quantification of proteomic and metabolic burdens predicts growth retardation and overflow metabolism in recombinant Escherichia coli

Escherichia coli has been the host organism most frequently investigated for efficient recombinant protein production. However, the production of a foreign protein in recombinant E. coli often leads to growth deterioration and elevated secretion of acetic acid. Such observed phenomena have been widely linked with cell stress responses and metabolic burdens originated particularly from the increased energy demand. In this study, flux balance analysis and dynamic flux balance analysis were applied to investigate the observed growth physiology of recombinant E. coli, incorporating the proteome allocation theory and an adjustable maintenance energy level (ATPM) to capture the proteomic and energetic burdens introduced by recombinant protein synthesis. Model predictions of biomass growth, substrate consumption, acetate excretion, and protein production with two different strains were in good agreement with the experimental data, indicating that the constraint on the available proteomic resource and the change in ATPM might be important contributors governing the growth physiology of recombinant strains. The modeling framework developed in this work, currently with several limitations to overcome, offers a starting point for the development of a practical, model-based tool to guide metabolic engineering decisions for boosting recombinant protein production.     

Understanding the Adaptive Growth Strategy of Lactobacillus plantarum by In Silico Optimisation

In the study of metabolic networks, optimization techniques are often used to predict flux distributions, and hence, metabolic phenotype. Flux balance analysis in particular has been successful in predicting metabolic phenotypes. However, an inherent limitation of a stoichiometric approach such as flux balance analysis is that it can predict only flux distributions that result in maximal yields. Hence, previous attempts to use FBA to predict metabolic fluxes in Lactobacillus plantarum failed, as this lactic acid bacterium produces lactate, even under glucose-limited chemostat conditions, where FBA predicted mixed acid fermentation as an alternative pathway leading to a higher yield. In this study we tested, however, whether long-term adaptation on an unusual and poor carbon source (for this bacterium) would select for mutants with optimal biomass yields. We have therefore adapted Lactobacillus plantarum to grow well on glycerol as its main growth substrate. After prolonged serial dilutions, the growth yield and corresponding fluxes were compared to in silico predictions. Surprisingly, the organism still produced mainly lactate, which was corroborated by FBA to indeed be optimal. To understand these results, constraint-based elementary flux mode analysis was developed that predicted 3 out of 2669 possible flux modes to be optimal under the experimental conditions. These optimal pathways corresponded very closely to the experimentally observed fluxes and explained lactate formation as the result of competition for oxygen by the other flux modes. Hence, these results provide thorough understanding of adaptive evolution, allowing in silico predictions of the resulting flux states, provided that the selective growth conditions favor yield optimization as the winning strategy.     

Inclusion of maintenance energy improves the intracellular flux predictions of CHO

Chinese hamster ovary (CHO) cells are the leading platform for the production of biopharmaceuticals with human-like glycosylation. The standard practice for cell line generation relies on trial and error approaches such as adaptive evolution and high-throughput screening, which typically take several months. Metabolic modeling could aid in designing better producer cell lines and thus shorten development times. The genome-scale metabolic model (GSMM) of CHO can accurately predict growth rates. However, in order to predict rational engineering strategies it also needs to accurately predict intracellular fluxes. In this work we evaluated the agreement between the fluxes predicted by parsimonious flux balance analysis (pFBA) using the CHO GSMM and a wide range of ¹³C metabolic flux data from literature. While glycolytic fluxes were predicted relatively well, the fluxes of tricarboxylic acid (TCA) cycle were vastly underestimated due to too low energy demand. Inclusion of computationally estimated maintenance energy significantly improved the overall accuracy of intracellular flux predictions. Maintenance energy was therefore determined experimentally by running continuous cultures at different growth rates and evaluating their respective energy consumption. The experimentally and computationally determined maintenance energy were in good agreement. Additionally, we compared alternative objective functions (minimization of uptake rates of seven nonessential metabolites) to the biomass objective. While the predictions of the uptake rates were quite inaccurate for most objectives, the predictions of the intracellular fluxes were comparable to the biomass objective function.     

Simulating Serial-Target Antibacterial Drug Synergies Using Flux Balance Analysis

Flux balance analysis (FBA) is an increasingly useful approach for modeling the behavior of metabolic systems. However, standard FBA modeling of genetic knockouts cannot predict drug combination synergies observed between serial metabolic targets, even though such synergies give rise to some of the most widely used antibiotic treatments. Here we extend FBA modeling to simulate responses to chemical inhibitors at varying concentrations, by diverting enzymatic flux to a waste reaction. This flux diversion yields very similar qualitative predictions to prior methods for single target activity. However, we find very different predictions for combinations, where flux diversion, which mimics the kinetics of competitive metabolic inhibitors, can explain serial target synergies between metabolic enzyme inhibitors that we confirmed in Escherichia coli cultures. FBA flux diversion opens the possibility for more accurate genome-scale predictions of drug synergies, which can be used to suggest treatments for infections and other diseases.     

L-缬氨酸生物合成中的代谢流量分析

应用流量平衡模型 ,通过物料衡算和MATLAB线性规划方法得到了发酵中后期L 缬氨酸合成过程的代谢流量分步。代谢流分析结果表明 ,在分批培养生成L 缬氨酸的过程中 ,有62 8%的葡萄糖进入糖酵解途径生成L 缬氨酸 ,38 2 %进入HMP途径 ,仅 9 2 %的碳架进入TCA循环。实验条件下的代谢流 (58)与理想代谢流 (92 31 )相比 ,仍应从遗传改造和发酵控制方面降低TCA循环的代谢流 ,减少副产氨基酸的生成来进一步提高缬氨酸的产率。     

Flux balance analysis in the production of clavulanic acid by Streptomyces clavuligerus

In this work, in silico flux balance analysis is used for predicting the metabolic behavior of Streptomyces clavuligerus during clavulanic acid production. To choose the best objective function for use in the analysis, three different optimization problems are evaluated inside the flux balance analysis formulation: (i) maximization of the specific growth rate, (ii) maximization of the ATP yield, and (iii) maximization of clavulanic acid production. Maximization of ATP yield showed the best predictions for the cellular behavior. Therefore, flux balance analysis using ATP as objective function was used for analyzing different scenarios of nutrient limitations toward establishing the effect of limiting the carbon, nitrogen, phosphorous, and oxygen sources on the growth and clavulanic acid production rates. Obtained results showed that ammonia and phosphate limitations are the ones most strongly affecting clavulanic acid biosynthesis. Furthermore, it was possible to identify the ornithine flux from the urea cycle and the α‐ketoglutarate flux from the TCA cycle as the most determinant internal fluxes for promoting clavulanic acid production. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:1226–1236, 2015     

The Ability of Flux Balance Analysis to Predict Evolution of Central Metabolism Scales with the Initial Distance to the Optimum

The most powerful genome-scale framework to model metabolism, flux balance analysis (FBA), is an evolutionary optimality model. It hypothesizes selection upon a proposed optimality criterion in order to predict the set of internal fluxes that would maximize fitness. Here we present a direct test of the optimality assumption underlying FBA by comparing the central metabolic fluxes predicted by multiple criteria to changes measurable by a ¹³C-labeling method for experimentally-evolved strains. We considered datasets for three Escherichia coli evolution experiments that varied in their length, consistency of environment, and initial optimality. For ten populations that were evolved for 50,000 generations in glucose minimal medium, we observed modest changes in relative fluxes that led to small, but significant decreases in optimality and increased the distance to the predicted optimal flux distribution. In contrast, seven populations evolved on the poor substrate lactate for 900 generations collectively became more optimal and had flux distributions that moved toward predictions. For three pairs of central metabolic knockouts evolved on glucose for 600–800 generations, there was a balance between cases where optimality and flux patterns moved toward or away from FBA predictions. Despite this variation in predictability of changes in central metabolism, two generalities emerged. First, improved growth largely derived from evolved increases in the rate of substrate use. Second, FBA predictions bore out well for the two experiments initiated with ancestors with relatively sub-optimal yield, whereas those begun already quite optimal tended to move somewhat away from predictions. These findings suggest that the tradeoff between rate and yield is surprisingly modest. The observed positive correlation between rate and yield when adaptation initiated further from the optimum resulted in the ability of FBA to use stoichiometric constraints to predict the evolution of metabolism despite selection for rate.     

Stoichiometric flux balance models quantitatively predict growth and metabolic by-product secretion in wild-type Escherichia coli W3110.

Flux balance models of metabolism use stoichiometry of metabolic pathways, metabolic demands of growth, and optimality principles to predict metabolic flux distribution and cellular growth under specified environmental conditions. These models have provided a mechanistic interpretation of systemic metabolic physiology, and they are also useful as a quantitative tool for metabolic pathway design. Quantitative predictions of cell growth and metabolic by-product secretion that are experimentally testable can be obtained from these models. In the present report, we used independent measurements to determine the model parameters for the wild-type Escherichia coli strain W3110. We experimentally determined the maximum oxygen utilization rate (15 mmol of O2 per g [dry weight] per h), the maximum aerobic glucose utilization rate (10.5 mmol of Glc per g [dry weight] per h), the maximum anaerobic glucose utilization rate (18.5 mmol of Glc per g [dry weight] per h), the non-growth-associated maintenance requirements (7.6 mmol of ATP per g [dry weight] per h), and the growth-associated maintenance requirements (13 mmol of ATP per g of biomass). The flux balance model specified by these parameters was found to quantitatively predict glucose and oxygen uptake rates as well as acetate secretion rates observed in chemostat experiments. We have formulated a predictive algorithm in order to apply the flux balance model to describe unsteady-state growth and by-product secretion in aerobic batch, fed-batch, and anaerobic batch cultures. In aerobic experiments we observed acetate secretion, accumulation in the culture medium, and reutilization from the culture medium. In fed-batch cultures acetate is cometabolized with glucose during the later part of the culture period.(ABSTRACT TRUNCATED AT 250 WORDS)     

Dynamic flux balance analysis of the metabolism of Saccharomyces cerevisiae during the shift from fully respirative or respirofermentative metabolic states to anaerobiosis

Dynamic flux balance analysis was utilized to simulate the metabolic behaviour of initially fully respirative and respirofermentative steady-state cultures of Saccharomyces?cerevisiae during sudden oxygen depletion. The hybrid model for the dynamic flux balance analysis included a stoichiometric genome-scale metabolic model as a static part and dynamic equations for the uptake of glucose and the cessation of respirative metabolism. The yeast consensus genome-scale metabolic model [Herrg?rd MJ et?al. (2008) Nat Biotechnol26, 1155-1160; Dobson PD et?al. (2010) BMC Syst Biol4, 145] was refined with respect to oxygen-dependent energy metabolism and further modified to reflect S.?cerevisiae anabolism in the absence of oxygen. Dynamic flux balance analysis captured well the essential features of the dynamic metabolic behaviour of S.?cerevisiae during adaptation to anaerobiosis. Modelling and simulation enabled the identification of short time-scale flux distribution dynamics under the transition to anaerobic metabolism, during which the specific growth rate was reduced, as well as longer time-scale process dynamics when the specific growth rate recovered. Expression of the metabolic genes was set into the context of the identified dynamics. Metabolic gene expression responses associated with the specific growth rate and with the cessation of respirative metabolism were distinguished.     

Genome-scale reconstruction and analysis of the Pseudomonas putida KT2440 metabolic network facilitates applications in biotechnology

A cornerstone of biotechnology is the use of microorganisms for the efficient production of chemicals and the elimination of harmful waste. Pseudomonas putida is an archetype of such microbes due to its metabolic versatility, stress resistance, amenability to genetic modifications, and vast potential for environmental and industrial applications. To address both the elucidation of the metabolic wiring in P. putida and its uses in biocatalysis, in particular for the production of non-growth-related biochemicals, we developed and present here a genome-scale constraint-based model of the metabolism of P. putida KT2440. Network reconstruction and flux balance analysis (FBA) enabled definition of the structure of the metabolic network, identification of knowledge gaps, and pin-pointing of essential metabolic functions, facilitating thereby the refinement of gene annotations. FBA and flux variability analysis were used to analyze the properties, potential, and limits of the model. These analyses allowed identification, under various conditions, of key features of metabolism such as growth yield, resource distribution, network robustness, and gene essentiality. The model was validated with data from continuous cell cultures, high-throughput phenotyping data, (13)C-measurement of internal flux distributions, and specifically generated knock-out mutants. Auxotrophy was correctly predicted in 75% of the cases. These systematic analyses revealed that the metabolic network structure is the main factor determining the accuracy of predictions, whereas biomass composition has negligible influence. Finally, we drew on the model to devise metabolic engineering strategies to improve production of polyhydroxyalkanoates, a class of biotechnologically useful compounds whose synthesis is not coupled to cell survival. The solidly validated model yields valuable insights into genotype-phenotype relationships and provides a sound framework to explore this versatile bacterium and to capitalize on its vast biotechnological potential.     

A possibilistic framework for constraint-based metabolic flux analysis

Background  

Constraint-based models allow the calculation of the metabolic flux states that can be exhibited by cells, standing out as a powerful analytical tool, but they do not determine which of these are likely to be existing under given circumstances. Typical methods to perform these predictions are (a) flux balance analysis, which is based on the assumption that cell behaviour is optimal, and (b) metabolic flux analysis, which combines the model with experimental measurements.     

柠檬酸钠对L-组氨酸发酵代谢流分布的影响

目的：建立谷氨酸棒杆菌TL1105生物合成L-组氨酸的代谢网络模型,并进行代谢网络计量分析。方法：通过所构建的L-组氨酸代谢网络模型,利用MATLAB软件计算出添加柠檬酸钠和不添加柠檬酸钠发酵中后期代谢网络的代谢流分布。结果：在L-组氨酸分批发酵过程中,在发酵初期未添加柠檬酸钠的条件下流向戊糖磷酸途径（HMP）的代谢流为9.59,合成组氨酸的代谢流为8.91;在发酵初期添加2g/L柠檬酸钠的条件下流向HMP的代谢流为12.74,合成组氨酸的代谢流为9.61。结论：在发酵初期添加柠檬酸钠能够改变L-组氨酸生物合成途径的关键节点6-磷酸葡萄糖、丙酮酸及乙酰辅酶A的代谢流分布,保持糖酵解途径、三羧酸循环与HMP之间代谢流量平衡,有利于提高L-组氨酸生物合成途径的代谢流量,最终使流向组氨酸的代谢流增加了7.86%。     

Selection of objective function in genome scale flux balance analysis for process feed development in antibiotic production

Using flux variability analysis of a genome scale metabolic network of Streptomyces coelicolor, a series of reactions were identified, from disparate pathways that could be combined into an actinorhodin-generating mini-network. Candidate process feed nutrients that might be expected to influence this network were used in process simulations and in silico predictions compared to experimental findings. Ranking potential process feeds by flux balance analysis optimisation, using either growth or antibiotic production as objective function, did not correlate with experimental actinorhodin yields in fed processes. However, the effect of the feeds on glucose assimilation rate (using glucose uptake as objective function) ranked them in the same order as in vivo antibiotic production efficiency, consistent with results of a robustness analysis of the effect of glucose assimilation on actinorhodin production.     

Flux Imbalance Analysis and the Sensitivity of Cellular Growth to Changes in Metabolite Pools

Stoichiometric models of metabolism, such as flux balance analysis (FBA), are classically applied to predicting steady state rates - or fluxes - of metabolic reactions in genome-scale metabolic networks. Here we revisit the central assumption of FBA, i.e. that intracellular metabolites are at steady state, and show that deviations from flux balance (i.e. flux imbalances) are informative of some features of in vivo metabolite concentrations. Mathematically, the sensitivity of FBA to these flux imbalances is captured by a native feature of linear optimization, the dual problem, and its corresponding variables, known as shadow prices. First, using recently published data on chemostat growth of Saccharomyces cerevisae under different nutrient limitations, we show that shadow prices anticorrelate with experimentally measured degrees of growth limitation of intracellular metabolites. We next hypothesize that metabolites which are limiting for growth (and thus have very negative shadow price) cannot vary dramatically in an uncontrolled way, and must respond rapidly to perturbations. Using a collection of published datasets monitoring the time-dependent metabolomic response of Escherichia coli to carbon and nitrogen perturbations, we test this hypothesis and find that metabolites with negative shadow price indeed show lower temporal variation following a perturbation than metabolites with zero shadow price. Finally, we illustrate the broader applicability of flux imbalance analysis to other constraint-based methods. In particular, we explore the biological significance of shadow prices in a constraint-based method for integrating gene expression data with a stoichiometric model. In this case, shadow prices point to metabolites that should rise or drop in concentration in order to increase consistency between flux predictions and gene expression data. In general, these results suggest that the sensitivity of metabolic optima to violations of the steady state constraints carries biologically significant information on the processes that control intracellular metabolites in the cell.     

Cellular Growth Arrest and Persistence from Enzyme Saturation

Metabolic efficiency depends on the balance between supply and demand of metabolites, which is sensitive to environmental and physiological fluctuations, or noise, causing shortages or surpluses in the metabolic pipeline. How cells can reliably optimize biomass production in the presence of metabolic fluctuations is a fundamental question that has not been fully answered. Here we use mathematical models to predict that enzyme saturation creates distinct regimes of cellular growth, including a phase of growth arrest resulting from toxicity of the metabolic process. Noise can drive entry of single cells into growth arrest while a fast-growing majority sustains the population. We confirmed these predictions by measuring the growth dynamics of Escherichia coli utilizing lactose as a sole carbon source. The predicted heterogeneous growth emerged at high lactose concentrations, and was associated with cell death and production of antibiotic-tolerant persister cells. These results suggest how metabolic networks may balance costs and benefits, with important implications for drug tolerance.     

L-色氨酸生物合成的代谢流量分析

建立了谷氨酸棒杆菌合成L-色氨酸(L-Try)的代谢流量平衡模型,应用该模型计算出发酵中后期的代谢流分布并通过MATLAB软件线性规划得到Try理想代谢流分布。结果表明75．15％的碳架进入糖酵解,24．85％的碳架进入HMP途径;但与理想代谢流相比,应从遗传改造和发酵控制方面降低TCA循环的代谢流,减少副产氨基酸的生成,摸索最适的溶氧控制对提高Try产率至关重要。     

辅因子在微生物细胞工厂中的代谢调控与应用*

生物体中大部分酶催化反应都需要辅因子参与,辅因子平衡对维持正常的细胞代谢至关重要,而辅因子失衡则会导致细胞生长和生产的紊乱。在微生物细胞工厂的构建中,通过调节辅因子代谢平衡来提高产物合成途径的效率,从而调控细胞生长与产物生产,使代谢流能够最大限度地流向目标产物,已经成为代谢调控的重要手段。目前常见的用于代谢调控的辅因子有NAD(P)H/NAD(P)⁺、辅酶、ATP/ADP等。围绕这几种辅因子的代谢途径及功能分类进行了综述,并总结了微生物中不同产物利用辅因子平衡策略进行合成调控的研究,以期为各类化合物的高效生物合成提供参考。     

Prediction of dynamic behavior of mutant strains from limited wild-type data

Metabolic engineering is the field of introducing genetic changes in organisms so as to modify their function towards synthesizing new products of high impact to society. However, engineered cells frequently have impaired growth rates thus seriously limiting the rate at which such products are made. The problem is attributable to inadequate understanding of how a metabolic network functions in a dynamic sense. Predictions of mutant strain behavior in the past have been based on steady state theories such as flux balance analysis (FBA), minimization of metabolic adjustment (MOMA), and regulatory on/off minimization (ROOM). Such predictions are restricted to product yields and cannot address productivity, which is of focal interest to applications. We demonstrate that our framework ( [Song and Ramkrishna, 2010] and [Song and Ramkrishna, 2011]), based on a “cybernetic” view of metabolic systems, makes predictions of the dynamic behavior of mutant strains of Escherichia coli from a limited amount of data obtained from the wild-type. Dynamic frameworks must necessarily address the issue of metabolic regulation, which the cybernetic approach does by postulating that metabolism is an optimal dynamic response of the organism to the environment in driving reactions towards ensuring survival. The predictions made in this paper are without parallel in the literature and lay the foundation for rational metabolic engineering.     

Predicting Growth Conditions from Internal Metabolic Fluxes in an In-Silico Model of E. coli

A widely studied problem in systems biology is to predict bacterial phenotype from growth conditions, using mechanistic models such as flux balance analysis (FBA). However, the inverse prediction of growth conditions from phenotype is rarely considered. Here we develop a computational framework to carry out this inverse prediction on a computational model of bacterial metabolism. We use FBA to calculate bacterial phenotypes from growth conditions in E. coli, and then we assess how accurately we can predict the original growth conditions from the phenotypes. Prediction is carried out via regularized multinomial regression. Our analysis provides several important physiological and statistical insights. First, we show that by analyzing metabolic end products we can consistently predict growth conditions. Second, prediction is reliable even in the presence of small amounts of impurities. Third, flux through a relatively small number of reactions per growth source (∼10) is sufficient for accurate prediction. Fourth, combining the predictions from two separate models, one trained only on carbon sources and one only on nitrogen sources, performs better than models trained to perform joint prediction. Finally, that separate predictions perform better than a more sophisticated joint prediction scheme suggests that carbon and nitrogen utilization pathways, despite jointly affecting cellular growth, may be fairly decoupled in terms of their dependence on specific assortments of molecular precursors.