On-line optimization of glutamate production based on balanced metabolic control by RQ

In glutamate fermentations by Corynebacterium glutamicum, higher glutamate concentration could be achieved by constantly controlling dissolved oxygen concentration (DO) at a lower level; however, by-product lactate also severely accumulated. The results of analyzing activities changes of the two key enzymes, glutamate and lactate dehydrogenases involved with the fermentation, and the entire metabolic network flux analysis showed that the lactate overproduction was because the metabolic flux in TCA cycle was too low to balance the glucose glycolysis rate. As a result, the respiratory quotient (RQ) adaptive control based “balanced metabolic control” (BMC) strategy was proposed and used to regulate the TCA metabolic flux rate at an appropriate level to achieve the metabolic balance among glycolysis, glutamate synthesis, and TCA metabolic flux. Compared with the best results of various DO constant controls, the BMC strategy increased the maximal glutamate concentration by about 15% and almost completely repressed the lactate accumulation with competitively high glutamate productivity.     

Studies of regulatory metabolism in Moniezia expansa: general considerations.

The activities of selected enzymes in the branched metabolic pathway to succinate or lactate were determined in cytosol and mitochondrial fractions. The enzymes of lowest activity in the cytosol, and thus possibly regulatory, are phosphofructokinase and pyruvate kinase. Malic enzyme activity could scarcely be detected in either compartment; phosphoenolpyruvate carboxykinase and malate dehydrogenase occur in both. The end products of metabolism are succinate and lactate; under anaerobic conditions lactate production increases whereas succinate production shows a small decrease. The presence of glucose in the medium does not influence the change, but causes an increase in total endproduct accumulation. Levels of metabolic intermediates in worms incubated aerobically and anaerobically are presented, and ‘cross-over’ plots and calculations of apparent equilibrium constants identify hexokinase, phosphofructokinase and pyruvate kinase as regulatory. Under aerobic conditions a large increase in the size of the malate pool is observed suggesting that the depression of lactate production is produced by its inhibitory effect on pyruvate kinase. Adenine nucleotide levels are maintained whether or not the worm is incubated under anaerobic conditions.     

Incorporating expression data in metabolic modeling: a case study of lactate dehydrogenase

Integrating biological information from different sources to understand cellular processes is an important problem in systems biology. We use data from mRNA expression arrays and chemical kinetics to formulate a metabolic model relevant to K562 erythroleukemia cells. MAP kinase pathway activation alters the expression of metabolic enzymes in K562 cells. Our array data show changes in expression of lactate dehydrogenase (LDH) isoforms after treatment with phorbol 12-myristate 13-acetate (PMA), which activates MAP kinase signaling. We model the change in lactate production which occurs when the MAP kinase pathway is activated, using a non-equilibrium, chemical-kinetic model of homolactic fermentation. In particular, we examine the role of LDH isoforms, which catalyse the conversion of pyruvate to lactate. Changes in the isoform ratio are not the primary determinant of the production of lactate. Rather, the total concentration of LDH controls the lactate concentration.     

Elevated pCO2 affects the lactate metabolic shift in CHO cell culture processes

The shift from lactate production to consumption in CHO cell metabolism is a key event during cell culture cultivations and is connected to increased culture longevity and final product titers. However, the mechanisms controlling this metabolic shift are not yet fully understood. Variations in lactate metabolism have been mainly reported to be induced by process pH and availability of substrates like glucose and glutamine. The aim of this study was to investigate the effects of elevated pCO₂ concentrations on the lactate metabolic shift phenomena in CHO cell culture processes. In this publication, we show that at elevated pCO₂ in batch and fed‐batch cultures, the lactate metabolic shift was absent in comparison to control cultures at lower pCO₂ values. Furthermore, through metabolic flux analysis we found a link between the lactate metabolic shift and the ratio of NADH producing and regenerating intracellular pathways. This ratio was mainly affected by a reduced oxidative capacity of cultures at elevated pCO₂. The presented results are especially interesting for large‐scale and perfusion processes where increased pCO₂ concentrations are likely to occur. Our results suggest, that so far unexplained metabolic changes may be connected to increased pCO₂ accumulation in larger scale fermentations. Finally, we propose several mechanisms through which increased pCO₂ might affect the cell metabolism and briefly discuss methods to enable the lactate metabolic shift during cell cultivations.     

Multiplicity of Steady States in Glycolysis and Shift of Metabolic State in Cultured Mammalian Cells

Cultured mammalian cells exhibit elevated glycolysis flux and high lactate production. In the industrial bioprocesses for biotherapeutic protein production, glucose is supplemented to the culture medium to sustain continued cell growth resulting in the accumulation of lactate to high levels. In such fed-batch cultures, sometimes a metabolic shift from a state of high glycolysis flux and high lactate production to a state of low glycolysis flux and low lactate production or even lactate consumption is observed. While in other cases with very similar culture conditions, the same cell line and medium, cells continue to produce lactate. A metabolic shift to lactate consumption has been correlated to the productivity of the process. Cultures that exhibited the metabolic shift to lactate consumption had higher titers than those which didn’t. However, the cues that trigger the metabolic shift to lactate consumption state (or low lactate production state) are yet to be identified. Metabolic control of cells is tightly linked to growth control through signaling pathways such as the AKT pathway. We have previously shown that the glycolysis of proliferating cells can exhibit bistability with well-segregated high flux and low flux states. Low lactate production (or lactate consumption) is possible only at a low glycolysis flux state. In this study, we use mathematical modeling to demonstrate that lactate inhibition together with AKT regulation on glycolysis enzymes can profoundly influence the bistable behavior, resulting in a complex steady-state topology. The transition from the high flux state to the low flux state can only occur in certain regions of the steady state topology, and therefore the metabolic fate of the cells depends on their metabolic trajectory encountering the region that allows such a metabolic state switch. Insights from such switch behavior present us with new means to control the metabolism of mammalian cells in fed-batch cultures.     

Kinetic modelling of lactate utilization and butyrate production by key human colonic bacterial species

Butyrate is the preferred energy source for colonocytes and has an important role in gut health; in contrast, accumulation of high concentrations of lactate is detrimental to gut health. The major butyrate-producing bacterial species in the human colon belong to the Firmicutes. Eubacterium hallii and a new species, Anaerostipes coli SS2/1, members of clostridial cluster XIVa, are able to utilize lactate and acetate via the butyryl CoA : acetate CoA transferase route, the main metabolic pathway for butyrate synthesis in the human colon. Here we provide a mathematical model to analyse the production of butyrate by lactate-utilizing bacteria from the human colon. The model is an aggregated representation of the fermentation pathway. The parameters of the model were estimated using total least squares and maximum likelihood, based on in vitro experimental data with E. hallii L2-7 and A. coli SS2/1. The findings of the mathematical model adequately match those from the bacterial batch culture experiments. Such an in silico approach should provide insight into carbohydrate fermentation and short-chain fatty acid cross-feeding by dominant species of the human colonic microbiota.     

Comparative genome-scale metabolic modeling of actinomycetes: the topology of essential core metabolism

Actinomycetes are highly important bacteria. On one hand, some of them cause severe human and plant diseases, on the other hand, many species are known for their ability to produce antibiotics. Here we report the results of a comparative analysis of genome-scale metabolic models of 37 species of actinomycetes. Based on in silico knockouts we generated topological and genomic maps for each organism. Combining the collection of genome-wide models, we constructed a global enzyme association network to identify both a conserved "core network" and an "essential core network" of the entire group. As has been reported for low-degree metabolites in several organisms, low-degree enzymes (in linear pathways) turn out to be generally more essential than high-degree enzymes (in metabolic hubs).     

Metabolic Engineering of Escherichia coli for Enhanced Production of Succinic Acid, Based on Genome Comparison and In Silico Gene Knockout Simulation

Comparative analysis of the genomes of mixed-acid-fermenting Escherichia coli and succinic acid-overproducing Mannheimia succiniciproducens was carried out to identify candidate genes to be manipulated for overproducing succinic acid in E. coli. This resulted in the identification of five genes or operons, including ptsG, pykF, sdhA, mqo, and aceBA, which may drive metabolic fluxes away from succinic acid formation in the central metabolic pathway of E. coli. However, combinatorial disruption of these rationally selected genes did not allow enhanced succinic acid production in E. coli. Therefore, in silico metabolic analysis based on linear programming was carried out to evaluate the correlation between the maximum biomass and succinic acid production for various combinatorial knockout strains. This in silico analysis predicted that disrupting the genes for three pyruvate forming enzymes, ptsG, pykF, and pykA, allows enhanced succinic acid production. Indeed, this triple mutation increased the succinic acid production by more than sevenfold and the ratio of succinic acid to fermentation products by ninefold. It could be concluded that reducing the metabolic flux to pyruvate is crucial to achieve efficient succinic acid production in E. coli. These results suggest that the comparative genome analysis combined with in silico metabolic analysis can be an efficient way of developing strategies for strain improvement.     

Metabolic flux analysis for calcium dependent antibiotic (CDA) production in Streptomyces coelicolor

The calcium dependent antibiotic (CDA) is a nonribosomal lipopeptide produced by Streptomyces coelicolor. We constructed a metabolic network of more than 400 reactions for the primary and secondary metabolism of S. coelicolor and used computational metabolic flux balancing to investigate some of the factors affecting growth and production of CDA. Computational results indicated that the CDA production was concomitant with growth. Computational specific growth rates were twice as high as the experimental specific growth rates. Metabolic flux distributions and sensitivity analyses computed for various phases of the batch culture indicated that the specific CDA production rate was affected by nitrogen assimilation, pentose phosphate pathway, shikimate biosynthesis, and oxoglutarate fluxes. Consequently, these metabolic targets were tested using genetic deletions in the model which increased the in silico specific CDA production rate.     

Identification of metabolic engineering targets for the enhancement of 1,4-butanediol production in recombinant E. coli using large-scale kinetic models

Rational metabolic engineering methods are increasingly employed in designing the commercially viable processes for the production of chemicals relevant to pharmaceutical, biotechnology, and food and beverage industries. With the growing availability of omics data and of methodologies capable to integrate the available data into models, mathematical modeling and computational analysis are becoming important in designing recombinant cellular organisms and optimizing cell performance with respect to desired criteria. In this contribution, we used the computational framework ORACLE (Optimization and Risk Analysis of Complex Living Entities) to analyze the physiology of recombinant Escherichia coli producing 1,4-butanediol (BDO) and to identify potential strategies for improved production of BDO. The framework allowed us to integrate data across multiple levels and to construct a population of large-scale kinetic models despite the lack of available information about kinetic properties of every enzyme in the metabolic pathways. We analyzed these models and we found that the enzymes that primarily control the fluxes leading to BDO production are part of central glycolysis, the lower branch of tricarboxylic acid (TCA) cycle and the novel BDO production route. Interestingly, among the enzymes between the glucose uptake and the BDO pathway, the enzymes belonging to the lower branch of TCA cycle have been identified as the most important for improving BDO production and yield. We also quantified the effects of changes of the target enzymes on other intracellular states like energy charge, cofactor levels, redox state, cellular growth, and byproduct formation. Independent earlier experiments on this strain confirmed that the computationally obtained conclusions are consistent with the experimentally tested designs, and the findings of the present studies can provide guidance for future work on strain improvement. Overall, these studies demonstrate the potential and effectiveness of ORACLE for the accelerated design of microbial cell factories.     

The impact of anti-apoptotic gene Bcl-2∆ expression on CHO central metabolism

Anti-apoptosis engineering is an established technique to prolong the viability of mammalian cell cultures used for industrial production of recombinant proteins. However, the effect of overexpressing anti-apoptotic proteins on central carbon metabolism has not been systematically studied. We transfected CHO-S cells to express Bcl-2∆, an engineered anti-apoptotic gene, and selected clones that differed in their Bcl-2∆ expression and caspase activity. ¹³C metabolic flux analysis (MFA) was then applied to elucidate the metabolic alterations induced by Bcl-2∆. Expression of Bcl-2Δ reduced lactate accumulation by redirecting the fate of intracellular pyruvate toward mitochondrial oxidation during the lactate-producing phase, and it significantly increased lactate re-uptake during the lactate-consuming phase. This flux redistribution was associated with significant increases in biomass yield, peak viable cell density (VCD), and integrated VCD. Additionally, Bcl-2∆ expression was associated with significant increases in isocitrate dehydrogenase and NADH oxidase activities, both rate-controlling mitochondrial enzymes. This is the first comprehensive ¹³C MFA study to demonstrate that expression of anti-apoptotic genes has a significant impact on intracellular metabolic fluxes, especially in controlling the fate of pyruvate carbon, which has important biotechnology applications for reducing lactate accumulation and enhancing productivity in mammalian cell cultures.     

Metabolic flux network and analysis of fermentative hydrogen production

Fermentative hydrogen production (FHP) has received a great R & D interest in recent decades, as it offers a potential means of producing H₂ from a variety of renewable resources, even wastewater via a low energy continuous process. Various extracellular metabolites including ethanol, acetate, butyrate and lactate can be produced during the fermentation, building a complex metabolic network of the FHP. Except for the recognition of its complexity, the metabolic flux network has not been well understood. Studies on biochemical reactions and metabolic flux network associated with the FHP in anaerobic fermentation system have only been drawn attention in recent years. This review summarizes the biochemical reactions taking place in the metabolic network of FHP. We discuss how the key operation factors influence metabolism in the FHP process. Recently developed and applied technologies for metabolic flux analysis have been described. Future studies on the metabolic network to enhance fermentative hydrogen production by strict anaerobes are recommended. It is expected that this review can provide useful information in terms of fundamental knowledge and update technology for scientists and research engineers in the field of biological hydrogen production.     

Metabolic profiling by 13C-NMR spectroscopy: [1,2-13C2]glucose reveals a heterogeneous metabolism in human leukemia T cells

Metabolic profiling is defined as the simultaneous assessment of substrate fluxes within and among the different pathways of metabolite synthesis and energy production under various physiological conditions. The use of stable-isotope tracers and the analysis of the distribution of labeled carbons in various intermediates, by both mass spectrometry and NMR spectroscopy, allow the role of several metabolic processes in cell growth and death to be defined. In the present paper we describe the metabolic profiling of Jurkat cells by isotopomer analysis using (13)C-NMR spectroscopy and [1,2-(13)C(2)]glucose as the stable-isotope tracer. The isotopomer analysis of the lactate, alanine, glutamate, proline, serine, glycine, malate and ribose-5-phosphate moiety of nucleotides has allowed original integrated information regarding the pentose phosphate pathway, TCA cycle, and amino acid metabolism in proliferating human leukemia T cells to be obtained. In particular, the contribution of the glucose-6-phosphate dehydrogenase and transketolase activities to phosphoribosyl-pyrophosphate synthesis was evaluated directly by the determination of isotopomers of the [1'-(13)C], [4',5'-(13)C(2)]ribosyl moiety of nucleotides. Furthermore, the relative contribution of the glycolysis and pentose cycle to lactate production was estimated via analysis of lactate isotopomers. Interestingly, pyruvate carboxylase and pyruvate dehydrogenase flux ratios measured by glutamate isotopomers and the production of isotopomers of several metabolites showed that the metabolic processes described could not take place simultaneously in the same macrocompartments (cells). Results revealed a heterogeneous metabolism in an asynchronous cell population that may be interpreted on the basis of different metabolic phenotypes of subpopulations in relation to different cell cycle phases.     

Bioinformatics and metabolic flux analysis highlight a new mechanism involved in lactate oxidation in Clostridium tyrobutyricum

Climate change and environmental issues compel us to find alternatives to the production of molecules of interest from petrochemistry. This study aims at understanding the production of butyrate, hydrogen, and CO₂ from the oxidation of lactate with acetate in Clostridium tyrobutyricum and thus proposes an alternative carbon source to glucose. This specie is known to produce more butyrate than the other butyrate-producing clostridia species due to a lack of solvent genesis phase. The recent discoveries on flavin-based electron bifurcation and confurcation mechanism as a mode of energy conservation led us to suggest a new metabolic scheme for the formation of butyrate from lactate-acetate co-metabolism. While searching for genes encoding for EtfAB complexes and neighboring genes in the genome of C. tyrobutyricum, we identified a cluster of genes involved in butyrate formation and another cluster involved in lactate oxidation homologous to Acetobacterium woodii. A phylogenetic approach encompassing other butyrate-producing and/or lactate-oxidizing species based on EtfAB complexes confirmed these results. A metabolic scheme on the production of butyrate, hydrogen, and CO₂ from the lactate-acetate co-metabolism in C. tyrobutyricum was constructed and then confirmed with data of steady-state continuous culture. This in silico metabolic carbon flux analysis model showed the coherence of the scheme from the carbon recovery, the cofactor ratio, and the ATP yield. This study improves our understanding of the lactate oxidation metabolic pathways and the role of acetate and intracellular redox balance, and paves the way for the production of molecules of interest as butyrate and hydrogen with C. tyrobutyricum.

     

Uncoupling of grwoth and acid production in Streptococcus cremoris

In lactose and leucine-limited continuous cultures of Streptococcus cremoris a linear relationship exists between specific rate of lactate production and specific growth rate. The rate of acid production in leucine-limited cultures is much higher than in lactose-limited cultures, indicating that under these conditions metabolic energy production is not coupled to growth and that metabolic energy has to be dissipated S. cremoris contains phosphofructokinase and fructose-1,6-diphosphatase, joint action of these two enzymes results in an ATP consuming futile cycle. Analyses of intracellular metabolite pools suggested that AMP and phosphoenolpyruvate play important roles in the regulation of the activity of this futile cycle.Abbreviations PEP Phosphoenolpyruvate - PFK phosphofructokinase (EC 2.7.1.11) - FBPase fructose-1,6-bisphosphatase (EC 3.1.3.11)     

Experimental determination of control of glycolysis in Lactococcus lactis

The understanding of control of metabolic processes requires quantitative studies of the importance of the different enzymatic steps for the magnitude of metabolic fluxes and metabolite concentrations. An important element in such studies is the modulation of enzyme activities in small steps above and below the wild-type level. We review a genetic approach that is well suited for both Metabolic Optimization and Metabolic Control Analysis and studies on the importance of a number of glycolytic enzymes for metabolic fluxes in Lactococcus lactis. The glycolytic enzymes phosphofructokinase (PFK), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), pyruvate kinase (PYK) and lactate dehydrogenase (LDH) are shown to have no significant control on the glycolytic flux in exponentially growing cells of L. lactis MG1363. Introduction of an uncoupled ATPase activity results in uncoupling of glycolysis from biomass production. With MG1363 growing in defined medium supplemented with glucose, the ATP demanding processes do not have a significant control on the glycolytic flux; it appears that glycolysis is running at maximal rate. It is likely that the flux control is distributed over many enzymes in L. lactis, but it cannot yet be excluded that one of the remaining glycolytic steps is a rate-limiting step for the glycolytic flux.     

Ensemble modeling for strain development of l-lysine-producing Escherichia coli

One of the main strategies to improve the production of relevant metabolites has been the manipulation of single or multiple key genes in the metabolic pathways. This kind of strategy requires several rounds of experiments to identify enzymes that impact either yield or productivity. The use of mathematical tools to facilitate this process is desirable. In this work, we apply the Ensemble Modeling (EM) framework, which uses phenotypic data (effects of enzyme overexpression or genetic knockouts on the steady-state production rate) to screen for potential models capable of describe existing data and thus gaining insight to improve strains for l-lysine production. Described herein is a strategy to generate a set of kinetic models that describe a set of enzyme overexpression phenotypes previously determined in an Escherichia coli strain that produces increased levels of l-lysine in an industrial laboratory. This final ensemble of models captures the kinetic characteristics of the cell through screening of phenotypes after sequential overexpression of enzymes. Furthermore, these models demonstrate some predictive capability, as starting from the reference producing strain (overexpressing desensitized dihydrodipicolinate synthetase (dapA*)) this set of models is able to predict that the desensitization of aspartate kinase (lysC*) is the next rate-controlling step in the l-lysine pathway. Moreover, this set of models allows for the generation of further targets for testing, for example, phosphoenolpyruvate (Ppc), aspartate aminotransferase (AspC), and glutamate dehydrogenase (GdhA). This work demonstrates the usefulness, applicability, and scope that the Ensemble Modeling framework offers to build production strains.     

Combining mechanistic and data‐driven approaches to gain process knowledge on the control of the metabolic shift to lactate uptake in a fed‐batch CHO process

A growing body of knowledge is available on the cellular regulation of overflow metabolism in mammalian hosts of recombinant protein production. However, to develop strategies to control the regulation of overflow metabolism in cell culture processes, the effect of process parameters on metabolism has to be well understood. In this study, we investigated the effect of pH and temperature shift timing on lactate metabolism in a fed‐batch Chinese hamster ovary (CHO) process by using a Design of Experiments (DoE) approach. The metabolic switch to lactate consumption was controlled in a broad range by the proper timing of pH and temperature shifts. To extract process knowledge from the large experimental dataset, we proposed a novel methodological concept and demonstrated its usefulness with the analysis of lactate metabolism. Time‐resolved metabolic flux analysis and PLS‐R VIP were combined to assess the correlation of lactate metabolism and the activity of the major intracellular pathways. Whereas the switch to lactate uptake was mainly triggered by the decrease in the glycolytic flux, lactate uptake was correlated to TCA activity in the last days of the cultivation. These metabolic interactions were visualized on simple mechanistic plots to facilitate the interpretation of the results. Taken together, the combination of knowledge‐based mechanistic modeling and data‐driven multivariate analysis delivered valuable insights into the metabolic control of lactate production and has proven to be a powerful tool for the analysis of large metabolic datasets. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:1657–1668, 2015     

Lethality and synthetic lethality in the genome-wide metabolic network of Escherichia coli

Recent genomic analyses on the cellular metabolic network show that reaction flux across enzymes are diverse and exhibit power-law behavior in its distribution. While intuition might suggest that the reactions with larger fluxes are more likely to be lethal under the blockade of its catalysing gene products or gene knockouts, we find, by in silico flux analysis, that the lethality rarely has correlations with the flux level owing to the widespread backup pathways innate in the genome-wide metabolism of Escherichia coli. Lethal reactions, of which the deletion generates cascading failure of following reactions up to the biomass reaction, are identified in terms of the Boolean network scheme as well as the flux balance analysis. The avalanche size of a reaction, defined as the number of subsequently blocked reactions after its removal, turns out to be a useful measure of lethality. As a means to elucidate phenotypic robustness to a single deletion, we investigate synthetic lethality in reaction level, where simultaneous deletion of a pair of nonlethal reactions leads to the failure of the biomass reaction. Synthetic lethals identified via flux balance and Boolean scheme are consistently shown to act in parallel pathways, working in such a way that the backup machinery is compromised.