Skip to main content

A fast approach to determine a fed batch feeding profile for recombinant Pichia pastoris strains



The microorganism Pichia pastoris is a commonly used microbial host for the expression of recombinant proteins in biotechnology and biopharmaceutical industry. To speed up process development, a fast methodology to determine strain characteristic parameters, which are needed to subsequently set up fed batch feeding profiles, is required.


Here, we show the general applicability of a novel approach to quantify a certain minimal set of bioprocess-relevant parameters, i.e. the adaptation time of the culture to methanol, the specific substrate uptake rate during the adaptation phase and the maximum specific substrate uptake rate, based on fast and easy-to-do batch cultivations with repeated methanol pulses in a batch culture. A detailed analysis of the adaptation of different P. pastoris strains to methanol was conducted and revealed that each strain showed very different characteristics during adaptation, illustrating the need of individual screenings for an optimal parameter definition during this phase. Based on the results obtained in batch cultivations, dynamic feeding profiles based on the specific substrate uptake rate were employed for different P. pastoris strains. In these experiments the maximum specific substrate uptake rate, which had been defined in batch experiments, also represented the upper limit of methanol uptake, underlining the validity of the determined process-relevant parameters and the overall experimental strategy.


In this study, we show that a fast approach to determine a minimal set of strain characteristic parameters based on easy-to-do batch cultivations with methanol pulses is generally applicable for different P. pastoris strains and that dynamic fed batch strategies can be designed on the specific substrate uptake rate without running the risk of methanol accumulation.


Advances in molecular biology, cloning techniques and strain improvement allowed an increasing use of recombinant organisms for the industrial production of a variety of substances like organic acids, antibiotics, enzymes and amino acids. In this context the methylotrophic yeast Pichia pastoris is one of the most important host organisms for the expression of recombinant proteins. To meet industrial demands, a fast and easy-to-do characterization of recombinant Pichia strains to extract bioprocess-relevant strain characteristic parameters for the subsequent set-up of production processes is essential to speed up process development. Normally, this strain characterization procedure is done by time-consuming experiments, which require complex and costly equipment, like continuous cultures [1, 2] or several, consecutive fed batch cultivations, operated at different conditions [3, 4]. Parameters, which have to be extracted out of these experiments, describe the best operating conditions for each strain as well as the optimal condition for the adaptation of the culture from the growth substrate (e.g. glucose or glycerol) to the inducer methanol. To date, different strategies are employed regarding the adaptation of Pichia to methanol; two prominent examples are: 1) after a fed batch on glycerol or glucose a certain low concentration or flow of methanol is applied to the culture which is then increased to a pre-defined maximum and constantly maintained throughout the whole cultivation time [5, 6], and 2) the glycerol flow in the reactor is decreased following a linear function during a concomitant addition of methanol, a period which is called transition phase, to slowly adapt the culture to methanol [710]. These methods are often based on specific experiences with a certain strain, but are nevertheless often used as a general approach for different Pichia strains in following studies, without taking into account the specific requirements of the single strains during adaptation.

After adaptation of the culture to methanol, different feeding strategies can be employed for recombinant protein production with P. pastoris. Besides the two common strategies of either using a feed forward regime or a controlled specific growth rate (μ) [6, 7, 914], a few studies have also described the importance of the specific substrate uptake rate (qs) on recombinant protein production [15, 16]. A direct correlation between qs and the specific productivity (qp) was shown [15, 16], and it was clearly stated, that qs was the most important induction parameter in these experiments [17]. Based on those findings and motivated by problems which occur, when more traditional feeding strategies are applied (e.g. possible accumulation of methanol caused by changing cell capacities during cultivation or the need of expensive monitoring equipment to allow μ-controlled feeding [9, 14]), we have focused our research on the specific substrate uptake rate (qs) and have recently shown optimization potential using dynamic feeding profiles based on this parameter [18].

In our previous study we also developed a fast approach based on batch experiments with methanol pulses to extract a minimal set of strain characteristic parameters (i.e. Δtimeadapt - time for adaptation, qs adapt - specific substrate uptake rate during adaptation, qs max - maximum specific substrate uptake rate), which are required to set up a subsequent feeding regime based on qs. However, our previous study dealt with the development and the application of this approach for only one recombinant P. pastoris MutS strain [18]. In the present work, we characterized various P. pastoris strains with different phenotypes (Muts and Mut+) expressing different target enzymes using the above mentioned strategy. We analyzed the required time for adaptation to methanol of each strain in detail and could reliably derive certain strain characteristic parameters from pulse experiments to fed batch cultivations. With the variety of used strains in this study, we demonstrate that this approach is generally applicable for different P. pastoris strains and is thus a valuable tool for fast process development, which is especially interesting in an industrial environment.

Materials and methods

The experiments conducted in the present study were performed according to our previous study [18], and are thus only described briefly here.

Microorganisms and recombinant proteins

Different P. pastoris strains with different phenotypes expressing different target enzymes were used in this study to prove the general applicability of our strategy. A list of the various strains is given in Table 1.

Table 1 Different P. pastoris strains used in this study

All recombinant genes used in this study were under the control of the AOX1 promoter. The expressed HRP gene coded for the isoenzyme HRP C1A. The strain KM71H PDI HRP concomitantly expressed HRP and the chaperone protein disulfide isomerase (PDI), which was under the control of a modified AOX1 promoter [19]. The strains KM71H, KM71H HRP, KM71H PDI HRP and CBS7435 HRP were gratefully provided by Prof. Anton Glieder (Graz University of Technology, Austria). The strain SMD1168H GalOX was constructed by Spadiut et al., as described elsewhere [20].

Culture Media

Precultures were performed in complex yeast nitrogen base media (YNBM), whereas batch and fed batch cultivations were done in defined basal salt media (BSM; [21]). The glucose feed was prepared with glucose (250 g·l-1), trace element solution PTM1 (12 ml·l-1) and antifoam Struktol J650 (0.3 ml/l). The methanol feed was composed of methanol (300 g·l-1), PTM1 (4 ml·l-1) and Struktol J650 (0.3 ml·l-1). The induction period for HRP expression was carried out in the presence of δ-Aminolevulinic acid (δ-ALA) in a final concentration of 1 mM. The concentration of the base NH4OH was determined by titration with 0.25 M potassium hydrogen phthalate (KHP).

Experimental Procedure


Frozen stocks (-80°C) were precultivated in 100 ml of YNBM in 1000 ml shake flasks at 28°C and 230 rpm for max. 24 hours.

Batch cultivation with methanol pulses

Batch cultivations were carried out in a 5 l working volume glass bioreactor (Infors, Switzerland) at 28°C and a fixed agitation speed of 1200 rpm. The culture was aerated with 1 vvm dried air and off-gas was measured by using an infrared cell for CO2 and a paramagnetic cell for O2 concentration (Servomex, Switzerland). Process parameters were recorded and logged in a process information management system (PIMS; Lucullus, Biospectra, Switzerland). After the complete consumption of glucose, which was indicated by an increase of dissolved oxygen and a drop in off-gas activity, the first methanol pulse (adaptation pulse) with a final concentration of 0.5% (v/v) was conducted with pure methanol (supplemented with PTM1, 12 ml·l-1 of methanol). Following pulses were performed with 1% (v/v) concentration of methanol. For all strains, several pulses were conducted after the adaptation pulse to generate consistent data for each strain. For each methanol-pulse, at least two samples were taken to determine the concentrations of substrate and product as well as dry cell weight and OD600 to calculate the specific substrate uptake rate qs.

Fed batch cultivations

Fed batch cultivations were carried out in a 5 l working volume glass bioreactor (Infors, Switzerland) in 2-fold concentrated BSM medium at 28°C and 1500 rpm. The culture was aerated with at least 1 vvm to keep dissolved oxygen levels > 30%. In case, air flow was limited, pure oxygen was added. The fed batch feed was measured and controlled using a gravimetrically based PID flow controller. At several time points during fed batch cultivations, samples were taken and analyzed for accumulated methanol, biomass concentration (dry cell weight and optical density OD600) and, if applicable, enzymatic activity. Based on the total biomass content, feeding rates were adjusted manually corresponding to the defined qs set point. All fed batches described in this study were conducted in the same way: after an adaptation period at qs of 0.5 mmol·g-1·h-1, a stepwise increase of qs up to qs max of the respective strain was carried out with step times of 24 hours.

Analysis of growth- and expression-parameters

Dry cell weight (DCW), OD600, substrate concentrations as well as the catalytic activity of HRP were determined as described before [18]. However, in this study the ABTS solution for HRP activity measurements was prepared in 50 mM KH2PO4-buffer at pH 6.5 and the calibration range was expanded to 2.0 U·ml-1. Also GalOX activity was measured with an ABTS assay; i.e. a sample of diluted enzyme (10 μl) was added to 990 μl of assay buffer containing horseradish peroxidase (222 U) (Type VI-A, Sigma-Aldrich, P6782), ABTS (17.7 mg), KH2PO4-buffer (50 mM, pH 6.5) and D-galactose (300 mM). The absorbance change at 420 nm (ε420 = 42.3 mM-1·cm-1) was recorded at 30°C for 180 seconds. One Unit of GalOX activity was defined as the amount of enzyme necessary for the oxidation of 2 μmol of ABTS per min, corresponding to the consumption of 1 μmol of O2 per min. An additional post-translational activation of GalOX by adding CuSO4 to the samples in the presence of oxygen before activity measurements, as described elsewhere [20], was not executed.

Specific rate calculations

Batch cultivations

To obtain the specific rates for substrate uptake and productivity, samples were taken at certain time points during the methanol pulses (i.e. beginning of pulse, maximum off-gas of pulse, end of pulse) and analyzed offline for biomass content, methanol concentration and, if applicable, enzymatic activity. Determined values at the beginning and the end of the respective pulse were used to calculate an average rate of the specific substrate uptake, which was corrected for stripping using Antoine's equation, and the specific productivity. Errors for specific rates were set to 10%, according to our previous study [18]. Online calculated carbon dioxide evolution rate (CER) was divided by actual biomass concentrations to obtain the specific carbon dioxide production rate (qCO2). In addition, a time derivative of the qCO2 signal (i.e. qCO2') was calculated using a time window of 30 minutes (15 minutes before and 15 minutes after the actual time point).

Fed batch cultures

During different cultivation periods, representing defined qs set points, several samples were taken and OD600 measurements were used to calculate the actual total biomass content, which allowed adjustments of the methanol feed flow to the actual qs set point. Specific rates were calculated using DCW and the amount of consumed methanol, which was determined gravimetrically. Presented results correspond to an average value over the respective qs set point period. Again, errors for specific rates were set to 10%.

Results and Discussion

A fast approach to derive a minimal set of strain characteristic parameters relevant for bioprocess development

Each P. pastoris strain was cultivated in an easy-to-do batch system with methanol pulses to obtain certain strain characteristic parameters during the adaptation period (Δtimeadapt - time for adaptation of the culture to methanol, qs adapt - specific substrate uptake rate during the adaptation pulse) and the maximum specific substrate uptake rate (qs max). These parameters were extracted and consecutively transformed into a feeding profile for fed batch operations based on qs.

Adaptation of the culture to methanol

After depletion of glucose in batch cultivation, a methanol adaptation pulse with a final concentration of 0.5% (v/v) was applied. The time required to develop a maximum in off-gas activity was used to define Δtimeadapt, according to our previous study [18], and is shown here as the specific carbon dioxide production rate (qCO2; Figure 1).

Figure 1
figure 1

Adaptation pulses with 0.5% (v/v) methanol after glucose depletion for different P. pastoris strains. Straight line, specific carbon dioxide production rate qCO2; dashed line, time-derivative of the specific carbon dioxide production rate qCO2' (d(qCO2)/dt); A, KM71H; B, KM71H HRPa; C, KM71H PDI HRP; D, CBS7435 HRP; E, SMD1168H GalOX. adata taken from [18]

The used strains showed very different metabolic characteristics during the adaptation to methanol. The shortest Δtimeadapt of 3.5 h was detected for the KM71H strain, which was a MutS strain, not carrying a recombinant gene for heterologous protein expression. The other MutS strains KM71H HRP, KM71H PDI HRP and CBS7435 HRP (Figure 1B, C and 1D, respectively) showed 2-fold longer adaptation times compared to the KM71H strain (see also Table 2). These results clearly show that recombinant Pichia strains, which heterologously produce proteins upon the presence of the inducer methanol, carry an additional biological burden which significantly changes their metabolism and slows down their adaptation to methanol.

Table 2 Batch experiments with methanol pulses to determine strain specific parameters of different P. pastoris strains.

Surprisingly, the P. pastoris Mut+ strain SMD1168H GalOX showed the longest adaptation time of 14 h (Figure 1E). However, this maximum was just detected as a result of methanol depletion and thus Monod kinetics, rather than representing the real time point of full methanol adaptation. As shown in Figure 1, the Mut+ strain SMD1168H GalOX showed a very different metabolic behaviour with a flat slope for qCO2 compared to the qCO2 curves of the MutS strains, which might be due to intracellular regulation and inhibition events, caused by produced H2O2 and the absence of sufficient catalases in the cells. Regulation events like this have been described in various systems before [2224]. Thus, the determined Δtimeadapt of 14 h for the Mut+ strain is questionable and the cells had probably been adapted to methanol before.

Consequently, the usefulness of the strain characteristic parameter Δtimeadapt, especially with regard to the Mut+ strain SMD1168H GalOX, was checked by introducing a time derivative of the qCO2 signal (qCO2'). Since qCO2 and the specific growth rate μ are linearly related, the maximum of qCO2' represents the time point of adaptation of the culture to the new substrate methanol. At the maximum qCO2' the cells should be already fit for further assimilation to methanol, and thus this point represents a good starting point for consecutive fed batch cultivations without running the risk of methanol accumulation. A good example for this is shown in Figure 1A for the strain KM71H, where a rising slope is shown with a single qCO2' maximum after 2.1 h. However, the other MutS strains tested (KM71H HRP, KM71H PDI HRP and CBS7435 HRP) were characterized by a more bumpy qCO2 curve, resulting in a qCO2' signal with several shoulders (Figure 1B, C and 1D, respectively). Apparently, the adaptation of these strains to methanol did not happen as straight-forward as for the strain KM71H, but with local minima and maxima of the metabolic capacity probably caused by regulatory events upon an excess of methanol. For the Mut+ strain SMD1168H GalOX the specific carbon dioxide production was rather low compared to the MutS strains. The maximum in qCO2' was determined already after 2.7 h, which was by far sooner than the observed Δtimeadapt of 14 h. The quite constant qCO2' signal over time (Figure 1E), further supports the fact that the cell capacity was already adapted to its maximum after this short time of 2.7 h and that the cells had been fully adapted to methanol much sooner than the Δtimeadapt of 14 h.

In general, the more detailed analysis of the different strains in their adaptation to methanol revealed three different patterns in the qCO2' signal: 1) with a single maximum, 2) with several shoulders and 3) a quite constant signal over time. This underlines the necessity for individual analyses of different strains in order to quantitatively characterize them during the adaptation phase in the presence of methanol excess. Of course, the observed maximum in off-gas activity is dependent on the affinity of the cells to the substrate methanol. The maximum in qCO2 could also be reached in terms of substrate limitation and a consequent drop in the qCO2 signal due to Monod kinetics rather than by the maximum metabolic adaptation to methanol, which in this study can clearly be seen for the Mut+ strain SMD1186H GalOX (Figure 1E). However, for all strains tested in this study, maximum values of qCO2' were reached before the maximum off-gas activity (Figure 1 and Table 2), demonstrating that the applied concentration of methanol in the adaptation pulse was high enough to guarantee that the maximum qCO2' was reached independent of Monod kinetic effects. The validity of qCO2' as a reliable signal to detect the adaptation of the culture to methanol is further underlined when analyzing the respiratory quotient (RQ) during the adaptation pulse, which is exemplarily shown for the strain KM71H in Figure 2. During the adaptation pulse, RQ fluctuates with local minima and maxima until the signal becomes rather constant indicating the adaptation of the culture to methanol, which actually coincides with the maximum of qCO2'. Bespoken fluctuations of RQ at the beginning of the adaptation pulse represent the differences in catabolic and anabolic activity of the adapting cells. Similar effects have been observed for Saccharomyces cerevisiae [22, 25, 26]. During the following pulses, RQ shows a rather constant signal indicating that the cells had already been adapted to the new substrate. These findings validate the parameter qCO2' as a reliable indicator for methanol adaptation.

Figure 2
figure 2

Batch cultivation of the strain KM71H with a 0.5% (v/v) methanol pulse for adaptation and 1% (v/v) methanol pulses for q s max determination. Straight line, calculated carbon dioxide evolution rate (CER); dotted line, respiratory quotient RQ (CER/OUR); black dot, specific substrate uptake rate (qs) for each pulse.

However, despite the advantage of describing the time point of adaptation of the culture to methanol more accurately, the use of qCO2' as a parameter to determine the starting point of the following fed batch could be risky because of the described fluctuations in qCO2 (Figure 1). On the other hand, the maximum off-gas activity (Δtimeadapt) is a parameter which safely describes methanol adaptation. Another advantage of using Δtimeadapt instead of qCO2' is the fact that due to no significant biomass increase during the adaptation pulses, the carbon dioxide evolution rate (CER), which can easily be derived in online mode, can be used to determine Δtimeadapt, as also shown in our previous study [18], and thus describes a valuable online tool for process monitoring and control. Consequently, we stuck to Δtimeadapt as a minimum and safe parameter for complete methanol adaptation, as we have done previously [18], while the maximum of qCO2' should be regarded as a possible minimum prerequisite to start the fed batch feed.

Determination of the specific substrate uptake rates (qs adapt and qs max) using batch cultivations with repeated methanol pulses

The frequent determination of biomass and methanol concentrations allowed specific rate calculations for methanol uptake (qs) during the methanol pulses. Adaptation pulses of 0.5% (v/v) methanol were used to determine the parameter qs adapt (specific substrate uptake rate during the adaptation pulse) for each strain (Table 2). The knowledge of qs adapt allows the operator to adjust a specific, optimal flow of methanol during the adaptation for each single strain, thus preventing methanol accumulation.

After the methanol of the adaptation pulse was depleted, several pulses with 1% (v/v) methanol were conducted to determine the maximum specific substrate uptake rate (qs max) for each P. pastoris strain. In Figure 2, this strategy is exemplarily shown for the MutS strain KM71H, for which the specific substrate uptake rate was calculated with 0.96 ± 0.09 mmol·g-1·h-1 in the adaptation period (qs adapt) and with 1.94 ± 0.19 mmol·g-1·h-1 as a maximum during pulses (Figure 2 and Table 2). All P. pastoris strains were characterized with the above mentioned strategy and results are summarized in Table 2. Very different values for the single strain characteristic parameters of the different P. pastoris strains were determined and it becomes evident, that P. pastoris strains require specific conditions for an optimal adaptation to methanol and that the maximum levels of methanol uptake differ significantly between the single strains.

Summing up, our results underline the importance of analyzing different P. pastoris strains in more detail during methanol adaptation and regarding the maximum substrate uptake rate. By a more detailed analysis of each strain, which can be done in a fast way by the strategy described here, the bioprocess-relevant strain characteristic parameters Δtimeadapt, qs adapt and qs max can be easily extracted and used for setting up a consecutive fed batch experiment based on qs. Additionally, the more detailed individual analysis of each strain during the adaptation phase delivers important information for an early process development.

Fed batch design based on batch cultivations with methanol pulses

After characterization experiments in batch cultivations, fed batch experiments using qs based feeding profiles were carried out, according to our previous study [18]. Besides proving the general applicability of this feeding strategy on different P. pastoris strains, we wanted to check for parameter consistency, i.e. whether qs max, which had been determined in batch experiments before, could be reached in fed batch cultivations without observable methanol accumulation.

After a batch phase on glucose as substrate (volume 1.5 l), an exponential fed batch cultivation with glucose yielded in biomass concentrations of up to 70 g/l in a final volume of 2.5 l. As soon as glucose was depleted, a sample was taken to determine the current biomass concentration by measuring the OD600 and the DCW. Afterwards, all cultures were induced with a flow corresponding to a qs below qs adapt to guarantee a certain safety margin. The observed adaptation times during the fed batch cultivations for all strains are summarized in Table 3. All the MutS strains were characterized by slightly longer adaptation times in fed-batches compared to the batch cultivations, which might be due to the mode of providing methanol, i.e. in batch experiments methanol was pulsed into the reactor, resulting in a temporal excess of methanol, whereas in fed batch cultivations methanol was slowly fed into the bioreactor according to the actual biomass content, which is why it took the cells longer to fully adapt to the new substrate methanol.

Table 3 Results of dynamic fed batch cultivations of different P. pastoris strains based on qs.

However, the Mut+ strain SMD1168H GalOX showed a completely different behaviour: the adaptation time in the fed batch was much shorter than the one observed in the batch experiment, which might be due to the fact that the sudden excess of methanol in the batch-pulse experiment resulted in inhibition events caused by produced H2O2 and thus Δtimeadapt, which had been determined in the batch experiment, did not really describe the time point of adaptation of the culture to methanol, but rather qCO2' (vide supra).

However, as soon as the maximum in off-gas activity was reached in the fed batch cultivations, the feed was increased to 1.0 mmol·g-1·h-1 and then stepwise (0.5 mmol·g-1·h-1/step) every 24 hours up to qs max. Since qs max for the strain KM71H HRP PDI had been determined with just 1.08 mmol·g-1·h-1, smaller increases of 0.25 mmol·g-1·h-1 in the qs steps were performed. This strategy, describing a stepwise increase of qs to qs max, was chosen to allow the detection of possible dependencies between qs and qp. At several time points during each step, samples were taken and, based on the apparent biomass content (estimated by OD600 measurements) feeding rates were adjusted manually corresponding to the defined qs set point. Since these regular adjustments of qs to the actual biomass content were performed, the feeding profile actually represented an accelerated exponential feeding profile, which has proven to result in higher specific productivities compared to other feeding profiles tested [18].

In Figure 3, this fed batch strategy, which was applied for all P. pastoris strains in this study, is exemplarily shown for the Muts strain KM71H and the Mut+ strain SMD1186H GalOX. The qs set points were increased stepwise to a value of 2.00 mmol·g-1·h-1 for the KM71H strain, and, as shown in Figure 3A, methanol accumulation was only observed when the feeding rate exceeded values above the respective qs max of 1.94 ± 0.19 mmol·g-1·h-1. When we stopped feeding after ~ 70 h, the accumulated methanol was consumed immediately. For the strain SMD1168H GalOX, the qs set point was stepwise increased to 2.5 mmol·g-1·h-1 (qs max = 2.62 mmol·g-1·h-1) and no methanol accumulation was detected when these qs steps were conducted (Figure 3B). The same dynamic feeding strategy was applied to the other P. pastoris strains and all essential results are summarized in Table 3.

Figure 3
figure 3

Fed batch cultivations of different P. pastoris strains on methanol with a stepwise increase of q s to q s max . A, P. pastoris Muts strain KM71H; B, P. pastoris Mut+ strain SMD1186H GalOX; straight line, set point for qs; black dot, calculated qs values; black triangle, methanol concentration in the supernatant.

As shown in Table 3, values for qs max, which had been determined in batch pulsing experiments before, were reached in fed batch experiments without methanol accumulation. However, when qs set points were further increased, methanol accumulation was observed. This proves that the values for qs max from batch experiments can be found again in consecutive fed batch experiments for different P. pastoris strains and thus shows the great potential of this approach: by determining qs max in fast and easy-to-do batch experiments, the operator does not run the risk of overfeeding methanol in consecutive fed batch cultivations. However, we recommend a maximum feed flow below qs max to guarantee a certain safety margin.

Regarding the productivity, for all the strains, except for SMD1168H GalOX, same or even higher specific productivities were obtained in fed batch cultures compared to the batch experiments, signifying that volumetric productivities were higher in fed batch cultivations compared to batch experiments (Table 2 and Table 3). This shows that our approach employing dynamic feeding profiles [18] can be successfully applied on different P. pastoris strains. The lower qp max of the strain SMD1168H GalOX in fed batch cultivation probably resulted from copper-limitation [20], but was not investigated any further, since it was not goal of this study to optimize the production of recombinant enzymes of the single P. pastoris strains.

Summing up, in this study we show the general applicability of a fast approach to determine certain strain characteristic parameters, which were extracted out of batch experiments and were verified in subsequent fed batch cultures of different P. pastoris strains, making this approach a valuable tool for fast bioprocess development.


In the present study we prove that the fast approach to determine bioprocess-relevant strain characteristic parameters and the novel dynamic feeding strategy based on qs, which we have described recently for one recombinant P. pastoris strain [18], are applicable for a variety of P. pastoris strains with different phenotypes producing different recombinant proteins. This underlines the great potential of this strategy as a fast and simple tool to quantify a minimal set of parameters needed to set up consecutive fed batch regimes, which is particularly important for industry, where a fast process development is essential.

Our strategy describes:

  1. 1.

    a batch experiment with

  • ♦ a 0.5% (v/v) methanol adaptation pulse to determine Δtimeadapt and qs adapt

  • ♦ at least 4 consecutive 1.0% (v/v) methanol pulses to determine qs max

  1. 2.

    a dynamic fed batch feeding strategy based on qs, where after the adaptation of the culture to methanol (described by a maximum in off-gas activity), qs set points can be increased to qs max without observable methanol accumulation

We further show that a detailed analysis of the adaptation to methanol reveals a variability of adaptation characteristics of the different strains, highlighting that an individual analysis of potentially new strains in this respect is required to allow quantitative strain characterization and to derive parameters necessary for a consecutive fed batch set-up. The parameter Δtimeadapt safely describes the transition condition during methanol adaptation. Since also the carbon dioxide evolution rate (CER), as well as the oxygen uptake rate (OUR), can be used to determine Δtimeadapt, it further describes an online available data source allowing real-time monitoring and controlling of bioprocesses, which is essential under the aspect of Process Analytical Technology (PAT).

In this study we show that easy-to-do batch experiments with methanol pulses delivered valid and safe strain characteristic parameters which were consistent and precise enough to set up fed batch feeding profiles based on the specific substrate uptake rate. Our strategy is faster than the usually used continuous cultures or consecutive fed batch cultivations, and therefore allows faster process development. Besides, the strategy described here can be carried out using standard equipment without the need of cost-intensive tools; only a standard bioreactor connected to an off-gas analysis system and an HPLC or GC to determine methanol concentrations are required to carry out the experiments.



mtime for adaptation of the culture to the new substrate (methanol) [h]


specific growth rate [h-1]


2, 2' azino bis 3-ethylbenzthiazoline-6-sulphonic acid


CER- carbon dioxide evolution rate [mmol·l-1·h-1]


derivative of qCO2


dry cell weight


high performance liquid chromatography


horseradish peroxidise


galactose 6-oxidase


methanol utilization phenotype




methanol utilization slow phenotype


optical density at 600 nm [AU]


proportional-integrative-derivative controller


specific carbon dioxide production rate [mmol·g-1·h-1]


derivative of qCO2


specific productivity [U·g-1·h-1]

qs adapt:

specific substrate uptake rate during adaptation [mmol·g-1·h-1]

qs max:

maximum specific substrate uptake rate [mmol·g-1·h-1]


specific substrate uptake rate [mmol·g-1·h-1]


rounds per minute


respiratory quotient


volume gas flow per volume liquid per minute


  1. d'Anjou MC, Daugulis AJ: A rational approach to improving productivity in recombinant Pichia pastoris fermentation. Biotechnology and bioengineering. 2001, 72: 1-11. 10.1002/1097-0290(20010105)72:1<1::AID-BIT1>3.0.CO;2-T 10.1002/1097-0290(20010105)72:1<1::AID-BIT1>3.0.CO;2-T

    Article  Google Scholar 

  2. Curvers S, Linnemann J, Klauser T, Wandrey C, Takors R: Recombinant protein production with Pichia pastoris in continuous fermentation - kinetic analysis of growth and product formation. Engineering in Life Sciences. 2002, 2: 229-235. 10.1002/1618-2863(20020806)2:8<229::AID-ELSC229>3.0.CO;2-9. 10.1002/1618-2863(20020806)2:8<229::AID-ELSC229>3.0.CO;2-9

    Article  Google Scholar 

  3. Ren H, Yuan J: Model-based specific growth rate control for Pichia pastoris to improve recombinant protein production. Journal of Chemical Technology and Biotechnology. 2005, 80: 1268-1272. 10.1002/jctb.1321. 10.1002/jctb.1321

    Article  Google Scholar 

  4. Ren HT, Yuan JQ, Bellgardt KH: Macrokinetic model for methylotrophic Pichia pastoris based on stoichiometric balance. Journal of Biotechnology. 2003, 106: 53-68. 10.1016/j.jbiotec.2003.08.003

    Article  Google Scholar 

  5. Invitrogen: Pichia Fermentation Process Guidelines. 2002,

    Google Scholar 

  6. Jahic M, Rotticci-Mulder JC, Martinelle M, Hult K, Enfors SO: Modeling of growth and energy metabolism of Pichia pastoris producing a fusion protein. Bioprocess and Biosystems Engineering. 2002, 24: 385-393. 10.1007/s00449-001-0274-5. 10.1007/s00449-001-0274-5

    Article  Google Scholar 

  7. Trinh LB, Phue JN, Shiloach J: Effect of methanol feeding strategies on production and yield of recombinant mouse endostatin from Pichia pastoris. Biotechnology and Bioengineering. 2003, 82: 438-444. 10.1002/bit.10587

    Article  Google Scholar 

  8. Zhang W, Bevins MA, Plantz BA, Smith LA, Meagher MM: Modeling Pichia pastoris growth on methanol and optimizing the production of a recombinant protein, the heavy-chain fragment C of botulinum neurotoxin, serotype A. Biotechnology and Bioengineering. 2000, 70: 1-8. 10.1002/1097-0290(20001005)70:1<1::AID-BIT1>3.0.CO;2-Y

    Article  Google Scholar 

  9. Zhang W, Sinha J, Smith LA, Inan M, Meagher MM: Maximization of Production of Secreted Recombinant Proteins in Pichia pastoris Fed-Batch Fermentation. Biotechnology Progress. 2005, 21: 386-393.

    Article  Google Scholar 

  10. Plantz BA, Sinha J, Villarete L, Nickerson KW, Schlegel VL: Pichia pastoris fermentation optimization: energy state and testing a growth-associated model. Applied Microbiology and Biotechnology. 2006, 72: 297-305. 10.1007/s00253-005-0271-8

    Article  Google Scholar 

  11. Potgieter TI, Kersey SD, Mallem MR, Nylen AC, d'Anjou M: Antibody expression kinetics in glycoengineered Pichia pastoris. Biotechnology and Bioengineering. 2010, 106: 918-927. 10.1002/bit.22756

    Article  Google Scholar 

  12. Kobayashi K, Kuwae S, Ohya T, Ohda T, Ohyama M, Tomomitsu K: High level secretion of recombinant human serum albumin by fed-batch fermentation of the methylotrophic yeast, Pichia pastoris, based on optimal methanol feeding strategy. Journal of Bioscience and Bioengineering. 2000, 90: 280-288.

    Article  Google Scholar 

  13. Jacobs P, Inan M, Festjens N, Haustraete J, Van Hecke A, Meagher M, Contreras R, Callewaert N: Fed-batch fermentation of GM-CSF-producing glycoengineered Pichia pastoris under controlled specific growth rate. Microbial Cell Factories. 2010, 9: 93- 10.1186/1475-2859-9-93

    Article  Google Scholar 

  14. Sinha J, Plantz BA, Zhang W, Gouthro M, Schlegel V, Liu C-P, Meagher MM: Improved production of recombinant ovine interferon-t by Mut+ strain of Pichia pastoris using an optimized methanol feed profile. Biotechnology Progress. 2003, 19: 794-802. 10.1021/bp025744q

    Article  Google Scholar 

  15. Khatri NK, Hoffmann F: Oxygen-limited control of methanol uptake for improved production of a single-chain antibody fragment with recombinant Pichia pastoris. Applied Microbiology and Biotechnology. 2006, 72: 492-498. 10.1007/s00253-005-0306-1

    Article  Google Scholar 

  16. Khatri NK, Hoffmann F: Impact of methanol concentration on secreted protein production in oxygen-limited cultures of recombinant Pichia pastoris. Biotechnology and Bioengineering. 2006, 93: 871-879. 10.1002/bit.20773

    Article  Google Scholar 

  17. Cunha AE, Clemente JJ, Gomes R, Pinto F, Thomaz M, Miranda S, Pinto R, Moosmayer D, Donner P, Carrondo MJT: Methanol induction optimization for scFv antibody fragment production in Pichia pastoris. Biotechnology and bioengineering. 2004, 86: 458-467. 10.1002/bit.20051

    Article  Google Scholar 

  18. 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. Microbial Cell Factories. 2011, 10: 14- 10.1186/1475-2859-10-14

    Article  Google Scholar 

  19. Hartner FS, Ruth C, Langenegger D, Johnson SN, Hyka P, Lin-Cereghino GP, Lin-Cereghino J, Kovar K, Cregg JM, Glieder A: Promoter library designed for fine-tuned gene expression in Pichia pastoris. Nucleic Acids Research. 2008, 36: e76/71-e76/15.

    Article  Google Scholar 

  20. Spadiut O, Olsson L, Brumer H: A comparative summary of expression systems for the recombinant production of galactose oxidase. Microb Cell Fact. 2010, 9: 68- 10.1186/1475-2859-9-68

    Article  Google Scholar 

  21. Cregg JM: Pichia Protocols, Second Edition. Methods Mol Biol. Edited by: Cregg JM. 2007, 268-Totowa, NJ, U. S.,

    Google Scholar 

  22. Herwig C, von Stockar U: A small metabolic flux model to identify transient metabolic regulations in Saccharomyces cerevisiae. Bioprocess and Biosystems Engineering. 2002, 24: 395-403. 10.1007/s00449-001-0277-2. 10.1007/s00449-001-0277-2

    Article  Google Scholar 

  23. Diderich JA, Schepper M, van Hoek P, Luttik MA, van Dijken JP, Pronk JT, Klaassen P, Boelens HF, de Mattos MJ, van Dam K, Kruckeberg AL: Glucose uptake kinetics and transcription of HXT genes in chemostat cultures of Saccharomyces cerevisiae. The Journal of biological chemistry. 1999, 274: 15350-15359. 10.1074/jbc.274.22.15350

    Article  Google Scholar 

  24. Rizzi M, Theobald U, Querfurth E, Rohrhirsch T, Baltes M, Reuss M: In vivo investigations of glucose transport in Saccharomyces cerevisiae. Biotechnology and Bioengineering. 1996, 49: 316-327.

    Article  Google Scholar 

  25. Duboc P, Von Stockar U, Villadsen J: Simple generic model for dynamic experiments with Saccharomyces cerevisiae in continuous culture: decoupling between anabolism and catabolism. Biotechnology and Bioengineering. 1998, 60: 180-189. 10.1002/(SICI)1097-0290(19981020)60:2<180::AID-BIT5>3.0.CO;2-J

    Article  Google Scholar 

  26. Duboc P, Cascao-Pereira LG, von Stockar U: Identification and control of oxidative metabolism in Saccharomyces cerevisiae during transient growth using calorimetric measurements. Biotechnology and Bioengineering. 1998, 57: 610-619. 10.1002/(SICI)1097-0290(19980305)57:5<610::AID-BIT13>3.0.CO;2-F

    Article  Google Scholar 

Download references


The authors are very grateful to Prof. Anton Glieder and his workgroup for providing all P. pastoris MutS strains expressing the enzyme horseradish peroxidase. Additionally, the authors express their gratitude to Josef Schulze, who supported the experimental work in the lab.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Christoph Herwig.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

CD designed and performed the majority of the experiments, analyzed the data and drafted the manuscript. OS designed and performed some experiments and finalized the manuscript. CH conceived the study and supervised research. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Dietzsch, C., Spadiut, O. & Herwig, C. A fast approach to determine a fed batch feeding profile for recombinant Pichia pastoris strains. Microb Cell Fact 10, 85 (2011).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Pichia pastoris
  • strain characterization
  • specific substrate uptake rate
  • batch cultivation
  • methanol pulse
  • dynamic feeding profile