Metabolic flux analysis of recombinant Pichia pastoris growing on different glycerol/methanol mixtures by iterative fitting of NMR-derived 13C-labelling data from proteinogenic amino acids

Continuous Process Verification 4.0 application in upstream: adaptiveness implementation managed by AI in the hypoxic bioprocess of the Pichia pastoris cell factory

The experimental approach developed in this research demonstrated how the cloud, the Internet of Things (IoT), edge computing, and Artificial Intelligence (AI), considered key technologies in Industry 4.0, provide the expected horizon for adaptive vision in Continued Process Verification (CPV), the final stage of Process Validation (PV). Pichia pastoris producing Candida rugosa lipase 1 under the regulation of the constitutive GAP promoter was selected as an experimental bioprocess. The bioprocess worked under hypoxic conditions in carbon-limited fed-batch cultures through a physiological control based on the respiratory quotient (RQ). In this novel bioprocess, a digital twin (DT) was built and successfully tested. The implementation of online sensors worked as a bridge between the microorganism and AI models, to provide predictions from the edge and the cloud. AI models emulated the metabolism of Pichia based on critical process parameters and actionable factors to achieve the expected quality attributes. This innovative AI-aided Adaptive-Proportional Control strategy (AI-APC) improved the reproducibility comparing to a Manual-Heuristic Control strategy (MHC), showing better performance than the Boolean-Logic-Controller (BLC) tested. The accuracy, indicated by the Mean Relative Error (MRE), was for the AI-APC lower than 4%, better than the obtained for MHC (10%) and BLC (5%). Moreover, in terms of precision, the same trend was observed when comparing the Root Mean Square Deviation (RMSD) values, becoming lower as the complexity of the controller increases. The successful automatic real time control of the bioprocess orchestrated by AI models proved the 4.0 capabilities brought by the adaptive concept and its validity in biopharmaceutical upstream operations.

Candida boidinii isolates from olive curation water: a promising platform for methanol-based biomanufacturing

Methanol is a promising feedstock for biomanufacturing, but the efficiency of methanol-based bioprocesses is limited by the low rate of methanol utilization pathways and methanol toxicity. Yeast diversity is an attractive biological resource to develop efficient bioprocesses since any effort with strain improvement is more deserving if applied to innate robust strains with relevant catabolic and biosynthetic potential. The present study is in line with such rational and describes the isolation and molecular identification of seven isolates of the methylotrophic species Candida boidinii from waters derived from the traditional curation of olives, in different years, and from contaminated superficial soil near fuel stations. The yeast microbiota from those habitats was also characterized. The four C. boidinii isolates obtained from the curation of olives' water exhibited significantly higher maximum specific growth rates (range 0.15-0.19 h-1), compared with the three isolates obtained from the fuel contaminated soils (range 0.05-0.06 h-1) when grown on methanol as the sole C-source (1% (v/v), in shake flasks, at 30°C). The isolates exhibit significant robustness towards methanol toxicity that increases as the cultivation temperature decreases from 30°C to 25°C. The better methanol-based growth performance exhibited by C. boidinii isolates from olives´ soaking waters could not be essentially attributed to higher methanol tolerance. These methanol-efficient catabolizing isolates are proposed as a promising platform to develop methanol-based bioprocesses.

Scalable protein production by Komagataella phaffii enabled by ARS plasmids and carbon source-based selection

BACKGROUND

Most recombinant Komagataella phaffii (Pichia pastoris) strains for protein production are generated by genomic integration of expression cassettes. The clonal variability in gene copy numbers, integration loci and consequently product titers limit the aptitude for high throughput applications in drug discovery, enzyme engineering or most comparative analyses of genetic elements such as promoters or secretion signals. Circular episomal plasmids with an autonomously replicating sequence (ARS), an alternative which would alleviate some of these limitations, are inherently unstable in K. phaffii. Permanent selection pressure, mostly enabled by antibiotic resistance or auxotrophy markers, is crucial for plasmid maintenance and hardly scalable for production. The establishment and use of extrachromosomal ARS plasmids with key genes of the glycerol metabolism (glycerol kinase 1, GUT1, and triosephosphate isomerase 1, TPI1) as selection markers was investigated to obtain a system with high transformation rates that can be directly used for scalable production processes in lab scale bioreactors.

RESULTS

In micro-scale deep-well plate experiments, ARS plasmids employing the Ashbya gossypii TEF1 (transcription elongation factor 1) promoter to regulate transcription of the marker gene were found to deliver high transformation efficiencies and the best performances with the reporter protein (CalB, lipase B of Candida antarctica) for both, the GUT1- and TPI1-based, marker systems. The GUT1 marker-bearing strain surpassed the reference strain with integrated expression cassette by 46% upon re-evaluation in shake flask cultures regarding CalB production, while the TPI1 system was slightly less productive compared to the control. In 5 L bioreactor methanol-free fed-batch cultivations, the episomal production system employing the GUT1 marker led to 100% increased CalB activity in the culture supernatant compared to integration construct.

CONCLUSIONS

For the first time, a scalable and methanol-independent expression system for recombinant protein production for K. phaffii using episomal expression vectors was demonstrated. Expression of the GUT1 selection marker gene of the new ARS plasmids was refined by employing the TEF1 promoter of A. gossypii. Additionally, the antibiotic-free marker toolbox for K. phaffii was expanded by the TPI1 marker system, which proved to be similarly suited for the use in episomal plasmids as well as integrative expression constructs for the purpose of recombinant protein production.

Protein production dynamics and physiological adaptation of recombinant Komagataella phaffii at near-zero growth rates

BACKGROUND

Specific productivity (qP) in yeast correlates with growth, typically peaking at intermediate or maximum specific growth rates (μ). Understanding the factors limiting productivity at extremely low μ might reveal decoupling strategies, but knowledge of production dynamics and physiology in such conditions is scarce. Retentostats, a type of continuous cultivation, enable the well-controlled transition to near-zero µ through the combined retention of biomass and limited substrate supply. Recombinant Komagataella phaffii (syn Pichia pastoris) secreting a bivalent single domain antibody (VHH) was cultivated in aerobic, glucose-limited retentostats to investigate recombinant protein production dynamics and broaden our understanding of relevant physiological adaptations at near-zero growth conditions.

RESULTS

By the end of the retentostat cultivation, doubling times of approx. two months were reached, corresponding to µ = 0.00047 h-1. Despite these extremely slow growth rates, the proportion of viable cells remained high, and de novo synthesis and secretion of the VHH were observed. The average qP at the end of the retentostat was estimated at 0.019 mg g-1 h-1. Transcriptomics indicated that genes involved in protein biosynthesis were only moderately downregulated towards zero growth, while secretory pathway genes were mostly regulated in a manner seemingly detrimental to protein secretion. Adaptation to near-zero growth conditions of recombinant K. phaffii resulted in significant changes in the total protein, RNA, DNA and lipid content, and lipidomics revealed a complex adaptation pattern regarding the lipid class composition. The higher abundance of storage lipids as well as storage carbohydrates indicates that the cells are preparing for long-term survival.

CONCLUSIONS

In conclusion, retentostat cultivation proved to be a valuable tool to identify potential engineering targets to decouple growth and protein production and gain important insights into the physiological adaptation of K. phaffii to near-zero growth conditions.

Methanol bioconversion into C3, C4, and C5 platform chemicals by the yeast Ogataea polymorpha

BACKGROUND

One carbon (C1) molecules such as methanol have the potential to become sustainable feedstocks for biotechnological processes, as they can be derived from CO2 and green hydrogen, without the need for arable land. Therefore, we investigated the suitability of the methylotrophic yeast Ogataea polymorpha as a potential production organism for platform chemicals derived from methanol. We selected acetone, malate, and isoprene as industrially relevant products to demonstrate the production of compounds with 3, 4, or 5 carbon atoms, respectively.

RESULTS

We successfully engineered O. polymorpha for the production of all three molecules and demonstrated their production using methanol as carbon source. We showed that the metabolism of O. polymorpha is well suited to produce malate as a product and demonstrated that the introduction of an efficient malate transporter is essential for malate production from methanol. Through optimization of the cultivation conditions in shake flasks, which included pH regulation and constant substrate feeding, we were able to achieve a maximum titer of 13 g/L malate with a production rate of 3.3 g/L/d using methanol as carbon source. We further demonstrated the production of acetone and isoprene as additional heterologous products in O. polymorpha, with maximum titers of 13.6 mg/L and 4.4 mg/L, respectively.

CONCLUSION

These findings highlight how O. polymorpha has the potential to be applied as a versatile cell factory and contribute to the limited knowledge on how methylotrophic yeasts can be used for the production of low molecular weight biochemicals from methanol. Thus, this study can serve as a point of reference for future metabolic engineering in O. polymorpha and process optimization efforts to boost the production of platform chemicals from renewable C1 carbon sources.

Design and Construction of Artificial Biological Systems for One-Carbon Utilization

The third-generation (3G) biorefinery aims to use microbial cell factories or enzymatic systems to synthesize value-added chemicals from one-carbon (C1) sources, such as CO2, formate, and methanol, fueled by renewable energies like light and electricity. This promising technology represents an important step toward sustainable development, which can help address some of the most pressing environmental challenges faced by modern society. However, to establish processes competitive with the petroleum industry, it is crucial to determine the most viable pathways for C1 utilization and productivity and yield of the target products. In this review, we discuss the progresses that have been made in constructing artificial biological systems for 3G biorefineries in the last 10 years. Specifically, we highlight the representative works on the engineering of artificial autotrophic microorganisms, tandem enzymatic systems, and chemo-bio hybrid systems for C1 utilization. We also prospect the revolutionary impact of these developments on biotechnology. By harnessing the power of 3G biorefinery, scientists are establishing a new frontier that could potentially revolutionize our approach to industrial production and pave the way for a more sustainable future.

Methanol-based biomanufacturing of fuels and chemicals using native and synthetic methylotrophs

Methanol has recently gained significant attention as a potential carbon substrate for the production of fuels and chemicals, owing to its high degree of reduction, abundance, and low price. Native methylotrophic yeasts and bacteria have been investigated for the production of fuels and chemicals. Alternatively, synthetic methylotrophic strains are also being developed by reconstructing methanol utilization pathways in model microorganisms, such as Escherichia coli. Owing to the complex metabolic pathways, limited availability of genetic tools, and methanol/formaldehyde toxicity, the high-level production of target products for industrial applications are still under development to satisfy commercial feasibility. This article reviews the production of biofuels and chemicals by native and synthetic methylotrophic microorganisms. It also highlights the advantages and limitations of both types of methylotrophs and provides an overview of ways to improve their efficiency for the production of fuels and chemicals from methanol.

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

Enabling growth-decoupled Komagataella phaffii recombinant protein production based on the methanol-free PDH promoter

The current transition towards the circular bioeconomy requires a rational development of biorefineries to sustainably fulfill the present demands. The use of Komagataella phaffii (Pichia pastoris) can meet this challenge, since it has the capability to use crude glycerol as a carbon-source, a by-product from the biodiesel industry, while producing high- and low-added value products. Recombinant protein production (RPP) using K. phaffii has often been driven either by the methanol induced AOX1 promoter (P_AOX1) and/or the constitutive GAP promoter (P_GAP). In the last years, strong efforts have been focused on developing novel expression systems that expand the toolbox variety of K. phaffii to efficiently produce diverse proteins that requires different strategies. In this work, a study was conducted towards the development of methanol-free expression system based on a heat-shock gene promoter (P_DH) using glycerol as sole carbon source. Using this promoter, the recombinant expression is strongly induced in carbon-starving conditions. The classical P_GAP was used as a benchmark, taking for both strains the lipase B from Candida antarctica (CalB) as model protein. Titer of CalB expressed under P_DH outperformed P_GAP controlled expression in shake-flask cultivations when using a slow-release continuous feeding technology, confirming that P_DH is induced under pseudo-starving conditions. This increase was also confirmed in fed-batch cultivations. Several optimization rounds were carried out for P_DH under different feeding and osmolarity conditions. In all of them the P_DH controlled process outperformed the P_GAP one in regard to CalB titer. The best P_DH approach reached 3.6-fold more specific productivity than P_GAP fed-batch at low μ. Compared to the optimum approach for P_GAP-based process, the best P_DH fed-batch strategy resulted in 2.3-fold higher titer, while the specific productivity was very similar. To summarize, P_DH is an inducible promoter that exhibited a non-coupled growth regulation showing high performance, which provides a methanol-free additional solution to the usual growth-coupled systems for RPP. Thus, this novel system emerges as a potential alternative for K. phaffii RPP bioprocess and for revaluing crude glycerol, promoting the transition towards a circular economy.

Metabolic engineering of Kluyveromyces marxianus for biomass-based applications

Kluyveromyces marxianus ATCC 26,548 was cultivated in aerobic chemostats with [1-13C] and [U-13C] glucose as carbon source under three different growth conditions (0.10, 0.25, and 0.5 h-1) to evaluate metabolic fluxes. Carbon balances closed always within 97-102%. Growth was carbon limited, and the cell yield on glucose was the same. The extracellular side-product formation was very low, totaling 0.0008 C-mol C-mol-1 substrate at 0.5 h-1. The intracellular flux ratios did not show significant variation for metabolic flux analysis from labelling and biomass composition and metabolic flux ratio analysis from labelling. The observed strictly oxidative metabolism and the stability of the metabolism in terms of fluxes even at high growth rates, without triggering out the synthesis of by-products, is an extremely desired condition that underlines the potential of K. marxianus for biotechnological biomass-related applications and the comprehension of the metabolic pools and pathways is an important step to engineering this organism.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03324-x.

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.

Innovative Bioprocess Strategies Combining Physiological Control and Strain Engineering of Pichia pastoris to Improve Recombinant Protein Production

The combination of strain and bioprocess engineering strategies should be considered to obtain the highest levels of recombinant protein production (RPP) while assuring product quality and process reproducibility of heterologous products. In this work, two complementary approaches were investigated to improve bioprocess efficiency based on the yeast P. pastoris. Firstly, the performance of two Candida rugosa lipase 1 producer clones with different gene dosage under the regulation of the constitutive P_GAP were compared in chemostat cultures with different oxygen-limiting conditions. Secondly, hypoxic conditions in carbon-limited fed-batch cultures were applied by means of a physiological control based on the respiratory quotient (RQ). Stirring rate was selected to maintain RQ between 1.4 and 1.6, since it was found to be the most favorable in chemostat. As the major outcome, between 2-fold and 4-fold higher specific production rate (q_P) values were observed when comparing multicopy clone (MCC) and single-copy clone (SCC), both in chemostat and fed-batch. Additionally, when applying oxygen limitation, between 1.5-fold and 3-fold higher q_P values were obtained compared with normoxic conditions. Thus, notable increases of up to 9-fold in the production rates were reached. Furthermore, transcriptional analysis of certain key genes related to RPP and central carbon metabolism were performed. Results seem to indicate the presence of a limitation in post-transcriptional protein processing steps and a possible transcription attenuation of the target gene in the strains with high gene dosage. The entire approach, including both strain and bioprocess engineering, represents a relevant novelty involving physiological control in Pichia cell factory and is of crucial interest in bioprocess optimization, boosting RPP, allowing bioproducts to be economically competitive in the market, and helping develop the bioeconomy.

Aerobic Utilization of Methanol for Microbial Growth and Production

Methanol is a reduced one-carbon (C1) compound. It supports growth of aerobic methylotrophs that gain ATP from reduced redox equivalents by respiratory phosphorylation in their electron transport chains. Notably, linear oxidation of methanol to carbon dioxide may yield three reduced redox equivalents if methanol oxidation is NAD-dependent as, e.g., in Bacillus methanolicus. Methanol has a higher degree of reduction per carbon than glucose (6 vs. 4), and thus, lends itself as an ideal carbon source for microbial production of reduced target compounds. However, C-C bond formation in the RuMP or serine cycle, a prerequisite for production of larger molecules, requires ATP and/or reduced redox equivalents. Moreover, heat dissipation and a high demand for oxygen during catabolic oxidation of methanol may pose challenges for fermentation processes. In this chapter, we summarize metabolic pathways for aerobic methanol utilization, aerobic methylotrophs as industrial production hosts, strain engineering, and methanol bioreactor processes. In addition, we provide technological and market outlooks.

Strategies and challenges with the microbial conversion of methanol to high-value chemicals

As alternatives to traditional fermentation substrates, methanol (CH₃ OH), carbon dioxide (CO₂ ) and methane (CH₄ ) represent promising one-carbon (C1) sources that are readily available at low-cost and share similar metabolic pathway. Of these C1 compounds, methanol is used as a carbon and energy source by native methylotrophs, and can be obtained from CO₂ and CH₄ by chemical catalysis. Therefore, constructing and rewiring methanol utilization pathways may enable the use of one-carbon sources for microbial fermentations. Recent bioengineering efforts have shown that both native and nonnative methylotrophic organisms can be engineered to convert methanol, together with other carbon sources, into biofuels and other commodity chemicals. However, many challenges remain and must be overcome before industrial-scale bioprocessing can be established using these engineered cell refineries. Here, we provide a comprehensive summary and comparison of methanol metabolic pathways from different methylotrophs, followed by a review of recent progress in engineering methanol metabolic pathways in vitro and in vivo to produce chemicals. We discuss the major challenges associated with establishing efficient methanol metabolic pathways in microbial cells, and propose improved designs for future engineering.

Established tools and emerging trends for the production of recombinant proteins and metabolites in Pichia pastoris

Besides bakers' yeast, the methylotrophic yeast Komagataella phaffii (also known as Pichia pastoris) has been developed into the most popular yeast cell factory for the production of heterologous proteins. Strong promoters, stable genetic constructs and a growing collection of freely available strains, tools and protocols have boosted this development equally as thorough genetic and cell biological characterization. This review provides an overview of state-of-the-art tools and techniques for working with P. pastoris, as well as guidelines for the production of recombinant proteins with a focus on small-scale production for biochemical studies and protein characterization. The growing applications of P. pastoris for in vivo biotransformation and metabolic pathway engineering for the production of bulk and specialty chemicals are highlighted as well.

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.

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.

Evaluation of Product Distribution in Chemostat and Batch Fermentation in Lactic Acid-Producing Komagataella phaffii Strains Utilizing Glycerol as Substrate

Lactic acid is the monomeric unit of polylactide (PLA), a bioplastic widely used in the packaging, automotive, food, and pharmaceutical industries. Previously, the yeast Komagataella phaffii was genetically modified for the production of lactate from glycerol. For this, the bovine L-lactate dehydrogenase- (LDH)-encoding gene was inserted and the gene encoding the pyruvate decarboxylase (PDC) was disrupted, resulting in the GLp strain. This showed a yield of 67% L-lactic acid and 20% arabitol as a by-product in batches with oxygen limitation. Following up on these results, the present work endeavored to perform a detailed study of the metabolism of this yeast, as well as perturbing arabitol synthesis in an attempt to increase lactic acid titers. The GLp strain was cultivated in a glycerol-limited chemostat at different dilution rates, confirming that the production of both lactic acid and arabitol is dependent on the specific growth rate (and consequently on the concentration of the limiting carbon source) as well as on the oxygen level. Moreover, disruption of the gene encoding arabitol dehydrogenase (ArDH) was carried out, resulting in an increase of 20% in lactic acid and a 50% reduction in arabitol. This study clarifies the underlying metabolic reasons for arabitol formation in K. phaffii and points to ways for improving production of lactic acid using K. phaffii as a biocatalyst.

CO2 to succinic acid - Estimating the potential of biocatalytic routes

Microbial carbon dioxide assimilation and conversion to chemical platform molecules has the potential to be developed as economic, sustainable processes. The carbon dioxide assimilation can proceed by a variety of natural pathways and recently even synthetic CO2 fixation routes have been designed. Early assessment of the performance of the different carbon fixation alternatives within biotechnological processes is desirable to evaluate their potential. Here we applied stoichiometric metabolic modeling based on physiological and process data to evaluate different process variants for the conversion of C1 carbon compounds to the industrial relevant platform chemical succinic acid. We computationally analyzed the performance of cyanobacteria, acetogens, methylotrophs, and synthetic CO2 fixation pathways in Saccharomyces cerevisiae in terms of production rates, product yields, and the optimization potential. This analysis provided insight into the economic feasibility and allowed to estimate the future industrial applicability by estimating overall production costs. With reported, or estimated data of engineered or wild type strains, none of the simulated microbial succinate production processes showed a performance allowing competitive production. The main limiting factors were identified as gas and photon transfer and metabolic activities whereas metabolic network structure was not restricting. In simulations with optimized parameters most process alternatives reached economically interesting values, hence, represent promising alternatives to sugar-based fermentations.

Comparison of Yeasts as Hosts for Recombinant Protein Production

Recombinant protein production emerged in the early 1980s with the development of genetic engineering tools, which represented a compelling alternative to protein extraction from natural sources. Over the years, a high level of heterologous protein was made possible in a variety of hosts ranging from the bacteria Escherichia coli to mammalian cells. Recombinant protein importance is represented by its market size, which reached $1654 million in 2016 and is expected to reach $2850.5 million by 2022. Among the available hosts, yeasts have been used for producing a great variety of proteins applied to chemicals, fuels, food, and pharmaceuticals, being one of the most used hosts for recombinant production nowadays. Historically, Saccharomyces cerevisiae was the dominant yeast host for heterologous protein production. Lately, other yeasts such as Komagataella sp., Kluyveromyces lactis, and Yarrowia lipolytica have emerged as advantageous hosts. In this review, a comparative analysis is done listing the advantages and disadvantages of using each host regarding the availability of genetic tools, strategies for cultivation in bioreactors, and the main techniques utilized for protein purification. Finally, examples of each host will be discussed regarding the total amount of protein recovered and its bioactivity due to correct folding and glycosylation patterns.

Recent Advances in Pichia pastoris as Host for Heterologous Expression System for Lipases: A Review

The production of heterologous lipases is one of the most promising strategies to increase the productivity of the bioprocesses and to reduce costs, with the final objective that more industrial lipase applications could be implemented.In this chapter, an overview of the new success in synthetic biology, with traditional molecular genetic techniques and bioprocess engineering in the last 5 years in the cell factory Pichia pastoris, the most promising host system for heterologous lipase production, is presented.The goals get on heterologous Candida antarctica, Rhizopus oryzae, and Candida rugosa lipases, three of the most common lipases used in biocatalysis, are showed. Finally, new cell factories producing heterologous lipases are presented.

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.

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.

Metabolic modeling to identify engineering targets for Komagataella phaffii: The effect of biomass composition on gene target identification

Genome‐scale metabolic models are valuable tools for the design of novel strains of industrial microorganisms, such as Komagataella phaffii (syn. Pichia pastoris). However, as is the case for many industrial microbes, there is no executable metabolic model for K. phaffiii that confirms to current standards by providing the metabolite and reactions IDs, to facilitate model extension and reuse, and gene‐reaction associations to enable identification of targets for genetic manipulation. In order to remedy this deficiency, we decided to reconstruct the genome‐scale metabolic model of K. phaffii by reconciling the extant models and performing extensive manual curation in order to construct an executable model (Kp.1.0) that conforms to current standards. We then used this model to study the effect of biomass composition on the predictive success of the model. Twelve different biomass compositions obtained from published empirical data obtained under a range of growth conditions were employed in this investigation. We found that the success of Kp1.0 in predicting both gene essentiality and growth characteristics was relatively unaffected by biomass composition. However, we found that biomass composition had a profound effect on the distribution of the fluxes involved in lipid, DNA, and steroid biosynthetic processes, cellular alcohol metabolic process, and oxidation‐reduction process. Furthermore, we investigated the effect of biomass composition on the identification of suitable target genes for strain development. The analyses revealed that around 40% of the predictions of the effect of gene overexpression or deletion changed depending on the representation of biomass composition in the model. Considering the robustness of the in silico flux distributions to the changing biomass representations enables better interpretation of experimental results, reduces the risk of wrong target identification, and so both speeds and improves the process of directed strain development.

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

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.

PFA toolbox: a MATLAB tool for Metabolic Flux Analysis

BACKGROUND

Metabolic Flux Analysis (MFA) is a methodology that has been successfully applied to estimate metabolic fluxes in living cells. However, traditional frameworks based on this approach have some limitations, particularly when measurements are scarce and imprecise. This is very common in industrial environments. The PFA Toolbox can be used to face those scenarios.

RESULTS

Here we present the PFA (Possibilistic Flux Analysis) Toolbox for MATLAB, which simplifies the use of Interval and Possibilistic Metabolic Flux Analysis. The main features of the PFA Toolbox are the following: (a) It provides reliable MFA estimations in scenarios where only a few fluxes can be measured or those available are imprecise. (b) It provides tools to easily plot the results as interval estimates or flux distributions. (c) It is composed of simple functions that MATLAB users can apply in flexible ways. (d) It includes a Graphical User Interface (GUI), which provides a visual representation of the measurements and their uncertainty. (e) It can use stoichiometric models in COBRA format. In addition, the PFA Toolbox includes a User's Guide with a thorough description of its functions and several examples.

CONCLUSIONS

The PFA Toolbox for MATLAB is a freely available Toolbox that is able to perform Interval and Possibilistic MFA estimations.

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.

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.

Systems-level organization of yeast methylotrophic lifestyle

Background

Some yeasts have evolved a methylotrophic lifestyle enabling them to utilize the single carbon compound methanol as a carbon and energy source. Among them, Pichia pastoris (syn. Komagataella sp.) is frequently used for the production of heterologous proteins and also serves as a model organism for organelle research. Our current knowledge of methylotrophic lifestyle mainly derives from sophisticated biochemical studies which identified many key methanol utilization enzymes such as alcohol oxidase and dihydroxyacetone synthase and their localization to the peroxisomes. C1 assimilation is supposed to involve the pentose phosphate pathway, but details of these reactions are not known to date.

Results

In this work we analyzed the regulation patterns of 5,354 genes, 575 proteins, 141 metabolites, and fluxes through 39 reactions of P. pastoris comparing growth on glucose and on a methanol/glycerol mixed medium, respectively. Contrary to previous assumptions, we found that the entire methanol assimilation pathway is localized to peroxisomes rather than employing part of the cytosolic pentose phosphate pathway for xylulose-5-phosphate regeneration. For this purpose, P. pastoris (and presumably also other methylotrophic yeasts) have evolved a duplicated methanol inducible enzyme set targeted to peroxisomes. This compartmentalized cyclic C1 assimilation process termed xylose-monophosphate cycle resembles the principle of the Calvin cycle and uses sedoheptulose-1,7-bisphosphate as intermediate. The strong induction of alcohol oxidase, dihydroxyacetone synthase, formaldehyde and formate dehydrogenase, and catalase leads to high demand of their cofactors riboflavin, thiamine, nicotinamide, and heme, respectively, which is reflected in strong up-regulation of the respective synthesis pathways on methanol. Methanol-grown cells have a higher protein but lower free amino acid content, which can be attributed to the high drain towards methanol metabolic enzymes and their cofactors. In context with up-regulation of many amino acid biosynthesis genes or proteins, this visualizes an increased flux towards amino acid and protein synthesis which is reflected also in increased levels of transcripts and/or proteins related to ribosome biogenesis and translation.

Conclusions

Taken together, our work illustrates how concerted interpretation of multiple levels of systems biology data can contribute to elucidation of yet unknown cellular pathways and revolutionize our understanding of cellular biology.

Electronic supplementary material

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

Engineering the supply chain for protein production/secretion in yeasts and mammalian cells

Metabolic bottlenecks play an increasing role in yeasts and mammalian cells applied for high-performance production of proteins, particularly of pharmaceutical ones that require complex posttranslational modifications. We review the present status and developments focusing on the rational metabolic engineering of such cells to optimize the supply chain for building blocks and energy. Methods comprise selection of beneficial genetic modifications, rational design of media and feeding strategies. Design of better producer cells based on whole genome-wide metabolic network analysis becomes increasingly possible. High-resolution methods of metabolic flux analysis for the complex networks in these compartmented cells are increasingly available. We discuss phenomena that are common to both types of organisms but also those that are different with respect to the supply chain for the production and secretion of pharmaceutical proteins.

Validation of an FBA model for Pichia pastoris in chemostat cultures

Background

Constraint-based metabolic models and flux balance analysis (FBA) have been extensively used in the last years to investigate the behavior of cells and also as basis for different industrial applications. In this context, this work provides a validation of a small-sized FBA model of the yeast Pichia pastoris. Our main objective is testing how accurate is the hypothesis of maximum growth to predict the behavior of P. pastoris in a range of experimental environments.

Results

A constraint-based model of P. pastoris was previously validated using metabolic flux analysis (MFA). In this paper we have verified the model ability to predict the cells behavior in different conditions without introducing measurements, experimental parameters, or any additional constraint, just by assuming that cells will make the best use of the available resources to maximize its growth. In particular, we have tested FBA model ability to: (a) predict growth yields over single substrates (glucose, glycerol, and methanol); (b) predict growth rate, substrate uptakes, respiration rates, and by-product formation in scenarios where different substrates are available (glucose, glycerol, methanol, or mixes of methanol and glycerol); (c) predict the different behaviors of P. pastoris cultures in aerobic and hypoxic conditions for each single substrate. In every case, experimental data from literature are used as validation.

Conclusions

We conclude that our predictions based on growth maximisation are reasonably accurate, but still far from perfect. The deviations are significant in scenarios where P. pastoris grows on methanol, suggesting that the hypothesis of maximum growth could be not dominating in these situations. However, predictions are much better when glycerol or glucose are used as substrates. In these scenarios, even if our FBA model is small and imposes a strong assumption regarding how cells will regulate their metabolic fluxes, it provides reasonably good predictions in terms of growth, substrate preference, product formation, and respiration rates.

Electronic supplementary material

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

Impacts of high β-galactosidase expression on central metabolism of recombinant Pichia pastoris GS115 using glucose as sole carbon source via (13)C metabolic flux analysis

The yeast Pichia pastoris GS115 is a widely used microbial cell factory for the production of heterologous protein. In order to reveal the impacts of high heterologous protein expression on the central metabolism of Pichia pastoris GS115 using glucose as sole carbon source, we engineered a high β-galactosidase expression strain P. pastoris G1HL and a low expression control strain P. pastoris GHL through controlling the initiation strength of constitutive promoter pGAP. The carbon flux distributions in these two strains were quantified via (13)C metabolic flux analysis. Compared to the control strain, G1HL showed a lower growth rate, a higher flux through glycolysis pathway, a higher flux through pentose phosphate pathway, and a lower flux through by-products secretion pathway. The metabolic flux redistribution in G1HL was thought to compensate the increased redox cofactors and energy demands caused by the high protein expression. Although the fluxes through Krebs cycle in two engineered strains were almost the same, they were significantly lower than those in wild strain. The enhanced expression of β-galactosidase by glutamate supplementation demonstrated the potential of P. pastoris GS115 to catabolize more carbon through the Krebs cycle for even higher protein expression. In conclusion, our work indicates that P. pastoris GS115 can readjusts the central metabolism for higher heterologous protein expression and provides strategies for strain development or process optimization for enhancing production of heterologous protein.

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