Physiological responses of Saccharomyces cerevisiae to industrially relevant conditions: Slow growth, low pH, and high CO2 levels

High throughput 13C-metabolic flux analysis of 3-hydroxypropionic acid producing Pichia pastoris reveals limited availability of acetyl-CoA and ATP due to tight control of the glycolytic flux

BACKGROUND

Production of 3-hydroxypropionic acid (3-HP) through the malonyl-CoA pathway has yielded promising results in Pichia pastoris (Komagataella phaffii), demonstrating the potential of this cell factory to produce this platform chemical and other acetyl-CoA-derived products using glycerol as a carbon source. However, further metabolic engineering of the original P. pastoris 3-HP-producing strains resulted in unexpected outcomes, e.g., significantly lower product yield and/or growth rate. To gain an understanding on the metabolic constraints underlying these observations, the fluxome (metabolic flux phenotype) of ten 3-HP-producing P. pastoris strains has been characterized using a high throughput ¹³C-metabolic flux analysis platform. Such platform enabled the operation of an optimised workflow to obtain comprehensive maps of the carbon flux distribution in the central carbon metabolism in a parallel-automated manner, thereby accelerating the time-consuming strain characterization step in the design-build-test-learn cycle for metabolic engineering of P. pastoris.

RESULTS

We generated detailed maps of the carbon fluxes in the central carbon metabolism of the 3-HP producing strain series, revealing the metabolic consequences of different metabolic engineering strategies aimed at improving NADPH regeneration, enhancing conversion of pyruvate into cytosolic acetyl-CoA, or eliminating by-product (arabitol) formation. Results indicate that the expression of the POS5 NADH kinase leads to a reduction in the fluxes of the pentose phosphate pathway reactions, whereas an increase in the pentose phosphate pathway fluxes was observed when the cytosolic acetyl-CoA synthesis pathway was overexpressed. Results also show that the tight control of the glycolytic flux hampers cell growth due to limited acetyl-CoA biosynthesis. When the cytosolic acetyl-CoA synthesis pathway was overexpressed, the cell growth increased, but the product yield decreased due to higher growth-associated ATP costs. Finally, the six most relevant strains were also cultured at pH 3.5 to assess the effect of a lower pH on their fluxome. Notably, similar metabolic fluxes were observed at pH 3.5 compared to the reference condition at pH 5.

CONCLUSIONS

This study shows that existing fluoxomics workflows for high-throughput analyses of metabolic phenotypes can be adapted to investigate P. pastoris, providing valuable information on the impact of genetic manipulations on the metabolic phenotype of this yeast. Specifically, our results highlight the metabolic robustness of P. pastoris's central carbon metabolism when genetic modifications are made to increase the availability of NADPH and cytosolic acetyl-CoA. Such knowledge can guide further metabolic engineering of these strains. Moreover, insights into the metabolic adaptation of P. pastoris to an acidic pH have also been obtained, showing the capability of the fluoxomics workflow to assess the metabolic impact of environmental changes.

Challenges and progress towards industrial recombinant protein production in yeasts: A review

Recombinant proteins (RP) are widely used as biopharmaceuticals, industrial enzymes, or sustainable food source. Yeasts, with their ability to produce complex proteins through a broad variety of cheap carbon sources, have emerged as promising eukaryotic production hosts. As such, the prevalence of yeasts as favourable production organisms in commercial RP production is expected to increase. Yet, with the selection of a robust production host on the one hand, successful scale-up is dependent on a thorough understanding of the challenging environment and limitations of large-scale bioreactors on the other hand. In the present work, several prominent yeast species, including Saccharomyces cerevisiae, Pichia pastoris, Yarrowia lipolytica, Kluyveromyces lactis and Kluyveromyces marxianus are reviewed for their current state and performance in commercial RP production. Thereafter, the impact of principal process control parameters, including dissolved oxygen, pH, substrate concentration, and temperature, on large-scale RP production are discussed. Finally, technical challenges of process scale-up are identified. To that end, process intensification strategies to enhance industrial feasibility are summarized, specifically highlighting fermentation strategies to ensure sufficient cooling capacity, overcome oxygen limitation, and increase protein quality and productivity. As such, this review aims to contribute to the pursuit of sustainable yeast-based RP production.

Genome-wide transcriptional regulation in Saccharomyces cerevisiae in response to carbon dioxide

Sugar metabolism by Saccharomyces cerevisiae produces ample amounts of CO2 under both aerobic and anaerobic conditions. High solubility of CO2 in fermentation media, contributing to enjoyable sensory properties of sparkling wine and beers by S. cerevisiae, might affect yeast metabolism. To elucidate the overlooked effects of CO2 on yeast metabolism, we examined glucose fermentation by S. cerevisiae under CO2 as compared to N2 and O2 limited conditions. While both CO2 and N2 conditions are considered anaerobic, less glycerol and acetate but more ethanol were produced under CO2 condition. Transcriptomic analysis revealed that significantly decreased mRNA levels of GPP1 coding for glycerol-3-phosphate phosphatase in glycerol synthesis explained the reduced glycerol production under CO2 condition. Besides, transcriptional regulations in signal transduction, carbohydrate synthesis, heme synthesis, membrane and cell wall metabolism, and respiration were detected in response to CO2. Interestingly, signal transduction was uniquely regulated under CO2 condition, where up-regulated genes (STE3, MSB2, WSC3, STE12 and TEC1) in the signal sensors and transcriptional factors suggested that MAPK signaling pathway plays a critical role in CO2 sensing and CO2-induced metabolisms in yeast. Our study identifies CO2 as an external stimulus for modulating metabolic activities in yeast and a transcriptional effector for diverse applications.

Robustness: linking strain design to viable bioprocesses

Microbial cell factories are becoming increasingly popular for the sustainable production of various chemicals. Metabolic engineering has led to the design of advanced cell factories; however, their long-term yield, titer, and productivity falter when scaled up and subjected to industrial conditions. This limitation arises from a lack of robustness - the ability to maintain a constant phenotype despite the perturbations of such processes. This review describes predictable and stochastic industrial perturbations as well as state-of-the-art technologies to counter process variability. Moreover, we distinguish robustness from tolerance and discuss the potential of single-cell studies for improving system robustness. Finally, we highlight ways of achieving consistent and comparable quantification of robustness that can guide the selection of strains for industrial bioprocesses.

Recombinant protein linker production as a basis for non-invasive determination of single-cell yeast age in heterogeneous yeast populations

The physiological and metabolic diversity of a yeast culture is the sum of individual cell phenotypes. As well as environmental conditions, genetics, and numbers of cell divisions, a major factor influencing cell characteristics is cell age. A postcytokinesis bud scar on the mother cell, a benchmark in the replicative life span, is a quantifiable indicator of cell age, characterized by significant amounts of chitin. We developed a binding process for visualizing the bud scars of Saccharomyces pastorianus var. carlsbergensis using a protein linker containing a polyhistidine tag, a superfolder green fluorescent protein (sfGFP), and a chitin-binding domain (His6-SUMO-sfGFP-ChBD). The binding did not affect yeast viability; thus, our method provides the basis for non-invasive cell age determination using flow cytometry. The His6-SUMO-sfGFP-ChBD protein was synthesized in Escherichia coli, purified using two-stage chromatography, and checked for monodispersity and purity. Linker-cell binding and the characteristics of the bound complex were determined using flow cytometry and confocal laser scanning microscopy (CLSM). Flow cytometry showed that protein binding increased to 60 455 ± 2706 fluorescence units per cell. The specific coupling of the linker to yeast cells was additionally verified by CLSM and adsorption isotherms using yeast cells, E. coli cells, and chitin resin. We found a relationship between the median bud scar number, the median of the fluorescence units, and the chitin content of yeast cells. A fast measurement of yeast population dynamics by flow cytometry is possible, using this protein binding technique. Rapid qualitative determination of yeast cell age distribution can therefore be performed.

Uncoupling growth and succinic acid production in an industrial Saccharomyces cerevisiae strain

This study explores the relation between biomass‐specific succinic acid (SA) production rate and specific growth rate of an engineered industrial strain of Saccharomyces cerevisiae, with the aim to investigate the extent to which growth and product formation can be uncoupled. Ammonium‐limited aerobic chemostat and retentostat cultures were grown at different specific growth rates under industrially relevant conditions, that is, at a culture pH of 3 and with sparging of a 1:1 CO₂–air mixture. Biomass‐specific SA production rates decreased asymptotically with decreasing growth rate. At near‐zero growth rates, the engineered strain maintained a stable biomass‐specific SA production rate for over 500 h, with a SA yield on glucose of 0.61 mol mol⁻¹. These results demonstrate that uncoupling of growth and SA production could indeed be achieved. A linear relation between the biomass‐specific SA production rate and glucose consumption rate indicated the coupling of SA production rate and the flux through primary metabolism. The low culture pH resulted in an increased death rate, which was lowest at near‐zero growth rates. Nevertheless, a significant amount of non‐viable biomass accumulated in the retentostat cultures, thus underlining the importance of improving low‐pH tolerance in further strain development for industrial SA production with S. cerevisiae.

Impact of Altered Trehalose Metabolism on Physiological Response of Penicillium chrysogenum Chemostat Cultures during Industrially Relevant Rapid Feast/Famine Conditions

Due to insufficient mass transfer and mixing issues, cells in the industrial-scale bioreactor are often forced to experience glucose feast/famine cycles, mostly resulting in reduced commercial metrics (titer, yield and productivity). Trehalose cycling has been confirmed as a double-edged sword in the Penicillium chrysogenum strain, which facilitates the maintenance of a metabolically balanced state, but it consumes extra amounts of the ATP responsible for the repeated breakdown and formation of trehalose molecules in response to extracellular glucose perturbations. This loss of ATP would be in competition with the high ATP-demanding penicillin biosynthesis. In this work, the role of trehalose metabolism was further explored under industrially relevant conditions by cultivating a high-yielding Penicillium chrysogenum strain, and the derived trehalose-null strains in the glucose-limited chemostat system where the glucose feast/famine condition was imposed. This dynamic feast/famine regime with a block-wise feed/no feed regime (36 s on, 324 s off) allows one to generate repetitive cycles of moderate changes in glucose availability. The results obtained using quantitative metabolomics and stoichiometric analysis revealed that the intact trehalose metabolism is vitally important for maintaining penicillin production capacity in the Penicillium chrysogenum strain under both steady state and dynamic conditions. Additionally, cells lacking such a key metabolic regulator would become more sensitive to industrially relevant conditions, and are more able to sustain metabolic rearrangements, which manifests in the shrinkage of the central metabolite pool size and the formation of ATP-consuming futile cycles.

Predicting By-Product Gradients of Baker’s Yeast Production at Industrial Scale: A Practical Simulation Approach

Scaling up bioprocesses is one of the most crucial steps in the commercialization of bioproducts. While it is known that concentration and shear rate gradients occur at larger scales, it is often too risky, if feasible at all, to conduct validation experiments at such scales. Using computational fluid dynamics equipped with mechanistic biochemical engineering knowledge of the process, it is possible to simulate such gradients. In this work, concentration profiles for the by-products of baker’s yeast production are investigated. By applying a mechanistic black-box model, concentration heterogeneities for oxygen, glucose, ethanol, and carbon dioxide are evaluated. The results suggest that, although at low concentrations, ethanol is consumed in more than 90% of the tank volume, which prevents cell starvation, even when glucose is virtually depleted. Moreover, long exposure to high dissolved carbon dioxide levels is predicted. Two biomass concentrations, i.e., 10 and 25 g/L, are considered where, in the former, ethanol production is solely because of overflow metabolism while, in the latter, 10% of the ethanol formation is due to dissolved oxygen limitation. This method facilitates the prediction of the living conditions of the microorganism and its utilization to address the limitations via change of strain or bioreactor design or operation conditions. The outcome can also be of value to design a representative scale-down reactor to facilitate strain studies.

Understanding gradients in industrial bioreactors

Gradients in industrial bioreactors have attracted substantial research attention since exposure to fluctuating environmental conditions has been shown to lead to changes in the metabolome, transcriptome as well as population heterogeneity in industrially relevant microorganisms. Such changes have also been found to impact key process parameters like the yield on substrate and the productivity. Hence, understanding gradients is important from both the academic and industrial perspectives. In this review the causes of gradients are outlined, along with their impact on microbial physiology. Quantifying the impact of gradients requires a detailed understanding of both fluid flow inside industrial equipment and microbial physiology. This review critically examines approaches used to investigate gradients including large-scale experimental work, computational methods and scale-down approaches. Avenues for future work have been highlighted, particularly the need for further coordinated development of both in silico and experimental tools which can be used to further the current understanding of gradients in industrial equipment.

Engineering Saccharomyces cerevisiae for Succinic Acid Production From Glycerol and Carbon Dioxide

Previously, our lab replaced the endogenous FAD-dependent pathway for glycerol catabolism in S. cerevisiae by the synthetic NAD-dependent dihydroxyacetone (DHA) pathway. The respective modifications allow the full exploitation of glycerol’s higher reducing power (compared to sugars) for the production of the platform chemical succinic acid (SA) via a reductive, carbon dioxide fixing and redox-neutral pathway in a production host robust for organic acid production. Expression cassettes for three enzymes converting oxaloacetate to SA in the cytosol (“SA module”) were integrated into the genome of UBR2_CBS-DHA, an optimized CEN.PK derivative. Together with the additional expression of the heterologous dicarboxylic acid transporter DCT-02 from Aspergillus niger, a maximum SA titer of 10.7 g/L and a yield of 0.22 ± 0.01 g/g glycerol was achieved in shake flask (batch) cultures. Characterization of the constructed strain under controlled conditions in a bioreactor supplying additional carbon dioxide revealed that the carbon balance was closed to 96%. Interestingly, the results of the current study indicate that the artificial “SA module” and endogenous pathways contribute to the SA production in a highly synergistic manner.

Physiological responses of Saccharomyces cerevisiae to industrially relevant conditions: Slow growth, low pH, and high CO2 levels

Engineered strains of Saccharomyces cerevisiae are used for industrial production of succinic acid. Optimal process conditions for dicarboxylic‐acid yield and recovery include slow growth, low pH, and high CO₂. To quantify and understand how these process parameters affect yeast physiology, this study investigates individual and combined impacts of low pH (3.0) and high CO₂ (50%) on slow‐growing chemostat and retentostat cultures of the reference strain S. cerevisiae CEN.PK113‐7D. Combined exposure to low pH and high CO₂ led to increased maintenance‐energy requirements and death rates in aerobic, glucose‐limited cultures. Further experiments showed that these effects were predominantly caused by low pH. Growth under ammonium‐limited, energy‐excess conditions did not aggravate or ameliorate these adverse impacts. Despite the absence of a synergistic effect of low pH and high CO₂ on physiology, high CO₂ strongly affected genome‐wide transcriptional responses to low pH. Interference of high CO₂ with low‐pH signaling is consistent with low‐pH and high‐CO₂ signals being relayed via common (MAPK) signaling pathways, notably the cell wall integrity, high‐osmolarity glycerol, and calcineurin pathways. This study highlights the need to further increase robustness of cell factories to low pH for carboxylic‐acid production, even in organisms that are already applied at industrial scale.