Comprehensive stable-isotope tracing of glucose and amino acids identifies metabolic by-products and their sources in CHO cell culture.

Mammalian cell culture processes are widely utilized for biotherapeutics production, disease diagnostics, and biosensors, and hence, should be optimized to support robust cell growth and viability. However, toxic by-products accumulate in cultures due to inefficiencies in metabolic activities and nutrient utilization. In this study, we applied comprehensive 13C stable-isotope tracing of amino acids and glucose to two Immunoglobulin G (IgG) producing Chinese Hamster Ovary (CHO) cell lines to identify secreted by-products and trace their origins. CHO cells were cultured in media formulations missing a single amino acid or glucose supplemented with a 13C-tracer of the missing substrate, followed by gas chromatography-mass spectrometry (GC-MS) analysis to track labeled carbon flows and identify by-products. We tracked the sources of all secreted by-products and verified the identity of 45 by-products, majority of which were derived from glucose, leucine, isoleucine, valine, tyrosine, tryptophan, methionine, and phenylalanine. In addition to by-products identified previously, we identified several metabolites including 2-hydroxyisovaleric acid, 2-aminobutyric acid, L-alloisoleucine, ketoisoleucine, 2-hydroxy-3-methylvaleric acid, desmeninol, and 2-aminobutyric acid. When added to CHO cell cultures at different concentrations, certain metabolites inhibited cell growth while others including 2-hydroxy acids, surprisingly, reduced lactate accumulation. In vitro enzymatic analysis indicated that 2-hydroxy acids were metabolized by lactate dehydrogenase suggesting a possible mechanism for lowered lactate accumulation, e.g., competitive substrate inhibition. The 13C-labeling assisted metabolomics pipeline developed and the metabolites identified will serve as a springboard to reduce undesirable by-products accumulation and alleviate inefficient substrate utilization in mammalian cultures used for biomanufacturing and other applications through altered media formulations and pathway engineering strategies.

Program for Integration and Rapid Analysis of Mass Isotopomer Distributions (PIRAMID).

SUMMARY: The analysis of stable isotope labeling experiments requires accurate, efficient, and reproducible quantification of mass isotopomer distributions (MIDs), which is not a core feature of general-purpose metabolomics software tools that are optimized to quantify metabolite abundance. Here, we present PIRAMID (Program for Integration and Rapid Analysis of Mass Isotopomer Distributions), a MATLAB-based tool that addresses this need by offering a user-friendly, graphical user interface-driven program to automate the extraction of isotopic information from mass spectrometry (MS) datasets. This tool can simultaneously extract ion chromatograms for various metabolites from multiple data files in common vendor-agnostic file formats, locate chromatographic peaks based on a targeted list of characteristic ions and retention times, and integrate MIDs for each target ion. These MIDs can be corrected for natural isotopic background based on the user-defined molecular formula of each ion. PIRAMID offers support for datasets acquired from low- or high-resolution MS, and single (MS) or tandem (MS/MS) instruments. It also enables the analysis of single or dual labeling experiments using a variety of isotopes (i.e. 2H, 13C, 15N, 18O, 34S). DATA AVAILABILITY AND IMPLEMENTATION: MATLAB p-code files are freely available for non-commercial use and can be downloaded from https://mfa.vueinnovations.com/. Commercial licenses are also available. All the data presented in this publication are available under the "Help_menu" folder of the PIRAMID software.

Ala-Cys-Cys-Ala dipeptide dimer alleviates problematic cysteine and cystine levels in media formulations and enhances CHO cell growth and metabolism.

Cysteine and cystine are essential amino acids present in mammalian cell cultures. While contributing to biomass synthesis, recombinant protein production, and antioxidant defense mechanisms, cysteine poses a major challenge in media formulations owing to its poor stability and oxidation to cystine, a cysteine dimer. Due to its poor solubility, cystine can cause precipitation of feed media, formation of undesired products, and consequently, reduce cysteine bioavailability. In this study, a highly soluble cysteine containing dipeptide dimer, Ala-Cys-Cys-Ala (ACCA), was evaluated as a suitable alternative to cysteine and cystine in CHO cell cultures. Replacing cysteine and cystine in basal medium with ACCA did not sustain cell growth. However, addition of ACCA at 4 mM and 8 mM to basal medium containing cysteine and cystine boosted cell growth up to 15% and 27% in CHO-GS and CHO-K1 batch cell cultures respectively and led to a proportionate increase in IgG titer. 13C-Metabolic flux analysis revealed that supplementation of ACCA reduced glycolytic fluxes by 20% leading to more efficient glucose metabolism in CHO-K1 cells. In fed-batch cultures, ACCA was able to replace cysteine and cystine in feed medium. Furthermore, supplementation of ACCA at high concentrations in basal medium eliminated the need for any cysteine equivalents in feed medium and increased cell densities and viabilities in fed-batch cultures without any significant impact on IgG charge variants. Taken together, this study demonstrates the potential of ACCA to improve CHO cell growth, productivity, and metabolism while also facilitating the formulation of cysteine- and cystine-free feed media. Such alternatives to cysteine and cystine will pave the way for enhanced biomanufacturing by increasing cell densities in culture and extending the storage of highly concentrated feed media as part of achieving intensified bioproduction processes.

Elucidating uptake and metabolic fate of dipeptides in CHO cell cultures using 13C labeling experiments and kinetic modeling.

The rapidly growing market of biologics including monoclonal antibodies has stimulated the need to improve biomanufacturing processes including mammalian host systems such as Chinese Hamster Ovary (CHO) cells. Cell culture media formulations continue to be enhanced to enable intensified cell culture processes and optimize cell culture performance. Amino acids, major components of cell culture media, are consumed in large amounts by CHO cells. Due to their low solubility and poor stability, certain amino acids including tyrosine, leucine, and phenylalanine can pose major challenges leading to suboptimal bioprocess performance. Dipeptides have the potential to replace amino acids in culture media. However, very little is known about the cleavage, uptake, and utilization kinetics of dipeptides in CHO cell cultures. In this study, replacing amino acids, including leucine and tyrosine by their respective dipeptides including but not limited to Ala-Leu and Gly-Tyr, supported similar cell growth, antibody production, and lactate profiles. Using 13C labeling techniques and spent media studies, dipeptides were shown to undergo both intracellular and extracellular cleavage in cultures. Extracellular cleavage increased with the culture duration, indicating cleavage by host cell proteins that are likely secreted and accumulate in cell culture over time. A kinetic model was built and for the first time, integrated with 13C labeling experiments to estimate dipeptide utilization rates, in CHO cell cultures. Dipeptides with alanine at the N-terminus had a higher utilization rate than dipeptides with alanine at the C-terminus and dipeptides with glycine instead of alanine at N-terminus. Simultaneous supplementation of more than one dipeptide in culture led to reduction in individual dipeptide utilization rates indicating that dipeptides compete for the same cleavage enzymes, transporters, or both. Dipeptide utilization rates in culture and cleavage rates in cell-free experiments appeared to follow Michaelis-Menten kinetics, reaching a maximum at higher dipeptide concentrations. Dipeptide utilization behavior was found to be similar in cell-free and cell culture environments, paving the way for future testing approaches for dipeptides in cell-free environments prior to use in large-scale bioreactors. Thus, this study provides a deeper understanding of the fate of dipeptides in CHO cell cultures through an integration of cell culture, 13C labeling, and kinetic modeling approaches providing insights in how to best use dipeptides in media formulations for robust and optimal mammalian cell culture performance.

13C-Metabolic flux analysis of 3T3-L1 adipocytes illuminates its core metabolism under hypoxia.

Hypoxia has been identified as a major factor in the pathogenesis of adipose tissue inflammation, which is a hallmark of obesity and obesity-linked type 2 diabetes mellitus. In this study, we have investigated the impact of hypoxia (1% oxygen) on the physiology and metabolism of 3T3-L1 adipocytes, a widely used cell culture model of adipose. Specifically, we applied parallel labeling experiments, isotopomer spectral analysis, and 13C-metabolic flux analysis to quantify the impact of hypoxia on adipogenesis, de novo lipogenesis and metabolic flux reprogramming in adipocytes. We found that 3T3-L1 cells can successfully differentiate into lipid-accumulating adipocytes under hypoxia, although the production of lipids was reduced by about 40%. Quantitative flux analysis demonstrated that short-term (1 day) and long-term (7 days) exposure to hypoxia resulted in similar reprogramming of cellular metabolism. Overall, we found that hypoxia: 1) reduced redox and energy generation by more than 2-fold and altered the patterns of metabolic pathway contributions to production and consumption of energy and redox cofactors; 2) redirected glucose metabolism from pentose phosphate pathway and citric acid cycle to lactate production; 3) rewired glutamine metabolism, from net glutamine production to net glutamine catabolism; 4) suppressed branched chain amino acid consumption; and 5) reduced biosynthesis of odd-chain fatty acids and mono-unsaturated fatty acids, while synthesis of saturated even-chain fatty acids was not affected. Together, these results highlight the profound impact of extracellular microenvironment on adipocyte metabolic activity and function.

Chemical inhibitors of hexokinase-2 enzyme reduce lactate accumulation, alter glycosylation processing, and produce altered glycoforms in CHO cell cultures.

Chinese hamster ovary (CHO) cells, predominant hosts for recombinant biotherapeutics production, generate lactate as a major glycolysis by-product. High lactate levels adversely impact cell growth and productivity. The goal of this study was to reduce lactate in CHO cell cultures by adding chemical inhibitors to hexokinase-2 (HK2), the enzyme catalyzing the conversion of glucose to glucose 6-phosphate, and examine their impact on lactate accumulation, cell growth, protein titers, and N-glycosylation. Five inhibitors of HK2 enzyme at different concentrations were evaluated, of which 2-deoxy- d-glucose (2DG) and 5-thio- d-glucose (5TG) successfully reduced lactate accumulation with only limited impacts on CHO cell growth. Individual 2DG and 5TG supplementation led to a 35%-45% decrease in peak lactate, while their combined supplementation resulted in a 60% decrease in peak lactate. Inhibitor supplementation led to at least 50% decrease in moles of lactate produced per mol of glucose consumed. Recombinant EPO-Fc titers peaked earlier relative to the end of culture duration in supplemented cultures leading to at least 11% and as high as 32% increase in final EPO-Fc titers. Asparagine, pyruvate, and serine consumption rates also increased in the exponential growth phase in 2DG and 5TG treated cultures, thus, rewiring central carbon metabolism due to low glycolytic fluxes. N-glycan analysis of EPO-Fc revealed an increase in high mannose glycans from 5% in control cultures to 25% and 37% in 2DG and 5TG-supplemented cultures, respectively. Inhibitor supplementation also led to a decrease in bi-, tri-, and tetra-antennary structures and up to 50% lower EPO-Fc sialylation. Interestingly, addition of 2DG led to the incorporation of 2-deoxy-hexose (2DH) on EPO-Fc N-glycans and addition of 5TG resulted in the first-ever observed N-glycan incorporation of 5-thio-hexose (5TH). Six percent to 23% of N-glycans included 5TH moieties, most likely 5-thio-mannose and/or 5-thio-galactose and/or possibly 5-thio-N-acetylglucosamine, and 14%-33% of N-glycans included 2DH moieties, most likely 2-deoxy-mannose and/or 2-deoxy-galactose, for cultures treated with different concentrations of 5TG and 2DG, respectively. Our study is the first to evaluate the impact of these glucose analogs on CHO cell growth, protein production, cell metabolism, N-glycosylation processing, and formation of alternative glycoforms.

Coordinated reprogramming of metabolism and cell function in adipocytes from proliferation to differentiation.

Adipose tissue plays a major role in regulating lipid and energy homeostasis by storing excess nutrients, releasing energetic substrates through lipolysis, and regulating metabolism of other tissues and organs through endocrine and paracrine signaling. Adipocytes within fat tissues store excess nutrients through increased cell number (hyperplasia), increased cell size (hypertrophy), or both. The differentiation of pre-adipocytes into mature lipid-accumulating adipocytes requires a complex interaction of metabolic pathways that is still incompletely understood. Here, we applied parallel labeling experiments and 13C-metabolic flux analysis to quantify precise metabolic fluxes in proliferating and differentiated 3T3-L1 cells, a widely used model to study adipogenesis. We found that morphological and biomass composition changes in adipocytes were accompanied by significant shifts in metabolic fluxes, encompassing all major metabolic pathways. In contrast to proliferating cells, differentiated adipocytes 1) increased glucose uptake and redirected glucose utilization from lactate production to lipogenesis and energy generation; 2) increased pathway fluxes through glycolysis, oxidative pentose phosphate pathway and citric acid cycle; 3) reduced lactate secretion, resulting in increased ATP generation via oxidative phosphorylation; 4) rewired glutamine metabolism, from glutaminolysis to de novo glutamine synthesis; 5) increased cytosolic NADPH production, driven mostly by increased cytosolic malic enzyme flux; 6) increased production of monounsaturated C16:1; and 7) activated a mitochondrial pyruvate cycle through simultaneous activity of pyruvate carboxylase, malate dehydrogenase and malic enzyme. Taken together, these results quantitatively highlight the complex interplay between pathway fluxes and cell function in adipocytes, and suggest a functional role for metabolic reprogramming in adipose differentiation and lipogenesis.

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 CO2 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 13C-metabolic flux analysis (13C-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 13C-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. 13C-MFA during mixotrophic growth using the parallel tracers [1-13C]fructose, [1,2-13C]fructose, [1,2,3-13C]fructose, and [U-13C]asparagine found that externally supplied CO2 contributed the majority of carbon consumed. All internally-produced CO2 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.

A guide to metabolic flux analysis in metabolic engineering: Methods, tools and applications.

The field of metabolic engineering is primarily concerned with improving the biological production of value-added chemicals, fuels and pharmaceuticals through the design, construction and optimization of metabolic pathways, redirection of intracellular fluxes, and refinement of cellular properties relevant for industrial bioprocess implementation. Metabolic network models and metabolic fluxes are central concepts in metabolic engineering, as was emphasized in the first paper published in this journal, "Metabolic fluxes and metabolic engineering" (Metabolic Engineering, 1: 1-11, 1999). In the past two decades, a wide range of computational, analytical and experimental approaches have been developed to interrogate the capabilities of biological systems through analysis of metabolic network models using techniques such as flux balance analysis (FBA), and quantify metabolic fluxes using constrained-based modeling approaches such as metabolic flux analysis (MFA) and more advanced experimental techniques based on the use of stable-isotope tracers, i.e. 13C-metabolic flux analysis (13C-MFA). In this review, we describe the basic principles of metabolic flux analysis, discuss current best practices in flux quantification, highlight potential pitfalls and alternative approaches in the application of these tools, and give a broad overview of pragmatic applications of flux analysis in metabolic engineering practice.

Regulatory interventions improve the biosynthesis of limiting amino acids from methanol carbon to improve synthetic methylotrophy in Escherichia coli.

Synthetic methylotrophy aims to engineer methane and methanol utilization pathways in platform hosts like Escherichia coli for industrial bioprocessing of natural gas and biogas. While recent attempts to engineer synthetic methylotrophs have proved successful, autonomous methylotrophy, that is, the ability to utilize methane or methanol as sole carbon and energy substrates, has not yet been realized. Here, we address an important limitation of autonomous methylotrophy in E. coli: the inability of the organism to synthesize several amino acids when grown on methanol. We targeted global and local amino acid regulatory networks. Those include removal of amino acid allosteric feedback inhibition (argAH15Y , ilvAL447F , hisGE271K , leuAG462D , proBD107N , thrAS345F , trpES40F ), knockouts of transcriptional repressors (ihfA, metJ); and overexpression of amino acid biosynthetic operons (hisGDCBHAFI, leuABCD, thrABC, trpEDCBA) and transcriptional regulators (crp, purR). Compared to the parent methylotrophic E. coli strain that was unable to synthesize these amino acids from methanol carbon, these strategies resulted in improved biosynthesis of limiting proteinogenic amino acids (histidine, leucine, lysine, methionine, phenylalanine, threonine, tyrosine) from methanol carbon. In several cases, improved amino acid biosynthesis from methanol carbon led to improvements in methylotrophic growth in methanol minimal medium supplemented with a small amount of yeast extract. This study addresses a key limitation currently preventing autonomous methylotrophy in E. coli and possibly other synthetic methylotrophs and provides insight as to how this limitation can be alleviated via global and local regulatory modifications.

Adaptive laboratory evolution of methylotrophic Escherichia coli enables synthesis of all amino acids from methanol-derived carbon.

Recent attempts to create synthetic Escherichia coli methylotrophs identified that de novo biosynthesis of amino acids, in the presence of methanol, presents significant challenges in achieving autonomous methylotrophic growth. Previously engineered methanol-dependent strains required co-utilization of stoichiometric amounts of co-substrates and methanol. As such, these strains could not be evolved to grow on methanol alone. In this work, we have explored an alternative approach to enable biosynthesis of all amino acids from methanol-derived carbon in minimal media without stoichiometric coupling. First, we identified that biosynthesis of threonine was limiting the growth of our methylotrophic E. coli. To address this, we performed adaptive laboratory evolution to generate a strain that grew efficiently in minimal medium with methanol and threonine. Methanol assimilation and growth of the evolved strain were analyzed, and, interestingly, we found that the evolved strain synthesized all amino acids, including threonine, from methanol-derived carbon. The evolved strain was then further engineered through overexpression of an optimized threonine biosynthetic pathway. We show that the resulting methylotrophic E. coli strain has a methanol-dependent growth phenotype with homoserine as co-substrate. In contrast to previous methanol-dependent strains, co-utilization of homoserine is not stoichiometrically linked to methanol assimilation. As such, future engineering of this strain and successive adaptive evolution could enable autonomous growth on methanol as the sole carbon source. KEY POINTS: • Adaptive evolution of E. coli enables biosynthesis of all amino acids from methanol. • Overexpression of threonine biosynthesis pathway improves methanol assimilation. • Methanol-dependent growth is seen in minimal media with homoserine as co-substrate.

Dissecting the genetic and metabolic mechanisms of adaptation to the knockout of a major metabolic enzyme in Escherichia coli.

Unraveling the mechanisms of microbial adaptive evolution following genetic or environmental challenges is of fundamental interest in biological science and engineering. When the challenge is the loss of a metabolic enzyme, adaptive responses can also shed significant insight into metabolic robustness, regulation, and areas of kinetic limitation. In this study, whole-genome sequencing and high-resolution 13C-metabolic flux analysis were performed on 10 adaptively evolved pgi knockouts of Escherichia coliPgi catalyzes the first reaction in glycolysis, and its loss results in major physiological and carbon catabolism pathway changes, including an 80% reduction in growth rate. Following adaptive laboratory evolution (ALE), the knockouts increase their growth rate by up to 3.6-fold. Through combined genomic-fluxomic analysis, we characterized the mutations and resulting metabolic fluxes that enabled this fitness recovery. Large increases in pyridine cofactor transhydrogenase flux, correcting imbalanced production of NADPH and NADH, were enabled by direct mutations to the transhydrogenase genes sthA and pntAB The phosphotransferase system component crr was also found to be frequently mutated, which corresponded to elevated flux from pyruvate to phosphoenolpyruvate. The overall energy metabolism was found to be strikingly robust, and what have been previously described as latently activated Entner-Doudoroff and glyoxylate shunt pathways are shown here to represent no real increases in absolute flux relative to the wild type. These results indicate that the dominant mechanism of adaptation was to relieve the rate-limiting steps in cofactor metabolism and substrate uptake and to modulate global transcriptional regulation from stress response to catabolism.

Improving synthetic methylotrophy via dynamic formaldehyde regulation of pentose phosphate pathway genes and redox perturbation.

Escherichia coli is an ideal choice for constructing synthetic methylotrophs capable of utilizing the non-native substrate methanol as a carbon and energy source. All current E. coli-based synthetic methylotrophs require co-substrates. They display variable levels of methanol-carbon incorporation due to a lack of native regulatory control of biosynthetic pathways, as E. coli does not recognize methanol as a proper substrate despite its ability to catabolize it. Here, using the E. coli formaldehyde-inducible promoter Pfrm, we implement dynamic expression control of select pentose-phosphate genes in response to the formaldehyde produced upon methanol oxidation. Genes under Pfrm control exhibited 8- to 30-fold transcriptional upregulation during growth on methanol. Formaldehyde-induced episomal expression of the B. methanolicus rpe and tkt genes involved in the regeneration of ribulose 5-phosphate required for formaldehyde fixation led to significantly improved methanol assimilation into intracellular metabolites, including a 2-fold increase of 13C-methanol into glutamate. Using a simple strategy for redox perturbation by deleting the E. coli NAD-dependent malate dehydrogenase gene maldh, we demonstrate 5-fold improved biomass formation of cells growing on methanol in the presence of a small concentration of yeast extract. Further improvements in methanol utilization are achieved via adaptive laboratory evolution and heterologous rpe and tkt expression. A short-term in vivo13C-methanol labeling assay was used to determine methanol assimilation activity for Δmaldh strains, and demonstrated dramatically higher labeling in intracellular metabolites, including a 6-fold and 1.8-fold increase in glycine labeling for the rpe/tkt and evolved strains, respectively. The combination of formaldehyde-controlled pentose phosphate pathway expression and redox perturbation with the maldh knock-out greatly improved both growth benefit with methanol and methanol carbon incorporation into intracellular metabolites.

Triggering the stringent response enhances synthetic methanol utilization in Escherichia coli.

Synthetic methylotrophy aims to engineer methane and methanol utilization pathways in platform hosts like Escherichia coli for industrial bioprocessing of natural gas and biogas. While recent attempts to engineer synthetic methylotrophs have proved successful, autonomous methylotrophy, i.e. the ability to utilize methane or methanol as sole carbon and energy substrates, has not yet been realized. Here, we address an important limitation of autonomous methylotrophy in E. coli: the inability of the organism to synthesize several amino acids when grown on methanol. By activating the stringent/stress response via ppGpp overproduction, or DksA and RpoS overexpression, we demonstrate improved biosynthesis of proteinogenic amino acids via endogenous upregulation of amino acid synthesis pathway genes. Thus, we were able to achieve biosynthesis of several limiting amino acids from methanol-derived carbon, in contrast to the control methylotrophic E. coli strain. This study addresses a key limitation currently preventing autonomous methylotrophy in E. coli and possibly other synthetic methylotrophs and provides insight as to how this limitation can be alleviated via stringent/stress response activation.

Engineering Escherichia coli for methanol-dependent growth on glucose for metabolite production.

Synthetic methylotrophy aims to engineer methane and methanol utilization pathways in platform hosts like Escherichia coli for industrial bioprocessing of natural gas and biogas. While recent attempts to engineer synthetic methanol auxotrophs have proved successful, these studies focused on scarce and expensive co-substrates. Here, we engineered E. coli for methanol-dependent growth on glucose, an abundant and inexpensive co-substrate, via deletion of glucose 6-phosphate isomerase (pgi), phosphogluconate dehydratase (edd), and ribose 5-phosphate isomerases (rpiAB). Since the parental strain did not exhibit methanol-dependent growth on glucose in minimal medium, we first achieved methanol-dependent growth via amino acid supplementation and used this medium to evolve the strain for methanol-dependent growth in glucose minimal medium. The evolved strain exhibited a maximum growth rate of 0.15 h-1 in glucose minimal medium with methanol, which is comparable to that of other synthetic methanol auxotrophs. Whole genome sequencing and 13C-metabolic flux analysis revealed the causative mutations in the evolved strain. A mutation in the phosphotransferase system enzyme I gene (ptsI) resulted in a reduced glucose uptake rate to maintain a one-to-one molar ratio of substrate utilization. Deletion of the e14 prophage DNA region resulted in two non-synonymous mutations in the isocitrate dehydrogenase (icd) gene, which reduced TCA cycle carbon flux to maintain the internal redox state. In high cell density glucose fed-batch fermentation, methanol-dependent acetone production resulted in 22% average carbon labeling of acetone from 13C-methanol, which far surpasses that of the previous best (2.4%) found with methylotrophic E. coli Δpgi. This study addresses the need to identify appropriate co-substrates for engineering synthetic methanol auxotrophs and provides a basis for the next steps toward industrial one-carbon bioprocessing.

From Escherichia coli mutant 13C labeling data to a core kinetic model: A kinetic model parameterization pipeline.

Kinetic models of metabolic networks offer the promise of quantitative phenotype prediction. The mechanistic characterization of enzyme catalyzed reactions allows for tracing the effect of perturbations in metabolite concentrations and reaction fluxes in response to genetic and environmental perturbation that are beyond the scope of stoichiometric models. In this study, we develop a two-step computational pipeline for the rapid parameterization of kinetic models of metabolic networks using a curated metabolic model and available 13C-labeling distributions under multiple genetic and environmental perturbations. The first step involves the elucidation of all intracellular fluxes in a core model of E. coli containing 74 reactions and 61 metabolites using 13C-Metabolic Flux Analysis (13C-MFA). Here, fluxes corresponding to the mid-exponential growth phase are elucidated for seven single gene deletion mutants from upper glycolysis, pentose phosphate pathway and the Entner-Doudoroff pathway. The computed flux ranges are then used to parameterize the same (i.e., k-ecoli74) core kinetic model for E. coli with 55 substrate-level regulations using the newly developed K-FIT parameterization algorithm. The K-FIT algorithm employs a combination of equation decomposition and iterative solution techniques to evaluate steady-state fluxes in response to genetic perturbations. k-ecoli74 predicted 86% of flux values for strains used during fitting within a single standard deviation of 13C-MFA estimated values. By performing both tasks using the same network, errors associated with lack of congruity between the two networks are avoided, allowing for seamless integration of data with model building. Product yield predictions and comparison with previously developed kinetic models indicate shifts in flux ranges and the presence or absence of mutant strains delivering flux towards pathways of interest from training data significantly impact predictive capabilities. Using this workflow, the impact of completeness of fluxomic datasets and the importance of specific genetic perturbations on uncertainties in kinetic parameter estimation are evaluated.

Metabolic flux responses to deletion of 20 core enzymes reveal flexibility and limits of E. coli metabolism.

Despite remarkable progress in mapping biochemistry and gene-protein-reaction relationships, quantitative systems-level understanding of microbial metabolism remains a persistent challenge. Here, 13C-metabolic flux analysis was applied to interrogate metabolic responses to 20 genetic perturbations in all viable Escherichia coli single gene knockouts in upper central metabolic pathways. Strains with severe growth defects displayed highly altered intracellular flux patterns and were the most difficult to predict using current constraint-based modeling approaches. In the ΔpfkA strain, an unexpected glucose-secretion phenotype was identified. The broad range of flux rewiring responses that were quantified suggest that some compensating pathways are more flexible than others, resulting in a more robust physiology. The fact that only 2 out of 20 strains displayed an increased net pathway-flux capacity points to a fundamental rate limitation of E. coli core metabolism. In cataloguing the various cellular responses, our results provide a critical resource for kinetic model development and efforts focused on genotype-to-phenotype predictions.

Deletion of four genes in Escherichia coli enables preferential consumption of xylose and secretion of glucose.

Overcoming carbon catabolite repression presents a significant challenge, largely due to the complex regulatory networks governing substrate catabolism, even in microbial cells. In this work, we have engineered an E. coli strain, which we have named X2G, that not only exhibits a reversed substrate preference for xylose over glucose, but also demonstrates an unusual ability to produce significant amounts of glucose. We obtained this non-intuitive phenotype by deleting four genes in upper central metabolism: ptsI, glk, pfkA, and zwf, which respectively encode Enzyme I of the phosphotransferase system, glucokinase, the dominant isozyme of phosphofructokinase, and glucose-6-phosphate dehydrogenase. The deletion of ptsI and glk blocks glucose uptake in E. coli, while the deletion of pfkA and zwf prevents the reassimilation of carbons through glycolysis and the oxidative pentose phosphate pathway, respectively. Our strain X2G is capable of converting 34% of the carbon it takes up as xylose into exported glucose. This corresponds to a glucose production rate of 1.4 ±â€¯0.3 mmol/gDW/h at a specific growth rate of 0.25 ±â€¯0.03 h-1, or about 1.8 ±â€¯0.1 mM of glucose accumulated for every unit increase in OD600. Despite a 22% decrease in xylose uptake rate, a 33% decrease in biomass yield, and a 52% decrease in acetate production rate relative to the wild-type, the intracellular flux profile and cofactor allocation of X2G remain largely unperturbed, as elucidated through 13C-metabolic flux analysis. Further quantification of the pool sizes of key intracellular metabolites revealed that glucose secretion by X2G is likely driven by the substantial accumulation of intracellular glucose 6-phosphate, fructose 6-phosphate, glucose and fructose at levels greater than 20x of that in wild-type E. coli. Combined, our results shed light on the flexibility of central metabolism, and the opportunities this affords for producing value-added pentose- and hexose-derived products from lignocellulosic feedstocks.

Predicting Dynamic Metabolic Demands in the Photosynthetic Eukaryote Chlorella vulgaris.

Phototrophic organisms exhibit a highly dynamic proteome, adapting their biomass composition in response to diurnal light/dark cycles and nutrient availability. Here, we used experimentally determined biomass compositions over the course of growth to determine and constrain the biomass objective function (BOF) in a genome-scale metabolic model of Chlorella vulgaris UTEX 395 over time. Changes in the BOF, which encompasses all metabolites necessary to produce biomass, influence the state of the metabolic network thus directly affecting predictions. Simulations using dynamic BOFs predicted distinct proteome demands during heterotrophic or photoautotrophic growth. Model-driven analysis of extracellular nitrogen concentrations and predicted nitrogen uptake rates revealed an intracellular nitrogen pool, which contains 38% of the total nitrogen provided in the medium for photoautotrophic and 13% for heterotrophic growth. Agreement between flux and gene expression trends was determined by statistical comparison. Accordance between predicted flux trends and gene expression trends was found for 65% of multisubunit enzymes and 75% of allosteric reactions. Reactions with the highest agreement between simulations and experimental data were associated with energy metabolism, terpenoid biosynthesis, fatty acids, nucleotides, and amino acid metabolism. Furthermore, predicted flux distributions at each time point were compared with gene expression data to gain new insights into intracellular compartmentalization, specifically for transporters. A total of 103 genes related to internal transport reactions were identified and added to the updated model of C. vulgaris, iCZ946, thus increasing our knowledgebase by 10% for this model green alga.

Metabolism in dense microbial colonies: 13C metabolic flux analysis of E. coli grown on agar identifies two distinct cell populations with acetate cross-feeding.

In this study, we have investigated for the first time the metabolism of E. coli grown on agar using 13C metabolic flux analysis (13C-MFA). To date, all 13C-MFA studies on microbes have been performed with cells grown in liquid culture. Here, we extend the scope of 13C-MFA to biological systems where cells are grown in dense microbial colonies. First, we identified new optimal 13C tracers to quantify fluxes in systems where the acetate yield cannot be easily measured. We determined that three parallel labeling experiments with the tracers [1,2-13C]glucose, [1,6-13C]glucose, and [4,5,6-13C]glucose permit precise estimation of not only intracellular fluxes, but also of the amount of acetate produced from glucose. Parallel labeling experiments were then performed with wild-type E. coli and E. coli ΔackA grown in liquid culture and on agar plates. Initial attempts to fit the labeling data from wild-type E. coli grown on agar did not produce a statistically acceptable fit. To resolve this issue, we employed the recently developed co-culture 13C-MFA approach, where two E. coli subpopulations were defined in the model that engaged in metabolite cross-feeding. The flux results identified two distinct E. coli cell populations, a dominant cell population (92% of cells) that metabolized glucose via conventional metabolic pathways and secreted a large amount of acetate (~40% of maximum theoretical yield), and a second smaller cell population (8% of cells) that consumed the secreted acetate without any glucose influx. These experimental results are in good agreement with recent theoretical simulations. Importantly, this study provides a solid foundation for future investigations of a wide range of problems involving microbial biofilms that are of great interest in biotechnology, ecology and medicine, where metabolite cross-feeding between cell populations is a core feature of the communities.