- Open Access
Development of a fed-batch process for a recombinant Pichia pastoris Δoch1 strain expressing a plant peroxidase
Microbial Cell Factoriesvolume 14, Article number: 1 (2015)
Pichia pastoris is a prominent host for recombinant protein production, amongst other things due to its capability of glycosylation. However, N-linked glycans on recombinant proteins get hypermannosylated, causing problems in subsequent unit operations and medical applications. Hypermannosylation is triggered by an α-1,6-mannosyltransferase called OCH1. In a recent study, we knocked out OCH1 in a recombinant P. pastoris CBS7435 MutS strain (Δoch1) expressing the biopharmaceutically relevant enzyme horseradish peroxidase. We characterized the strain in the controlled environment of a bioreactor in dynamic batch cultivations and identified the strain to be physiologically impaired. We faced cell cluster formation, cell lysis and uncontrollable foam formation.
In the present study, we investigated the effects of the 3 process parameters temperature, pH and dissolved oxygen concentration on 1) cell physiology, 2) cell morphology, 3) cell lysis, 4) productivity and 5) product purity of the recombinant Δoch1 strain in a multivariate manner. Cultivation at 30°C resulted in low specific methanol uptake during adaptation and the risk of methanol accumulation during cultivation. Cell cluster formation was a function of the C-source rather than process parameters and went along with cell lysis. In terms of productivity and product purity a temperature of 20°C was highly beneficial. In summary, we determined cultivation conditions for a recombinant P. pastoris Δoch1 strain allowing high productivity and product purity.
The methylotrophic yeast Pichia pastoris is an attractive host for the recombinant production of proteins and biopharmaceuticals (e.g. [1-3]). It can grow on inexpensive media to high cell densities , numerous molecular manipulation tools are available  and high production titers are possible [5,6]. Due to the capacity of performing posttranslational modifications, like glycosylation, P. pastoris is attractive for the production of eukaryotic proteins (e.g. [3,7-10]). However, the glycosylation capacity of this yeast also is a curse: native glycosyltransferases recognize the aminoacid motif N-X-S/T and link N-glycans to the asparagine [11,12]. In contrast to mammalians, however, no trimming reactions of the attached glycans happen, but the glycans are further extended, a phenomenon known as hyperglycosylation . The first reaction of this cascade is catalyzed by an α-1,6-mannosyltransferase (OCH1) localized in the Golgi apparatus [14,15]. Hyperglycosylation describes a huge problem since not only the physico-chemical properties of the target protein get masked leading to difficulties in the downstream process , but also yeast derived glycans are not compatible with the human organism and can cause immunogenic reactions . Consequently, there have been numerous attempts to manipulate the native glycosylation machinery of P. pastoris (e.g. [18-22]). In a recent study, we deleted OCH1 in a recombinant P. pastoris strain (Δoch1) and physiologically characterized the strain in the controlled environment of a bioreactor . We purified the recombinant product horseradish peroxidase (EC 18.104.22.168; HRP; e.g. ) and analyzed catalytic constants, thermal stability as well as protein glycosylation. Although the Δoch1 strain produced the recombinant protein with shorter glycans of considerably increased homogeneity, the strain was physiologically impaired and thus hard to cultivate. We faced cell cluster formation, cell lysis and uncontrollable foam formation [25,26].
In the present study, we investigated the effects of the 3 process parameters temperature, pH and dissolved oxygen concentration (dO2) on 1) cell physiology, 2) cell morphology, 3) cell lysis, 4) productivity and 5) product purity in a multivariate manner to identify fed-batch operating conditions for the recombinant Δoch1 strain which give both high productivity and product purity, and hamper methanol accumulation as well as cell lysis and consequent foam formation.
Material and methods
A P. pastoris CBS7435 MutS Δoch1 strain carrying the gene coding for the HRP isoenzyme A2A was provided by Prof. Anton Glieder (University of Technology, Graz, Austria). Strain generation and isoenzyme characteristics were described previously [23,27]. A recombinant P. pastoris CBS7435 MutS strain with intact OCH1 expressing HRP A2A, hereafter called wildtype OCH1 strain, was included as reference.
Design of experiments
A 23-level full factorial screening approach with 2 centre points was set up with the program MODDE (Umetrics, Sweden) to explore the influence of the 3 factors temperature (20-30°C), pH (5.0-7.0) and dO2 (10–30%) as well as their linear interactions on different response parameters resulting in a total of 10 fed-batch cultivations (Table 1). We chose the limits for temperature with 20-30°C, since this temperature range is reported for yeasts (e.g. [28-31]). For pH we investigated values between pH 5.0 and 7.0 (e.g. ), since P. pastoris does not grow well at more acidic or alkaline conditions and also HRP exhibits high stability in this pH range . Finally, we investigated dO2 levels between 10–30%, which is again a range which had been used for P. pastoris before (e.g. [30,32,33]).
We analyzed the effects of the 3 factors on 1) strain physiology (specific substrate uptake rate during methanol adaptation (qs MeOH adapt), methanol accumulation during fed-batch (MeOHaccum), biomass yield (YX/S) and CO2 yield (YCO2/S)), 2) strain morphology (cell size distribution), 3) cell lysis (extracellular DNA content), 4) productivity (space-time-yield STY [U∙L−1∙h−1], specific productivity qp [U∙g−1∙h−1]) and 5) product purity (defined as the amount of target protein in relation to the amount of total extracellular protein [U∙mg−1]).
The P. pastoris Δoch1 strain expressing HRP isoenzyme A2A was cultivated in the controlled environment of a bioreactor. Batch and fed-batch phase were performed on glycerol, followed by a methanol adaptation pulse. Afterwards, a methanol fed-batch with a controlled feed rate corresponding to a certain specific substrate uptake rate of methanol (qs MeOH) was done.
Precultures were done in yeast nitrogen base medium (YNBM; 0.1 M potassium phosphate buffer pH 6.0, 3.4 g · L−1 YNB w/o amino acids and ammonia sulfate, 10 g · L−1 (NH4)2SO4, 400 mg · L−1 biotin, 20 g · L−1 glucose). Zeocine was added at a concentration of 100 μg∙L−1.
Batch and fed-batch cultivations were performed in 2-fold concentrated basal salt medium (BSM; 21.6 mL · L−1 85% phosphoric acid, 0.36 g∙L−1 CaSO4 · 2H2O, 27.24 g∙L−1 K2SO4, 4.48 g∙L−1 MgSO4 · 7H2O, 8.26 g∙L−1 KOH, 0.3 mL∙L−1 Antifoam Struktol J650, 4.35 mL∙L−1 PTM1, NH4OH as N-source). Trace element solution (PTM1) was made of 6.0 g · L−1 CuSO4 · 5H2O, 0.08 g · L−1 NaI, 3.0 g · L−1 MnSO4 · H2O, 0.2 g · L−1 Na2MoO4 · 2H2O, 0.02 g · L−1 H3BO3, 0.5 g · L−1 CoCl2, 20.0 g · L−1 ZnCl2, 65.0 g · L−1 FeSO4 · 7H2O, 0.2 g · L−1 biotin, 5 mL · L−1 H2SO4. Induction was carried out in presence of the heme-precursor Δ-aminolevulinic acid at a final concentration of 1 mM. The concentration of the base NH4OH was determined by titration with 0.25 M potassium hydrogen phthalate.
Frozen stocks (−80°C) were cultivated in 100 mL YNBM-Zeocine in 1,000 mL shake flasks at 30°C and 230 rpm for 48 h. Then, the preculture was transferred aseptically to the culture vessel. The inoculum volume was 10% of the final starting volume.
Batch and non-induced fed-batch
Batch cultivations were carried out in a 3 L working volume glass bioreactor (Infors, Switzerland). Basal salt medium was sterilized in the bioreactor and pH was adjusted to pH 6.0, which corresponds to the centre point of the subsequent DoE, by concentrated NH4OH solution after autoclaving. Sterile filtered trace elements were transferred to the reactor. Dissolved oxygen (dO2) was measured with a sterilizable fluorescence dissolved oxygen electrode (Visiferm DO425, Hamilton, Germany). The pH was measured with a sterilizable electrode (Easyferm™, Hamilton, Switzerland) and maintained constant with a PID controller using NH4OH solution (2 to 3 M). Base consumption was determined gravimetrically. Cultivation temperature was set to 30°C and agitation was fixed to 900 rpm. The culture was aerated with 1.0 vvm dried air and off-gas of the culture was measured by using an infrared cell for CO2 and a paramagnetic cell for O2 concentration (Servomax, Switzerland). Temperature, pH, dO2, agitation as well as CO2 and O2 in the off-gas were measured online and logged in a process information management system (PIMS; Lucullus, Biospectra, Switzerland). After the complete consumption of the substrate glycerol, indicated by an increase of dO2 and a drop in off-gas activity, an exponential fed-batch phase on glycerol with a specific growth rate of μ = 0.08 h−1 was performed. Based on the amount of glycerol used in the batch and the biomass yield on glycerol of YX/S = 0.47 g∙g−1, we calculated the biomass concentration after the batch phase. The feed rate was then determined by equations 1 and 2 and controlled by the PIMS. The fed-batch on glycerol was stopped when the reactor volume reached 2.2 L.
F0 = initial feed rate [g∙h−1]; X = calculated biomass concentration [g∙L−1]; V = volume in the bioreactor [L]; μ = specific growth rate [h−1]; δFeed = density glycerol feed [g∙L−1]; YX/S = biomass yield; cFeed = concentration feed [g∙L−1]; F = calculated feed rate [g∙h−1]; e, Euler constant; t = time [h]
Before fed-batch experiments, a single batch cultivation with dynamic methanol pulses was performed to determine qsmax MeOH at 20°C, 25°C and 30°C [25,26,34]. After complete consumption of glycerol, a methanol adaptation pulse (supplemented with 12 mL · L−1 PTM1) of a final concentration of 0.5% (v/v) was conducted. Following pulses were performed with 1% methanol (v/v). At least 3 consecutive methanol pulses were analyzed at each temperature. For each pulse, at least two samples were taken to determine the concentrations of substrate and product, as well as dry cell weight to calculate specific rates and yields. Furthermore, biomass was used to determine the correlation between dry cell weight (DCW) and optical density at 600 nm (OD600). The DCW of diluted samples (1:2, 1:4, 1:6, 1:8 and 1:10) was plotted against the corresponding OD600 values (Genesys 20; Thermo Scientific, Austria) resulting in a so-called α-factor, which was used for biomass calculation in subsequent fed-batch cultivations (equation 3).
X = calculated biomass [g∙L−1]; OD600 = optical density at 600 nm; α = correlation factor between DCW and OD600 value.
After glycerol depletion, temperature, pH and dO2 were changed according to Table 1. To control dO2, a PID controller was implemented regulating dO2 by aeration and not by stirrer speed. Thus, potential influences on cell cluster formation were omitted. After parameters reached the target values, a 0.5% (v/v) methanol adaptation pulse was applied. Concomitantly, the heme-precursor Δ-aminolevulinic acid was aseptically added at a final concentration of 1 mM. When the culture was adapted to methanol, indicated by an increase of dO2 and a drop in off-gas activity, a methanol fed-batch at a feed rate corresponding to a qs MeOH of 0.2 mmol∙g−1∙h−1 was performed for around 100 hours. The feed rate was determined by equation 4 and controlled by the PIMS.
F = calculated feed rate [g∙h−1]; X = biomass concentration calculated from OD600 [g∙L−1]; V = volume in the bioreactor [L]; qs MeOH = specific methanol uptake rate [mmol∙g−1∙h−1]; CFeed = density of methanol feed [g∙L−1]; cFeed = concentration of methanol feed [g∙L−1].
Dry cell weight was determined by centrifugation of 5 mL culture broth (5,000 rpm, 4°C, 10 min) in a laboratory centrifuge (Sigma 4K15, rotor 11156), washing the pellet with 5 mL deionized water and subsequent drying at 105°C to a constant weight in an oven. The enzymatic activity of HRP was measured using an ABTS assay in a CuBiAn XC enzymatic robot (Innovatis, Germany). Ten μl of sample were mixed with 140 μl 1 mM ABTS solution (50 mM KH2PO4, pH 6.5). The reaction mixture was incubated at 37°C for 5 min before the reaction was started by the addition of 20 μl 0.078% H2O2 (v/v). Changes in absorbance at 415 nm were measured for 80 seconds and rates were calculated. The standard curve was prepared using a commercially available HRP preparation (Type VI-A, Sigma-Aldrich, USA) in the range from 0.02 to 2.0 U · ml−1. Protein concentrations were determined at 595 nm using the Bradford Protein Assay Kit (Bio-Rad Laboratories GmbH, Austria) with bovine serum albumin as standard. Extracellular DNA content was measured by the Nanodrop 1000 device (ThermoScientific, Austria). Concentrations of methanol and potential metabolites were determined in cell-free samples by HPLC (Agilent Technologies, USA) equipped with an ion-exchange column (Supelcogel C-610H Sigma-Aldrich, USA) and a refractive index detector (Agilent Technologies, USA). The mobile phase was 0.1% H3PO4 with a constant flow rate of 0.5 mL · min−1 and the system was run isocratically at 30°C. All measurements were done in duplicates.
Changes in the morphology of the recombinant P. pastoris Δoch1 strain during the bioprocess were monitored by a Malvern Mastersizer 2000, measuring the cell size distribution in all samples. Frozen cell pellets (stored at −20°C) with a biomass concentration of 10–20 g∙L−1 were thawed and resuspended in 10 mL deionized water. After dispersing the samples with module Hydro2000S, 6 to 20 drops of cell suspension were dripped into the water tank of the Malvern Mastersizer 2000 until the laser diffraction reached a value between 10–12%. Samples were analyzed between 0 – 10,000 μm. Furthermore, all samples were analyzed microscopically using a Zeiss Epifluorescence Axio Observer Z1 deconvolution microscope (Carl Zeiss, Germany) equipped with a LD Plan-Neofluar 63x objective (+10x ocular) and the LED illumination system Colibri.
Electrophoresis was done with aliquots of supernatants obtained at the end of the cultivation. SDS-PAGE was performed using a 6% stacking gel and a 12% separating gel in 1x Tris-glycine buffer. Gels were run in the vertical electrophoresis Mini-PROTEAN Tetra Cell apparatus (Biorad; Austria) at 150 V for about 2 h. Gels were stained with Coomassie blue. The protein mass standard used was the PageRuler Prestained Ladder (Fermentas; Austria).
Protein identification and glycopeptide analysis by LC-ESI-MS
Relevant protein bands were cut out and digested in gel. S-alkylation with iodoacetamide and digestion with sequencing grade modified trypsin (Promega) were performed. The peptide mixture was analysed using a Dionex Ultimate 3000 system directly linked to a QTOF instrument (maXis 4G ETD, Bruker) equipped with the standard ESI source in the positive ion, DDA mode (= switching to MSMS mode for eluting peaks). MS-scans were recorded (range: 150–2200 Da) and the 6 highest peaks were selected for fragmentation. Instrument calibration was performed using ESIcalibration mixture (Agilent). For separation of the peptides a Thermo BioBasic C18 separation column (5 μm particle size, 150°0.360 mm) was used. A gradient from 95% solvent A and 5% solvent B (Solvent A: 65 mM ammonium formiate buffer, B: 100% ACCN) to 32% B in 45 min was applied, followed by a 15 min gradient from 32% B to 75% B, at a flow rate of 6 μL∙min−1. The analysis files were converted using Data Analysis 4.0 (Bruker) to XML files, which are suitable to perform MS/MS ion searches with MASCOT (embedded in ProteinScape 3.0, Bruker) for protein identification. Only proteins identified with at least 2 peptides with a protein score higher than 80 were accepted. For searches the SwissProt database was used.
Results and discussion
In the present study we developed a fed-batch bioprocess for a recombinant P. pastoris Δoch1 strain. We analyzed the effects of 3 process parameters (temperature, pH and dO2) on strain physiology and morphology , as well as on productivity and product purity in a multivariate manner.
Dynamic batch cultivation to determine qs max MeOH
We performed a dynamic batch cultivation with methanol pulses to determine the maximum specific methanol uptake rate (qs max MeOH) of the Δoch1 strain at 20°C, 25°C and 30°C [25,26]. The average qs max MeOH determined for at least 3 consecutive 1% (v/v) methanol pulses were 0.43 mmol∙g−1∙h−1 at 30°C, 0.74 mmol∙g−1∙h−1 at 25°C and 0.85 mmol∙g−1∙h−1 at 20°C. Interestingly, qs max MeOH was strongly temperature dependent. At 20°C, cells specifically consumed the double amount of methanol per time compared to 30°C. In previous studies with different microorganisms it was shown that substrate uptake usually declines with lower temperature [35-37]. Thus, we also analyzed a reference strain with intact OCH1, hereafter referred to as wildtype OCH1 strain, and indeed observed decreasing qs MeOH with decreasing temperature, namely 1.30 mmol∙g−1∙h−1 at 28°C, 1.20 mmol∙g−1∙h−1 at 24°C and 0.95 mmol∙g−1∙h−1 at 20°C. So far, we do not have a physiological explanation for the opposite behaviour of the recombinant Δoch1 strain. However, since the subsequent DoE covered a range from 20-30°C we designed the fed-batch strategy for the recombinant Δoch1 strain in a way to constantly feed at a rate corresponding to a rather low qs MeOH = 0.2 mmol∙g−1∙h−1. This was less than half of the lowest qs max MeOH determined in the dynamic batch experiment, and thus guaranteed a certain safety margin to avoid methanol accumulation.
DoE fed-batch cultivations
We analyzed the effects of the 3 process parameters temperature, pH and dO2 on 1) strain physiology, 2) strain morphology, 3) cell lysis, 4) productivity and 5) product purity of the recombinant Δoch1 strain. The goal of this multivariate study was to find operating conditions in fed-batch mode which give high productivity and product purity and concomitantly hamper methanol accumulation as well as cell lysis and consequent foam formation. Therefore, we conducted 10 fed-batch experiments in the controlled environment of a bioreactor (Table 1). The batch and the non-induced fed-batch were always performed under the same conditions. Only prior methanol adaptation (i.e. the phase during the Pichia pastoris cells get adapted to methanol) process parameters were changed. The carbon dioxide evolution rate (CER) depicting such a bioprocess is exemplarily shown for cultivation DoE6 in Figure 1.
Effects on strain physiology
By frequent sampling and subsequent analyses, we physiologically characterized the recombinant Δoch1 strain. The most relevant strain characteristic parameters are summarized in Table 2. The specific glycerol uptake rate in the batch phase (qs Gly batch) was very similar in all cultivations. Only during the induction phase strain characteristic parameters changed. The qs MeOH in the different fed-batches was close to the set value of 0.2 mmol∙g−1∙h−1, underlining the validity of the feeding strategy. Furthermore, closing C-balances underline the reliability of the data (Table 2).
Specific substrate uptake rate during methanol adaptation (qs MeOH adapt)
As shown in Figure 2, all 3 factors significantly affected qs MeOH adapt with temperature being the most significant one (p-values: temperature = 0.0011, pH = 0.0406, dO2 = 0.0055). The summary of fit plot indicated a valid model (R2 = 0.99, Q2 = 0.90, Model Validity = 0.95, Reproducibility = 0.95). At 30°C and a dO2 of 30% the recombinant Δoch1 strain only consumed 0.1 mmol∙g−1∙h−1, whereas this value was more than doubled at 20°C. Interestingly, the wildtype OCH1 strain showed the opposite as at 30°C 0.37 mmol∙g−1∙h−1 methanol were consumed but only half was consumed at 20°C (0.18 mmol∙g−1∙h−1). Summarizing, the lower the temperature, pH and dO2, the higher was qs MeOH adapt for the recombinant Δoch1 strain (Figure 2).
Methanol accumulation during fed-batch (MeOHaccum)
During induction the feed rate was regulated to correspond to a constant qs MeOH = 0.2 mmol∙g−1∙h−1, which was 2-fold lower than the lowest qs max MeOH at 30°C determined in the dynamic batch experiment. However, methanol still accumulated in cultivations DoE3, DoE4 and DoE10. For DoE3 methanol accumulation set on after around 80 h and for DoE10 after around 70 h of induction, whereas for DoE4 methanol accumulated right from the start. In Table 2 methanol concentrations at the end of cultivation are shown. Methanol only accumulated in cultivations at 30°C which is in agreement with both the results from the dynamic batch cultivation at 30°C where qs max for methanol was the lowest, as well as with the low qs MeOH adapt in fed batch cultivations at 30°C (Figure 2). However, we could not identify a significant factor causing this phenomenon since in DoE5 (30°C, pH 7.0, 10% dO2) no methanol accumulated (Table 2). Apparently, induction at 30°C favours methanol accumulation over time, a phenomenon which was observed before . However, using the current DoE screening approach based on only 3 factors and the underlying linear regression model, we could not identify significant factors and their interactions causing methanol accumulation. Apparently, a temperature of 30°C favours changes in cell physiology and a reduction of methanol uptake over time which is why a lower cultivation temperature of 20°C for fed-batch production processes with this Δoch1 strain is highly recommended. In contrast, the wildtype OCH1 strain can be cultivated at 30°C without experiencing changes in methanol uptake.
We investigated the influence of the 3 process parameters on both biomass yield (YX/S) and CO2 yield (YCO2/S). However, yields were not affected by varying the 3 process parameters in the design space (Table 2).
Effects on strain morphology
In a previous study, we had observed cluster formation of the recombinant Δoch1 strain during batch cultivation . To analyze this phenomenon in more detail and monitor it in the different phases of a fed-batch, samples taken during the different fed-batch cultivations were analyzed for cell size distribution in a Malvern Mastersizer. The cell size of an average P. pastoris cell lies between 4–6 μm . Of course, budding cells, where two or more daughter cells are still attached to the mother cell, are larger, which is why we categorized all signals < 15 μm as single and budding cells and all signals > 15 μm as cell clusters. We were especially interested in comparing the samples taken after batch and after non-induced fed-batch, which were both performed on glycerol, with the sample taken at the end of the induction phase on methanol (Table 3). We performed these measurements for all bioreactor cultivations, calculated mean values and standard deviations and compared the results with a wildtype OCH1 strain which was cultivated in fed-batch mode at 25°C, pH 6.0 and 20% dO2, corresponding to the DoE centre point (Table 1).
As shown in Table 3, the size distribution of the cells/cell clusters at the chosen time points of the different DoE cultivations was quite similar. In fact, all 3 factors were not significant for cluster formation. In Figure 3 the average amount of signals < 15 μm of the Δoch1 strain were compared to the signals of the wildtype OCH1 strain.
More than 80% of the signals for the wildtype OCH1 strain were smaller than 15 μm during the whole bioprocess (grey bars in Figure 3), whereas the Δoch1 strain showed different size distributions dependent on the cultivation phase. After the batch on glycerol more than 50% of the signals indicated cell clusters. These clusters were still prominent after the glycerol fed-batch (40%). Only after switching the substrate to methanol, the clusters slowly disappeared. After 100 h of induction, around 80% of the signals were < 15 μm, which was similar to the wildtype OCH1 strain. Cell cluster formation was also followed by light microscopy. In Figure 4, typical Δoch1 cell clusters  and budding cells for the control strain after the glycerol fed-batch are shown.
Apparently, cell cluster formation of the recombinant Δoch1 strain was a function of the C-source and not of any of the 3 investigated process parameters. When grown on methanol, cell clusters disappeared over time. However, Δoch1 cells stayed much smaller than wildtype OCH1 cells (Figure 4) which might result from the altered cell wall composition with enhanced chitin deposition .
Effects on cell lysis
Since we had observed intensive foam formation for the recombinant Δoch1 strain before , which describes a huge problem for potential scale up, we monitored cell lysis by measuring extracellular DNA content over time. In Table 4 the extracellular DNA contents at the end of the different cultivations are shown. Although the extracellular DNA content constantly increased over time, more than 80% of the final DNA amount was already present after the fed-batch on glycerol for all cultivations. This strongly indicated that cell lysis went along with cell cluster formation and was not affected by any of the 3 process parameters. We speculate that cells in the center of the clusters get limited in either nutrients or oxygen and thus lyse. Once these clusters disappear due to the switch from glycerol to methanol, cell lysis is diminished. Since none of the 3 process parameters affected lysis, we had to adapt the cultivation strategy during glycerol fed-batch elsewise to avoid extensive foam formation. Initially, agitation was fixed to 900 rpm and the culture was aerated with 1.0 vvm dried air to guarantee dO2 values of > 30%. However, to avoid intensive foam formation, we decreased agitation to 600 rpm and aeration to 0.5 vvm, but added pure oxygen. Thus, although cells still lysed during glycerol fed-batch, intensive foam formation was avoided.
Effects on productivity and product purity
The effects of temperature, pH and dO2 on qp (Additional file 1: Figure S1A), STY (Additional file 1: Figure S1B) and product purity (Additional file 1: Figure S1C) were analyzed and the respective data at the end of cultivation are summarized in Table 5.
The only significant factor for all 3 responses was the temperature, whereas pH and dO2 had no effect (p-values for qp: temperature = 0.000072, pH = 0.77, dO2 = 0.33; p-values for STY: temperature = 0.000026, pH = 0.37, dO2 = 0.90; p-values for specific activity: temperature = 0.000024, pH = 0.49, dO2 = 0.45). The summary of fit plots indicated valid models (for qp: R2 = 0.94, Q2 = 0.88, Model Validity = 0.77, Reproducibility = 0.97; for STY: R2 = 0.95, Q2 = 0.81, Model Validity = 0.55, Reproducibility = 0.99; for specific activity: R2 = 0.95, Q2 = 0.93, Model Validity = 0.87, Reproducibility = 0.95). To visualize the effects, we plotted the responses versus temperature and dO2 at pH 5.0 (Additional file 1: Figure S1). Summarizing, both highest productivity and product purity were obtained at 20°C and a low dO2 at pH 5.0. Apparently, at lower temperature the recombinant Δoch1 strain secreted more active HRP and less contaminating proteins compared to higher temperatures. This is not only important to increase the total amount of active product over time, but also for the subsequent downstream process. We observed a similar trend for the wildtype OCH1 strain where qp and product purity were both higher at lower temperature: (20°C: qp = 15 U∙g−1∙h−1, specific activity = 95 U∙mg−1; 24°C: qp = 11 U∙g−1∙h−1, specific activity = 50 U∙mg−1; 28°C: qp = 6 U∙g−1∙h−1, specific activity = 35 U∙mg−1).
Electrophoresis and glycopeptides analysis by LC-ESI-MS
Cell-free cultivation broths of the single fed-batch cultivations were analyzed by SDS-PAGE. Cell-free cultivation broth of a wildtype OCH1 strain as well as commercially available enzyme preparation from plant were included for comparison (Figure 5).
To identify the different HRP proteins, respective bands were excised (indicated in Figure 5) and glycopeptides were analyzed by HPLC-ESI-MS. In fact, we were able to identify the different HRP preparations (Additional file 2: Table S1). We observed apparent size differences between the 2 HRP A2A preparations from the Δoch1 strain and the wildtype OCH1 strain, which might result from the different degree of glycosylation (lanes 1 and 2), and between HRP A2A and HRP C1A from the wildtype OCH1 strains (lanes 2 and 5; Additional file 2: Table S1).
Verification runs at 15°C
The results of the DoE revealed a low cultivation temperature to be highly beneficial for productivity and product purity. To check if both responses could be further increased, 2 additional fed-batch cultivations at 15°C, pH 5.0 and dO2 of 10% were done. The average qs MeOH adapt was 0.18 mmol∙g−1∙h−1, which was comparable to the results at 20°C and 25°C (Table 2). Apparently, significant changes in qs MeOH adapt of the recombinant Δoch1 strain only happened between 25°C and 30°C. Surprisingly, at 15°C productivity and product purity were only half compared to 20°C. We determined a qp of 1.50 [U∙g−1∙h−1], a STY of 66 [U∙L−1∙h−1] and a specific activity of 23.3 [U∙mg−1] at 15°C, whereas in DoE2 (20°C, pH 5.0, dO2 10%), qp of 3.46 [U∙g−1∙h−1], a STY of 153 [U∙L−1∙h−1] and a specific activity of 56.5 [U∙mg−1] were measured. Thus, we concluded that the highest values for productivity and product purity, which are the main target variables in production processes, might be obtained between 15°C and 20°C, a pH of 5.0 and a dO2 of 10%.
In the present study we developed a fed-batch bioprocess for a recombinant P. pastoris Δoch1 strain. High productivity and product purity were reached when this strain was cultivated at pH 5.0, dO2 of 10% and a temperature of 20°C. At 30°C methanol accumulated over time due to apparent changes in cell metabolism. Cell cluster formation, which was accompanied by cell lysis, was dependent on the C-source. To avoid intensive foam formation during glycerol batch and fed-batch, aeration and stirrer speed had to be reduced.
Currently, we are investigating cell cluster formation on other C-sources, like glucose and sorbitol, and analyze the temperature range between 15-20°C in more detail to find the true optimum for the recombinant Δoch1 strain. However, the present study already describes a good basis for bioprocess engineers working with glyco-engineered Δoch1 yeast strains.
Zhou X, Yu Y, Tao J, Yu L. Production of LYZL6, a novel human c-type lysozyme, in recombinant Pichia pastoris employing high cell density fed-batch fermentation. J Biosci Bioeng. 2014;118:420–5.
Ahmad M, Hirz M, Pichler H, Schwab H. Protein expression in Pichia pastoris: recent achievements and perspectives for heterologous protein production. Appl Microbiol Biotechnol. 2014;98:5301–17.
Weinacker D, Rabert C, Zepeda AB, Figueroa CA, Pessoa A, Farias JG. Applications of recombinant Pichia pastoris in the healthcare industry. Braz J Microbiol. 2013;44:1043–8.
Felber M, Pichler H, Ruth C. Strains and molecular tools for recombinant protein production in Pichia pastoris. Methods Mol Biol. 2014;1152:87–111.
Hasslacher M, Schall M, Hayn M, Bona R, Rumbold K, Luckl J, et al. High-level intracellular expression of hydroxynitrile lyase from the tropical rubber tree Hevea brasiliensis in microbial hosts. Protein Expr Purif. 1997;11:61–71.
Werten MW, van den Bosch TJ, Wind RD, Mooibroek H, de Wolf FA. High-yield secretion of recombinant gelatins by Pichia pastoris. Yeast. 1999;15:1087–96.
Cereghino GPL, Cereghino JL, Ilgen C, Cregg JM. Production of recombinant proteins in fermenter cultures of the yeast Pichia pastoris. Curr Opin Biotechnol. 2002;13:329–32.
Li H, Hui X, Yang S, Hu X, Tang X, Li P, et al. High level expression, efficient purification and bioactivity assay of recombinant human platelet-derived growth factor AA dimer (PDGF-AA) from methylotrophic yeast Pichia pastoris. Protein Expr Purif. 2013;91:221–7.
Spadiut O, Capone S, Krainer F, Glieder A, Herwig C. Microbials for the production of monoclonal antibodies and antibody fragments. Trends Biotechnol. 2014;32:54–60.
Hamilton SR, Bobrowicz P, Bobrowicz B, Davidson RC, Li HJ, Mitchell T, et al. Production of complex human glycoproteins in yeast. Science. 2003;301:1244–6.
Gong B, Cukan M, Fisher R, Li H, Stadheim TA, Gerngross T. Characterization of N-linked glycosylation on recombinant glycoproteins produced in Pichia pastoris using ESI-MS and MALDI-TOF. Methods Mol Biol. 2009;534:213–23.
Li H, Miele RG, Mitchell TI, Gerngross TU. N-linked glycan characterization of heterologous proteins. Methods Mol Biol. 2007;389:139–50.
Montesino R, Garcia R, Quintero O, Cremata JA. Variation in N-linked oligosaccharide structures on heterologous proteins secreted by the methylotrophic yeast Pichia pastoris. Protein Expr Purif. 1998;14:197–207.
Nagasu T, Shimma Y, Nakanishi Y, Kuromitsu J, Iwama K, Nakayama K, et al. Isolation of new temperature-sensitive mutants of Saccharomyces cerevisiae deficient in mannose outer chain elongation. Yeast. 1992;8:535–47.
Nakayama K, Nagasu T, Shimma Y, Kuromitsu J, Jigami Y. OCH1 encodes a novel membrane bound mannosyltransferase: outer chain elongation of asparagine-linked oligosaccharides. EMBO J. 1992;11:2511–9.
Spadiut O, Rossetti L, Dietzsch C, Herwig C. Purification of a recombinant plant peroxidase produced in Pichia pastoris by a simple 2-step strategy. Protein Expr Purif. 2012;86:89–97.
Yang YL, Chang SH, Gong X, Wu J, Liu B. Expression, purification and characterization of low-glycosylation influenza neuraminidase in alpha-1,6-mannosyltransferase defective Pichia pastoris. Mol Biol Rep. 2012;39:857–64.
Wildt S, Gerngross TU. The humanization of N-glycosylation pathways in yeast. Nat Rev Microbiol. 2005;3:119–28.
Hamilton SR, Davidson RC, Sethuraman N, Nett JH, Jiang Y, Rios S, et al. Humanization of yeast to produce complex terminally sialylated glycoproteins. Science. 2006;313:1441–3.
Hamilton SR, Gerngross TU. Glycosylation engineering in yeast: the advent of fully humanized yeast. Curr Opin Biotechnol. 2007;18:387–92.
De Pourcq K, De Schutter K, Callewaert N. Engineering of glycosylation in yeast and other fungi: current state and perspectives. Appl Microbiol Biotechnol. 2010;87:1617–31.
Vervecken W, Kaigorodov V, Callewaert N, Geysens S, De Vusser K, Contreras R. In vivo synthesis of mammalian-like, hybrid-type N-glycans in Pichia pastoris. Appl Environ Microbiol. 2004;70:2639–46.
Krainer FW, Gmeiner C, Neutsch L, Windwarder M, Pletzenauer R, Herwig C, et al. Knockout of an endogenous mannosyltransferase increases the homogeneity of glycoproteins produced in Pichia pastoris. Sci Rep. 2013;3:3279.
Spadiut O, Herwig C. Production and purification of the multifunctional enzyme horseradish peroxidase: a review. Pharm Bioproc. 2013;1:283–95.
Dietzsch C, Spadiut O, Herwig C. A fast approach to determine a fed batch feeding profile for recombinant Pichia pastoris strains. Microb Cell Fact. 2011;10:85.
Dietzsch C, Spadiut O, Herwig C. A dynamic method based on the specific substrate uptake rate to set up a feeding strategy for Pichia pastoris. Microb Cell Fact. 2011;10:14.
Krainer FW, Pletzenauer R, Rossetti L, Herwig C, Glieder A, Spadiut O. Purification and basic biochemical characterization of 19 recombinant plant peroxidase isoenzymes produced in Pichia pastoris. Protein Expr Purif. 2013;95C:104–12.
Jeong GM, Lee YJ, Kim YS, Jeong KJ. High-level production of Fc-fused kringle domain in Pichia pastoris. J Ind Microbiol Biotechnol. 2014;41:989–96.
Peng Z, Wang A, Feng Q, Wang Z, Ivanova IV, He X, et al. High-level expression, purification and characterisation of porcine beta-defensin 2 in Pichia pastoris and its potential as a cost-efficient growth promoter in porcine feed. Appl Microbiol Biotechnol. 2014;98:5487–97.
Berdichevsky M, d'Anjou M, Mallem MR, Shaikh SS, Potgieter TI. Improved production of monoclonal antibodies through oxygen-limited cultivation of glycoengineered yeast. J Biotechnol. 2011;155:217–24.
Dragosits M, Frascotti G, Bernard-Granger L, Vazquez F, Giuliani M, Baumann K, et al. Influence of growth temperature on the production of antibody Fab fragments in different microbes: a host comparative analysis. Biotechnol Prog. 2011;27:38–46.
Potgieter TI, Cukan M, Drummond JE, Houston-Cummings NR, Jiang Y, Li F, et al. Production of monoclonal antibodies by glycoengineered Pichia pastoris. J Biotechnol. 2009;139:318–25.
Minning S, Serrano A, Ferrer P, Sola C, Schmid RD, Valero F. Optimization of the high-level production of Rhizopus oryzae lipase in Pichia pastoris. J Biotechnol. 2001;86:59–70.
Krainer FW, Dietzsch C, Hajek T, Herwig C, Spadiut O, Glieder A. Recombinant protein expression in Pichia pastoris strains with an engineered methanol utilization pathway. Microb Cell Fact. 2012;11:22.
Reay DS, Nedwell DB, Priddle J, Ellis-Evans JC. Temperature dependence of inorganic nitrogen uptake: reduced affinity for nitrate at suboptimal temperatures in both algae and bacteria. Appl Environ Microbiol. 1999;65:2577–84.
Pomeroy LR, Wiebe WJ. Temperature and substrates as interactive limiting factors for marine heterotrophic bacteria. Aquat Microb Ecol. 2001;23:187–204.
Nedwell DB. Effect of low temperature on microbial growth: lowered affinity for substrates limits growth at low temperature. Fems Microbiol Ecol. 1999;30:101–11.
Rambourg A, Clermont Y, Ovtracht L, Kepes F. Three-dimensional structure of tubular networks, presumably Golgi in nature, in various yeast strains: a comparative study. Anat Rec. 1995;243:283–93.
The authors thank the Austrian Science Fund FWF (project: P24861-B19) for financial support and Prof. Anton Glieder (University of Technology, Graz, Austria) for providing the recombinant P. pastoris Δoch1 strain. Furthermore we want to thank Simona Capone for great assistance in the lab.
The authors declare that they have no competing interests.
CG, AS and MF conducted bioreactor runs and biochemical experiments. SK helped with Malvern Mastersizer analyses. DM and FA conducted mass spectrometry analyses. CH helped in planning the study. CG, FA and OS analyzed the data. OS planned the study, supervised research and wrote the manuscript. All authors read and approved the final manuscript.
Contour plot showing A, specific productivity; B, space-time-yield, and C, specific activity as an indicator for product purity at the end of cultivation in dependence on temperature and dO2.
Results of the protein search with MASCOT (embedded in ProteinScape 3.0, Bruker) for protein identification of the excised bands. Only proteins identified with at least 2 peptides and a protein score higher than 80 were accepted. For searches the SwissProt database was used.