Skip to main content

Modeling and measuring intracellular fluxes of secreted recombinant protein in Pichia pastoris with a novel 34S labeling procedure

Abstract

Background

The budding yeast Pichia pastoris is widely used for protein production. To determine the best suitable strategy for strain improvement, especially for high secretion, quantitative data of intracellular fluxes of recombinant protein are very important. Especially the balance between intracellular protein formation, degradation and secretion defines the major bottleneck of the production system. Because these parameters are different for unlimited growth (shake flask) and carbon-limited growth (bioreactor) conditions, they should be determined under "production like" conditions. Thus labeling procedures must be compatible with minimal production media and the usage of bioreactors. The inorganic and non-radioactive 34S labeled sodium sulfate meets both demands.

Results

We used a novel labeling method with the stable sulfur isotope 34S, administered as sodium sulfate, which is performed during chemostat culivations. The intra- and extracellular sulfur 32 to 34 ratios of purified recombinant protein, the antibody fragment Fab3H6, are measured by HPLC-ICP-MS. The kinetic model described here is necessary to calculate the kinetic parameters from sulfur ratios of consecutive samples as well as for sensitivity analysis. From the total amount of protein produced intracellularly (143.1 μg g-1 h-1 protein per yeast dry mass and time) about 58% are degraded within the cell, 35% are secreted to the exterior and 7% are inherited to the daughter cells.

Conclusions

A novel 34S labeling procedure that enables in vivo quantification of intracellular fluxes of recombinant protein under "production like" conditions is described. Subsequent sensitivity analysis of the fluxes by using MATLAB, indicate the most promising approaches for strain improvement towards increased secretion.

Background

The production of recombinant proteins in yeast has to compete with other host organisms, mainly bacteria and mammalian cell lines. Strain improvement therefore is an essential step between the discovery of a new protein and its large scale production. Yeasts like Pichia pastoris grow faster and to a higher cell density compared to mammalian cells, however the low specific productivity (the amount of secreted protein per unit biomass and time) is their major drawback [1]. A lot of efforts have already been made to find and overcome specific bottlenecks in the cellular protein production and secretory system [reviewed by [2]].

At genomic level increasing the gene copy number as well as the promoter strength leads to higher productivities [35]. The overload of the endoplasmic reticulum (ER) with recombinant protein may induce the unfolded protein response (UPR) [68] followed by enhanced ER-associated degradation (ERAD) [9, 10]. Among many other things, UPR reduces overall translation speed [11] and enforces ERAD via the Ire1 signaling cascade [12]. ERAD causes proteolytic digestion of malfolded protein in the cytosolic proteasome [13]. Thus, reduced ER-stress can be beneficial for recombinant protein production. Therefore, many attempts have been made to improve the complex process of protein maturation, mainly by co-overexpressing ER resident chaperons or foldases like BiP / Kar2, Pdi1 or calnexin [1416]. Furthermore the transport from the ER to the Golgi and finally into the exterior can be improved by co-overexpression of proteins involved in this pathways. Examples are Sso1 and Sso2, both coding for plasma membrane t-SNARE proteins [17] or Cog6, Coy1 and Bmh2, all coding for proteins involved in vesicular transport [18].

In the strain improvement process by cell engineering it is required to achieve high yields in short time. A focused and systematic approach therefore would be to identify the most important bottleneck in recombinant protein synthesis being the one which modification has the highest impact on protein titers.

Kinetic models are a valuable tool in this regard, as they give insights into intracellular fluxes. The formal kinetic description of the processing and transport of secreted proteins are already known for quite a while [19, 20]. However, the challenge is the experimental determination of the parameters needed in those models. Furthermore it is necessary to make as few assumptions as possible so that a production process can still be described. In this regard the experiments have to be done under carbon limited, production "similar", growth in bioreactors under defined and controlled conditions instead of using shake flask cultivations. This is usually not possible when labeling is performed with radioactive isotopes or when protein kinetics is measured with microscopic tools, like fluorescence microscope imaging. Handling of large volumes of radioactive material is not feasible for risk of contamination. Microscopic imaging on the other side quantifies the protein fluxes by comparing images of living cells over time [21]. The advantage is that single cells are analyzed instead of an average. However cells are exposed to non-defined conditions which are likely to be different to the bioreactor. Furthermore, this method is limited to fluorescent or fluorescent-tagged proteins.

It is important that the model of choice accounts for intracellular degradation as well, because a substantial amount of recombinant protein may be degraded via ERAD or other pathways. Also the dilution by growth, especially in the fast growing prokaryotes and yeast cells, has to be taken into account. For example at a specific growth rate μ = 0.1 h-1, 7% of the intracellular protein is inherited to the daughter cells (as described in this study). Several kinetic model studies for antibody production have been developed [2225], but intracellular protein degradation as well as protein inheritance to the daughter cells was not taken into account in these cases.

In this work we present a novel 34S labeling method during chemostat cultivation, providing data to consider intracellular protein formation, intracellular degradation, secretion and dilution by growth in the kinetic model. Further sensitivity analysis of the intracellular protein flux enables the estimation of their impact and serves as a decisions basis for further strain improvement.

Results

Structured Kinetic Model

The structured kinetic model describes the intra- and extracellular recombinant protein pools. The dynamics of the intra- and extracellular Fab3H6 pools are expressed in two separate differential equations (1) and (3), which are modified from Noe and Delenick [19] and Batt and Kompala [20].

The amount of intracellular protein (Pi) depends on the fluxes of intracellular protein formation (qPi), protein secretion (qSec), intracellular protein degradation (qDeg) and protein dilution into the daughter cells (qDil). Due to the rapid growth of yeast cells, the inheritance of protein to daughter cells has to be taken into account. The dynamics of protein secretion, degradation and dilution can be assumed to be a first-order kinetic [22, 23, 26], with specific time constants (K). Equation (1) can therefore be transformed into equation (2).

(1)
(2)

In chemostat culture, the extracellular product pool (Pe) is a function of the secreted protein (qSec) and the amount of protein in the media that is harvested (qHar) (equation 3). Also here first-order kinetic can be assumed and equation (3) transformed into equation (4).

(3)
(4)

The observed viability of more than 99.0% ensures that the dilution rate (D) and the specific growth rate (μ) are equal in these continuous cultivations and thus also the time constants of dilution KDil and harvest KHar are the same.

Calculation of the time constants for secretion KSec and degradation KDeg

Due to steady state conditions during continuous cultivation no change in the total amount of Pi or Pe occurred over time. Therefore equation (4) is zero and can be transformed into equation (5).

(5)

To determine the time constant for protein degradation KDeg, the decrease of the amount of 32S-containing protein (P32S) during the continuous labeling experiment with 34S is used. The fraction of non-labeled sodium sulfate during labeling remains always below 1%. It is therefore assumed that the intracellular 32S protein formation during labeling is negligible. Thus equation (2) can be converted to equation (6) and further simplified for q(Pi32S) (equation 7). The solution of this new differential equation can be transformed into a linear system, shown in equation (8).

(6)
(7)
(8)

Linear regression analysis was performed with the statistics software package R. The slope represents the sum of the three time constants for degradation, secretion and dilution. KDeg is calculated by subtracting the known KSec and KDil from the absolute value of the slope.

Calculation of the rate of intracellular recombinant protein formation (qPi)

Due to steady state considerations, equation (1) can be set to zero. The result is shown in equation (9), were qPi is depicted. This equation expresses that all protein that is secreted, degraded or transferred into daughter cells is produced in the cells before.

(9)

Chemostat cultivation

Extracellular and intracellular Fab3H6 concentrations remained constant over time (data not shown). The average concentrations were 91.8 μg per gram YDM (STD: 9.6%) intracellular and 507.5 μg per gram YDM extracellular (STD: 12.1%). The standard deviation (STD) derived from two independent chemostat cultivations.

Immunoprecipitation

The antibody fragment Fab3H6 was immunopurified during the continuous labeling from the fermentation supernatants and the cell biomass respectively (figure 1). In the supernatant almost only the dimer was present, whereas the cell lysate contained a substantial amount of monomers. Some impurities were detected in the immunoprecipitates (IPs) from the cell lysates. This might be a result of covalent and hydrophobic interactions of the nascent peptide chains during folding, assembly and transport. Therefore the decrease of the intracellular 32S / 34S ratio has to be corrected by subtracting the background levels derived from impurities. According to the observed half times of secretion and degradation (table 1), all intracellular recombinant protein should contain less than 0.1% 32S 140 min post labeling. Therefore it is assumed that from this time point on the measured 32S signal derived from impurities. The background was calculated by extrapolating the 32S / 34S ratios from 140 min post labeling on towards the beginning of the labeling.

Figure 1
figure 1

Immunoprecipitation of Fab3H6. Immunoprecipitation of extracellular and intracellular Fab3H6 during 34S-labeling. Lanes 2-6: consecutive extracellular samples, lanes 7-11: consecutive intracellular samples. Immunoprecipitates were separated by non-reduced SDS-PAGE and silver stained. Dimers and monomers of the Fab fragment are indicated by arrows. The protein ladder (Page Ruler, Invitrogen) is seen in lane 1.

Table 1 Intracellular fluxes of the antibody fragment Fab3H6

Sulfur isotopic distribution during labeling

The 32S amount of compounds smaller than 3 kDa, mainly sulfate, in the media is shown in figure 2A (open circles). The values are always around or below 1%, which is a prerequisite for the kinetic model (equation 7). The decrease of the unlabeled antibody fragment Fab3H6 from two replicative, independent chemostat cultivations is shown in figure 2A and 2B. The amount of intracellular non-labeled 32S recombinant protein is always below the extracellular value. This is not surprising because labeled protein has to be produced within the cell before being released to the media. In some samples the signal to noise ratio was not sufficient for ICP-MS measurement so that these were not taken into account. Therefore, the number of data points in figure 2A and figure 2B are not identical.

Figure 2
figure 2

34S labeling in chemostat culture. (A + C) and (B + D) respectively derive from two independent, replicate fermentations. In (A) and (B) the decrease of the 32S content of intra- and extracellular Fab3H6 during continuous 34S labeling in chemostat cultivation is shown. Additionally, in (A) the 32S content of small compounds, mainly sulfate, in the media was measured (open circles). The linear regression analysis of the normalized 32S contents of the intracellular Fab3H6 (equation 8) is presented in (C) and (D). The dashed lines illustrate the 95% confidence intervals. The slope represents the sum of the time constants of protein secretion, degradation and dilution by growth.

Intracellular Fab3H6 fluxes

The dilution rate of the chemostat defines the growth rate of the culture. The specific growth rate expressed in minutes equates the time constant of the dilution by growth KDil = 0.00167 min-1 (table 1). The corresponding half time t1/2 Dil = 415.9 min is calculated by equation (10):

(10)

The relative standard deviation (STD) of 0.05% results from the discontinuous harvest (see material an methods). The rate by which the protein is inherited to the daughter cells qDil = 9.18 μg g-1 h-1 (STD: 9.6%).

For the determination of KSec the quotient from extra- and intracellular protein concentration is multiplied with KDil (equation 5). KSec was calculated to be 0.0092 min-1 (STD: 2.5%) with the corresponding half time t1/2 Sec = 75.3 minutes (STD: 2.5%) and the secretion rate qSec = 50.7 μg g-1 h-1 (STD: 12.1%) (table 1).

Before the degradation rate can be calculated, the intracellular values of the unlabeled 32S Fab3H6 content (see figure 2B) have to be transformed according to equation (8). This is done by taking the logarithm of the normalized data points. The linear regression analyses of the two biological replicates are shown in figure 2C and figure 2D respectively, with the transformed values on the y-axis and the time on the x-axis. KDeg is calculated by subtracting KSec and KDil from the positive slopes 0.0240 min-1 and 0.0258 min-1. The time constant for intracellular degradation KDeg = 0.0151 min-1 (STD 2.8%) with the corresponding t1/2 Deg = 45.8 min (STD 2.8%) and qDeg = 83.3 μg g-1 h-1 (STD: 9.3%). Furthermore, the 95% confidence interval is represented by the dashed line in figure 2D.

The rate of intracellular protein formation qPi is simply the sum of qDil, qSec and qDeg and is 143.1 μg g-1 h-1 (STD: 8.8%).

Figure 3 gives an overview of the rates, where qPi of 143.1 μg g-1 h-1 represent 100%. The other rates (qDeg, qSec and qDil) are shown relative to the intracellular protein formation. From 100% being produced within the cell, 58% are degraded intracellularly, 35% are secreted to the exterior and 7% are inherited to daughter cells.

Figure 3
figure 3

Overview of the intracellular fluxes of the recombinant antibody fragment Fab3H6 in P. pastoris.

Sensitivity analysis

The sensitivity analysis reflecting the impact of the different model parameters to the extracellular Fab3H6 concentration is shown in figure 4. Parameters, describing either protein formation (qPi), secretion (t1/2 Sec) or degradation (t1/2 Deg) are continuously changed in MATLAB from their starting point up to +/- 10 fold. Each fold increase or decrease is pictured as a solid line. In the three dimensional plots, two of the model parameters are plotted against the corresponding Fab3H6 titers, resulting in the presented areas (figure 4A-D). The 10 fold improvement of each parameter is shown as a white arrow.

Figure 4
figure 4

Sensitivity analysis of the extracellular Fab3H6 concentration. In each picture the effect of two different input variables on the extracellular Fab3H6 concentration is shown. The red dot in the middle represents the measured values in the chemostat. These initial values were changed from -10 to +10 fold resulting in depicted areas. In (A) t1/2 Deg and t1/2 Sec, in (B) t1/2 Sec and qPi, in (C) and (D) t1/2 Deg and qPi were analyzed. Furthermore, in (D) the maximum secretion capacity was set to 101.4 μg g-1 h-1 (2 fold of the measured value).

A 10 fold "faster" secretion, t1/2 Sec = 7.5 min, would increase extracellular Fab3H6 concentration by the factor 2.5. A 10 fold "improved" degradation, t1/2 Deg = 458 min, results in a 2.0 fold increase. By combining both effects the titers are estimated to be 3.5 fold higher (figure 4A).

The intracellular protein formation rate qpi has an almost linear effect on the extracellular protein concentration (figure 4B and figure 4C). By increasing qPi by the factor of 10, from 143.1 μg g-1 h-1 to 1.43 mg g-1 h-1, the predicted Fab3H6 concentrations are also enhanced by approximately the same factor to about 5000 μg g-1.

In figure 4A-C no limit in the secretion pathway was anticipated in the underling kinetic model. This assumption however is not realistic, especially in yeast cells. It has been reported many times that the secretory capacity can be a bottleneck during recombinant protein production [2, 4, 14, 18]. Therefore, an additional model constraint was implemented in figure 4D, where the maximum secretion rate was set to be twice the default value. As a result the Fab3H6 titers reach a plateau that can be achieved by improving either t1/2 Deg or qPi. However the efforts by which the two parameters have to be varied are different, the degradation has to be changed by the factor of 4-6 and the intracellular protein formation by the factor 2.

Discussion

Kinetic models describing the cellular fluxes of recombinant secreted proteins are available as well as their mathematic solutions [19, 20]. However, in their practical application these models have usually been simplified by omission of parameters like dilution by growth as well as intracellular degradation. In our experiments we showed that qDeg and qDil are of substantial dimensions, together being 65% of all intracellular synthesized protein. Therefore, especially if a model should be applied for further strain improvement applications, it appears to be essential that those parameters are considered in the model at least in microbial expression systems.

Growth conditions like specific growth rate, media composition or oxygen limitations, just to name some, affect recombinant protein production and thus alter its intracellular flux [27, 28] So, it is necessary that the fluxes are measured under defined and reproducible conditions. In our case this is done in carbon limited chemostat cultivations using chemically defined minimal medium. Labeling with the non-radioactive and inorganic sulfur-isotope 34S, administered as sodium sulfate, enables the usage of bioreactors and a standard minimal production medium without amino acids. The common protein labeling procedures with radioactively labeled amino acids like 35S and 3H do not meet these demands.

Continuous labeling does not require any washing and centrifugation steps of the biomass and thus avoids the physiologically undefined conditions during this procedure. Non-labeled 32S and labeled 34S are measured simultaneously and thus the continuous labeling can be also interpreted as a chase of 32S. Therefore the same amount of information like with a pulse-chase experiment using common 35S amino acid labeling can be obtained. However, the detection limits, mainly depending on the resolution of the mass spectrometry, might be higher.

For efficient labeling, all sources of sulfur "contaminations" have to be avoided. This means that the sulfur in the media, and the trace salt solution within, has to be reduced to only one source, the labeled sodium sulfate. In our case the necessary changes were minimal, just the anions of the salts had to be adapted. The remaining sulfate concentration in the culture broth was quantified by nephelometric analysis with BaCl2. It contained about 10% of the sulfate of the feed medium. The amount of pulsed sodium sulfate has to be chosen due to labeling and osmolaric considerations. The non-labeled sulfate should be low and the change in osmolarity should not affect the protein fluxes. With the addition of 600 mg labeled sodium sulfate, the 32S concentration was below 1% throughout the experiment and the osmolarity increase was not detectable.

With this experimental setup the time constants and the corresponding half times of intracellular formation, intracellular degradation, dilution by growth as well as for secretion of the product can be determined. In the case of the recombinant secreted Fab3H6 35% of the protein is secreted wereas the majority (58%) is intracellularly degraded. The half time of protein degradation is 45.8 minutes, which is similar to the average half time of S. cerevisiae proteins, being 43 minutes [26]. The half time of secretion of 75.3 minutes seems to be rather slow, compared to e.g. carboxypeptidase in yeast [29]. However, the antibody fragment Fab3H6 is a heterodimer and thus could be more difficult to assembly and to secrete. In NS0 cells, Yee et al. [30] reported a half-time of one hour for IgG secretion for the major fraction of (70%) of IgG molecules. The remaining amount had a half-time close to the doubling time of NS0.

The major application of the model should be the support of a systematic strain improvement process. Therefore three dimensional sensitivity analyses were performed (figure 4). The impact of each parameter on overall expression titers as well as their synergistic effect was estimated. A ten fold improvement of the degradation or secretion half time resulted in two to three fold increased protein titers, respectively. Similar fold changes are often achieved when protein folding is engineered, like via PDI1 or BiP / KAR2 over-expression [14, 15]. The protein synthesis rate, however, affects the secreted protein titer in a linear manner. By increasing the gene copy numbers, ten fold or higher protein production rates have already been achieved [3, 31]. The linear correlation between gene copy numbers and expression levels may only be valid as long as no bottlenecks in folding or secretion occur [4]. Figure 4D takes a limitation in the secretory pathway into account. As a result protein titers reach a plateau, where further increase in the protein formation rate has no more effect. In such situations cell engineering of the secretory pathway may be the only way to break through the plateau by opening the secretion bottleneck.

Sensitivity analysis gives hints which parameter is most worthy to be modified in this regard. However strain development is an iterative process and it is likely that the change of one parameter also varies the others. Therefore after each engineering step the new fluxes should be determined. So we believe that the continuous 34S labeling described in this work may be a valuable tool for systematic strain improvement processes and deeper understanding of large scale production in bioreactors.

Methods

MATLAB implementation of the model and sensitivity analysis

The implementation of the model was done in MATLAB (additional file 1: MATLAB implementation of the kinetic model), and in its framework it is based on the work by Bibila et al. [22, 32]. The model is formulated as set of ordinary differential equations, which are then solved over time using an ODE solver. For our secretion model we used an implicit linear multistep solver (MATLAB ode15s), because it is more appropriate for chemical or biochemical problems than an explicit Runge-Kutta pair solver [33].

To evaluate the behavior of the model concerning the input parameters, namely degradation, secretion and intracellular protein formation, always two of these parameters were varied against each other. The values were then plotted against the extracellular Fab3H6 concentrations in a three dimensional representation. The three dimensional plots were created using the MATLAB functions meshgrid and surf (additional file 2: Sensitivity analysis of the kinetic parameters).

Yeast strain

The P. pastoris strain X-33 used in this study expressed the antibody fragment Fab3H6, previously described by Baumann et al. [27] and Dragosits et al. [34]. Both antibody chains are under the control of the constitutive GAP-promoter and are secreted via the Saccharomyes cerevisiae α-mating factor secretion signal. Fab3H6 is the anti-idiotypic antibody of the HIV neutralizing antibody 2F5. It has a molecular weight of 47.38 kDa and has 16 sulfur containing amino acids, 11 cysteines and 5 methionines [35, 36].

Chemostat cultivation

A preculture was incubated at 28°C for 24 h and 180 rpm on YPG (per liter: 10 g yeast extract, 20 g peptone, 10 g glycerol). The culture was harvested by centrifugation, resuspended in 50 ml sterile batch medium and used to inoculate 1.0 L batch medium in the bioreactor (Minifors, Infors, Switzerland) to a starting optical density (OD600) of 1.0.

After a batch phase of approximately 24 hours the continuous culture was started at a dilution rate of D = 0.1 h-1 with a corresponding feed medium and harvest flow rate of 100 g h-1. Cultivation conditions were controlled constantly, the temperature at 25°C, pH at 5.0 with 25% ammonium hydroxide and pO2 at 20% by controlling the stirrer speed between 600 and 1200 rpm. Air flow was kept constant at 1.5 vvm (volume gas per volume medium and minute).

The batch medium contained per liter: 40 g glycerol, 2.0 g citric acid, 12.6 g (NH)2HPO4, 0.5 g MgSO4 • 7 H2O, 0.9 g KCl, 0.022 g CaCl2 • 2 H2O, 2 ml biotin stock solution (0.2 g L-1) and 4.6 ml PTM1 trace salt stock solution. The pH was set to 5.0 with 25% HCl. The PTM1 trace salt stock solution contained per liter: 65.0 g FeSO4 • 7 H2O, 20.0 g ZnCl2, 6.0 g CuSO4 • 5 H2O, 3.36 g MnSO4 • 1 H2O, 0.82 g CoCl2 • 6 H2O, 0.2 g Na2MoO4 • 2 H2O, 0.08 g NaI, 0.02 g H3BO3 and 5 ml H2SO4 (95 - 98%).

In the chemostat medium sodium sulfate was used as the only sulfur source. Per liter this medium contained 1.0 g citric acid monohydrate, 55 g glucose • 1 H2O, 9.83 g (NH4)2HPO4, 0.41 g MgCl2 • 6 H2O, 0.29 g Na2SO4, 1.7 g KCl, 0.01 g CaCl2 • 2 H2O, 2.0 ml biotin stock solution (0.2 g L-1) and 1.6 g PTM2 trace salt stock solution. The trace salt stock solution PTM2 contained (per liter): 63.3 g FeCl2 • 6 H2O, 20.0 g ZnCl2, 5.77 g CuCl2 2H2O, 3.94 g MnCl2 • 4 H2O, 0.82 g CoCl2 • 6 H2O, 0.2 g Na2MoO4 • 2 H2O, 0.08 g NaI, 0.02 g H3BO3 and 5 ml HCl (32%).

For the continuous labeling enriched 34S sodium sulfate (isotopic distribution: < 0.1% 32S, 1.1% 33S, 98.8% 34S and < 0.05% 36S) from Isoflex USA was used in the chemostat medium.

34S labeling

Cells were grown for at least 5 resident times in chemostat to ensure steady state conditions. Continuous 34S labeling was started by changing the feed to the 34S enriched medium. In addition, at the same time, a labeled sodium sulfate pulse, 5 mL sterile solution containing in total 600 mg of 34S-labeled sodium sulfate, was administered into the bioreactor. The change in the 32S to 34S ratio of the intra- and extracellular Fab3H6 was followed for 8 hours.

Sampling

Samples were taken to determine the yeast cell dry mass, the extra- and intracellular Fab3H6 concentration and for immunoprecipitation of extra- and intracellular Fab3H6. During the 34S labeling samples were taken up to 6 times per hour. In this cases sample volume has to be kept small and therefore no biomass analysis was performed thereof.

The dead volume of the harvest port is 5 mL and therefore the first 5 mL culture broth were withdrawn. 10 mL culture broth were used for yeast dry mass (YDM) determination. For intracellular Fab3H6 measurements four cell pellets of 2 mL culture (0.05 g YDM each in capped screw tube, Biozym) were collected (1 min centrifugation at 4°C and 13.000 rpm, followed by quick freezing in liquid nitrogen). The supernatant was used for all extracellular measurements.

To enable the necessary sampling volume of at least 13 mL, the continuous harvest was replaced by discontinuous sampling. The sample volumes taken out from the bioreactor were exactly the volumes that should have been harvested by the pump. This caused a slight variation of the preset dilution rate of 0.1 h-1. However, the calculated dilution rates are in the range of minimum 0.0998 h-1 and maximum 0.1013 h-1 (calculation not shown) and it can be assumed that these changes have no significant influence.

Mechanical cell lysis of P. pastoris

The cell pellets, containing 0.05 g YDM in capped screwing tubes, were washed with 0.5 mL PBS (per liter: 8.0 g NaCl, 0.2 g KCl, 1.8 g Na2HPO4 • 2 H2O, 0.24 g KH2PO4) and further resuspended in 0.5 mL lysis buffer. 0.5 mL of glass beads (acid washed, 0.4-0.6 mm, Satorius) were added. Cells were mechanically disrupted by using the FastPrep system (MP Biomedicals; settings: 3 times 20 s shaking at 6.5 m s-1). At the bottom of the tubes small holes were pierced with a hot needle and the tubes were put onto a 2 mL eppendorf tube. Cell lysates were collected by moderate centrifugation (1 min with 1.000 g). Cells debris and glass beads were washed with 0.5 mL lysis buffer followed by the same moderate centrifugation step. Lysates were cleared by centrifugation (15 min at 13.000 g) and the supernatants were taken for further analysis (Fab immunoprecipitation or quantification). During the whole procedure samples were kept at maximum 4°C.

Lysis buffer: 1% (w/v) triton X-100, 50 mM Tris·HCl pH = 7.4, 300 mM NaCl, 5 mM EDTA and 0.02% (w/v) sodium azide. Immediately before use inhibitors were added: 1 tablet of protease inhibitor (Sigma, S8820) per 20 mL buffer, proteasome inhibitor MG-132 to a final concentration of 5 μM and lysosomal inhibitor chloroquine to a final concentration of 50 μM.

Fab3H6 immunoprecipitation

The procedure was adapted from the Current Protocols in Molecular Biology, chapter 10.16 "Immunoprecipitation" [37]. 80 μL of the anti-human IgG agarose suspension (Sigma-Aldrich) were added to 1 mL cleared cell lysate or 1 mL of culture supernatant and incubated for 2 hours at 4°C in a tube rotator. The agarose slurry was transferred onto PVDF membranes (Ultrafree®-MC Centrifugal Filter Units, 0.45 μm, Millipore) and washed four times with 0.5 mL ice-cold washing buffer and twice with 0.5 mL ice-cold PBS (centrifuged each time 5 sec at 3.000 g). To disaggregate the proteins from the matrix, the agarose suspension was incubated two times for 10 minutes with 120 μL elution buffer at room temperature. The eluates were collected by centrifugation (5 sec at 3.000 g). The buffer was changed to PBS and the solution was concentrated by the factor 10 by using the 10 kDa Amicon® Ultra-0.5 centrifugal filter devices (Millipore). The immunopurified Fab3H6 solutions were used for sulfur isotope determination (see below).

Wash buffer: 0.1% (w/v) triton X-100, 50 mM Tris·HCl pH = 7.4, 300 mM NaCl, 5 mM EDTA and 0.02% (w/v) sodium azide. Immediately before use 0.1% sodium deoxycholate was added.

Elution buffer: 6 M guanidine hydrochloride, 100 mM Tris·HCl pH = 8.5, 5 mM EDTA and 0.02% (w/v) sodium azide.

Sulfur isotope ratio determination

An inert titanium HPLC gradient system (Rheos 2000, Flux Instruments AG, Basel, Switzerland) with a metal-free autosampler (HTC PAL Autosampler, Thermo Fisher Scientific Inc., Waltham, USA) was used in combination with a high resolution inductively coupled plasma sector field mass spectrometer, ICP-SFMS (Element 2, Thermo Scientific Inc., Bremen, Germany). Sample introduction system consisted of a nebulizer (PFA-ST, Elemental Scientific Inc., Omaha, Nebrasca, USA) and a cooled (5°C) cyclonic silica glass spray chamber (PC3, ISA Elemental Scientific). Measured isotopes were 34S and 32S. Mass resolution was set to 4000. 32S was used as lock mass during measurement for instrumental mass drift correction. Dwell time per isotope was 0.1 sec. For separation by size exclusion chromatography a KW402.51E column (Shodex, Showa Denko K. K., Kawasaki, Japan) was used. Column dimensions were 1 × 150 mm. Separation was carried out under native conditions, with a 50 mM ammonium acetate, pH 6 eluent. The SEC flow was 50 μL min-1; injection volume was set to 2 μL. The integration of all chromatographic data from SEC-ICP-SFMS analysis was carried out using Chromeleon software (Version 6.7, Dionex, Sunnyvale, California, USA). Fab monomer and heterodimer dimer coeluted from the SEC column under the selected conditions. The 34S/32S ratio in the Fab samples was determined with a long term repeatability of 5% (12 hours, N = 5).

Sulfur isotope determination of the inorganic compounds in the supernatant

The supernatant of the culture broth has been separated from all molecules with a size larger than 3 kDa by ultrafiltration (3 kDa Amicon®, Millipore). The sulfur in the remaining solution, mainly from inorganic sulfate, was analyzed by ICP-SFMS for its sulfur ratio (see above).

Biomass determination by dry cell mass

Two times 5 mL culture broth were centrifuged. The pellets were resuspended in 10 mL RO-H2O (reverse osmosis water) and recentrifuged. The washed pellets were again resuspended in RO-H2O, transferred to weighed beakers and dried at 105°C until constant weight.

SDS-PAGE and silver staining

10 μL of the IP concentrate were run on a non-reducing sodium dodecyl sulfate (SDS) 4 to 12% polyacrylamide gel (Invitrogen) with MOPS buffer (morpholinepropanesulfonic acid) at 200 V for one hour and silver stained according to the protocol in the Current Protocols in Molecular Biology, chapter 10.6 [38].

Fab3H6 quantification

Quantification was done by sandwich ELISA as described in a previous study [27].

Abbreviations

KDeg:

time constant for intracellular protein degradation

KDil:

time constant for protein dilution

KHar:

time constant for protein harvest

KSec:

time constant for protein secretion

Pe:

extracellular protein concentration

Pi:

intracellular protein concentration

P32S:

concentration of the 32S isotope containing protein

qDeg:

intracellular protein degradation rate

qDil:

protein dilution rate

qHar:

protein harvest rate

qPi:

intracellular protein formation rate

qPi32S:

intracellular protein formation rate of the 32S isotope containing protein

qSec:

protein secretion rate

t1/2 Deg:

half time of intracellular protein degradation

t1/2 Dil:

half time of protein dilution

t1/2 Sec:

half time of protein secretion.

References

  1. Porro D, Gasser B, Fossati T, Maurer M, Branduardi P, Sauer M, Mattanovich D: Production of recombinant proteins and metabolites in yeasts: when are these systems better than bacterial production systems?. Appl Microbiol Biotechnol. 2011, 89 (4): 939-948. 10.1007/s00253-010-3019-z.

    Article  CAS  Google Scholar 

  2. Idiris A, Tohda H, Kumagai H, Takegawa K: Engineering of protein secretion in yeast: strategies and impact on protein production. Appl Microbiol Biotechnol. 2010, 86 (2): 403-417. 10.1007/s00253-010-2447-0.

    Article  CAS  Google Scholar 

  3. Lopes TS, Hakkaart GJ, Koerts BL, Raué HA, Planta RJ: Mechanism of high-copy-number integration of pMIRY-type vectors into the ribosomal DNA of Saccharomyces cerevisiae. Gene. 1991, 105 (1): 83-90. 10.1016/0378-1119(91)90516-E.

    Article  CAS  Google Scholar 

  4. Marx H, Mecklenbräuker A, Gasser B, Sauer M, Mattanovich D: Directed gene copy number amplification in Pichia pastoris by vector integration into the ribosomal DNA locus. FEMS Yeast Res. 2009, 9 (8): 1260-1270. 10.1111/j.1567-1364.2009.00561.x.

    Article  CAS  Google Scholar 

  5. 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 Res. 2008, 36 (12): e76-10.1093/nar/gkn369.

    Article  Google Scholar 

  6. Gasser B, Maurer M, Rautio J, Sauer M, Bhattacharyya A, Saloheimo M, Penttilä M, Mattanovich D: Monitoring of transcriptional regulation in Pichia pastoris under protein production conditions. BMC Genomics. 2007, 8: 179-10.1186/1471-2164-8-179.

    Article  Google Scholar 

  7. Hohenblum H, Gasser B, Maurer M, Borth N, Mattanovich D: Effects of gene dosage, promoters, and substrates on unfolded protein stress of recombinant Pichia pastoris. Biotechnol Bioeng. 2004, 85 (4): 367-375. 10.1002/bit.10904.

    Article  CAS  Google Scholar 

  8. Kauffman KJ, Pridgen EM, Doyle FJ, Dhurjati PS, Robinson AS: Decreased protein expression and intermittent recoveries in BiP levels result from cellular stress during heterologous protein expression in Saccharomyces cerevisiae. Biotechnol Prog. 2002, 18 (5): 942-950. 10.1021/bp025518g.

    Article  CAS  Google Scholar 

  9. Ng DT, Spear ED, Walter P: The unfolded protein response regulates multiple aspects of secretory and membrane protein biogenesis and endoplasmic reticulum quality control. J Cell Biol. 2000, 150 (1): 77-88. 10.1083/jcb.150.1.77.

    Article  CAS  Google Scholar 

  10. Travers KJ, Patil CK, Wodicka L, Lockhart DJ, Weissman JS, Walter P: Functional and genomic analyses reveal an essential coordination between the unfolded protein response and ER-associated degradation. Cell. 2000, 101 (3): 249-258. 10.1016/S0092-8674(00)80835-1.

    Article  CAS  Google Scholar 

  11. Harding HP, Zhang Y, Ron D: Protein translation and folding are coupled by an endoplasmic-reticulum-resident kinase. Nature. 1999, 397 (6716): 271-274. 10.1038/16729.

    Article  CAS  Google Scholar 

  12. Schröder M, Kaufman RJ: ER stress and the unfolded protein response. Mutat Res. 2005, 569 (1-2): 29-63.

    Article  Google Scholar 

  13. Hirsch C, Gauss R, Horn S, Neuber O, Sommer T: The ubiquitylation machinery of the endoplasmic reticulum. Nature. 2009, 458 (7237): 453-460. 10.1038/nature07962.

    Article  CAS  Google Scholar 

  14. Zhang W, Zhao HL, Xue C, Xiong XH, Yao XQ, Li XY, Chen HP, Liu ZM: Enhanced secretion of heterologous proteins in Pichia pastoris following overexpression of Saccharomyces cerevisiae chaperone proteins. Biotechnol Prog. 2006, 22 (4): 1090-1095. 10.1021/bp060019r.

    Article  CAS  Google Scholar 

  15. Damasceno LM, Anderson KA, Ritter G, Cregg JM, Old LJ, Batt CA: Cooverexpression of chaperones for enhanced secretion of a single-chain antibody fragment in Pichia pastoris. Appl Microbiol Biotechnol. 2007, 74 (2): 381-389. 10.1007/s00253-006-0652-7.

    Article  CAS  Google Scholar 

  16. Klabunde J, Kleebank S, Piontek M, Hollenberg CP, Hellwig S, Degelmann A: Increase of calnexin gene dosage boosts the secretion of heterologous proteins by Hansenula polymorpha. FEMS Yeast Res. 2007, 7 (7): 1168-1180. 10.1111/j.1567-1364.2007.00271.x.

    Article  CAS  Google Scholar 

  17. Ruohonen L, Toikkanen J, Tieaho V, Outola M, Soderlund H, Keranen S: Enhancement of protein secretion in Saccharomyces cerevisiae by overproduction of Sso protein, a late-acting component of the secretory machinery. Yeast. 1997, 13 (4): 337-351. 10.1002/(SICI)1097-0061(19970330)13:4<337::AID-YEA98>3.0.CO;2-K.

    Article  CAS  Google Scholar 

  18. Gasser B, Sauer M, Maurer M, Stadlmayr G, Mattanovich D: Transcriptomics-based identification of novel factors enhancing heterologous protein secretion in yeasts. Appl Environ Microbiol. 2007, 73 (20): 6499-6507. 10.1128/AEM.01196-07.

    Article  CAS  Google Scholar 

  19. Noe D, Delenick J: Quantitative analysis of membrane and secretory protein processing and intracellular transport. J Cell Sci. 1989, 92 (Pt 3): 449-459.

    CAS  Google Scholar 

  20. Batt BC, Kompala DS: A structured kinetic modeling framework for the dynamics of hybridoma growth and monoclonal antibody production in continuous suspension cultures. Biotechnol Bioeng. 1989, 34 (4): 515-531. 10.1002/bit.260340412.

    Article  CAS  Google Scholar 

  21. Hirschberg K, Lippincott-Schwartz J: Secretory pathway kinetics and in vivo analysis of protein traffic from the Golgi complex to the cell surface. FASEB J. 1999, 13 (Suppl 2): S251-256.

    CAS  Google Scholar 

  22. Bibila T, Flickinger M: A model of interorganelle monoclonal antibody transport and secretion in mouse hybridoma cells. Biotechnol Bioeng. 1991, 38 (7): 767-780. 10.1002/bit.260380711.

    Article  CAS  Google Scholar 

  23. Whiteley E, Hsu T, Betenbaugh M: Modeling assembly, aggregation, and chaperoning of immunoglobulin G production in insect cells. Biotechnol Bioeng. 1997, 56 (1): 106-116. 10.1002/(SICI)1097-0290(19971005)56:1<106::AID-BIT12>3.0.CO;2-I.

    Article  CAS  Google Scholar 

  24. Gonzalez R, Andrews BA, Asenjo JA: Metabolic control analysis of monoclonal antibody synthesis. Biotechnol Prog. 2001, 17 (2): 217-226. 10.1021/bp000165b.

    Article  CAS  Google Scholar 

  25. Ho Y, Varley J, Mantalaris A: Development and analysis of a mathematical model for antibody-producing GS-NS0 cells under normal and hyperosmotic culture conditions. Biotechnol Prog. 2006, 22 (6): 1560-1569.

    Article  CAS  Google Scholar 

  26. Belle A, Tanay A, Bitincka L, Shamir R, O'Shea EK: Quantification of protein half-lives in the budding yeast proteome. Proc Natl Acad Sci USA. 2006, 103 (35): 13004-13009. 10.1073/pnas.0605420103.

    Article  CAS  Google Scholar 

  27. Baumann K, Maurer M, Dragosits M, Cos O, Ferrer P, Mattanovich D: Hypoxic fed-batch cultivation of Pichia pastoris increases specific and volumetric productivity of recombinant proteins. Biotechnol Bioeng. 2008, 100 (1): 177-183. 10.1002/bit.21763.

    Article  CAS  Google Scholar 

  28. Maurer M, Kühleitner M, Gasser B, Mattanovich D: Versatile modeling and optimization of fed batch processes for the production of secreted heterologous proteins with Pichia pastoris. Microb Cell Fact. 2006, 5: 37-10.1186/1475-2859-5-37.

    Article  Google Scholar 

  29. Losev E, Reinke C, Jellen J, Strongin D, Bevis B, Glick B: Golgi maturation visualized in living yeast. Nature. 2006, 441 (7096): 1002-1006. 10.1038/nature04717.

    Article  CAS  Google Scholar 

  30. Yee JC, Jacob NM, Jayapal KP, Kok YJ, Philp R, Griffin TJ, Hu WS: Global assessment of protein turnover in recombinant antibody producing myeloma cells. J Biotechnol. 2010, 148 (4): 182-193. 10.1016/j.jbiotec.2010.06.005.

    Article  CAS  Google Scholar 

  31. Pignède G, Wang HJ, Fudalej F, Seman M, Gaillardin C, Nicaud JM: Autocloning and amplification of LIP2 in Yarrowia lipolytica. Appl Environ Microbiol. 2000, 66 (8): 3283-3289. 10.1128/AEM.66.8.3283-3289.2000.

    Article  Google Scholar 

  32. Bibila TA, Flickinger MC: Use of a structured kinetic model of antibody synthesis and secretion for optimization of antibody production systems: I. Steady-state analysis. Biotechnol Bioeng. 1992, 39 (3): 251-261. 10.1002/bit.260390302.

    Article  CAS  Google Scholar 

  33. Shampine LF, Reichelt MW, Kierzenka JA: Solving index-I DAEs in MATLAB and Simulink. Siam Review. 1999, 41 (3): 538-552. 10.1137/S003614459933425X.

    Article  Google Scholar 

  34. Dragosits M, Stadlmann J, Albiol J, Baumann K, Maurer M, Gasser B, Sauer M, Altmann F, Ferrer P, Mattanovich D: The effect of temperature on the proteome of recombinant Pichia pastoris. J Proteome Res. 2009, 8 (3): 1380-1392. 10.1021/pr8007623.

    Article  CAS  Google Scholar 

  35. Gach J, Quendler H, Strobach S, Katinger H, Kunert R: Structural analysis and in vivo administration of an anti-idiotypic antibody against mAb 2F5. Mol Immunol. 2008, 45 (4): 1027-1034. 10.1016/j.molimm.2007.07.030.

    Article  CAS  Google Scholar 

  36. Kunert R, Weik R, Siegler G, Katinger H: Anti-idiotypic antibody inducing hiv-1 neutralizing antibodies. US Patent. 2005, 2005080240

    Google Scholar 

  37. Bonifacino J, Dell'Angelica E, Springer T: Immunoprecipitation. Curr Protoc Mol Biol. 2001, 10 (10): 16.

    Google Scholar 

  38. Sasse J, Gallagher SR: Staining proteins in gels. Curr Protoc Mol Biol. 2009, 10 (10): 16.

    Google Scholar 

Download references

Acknowledgements

We want to thank Katrin Benakovitsch for support in method development, especially with respect to immunoprecipitation, and Marizela Delic and Markus Buchetics for fruitful discussions. This project was financially supported by Boehringer Ingelheim RCV (Vienna, Austria). ABG was supported by the project "Competence team Molecular Biology" by the City of Vienna, MA27.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Michael Maurer.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MP performed most of the experimental work, formulated and solved the kinetic model, and drafted the manuscript. MM contributed to study design, data analysis and interpretation. GK and SH designed and performed ICP-MS analysis of sulfur isotopes. ABG implemented the kinetic model in MATLAB and set up a MATLAB script for sensitivity analysis. DM conceived of the study and contributed to data analysis, interpretation and manuscript writing. All authors read and approved the final manuscript.

Electronic supplementary material

12934_2011_542_MOESM1_ESM.M

Additional file 1: MATLAB implementation of the kinetic model. Time courses of the intracellular and extracellular labeled Fab concentrations after start of labeling are modeled dependent of the kinetic parameters of protein formation, degradation and secretion. Continuous cultivation (constant YDM and dilution rate) are prerequisits for numerically solving the differential equations. (M 1 KB)

12934_2011_542_MOESM2_ESM.M

Additional file 2: Sensitivity analysis of the kinetic parameters. At a time two of the parameters, namely degradation, secretion and intracellular protein formation are varied against each other to determine the effect on extracellular Fab3H6 concentrations. Three dimensional plots are automatically created using the MATLAB functions meshgrid and surf. (M 5 KB)

Authors’ original submitted files for images

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 ( https://creativecommons.org/licenses/by/2.0 ), 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

Pfeffer, M., Maurer, M., Köllensperger, G. et al. Modeling and measuring intracellular fluxes of secreted recombinant protein in Pichia pastoris with a novel 34S labeling procedure. Microb Cell Fact 10, 47 (2011). https://doi.org/10.1186/1475-2859-10-47

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1475-2859-10-47

Keywords