Metabolic flux analysis for recombinant protein production by Pichia pastoris using dual carbon sources: Effects of methanol feeding rate

Bioconversion of C1 feedstocks for chemical production using Pichia pastoris

Bioconversion of C1 feedstocks for chemical production offers a promising solution to global challenges such as the energy and food crises and climate change. The methylotroph Pichia pastoris is an attractive host system for the production of both recombinant proteins and chemicals from methanol. Recent studies have also demonstrated its potential for utilizing CO2 through metabolic engineering or coupling with electrocatalysis. This review focuses on the bioconversion of C1 feedstocks for chemical production using P. pastoris. Herein the challenges and feasible strategies for chemical production in P. pastoris are discussed. The potential of P. pastoris to utilize other C1 feedstocks - including CO2 and formate - is highlighted, and new insights from the perspectives of synthetic biology and material science are proposed.

Developing a metabolic model-based fed-batch feeding strategy for Pichia pastoris fermentation through fine-tuning of the methanol utilization pathway

Pichia pastoris is a commonly used microbial host for recombinant protein production. It is mostly cultivated in fed-batch mode, in which the environment of the cell is continuously changing. Hence, it is vital to understand the influence of feeding strategy parameters on the intracellular reaction network to fine-tune bioreactor performance. This study used dynamic flux balance analysis (DFBA) integrated with transcriptomics data to simulate the recombinant P. pastoris (Mut^s ) growth during the induction phase for three fed-batch strategies, conducted at constant specific growth rates (μ-stat). The induction phase was split into equal time intervals, and the correlated reactions with protein yield were identified in the three fed-batch strategies using the Pearson correlation coefficient. Subsequently, principal component analysis (PCA) was applied to cluster induction phase time intervals and identify the role of correlated reactions on metabolic differentiation of time intervals. It was found that increasing fluxes through the methanol dissimilation pathway increased protein yield. By adding a methanol assimilation pathway inhibitor (HgCl₂ ) to the shake flask medium growing on glycerol: methanol mixture (10%: 90%, v/v), the protein titre increased by 60%. As per DFBA, the higher the methanol to biomass flux ratio (R_meoh/Δx ), the higher the protein yield. Finally, a novel feeding strategy was developed to increase the amount of R_meoh/Δx compared to the three feeding strategies. The concentration of recombinant human growth hormone (rhGH), used as the model protein, increased by 16% compared to the optimal culture result obtained previously (800 mg L^-1 to 928 mg L^-1 ), while production yield improved by 85% (24.8 mg g_DCW ^-1 to 46 mg g_DCW ^-1 ).

Role of Dissimilative Pathway of Komagataella phaffii (Pichia pastoris): Formaldehyde Toxicity and Energy Metabolism

Komagataella phaffii (aka Pichia pastoris) is a yeast able to grow in methanol as the sole carbon and energy source. This substrate is converted into formaldehyde, a toxic intermediary that can either be assimilated to biomass or dissimilated to CO₂ through the enzymes formaldehyde dehydrogenase (FLD) and formate dehydrogenase, also producing energy in the form of NADH. The dissimilative pathway has been described as an energy producing and a detoxifying route, but conclusive evidence has not been provided for this. In order to elucidate this theory, we generated mutants lacking the FLD activity (Δfld1) and used flux analysis to evaluate the metabolic impact of this disrupted pathway. Unexpectedly, we found that the specific growth rate of the Δfld1 strain was only slightly lower (92%) than the control. In contrast, the sensitivity to formaldehyde pulses (up to 8mM) was significantly higher in the Δfld1 mutant strain and was associated with a higher maintenance energy. In addition, the intracellular flux estimation revealed a high metabolic flexibility of K. phaffii in response to the disrupted pathway. Our results suggest that the role of the dissimilative pathway is mainly to protect the cells from the harmful effect of formaldehyde, as they were able to compensate for the energy provided from this pathway when disrupted.

Comparative transcriptome and metabolome analyses reveal the methanol dissimilation pathway of Pichia pastoris

BACKGROUND

Pichia pastoris (Komagataella phaffii) is a model organism widely used for the recombinant expression of eukaryotic proteins, and it can metabolize methanol as its sole carbon and energy source. Methanol is oxidized to formaldehyde by alcohol oxidase (AOX). In the dissimilation pathway, formaldehyde is oxidized to CO₂ by formaldehyde dehydrogenase (FLD), S-hydroxymethyl glutathione hydrolase (FGH) and formate dehydrogenase (FDH).

RESULTS

The transcriptome and metabolome of P. pastoris were determined under methanol cultivation when its dissimilation pathway cut off. Firstly, Δfld and Δfgh were significantly different compared to the wild type (GS115), with a 60.98% and 23.66% reduction in biomass, respectively. The differential metabolites between GS115 and Δfld were mainly enriched in ABC transporters, amino acid biosynthesis, and protein digestion and absorption. Secondly, comparative transcriptome between knockout and wild type strains showed that oxidative phosphorylation, glycolysis and the TCA cycle were downregulated, while alcohol metabolism, proteasomes, autophagy and peroxisomes were upregulated. Interestingly, the down-regulation of the oxidative phosphorylation pathway was positively correlated with the gene order of dissimilation pathway knockdown. In addition, there were significant differences in amino acid metabolism and glutathione redox cycling that raised our concerns about formaldehyde sorption in cells.

CONCLUSIONS

This is the first time that integrity of dissimilation pathway analysis based on transcriptomics and metabolomics was carried out in Pichia pastoris. The blockage of dissimilation pathway significantly down-regulates the level of oxidative phosphorylation and weakens the methanol assimilation pathway to the point where deficiencies in energy supply and carbon fixation result in inefficient biomass accumulation and genetic replication. In addition, transcriptional upregulation of the proteasome and autophagy may be a stress response to resolve formaldehyde-induced DNA-protein crosslinking.

Metabolic engineering of Pichia pastoris for malic acid production from methanol

The application of rational design in reallocating metabolic flux to accumulate desired chemicals is always restricted by the native regulatory network. In this study, recombinant Pichia pastoris was constructed for malic acid production from sole methanol through rational redistribution of metabolic flux. Different malic acid accumulation modules were systematically evaluated and optimized in P. pastoris. The recombinant PP-CM301 could produce 8.55 g/L malic acid from glucose, which showed a 3.45-fold increase compared to the parent strain. To improve the efficiency of site-directed gene knockout, NHEJ-related protein Ku70 was destroyed, whereas leading to the silencing of heterogenous genes. Hence, genes related to by-product generation were deleted via a specially designed FRT/FLP system, which successfully reduced succinic acid and ethanol production. Furthermore, a key node in the methanol assimilation pathway, glucose-6-phosphate isomerase was knocked out to liberate metabolic fluxes trapped in the XuMP cycle, which finally enabled 2.79 g/L malic acid accumulation from sole methanol feeding with nitrogen source optimization. These results will provide guidance and reference for the metabolic engineering of P. pastoris to produce value-added chemicals from methanol.

Combining metabolic flux analysis and adaptive evolution to enhance lipase production in Bacillus subtilis

Metabolic fluxes during lipase production by Bacillus subtilis CICC 20034 in synthetic medium were studied using metabolic flux analysis (MFA). The MFA showed that lipase production was dependent on, and coupled to the tributyrin uptake rate, formation of biomass, lactate, ATP, as well as amino acids from the aspartate and glutamate family. Using tributyrin as the sole carbon source, an adaptive evolution strategy was applied to increase the tributyrin uptake rate. B. subtilis SPZ1 was obtained from CICC 20034 by adaptive evolution over 1000 generations of growth-based selection. The tributyrin consumption rate of strain SPZ1 reached 0.89 g/(L·h) which was 1.9-fold higher than that of the original strain. The MFA indicated that the 212% increase of tributyrin uptake flux contributed to the 556% increase of lipase flux. Consequently, the lipase activity (0.65 U/mL) of strain SPZ1 was 1.9-fold higher than that of the original strain. This was the highest lipase activity obtained by fermentation in synthetic medium reported for Bacillus strains. In complex culture medium, lipase activity of SPZ1 reached 3.3 U/mL.

Yeast Expression Systems: Overview and Recent Advances

Yeasts are outstanding hosts for the production of functional recombinant proteins with industrial or medical applications. Great attention has been emerged on yeast due to the inherent advantages and new developments in this host cell. For the production of each specific product, the most appropriate expression system should be identified and optimized both on the genetic and fermentation levels, considering the features of the host, vector and expression strategies. Currently, several new systems are commercially available; some of them are private and need licensing. The potential for secretory expression of heterologous proteins in yeast proposed this system as a candidate for the production of complex eukaryotic proteins. The common yeast expression hosts used for recombinant proteins' expression include Saccharomyces cerevisiae, Pichia pastoris, Hansenula polymorpha, Yarrowia lipolytica, Arxula adeninivorans, Kluyveromyces lactis, and Schizosaccharomyces pombe. This review is dedicated to discuss on significant characteristics of the most common methylotrophic and non-methylotrophic yeast expression systems with an emphasis on their advantages and new developments.

High Expression of Human Cathepsin S by Recombinant Pichia pastoris with Cod Skin as an Organic Co-Nitrogen Source

Human cathepsin S production by recombinant Pichia pastoris using cod skin as the co-nitrogen source was investigated in this study. The addition of carbon sources of glycerol in the fed-batch phase and of methanol in the induction stage was also investigated. A new approach to the highly expression of human cathepsin S was developed using 90 g/L of cod skin (wet weight). After 24 h of the initial fermentation, 4% glycerol (v/v, glycerol/culture) was added once to enhance the cell density (OD600) in the cultivation. Then, adding and maintaining methanol at 0.5% (v/v, methanol/cultivation) after about 48 h of fermentation achieved a high expression of human cathepsin S in a 5-L bioreactor. The results demonstrate that the maximum activity of human cathepsin S in the fermentation supernatant reached 7,152 U/L after 96 h of methanol induction. The methylotrophic yeast P. pastoris grown in the medium containing cod skin (90 g/L) as the co-nitrogen source provided a 21% higher cell density (OD600) and 18.3% higher human cathepsin S yield than P. pastoris grown in BMGY medium. For the first time, human cathepsin S was successfully expressed by P. pastoris with cod skin as the co-nitrogen source. The glycerol fed-batch controlling strategy and method of maintaining methanol at a constant concentration of 0.5% (v/v, methanol/cultivation) in the induction stage was efficient for P. pastoris growth and the expression of human cathepsin S.

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.

Engineering strategies for enhanced production of protein and bio-products in Pichia pastoris: A review

Pichia pastoris has been recognized as one of the most industrially important hosts for heterologous protein production. Despite its high protein productivity, the optimization of P. pastoris cultivation is still imperative due to strain- and product-specific challenges such as promoter strength, methanol utilization type and oxygen demand. To address the issues, strategies involving genetic and process engineering have been employed. Optimization of codon usage and gene dosage, as well as engineering of promoters, protein secretion pathways and methanol metabolic pathways have proved beneficial to innate protein expression levels. Large-scale production of proteins via high cell density fermentation additionally relies on the optimization of process parameters including methanol feed rate, induction temperature and specific growth rate. Recent progress related to the enhanced production of proteins in P. pastoris via various genetic engineering and cultivation strategies are reviewed. Insight into the regulation of the P. pastoris alcohol oxidase 1 (AOX1) promoter and the development of methanol-free systems are highlighted. Novel cultivation strategies such as mixed substrate feeding are discussed. Recent advances regarding substrate and product monitoring techniques are also summarized. Application of P. pastoris to the production of biodiesel and other value-added products via metabolic engineering are also reviewed. P. pastoris is becoming an indispensable platform through the use of these combined engineering strategies.

Enhancing monellin production by Pichia pastoris at low cell induction concentration via effectively regulating methanol metabolism patterns and energy utilization efficiency

In heterologous protein productions by P. pastoris, methanol induction is generally initiated when cell concentration reaches very high density. The alternative strategy by initiating methanol induction at lower cells concentration was also reported to be effective in easing DO control, reducing toxic by-metabolites accumulation and increasing targeted proteins titers. However, the methanol/energy regulation mechanisms are seldom reported. We theoretically analyzed the methanol/energy metabolisms in protein expression process with the strategies of initiating induction at higher or lower cells concentrations, using monellin production as a prototype. When initiating induction at lower cells concentration and controlling induction temperature at 30°C, monellin concentration reached the highest levels of 2.62~2.71 g/L, which was 2.5~4.9 fold of those obtained with the strategy of initiating induction at higher cells concentration. With the desired induction strategy, 1) carbon metabolism ratio directing into the precursors synthesis route for monellin production reached the highest level of 65%, carbon metabolism ratios towards to precursors synthesis and ATP regeneration routes were regulated at relatively balanced levels; 2) monellin synthesis was completely cell growth associated, with the largest associated coefficient and higher specific growth rate; 3) theoretical NADH (energy) utilization efficiency η was the highest, and η stayed high levels (≥0.8) during most period (89%) within induction phase to supply sufficient energy in supporting monellin synthesis.

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

Investigating the effect of carbon source on rabies virus glycoprotein production in Pichia pastoris by a transcriptomic approach

Several factors affect protein expression in Pichia pastoris, one among them is the carbon source. In this work, we studied the effect of this factor on the expression level of rabies virus glycoprotein (RABV-G) in two recombinant clones harboring seven copies of the gene of interest. The expression was driven either by the constitutive glyceraldehyde-3-phosphate dehydrogenase (GAP) promoter or the inducible alcohol oxidase1 (AOX1) promoter. Clones were compared in terms of cell physiology and carbon source metabolism. The transcription levels of 16 key genes involved in the central metabolic pathway, the methanol catabolism, and the oxidative stress were investigated in both clones. Cell size, as a parameter reflecting cell physiological changes, was also monitored. Our results showed that when glucose was used as the sole carbon source, large cells were obtained. Transcript levels of the genes of the central metabolic pathway were also upregulated, whereas antioxidative gene transcript levels were low. By contrast, the use of methanol as a carbon source generated small cells and a shift in carbon metabolism toward the dissimilatory pathway by the upregulation of formaldehyde dehydrogenase gene and the downregulation of those of the central metabolic. These observations are in favor of the use of glucose to enhance the expression of RABV-G in P. pastoris.

Cultivation of Pichia pastoris carrying the scFv anti LDL (-) antibody fragment. Effect of preculture carbon source

Antibodies and antibody fragments are nowadays among the most important biotechnological products, and Pichia pastoris is one of the most important vectors to produce them as well as other recombinant proteins. The conditions to effectively cultivate a P. pastoris strain previously genetically modified to produce the single-chain variable fragment anti low density lipoprotein (−) under the control of the alcohol oxidase promoter have been investigated in this study. In particular, it was evaluated if, and eventually how, the carbon source (glucose or glycerol) used in the preculture preceding cryopreservation in 20% glycerol influences both cell and antibody fragment productions either in flasks or in bioreactor. Although in flasks the volumetric productivity of the antibody fragment secreted by cells precultured, cryopreserved and reactivated in glycerol was 42.9% higher compared with cells precultured in glucose, the use of glycerol in bioreactor led to a remarkable shortening of the lag phase, thereby increasing it by no less than thrice compared to flasks. These results are quite promising in comparison with those reported in the literature for possible future industrial applications of this cultivation, taking into account that the overall process time was reduced by around 8 h.

Overcoming nutrient limitations for cell-based production of influenza vaccine

Metabolic analysis for medium optimization represents a very useful strategy in the process development of production of vaccines in cells. During influenza vaccine production, viruses hijack host cells and take advantage of host's metabolism. As a consequence, the nutritional demand of host cells should undergo a profound change, and usually more nutrients such as glucose and amino acids should be consumed. As such, the maintaining media used in virus production processes often cannot provide sufficient nutrients, and novel methods are urged to be established to address this severe issue of nutritional limitation. A detailed study on impacts of influenza virus on cell death and metabolism, with a profound analysis of nutritional requirements during virus production process, followed by a rational medium optimization is expected to be the most straightfoward and effective strategy. This would ensure a balanced and adequate nutritional supply, which should minimize cell death and improve both cell-specific virus yield and total influenza virus production. Such a metabolic analysis-based medium optimization would lay a solid foundation for the development of cell culture technology in influenza vaccine production.

Influence of methanol/sorbitol co-feeding rate on pAOX1 induction in a Pichia pastoris Mut+ strain in bioreactor with limited oxygen transfer rate

High Pichia pastoris biomass density could be obtained using high co-feeding rate of methanol and sorbitol in a fed-batch or continuous culture, while further higher feeding rate finally leads to oxygen limitation in bioreactor. In the literature, there is lack of report about AOX1 promoter regulation with regard to dissolved oxygen level (DO). Therefore, in this work, chemostat cultures were performed to investigate the cell growth, metabolism and regulation of the AOX1 promoter (pAOX1) regarding co-feeding rate of optimized methanol/sorbitol mixture (methanol fraction 0.60 C-mol/C-mol) using a P. pastoris Mut+/pAOX1-lacZ strain. The oxygen transfer rates (OTR) in bioreactor were kept in the range of typical values of large bioreactor, i.e., 4–8 g/(L h) if DO equals 30 % saturation or 5–10 g/(L h) if DO nears zero. For DO >0, an increase of the carbon fed led to an increase of pAOX1 induction. By contrast, when dissolved oxygen was completely depleted, methanol accumulated, causing a 30 % decrease of pAOX1 induction. However, this decrease is more likely to be lined to methanol accumulation than to low level of dissolved oxygen (<4 % DO). Methanol/sorbitol co-feeding allowed cells to adapt to oxygen transient limitations that often occur at industrial scale with reduced effect on pAOX1 induction. The optimal feeding rate tested here was 6.6 mmol C (DCW h)−1 at an OTR of 8.28 g O2(L h)−1 with over fivefold pAOX1 induction (probably directly associated with target protein productivity) compared with previous work.

Optimization of lipid production with a genome-scale model of Yarrowia lipolytica

Background

Yarrowia lipolytica is a non-conventional yeast that is extensively investigated for its ability to excrete citrate or to accumulate large amounts of storage lipids, which is of great significance for single cell oil production. Both traits are thus of interest for basic research as well as for biotechnological applications but they typically occur simultaneously thus lowering the respective yields. Therefore, engineering of strains with high lipid content relies on novel concepts such as computational simulation to better understand the two competing processes and to eliminate citrate excretion.

Results

Using a genome-scale model (GSM) of baker's yeast as a scaffold, we reconstructed the metabolic network of Y. lipolytica and optimized it for use in flux balance analysis (FBA), with the aim to simulate growth and lipid production phases of this yeast. We validated our model and found the predictions of the growth behavior of Y. lipolytica in excellent agreement with experimental data. Based on these data, we successfully designed a fed-batch strategy to avoid citrate excretion during the lipid production phase. Further analysis of the network suggested that the oxygen demand of Y. lipolytica is reduced upon induction of lipid synthesis. According to this finding we hypothesized that a reduced aeration rate might induce lipid accumulation. This prediction was indeed confirmed experimentally. In a fermentation combining these two strategies lipid content of the biomass was increased by 80 %, and lipid yield was improved more than four-fold, compared to standard conditions.

Conclusions

Genome scale network reconstructions provide a powerful tool to predict the effects of genetic modifications and the metabolic response to environmental conditions. The high accuracy and the predictive value of a newly reconstructed GSM of Y. lipolytica to optimize growth conditions for lipid accumulation are demonstrated. Based on these findings, further strategies for engineering Y. lipolytica towards higher efficiency in single cell oil production are discussed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-015-0217-4) contains supplementary material, which is available to authorized users.

Functional recombinant protein is present in the pre-induction phases of Pichia pastoris cultures when grown in bioreactors, but not shake-flasks

Background

Pichia pastoris is a widely-used host for recombinant protein production; expression is typically driven by methanol-inducible alcohol oxidase (AOX) promoters. Recently this system has become an important source of recombinant G protein-coupled receptors (GPCRs) for structural biology and drug discovery. The influence of diverse culture parameters (such as pH, dissolved oxygen concentration, medium composition, antifoam concentration and culture temperature) on productivity has been investigated for a wide range of recombinant proteins in P. pastoris. In contrast, the impact of the pre-induction phases on yield has not been as closely studied. In this study, we examined the pre-induction phases of P. pastoris bioreactor cultivations producing three different recombinant proteins: the GPCR, human A_2a adenosine receptor (hA_2aR), green fluorescent protein (GFP) and human calcitonin gene-related peptide receptor component protein (as a GFP fusion protein; hCGRP-RCP-GFP).

Results

Functional hA_2aR was detected in the pre-induction phases of a 1 L bioreactor cultivation of glycerol-grown P. pastoris. In a separate experiment, a glycerol-grown P. pastoris strain secreted soluble GFP prior to methanol addition. When glucose, which has been shown to repress AOX expression, was the pre-induction carbon source, hA_2aR and GFP were still produced in the pre-induction phases. Both hA_2aR and GFP were also produced in methanol-free cultivations; functional protein yields were maintained or increased after depletion of the carbon source. Analysis of the pre-induction phases of 10 L pilot scale cultivations also demonstrated that pre-induction yields were at least maintained after methanol induction, even in the presence of cytotoxic concentrations of methanol. Additional bioreactor data for hCGRP-RCP-GFP and shake-flask data for GFP, horseradish peroxidase (HRP), the human tetraspanins hCD81 and CD82, and the tight-junction protein human claudin-1, demonstrated that bioreactor but not shake-flask cultivations exhibit recombinant protein production in the pre-induction phases of P. pastoris cultures.

Conclusions

The production of recombinant hA_2aR, GFP and hCGRP-RCP-GFP can be detected in bioreactor cultivations prior to methanol induction, while this is not the case for shake-flask cultivations of GFP, HRP, hCD81, hCD82 and human claudin-1. This confirms earlier suggestions of leaky expression from AOX promoters, which we report here for both glycerol- and glucose-grown cells in bioreactor cultivations. These findings suggest that the productivity of AOX-dependent bioprocesses is not solely dependent on induction by methanol. We conclude that in order to maximize total yields, pre-induction phase cultivation conditions should be optimized, and that increased specific productivity may result in decreased biomass yields.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-014-0127-y) contains supplementary material, which is available to authorized users.

Improving functional annotation for industrial microbes: a case study with Pichia pastoris

•

The current status of the Pichia pastoris genome is shown to lack extensive functional annotation.

•

GO annotation transfer and literature curation pipelines improve the functional annotation of genomes.

•

Pipelines and tools that can improve the annotation status of the genomes of Pichia pastoris and many industrial microbes are considered.

•

Well-annotated genome sequences will facilitate the utilization of these microbes in a broader range of synthetic biology applications.

The research communities studying microbial model organisms, such as Escherichia coli or Saccharomyces cerevisiae, are well served by model organism databases that have extensive functional annotation. However, this is not true of many industrial microbes that are used widely in biotechnology. In this Opinion piece, we use Pichia (Komagataella) pastoris to illustrate the limitations of the available annotation. We consider the resources that can be implemented in the short term both to improve Gene Ontology (GO) annotation coverage based on annotation transfer, and to establish curation pipelines for the literature corpus of this organism.

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

A quantitative study of methanol/sorbitol co-feeding process of a Pichia pastoris Mut⁺/pAOX1-lacZ strain

Background

One of the main challenges for heterologous protein production by the methylotrophic yeast Pichia pastoris at large-scale is related to its high oxygen demand. A promising solution is a co-feeding strategy based on a methanol/sorbitol mixture during the induction phase. Nonetheless, a deep understanding of the cellular physiology and the regulation of the AOX1 promoter, used to govern heterologous protein production, during this co-feeding strategy is still scarce.

Results

Transient continuous cultures with a dilution rate of 0.023 h^-1 at 25°C were performed to quantitatively assess the benefits of a methanol/sorbitol co-feeding process with a Mut⁺ strain in which the pAOX1-lacZ construct served as a reporter gene. Cell growth and metabolism, including O₂ consumption together with CO₂ and heat production were analyzed with regard to a linear change of methanol fraction in the mixed feeding media. In addition, the regulation of the promoter AOX1 was investigated by means of β-galactosidase measurements. Our results demonstrated that the cell-specific oxygen consumption (qO₂) could be reduced by decreasing the methanol fraction in the feeding media. More interestingly, maximal β-galactosidase cell-specific activity (>7500 Miller unit) and thus, optimal pAOX1 induction, was achieved and maintained in the range of 0.45 ~ 0.75 C-mol/C-mol of methanol fraction. In addition, the qO₂ was reduced by 30% at most in those conditions. Based on a simplified metabolic network, metabolic flux analysis (MFA) was performed to quantify intracellular metabolic flux distributions during the transient continuous cultures, which further shed light on the advantages of methanol/sorbitol co-feeding process. Finally, our observations were further validated in fed-batch cultures.

Conclusion

This study brings quantitative insight into the co-feeding process, which provides valuable data for the control of methanol/sorbitol co-feeding, aiming at enhancing biomass and heterologous protein productivities under given oxygen supply. According to our results, β-galactosidase productivity could be improved about 40% using the optimally mixed feed.

Metabolic reconstruction and flux analysis of industrial Pichia yeasts

Pichia yeasts have been recognized as important microbial cell factories in the biotechnological industry. Notably, the Pichia pastoris and Pichia stipitis species have attracted much research interest due to their unique cellular physiology and metabolic capability: P. pastoris has the ability to utilize methanol for cell growth and recombinant protein production, while P. stipitis is capable of assimilating xylose to produce ethanol under oxygen-limited conditions. To harness these characteristics for biotechnological applications, it is highly required to characterize their metabolic behavior. Recently, following the genome sequencing of these two Pichia species, genome-scale metabolic networks have been reconstructed to model the yeasts' metabolism from a systems perspective. To date, there are three genome-scale models available for each of P. pastoris and P. stipitis. In this mini-review, we provide an overview of the models, discuss certain limitations of previous studies, and propose potential future works that can be conducted to better understand and engineer Pichia yeasts for industrial applications.

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.

In silico aided metabolic engineering of Streptomyces roseosporus for daptomycin yield improvement

In silico metabolic network models are valuable tools for strain improvement with desired properties. In this work, based on the comparisons of each pathway flux under two different objective functions for the reconstructed metabolic network of Streptomyces roseosporus, three potential targets of zwf2 (code for glucose-6-phosphate hydrogenase), dptI (code for α-ketoglutarate methyltransferase), and dptJ (code for tryptophan oxygenase) were identified and selected for the genetic modifications. Overexpression of zwf2, dptI, and dptJ genes increased the daptomycin concentration up to 473.2, 452.5, and 489.1 mg/L, respectively. Furthermore, co-overexpression of three genes in series resulted in a 34.4% higher daptomycin concentration compared with the parental strain, which ascribed to the synergistic effect of the enzymes responsible for daptomycin biosynthesis. Finally, the engineered strain enhanced the yield of daptomycin up to 581.5 mg/L in the fed-batch culture, which was approximately 43.2% higher than that of the parental strain. These results demonstrated that the metabolic network based on in silico prediction would be accurate, reasonable, and practical for target gene identification and strain improvement.

Metabolic flux analysis and principal nodes identification for daptomycin production improvement by Streptomyces roseosporus

In the present work, a comprehensive metabolic network of Streptomyces roseosporus LC-54-20 was proposed for daptomycin production. The analysis of extracellular metabolites throughout the batch fermentation was evaluated in addition to daptomycin and biomass production. Metabolic flux distributions were based on stoichiometrical reaction as well as the extracellular metabolites fluxes. Experimental and calculated values for both the specific growth rate and daptomycin production rate indicated that the in silico model proved a powerful tool to analyze the metabolic behaviors based on the analysis under different initial glucose concentrations throughout the fermentation. Through manipulating different pH values, the production rates of various extracellular metabolites were also presented in this paper. Flux distribution variations revealed that the daptomycin production could be significantly influenced by the branch points of glucose 6-phosphate, 3-phosphoglycerate, phosphoenolpyruvate, pyruvate, and oxaloacetate. The five principal metabolites were certified as the flexible nodes and could form potential bottlenecks for a further enhancement of daptomycin production. Furthermore, various concentrations of the five precursors were added into the batch fermentation and led to the enhancement of daptomycin concentration and production rate.

Quantification of metabolic limitations during recombinant protein production in Escherichia coli

Escherichia coli is one of the major microorganisms for recombinant protein production because it has been best characterized in terms of molecular genetics and physiology, and because of the availability of various expression vectors and strains. The synthesis of proteins is one of the most energy consuming processes in the cell, with the result that cellular energy supply may become critical. Indeed, the so called metabolic burden of recombinant protein synthesis was reported to cause alterations in the operation of the host's central carbon metabolism. To quantify these alterations in E. coli metabolism in dependence of the rate of recombinant protein production, (13)C-tracer-based metabolic flux analysis in differently induced cultures was used. To avoid dilution of the (13)C-tracer signal by the culture history, the recombinant protein produced was used as a flux probe, i.e., as a read out of intracellular flux distributions. In detail, an increase in the generation rate rising from 36 mmol(ATP)g(CDW)(-1)h(-1) for the reference strain to 45 mmol(ATP)g(CDW)(-1)h(-1) for the highest yielding strain was observed during batch cultivation. Notably, the flux through the TCA cycle was rather constant at 2.5±0.1 mmol g(CDW)(-1)h(-1), hence was independent of the induced strength for gene expression. E. coli compensated for the additional energy demand of recombinant protein synthesis by reducing the biomass formation to almost 60%, resulting in excess NADPH. Speculative, this excess NADPH was converted to NADH via the soluble transhydrogenase and subsequently used for ATP generation in the electron transport chain. In this study, the metabolic burden was quantified by the biomass yield on ATP, which constantly decreased from 11.7g(CDW)mmol(ATP)(-1) for the reference strain to 4.9g(CDW)mmol(ATP)(-1) for the highest yielding strain. The insights into the operation of the metabolism of E. coli during recombinant protein production might guide the optimization of microbial hosts and fermentation conditions.

Carbon metabolism limits recombinant protein production in Pichia pastoris

The yeast Pichia pastoris enables efficient (high titer) recombinant protein production. As the molecular tools required are well established and gene specific optimizations of transcription and translation are becoming available, metabolism moves into focus as possible limiting factor of recombinant protein production in P. pastoris. To investigate the impact of recombinant protein production on metabolism systematically, we constructed strains that produced the model protein β-aminopeptidase BapA of Sphingosinicella xenopeptidilytica at different production yields. The impact of low to high BapA production on cell physiology was quantified. The data suggest that P. pastoris compensates for the additional resources required for recombinant protein synthesis by reducing by-product formation and by increasing energy generation via the TCA cycle. Notably, the activity of the TCA cycle was constant with a rate of 2.1 ± 0.1 mmol g CDW-1 h(-1) irrespective of significantly reduced growth rates in high BapA producing strains, suggesting an upper limit of TCA cycle activity. The reduced growth rate could partially be restored by providing all 20 proteinogenic amino acids in the fermentation medium. Under these conditions, the rate of BapA synthesis increased twofold. The successful supplementation of the growth medium by amino acids to unburden cellular metabolism during recombinant protein production suggests that the metabolic network is a valid target for future optimization of protein production by P. pastoris.

Genome-scale metabolic reconstruction and in silico analysis of methylotrophic yeast Pichia pastoris for strain improvement

Background

Pichia pastoris has been recognized as an effective host for recombinant protein production. A number of studies have been reported for improving this expression system. However, its physiology and cellular metabolism still remained largely uncharacterized. Thus, it is highly desirable to establish a systems biotechnological framework, in which a comprehensive in silico model of P. pastoris can be employed together with high throughput experimental data analysis, for better understanding of the methylotrophic yeast's metabolism.

Results

A fully compartmentalized metabolic model of P. pastoris (iPP668), composed of 1,361 reactions and 1,177 metabolites, was reconstructed based on its genome annotation and biochemical information. The constraints-based flux analysis was then used to predict achievable growth rate which is consistent with the cellular phenotype of P. pastoris observed during chemostat experiments. Subsequent in silico analysis further explored the effect of various carbon sources on cell growth, revealing sorbitol as a promising candidate for culturing recombinant P. pastoris strains producing heterologous proteins. Interestingly, methanol consumption yields a high regeneration rate of reducing equivalents which is substantial for the synthesis of valuable pharmaceutical precursors. Hence, as a case study, we examined the applicability of P. pastoris system to whole-cell biotransformation and also identified relevant metabolic engineering targets that have been experimentally verified.

Conclusion

The genome-scale metabolic model characterizes the cellular physiology of P. pastoris, thus allowing us to gain valuable insights into the metabolism of methylotrophic yeast and devise possible strategies for strain improvement through in silico simulations. This computational approach, combined with synthetic biology techniques, potentially forms a basis for rational analysis and design of P. pastoris metabolic network to enhance humanized glycoprotein production.

Metabolic flux analysis in eukaryotes

Metabolic flux analysis (MFA) represents a powerful tool for systems biology research on eukaryotic cells. This review describes recent advances, the challenges as well as applications of metabolic flux analysis comprising fungi, mammalian cells and plants. While MFA is widely established and applied in microorganisms, it remains still a challenge to adapt these methods to eukaryotic cell systems having a higher complexity particularly concerning compartmentation or media composition. In fungi MFA was used in the past few years to analyze a variety of conditions and factors and their effects on cellular metabolism. In mammalian cells MFA was applied mainly in cell culture technology and in medical and toxicological research. (13)C metabolic studies on native whole plants are additionally challenging by the fact that CO(2) is usually the only carbon source.