Integration and Validation of the GenomeScale Metabolic Models of Pichia pastoris : a Comprehensive Update of Protein Glycosylation Pathways, Lipid and Energy Metabolism

Motivation: Genome-scale metabolic models (GEMs) are tools that allow predicting a phenotype from a genotype under certain environmental conditions. GEMs have been developed in the last ten years for a broad range of organisms, and are used for multiple purposes such as discovering new properties of metabolic networks, predicting new targets for metabolic engineering, as well as optimizing the cultivation conditions for biochemicals or recombinant protein production. Pichia pastoris is one of the most widely used organisms for heterologous protein expression. There are different GEMs for this methylotrophic yeast of which the most relevant and complete in the published literature are iPP668, PpaMBEL1254 and iLC915. However, these three models differ regarding certain pathways, terminology for metabolites and reactions and annotations. Moreover, GEMs for some species are typically built based on the reconstructed models of related model organisms. In these cases, some organism-specific pathways could be missing or misrepresented. Results: In order to provide an updated and more comprehensive GEM for P. pastoris, we have reconstructed and validated a consensus model integrating and merging all three existing models. In this step a comprehensive review and integration of the metabolic pathways included in each one of these three versions was performed. In addition, the resulting iMT1026 model includes a new description of some metabolic processes. Particularly new information described in recently published literature is included, mainly related to fatty acid and sphingolipid metabolism, glycosylation and cell energetics. Finally the reconstructed model was tested and validated, by comparing the results of the simulations with available empirical physiological datasets results obtained from a wide range of experimental conditions, such as different carbon sources, distinct oxygen availability conditions, as well as producing of two different recombinant proteins. In these simulations, the iMT1026 model has shown a better performance than the previous existing models

Contextualized genome-scale model unveils high-order metabolic effects of the specific growth rate and oxygenation level in recombinant Pichia pastoris

Pichia pastoris is recognized as a biotechnological workhorse for recombinant protein expression. The metabolic performance of this microorganism depends on genetic makeup and culture conditions, amongst which the specific growth rate and oxygenation level are critical. Despite their importance, only their individual effects have been assessed so far, and thus their combined effects and metabolic consequences still remain to be elucidated. In this work, we present a comprehensive framework for revealing high-order (i.e., individual and combined) metabolic effects of the above parameters in glucose-limited continuous cultures of P. pastoris, using thaumatin production as a case study. Specifically, we employed a rational experimental design to calculate statistically significant metabolic effects from multiple chemostat data, which were later contextualized using a refined and highly predictive genome-scale metabolic model of this yeast under the simulated conditions. Our results revealed a negative effect of the oxygenation on the specific product formation rate (thaumatin), and a positive effect on the biomass yield. Notably, we identified a novel positive combined effect of both the specific growth rate and oxygenation level on the specific product formation rate. Finally, model predictions indicated an opposite relationship between the oxygenation level and the growth-associated maintenance energy (GAME) requirement, suggesting a linear GAME decrease of 0.56 mmol ATP/g per each 1% increase in oxygenation level, which translated into a 44% higher metabolic cost under low oxygenation compared to high oxygenation. Overall, this work provides a systematic framework for mapping high-order metabolic effects of different culture parameters on the performance of a microbial cell factory. Particularly in this case, it provided valuable insights about optimal operational conditions for protein production in P. pastoris

Quantitative Metabolomics and Instationary 13C-Metabolic Flux Analysis Reveals Impact of Recombinant Protein Production on Trehalose and Energy Metabolism in Pichia pastoris

Pichia pastoris has been recognized as an effective host for recombinant protein production. In this work, we combine metabolomics and instationary 13C metabolic flux analysis (INST 13C-MFA) using GC-MS and LC-MS/MS to evaluate the potential impact of the production of a Rhizopus oryzae lipase (Rol) on P. pastoris central carbon metabolism. Higher oxygen uptake and CO2 production rates and slightly reduced biomass yield suggest an increased energy demand for the producing strain. This observation is further confirmed by 13C-based metabolic flux analysis. In particular, the flux through the methanol oxidation pathway and the TCA cycle was increased in the Rol-producing strain compared to the reference strain. Next to changes in the flux distribution, significant variations in intracellular metabolite concentrations were observed. Most notably, the pools of trehalose, which is related to cellular stress response, and xylose, which is linked to methanol assimilation, were significantly increased in the recombinant strain

Continuous Cultivation as a Tool Toward the Rational Bioprocess Development With Pichia Pastoris Cell Factory

The methylotrophic yeast Pichia pastoris (Komagataella phaffii) is currently considered one of the most promising hosts for recombinant protein production (RPP) and metabolites due to the availability of several tools to efficiently regulate the recombinant expression, its ability to perform eukaryotic post-translational modifications and to secrete the product in the extracellular media. The challenge of improving the bioprocess efficiency can be faced from two main approaches: the strain engineering, which includes enhancements in the recombinant expression regulation as well as overcoming potential cell capacity bottlenecks; and the bioprocess engineering, focused on the development of rational-based efficient operational strategies. Understanding the effect of strain and operational improvements in bioprocess efficiency requires to attain a robust knowledge about the metabolic and physiological changes triggered into the cells. For this purpose, a number of studies have revealed chemostat cultures to provide a robust tool for accurate, reliable, and reproducible bioprocess characterization. It should involve the determination of key specific rates, productivities, and yields for different C and N sources, as well as optimizing media formulation and operating conditions. Furthermore, studies along the different levels of systems biology are usually performed also in chemostat cultures. Transcriptomic, proteomic and metabolic flux analysis, using different techniques like differential target gene expression, protein description and 13 C-based metabolic flux analysis, are widely described as valued examples in the literature. In this scenario, the main advantage of a continuous operation relies on the quality of the homogeneous samples obtained under steady-state conditions, where both the metabolic and physiological status of the cells remain unaltered in an all-encompassing picture of the cell environment. This contribution aims to provide the state of the art of the different approaches that allow the design of rational strain and bioprocess engineering improvements in Pichia pastoris toward optimizing bioprocesses based on the results obtained in chemostat cultures. Interestingly, continuous cultivation is also currently emerging as an alternative operational mode in industrial biotechnology for implementing continuous process operations

Rational development of bioprocess engineering strategies for recombinant protein production in Pichia pastoris (Komagataella phaffii) using the methanol-free GAP promoter. Where do we stand?

The increasing demand for recombinant proteins for a wide range of applications, from biopharmaceutical protein complexes to industrial enzymes, is leading to important growth in this market. Among the different efficient host organism alternatives commonly used for protein production, the yeast Pichia pastoris (Komagataella phaffii) is currently considered to be one of the most effective and versatile expression platforms. The promising features of this cell factory are giving rise to interesting studies covering the different aspects that contribute to improving the bioprocess efficiency, from strain engineering to bioprocess engineering. The numerous drawbacks of using methanol in industrial processes are driving interest towards methanol-free alternatives, among which the GAP promoter-based systems stand out. The aim of this work is to present the most promising innovative developments in operational strategies based on rational approaches through bioprocess engineering tools. This rational design should be based on physiological characterization of the producing strains under bioprocess conditions and its interrelation with specific rates. This review focuses on understanding the key factors that can enhance recombinant protein production in Pichia pastoris; they are the basis for a further discussion on future industrial applications with the aim of developing scalable alternative strategies that maximize yields and productivity

Increased dosage of AOX1 promoter-regulated expression cassettes leads to transcription attenuation of the methanol metabolism in Pichia pastoris

The methanol-regulated alcohol oxidase promoter (PAOX1) of Pichia pastoris is one of the strongest promoters for heterologous gene expression in this methylotrophic yeast. Although increasing gene dosage is one of the most common strategies to increase recombinant protein productivities, the increase of gene dosage of Rhizopus oryzae lipase (ROL) in P. pastoris has been previously shown to reduce cell growth, lipase production and substrate consumption in high-copy strains. To better assess that physiological response, transcriptomics analysis was performed of a subset of strains with 1 to 15 ROL copies. The macroscopic physiological parameters confirm that growth yield and carbon uptake rate are gene dosage dependent, and were supported by the transcriptomic data, showing the impact of increased dosage of AOX1 promoter-regulated expression cassettes on P. pastoris physiology under steady methanolic growth conditions. Remarkably, increased number of cassettes led to transcription attenuation of the methanol metabolism and peroxisome biogenesis in P. pastoris, concomitant with reduced secretion levels of the heterologous product. Moreover, our data also point to a block in ROL mRNA translation in the higher ROL-copies constructs, while the low productivities of multi-copy strains under steady growth conditions do not appear to be directly related to UPR and ERAD induction.

Quantitative metabolomics analysis of amino acid metabolism in recombinant Pichia pastoris under different oxygen availability conditions

Background: Environmental and intrinsic stress factors can result in the global alteration of yeast physiology, as evidenced by several transcriptional studies. Hypoxia has been shown to have a beneficial effect on the expression of recombinant proteins in Pichia pastoris growing on glucose. Furthermore, transcriptional profiling analyses revealed that oxygen availability was strongly affecting ergosterol biosynthesis, central carbon metabolism and stress responses, in particular the unfolded protein response. To contribute to the better understanding of the effect and interplay of oxygen availability and foreign protein secretion on central metabolism, a first quantitative metabolomic analysis of free amino acids pools in a recombinant P. pastoris strain growing under different oxygen availability conditions has been performed. Results: The values obtained indicate significant variations in the intracellular amino acid pools due to different oxygen availability conditions, showing an overall increase of their size under oxygen limitation. Notably, even while foreign protein productivities were relatively low (about 40-80 μg Fab/gDCW·h), recombinant protein production was found to have a limited but significant impact on the intracellular amino acid pools, which were generally decreased in the producing strain compared with the reference strain. However, observed changes in individual amino acids pools were not correlated with their corresponding relative abundance in the recombinant protein sequence, but to the overall cell protein amino acid compositional variations. Conclusions: Overall, the results obtained, combined with previous transcriptomic and proteomic analyses provide a systematic metabolic fingerprint of the oxygen availability impact on recombinant protein production in P. pastoris

Benchmarking recombinant Pichia pastoris for 3-hydroxypropionic acid production from glycerol

The use of the methylotrophic yeast Pichia pastoris (Komagataella phaffi) to produce heterologous proteins has been largely reported. However, investigations addressing the potential of this yeast to produce bulk chemicals are still scarce. In this study, we have studied the use of P. pastoris as a cell factory to produce the commodity chemical 3-hydroxypropionic acid (3-HP) from glycerol. 3-HP is a chemical platform which can be converted into acrylic acid and to other alternatives to petroleum-based products. To this end, the mcr gene from Chloroflexus aurantiacus was introduced into P. pastoris. This single modification allowed the production of 3-HP from glycerol through the malonyl-CoA pathway. Further enzyme and metabolic engineering modifications aimed at increasing cofactor and metabolic precursors availability allowed a 14-fold increase in the production of 3-HP compared to the initial strain. The best strain (PpHP6) was tested in a fed-batch culture, achieving a final concentration of 3-HP of 24.75 g l−1, a product yield of 0.13 g g−1 and a volumetric productivity of 0.54 g l−1 h−1, which, to our knowledge, is the highest volumetric productivity reported in yeast. These results benchmark P. pastoris as a promising platform to produce bulk chemicals for the revalorization of crude glycerol and, in particular, to produce 3-HP

A multi-level study of recombinant Pichia pastoris in different oxygen conditions

Background: Yeasts are attractive expression platforms for many recombinant proteins, and there is evidence for an important interrelation between the protein secretion machinery and environmental stresses. While adaptive responses to such stresses are extensively studied in Saccharomyces cerevisiae, little is known about their impact on the physiology of Pichia pastoris. We have recently reported a beneficial effect of hypoxia on recombinant Fab secretion in P. pastoris chemostat cultivations. As a consequence, a systems biology approach was used to comprehensively identify cellular adaptations to low oxygen availability and the additional burden of protein production. Gene expression profiling was combined with proteomic analyses and the 13C isotope labelling based experimental determination of metabolic fluxes in the central carbon metabolism. Results: The physiological adaptation of P. pastoris to hypoxia showed distinct traits in relation to the model yeast S. cerevisiae. There was a positive correlation between the transcriptomic, proteomic and metabolic fluxes adaptation of P. pastoris core metabolism to hypoxia, yielding clear evidence of a strong transcriptional regulation component of key pathways such as glycolysis, pentose phosphate pathway and TCA cycle. In addition, the adaptation to reduced oxygen revealed important changes in lipid metabolism, stress responses, as well as protein folding and trafficking. Conclusions: This systems level study helped to understand the physiological adaptations of cellular mechanisms to low oxygen availability in a recombinant P. pastoris strain. Remarkably, the integration of data from three different levels allowed for the identification of differences in the regulation of the core metabolism between P. pastoris and S. cerevisiae. Detailed comparative analysis of the transcriptomic data also led to new insights into the gene expression profiles of several cellular processes that are not only susceptible to low oxygen concentrations, but might also contribute to enhanced protein secretion

Integration and Validation of the GenomeScale Metabolic Models of Pichia pastoris : a Comprehensive Update of Protein Glycosylation Pathways, Lipid and Energy Metabolism

Motivation: Genome-scale metabolic models (GEMs) are tools that allow predicting a phenotype from a genotype under certain environmental conditions. GEMs have been developed in the last ten years for a broad range of organisms, and are used for multiple purposes such as discovering new properties of metabolic networks, predicting new targets for metabolic engineering, as well as optimizing the cultivation conditions for biochemicals or recombinant protein production. Pichia pastoris is one of the most widely used organisms for heterologous protein expression. There are different GEMs for this methylotrophic yeast of which the most relevant and complete in the published literature are iPP668, PpaMBEL1254 and iLC915. However, these three models differ regarding certain pathways, terminology for metabolites and reactions and annotations. Moreover, GEMs for some species are typically built based on the reconstructed models of related model organisms. In these cases, some organism-specific pathways could be missing or misrepresented. Results: In order to provide an updated and more comprehensive GEM for P. pastoris, we have reconstructed and validated a consensus model integrating and merging all three existing models. In this step a comprehensive review and integration of the metabolic pathways included in each one of these three versions was performed. In addition, the resulting iMT1026 model includes a new description of some metabolic processes. Particularly new information described in recently published literature is included, mainly related to fatty acid and sphingolipid metabolism, glycosylation and cell energetics. Finally the reconstructed model was tested and validated, by comparing the results of the simulations with available empirical physiological datasets results obtained from a wide range of experimental conditions, such as different carbon sources, distinct oxygen availability conditions, as well as producing of two different recombinant proteins. In these simulations, the iMT1026 model has shown a better performance than the previous existing models