Metabolic reconstruction and flux analysis of industrial Pichia yeasts

Komagataella phaffii as Emerging Model Organism in Fundamental Research

Komagataella phaffii (Pichia pastoris) is one of the most extensively applied yeast species in pharmaceutical and biotechnological industries, and, therefore, also called the biotech yeast. However, thanks to more advanced strain engineering techniques, it recently started to gain attention as model organism in fundamental research. So far, the most studied model yeast is its distant cousin, Saccharomyces cerevisiae. While these data are of great importance, they limit our knowledge to one organism only. Since the divergence of the two species 250 million years ago, K. phaffii appears to have evolved less rapidly than S. cerevisiae, which is why it remains more characteristic of the common ancient yeast ancestors and shares more features with metazoan cells. This makes K. phaffii a valuable model organism for research on eukaryotic molecular cell biology, a potential we are only beginning to fully exploit. As methylotrophic yeast, K. phaffii has the intriguing property of being able to efficiently assimilate methanol as a sole source of carbon and energy. Therefore, major efforts have been made using K. phaffii as model organism to study methanol assimilation, peroxisome biogenesis and pexophagy. Other research topics covered in this review range from yeast genetics including mating and sporulation behavior to other cellular processes such as protein secretion, lipid biosynthesis and cell wall biogenesis. In this review article, we compare data obtained from K. phaffii with S. cerevisiae and other yeasts whenever relevant, elucidate major differences, and, most importantly, highlight the big potential of using K. phaffii in fundamental research.

A Combinational Strategy for Effective Heterologous Production of Functional Human Lysozyme in Pichia pastoris

Human lysozyme (hLYZ), known for its bacteriolytic activity, is widely applied in the food and pharmaceutical industries as an antimicrobial agent. However, its extensive application was limited by its low large-scale production efficiency. In this study, a combinational method of integrating codon optimization, multiple gene copies, and ER molecular chaperone co-expression was developed to improve the heterologous production of hLYZ in Pichia pastoris GS115. Our results showed that increasing the copy number of the optimized hLYZ gene in P. pastoris could enhance its secretory production level up to 1.57-fold. The recombinant opt-hLYZ-6C strain that contains six copies of opt-hLYZ gene exhibited the highest mRNA transcription levels, giving the highest production of 0.22 ± 0.02 mg/mL of hLYZ in the medium supernatant with a bacteriolytic activity of 14,680 ± 300 U/mL against Micrococcus lysodeikticus in the shaking flask experiment. Moreover, co-overexpression of ER retention molecular chaperones, such as Pdi1 or Ero1, in the recombinant opt-hLYZ-6C strain both presented positive effects on the secretory production of hLYZ. Our further characterization indicated that tandem co-expression of Ero1 and Pdi1 together presented an added-up effect. The secretory production of hLYZ in the medium supernatant reached 0.34 ± 0.02 mg/mL of the recombinant opt-hLYZ-6C-EP strain in the shaking flask experiment, with a bacteriolytic activity of 21,200 ± 400 U/mL. Compared to the recombinant opt-hLYZ-1C strain, these final improvements were calculated as 2.43-fold and 2.30-fold on secretory protein levels and antibacterial activity, respectively. Finally, the recombinant opt-hLYZ-6C-EP strain was applied for high-density cultivation in 5 L of fermenter, in which the secretory yield of hLYZ reached 2.34 ± 0.02 mg/mL in the medium supernatant, with a bacteriolytic activity of 1.76 ± 0.02 × 10⁵ U/mL against M. lysodeikticus. All these numbers presented the highest heterologous production levels of hLYZ in microbial systems.

Integrating metabolic modeling and population heterogeneity analysis into optimizing recombinant protein production by Komagataella (Pichia) pastoris

The methylotrophic yeast Komagataella (Pichia) pastoris has become one of the most utilized cell factories for the production of recombinant proteins over the last three decades. This success story is linked to its specific physiological traits, i.e., the ability to grow at high cell density in inexpensive culture medium and to secrete proteins at high yield. Exploiting methanol metabolism is at the core of most P. pastoris-based processes but comes with its own challenges. Co-feeding cultures with glycerol/sorbitol and methanol is a promising approach, which can benefit from improved understanding and prediction of metabolic response. The development of profitable processes relies on the construction and selection of efficient producing strains from less efficient ones but also depends on the ability to master the bioreactor process itself. More specifically, how a bioreactor processes could be monitored and controlled to obtain high yield of production. In this review, new perspectives are detailed regarding a multi-faceted approach to recombinant protein production processes by P. pastoris; including gaining improved understanding of the metabolic pathways involved, accounting for variations in transcriptional and translational efficiency at the single cell level and efficient monitoring and control of methanol levels at the bioreactor level.

13 C metabolic flux profiling of Pichia pastoris grown in aerobic batch cultures on glucose revealed high relative anabolic use of TCA cycle and limited incorporation of provided precursors of branched-chain amino acids

Carbon metabolism of Crabtree-negative yeast Pichia pastoris was profiled using ¹³ C nuclear magnetic resonance (NMR) to delineate regulation during exponential growth and to study the import of two precursors for branched-chain amino acid biosynthesis, α-ketoisovalerate and α-ketobutyrate. Cells were grown in aerobic batch cultures containing (a) only glucose, (b) glucose along with the precursors, or (c) glucose and Val. The study provided the following new insights. First, ¹³ C flux ratio analyses of central metabolism reveal an unexpectedly high anaplerotic supply of the tricarboxylic acid cycle for a Crabtree-negative yeast, and show that a substantial fraction of glucose catabolism proceeds through the pentose phosphate pathway. A comparison with previous flux ratio analyses for batch cultures of Crabtree-negative Pichia stipitis and Crabtree-positive Saccharomyces cerevisiae indicate that the overall regulation of central carbon metabolism in P. pastoris is intermediate in between P. stipitis and S. cerevisiae. Second, excess α-ketoisovalerate in the medium is not transported into the cytoplasm indicating that P. pastoris lacks a suitable transporter. In contrast, excess Val is efficiently taken up and largely fulfills demands for both Val and Leu for protein synthesis. Third, excess α-ketobutyrate is transported into the mitochondria for Ile biosynthesis. However, the import does not efficiently inhibit the synthesis of α-ketobutyrate from pyruvate indicating that P. pastoris has not been optimized evolutionarily to take full advantage of this carbon source. These findings have direct implications for preparing uniformly ² H,¹³ C,¹⁵ N-labeled proteins containing protonated Ile, Val, and Leu methyl groups in P. pastoris for NMR-based structural biology.

ENZYMES

Acetohydroxy acid isomeroreductase (EC 1.1.1.86), branched-chain amino acid aminotransferase (BCAT, EC 2.6.1.42), fumarase (EC 4.2.1.2), malic enzyme (EC 1.1.1.39/1.1.1.40), phosphoenolpyruvate carboxykinase (EC 4.1.1.49), pyruvate carboxylase (EC 6.4.1.1), pyruvate kinase (EC 2.7.1.40), l-serine hydroxymethyltransferase (EC 2.1.2.1), threonine aldolase (EC 4.1.2.5), threonine dehydratase (EC 4.3.1.19); transketolase (EC 2.2.1.1), transaldolase (EC 2.2.1.2).

Characterization of a panARS-based episomal vector in the methylotrophic yeast Pichia pastoris for recombinant protein production and synthetic biology applications

Background

Recombinant protein production in the methylotrophic yeast Pichia pastoris largely relies on integrative vectors. Although the stability of integrated expression cassettes is well appreciated for most applications, the availability of reliable episomal vectors for this host would represent a useful tool to expedite cloning and high-throughput screening, ameliorating also the relatively high clonal variability reported in transformants from integrative vectors caused by off-target integration in the P. pastoris genome. Recently, heterologous and endogenous autonomously replicating sequences (ARS) were identified in P. pastoris by genome mining, opening the possibility of expanding the available toolbox to include efficient episomal plasmids. The aim of this technical report is to validate a 452-bp sequence (“panARS”) in context of P. pastoris expression vectors, and to compare their performance to classical integrative plasmids. Moreover, we aimed to test if such episomal vectors would be suitable to sustain in vivo recombination, using fragments for transformation, directly in P. pastoris cells.

Results

A panARS-based episomal vector was evaluated using blue fluorescent protein (BFP) as a reporter gene. Normalized fluorescence from colonies carrying panARS-BFP outperformed the level of signal obtained from integrative controls by several-fold, whereas endogenous sequences, identified from the P. pastoris genome, were not as efficient in terms of protein production. At the single cell level, panARS-BFP clones showed lower interclonal variability but higher intraclonal variation compared to their integrative counterparts, supporting the idea that heterologous protein production could benefit from episomal plasmids. Finally, efficiency of 2-fragment and 3-fragment in vivo recombination was tested using varying lengths of overlapping regions and molar ratios between fragments. Upon optimization, minimal background was obtained for in vivo assembled vectors, suggesting this could be a quick and efficient method to generate of episomal plasmids of interest.

Conclusions

An expression vector based on the panARS sequence was shown to outperform its integrative counterparts in terms of protein productivity and interclonal variability, facilitating recombinant protein expression and screening. Using optimized fragment lengths and ratios, it was possible to perform reliable in vivo recombination of fragments in P. pastoris. Taken together, these results support the applicability of panARS episomal vectors for synthetic biology approaches.

Integration and Validation of the Genome-Scale 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.

A blueprint of the amino acid biosynthesis network of hemiascomycetes

The structure and regulation of biosynthesis pathways in Saccharomyces cerevisiae have been detailed extensively. For other hemiascomycetes, genomic sequences are primarily available, whereas biochemical information on them is scarce. The resulting biochemical networks that are used for research in basic science and biotechnology are often biased by data from S. cerevisiae, assuming that there are often implicitly conserved structures between species. We examined the structure of the amino acid biosynthesis network in nine hemiascomycetes, spanning the phylogenetic clade. Differences in the genetic inventory included the presence and absence of isoenzymes and compartmentation of the pathways. Notably, no two hemiascomycetes had identical genetic inventories. For example, the lack of the mitochondrial αIPMS isoenzyme and presence of only one copy of the BCAA aminotransferase in Pichia pastoris indicate a disparately compartmented leucine biosynthesis pathway. Our findings suggest that αIPMS and BCAA aminotransferase are solely located in the cytosol of P. pastoris, requiring correction of the leucine biosynthesis pathway layout in this species. Our results argue for careful use of information from S. cerevisiae and for joint efforts to fill the knowledge gaps in other species. Such analysis will lead to contributions in biotechnology disciplines, such as protein production and compartment engineering.

(Post-)genomics approaches in fungal research

To date, hundreds of fungal genomes have been sequenced and many more are in progress. This wealth of genomic information has provided new directions to study fungal biodiversity. However, to further dissect and understand the complicated biological mechanisms involved in fungal life styles, functional studies beyond genomes are required. Thanks to the developments of current -omics techniques, it is possible to produce large amounts of fungal functional data in a high-throughput fashion (e.g. transcriptome, proteome, etc.). The increasing ease of creating -omics data has also created a major challenge for downstream data handling and analysis. Numerous databases, tools and software have been created to meet this challenge. Facing such a richness of techniques and information, hereby we provide a brief roadmap on current wet-lab and bioinformatics approaches to study functional genomics in fungi.

¹³C-based metabolic flux analysis of recombinant Pichia pastoris

Overexpression of a foreign protein may negatively affect several cell growth parameters, as well as cause cellular stress. Central (or core) metabolism plays a crucial role since it supplies energy, reduction equivalents, and precursor molecules for the recombinant product, cell's maintenance, and growth needs. However, the number of quantitative physiology studies of the impact of recombinant protein production on the central metabolic pathways of yeast cell factories has been traditionally rather limited, thereby hampering the application of rational strain engineering strategies targeting central metabolism.The development and application of quantitative physiology and modelling tools and methodologies is allowing for a systems-level understanding of the effect of bioprocess parameters such as growth rate, temperature, oxygen availability, and substrate(s) choice on metabolism, and its subsequent interactions with recombinant protein synthesis, folding, and secretion.Here, we review the recent developments and applications of (13)C-based metabolic flux analysis ((13)C-MFA) of Pichia pastoris and the gained understanding of the metabolic behavior of this yeast in recombinant protein production bioprocesses. We also discuss the potential of multilevel studies integrating (13)C-MFA with other omics analyses, as well as future perspectives on the metabolic modelling approaches to study and design metabolic engineering strategies for improved protein production.