Engineering the synthetic β-alanine pathway in Komagataella phaffii for conversion of methanol into 3-hydroxypropionic acid

BACKGROUND

Methanol is increasingly gaining attraction as renewable carbon source to produce specialty and commodity chemicals, as it can be generated from renewable sources such as carbon dioxide (CO2). In this context, native methylotrophs such as the yeast Komagataella phaffii (syn Pichia pastoris) are potentially attractive cell factories to produce a wide range of products from this highly reduced substrate. However, studies addressing the potential of this yeast to produce bulk chemicals from methanol are still scarce. 3-Hydroxypropionic acid (3-HP) is a platform chemical which can be converted into acrylic acid and other commodity chemicals and biopolymers. 3-HP can be naturally produced by several bacteria through different metabolic pathways.

RESULTS

In this study, production of 3-HP via the synthetic β-alanine pathway has been established in K. phaffii for the first time by expressing three heterologous genes, namely panD from Tribolium castaneum, yhxA from Bacillus cereus, and ydfG from Escherichia coli K-12. The expression of these key enzymes allowed a production of 1.0 g l-1 of 3-HP in small-scale cultivations using methanol as substrate. The addition of a second copy of the panD gene and selection of a weak promoter to drive expression of the ydfG gene in the PpCβ21 strain resulted in an additional increase in the final 3-HP titer (1.2 g l-1). The 3-HP-producing strains were further tested in fed-batch cultures. The best strain (PpCβ21) achieved a final 3-HP concentration of 21.4 g l-1 after 39 h of methanol feeding, a product yield of 0.15 g g-1, and a volumetric productivity of 0.48 g l-1 h-1. Further engineering of this strain aiming at increasing NADPH availability led to a 16% increase in the methanol consumption rate and 10% higher specific productivity compared to the reference strain PpCβ21.

CONCLUSIONS

Our results show the potential of K. phaffii as platform cell factory to produce organic acids such as 3-HP from renewable one-carbon feedstocks, achieving the highest volumetric productivities reported so far for a 3-HP production process through the β-alanine pathway.

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.

Combining Metabolic Engineering and Multiplexed Screening Methods for 3-Hydroxyprionic Acid Production in Pichia pastoris

Production of 3-hydroxypropionic acid (3-HP) in Pichia pastoris (syn. Komagataella phaffii) via the malonyl-CoA pathway has been recently demonstrated using glycerol as a carbon source, but the reported metrics were not commercially relevant. The flux through the heterologous pathway from malonyl-CoA to 3-HP was hypothesized as the main bottleneck. In the present study, different metabolic engineering approaches have been combined to improve the productivity of the original 3-HP producing strains. To do so, an additional copy of the gene encoding for the potential rate-limiting step of the pathway, i.e., the C-terminal domain of the malonyl-CoA reductase, was introduced. In addition, a variant of the endogenous acetyl-CoA carboxylase (ACC1S1132A) was overexpressed with the aim to increase the delivery of malonyl-CoA. Furthermore, the genes encoding for the pyruvate decarboxylase, aldehyde dehydrogenase and acetyl-CoA synthase, respectively, were overexpressed to enhance conversion of pyruvate into cytosolic acetyl-CoA, and the main gene responsible for the production of the by-product D-arabitol was deleted. Three different screening conditions were used to classify the performance of the different strains: 24-deep-well plates batch cultures, small-scale cultures in falcon tubes using FeedBeads® (i.e., slow release of glycerol over time), and mini bioreactor batch cultures. The best two strains from the FeedBeads® screening, PpHP8 and PpHP18, were tested in bioreactor fed-batch cultures using a pre-fixed exponentially increasing feeding rate. The strain PpHP18 produced up to 37.05 g L−1 of 3-HP at 0.712 g L−1 h−1 with a final product yield on glycerol of 0.194 Cmol−1 in fed-batch cultures. Remarkably, PpHP18 did not rank among the 2-top producer strains in small scale batch cultivations in deep-well plates and mini bioreactors, highlighting the importance of multiplexed screening conditions for adequate assessment of metabolic engineering strategies. These results represent a 50% increase in the product yield and final concentration, as well as over 30% increase in volumetric productivity compared to the previously obtained metrics for P. pastoris. Overall, the combination of glycerol as carbon source and a metabolically engineered P. pastoris strain resulted in the highest 3-HP concentration and productivity reported so far in yeast.

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.

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 ¹³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.

Redox Engineering by Ectopic Overexpression of NADH Kinase in Recombinant Pichia pastoris (Komagataella phaffii): Impact on Cell Physiology and Recombinant Production of Secreted Proteins

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 cerevisiaePOS5-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.

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_DCW 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.

Deregulation of methanol metabolism reverts transcriptional limitations of recombinant Pichia pastoris (Komagataella spp) with multiple expression cassettes under control of the AOX1 promoter

The methanol-regulated alcohol oxidase promoter (P_AOX1 ) of Pichia pastoris (syn. Komagataella spp. ) is one of the strongest promoters for heterologous gene expression. Although increasing the gene dosage is a common strategy to improve recombinant protein productivities, P. pastoris strains harboring more than two copies of a Rhizopus oryzae lipase gene (ROL) have previously shown a decrease in cell growth, lipase production, and substrate consumption, as well as a significant transcriptional downregulation of methanol metabolism. This pointed to a potential titration effect of key transcriptional factors methanol expression regulator 1 (Mxr1) and methanol-induced transcription factor (Mit1) regulating methanol metabolism caused by the insertion of multiple expression vectors. To prove this hypothesis, a set of strains carrying one and four copies of ROL (1C and 4C, respectively) were engineered to coexpress one or two copies of MXR1*, coding for an Mxr1 variant insensitive to repression by 14-3-3 regulatory proteins, or one copy of MIT1. Small-scale cultures revealed that growth, Rol productivity, and methanol consumption were improved in the 4C-MXR1* and 4C-MIT1, strains growing on methanol as a sole carbon source, whereas only a slight increase in productivity was observed for re-engineered 1C strains. We further verified the improved performance of these strains in glycerol-/methanol-limited chemostat cultures.

Metabolic flux balance analysis during lactate and glucose concomitant consumption in HEK293 cell cultures

At early stages of the exponential growth phase in HEK293 cell cultures, the tricarboxylic acid cycle is unable to process all the amount of NADH generated in the glycolysis pathway, being lactate the main by-product. However, HEK293 cells are also able to metabolize lactate depending on the environmental conditions. It has been recently observed that one of the most important modes of lactate metabolization is the cometabolism of lactate and glucose, observed even during the exponential growth phase. Extracellular lactate concentration and pH appear to be the key factors triggering the metabolic shift from glucose consumption and lactate production to lactate and glucose concomitant consumption. The hypothesis proposed for triggering this metabolic shift to lactate and glucose concomitant consumption is that HEK293 cells metabolize extracellular lactate as a response to both extracellular protons and lactate accumulation, by means of cotransporting them (extracellular protons and lactate) into the cytosol. At this point, there exists a considerable controversy about how lactate reaches the mitochondrial matrix: the first hypothesis proposes that lactate is converted into pyruvate in the cytosol, and afterward, pyruvate enters into the mitochondria; the second alternative considers that lactate enters first into the mitochondria, and then, is converted into pyruvate. In this study, lactate transport and metabolization into mitochondria is shown to be feasible, as evidenced by means of respirometry tests with isolated active mitochondria, including the depletion of lactate concentration of the respirometry assay. Although the capability of lactate metabolization by isolated mitochondria is demonstrated, the possibility of lactate being converted into pyruvate in the cytosol cannot be excluded from the discussion. For this reason, the calculation of the metabolic fluxes for an HEK293 cell line was performed for the different metabolic phases observed in batch cultures under pH controlled and noncontrolled conditions, considering both hypotheses. The main objective of this study is to evaluate the redistribution of cellular metabolism and compare the differences or similarities between the phases before and after the metabolic shift of HEK293 cells (shift observed when pH is not controlled). That is from a glucose consumption/lactate production phase to a glucose-lactate coconsumption phase. Interestingly, switching to a glucose and lactate cometabolization results in a better-balanced cell metabolism, with decreased glucose and amino acids uptake rates, affecting minimally cell growth. This behavior could be applied to further develop new approaches in terms of cell engineering and to develop improved cell culture strategies in the field of animal cell technology.

Metabolic flux analysis and the NAD(P)H/NAD(P)+ ratios in chemostat cultures of Azotobacter vinelandii

BACKGROUND

Azotobacter vinelandii is a bacterium that produces alginate and polyhydroxybutyrate (P3HB); however, the role of NAD(P)H/NAD(P)+ ratios on the metabolic fluxes through biosynthesis pathways of these biopolymers remains unknown. The aim of this study was to evaluate the NAD(P)H/NAD(P) + ratios and the metabolic fluxes involved in alginate and P3HB biosynthesis, under oxygen-limiting and non-limiting oxygen conditions.

RESULTS

The results reveal that changes in the oxygen availability have an important effect on the metabolic fluxes and intracellular NADPH/NADP+ ratio, showing that at the lowest OTR (2.4 mmol L-1 h-1), the flux through the tricarboxylic acid (TCA) cycle decreased 27.6-fold, but the flux through the P3HB biosynthesis increased 6.6-fold in contrast to the cultures without oxygen limitation (OTR = 14.6 mmol L-1 h-1). This was consistent with the increase in the level of transcription of phbB and the P3HB biosynthesis. In addition, under conditions without oxygen limitation, there was an increase in the carbon uptake rate (twofold), as well as in the flux through the pentose phosphate (PP) pathway (4.8-fold), compared to the condition of 2.4 mmol L-1 h-1. At the highest OTR condition, a decrease in the NADPH/NADP+ ratio of threefold was observed, probably as a response to the high respiration rate induced by the respiratory protection of the nitrogenase under diazotrophic conditions, correlating with a high expression of the uncoupled respiratory chain genes (ndhII and cydA) and induction of the expression of the genes encoding the nitrogenase complex (nifH).

CONCLUSIONS

We have demonstrated that changes in oxygen availability affect the internal redox state of the cell and carbon metabolic fluxes. This also has a strong impact on the TCA cycle and PP pathway as well as on alginate and P3HB biosynthetic fluxes.

Fine-tuning the P. pastoris iMT1026 genome-scale metabolic model for improved prediction of growth on methanol or glycerol as sole carbon sources

The 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.

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 (P_AOX1) 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.

Genome-scale metabolic reconstruction for the insidious bacterium in aquaculture Piscirickettsia salmonis

Piscirickettsia salmonis is a fish bacterium that causes the disease piscirickettsiosis in salmonids. This pathology is partially controlled by vaccines. The lack of knowledge has hindered its culture on laboratory and industrial scale. The study describes the metabolic phenotype of P. salmonis in culture. This study presents the first genome-scale model (iPF215) of the LF-89 strain of P. salmonis, describing the central metabolic pathway, biosynthesis and molecule degradation and transport mechanisms. The model was adjusted with experiment data, allowing the identification of the capacities that were not predicted by the automatic annotation of the genome sequences. The iPF215 model is comprised of 417 metabolites, 445 reactions and 215 genes, was used to reproduce the growth of P. salmonis (μ_max 0.052±0.005h^-1). The metabolic reconstruction of the P. salmonis LF-89 strain obtained in this research provides a baseline that describes the metabolic capacities of the bacterium and is the basis for developing improvements to its cultivation for vaccine formulation.

Droplet digital PCR-aided screening and characterization of Pichia pastoris multiple gene copy strains

Pichia (syn. Komagataella) pastoris is a widely used yeast platform for heterologous protein production. Expression cassettes are usually stably integrated into the genome of this host via homologous recombination. Although increasing gene dosage is a powerful strategy to improve recombinant protein production, an excess in the number of gene copies often leads to decreased product yields and increased metabolic burden, particularly for secreted proteins. We have constructed a series of strains harboring different copy numbers of a Rhizopus oryzae lipase gene (ROL), aiming to find the optimum gene dosage for secreted Rol production. In order to accurately determine ROL gene dosage, we implemented a novel protocol based on droplet digital PCR (ddPCR), and cross validated it with conventional real-time PCR. Gene copy number determination based on ddPCR allowed for an accurate ranking of transformants according to their ROL gene dosage. Results indicated that ddPCR was particularly superior at lower gene dosages (one to five copies) over quantitative real-time PCR (qPCR). This facilitated the determination of the optimal ROL gene dosage as low as two copies. The ranking of ROL gene dosage versus Rol yield was consistent at both small scale and bioreactor chemostat cultures, thereby easing clone characterization in terms of gene dosage dependent physiological effects, which could be discriminated even among strains differing by only one ROL copy. A selected two-copy strain showed twofold increase in Rol specific production in a chemostat culture over the single copy strain. Conversely, strains harboring more than two copies of the ROL gene showed decreased product and biomass yields, as well as altered substrate consumption specific rates, compared to the reference (one-copy) strain. Biotechnol. Bioeng. 2016;113: 1542-1551. © 2015 Wiley Periodicals, Inc.

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.

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.

¹³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.

¹³C-based metabolic flux analysis in yeast: the Pichia pastoris case

Metabolic flux analysis based on tracing patterns of stable isotopes, particularly (13)C, comprises a set of methodologies to experimentally quantify intracellular biochemical reaction rates, i.e., to measure carbon flux distributions through a metabolic network. This allows quantifying the response of a metabolic network to an environmental or genetic perturbation (i.e., the metabolic phenotype). Here, we describe a protocol based on growing yeast on a (13)C-labelled substrate and subsequent NMR detection of (13)C-patterns in proteinogenic amino acids. To calculate metabolic fluxes, we describe two complementary mathematical approaches using available software; namely, an approach based on the estimation of local ratios in network nodes, and a method based on a global iterative fitting approach. Furthermore, we consider specificities of these protocols for their application to the yeast Pichia pastoris growing on multicarbon substrates other than glucose (glycerol), as well as the case when methanol is used as co-substrate in combination with glucose or glycerol.

Metabolic flux analysis during the exponential growth phase of Saccharomyces cerevisiae in wine fermentations

As a consequence of the increase in global average temperature, grapes with the adequate phenolic and aromatic maturity tend to be overripe by the time of harvest, resulting in increased sugar concentrations and imbalanced C/N ratios in fermenting musts. This fact sets obvious additional hurdles in the challenge of obtaining wines with reduced alcohols levels, a new trend in consumer demands. It would therefore be interesting to understand Saccharomyces cerevisiae physiology during the fermentation of must with these altered characteristics. The present study aims to determine the distribution of metabolic fluxes during the yeast exponential growth phase, when both carbon and nitrogen sources are in excess, using continuous cultures. Two different sugar concentrations were studied under two different winemaking temperature conditions. Although consumption and production rates for key metabolites were severely affected by the different experimental conditions studied, the general distribution of fluxes in central carbon metabolism was basically conserved in all cases. It was also observed that temperature and sugar concentration exerted a higher effect on the pentose phosphate pathway and glycerol formation than on glycolysis and ethanol production. Additionally, nitrogen uptake, both quantitatively and qualitatively, was strongly influenced by environmental conditions. This work provides the most complete stoichiometric model used for Metabolic Flux Analysis of S. cerevisiae in wine fermentations employed so far, including the synthesis and release of relevant aroma compounds and could be used in the design of optimal nitrogen supplementation of wine fermentations.

Glucose-methanol co-utilization in Pichia pastoris studied by metabolomics and instationary ¹³C flux analysis

Background

Several studies have shown that the utilization of mixed carbon feeds instead of methanol as sole carbon source is beneficial for protein production with the methylotrophic yeast Pichia pastoris. In particular, growth under mixed feed conditions appears to alleviate the metabolic burden related to stress responses triggered by protein overproduction and secretion. Yet, detailed analysis of the metabolome and fluxome under mixed carbon source metabolizing conditions are missing. To obtain a detailed flux distribution of central carbon metabolism, including the pentose phosphate pathway under methanol-glucose conditions, we have applied metabolomics and instationary ¹³C flux analysis in chemostat cultivations.

Results

Instationary ¹³C-based metabolic flux analysis using GC-MS and LC-MS measurements in time allowed for an accurate mapping of metabolic fluxes of glycolysis, pentose phosphate and methanol assimilation pathways. Compared to previous results from NMR-derived stationary state labelling data (proteinogenic amino acids, METAFoR) more fluxes could be determined with higher accuracy. Furthermore, using a thermodynamic metabolic network analysis the metabolite measurements and metabolic flux directions were validated. Notably, the concentration of several metabolites of the upper glycolysis and pentose phosphate pathway increased under glucose-methanol feeding compared to the reference glucose conditions, indicating a shift in the thermodynamic driving forces. Conversely, the extracellular concentrations of all measured metabolites were lower compared with the corresponding exometabolome of glucose-grown P. pastoris cells.

The instationary ¹³C flux analysis resulted in fluxes comparable to previously obtained from NMR datasets of proteinogenic amino acids, but allowed several additional insights. Specifically, i) in vivo metabolic flux estimations were expanded to a larger metabolic network e.g. by including trehalose recycling, which accounted for about 1.5% of the glucose uptake rate; ii) the reversibility of glycolytic/gluconeogenesis, TCA cycle and pentose phosphate pathways reactions was estimated, revealing a significant gluconeogenic flux from the dihydroxyacetone phosphate/glyceraldehydes phosphate pool to glucose-6P. The origin of this finding could be carbon recycling from the methanol assimilatory pathway to the pentose phosphate pool. Additionally, high exchange fluxes of oxaloacetate with aspartate as well as malate indicated amino acid pool buffering and the activity of the malate/Asp shuttle; iii) the ratio of methanol oxidation vs utilization appeared to be lower (54 vs 79% assimilated methanol directly oxidized to CO₂).

Conclusions

In summary, the application of instationary ¹³C-based metabolic flux analysis to P. pastoris provides an experimental framework with improved capabilities to explore the regulation of the carbon and energy metabolism of this yeast, particularly for the case of methanol and multicarbon source metabolism.

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/g_DCW·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.

Metabolic flux profiling of recombinant protein secreting Pichia pastoris growing on glucose:methanol mixtures

Background

The methylotrophic yeast Pichia pastoris has emerged as one of the most promising yeast hosts for the production of heterologous proteins. Mixed feeds of methanol and a multicarbon source instead of methanol as sole carbon source have been shown to improve product productivities and alleviate metabolic burden derived from protein production. Nevertheless, systematic quantitative studies on the relationships between the central metabolism and recombinant protein production in P. pastoris are still rather limited, particularly when growing this yeast on mixed carbon sources, thus hampering future metabolic network engineering strategies for improved protein production.

Results

The metabolic flux distribution in the central metabolism of P. pastoris growing on a mixed feed of glucose and methanol was analyzed by Metabolic Flux Analysis (MFA) using ¹³C-NMR-derived constraints. For this purpose, we defined new flux ratios for methanol assimilation pathways in P. pastoris cells growing on glucose:methanol mixtures. By using this experimental approach, the metabolic burden caused by the overexpression and secretion of a Rhizopus oryzae lipase (Rol) in P. pastoris was further analyzed. This protein has been previously shown to trigger the unfolded protein response in P. pastoris. A series of ¹³C-tracer experiments were performed on aerobic chemostat cultivations with a control and two different Rol producing strains growing at a dilution rate of 0.09 h⁻¹ using a glucose:methanol 80:20 (w/w) mix as carbon source.

The MFA performed in this study reveals a significant redistristribution of carbon fluxes in the central carbon metabolism when comparing the two recombinant strains vs the control strain, reflected in increased glycolytic, TCA cycle and NADH regeneration fluxes, as well as higher methanol dissimilation rates.

Conclusions

Overall, a further ¹³C-based MFA development to characterise the central metabolism of methylotrophic yeasts when growing on mixed methanol:multicarbon sources has been implemented, thus providing a new tool for the investigation of the relationships between central metabolism and protein production. Specifically, the study points at a limited but significant impact of the conformational stress associated to secretion of recombinant proteins on the central metabolism, occurring even at modest production levels.

Development of quantitative metabolomics for Pichia pastoris

Accurate, reliable and reproducible measurement of intracellular metabolite levels has become important for metabolic studies of microbial cell factories. A first critical step for metabolomic studies is the establishment of an adequate quenching and washing protocol, which ensures effective arrest of all metabolic activity and removal of extracellular metabolites, without causing leakage of metabolites from the cells. Five different procedures based on cold methanol quenching and cell separation by filtration were tested for metabolomics of Pichia pastoris regarding methanol content and temperature of the quenching solution as key parameters. Quantitative evaluation of these protocols was carried out through mass balance analysis, based on metabolite measurements in all sample fractions, those are whole broth, quenched and washed cells, culture filtrate and quenching and washing solution. Finally, the optimal method was used to study the time profiles of free amino acid and central carbon metabolism intermediates in glucose-limited chemostat cultures. Acceptable recoveries (>90%) were obtained for all quenching procedures tested. However, quenching at -27°C in 60% v/v methanol performed slightly better in terms of leakage minimization. We could demonstrate that five residence times under glucose limitation are enough to reach stable intracellular metabolite pools. Moreover, when comparing P. pastoris and S. cerevisiae metabolomes, under the same cultivation conditions, similar metabolite fingerprints were found in both yeasts, except for the lower glycolysis, where the levels of these metabolites in P. pastoris suggested an enzymatic capacity limitation in that part of the metabolism. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11306-011-0308-1) contains supplementary material, which is available to authorized users.

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 ¹³C 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.

Macromolecular and elemental composition analysis and extracellular metabolite balances of Pichia pastoris growing at different oxygen levels

BACKGROUND

Analysis of the cell operation at the metabolic level requires collecting data of different types and to determine their confidence level. In addition, the acquired information has to be combined in order to obtain a consistent operational view. In the case of Pichia pastoris, information of its biomass composition at macromolecular and elemental level is scarce particularly when different environmental conditions, such as oxygen availability or, genetic backgrounds (e.g. recombinant protein production vs. non production conditions) are compared.

RESULTS

P. pastoris cells growing in carbon-limited chemostat cultures under different oxygenation conditions (% O2 in the bioreactor inlet gas: 21%, 11% and 8%, corresponding to normoxic, oxygen-limiting and hypoxic conditions, respectively), as well as under recombinant protein (antibody fragment, Fab) producing and non-producing conditions, were analyzed from different points of view. On the one hand, the macromolecular and elemental composition of the biomass was measured using different techniques at the different experimental conditions and proper reconciliation techniques were applied for gross error detection of the measured substrates and products conversion rates. On the other hand, fermentation data was analyzed applying elemental mass balances. This allowed detecting a previously missed by-product secreted under hypoxic conditions, identified as arabinitol (aka. arabitol). After identification of this C5 sugar alcohol as a fermentation by-product, the mass balances of the fermentation experiments were validated.

CONCLUSIONS

After application of a range of analytical and statistical techniques, a consistent view of growth parameters and compositional data of P. pastoris cells growing under different oxygenation conditions was obtained. The obtained data provides a first view of the effects of oxygen limitation on the physiology of this microorganism, while recombinant Fab production seems to have little or no impact at this level of analysis. Furthermore, the results will be highly useful in other complementary quantitative studies of P. pastoris physiology, such as metabolic flux analysis.

The effect of temperature on the proteome of recombinant Pichia pastoris

The impact of environmental factors on the productivity of yeast cells is poorly investigated so far. Therefore, it is a major concern to improve the understanding of cellular physiology of microbial protein production hosts, including the methylotrophic yeast Pichia pastoris. Two-Dimensional Fluorescence Difference Gel electrophoresis and protein identification via mass spectrometry were applied to analyze the impact of cultivation temperature on the physiology of a heterologous protein secreting P. pastoris strain. Furthermore, specific productivity was monitored and fluxes through the central carbon metabolism were calculated. Chemostat culture conditions were applied to assess the adaption to different growth temperatures (20, 25, 30 degrees C) at steady-state conditions. Many important cellular processes, including the central carbon metabolism, stress response and protein folding are affected by changing the growth temperature. A 3-fold increased specific productivity at lower cultivation temperature for an antibody Fab fragment was accompanied by a reduced flux through the TCA-cycle, reduced levels of proteins involved in oxidative stress response and lower cellular levels of molecular chaperones. These data indicate that folding stress is generally decreased at lower cultivation temperatures, enabling more efficient heterologous protein secretion in P. pastoris host cells.

[Primary angiosarcoma of the spleen with liver metastases. Case report and literature review]

Angiosarcomas of the spleen are infrequent tumors that often present with liver metastases. These tumors are usually multiple and have a poor prognosis. Immunohistochemical markers can aid preoperative diagnosis since radiological and ultrasonographic images are fairly nonspecific. The recommended treatment is splenectomy. The benefits of chemotherapy and radiotherapy have not been clearly demonstrated.

The MELISSA pilot plant facility as as integration test-bed for advanced life support systems

The different advances in the Micro Ecological Life Support System Alternative project (MELISSA), fostered and coordinated by the European Space Agency, as well as in other associated technologies, are integrated and demonstrated in the MELISSA Pilot Plant laboratory. During the first period of operation, the definition of the different compartments at an individual basis has been achieved, and the complete facility is being re-designed to face a new period of integration of all these compartments. The final objective is to demonstrate the potentiality of biological systems such as MELISSA as life support systems. The facility will also serve as a test bed to study the robustness and stability of the continuous operation of a complex biological system. This includes testing of the associated instrumentation and control for a safe operation, characterization of the chemical and microbial safety of the system, as well as tracking the genetic stability of the microbial strains used. The new period is envisaged as a contribution to the further development of more complete biological life support systems for long-term manned missions, that should be better defined from the knowledge to be gained from this integration phase. This contribution summarizes the current status of the Pilot Plant and the planned steps for the new period.

MELISSA: a loop of interconnected bioreactors to develop life support in space

The development of a loop of interconnected continuous bioreactors, aimed to provide life support in space, is reported. The complete loop concept consists of four bioreactors and one higher plant compartment. For its realization the continuous and controlled operation of the bioreactors is characterized, up to the pilot scale level, first for each individual reactor, second for the interconnected reactor operation. The results obtained with the two more advanced bioreactors in the Micro Ecological Life Support System Alternative (MELISSA) loop are described more specifically. These reactors consist of a packed-bed reactor working with an immobilized co-culture of Nitrosomonas and Nitrobacter cells, and an external loop gas-lift photobioreactor for the culture of the cyanobacteria Spirulina platensis. Their individual operation for long duration runs has been achieved and characterized, and their interconnected operation at pilot scale is reported.

Scale-up and design of a pilot-plant photobioreactor for the continuous culture of Spirulina platensis

Scale-up of bioreactors has the intrinsic difficulty of establishing a reliable relationship among physical parameters involved in the design of the new bioreactor and the physiology of the cultured cells. This is more critical in those cases where a more complex operation of the bioreactor is needed, such as in photobioreactors. A key issue in the operation of photobioreactors is establishing a quantification for the interaction between external illumination, internal light distribution and cell growth. In this paper an approach to the scale-up of a photobioreactor for the culture of Spirulina platensis, based on a mathematical model describing this interaction, and the operation of a previous reactor 10 times smaller is presented. The paper describes the approach followed in the scale-up, the analysis of different design constraints, the physical realization of the new bioreactor design, innovative use of plastic material walls to improve reactor safety, and finally the corroboration of its satisfactory operation.

Modeling photoheterotrophic growth kinetics of Rhodospirillum rubrum inrectangular photobioreactors

Based on a previously established model for radiant light transfer in photobioreactors (PBR), taking into account absorption and scattering of light, a new knowledge model for coupling radiant light energy available and local growth kinetics in PBRs for the photoheterotrophic bacteria Rhodospirillum rubrum is discussed. A revised method is presented for the calculation of the absorption and scattering coefficients. The specific characteristics of the electron-transfer chains in such microorganisms leads to definition of three different metabolic zones in the PBR, explaining the behavior of mean kinetics observed in a wide range of incident light fluxes. The model is validated in rectangular PBRs for five different carbon sources and proved robust and fully predictive. This approach can be considered for simulation and model-based predictive control of PBRs cultivating photoheterotrophic microorganisms.