SummaryThe methylotrophic yeast Pichia pastoris (Komagataella spp.) is widely used as cell factory for recombinant protein production. In the past recent years, important breakthroughs in the systems‐level quantitative analysis of its physiology have been achieved. This wealth of information has allowed the development of genome‐scale metabolic models, which make new approaches possible for host cell and bioprocess engineering. Nevertheless, the predictive accuracy of the previous consensus model required to be upgraded and validated with new experimental data sets for P. pastoris growing on glycerol or methanol as sole carbon sources, two of the most relevant substrates for this cell factory. In this study, we have characterized P. pastoris growing in chemostat cultures using glycerol or methanol as sole carbon sources over a wide range of growth rates, thereby providing physiological data on the effect of growth rate and culture conditions on biomass macromolecular and elemental composition. In addition, these data sets were used to improve the performance of the P. pastoris consensus genomic‐scale metabolic model iMT1026. Thereupon, new experimentally determined bounds, including the representation of biomass composition for these growth conditions, have been incorporated. As a result, here, we present version 3 (v3.0) of the consensus P. pastoris genome‐scale metabolic model as an update of the iMT1026 model. The v3.0 model was validated for growth on glycerol and methanol as sole carbon sources, demonstrating improved prediction capabilities over an extended substrate range including two biotechnologically relevant carbon sources.

MotivationGenome-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.ResultsIn 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.

High-level expression and secretion of heterologous proteins in yeast cause an increased energy demand, which may result in altered metabolic flux distributions. Moreover, recombinant protein overproduction often results in endoplasmic reticulum (ER) stress and oxidative stress, causing deviations from the optimal NAD(P)H regeneration balance. In this context, overexpression of genes encoding enzymes catalyzing endogenous NADPH-producing reactions, such as the oxidative branch of the pentose phosphate pathway, has been previously shown to improve protein production in Pichia pastoris (syn. Komagataella spp.). In this study, we evaluate the overexpression of the Saccharomyces cerevisiae POS5-encoded NADH kinase in a recombinant P. pastoris strain as an alternative approach to overcome such redox constraints. Specifically, POS5 was cooverexpressed in a strain secreting an antibody fragment, either by directing Pos5 to the cytosol or to the mitochondria. The physiology of the resulting strains was evaluated in continuous cultivations with glycerol or glucose as the sole carbon source, as well as under hypoxia (on glucose). Cytosolic targeting of Pos5 NADH kinase resulted in lower biomass-substrate yields but allowed for a 2-fold increase in product specific productivity. In contrast, Pos5 NADH kinase targeting to the mitochondria did not affect growth physiology and recombinant protein production significantly. Growth physiological parameters were in silico evaluated using the recent upgraded version (v3.0) of the P. pastoris consensus genome-scale metabolic model iMT1026, providing insights on the impact of POS5 overexpression on metabolic flux distributions. IMPORTANCE Recombinant protein overproduction often results in oxidative stress, causing deviations from the optimal redox cofactor regeneration balance. This becomes one of the limiting factors in obtaining high levels of heterologous protein production. Overexpression of redox-affecting enzymes has been explored in other organisms, such as Saccharomyces cerevisiae, as a means to fine tune the cofactor regeneration balance in order to obtain higher protein titers. In the present work, this strategy is explored in P. pastoris. In particular, one NADH kinase enzyme from S. cerevisiae (Pos5) is used, either in the cytosol or in mitochondria of P. pastoris, and its impact on the production of a model protein (antibody fragment) is evaluated. A significant improvement in the production of the model protein is observed when the kinase is directed to the cytosol. These results are significant in the field of heterologous protein production in general and in particular in the development of improved metabolic engineering strategies for P. pastoris.

No abstract

A full-length cDNA encoding digestive lipase (SmDL) was cloned from the pancreas of the smooth-hound (Mustelus mustelus). The obtained cDNA was 1350 bp long encoding 451 amino acids. The deduced amino acid sequence has high similarity with known pancreatic lipases. Catalytic triad and disulphide bond positions are also conserved. According to the established phylogeny, the SmDL was grouped with those of tuna and Sparidae lipases into one fish digestive lipase cluster. The recently purified enzyme shows no dependence for bile salts and colipase. For this, the residue-level interactions between lipase-colipase are yet to be clearly understood. The structural model of the SmDL was built, and several dissimilarities were noticed when analyzing the SmDL amino acids corresponding to those involved in HPL binding to colipase. Interestingly, the C-terminal domain of SmDL which holds the colipase shows a significant role for colipase interaction. This is apt to prevent the interaction between fish lipase and the pancreatic colipase which and can provide more explanation on the fact that the classical colipase is unable to activate the SmDL.