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.

A citric acid cycle-deficient Escherichia coli as an efficient chassis for aerobic fermentations.

Tricarboxylic acid cycle (TCA cycle) plays an important role for aerobic growth of heterotrophic bacteria. Theoretically, eliminating TCA cycle would decrease carbon dissipation and facilitate chemicals biosynthesis. Here, we construct an E. coli strain without a functional TCA cycle that can serve as a versatile chassis for chemicals biosynthesis. We first use adaptive laboratory evolution to recover aerobic growth in minimal medium of TCA cycle-deficient E. coli. Inactivation of succinate dehydrogenase is a key event in the evolutionary trajectory. Supply of succinyl-CoA is identified as the growth limiting factor. By replacing endogenous succinyl-CoA dependent enzymes, we obtain an optimized TCA cycle-deficient E. coli strain. As a proof of concept, the strain is engineered for high-yield production of four separate products. This work enhances our understanding of the role of the TCA cycle in E. coli metabolism and demonstrates the advantages of using TCA cycle-deficient E. coli strain for biotechnological applications.

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.

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.

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.

Experimental Evolution Reveals Unifying Systems-Level Adaptations but Diversity in Driving Genotypes.

Genotype-fitness maps of evolution have been well characterized for biological components, such as RNA and proteins, but remain less clear for systems-level properties, such as those of metabolic and transcriptional regulatory networks. Here, we take multi-omics measurements of 6 different E. coli strains throughout adaptive laboratory evolution (ALE) to maximal growth fitness. The results show the following: (i) convergence in most overall phenotypic measures across all strains, with the notable exception of divergence in NADPH production mechanisms; (ii) conserved transcriptomic adaptations, describing increased expression of growth promoting genes but decreased expression of stress response and structural components; (iii) four groups of regulatory trade-offs underlying the adjustment of transcriptome composition; and (iv) correlates that link causal mutations to systems-level adaptations, including mutation-pathway flux correlates and mutation-transcriptome composition correlates. We thus show that fitness landscapes for ALE can be described with two layers of causation: one based on system-level properties (continuous variables) and the other based on mutations (discrete variables). IMPORTANCE Understanding the mechanisms of microbial adaptation will help combat the evolution of drug-resistant microbes and enable predictive genome design. Although experimental evolution allows us to identify the causal mutations underlying microbial adaptation, it remains unclear how causal mutations enable increased fitness and is often explained in terms of individual components (i.e., enzyme rate) as opposed to biological systems (i.e., pathways). Here, we find that causal mutations in E. coli are linked to systems-level changes in NADPH balance and expression of stress response genes. These systems-level adaptation patterns are conserved across diverse E. coli strains and thus identify cofactor balance and proteome reallocation as dominant constraints governing microbial adaptation.

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.

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.

Improving the Methanol Tolerance of an Escherichia coli Methylotroph via Adaptive Laboratory Evolution Enhances Synthetic Methanol Utilization.

There is great interest in developing synthetic methylotrophs that harbor methane and methanol utilization pathways in heterologous hosts such as Escherichia coli for industrial bioconversion of one-carbon compounds. While there are recent reports that describe the successful engineering of synthetic methylotrophs, additional efforts are required to achieve the robust methylotrophic phenotypes required for industrial realization. Here, we address an important issue of synthetic methylotrophy in E. coli: methanol toxicity. Both methanol, and its oxidation product, formaldehyde, are cytotoxic to cells. Methanol alters the fluidity and biological properties of cellular membranes while formaldehyde reacts readily with proteins and nucleic acids. Thus, efforts to enhance the methanol tolerance of synthetic methylotrophs are important. Here, adaptive laboratory evolution was performed to improve the methanol tolerance of several E. coli strains, both methylotrophic and non-methylotrophic. Serial batch passaging in rich medium containing toxic methanol concentrations yielded clones exhibiting improved methanol tolerance. In several cases, these evolved clones exhibited a > 50% improvement in growth rate and biomass yield in the presence of high methanol concentrations compared to the respective parental strains. Importantly, one evolved clone exhibited a two to threefold improvement in the methanol utilization phenotype, as determined via 13C-labeling, at non-toxic, industrially relevant methanol concentrations compared to the respective parental strain. Whole genome sequencing was performed to identify causative mutations contributing to methanol tolerance. Common mutations were identified in 30S ribosomal subunit proteins, which increased translational accuracy and provided insight into a novel methanol tolerance mechanism. This study addresses an important issue of synthetic methylotrophy in E. coli and provides insight as to how methanol toxicity can be alleviated via enhancing methanol tolerance. Coupled improvement of methanol tolerance and synthetic methanol utilization is an important advancement for the field of synthetic methylotrophy.

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.

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.

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.

A guide to deciphering microbial interactions and metabolic fluxes in microbiome communities.

Microbiomes occupy nearly all environments on Earth. These communities of interacting microorganisms are highly complex, dynamic biological systems that impact and reshape the molecular composition of their habitats by performing complex biochemical transformations. The structure and function of microbiomes are influenced by local environmental stimuli and spatiotemporal changes. In order to control the dynamics and ultimately the function of microbiomes, we need to develop a mechanistic and quantitative understanding of the ecological, molecular, and evolutionary driving forces that govern these systems. Here, we describe recent advances in developing computational and experimental approaches that can promote a more fundamental understanding of microbial communities through comprehensive model-based analysis of heterogeneous data types across multiple scales, from intracellular metabolism, to metabolite cross-feeding interactions, to the emergent macroscopic behaviors. Ultimately, harnessing the full potential of microbiomes for practical applications will require developing new predictive modeling approaches and better tools to manipulate microbiome interactions.

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.

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.

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.

High-resolution 13C metabolic flux analysis.

Precise quantification of metabolic pathway fluxes in biological systems is of major importance in guiding efforts in metabolic engineering, biotechnology, microbiology, human health, and cell culture. 13C metabolic flux analysis (13C-MFA) is the predominant technique used for determining intracellular fluxes. Here, we present a protocol for 13C-MFA that incorporates recent advances in parallel labeling experiments, isotopic labeling measurements, and statistical analysis, as well as best practices developed through decades of experience. Experimental design to ensure that fluxes are estimated with the highest precision is an integral part of the protocol. The protocol is based on growing microbes in two (or more) parallel cultures with 13C-labeled glucose tracers, followed by gas chromatography-mass spectrometry (GC-MS) measurements of isotopic labeling of protein-bound amino acids, glycogen-bound glucose, and RNA-bound ribose. Fluxes are then estimated using software for 13C-MFA, such as Metran, followed by comprehensive statistical analysis to determine the goodness of fit and calculate confidence intervals of fluxes. The presented protocol can be completed in 4 d and quantifies metabolic fluxes with a standard deviation of ≤2%, a substantial improvement over previous implementations. The presented protocol is exemplified using an Escherichia coli ΔtpiA case study with full supporting data, providing a hands-on opportunity to step through a complex troubleshooting scenario. Although applications to prokaryotic microbial systems are emphasized, this protocol can be easily adjusted for application to eukaryotic organisms.

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.

Synthetic methylotrophy: Strategies to assimilate methanol for growth and chemicals production.

Methanol is an attractive and broadly available substrate for large-scale bioproduction of fuels and chemicals. It contains more energy and electrons per carbon than carbohydrates and can be cheaply produced from natural gas. Synthetic methylotrophy refers to the development of non-native methylotrophs such as Escherichia coli and Corynebacterium glutamicum to utilize methanol as a carbon source. Here, we discuss recent advances in engineering these industrial hosts to assimilate methanol for growth and chemicals production through the introduction of the ribulose monophosphate (RuMP) cycle. In addition, we present novel strategies based on flux coupling and adaptive laboratory evolution to engineer new strains that can grow exclusively on methanol.