¹³C-based metabolic flux analysis of Saccharomyces cerevisiae with a reduced Crabtree effect

13C-metabolic flux analysis of Saccharomyces cerevisiae in complex media

Saccharomyces cerevisiae is often cultivated in complex media for applications in food and other biochemical production. However, 13C-metabolic flux analysis (13C-MFA) has been conducted for S. cerevisiae cultivated in synthetic media, resulting in a limited understanding of the metabolic flux distributions under the complex media. In this study, 13C-MFA was applied to S. cerevisiae cultivated in complex media to quantify the metabolic fluxes in the central metabolic network. S. cerevisiae was cultivated in a synthetic dextrose (SD) medium supplemented with 20 amino acids (SD + AA) and yeast extract peptone dextrose (YPD) medium. The results revealed that glutamic acid, glutamine, aspartic acid, and asparagine are incorporated into the TCA cycle as carbon sources in parallel with glucose consumption. Based on these findings, we successfully conducted 13C-MFA of S. cerevisiae cultivated in SD + AA and YPD media using parallel labeling and measured amino acid uptake rates. Furthermore, we applied the developed approach to 13C-MFA of yeast cultivated in malt extract medium. The analysis revealed that the metabolic flux through the anaplerotic and oxidative pentose phosphate pathways was lower in complex media than in synthetic media. Owing to the reduced carbon loss by the branching pathways, carbon flow toward ethanol production via glycolysis could be elevated. 13C-MFA of S. cerevisiae cultured in complex media provides valuable insights for metabolic engineering and process optimization in industrial yeast fermentation.

Comparative study of two Saccharomyces cerevisiae strains with kinetic models at genome-scale

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by ¹³C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

Glucose assimilation rate determines the partition of flux at pyruvate between lactic acid and ethanol in Saccharomyces cerevisiae

Engineered Saccharomyces cerevisiae expressing a lactic acid dehydrogenase can metabolize pyruvate into lactic acid. However, three pyruvate decarboxylase (PDC) isozymes drive most carbon flux toward ethanol rather than lactic acid. Deletion of endogenous PDCs will eliminate ethanol production, but the resulting strain suffers from C₂ auxotrophy and struggles to complete a fermentation. Engineered yeast assimilating xylose or cellobiose produce lactic acid rather than ethanol as a major product without the deletion of any PDC genes. We report here that sugar flux, but not sensing, contributes to the partition of flux at the pyruvate branch point in S. cerevisiae expressing the Rhizopus oryzae lactic acid dehydrogenase (LdhA). While the membrane glucose sensors Snf3 and Rgt2 did not play any direct role in the option of predominant product, the sugar assimilation rate was strongly correlated to the partition of flux at pyruvate: fast sugar assimilation favors ethanol production while slow sugar assimilation favors lactic acid. Applying this knowledge, we created an engineered yeast capable of simultaneously converting glucose and xylose into lactic acid, increasing lactic acid production to approximately 17 g L^-1 from the 12 g L^-1 observed during sequential consumption of sugars. This work elucidates the carbon source-dependent effects on product selection in engineered yeast.

Comparative 13 C-metabolic flux analysis indicates elevation of ATP regeneration, carbon dioxide, and heat production in industrial Saccharomyces cerevisiae strains

BACKGROUND

Various industrial Saccharomyces cerevisiae strains are used for specific processes, such as sake, wine brewing and bread making. Understanding mechanisms underlying the fermentation performance of these strains would be useful for further engineering of the S. cerevisiae metabolism. However, the relationship between the fermentation performance, intra-cellular metabolic states, and other phenotypic characteristics of industrial yeasts is still unclear. In this study, ¹³ C-metabolic flux analysis of four diploid yeast strains-laboratory, sake, bread, and wine yeasts-was conducted.

RESULTS

While the Crabtree effect was observed for all strains, the metabolic flux level of glycolysis was elevated in bread and sake yeast. Furthermore, increased flux levels of the TCA cycle were commonly observed in the three industrial strains. The specific rates of CO₂ production, net ATP regeneration, and metabolic heat generation estimated from the metabolic flux distribution were two to three times greater than those of the laboratory strain. The elevation in metabolic heat generation was correlated with the tolerance to low-temperature stress.

CONCLUSION

These results indicate that the metabolic flux distribution of sake and bread yeast strains contributes to faster production of ethanol and CO₂ . It is also suggested that the generation of metabolic heat is preferable under the actual industrial fermentation conditions.

13C-metabolic flux analysis in glycerol-assimilating strains of Saccharomyces cerevisiae

Glycerol is an attractive raw material for the production of useful chemicals using microbial cells. We previously identified metabolic engineering targets for the improvement of glycerol assimilation ability in Saccharomyces cerevisiae based on adaptive laboratory evolution (ALE) and transcriptome analysis of the evolved cells. We also successfully improved glycerol assimilation ability by the disruption of the RIM15 gene encoding a Greatwall protein kinase together with overexpression of the STL1 gene encoding the glycerol/H⁺ symporter. To understand glycerol assimilation metabolism in the evolved glycerol-assimilating strains and STL1-overexpressing RIM15 disruptant, we performed metabolic flux analysis using ¹³C-labeled glycerol. Significant differences in metabolic flux distributions between the strains obtained from the culture after 35 and 85 generations in ALE were not found, indicating that metabolic flux changes might occur in the early phase of ALE (i.e., before 35 generations at least). Similarly, metabolic flux distribution was not significantly changed by RIM15 gene disruption. However, fluxes for the lower part of glycolysis and the TCA cycle were larger and, as a result, flux for the pentose phosphate pathway was smaller in the STL1-overexpressing RIM15 disruptant than in the strain obtained from the culture after 85 generations in ALE. It could be effective to increase flux for the pentose phosphate pathway to improve the glycerol assimilation ability in S. cerevisiae.

Repression of mitochondrial metabolism for cytosolic pyruvate-derived chemical production in Saccharomyces cerevisiae

Background

Saccharomyces cerevisiae is a suitable host for the industrial production of pyruvate-derived chemicals such as ethanol and 2,3-butanediol (23BD). For the improvement of the productivity of these chemicals, it is essential to suppress the unnecessary pyruvate consumption in S. cerevisiae to redirect the metabolic flux toward the target chemical production. In this study, mitochondrial pyruvate transporter gene (MPC1) or the essential gene for mitophagy (ATG32) was knocked-out to repress the mitochondrial metabolism and improve the production of pyruvate-derived chemical in S. cerevisiae.

Results

The growth rates of both aforementioned strains were 1.6-fold higher than that of the control strain. ¹³C-metabolic flux analysis revealed that both strains presented similar flux distributions and successfully decreased the tricarboxylic acid cycle fluxes by 50% compared to the control strain. Nevertheless, the intracellular metabolite pool sizes were completely different, suggesting distinct metabolic effects of gene knockouts in both strains. This difference was also observed in the test-tube culture for 23BD production. Knockout of ATG32 revealed a 23.6-fold increase in 23BD titer (557.0 ± 20.6 mg/L) compared to the control strain (23.5 ± 12.8 mg/L), whereas the knockout of MPC1 revealed only 14.3-fold increase (336.4 ± 113.5 mg/L). Further investigation using the anaerobic high-density fermentation test revealed that the MPC1 knockout was more effective for ethanol production than the 23BD production.

Conclusion

These results suggest that the engineering of the mitochondrial transporters and membrane dynamics were effective in controlling the mitochondrial metabolism to improve the productivities of chemicals in yeast cytosol.

Genome-Scale 13C Fluxomics Modeling for Metabolic Engineering of Saccharomyces cerevisiae

Synthetic biology is a rapidly developing field that pursues the application of engineering principles and development approaches to biological engineering. Synthetic biology is poised to change the way biology is practiced, and has important practical applications: for example, building genetically engineered organisms to produce biofuels, medicines, and other chemicals. Traditionally, synthetic biology has focused on manipulating a few genes (e.g., in a single pathway or genetic circuit), but its combination with systems biology holds the promise of creating new cellular architectures and constructing complex biological systems from the ground up. Enabling this merge of synthetic and systems biology will require greater predictive capability for modeling the behavior of cellular systems, and more comprehensive data sets for building and calibrating these models. The so-called "-omics" data sets can now be generated via high throughput techniques in the form of genomic, proteomic, transcriptomic, and metabolomic information on the engineered biological system. Of particular interest with respect to the engineering of microbes capable of producing biofuels and other chemicals economically and at scale are metabolomic datasets, and their insights into intracellular metabolic fluxes. Metabolic fluxes provide a rapid and easy to understand picture of how carbon and energy flow throughout the cell. Here, we present a detailed guide to performing metabolic flux analysis and modeling using the open source JBEI Quantitative Metabolic Modeling (jQMM) library. This library allows the user to transform metabolomics data in the form of isotope labeling data from a ¹³C labeling experiment into a determination of cellular fluxes that can be used to develop genetic engineering strategies for metabolic engineering.The jQMM library presents a complete toolbox for performing a range of different tasks of interest in metabolic engineering. Various different types of flux analysis and modeling can be performed such as flux balance analysis, ¹³C metabolic flux analysis, and two-scale ¹³C metabolic flux analysis (2S-¹³C MFA). 2S-¹³C MFA is a novel method that determines genome-scale fluxes without the need of every single carbon transition in the metabolic network. In addition to several other capabilities, the jQMM library can make model based predictions for how various genetic engineering strategies can be incorporated toward bioengineering goals: it can predict the effects of reaction knockouts on metabolism using both the MoMA and ROOM methodologies. In this chapter, we will illustrate the use of the jQMM library through a step-by-step demonstration of flux determination and knockout prediction in a complex eukaryotic model organism: Saccharomyces cerevisiae (S. cerevisiae). Included with this chapter is a digital Jupyter Notebook file that provides a computable appendix showing a self-contained example of jQMM usage, which can be changed to fit the user's specific needs. As an open source software project, users can modify and extend the code base to make improvements at will, allowing them to share their development work and contribute back to the jQMM modeling community.

Constraining Genome-Scale Models to Represent the Bow Tie Structure of Metabolism for 13C Metabolic Flux Analysis

Determination of internal metabolic fluxes is crucial for fundamental and applied biology because they map how carbon and electrons flow through metabolism to enable cell function. 13C Metabolic Flux Analysis (13C MFA) and Two-Scale 13C Metabolic Flux Analysis (2S-13C MFA) are two techniques used to determine such fluxes. Both operate on the simplifying approximation that metabolic flux from peripheral metabolism into central “core” carbon metabolism is minimal, and can be omitted when modeling isotopic labeling in core metabolism. The validity of this “two-scale” or “bow tie” approximation is supported both by the ability to accurately model experimental isotopic labeling data, and by experimentally verified metabolic engineering predictions using these methods. However, the boundaries of core metabolism that satisfy this approximation can vary across species, and across cell culture conditions. Here, we present a set of algorithms that (1) systematically calculate flux bounds for any specified “core” of a genome-scale model so as to satisfy the bow tie approximation and (2) automatically identify an updated set of core reactions that can satisfy this approximation more efficiently. First, we leverage linear programming to simultaneously identify the lowest fluxes from peripheral metabolism into core metabolism compatible with the observed growth rate and extracellular metabolite exchange fluxes. Second, we use Simulated Annealing to identify an updated set of core reactions that allow for a minimum of fluxes into core metabolism to satisfy these experimental constraints. Together, these methods accelerate and automate the identification of a biologically reasonable set of core reactions for use with 13C MFA or 2S-13C MFA, as well as provide for a substantially lower set of flux bounds for fluxes into the core as compared with previous methods. We provide an open source Python implementation of these algorithms at https://github.com/JBEI/limitfluxtocore.

A pyruvate carbon flux tugging strategy for increasing 2,3-butanediol production and reducing ethanol subgeneration in the yeast Saccharomyces cerevisiae

BACKGROUND

The yeast Saccharomyces cerevisiae is a promising host cell for producing a wide range of chemicals. However, attempts to metabolically engineer Crabtree-positive S. cerevisiae invariably face a common issue: how to reduce dominant ethanol production. Here, we propose a yeast metabolic engineering strategy for decreasing ethanol subgeneration involving tugging the carbon flux at an important hub branching point (e.g., pyruvate). Tugging flux at a central glycolytic overflow metabolism point arising from high glycolytic activity may substantially increase higher alcohol production in S. cerevisiae. We validated this possibility by testing 2,3-butanediol (2,3-BDO) production, which is routed via pyruvate as the important hub compound.

RESULTS

By searching for high-activity acetolactate synthase (ALS) enzymes that catalyze the important first-step reaction in 2,3-BDO biosynthesis, and tuning several fermentation conditions, we demonstrated that a stronger pyruvate pulling effect (tugging of pyruvate carbon flux) is very effective for increasing 2,3-BDO production and reducing ethanol subgeneration by S. cerevisiae. To further confirm the validity of the pyruvate carbon flux tugging strategy, we constructed an evolved pyruvate decarboxylase (PDC)-deficient yeast (PDCΔ) strain that lacked three isozymes of PDC. In parallel with re-sequencing to identify genomic mutations, liquid chromatography-tandem mass spectrometry analysis of intermediate metabolites revealed significant accumulation of pyruvate and NADH in the evolved PDCΔ strain. Harnessing the high-activity ALS and additional downstream enzymes in the evolved PDCΔ strain resulted in a high yield of 2,3-BDO (a maximum of 0.41 g g^-1 glucose consumed) and no ethanol subgeneration, thereby confirming the utility of our strategy. Using this engineered strain, we demonstrated a high 2,3-BDO titer (81.0 g L^-1) in a fed-batch fermentation using a high concentration of glucose as the sole carbon source.

CONCLUSIONS

We demonstrated that the pyruvate carbon flux tugging strategy is very effective for increasing 2,3-BDO production and decreasing ethanol subgeneration in Crabtree-positive S. cerevisiae. High activity of the common first-step enzyme for the conversion of pyruvate, which links to both the TCA cycle and amino acid biosynthesis, is likely important for the production of various chemicals by S. cerevisiae.

Rerouting of carbon flux in a glycogen mutant of cyanobacteria assessed via isotopically non-stationary 13 C metabolic flux analysis

Cyanobacteria, which constitute a quantitatively dominant phylum, have attracted attention in biofuel applications due to favorable physiological characteristics, high photosynthetic efficiency and amenability to genetic manipulations. However, quantitative aspects of cyanobacterial metabolism have received limited attention. In the present study, we have performed isotopically non-stationary ¹³ C metabolic flux analysis (INST-¹³ C-MFA) to analyze rerouting of carbon in a glycogen synthase deficient mutant strain (glgA-I glgA-II) of the model cyanobacterium Synechococcus sp. PCC 7002. During balanced photoautotrophic growth, 10-20% of the fixed carbon is stored in the form of glycogen via a pathway that is conserved across the cyanobacterial phylum. Our results show that deletion of glycogen synthase gene orchestrates cascading effects on carbon distribution in various parts of the metabolic network. Carbon that was originally destined to be incorporated into glycogen gets partially diverted toward alternate storage molecules such as glucosylglycerol and sucrose. The rest is partitioned within the metabolic network, primarily via glycolysis and tricarboxylic acid cycle. A lowered flux toward carbohydrate synthesis and an altered distribution at the glucose-1-phosphate node indicate flexibility in the network. Further, reversibility of glycogen biosynthesis reactions points toward the presence of futile cycles. Similar redistribution of carbon was also predicted by Flux Balance Analysis. The results are significant to metabolic engineering efforts with cyanobacteria where fixed carbon needs to be re-routed to products of interest. Biotechnol. Bioeng. 2017;114: 2298-2308. © 2017 Wiley Periodicals, Inc.

13C Metabolic Flux Analysis for Systematic Metabolic Engineering of S. cerevisiae for Overproduction of Fatty Acids

Efficient redirection of microbial metabolism into the abundant production of desired bioproducts remains non-trivial. Here, we used flux-based modeling approaches to improve yields of fatty acids in Saccharomyces cerevisiae. We combined ¹³C labeling data with comprehensive genome-scale models to shed light onto microbial metabolism and improve metabolic engineering efforts. We concentrated on studying the balance of acetyl-CoA, a precursor metabolite for the biosynthesis of fatty acids. A genome-wide acetyl-CoA balance study showed ATP citrate lyase from Yarrowia lipolytica as a robust source of cytoplasmic acetyl-CoA and malate synthase as a desirable target for downregulation in terms of acetyl-CoA consumption. These genetic modifications were applied to S. cerevisiae WRY2, a strain that is capable of producing 460 mg/L of free fatty acids. With the addition of ATP citrate lyase and downregulation of malate synthase, the engineered strain produced 26% more free fatty acids. Further increases in free fatty acid production of 33% were obtained by knocking out the cytoplasmic glycerol-3-phosphate dehydrogenase, which flux analysis had shown was competing for carbon flux upstream with the carbon flux through the acetyl-CoA production pathway in the cytoplasm. In total, the genetic interventions applied in this work increased fatty acid production by ~70%.

Achieving Metabolic Flux Analysis for S. cerevisiae at a Genome-Scale: Challenges, Requirements, and Considerations

Recent advances in 13C-Metabolic flux analysis (13C-MFA) have increased its capability to accurately resolve fluxes using a genome-scale model with narrow confidence intervals without pre-judging the activity or inactivity of alternate metabolic pathways. However, the necessary precautions, computational challenges, and minimum data requirements for successful analysis remain poorly established. This review aims to establish the necessary guidelines for performing 13C-MFA at the genome-scale for a compartmentalized eukaryotic system such as yeast in terms of model and data requirements, while addressing key issues such as statistical analysis and network complexity. We describe the various approaches used to simplify the genome-scale model in the absence of sufficient experimental flux measurements, the availability and generation of reaction atom mapping information, and the experimental flux and metabolite labeling distribution measurements to ensure statistical validity of the obtained flux distribution. Organism-specific challenges such as the impact of compartmentalization of metabolism, variability of biomass composition, and the cell-cycle dependence of metabolism are discussed. Identification of errors arising from incorrect gene annotation and suggested alternate routes using MFA are also highlighted.

A Method to Constrain Genome-Scale Models with 13C Labeling Data

Current limitations in quantitatively predicting biological behavior hinder our efforts to engineer biological systems to produce biofuels and other desired chemicals. Here, we present a new method for calculating metabolic fluxes, key targets in metabolic engineering, that incorporates data from ¹³C labeling experiments and genome-scale models. The data from ¹³C labeling experiments provide strong flux constraints that eliminate the need to assume an evolutionary optimization principle such as the growth rate optimization assumption used in Flux Balance Analysis (FBA). This effective constraining is achieved by making the simple but biologically relevant assumption that flux flows from core to peripheral metabolism and does not flow back. The new method is significantly more robust than FBA with respect to errors in genome-scale model reconstruction. Furthermore, it can provide a comprehensive picture of metabolite balancing and predictions for unmeasured extracellular fluxes as constrained by ¹³C labeling data. A comparison shows that the results of this new method are similar to those found through ¹³C Metabolic Flux Analysis (¹³C MFA) for central carbon metabolism but, additionally, it provides flux estimates for peripheral metabolism. The extra validation gained by matching 48 relative labeling measurements is used to identify where and why several existing COnstraint Based Reconstruction and Analysis (COBRA) flux prediction algorithms fail. We demonstrate how to use this knowledge to refine these methods and improve their predictive capabilities. This method provides a reliable base upon which to improve the design of biological systems.

While metabolic fluxes constitute the most direct window into a cell’s metabolism, their accurate measurement is non trivial. The gold standard for flux measurement involves providing a labeled feed where some of the carbon atoms have been substituted by isotopes with higher atomic mass (¹³C instead of ¹²C). The ensuing labeling found in intracellular metabolites is then used to computationally infer the metabolic fluxes that produced the observed pattern. However, this procedure is typically performed with small metabolic models encompassing only central carbon metabolism. The genomic revolution has afforded us easily available genomes and, with them, comprehensive genome-scale models of cellular metabolism. It would be desirable to use the ¹³C labeling experimental data to constrain genome-scale models: these data constrain fluxes very effectively and provide in the labeling data fit an obvious proof that the underlying model correctly explains measured quantities. Here, we introduce a rigorous, self-consistent method that uses the full amount of information contained in ¹³C labeling data to constrain fluxes for a genome-scale model where underlying assumptions are explicitly stated.