Adaptation of the dynamic method for measuring the specific respiration
Adaptation of the dynamic method for measuring the specific respiration
com
Technical Report
Marilena Martins Pamboukian1, Carlos Augusto Pereira2, Elisabeth de Fatima Pires Augusto3 and Aldo Tonso1
1 Departamento de Engenharia Química, Escola Politécnica, USP, Brazil
2 Laboratório de Imunologia Viral, Instituto Butantan, Brazil
3 Laboratório de Biotecnologia Industrial, Instituto de Pesquisas Tecnológicas – IPT, Brazil
Monitoring the specific respiration rate (QO2) is a valuable tool to evaluate cell growth and physi- Received 30 July 2010
ology. However, for low QO2 values the accuracy may depend on the measurement methodology, Revised 23 December 2010
as it is the case in animal cell culture. The widely used “Dynamic Method” imposes serious diffi- Accepted 26 April 2011
culties concerning oxygen transfer cancellation, especially through membrane oxygenation. This
paper presents an improved procedure to this method, through an automated control of the gas
inlet composition that can minimize the residual oxygen transfer driving force during the QO2
measurement phase. The improved technique was applied to animal cell cultivation, particularly
three recombinant S2 (Drosophila melanogaster) insect cell lines grown in a membrane aeration
bioreactor. The average measurements of the proposed method reached 98% of stationary liquid
phase balance method, taken as a reference, compared to 21% when the traditional method was
used. Furthermore, this methodology does not require knowledge of the volumetric transfer coef-
ficient kLa, which may vary during growth.
Keywords: Animal cell culture · Biochemical engineering · Diffusion membrane · Oxygen transfer · Specific respiration rate
According to Henry’s law, the concentration of The application of the dynamic method to mem-
DO relates to the oxygen fraction in the gas as: brane aeration bioreactors presents an additional
pO2 P ⋅ yO2 complication compared to gas sparging system: be-
cS = = (2) sides surface transfer that is not cancelled, the
H H
residual gas in the membrane after halting the flow
where H (L.atm/gO2) is the Henry’s law, pO2 (atm) is creates a positive gradient between this gas and the
the partial pressure of oxygen in the gas phase, P DO in the liquid, which produces an extra transfer.
(atm) is the total pressure and yO2 is the mole frac- This latter transfer can be significant since the tub-
tion of oxygen in the gas phase. ing surface can be many times higher than head
QO2 depends on the estimation of OUR and cell space surface.
concentration, the latter obtained on-line through This paper proposes a modification to the “Dy-
specific sensors (turbidimeter), or off-line through namic Method“ to reduce this O2 transfer. The re-
different techniques applied to samples taken from sults obtained with this new proposed dynamic
the bioreactor (e.g., for cell count, gravimetric, tur- method are compared to the traditional procedure
bidity). Two methods have been widely exploited to and that proposed by Kamen et al. [6]. These tech-
determine OUR in bioprocesses [7]: the gas balance niques were applied to three different recombinant
and the dynamic method. The gas balance is based S2 (Drosophila melanogaster) insect cell strains
on the measurement of the difference between grown in a membrane aeration bioreactor.
oxygen mole flow rates in the bioreactor inlet and
outlet gases. The dynamic method requires a total
cancellation of the oxygen transfer to the bioreac- 2 Materials and methods
tor for a short time interval, during which the re-
duction of the DO concentration, due exclusively to 2.1 Cells and culture medium
cell respiration, is monitored [8]. OUR is deter-
mined as the slope of the curve C = f (t) (Eq. 1). Three experiments were performed using three
Gas balance method allows on-line OUR moni- genetically modified cell lines derived from D.
toring, but, on the other hand, has disadvantages melanogaster – S2 (Invitrogen, USA): Run 1 was
that may hinder its application. This is a technique carried out with a cell strain (S2AcRVGP) constitu-
that requires additional equipment; more impor- tively expressing the rabies virus glycoprotein
tant is that small differences between inlet and out- (RVGP) [10]; Run 2 was carried out with a cell lin-
let gas composition may result in large inaccuracies eage [S2MtEGFP] expressing the enhanced green
as in animal cell cultivation that presents low oxy- fluorescent protein (EGFP) after induction with
gen consumption rates. CuSO4 [11]; Run 3 was carried out with a cell strain
In turn, the dynamic method requires a com- (S2AcHBsAgHy) constitutively synthesizing hepa-
plete cancellation of oxygen, which is usually ac- titis B surface antigens (HBsAg) [12].
complished, in sparged transfer systems, by stop- Cells were cultivated in batch mode, using
ping gas inlet. In effect, this does not cancel super- serum-free medium SF900 II® (Invitrogen, USA),
ficial O2 transfer at the bioreactor headspace, and were inoculated with 5 × 105 cells/mL [13].
which is not significant for high O2 transfer condi- In Run 2, 700 µM CuSO4 was added after
tions, but may result in an important error in OUR 2.5 days to induce protein synthesis, and in Run 3,
measuring for low O2 transfer situations. Neverthe- 5 mM sodium butyrate was added after 5 days to
less, in most animal cell cultures, sparging is not enhance protein expression.
commonly used, since bubbles can easily damage Total cell density was determined with a hema-
cells due to their lack of wall. The degree of sus- cytometer and the viability by the trypan blue ex-
ceptibility is strain dependent and insect cells are clusion method.
among the most susceptible. To minimize these
harmful effects, it is common to use bubble-free 2.2 System description
aeration systems; the most successful alternative is
the diffusion of gases through membranes. These The tests were carried out in a 2-L bench bioreac-
membranes are made of materials with high gas tor Discovery 100 (Inceltech, France) with 1-L
permeability such as silicone or teflon, leading to working volume, with temperature controlled at
an acceptable transfer rate between the gas inside 28°C and agitation at 100 rpm.
the membrane and the liquid phase in the bioreac- The oxygen supply was made through silicone
tor. Several meters of tubing can be installed inside membrane with the following characteristics:
the reactor, resulting in a large area of gas ex- Silastic® silicone tubing (Dow Corning Corpora-
change [9]. tion, USA) with internal diameter of 2 mm, outer di-
Figure 1. Scheme of the oxygen transfer control system through a diffu- OUR = kLa · (CS – C) (3)
sion membrane. (1) Oxygen mass flow controller, (2) mixed gas mass
flow controller, (3) control system, (4) headspace valve, (5) DO probe,
(6) exhaust valve (back pressure regulator), (7) diffusion membrane. Considering kLa, C and H constants and known
throughout the cultivation, OUR can be calculated
in real time based on CS (Eq. 5). In our system,
ameter of 3.2 mm and 8.3 m length (interface area kLa was determined as described by Bartholomew
of 0.052 m2). et al. [14] in the absence of cells, in water, under the
To maintain DO at 40% air saturation, a polaro- same experimental conditions: geometry of the
graphic sensor InPro 6110 (Mettler Toledo, bioreactor, gas flow rate and mean membrane pres-
Switzerland) was used to inform a LabVIEW (Na- sure. An average value of 1.18 ± 0.06 h-1 was ob-
tional Instruments, USA) control routine of any de- tained.
viation of C. In order to re-establish the desired C The PDM methodology applies the same princi-
value in liquid, the routine determined a new CS ple of the dynamic method described above, but
(Eq. 1) and the corresponding yO2 (Eq. 2) that was uses a different strategy to cancel oxygen transfer.
implemented through two mass flow controllers Instead of zeroing gas flow, a gas mixture with a yO2
(MFC) 5850EM (Brooks Instruments, Holland) that is in equilibrium with the O2 concentration of
used to blend O2 and N2 gas lines (Fig. 1). Total gas the liquid is circulated through the membrane,
flow rate was kept constant at 52 mL/min. seeking CS = C, and then cancelling the transfer
In the proposed configuration, oxygen transfer driving force during the OUR measurement. For
varies with oxygen partial pressure (pO2) within the this, the control system sends continuously a signal
silicone tubing, which depends on the total pres- to the MFC to adjust O2 flow so that the gas mixture
sure inside the membrane (P) and the mass frac- has the desired composition, i.e., Cs = C. Also, to en-
tion of oxygen in the aeration gas (yO2). As the total sure cancellation of surface transfer, valve 4 (Fig. 1)
flow rate is kept constant, the total pressure is also is opened, allowing the same gas mixture flow si-
constant. The total pressure of the system is ap- multaneously to the headspace.
proximately equal to the mean logarithmic pres-
sure at the silicone tube inlet and outlet [9]. The
outlet pressure is regulated by a valve (10BP Back 3 Results and discussion
Pressure Regulator; Fairchild, USA), which can be
adjusted from 0.03 to 2.0 bar. Growth profiles before induction (Run 2) and pro-
tein expression enhancement (Run 3) were very
2.3 Measurements of specific oxygen consumption similar for the three runs. As an illustration, Fig. 2
rate presents results of Run 2.Table 1 shows data for the
three experiments.
QO2 was calculated based on OUR and cell concen- In Fig. 2A, it can be seen that the DO was pre-
tration in each sample (Eq. 1). Three methods for cisely controlled. The peaks observed relate to DO
measuring OUR were used: (i) the traditional dy- concentration changes, and result from the appli-
namic method (TDM), (ii) continuous determina- cation of TDM and PDM. As expected, QO2 (Fig. 2A)
tion by stationary liquid phase balance (SLPB), and and specific growth rate (data not shown) had
(iii) determination using the proposed dynamic very similar profiles, showing a plateau of maxi-
method (PDM), with cancellation of the oxygen mum values, up to the CuSO4 induction (3.25 days).
transfer driving force. After this, QO2 diminished drastically, although
growth was still significant, indicating a stronger The explanation for this difference between the
effect of induction on oxygen uptake. Values of TDM and PDM derives from the fact that the TDM
4.16 × 10–10 mg O2/cell.h for wild S2 cells [13], and is unable to completely cancel the O2 transfer
11.8 × 10–10 to 23.7 × 10–10 mg O2/cell.h for modified through diffusion membranes. The PDM cancels
cells [15] has been reported, and are very similar to the transfer using a flow of gas through the mem-
those shown in this paper. Those QO2 values are in- brane and headspace, the composition of which is
ferior to that described in literature for other insect in equilibrium with the DO in the liquid, zeroing
cells: 48.0 × 1010 mg O2/cell.h [4, 13, 15, 16] or for the driving force.
mammalian cell lines (CHO, HeLa and NS0): Figure 2B clearly shows the differences be-
68.0 × 10–10 to 131.2 × 10–10 mg O2/cell.h [16–18]. tween the DO concentration profiles in the two
The QO2 measurements obtained using TDM methods, and explains the differences in QO2 val-
and PDM present very different results, with PDM ues. In the TDM, the membrane retains the gas cir-
always showing higher values than those of TDM. culation just before procedure initialization and
Since it is impossible to apply TDM and PDM si- maintains the transfer of oxygen to the system dur-
multaneously, individual results were compared ing the test. In this case, the DO concentration de-
with those of a third method (SLPB), to evaluate caying due to cellular consumption is always small-
their accuracy. er than expected (Fig. 2B), since a residual oxygen
Table 1 presents the ratio between the PDM and transfer exists, contributing to its increase. The
SLPB, as well as between TDM and SLPB, calculat- variation observed results from the balance be-
ed for all available data, and their average, SD and tween consumption and residual transfer (Eq. 1).
CV(%) (100 × SD/average) values. The average of Depending on the composition of gas trapped in
PDM/SLPB ratio was 98%, whereas that of TDM/ the system, differences can be very significant.
SLPB was just 21%. These average values are sta- To illustrate the errors involved in the TDM, we
tistically different (t-test, for α = 0.05) and attest estimated the oxygen transfer through the mem-
that the proposed method gives results that are in brane at the beginning of the method, when gas
most cases much closer to SLPB. flow rate is interrupted.
5.95 15.60 1.54 1.80 1.17 25.1 5.20 10.68 1.18 0.14 0.12 21.3
6.77 22.35 1.32 1.40 1.06 26.4 6.29 19.60 1.24 0.11 0.09 27.6
7.23 24.86 1.10 1.15 1.05 26.4 7.79 26.60 1.01 0.06 0.06 26.4
8.25 24.72 1.12 1.21 1.08 23.9
Run 1 – S2AcGPV
9.2 19.30 1.40 1.54 1.10 25.1
2.2 3.36 2.79 2.97 1.06 11.3 4.13 8.46 1.99 0.11 0.06 17.6
2.87 4.90 2.74 2.30 0.84 13.8 5.89 17.80 0.61 0.05 0.08 11.3
Run 2 –
4.68 10.15 2.17 1.50 0.69 21.3
S2MtEGFP
5.34 13.80 0.98 0.37 0.38 10.0
7.03 4.10 3.26 3.78 1.16 11.3 4.95 1.20 6.05 2.00 0.33 8.8
8.87 10.40 1.98 2.13 1.08 16.3 6.03 2.60 4.44 1.23 0.28 11.3
10.89 16.60 1.14 1.56 1.37 20.1 7.97 6.60 1.84 0.78 0.42 13.8
Run 3 –
9.91 14.20 1.23 0.46 0.37 18.8
S2AcHBsAg
average 0.98 average 0.21
SD 0.24 SD 0.15
CV 24.6 CV 71.7
a) X: viable cell concentration, CS residual: CS value just before starting the dynamic method.
www.biotechnology-journal.com
1501
Biotechnology Biotechnol. J. 2011, 6, 1497–1503
Journal
OTRresidual = kLa . (CS residual – C) (4) area. In contrast, the addition of a surfactant such
as Pluronic reduced the kLa up to 35% [21–23].
Where: OTRresidual is the oxygen transfer rate; Hence changes in the composition of the culture
CS residual is the oxygen saturation concentration in medium must change kLa, making constant kLa a
equilibrium with the gas inside the membrane (see hypothesis difficult to sustain. Moreover, there is
Table 1); C is the oxygen concentration in the cul- evidence that the presence of cells as suspended
ture medium, controlled at 40% of saturation with solids interfere with the measurement of kLa [24].
air. Thus, for Run 2, at 3.25 days (Fig. 2B; Table 1), Therefore, the assumption of a constant kLa, re-
we obtain: quired for the SLPB methodology, may compromise
the measurements.
OTRresidual = 1.01 × (17.6 – 0.4 × 6.99) = In this sense, the PDM, a more laborious tech-
15.0 mg O2/L.h nique that requires continual gas formulation,
adapts more easily to the common variations that
where kLa = 1.01 h–1 calculated through this PDM occur throughout the culture as a function of me-
(considered more precise) and C = 40% of oxygen tabolism and cell growth.Thus, it constitutes an ex-
solubility in water at 28°C, according to U.S. Geo- cellent choice for measuring specific respiration
logical Survey Techniques of Water-Resources In- rate in oxygen transfer through membranes or un-
vestigations (available at http://pubs.water.usgs. der low transfer conditions, when surface transfer
gov/twri9A6). contribution may significantly interfere in the
OTR calculated from the TDM values (Fig. 2B; TDM.
Table 1; 4.13 days; no oxygen accumulation) gives: