Evaluation of 1,2-diacyl-3-acetyl triacylglycerol production in Yarrowia lipolytica.

Plants produce many high-value oleochemical molecules. While oil-crop agriculture is performed at industrial scales, suitable land is not available to meet global oleochemical demand. Worse, establishing new oil-crop farms often comes with the environmental cost of tropical deforestation. The field of metabolic engineering offers tools to transplant oleochemical metabolism into tractable hosts while simultaneously providing access to molecules produced by non-agricultural plants. Here, we evaluate strategies for rewiring metabolism in the oleaginous yeast Yarrowia lipolytica to synthesize a foreign lipid, 3-acetyl-1,2-diacyl-sn-glycerol (acTAG). Oils made up of acTAG have a reduced viscosity and melting point relative to traditional triacylglycerol oils making them attractive as low-grade diesels, lubricants, and emulsifiers. This manuscript describes a metabolic engineering study that established acTAG production at g/L scale, exploration of the impact of lipid bodies on acTAG titer, and a techno-economic analysis that establishes the performance benchmarks required for microbial acTAG production to be economically feasible.

Comparative functional genomics identifies an iron-limited bottleneck in a Saccharomyces cerevisiae strain with a cytosolic-localized isobutanol pathway.

Metabolic engineering strategies have been successfully implemented to improve the production of isobutanol, a next-generation biofuel, in Saccharomyces cerevisiae. Here, we explore how two of these strategies, pathway re-localization and redox cofactor-balancing, affect the performance and physiology of isobutanol producing strains. We equipped yeast with isobutanol cassettes which had either a mitochondrial or cytosolic localized isobutanol pathway and used either a redox-imbalanced (NADPH-dependent) or redox-balanced (NADH-dependent) ketol-acid reductoisomerase enzyme. We then conducted transcriptomic, proteomic and metabolomic analyses to elucidate molecular differences between the engineered strains. Pathway localization had a large effect on isobutanol production with the strain expressing the mitochondrial-localized enzymes producing 3.8-fold more isobutanol than strains expressing the cytosolic enzymes. Cofactor-balancing did not improve isobutanol titers and instead the strain with the redox-imbalanced pathway produced 1.5-fold more isobutanol than the balanced version, albeit at low overall pathway flux. Functional genomic analyses suggested that the poor performances of the cytosolic pathway strains were in part due to a shortage in cytosolic Fe-S clusters, which are required cofactors for the dihydroxyacid dehydratase enzyme. We then demonstrated that this cofactor limitation may be partially recovered by disrupting iron homeostasis with a fra2 mutation, thereby increasing cellular iron levels. The resulting isobutanol titer of the fra2 null strain harboring a cytosolic-localized isobutanol pathway outperformed the strain with the mitochondrial-localized pathway by 1.3-fold, demonstrating that both localizations can support flux to isobutanol.

Assessing the impact of substrate-level enzyme regulations limiting ethanol titer in Clostridium thermocellum using a core kinetic model.

Clostridium thermocellum is a promising candidate for consolidated bioprocessing because it can directly ferment cellulose to ethanol. Despite significant efforts, achieved yields and titers fall below industrially relevant targets. This implies that there still exist unknown enzymatic, regulatory, and/or possibly thermodynamic bottlenecks that can throttle back metabolic flow. By (i) elucidating internal metabolic fluxes in wild-type C. thermocellum grown on cellobiose via 13C-metabolic flux analysis (13C-MFA), (ii) parameterizing a core kinetic model, and (iii) subsequently deploying an ensemble-docking workflow for discovering substrate-level regulations, this paper aims to reveal some of these factors and expand our knowledgebase governing C. thermocellum metabolism. Generated 13C labeling data were used with 13C-MFA to generate a wild-type flux distribution for the metabolic network. Notably, flux elucidation through MFA alluded to serine generation via the mercaptopyruvate pathway. Using the elucidated flux distributions in conjunction with batch fermentation process yield data for various mutant strains, we constructed a kinetic model of C. thermocellum core metabolism (i.e. k-ctherm138). Subsequently, we used the parameterized kinetic model to explore the effect of removing substrate-level regulations on ethanol yield and titer. Upon exploring all possible simultaneous (up to four) regulation removals we identified combinations that lead to many-fold model predicted improvement in ethanol titer. In addition, by coupling a systematic method for identifying putative competitive inhibitory mechanisms using K-FIT kinetic parameterization with the ensemble-docking workflow, we flagged 67 putative substrate-level inhibition mechanisms across central carbon metabolism supported by both kinetic formalism and docking analysis.

Metabolic Remodeling during Nitrogen Fixation in Zymomonas mobilis.

Zymomonas mobilis is an ethanologenic bacterium currently being developed for production of advanced biofuels. Recent studies have shown that Z. mobilis can fix dinitrogen gas (N2) as a sole nitrogen source. During N2 fixation, Z. mobilis exhibits increased biomass-specific rates of ethanol production. In order to better understand the physiology of Z. mobilis during N2 fixation and during changes in ammonium (NH4+) availability, we performed liquid chromatography-mass spectrometry (LC-MS)-based targeted metabolomics and shotgun proteomics under three regimes of nitrogen availability: continuous N2 fixation, gradual NH4+ depletion, and acute NH4+ addition to N2-fixing cells. We report dynamic changes in abundance of proteins and metabolites related to nitrogen fixation, motility, ammonium assimilation, amino acid biosynthesis, nucleotide biosynthesis, isoprenoid biosynthesis, and Entner-Doudoroff (ED) glycolysis, providing insight into the regulatory mechanisms that control these processes in Z. mobilis. Our analysis identified potential physiological mechanisms that may contribute to increased specific ethanol production during N2 fixation, including decreased activity of biosynthetic pathways, increased protein abundance of alcohol dehydrogenase (ADHI), and increased thermodynamic favorability of the ED pathway. Of particular relevance to advanced biofuel production, we found that intermediates in the methylerythritol phosphate (MEP) pathway for isoprenoid biosynthesis were depleted during N2 fixation, coinciding with decreased protein abundance of deoxyxylulose 5-phosphate synthase (DXS), the first enzyme in the pathway. This implies that DXS protein abundance serves as a native control point in regulating MEP pathway activity in Z. mobilis. The results of this study will inform metabolic engineering to further develop Z. mobilis as a platform organism for biofuel production. IMPORTANCE Biofuels and bioproducts have the potential to serve as environmentally sustainable replacements for petroleum-derived fuels and commodity molecules. Advanced fuels such as higher alcohols and isoprenoids are more suitable gasoline replacements than bioethanol. Developing microbial systems to generate advanced biofuels requires metabolic engineering to reroute carbon away from ethanol and other native products and toward desired pathways, such as the MEP pathway for isoprenoid biosynthesis. However, rational engineering of microbial metabolism relies on understanding metabolic control points, in terms of both enzyme activity and thermodynamic favorability. In Z. mobilis, the factors that control glycolytic rates, ethanol production, and isoprenoid production are still not fully understood. In this study, we performed metabolomic, proteomic, and thermodynamic analysis of Z. mobilis during N2 fixation. This analysis identified key changes in metabolite levels, enzyme abundance, and glycolytic thermodynamic favorability that occurred during changes in NH4+ availability, helping to inform future efforts in metabolic engineering.

Negative interactions determine Clostridioides difficile growth in synthetic human gut communities.

Understanding the principles of colonization resistance of the gut microbiome to the pathogen Clostridioides difficile will enable the design of defined bacterial therapeutics. We investigate the ecological principles of community resistance to C. difficile using a synthetic human gut microbiome. Using a dynamic computational model, we demonstrate that C. difficile receives the largest number and magnitude of incoming negative interactions. Our results show that C. difficile is in a unique class of species that display a strong negative dependence between growth and species richness. We identify molecular mechanisms of inhibition including acidification of the environment and competition over resources. We demonstrate that Clostridium hiranonis strongly inhibits C. difficile partially via resource competition. Increasing the initial density of C. difficile can increase its abundance in the assembled community, but community context determines the maximum achievable C. difficile abundance. Our work suggests that the C. difficile inhibitory potential of defined bacterial therapeutics can be optimized by designing communities featuring a combination of mechanisms including species richness, environment acidification, and resource competition.

Hostile Takeover: How Viruses Reprogram Prokaryotic Metabolism.

To reproduce, prokaryotic viruses must hijack the cellular machinery of their hosts and redirect it toward the production of viral particles. While takeover of the host replication and protein synthesis apparatus has long been considered an essential feature of infection, recent studies indicate that extensive reprogramming of host primary metabolism is a widespread phenomenon among prokaryotic viruses that is required to fulfill the biosynthetic needs of virion production. In this review we provide an overview of the most significant recent findings regarding virus-induced reprogramming of prokaryotic metabolism and suggest how quantitative systems biology approaches may be used to provide a holistic understanding of metabolic remodeling during lytic viral infection.

Investigating the Chemolithoautotrophic and Formate Metabolism of Nitrospira moscoviensis by Constraint-Based Metabolic Modeling and 13C-Tracer Analysis.

Nitrite-oxidizing bacteria belonging to the genus Nitrospira mediate a key step in nitrification and play important roles in the biogeochemical nitrogen cycle and wastewater treatment. While these organisms have recently been shown to exhibit metabolic flexibility beyond their chemolithoautotrophic lifestyle, including the use of simple organic compounds to fuel their energy metabolism, the metabolic networks controlling their autotrophic and mixotrophic growth remain poorly understood. Here, we reconstructed a genome-scale metabolic model for Nitrospira moscoviensis (iNmo686) and used flux balance analysis to evaluate the metabolic networks controlling autotrophic and formatotrophic growth on nitrite and formate, respectively. Subsequently, proteomic analysis and [13C]bicarbonate and [13C]formate tracer experiments coupled to metabolomic analysis were performed to experimentally validate model predictions. Our findings corroborate that N. moscoviensis uses the reductive tricarboxylic acid cycle for CO2 fixation, and we also show that N. moscoviensis can indirectly use formate as a carbon source by oxidizing it first to CO2 followed by reassimilation, rather than direct incorporation via the reductive glycine pathway. Our study offers the first measurements of Nitrospira's in vivo central carbon metabolism and provides a quantitative tool that can be used for understanding and predicting their metabolic processes. IMPORTANCE Nitrospira spp. are globally abundant nitrifying bacteria in soil and aquatic ecosystems and in wastewater treatment plants, where they control the oxidation of nitrite to nitrate. Despite their critical contribution to nitrogen cycling across diverse environments, detailed understanding of their metabolic network and prediction of their function under different environmental conditions remains a major challenge. Here, we provide the first constraint-based metabolic model of Nitrospira moscoviensis representing the ubiquitous Nitrospira lineage II and subsequently validate this model using proteomics and 13C-tracers combined with intracellular metabolomic analysis. The resulting genome-scale model will serve as a knowledge base of Nitrospira metabolism and lays the foundation for quantitative systems biology studies of these globally important nitrite-oxidizing bacteria.

Stepwise genetic engineering of Pseudomonas putida enables robust heterologous production of prodigiosin and glidobactin A.

Polyketide synthases (PKS) and nonribosomal peptide synthetases (NRPS) comprise biosynthetic pathways that provide access to diverse, often bioactive natural products. Metabolic engineering can improve production metrics to support characterization and drug-development studies, but often native hosts are difficult to genetically manipulate and/or culture. For this reason, heterologous expression is a common strategy for natural product discovery and characterization. Many bacteria have been developed to express heterologous biosynthetic gene clusters (BGCs) for producing polyketides and nonribosomal peptides. In this article, we describe tools for using Pseudomonas putida, a Gram-negative soil bacterium, as a heterologous host for producing natural products. Pseudomonads are known to produce many natural products, but P. putida production titers have been inconsistent in the literature and often low compared to other hosts. In recent years, synthetic biology tools for engineering P. putida have greatly improved, but their application towards production of natural products is limited. To demonstrate the potential of P. putida as a heterologous host, we introduced BGCs encoding the synthesis of prodigiosin and glidobactin A, two bioactive natural products synthesized from a combination of PKS and NRPS enzymology. Engineered strains exhibited robust production of both compounds after a single chromosomal integration of the corresponding BGC. Next, we took advantage of a set of genome-editing tools to increase titers by modifying transcription and translation of the BGCs and increasing the availability of auxiliary proteins required for PKS and NRPS activity. Lastly, we discovered genetic modifications to P. putida that affect natural product synthesis, including a strategy for removing a carbon sink that improves product titers. These efforts resulted in production strains capable of producing 1.1 g/L prodigiosin and 470 mg/L glidobactin A.

Autotrophic and mixotrophic metabolism of an anammox bacterium revealed by in vivo 13C and 2H metabolic network mapping.

Anaerobic ammonium-oxidizing (anammox) bacteria mediate a key step in the biogeochemical nitrogen cycle and have been applied worldwide for the energy-efficient removal of nitrogen from wastewater. However, outside their core energy metabolism, little is known about the metabolic networks driving anammox bacterial anabolism and use of different carbon and energy substrates beyond genome-based predictions. Here, we experimentally resolved the central carbon metabolism of the anammox bacterium Candidatus 'Kuenenia stuttgartiensis' using time-series 13C and 2H isotope tracing, metabolomics, and isotopically nonstationary metabolic flux analysis. Our findings confirm predicted metabolic pathways used for CO2 fixation, central metabolism, and amino acid biosynthesis in K. stuttgartiensis, and reveal several instances where genomic predictions are not supported by in vivo metabolic fluxes. This includes the use of the oxidative branch of an incomplete tricarboxylic acid cycle for alpha-ketoglutarate biosynthesis, despite the genome not having an annotated citrate synthase. We also demonstrate that K. stuttgartiensis is able to directly assimilate extracellular formate via the Wood-Ljungdahl pathway instead of oxidizing it completely to CO2 followed by reassimilation. In contrast, our data suggest that K. stuttgartiensis is not capable of using acetate as a carbon or energy source in situ and that acetate oxidation occurred via the metabolic activity of a low-abundance microorganism in the bioreactor's side population. Together, these findings provide a foundation for understanding the carbon metabolism of anammox bacteria at a systems-level and will inform future studies aimed at elucidating factors governing their function and niche differentiation in natural and engineered ecosystems.

In Vivo Thermodynamic Analysis of Glycolysis in Clostridium thermocellum and Thermoanaerobacterium saccharolyticum Using 13C and 2H Tracers.

Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are thermophilic anaerobic bacteria with complementary metabolic capabilities that utilize distinct glycolytic pathways for the conversion of cellulosic sugars to biofuels. We integrated quantitative metabolomics with 2H and 13C metabolic flux analysis to investigate the in vivo reversibility and thermodynamics of the central metabolic networks of these two microbes. We found that the glycolytic pathway in C. thermocellum operates remarkably close to thermodynamic equilibrium, with an overall drop in Gibbs free energy 5-fold lower than that of T. saccharolyticum or anaerobically grown Escherichia coli The limited thermodynamic driving force of glycolysis in C. thermocellum could be attributed in large part to the small free energy of the phosphofructokinase reaction producing fructose bisphosphate. The ethanol fermentation pathway was also substantially more reversible in C. thermocellum than in T. saccharolyticum These observations help explain the comparatively low ethanol titers of C. thermocellum and suggest engineering interventions that can be used to increase its ethanol productivity and glycolytic rate. In addition to thermodynamic analysis, we used our isotope tracer data to reconstruct the T. saccharolyticum central metabolic network, revealing exclusive use of the Embden-Meyerhof-Parnas (EMP) pathway for glycolysis, a bifurcated tricarboxylic acid (TCA) cycle, and a sedoheptulose bisphosphate bypass active within the pentose phosphate pathway.IMPORTANCE Thermodynamics constitutes a key determinant of flux and enzyme efficiency in metabolic networks. Here, we provide new insights into the divergent thermodynamics of the glycolytic pathways of C. thermocellum and T. saccharolyticum, two industrially relevant thermophilic bacteria whose metabolism still is not well understood. We report that while the glycolytic pathway in T. saccharolyticum is as thermodynamically favorable as that found in model organisms, such as E. coli or Saccharomyces cerevisiae, the glycolytic pathway of C. thermocellum operates near equilibrium. The use of a near-equilibrium glycolytic pathway, with potentially increased ATP yield, by this cellulolytic microbe may represent an evolutionary adaptation to growth on cellulose, but it has the drawback of being highly susceptible to product feedback inhibition. The results of this study will facilitate future engineering of high-performance strains capable of transforming cellulosic biomass to biofuels at high yields and titers.

Expression of Phosphofructokinase Is Not Sufficient to Enable Embden-Meyerhof-Parnas Glycolysis in Zymomonas mobilis ZM4.

Zymomonas mobilis is a bacterium that produces ethanol from glucose at up to 97% of theoretical efficiency on a carbon basis. One factor contributing to the high efficiency of ethanol production is that Z. mobilis has a low biomass yield. The low biomass yield may be caused partly by the low ATP yield of the Entner-Doudoroff (ED) glycolytic pathway used by Z. mobilis, which produces only one ATP per glucose consumed. To test the hypothesis that ATP yield limits biomass yield in Z. mobilis, we attempted to introduce the Embden-Meyerhof-Parnas (EMP) glycolytic pathway (with double the ATP yield) by expressing phosphofructokinase (Pfk I) from Escherichia coli. Expression of Pfk I caused growth inhibition and resulted in accumulation of mutations in the pfkA gene. Co-expression of additional EMP enzymes, fructose bisphosphate aldolase (Fba) and triose phosphate isomerase (Tpi), with Pfk I did not enable EMP flux, and resulted in production of glycerol as a side product. Further analysis indicated that heterologous reactions may have operated in the reverse direction because of native metabolite concentrations. This study reveals how the metabolomic context of a chassis organism influences the range of pathways that can be added by heterologous expression.

Near-equilibrium glycolysis supports metabolic homeostasis and energy yield.

Glycolysis plays a central role in producing ATP and biomass. Its control principles, however, remain incompletely understood. Here, we develop a method that combines 2H and 13C tracers to determine glycolytic thermodynamics. Using this method, we show that, in conditions and organisms with relatively slow fluxes, multiple steps in glycolysis are near to equilibrium, reflecting spare enzyme capacity. In Escherichia coli, nitrogen or phosphorus upshift rapidly increases the thermodynamic driving force, deploying the spare enzyme capacity to increase flux. Similarly, respiration inhibition in mammalian cells rapidly increases both glycolytic flux and the thermodynamic driving force. The thermodynamic shift allows flux to increase with only small metabolite concentration changes. Finally, we find that the cellulose-degrading anaerobe Clostridium cellulolyticum exhibits slow, near-equilibrium glycolysis due to the use of pyrophosphate rather than ATP for fructose-bisphosphate production, resulting in enhanced per-glucose ATP yield. Thus, near-equilibrium steps of glycolysis promote both rapid flux adaptation and energy efficiency.

2H and 13C metabolic flux analysis elucidates in vivo thermodynamics of the ED pathway in Zymomonas mobilis.

Zymomonas mobilis is an industrially relevant bacterium notable for its ability to rapidly ferment simple sugars to ethanol using the Entner-Doudoroff (ED) glycolytic pathway, an alternative to the well-known Embden-Meyerhof-Parnas (EMP) pathway used by most organisms. Recent computational studies have predicted that the ED pathway is substantially more thermodynamically favorable than the EMP pathway, a potential factor explaining the high glycolytic rate in Z. mobilis. Here, to investigate the in vivo thermodynamics of the ED pathway and central carbon metabolism in Z. mobilis, we implemented a network-level approach that integrates quantitative metabolomics with 2H and 13C metabolic flux analysis to estimate reversibility and Gibbs free energy (ΔG) of metabolic reactions. This analysis revealed a strongly thermodynamically favorable ED pathway in Z. mobilis that is nearly twice as favorable as the EMP pathway in E. coli or S. cerevisiae. The in vivo step-by-step thermodynamic profile of the ED pathway was highly similar to previous in silico predictions, indicating that maximizing ΔG for each pathway step likely constitutes a cellular objective in Z. mobilis. Our analysis also revealed novel features of Z. mobilis metabolism, including phosphofructokinase-like enzyme activity, tricarboxylic acid cycle anaplerosis via PEP carboxylase, and a metabolic imbalance in the pentose phosphate pathway resulting in excretion of shikimate pathway intermediates. The integrated approach we present here for in vivo ΔG quantitation may be applied to the thermodynamic profiling of pathways and metabolic networks in other microorganisms and will contribute to the development of quantitative models of metabolism.

Phenolic Amides Are Potent Inhibitors of De Novo Nucleotide Biosynthesis.

An outstanding challenge toward efficient production of biofuels and value-added chemicals from plant biomass is the impact that lignocellulose-derived inhibitors have on microbial fermentations. Elucidating the mechanisms that underlie their toxicity is critical for developing strategies to overcome them. Here, using Escherichia coli as a model system, we investigated the metabolic effects and toxicity mechanisms of feruloyl amide and coumaroyl amide, the predominant phenolic compounds in ammonia-pretreated biomass hydrolysates. Using metabolomics, isotope tracers, and biochemical assays, we showed that these two phenolic amides act as potent and fast-acting inhibitors of purine and pyrimidine biosynthetic pathways. Feruloyl or coumaroyl amide exposure leads to (i) a rapid buildup of 5-phosphoribosyl-1-pyrophosphate (PRPP), a key precursor in nucleotide biosynthesis, (ii) a rapid decrease in the levels of pyrimidine biosynthetic intermediates, and (iii) a long-term generalized decrease in nucleotide and deoxynucleotide levels. Tracer experiments using (13)C-labeled sugars and [(15)N]ammonia demonstrated that carbon and nitrogen fluxes into nucleotides and deoxynucleotides are inhibited by these phenolic amides. We found that these effects are mediated via direct inhibition of glutamine amidotransferases that participate in nucleotide biosynthetic pathways. In particular, feruloyl amide is a competitive inhibitor of glutamine PRPP amidotransferase (PurF), which catalyzes the first committed step in de novo purine biosynthesis. Finally, external nucleoside supplementation prevents phenolic amide-mediated growth inhibition by allowing nucleotide biosynthesis via salvage pathways. The results presented here will help in the development of strategies to overcome toxicity of phenolic compounds and facilitate engineering of more efficient microbial producers of biofuels and chemicals.