# Correlation between mass transfer coefficient k_{L}a and relevant operating parameters in cylindrical disposable shaken bioreactors on a bench-to-pilot scale

- Wolf Klöckner
^{1}, - Riad Gacem
^{1}, - Tibor Anderlei
^{2}, - Nicole Raven
^{3}, - Stefan Schillberg
^{3}, - Clemens Lattermann
^{1}and - Jochen Büchs
^{1}Email author

**7**:28

https://doi.org/10.1186/1754-1611-7-28

© Klöckner et al.; licensee BioMed Central Ltd. 2013

**Received: **23 May 2013

**Accepted: **4 November 2013

**Published: **2 December 2013

## Abstract

### Background

Among disposable bioreactor systems, cylindrical orbitally shaken bioreactors show important advantages. They provide a well-defined hydrodynamic flow combined with excellent mixing and oxygen transfer for mammalian and plant cell cultivations. Since there is no known universal correlation between the volumetric mass transfer coefficient for oxygen k_{L}a and relevant operating parameters in such bioreactor systems, the aim of this current study is to experimentally determine a universal k_{L}a correlation.

### Results

A Respiration Activity Monitoring System (RAMOS) was used to measure k_{L}a values in cylindrical disposable shaken bioreactors and Buckingham’s π-Theorem was applied to define a dimensionless equation for k_{L}a. In this way, a scale- and volume-independent k_{L}a correlation was developed and validated in bioreactors with volumes from 2 L to 200 L. The final correlation was used to calculate cultivation parameters at different scales to allow a sufficient oxygen supply of tobacco BY-2 cell suspension cultures.

### Conclusion

The resulting equation can be universally applied to calculate the mass transfer coefficient for any of seven relevant cultivation parameters such as the reactor diameter, the shaking frequency, the filling volume, the viscosity, the oxygen diffusion coefficient, the gravitational acceleration or the shaking diameter within an accuracy range of +/− 30%. To our knowledge, this is the first k_{L}a correlation that has been defined and validated for the cited bioreactor system on a bench-to-pilot scale.

## Keywords

## Background

The success of new biopharmaceuticals highly depends on their potential to compete with established products. Parameters such as a fast time to market, cost effectiveness and manufacturing flexibility are key issues that need to be considered while maintaining product quality [1]. Using disposable equipment can help reduce investment costs and increase flexibility and process safety [2] ultimately leading to enhanced competitiveness of the products. Consequently, disposable bioreactors are increasingly being used in biotechnological production processes within the past ten years [3].

Today, one of the most popular disposable bioreactors is the WAVE cultivation system that was first described in 1976 [4, 5]. Since then, various other types of disposable cultivation systems for different applications have been developed. A wide range of different reactor sizes, starting with disposable screening systems in microliter scale, up to bioreactors with capacities of several cubic meters are on the market [3]. Available cultivation systems and their applications are described in several review articles [2, 3, 6, 7].

Despite the increasing usage of disposable cultivation systems, only few scientific publications have reported about their characterization with respect to power input, oxygen transfer and mixing performance [8]. By contrast, stirred stainless steel bioreactors have been extensively characterized and optimized over the past 60 years [3]. A direct design transfer from stainless steel to stirred disposable bioreactor systems is not feasible because of different properties of the applied materials. For instance, it is not possible to keep geometric similarity between a stirrer made of stainless steel and a disposable stirrer, due to the reduced strength and rigidity of the disposable polymer material. New agitation concepts for disposable bioreactors are desired to simplify their design and make them more cost-efficient [6]. Surface aerated reactors without complex built-in components fulfil the requirement for a cost-efficient reactor design [9]. In these systems, oxygen transfer and power input are either introduced by a wave, rocking or shaking motion of the bioreactor [3]. Orbitally shaken bioreactors are advantageous due to their well-defined liquid distribution that allows a precise characterization of power input and oxygen transfer. The magnitude of volumetric power input in these systems is comparable to that of conventional stirred tank reactors [10, 11], indicating excellent mixing characteristics also for liquids at elevated viscosity.

In the 50 mL scale, the orbitally shaken TubeSpin Bioreactor was developed and tested as a high throughput system for mammalian cell cultivation [12, 13]. On a larger scale, different studies reported about power input [10, 11, 14, 15], mixing properties [16] and scale-up performance [15, 17–20] of disposable orbitally shaken reactors. In several studies the volumetric mass transfer coefficient k_{L}a was identified as one of the main parameters for successful scale-up and transfer of cultivation conditions from stirred tank to disposable shaken reactors [13, 18, 21, 22]. Zhang et al. [23] also described a helical track attached on the inner reactor wall as a potential means to increase the mass transfer coefficient. However, a scale- and volume-independent k_{L}a correlation for cylindrical orbitally shaken bioreactors has not been described so far. Thus, the aim of this current study is to experimentally determine such a k_{L}a correlation in cylindrical disposable shaken bioreactors with volumes from 2 L up to 200 L and apply it to the cultivation of suspended plant cells.

## Results and discussion

### Dimensional analysis

_{L}a correlation [24, 25]. According to the rules of Buckingham’s theory, the resulting correlation is restricted to systems with geometric similarity. Influencing variables for k

_{L}a in cylindrical orbitally shaken bioreactors and their corresponding units are the volumetric mass transfer coefficient for oxygen ${\mathrm{k}}_{\mathrm{L}}\mathrm{a}\phantom{\rule{0.25em}{0ex}}\left[\frac{1}{\mathrm{s}}\right]$, the reactor diameter d [m], the shaking diameter d

_{0}[m], the shaking frequency $\mathrm{n}\phantom{\rule{0.25em}{0ex}}\left[\frac{1}{\mathrm{s}}\right]$, the liquid volume V

_{L}[m

^{3}], the diffusion coefficient for oxygen ${\mathbb{D}}_{\mathrm{O}2}\phantom{\rule{0.25em}{0ex}}\left[\frac{{\mathrm{m}}^{2}}{\mathrm{S}}\right]$, the kinematic viscosity $\mathrm{\nu}\phantom{\rule{0.25em}{0ex}}\left[\frac{{\mathrm{m}}^{2}}{\mathrm{s}}\right]$ and the gravitational acceleration $\mathrm{g}\phantom{\rule{0.25em}{0ex}}\left[\frac{\mathrm{m}}{{\mathrm{s}}^{2}}\right]$. The following dimensionless numbers were formed with the influencing variables according to the rules of Buckingham’s theory:Mass transfer number

This set of dimensionless numbers was also proposed by Henzler and Schedel [26] to develop a scale-independent correlation for the mass transfer coefficient in shake flasks. The applied numbers for surface aerated bioreactors differ from dimensionless numbers that are commonly used to describe the mass transfer in bubble aerated bioreactors because the principle of mass transfer differs fundamentally in both systems. In contrast to surface aerated orbitally shaken bioreactors, mass transfer in bubble aerated bioreactors depends on the amount, size, break-up and coalescence of gas bubbles. The characteristics of gas bubbles are strongly influenced by the volumetric power input, which is therefore commonly used as a parameter for k_{L}a correlations for bubble aerated bioreactors. These aspects are not relevant for the utilized surface aerated bioreactors that are operated without bubble aeration.

Each exponent in Eq. 1 has to be determined separately by empirically varying the corresponding dimensionless number.

### Influence of the critical circulation frequency

_{c}has to be exceeded to induce a rotating liquid motion in cylindrical orbitally shaken vessels [10, 29]. Hardly any movement was observed between liquid bulk and reactor wall below the critical circulation frequency, resulting in very low power input and oxygen transfer values. Consequently, only k

_{L}a values that were measured at shaking frequencies above the critical circulation frequency were considered for determining the k

_{L}a correlation. The following equation was used to calculate the critical circulation frequency [10]:

According to Eq. 2, the critical circulating frequency N_{c} is a function of the inner reactor diameter d, the filling volume V_{L} and the acceleration of gravity g.

### Out-of-phase operation in orbitally shaken bioreactors

The transition from in-phase to out-of-phase operation in orbitally shaken bioreactors is accompanied by a strong decrease in power input and oxygen transfer [30–32]. Liquids with high viscosities in systems with low shaking diameters are prone to out-of-phase operation as described in several studies for shake flasks [32, 33]. Thus, a minimum value for the ratio of shaking diameter (d_{0}) to reactor diameter (d), expressed by the Geometric number, is required to avoid an undesired out-of-phase operation in cylindrical orbitally shaken bioreactors.

### Defining the k_{L}a correlation

_{L}a in cylindrical bioreactors with nominal volumes from 2 L to 50 L were determined according to Eq. 11 by using a RAMOS device in combination with a 0.5 M or 1 M sulfite system for OTR

_{max}measurements. Values for k

_{L}a in cylindrical reactors of different size are presented in Figure 1. A decrease in k

_{L}a with increasing filling volume was observed. However, similar k

_{L}a values were obtained for the same relative filling volume in the applied scales (10 L - 50 L). This result is, on the first glance, astonishing as the volumetric transfer area (a) decreases with increasing vessel size. Similar k

_{L}a values in different scales for the same relative filling volume can be explained by an increasing value for the coefficient k

_{L}with increasing reactor size because of a higher power input [10]. A precise characterization of the mass transfer coefficient using dimensional analysis is presented in the following section.

_{L}a number was determined by varying the shaking frequency (n) as shown in Figure 2. A variation of the shaking frequency only affects the Froude number while all remaining dimensionless numbers are kept constant, allowing the influence of the Froude number to be separated from the influence of the other dimensionless numbers. An average exponent of α = 1.06 was determined as shown in Figure 2. The validity of the exponent α is restricted to a range for the Froude number of between 0.013 and 0.097 as depicted in Figure 2.

_{L}a number was investigated by altering the shaking diameter d

_{0}. A modification of the shaking diameter in Eq. 1 leads to a variation of the Geometric number and the Froude number. Therefore, the influence of the Froude number on the k

_{L}a number had to be considered with the exponent α = 1.06. This calculation was conducted by multiplying the k

_{L}a number with the inverse Froude number as stated in Figure 4. An average exponent of γ = -1.06 was determined within a range of variation for the Geometric number from 0.04 to 0.42 as shown in Figure 4. A low value for the ratio of shaking diameter (d

_{o}) to reactor diameter (d), which is expressed with the Geometric number, can lead to an undesired out-of-phase operation. Thus, the minimum value for the Geometric number was set to 0.06 to prevent an out-of-phase operation during the measurements in the present work. The final range of applicability for the exponent γ of the Geometric number was defined between 0.06 and 0.42.

_{L}a number as stated in Figure 5. An average exponent of δ = -0.12 was determined with a validity range for the Galilei number of between 6.2∙10

^{10}and 1.88∙10

^{12}.

_{2}SO

_{3}) solutions with concentrations of 0.5 mol/L and 1 mol/L and measurements with deionized water, using the dynamic gassing-out method, were applied to vary ${\mathbb{D}}_{\mathrm{O}2}$ as stated in Table 1. A variation of the Na

_{2}SO

_{3}concentration leads to a change in the liquid viscosity which results in a change in the Galilei number and Schmidt number in Eq. 1. Thus, the influence of the Galilei number was considered during the determination of the exponent ϵ of the Schmidt number (Figure 6). The average value of ϵ = -0.12 can only be considered as a rough estimation because of the high experimental deviation for ϵ between values of −0.227 and 0.022 and the restricted range of variation for the Schmidt number between 417 and 878. However, the exponents of the Galilei number and Schmidt number are notably smaller than the exponents of the remaining numbers which implies that the influence of inaccuracies in these exponents on the k

_{L}a number is also notably smaller.

**Oxygen solubility, diffusion coefficient and kinematic viscosity for the applied solutions at 25°C**

Composition of the solution | Solubility (L | Diffusion coefficient for oxygen (${\mathbb{D}}_{\mathbf{O}\mathbf{2}}$) [m | Kinematic viscosity (ν) [m |
---|---|---|---|

Deionized water | 1.227 | 2.14 · 10 | 0.89 · 10 |

0.5 mol/L Na | 0.844 | 1.897 · 10 | 1.10 · 10 |

1 · 10 | |||

0.012 mol/L phosphate buffer (pH 8) | |||

1 mol/L Na | 0.561 | 1.688 · 10 | 1.48 · 10 |

1 · 10 | |||

0.012 mol/L phosphate buffer (pH 8) |

_{L}a correlation can be written as:

^{-3}in Eq. 3 was fitted by comparing all measured and calculated values using the method of least squares. Equation 3 can be simplified to the following final k

_{L}a correlation:

_{0}in Eq. 4. Equations 3 and 4 are only applicable for shaking frequencies that are higher than the critical frequency N

_{C}according to Eq. 2 and are restricted to the following variation ranges for the Froude number: 0.013 < $\frac{{\mathrm{n}}^{2}\cdot {\mathrm{d}}_{0}}{\mathrm{g}}$ < 0.097, the Volume number: 0.18 < $\frac{{\mathrm{V}}_{\mathrm{L}}}{{\mathrm{d}}^{3}}$ < 0.65, the Geometric number: 0.06 < $\frac{{\mathrm{d}}_{0}}{\mathrm{d}}$ < 0.42, the Galilei number: 6.2∙10

^{10}< $\frac{{\mathrm{d}}^{3}\cdot \mathrm{g}}{{\mathrm{\nu}}^{2}}$ < 1.88∙10

^{12}and the Schmidt number: 417 < $\frac{\mathrm{\nu}}{{\mathbb{D}}_{\mathrm{O}2}}$ < 878. Exponents of the reactor diameter d and the shaking frequency n in Eq. 4 are more than double those of respective exponents of correlations for shake flasks as summarized by Klöckner and Büchs [34]. For instance, Henzler and Schedel [26] proposed the following k

_{L}a correlation for shake flasks:

An exponent of 2.03 for the influence of the reactor diameter d and 1 for the influence of the shaking frequency n were specified in Eq. 5 for shake flasks. Different characteristics with respect to the oxygen transfer are based on differences in the shape of the reactor systems. The conical shape of the shake flask wall prevents a strong expansion of the liquid surface with increasing shaking frequency and vessel size. This is not the case in cylindrical orbitally shaken bioreactors, resulting in higher exponents for the corresponding variables n and d.

_{L}a values in this study were compared with calculated values using Eq. 4. A comparison between measured and calculated k

_{L}a values in scales from 2 L to 200 L is presented in Figure 7. Filled symbols indicate measuring values where the shaking frequency was higher than the critical frequency N

_{C}according to Eq. 2. A range for the accuracy of +/− 30% was determined for bioreactors with volumes from 2 L to 50 L as depicted in Figure 7. Open symbols that fall within the +/− 30% accuracy range are still in-phase but close to out-of-phase operation. It was necessary to exclude these data points to define a precise validity range for the correlation. Three measurement values of the 200 L reactor system were slightly above the 30% range (Figure 7). This deviation might be caused by the oxygen-transfer enhancing influence of welded plastic seams in the 200 L bag system. This effect was not considered in Eq. 4, because most of the measurements for determining the exponents in Eq. 4 were conducted with disposable systems with volumes ranging from 2 L to 50 L that provide a smooth and homogeneous surface of the inner reactor wall.

The application of dimensionless numbers requires geometrical similarity. The utilized reactor systems for determining the exponents in Eq. 4 had a cylindrical reactor wall with flat bottom geometry. Thus, Eq. 4 is only applicable for such cylindrical bioreactors. The developed k_{L}a correlation can be used to calculate suitable shaking parameters for a sufficient oxygen supply by calculating the maximum oxygen transfer capacity (OTR_{max}) according to Eq. 10. To avoid an oxygen limitation during scale-up it is essential to select shaking parameters where the OTR_{max} is higher than the oxygen uptake of the culture.

_{L}a correlation was used to determine suitable cultivation conditions for a sufficient oxygen supply of

*Nicotiana tabacum*cv. Bright Yellow-2 (BY-2) suspension cells. A maximum OTR value of 6–8 mmol/(L∙h) was measured during the cultivation of tobacco BY-2 cells in 250 mL shake flasks using Murashige and Skoog medium supplemented with 30 g/L sucrose [35]. Thus, a maximum oxygen transfer capacity of OTR

_{max}= 8 mmol/(L∙h) is required to ensure a sufficient oxygen supply during cultivation of BY-2 cells in orbitally shaken bioreactors. As specified with Eq. 10, the OTR

_{max}is a function of the k

_{L}a, the dissolved oxygen concentration at equilibrium conditions (${\mathrm{L}}_{{\mathrm{O}}_{2}}$) and the oxygen partial pressure in the gas phase (${\mathrm{p}}_{{\mathrm{O}}_{2}}$). The dissolved oxygen concentration at equilibrium conditions in Murashige and Skoog medium supplemented with 30 g/L sucrose at ambient air pressure and 25°C was specified by Curtis [36] as: ${\mathrm{L}}_{{\mathrm{O}}_{2}}\cdot {\mathrm{p}}_{{\mathrm{O}}_{2}}$ = 7.92 mg/L : 32 g/mol = 0.248 mmol/L. For the minimum required k

_{L}a value it follows:

_{L}a value of 32 1/h is required to avoid an oxygen limitation of tobacco BY-2 cells cultivated in Murashige and Skoog medium containing 30 g/L sucrose. The required shaking parameters were determined according to Eq. (4) with a measured kinematic viscosity of ν = 1.18 ∙ 10

^{−6}m

^{2}/s, a diffusion coefficient for oxygen of D

_{O2}= 1.72 ∙ 10

^{−9}m

^{2}/s estimated according to Jamnongwong et al. [37] for Murashige and Skoog medium and the gravitational acceleration of g = 9.81 m/s

^{2}. The required shaking frequency for a 10 L reactor system (d = 0.24 m) with a filling volume of V

_{L}= 3 L results from Eq. (4):

A similar trend of the OTR signal was achieved at a shaking frequency of n = 160 rpm in the 10 L reactor system compared to the 250 mL shake flask. A plateau of the OTR signal, as it usually appears during cultivations with oxygen limitation [38], was not observed.

_{L}= 5 L is determined analogously to Eq. (7) to: n = 146 rpm. Thus, a shaking frequency of n = 160 rpm leads to a sufficient oxygen transfer during the cultivation of BY-2 cells. Figure 9 shows a comparison of the oxygen transfer and growth parameters during the cultivation of BY-2 cells in a 20 L bioreactor system compared to a 250 mL shake flask culture. The dissolved oxygen tension (DOT) was additionally monitored in the 20 L bioreactor system (Figure 9). The DOT signal was always above 20% during the cultivation, indicating a sufficient oxygen supply. As depicted in Figure 9, a comparable trend of growth parameters between the 20 L bioreactor system and 250 mL shake flasks was observed. An increase of the filling volume in the 20 L reactor system from 5 L to 10 L at otherwise equal conditions leads to a reduced k

_{L}a value according to Eq. (7) of k

_{L}a =17 1/h. This value is significantly lower than the required k

_{L}a value of 32 1/h for the cultivation of BY-2 cells. Consequently, a filling volume of 10 L leads to an oxygen limitation at the specified conditions as shown in Figure 10. The oxygen limitation is indicated by the plateau of the OTR signal after a cultivation time of 88 h. At the same time the dissolved oxygen tension was close to zero as an additional proof for an oxygen limitation. The results demonstrated that the developed k

_{L}a correlation can be used to determine suitable cultivation conditions for a sufficient oxygen transfer at different scales.

## Conclusions

A universally applicable equation was defined for calculating the mass transfer coefficient k_{L}a in disposable cylindrical bioreactors with volumes ranging from 2 L to 200 L. Important parameters as the critical circulation frequency for the induction of liquid motion as well as in-phase operation conditions were considered during the experiments and the data evaluation. It was demonstrated that the final k_{L}a equation can be applied to determine the OTR_{max} of cylindrical bioreactors and, in this way, enables the identification of suitable cultivation conditions such as for plant cell suspension cultures. As a result, the derived k_{L}a equation is an essential tool for the correct application of cylindrical disposable shaken bioreactors on a bench- to pilot-scale.

## Materials and methods

### Cultivation and agitation systems

**Properties of the utilized disposable reactor systems**

Product name | Nominal volume | Reactor diameter (d) | Material | Shaking platform | N |
---|---|---|---|---|---|

Corning roller bottle, 850 cm | 2 L | 116.3 mm | PS | Kühner ISF1-X | 147 rpm |

Nalgene 10 L Clearboy | 10 L | 240 mm | PC | Kühner SR 200-X Pilot-Shaker | 77 rpm |

Nalgene 20 L Carboy | 20 L | 286 mm | PP | Kühner SR 200-X Pilot-Shaker | 76 rpm |

Nalgene 50 L Carboy | 50 L | 379 mm | PP | Kühner SR 200-X Pilot-Shaker | 69 rpm |

Sartorius 200 L CultiBag ORB | 200 L | 750 mm | LDPE | Kühner SB 200-X OrbShake | 35 rpm |

### Model for the gas/liquid oxygen transfer

_{L}a) and the concentration difference between oxygen at the saturated gas/liquid interface $\left({\mathrm{C}}_{\mathrm{L}}^{*}\right)$ and oxygen in the bulk phase (C

_{L}) as:

_{max}) is reached at a dissolved oxygen concentration in the liquid bulk phase close to zero (C

_{L}≈ 0 mol/L) [28, 41]. For the OTR

_{max}follows from Eq. 9:

Values for the mass transfer coefficient k_{L}a were calculated from the measured OTR_{max} signal according to Eq. 11.

### The sulfite reaction system

_{2}SO

_{3}) reaction can be used to reduce the dissolved oxygen in the bulk phase (C

_{L}), and, in this way, to simulate oxygen consumption by a biological culture [28, 42–44]. The main advantages of a chemical reaction compared to a biological culture are the constant and reproducible oxygen consumption without sterility problems. The cobalt catalyzed sulfite oxidation is described by the following stoichiometric equation [44]:

Different reaction rates of the sodium sulfite oxidation reaction can be adjusted by varying the cobalt concentration. A non-accelerated reaction rate with a Hatta number (Ha) of Ha < 0.3 is required for k_{L}a measurements [44]. In this range, the sulfite reaction is able to reduce the dissolved oxygen concentration (C_{L}) to values close to zero. An increase in the oxygen transfer rate leads to a slight increase in C_{L} as described by Maier et al. [42]. Thus, C_{L} needs to be considered during measurements with high oxygen transfer rates, as they are usually reached in bubble aerated stirred tank reactors [42]. As comparatively low oxygen transfer rates of less than 16 mmol/(L∙h) were achieved with the disposable reactor systems in the present work, the assumption of (C_{L} ≈ 0 mol/L) is, in this case, applicable for the determination of k_{L}a values.

Sodium sulfite (Roth, Karlsruhe, Germany, purity < 98%) dissolved in deionized water and catalyzed with cobalt sulfate (Fluka, Neu-Ulm, Germany) was used for the oxidation reaction in the liquid phase. Two different sulfite concentrations of 0.5 mol/L and 1 mol/L were used to vary the diffusion coefficient for oxygen $\left({\mathbb{D}}_{02}\right)$ in the liquid phase. Both solutions were buffered with 12 mmol Na_{2}HPO_{4}/NaH_{2}PO_{4} buffer, and a pH value of 8 was adjusted with 30% (w/w) sulfuric acid. The oxygen solubility of the solutions were calculated according to Weisenberger and Schumpe [39] and oxygen diffusion coefficients according to Akita [45]. A ratio of ${\mathbb{D}}_{\mathrm{O}2,\mathrm{Sulfite}}/{\mathbb{D}}_{\mathrm{O}2,\mathrm{Water}}$ = 0.886 was determined for the 0.5 mol/L sulfite solution and a ratio of ${\mathbb{D}}_{\mathrm{O}2,\mathrm{Sulfite}}/{\mathbb{D}}_{\mathrm{O}2,\mathrm{Water}}$ = 0.788 was found for the 1 mol/L sulfite solution in agreement with values calculated by Linek and Vacek [46] using the same model proposed by Akita [45]. The viscosity was measured with an Anton Paar MCR 301 rheometer (Anton Paar GmbH, Graz, Austria). Values for the oxygen solubility, diffusion coefficients and liquid viscosity are listed in Table 1.

### Adaption of the Respiration Activity Monitoring System (RAMOS)

_{max}in combination with the sulfite reaction system. Values in the 200 L bioreactor system were determined with the dynamic gassing out method. A detailed description of the RAMOS device for shake flasks is given by Anderlei and Büchs [38] and Anderlei et al. [47]. A scheme of the modified RAMOS for cylindrical reactors is shown in Figure 11. An electrochemical oxygen partial pressure sensor (MAX-250B, Maxtec, Salt Lake City, UT, USA) and a total pressure sensor (No. 26PCAFA6D, Honeywell Sensing and Control, Golden Valley, MN, USA) were integrated in the headspace of the reactor systems. Magnetic flipper valves (Type 0332 E, Bürkert, Ingelfingen, Germany) were used to switch between measuring and rinsing phase and a mass flow controller was used to control the aeration rate during the rinsing phase (see section about aeration). Otherwise, the setup of the RAMOS device and the OTR calculation were like those for shake flasks described by Anderlei and Büchs [38]. Measurement values for OTR

_{max}were used to determined k

_{L}a values according to Eq. 11.

### Aeration of the 2 L Corning roller bottle reactor

The Corning roller bottle reactor (Corning Inc., Acton, MA, USA) with vent cap is equipped with a 0.2 μm membrane that is integrated in the cap to allow a sufficient aeration during cultivation. A method developed by Anderlei et al. [48] for determining the mass transfer resistance of shake flask closures was applied to measure the diffusive mass transfer through the membrane of the cultivation system. A volume flow rate of 9 mL/min was determined with open ports at 25°C. This value was used for the active aeration with a Brooks 5850 TR mass flow controller (Brooks Instrument, Ede, The Netherlands) in the RAMOS as pictured in Figure 11. An oxygen sensor was placed on top of the reactor systems to measure the oxygen partial pressure in the gas phase. No further changes of the RAMOS for shake flasks were required to measure the OTR_{max} of the Corning roller bottle cultivation system.

### Aeration of cylindrical orbitally shaken bioreactors with volumes from 10 L to 200 L

The settings of the air flow rate for the active headspace aeration of the cylindrical shaken bioreactors was transferred from the RAMOS device for shake flasks and adjusted according to the nominal reactor volume. A volume flow rate of 10 mL/min is usually applied in the RAMOS device for 250 mL shake flasks to mimic the oxygen transfer through a conventional cotton plug of shake flasks with a narrow neck. Consequently, an air flow rate of 400 mL/min was used for a volume of 10 L, 800 mL/min for 20 L, 2 L/min for 50 L and 8 L/min for 200 L bioreactors, respectively. For reactor volumes from 2–50 L, the flow rate was adjusted with a Brooks 5850-E mass flow controller (Brooks Instrument) in the RAMOS as illustrated in Figure 11. A manual air flow meter was used for the dynamic gassing-out method in the 200 L reactor system. The air flow rates in different scales were high enough to keep the absolute headspace concentration of oxygen, measured with the oxygen partial pressure sensor, above 20% during all OTR measurements.

### Application of the dynamic gassing-out method for k_{L}a measurements with water

A defined and constant headspace volume is required for OTR_{max} measurements with RAMOS [38]. This could not be provided in the flexible 200 L bag reactor system. Thus, the dynamic gassing-out method, first described by Bandyopadhyay et al. [49], was used for k_{L}a measurements in the 200 L scale. Nitrogen was used to replace the dissolved oxygen in the liquid phase, and deionized water was used as medium. Oxygen-sensitive spots (type SP-PSt3-YAU-D5-YOP) from PreSens (PreSens GmbH, Regensburg, Germany) were applied to measure the dissolved oxygen tension (DOT) in the liquid phase. An electrochemical sensor (MAX-250B, Maxtec, Salt Lake City, UT, USA) was additionally used to measure the oxygen concentration in the headspace of the reactor. The reactor system was filled with deionized water according to the desired filling volume. Then, the headspace of the system was filled with nitrogen, and the reactor was shaken at 80 rpm until the DOT reached a value below 1%. Subsequently, the shaker was stopped, and the gas volume in the headspace was replaced with air until the relative oxygen concentration in the headspace reached a value above 98%. The shaker was then immediately started with the designated shaking frequency, and the DOT was recorded over time. During the required time to replace the nitrogen in the headspace with air, the diffusion of oxygen from the gas phase to the liquid phase could cause small oxygen concentration differences in the liquid phase. Therefore, measured DOT values during the first 60 s were not considered for the calculation of the mass transfer coefficient (k_{L}a) to ensure a sufficient mixing of the liquid phase. According to Tissot et al. [16], a mixing time of 60 s can be regarded as sufficient for a 200 L reactor system using liquids with water-like viscosity and shaking frequencies above 60 rpm. During all measurements the DOT signal after the 60 s mixing step was below 4%, indicating that only small amounts of oxygen were transferred to the liquid phase during the time needed to replace the gas in the headspace. Values for k_{L}a were calculated from the recorded DOT signal over time as described by Van Suijdam et al. [50] and recently summarized by Suresh et al. [51].

### Cultivation of *Nicotiana tabacum* cv. BY-2 plant cell suspension cultures

For evaluation of the k_{L}a correlation the transgenic *N. tabacum* cv. Bright Yellow-2 (BY-2) MTAD cell line producing the human antibody M12 was used. The generation of the transgenic BY-2 cell line is described by Raven et al. [52]. BY-2 plant cell suspension cultures were cultivated at 26°C in the dark using Murashige and Skoog media [53] with minimal organics (# M6899, Sigma Aldrich, Saint Louis, MO, USA) supplemented with 30 g/L sucrose, 0.2 g/L KH_{2}PO_{4}, 0.6 mg/L thiamine-HCl, 0.2 mg/L 2.4-Dichlorophenoxyacetic acid (2.4-D) and, where stated, additional 100 mmol/L KNO_{3}. The pH of the culture medium was adjusted to 5.8 with 1 mol/L KOH before autoclaving for 21 min at 121°C. Plant cell suspensions cultures were sub-cultivated weekly for cell maintenance in 250 mL Erlenmeyer flasks, filled with 50 ml cell suspension, sealed with a cotton plug and shaken at 180 rpm and 26°C in the dark. Inoculation was conducted by adding 5% (V/V) of a seven day old culture to fresh medium.

### Oxygen transfer rate measurements in shake flasks

The OTR signals in 250 mL Erlenmeyer flasks were measured with a Respiration Activity Monitoring System (RAMOS). A detailed description of the device and its applications is provided by Anderlei et al. 2004 [47] and Anderlei and Büchs 2001 [38]. The non-invasive measuring system allows online monitoring of the OTR, CTR and RQ without changing the culture conditions compared to conventional 250 mL shake flasks [47]. Conventional Erlenmeyer flasks were used in addition to the RAMOS flasks to take samples during the cultivations.

### Determination of fresh and dry cell weight

An electronic precision balance (SBC 31, Scaltec, Göttingen, Germany) was used to determine fresh and dry cell weight. Fresh cell weight was determined by vacuum filtration of 10 ml cell suspension for 3 min using Whatman filter paper grade 3 (# 1003–055, 55 mm diameter, Fisher Scientific GmbH, Schwerte, Germany). Prior to the filtration step the filter paper was weighted dry and after wetting with purified water. The fresh cell weight was determined from the difference in weight of the membrane with cells directly after filtration and the wet membrane without cells. The membrane with cells was dried at 105°C until the mass remained constant. The dry cell weight was determined from the difference in weight of the dried membrane with and without cells.

### Determination of osmolality, dissolved oxygen tension and pH

Osmolality was determined in the supernatant with a Gonotec Osmomat 030 (Gonotec GmbH, Berlin, Germany). The device was calibrated with a two point calibration prior to each measurement. The pH value was determined with a pH510 pH meter (Eutech, Fisher Scientific GmbH, Schwerte, Germany). The dissolved oxygen tension (DOT) was measured by using oxygen sensitive sensor spots (type SP-PSt3-YAU-D5-YOP) from PreSens (PreSens GmbH, Regensburg, Germany).

## Declarations

### Acknowledgments

The authors thank the company Kühner AG for their helpful support with equipment and materials. This work has been financed within the EU Project CoMoFarm (grant agreement no.: 227420).

## Authors’ Affiliations

## References

- Oosterhuis NMG, Van den Berg HJ:
**How multipurpose is a disposable bioreactor?***Biopharm International*2011,**24:**51-56.Google Scholar - Eibl R, Kaiser S, Lombriser R, Eibl D:
**Disposable bioreactors: the current state-of-the-art and recommended applications in biotechnology.***Appl Microbiol Biotechnol*2010,**86:**41-49. 10.1007/s00253-009-2422-9View ArticleGoogle Scholar - Eibl R, Löffelholz C, Eibl D:
**Single-use bioreactors—an overview.**In*Single-Use Technology in Biopharmaceutical Manufacture*. Edited by: Eibl R, Eibl D. Hoboken, New Jersey: John Wiley & Sons, Inc; 2010:33-51.View ArticleGoogle Scholar - Kybal J, Vlcek V:
**Simple device for stationary cultivation of microorganisms.***Biotechnol Bioeng*1976,**18:**1713-1718. 10.1002/bit.260181206View ArticleGoogle Scholar - Kybal J, Sikyta B:
**A device for cultivation of plant and animal cells.***Biotechnol Lett*1985,**7:**467-470. 10.1007/BF01199860View ArticleGoogle Scholar - Brod H, Vester A, Kauling J:
**Opportunities and limitations of disposable technologies in biopharmaceutical processes.***Chem Ing Tech*2012,**84:**633-645. 10.1002/cite.201100229View ArticleGoogle Scholar - Eibl R, Eibl D:
**Disposable bioreactors in cell culture-based upstream processing.***Bioprocess Int*2009,**7:**18-23.Google Scholar - Rao G, Moreira A, Brorson K:
**Disposable bioprocessing: the future has arrived.***Biotechnol Bioeng*2009,**102:**348-356. 10.1002/bit.22192View ArticleGoogle Scholar - Kato I, Tanaka H:
**Development of a novel box-shaped shake flask with efficient gas exchange capacity.***J Ferment Bioeng*1998,**85:**404-409. 10.1016/S0922-338X(98)80084-6View ArticleGoogle Scholar - Klöckner W, Tissot S, Wurm F, Büchs J:
**Power input correlation to characterize the hydrodynamics of cylindrical orbitally shaken bioreactors.***Biochem Eng J*2012,**65:**63-69.View ArticleGoogle Scholar - Kato Y, Peter CP, Akgün A, Büchs J:
**Power consumption and heat transfer resistance in large rotary shaking vessels.***Biochem Eng J*2004,**21:**83-91. 10.1016/j.bej.2004.04.011View ArticleGoogle Scholar - De Jesus MJ, Girard P, Bourgeois M, Baumgartner G, Jacko B, Amstutz H, Wurm FM:
**TubeSpin satellites: a fast track approach for process development with animal cells using shaking technology.***Biochem Eng J*2004,**17:**217-223. 10.1016/S1369-703X(03)00180-3View ArticleGoogle Scholar - Zhang X, Stettler M, De Sanctis D, Perrone M, Parolini N, Discacciati M, De Jesus M, Hacker D, Quarteroni A, Wurm F:
**Use of orbital shaken disposable bioreactors for mammalian cell cultures from the milliliter-scale to the 1,000-liter scale.***Adv Biochem Eng/Biotechnol*2010,**115:**33-53.View ArticleGoogle Scholar - Raval K, Kato Y, Büchs J:
**Comparison of torque method and temperature method for determination of power consumption in disposable shaken bioreactors.***Biochem Eng J*2007,**34:**224-227. 10.1016/j.bej.2006.12.017View ArticleGoogle Scholar - Raval KN, Liu C, Büchs J:
**Large-scale disposable shaking bioreactors: a promising choice.***Bioprocess Int*2006,**4:**46. 48,50Google Scholar - Tissot S, Farhat M, Hacker DL, Anderlei T, Kühner M, Comninellis C, Wurm F:
**Determination of a scale-up factor from mixing time studies in orbitally shaken bioreactors.***Biochem Eng J*2010,**52:**181-186. 10.1016/j.bej.2010.08.005View ArticleGoogle Scholar - Liu CM, Hong LN:
**Development of a shaking bioreactor system for animal cell cultures.***Biochem Eng J*2001,**7:**121-125. 10.1016/S1369-703X(00)00111-XView ArticleGoogle Scholar - Zhang X, Bürki C-A, Stettler M, De Sanctis D, Perrone M, Discacciati M, Parolini N, DeJesus M, Hacker DL, Quarteroni A, Wurm FM:
**Efficient oxygen transfer by surface aeration in shaken cylindrical containers for mammalian cell cultivation at volumetric scales up to 1000 L.***Biochem Eng J*2009,**45:**41-47. 10.1016/j.bej.2009.02.003View ArticleGoogle Scholar - Werner S, Eibl R, Lettenbauer C, Roll M, Eibl D, De Jesus M, Zhang X, Stettler M, Tissot S, Bürki C,
*et al*.:**Innovative, non-stirred bioreactors in scales from milliliters up to 1000 Liters for suspension cultures of cells using disposable bags and containers - a swiss contribution.***Chimia*2010,**64:**819-823. 10.2533/chimia.2010.819View ArticleGoogle Scholar - Stettler M, Zhang X, Hacker DL, De Jesus M, Wurm FM:
**Novel orbital shake bioreactors for transient production of CHO derived IgGs.***Biotechnol Progr*2007,**23:**1340-1346. 10.1021/bp070219iView ArticleGoogle Scholar - Tissot S, Michel PO, Hacker DL, Baldi L, De Jesus M, Wurm FM:
**k(L)a as a predictor for successful probe-independent mammalian cell bioprocesses in orbitally shaken bioreactors.***New Biotechnol*2012,**29:**387-394. 10.1016/j.nbt.2011.10.010View ArticleGoogle Scholar - Tissot S, Oberbek A, Reclari M, Dreyer M, Hacker DL, Baldi L, Farhat M, Wurm FM:
**Efficient and reproducible mammalian cell bioprocesses without probes and controllers?***New Biotechnol*2011,**28:**382-390. 10.1016/j.nbt.2011.02.004View ArticleGoogle Scholar - Zhang X, Stettler M, Reif O, Kocourek A, DeJesus M, Hacker DL, Wurm FM:
**Shaken helical track bioreactors: Providing oxygen to high-density cultures of mammalian cells at volumes up to 1000 L by surface aeration with air.***New Biotechnol*2008,**25:**68-75. 10.1016/j.nbt.2008.03.001View ArticleGoogle Scholar - Grimvall G:
**Characteristic quantities and dimensional analysis.***Sci Model Simul*2008,**15:**21-39. 10.1007/s10820-008-9102-4View ArticleGoogle Scholar - Buckingham E:
**On physically similar systems; illustrations of the use of dimensional equations.***Phys Rev*1914,**4:**345-376. 10.1103/PhysRev.4.345View ArticleGoogle Scholar - Henzler HJ, Schedel M:
**Suitability of the shaking flask for oxygen supply to microbial cultures.***Bioprocess Eng*1991,**7:**123-131. 10.1007/BF00369423View ArticleGoogle Scholar - Doig SD, Pickering SCR, Lye GJ, Baganz F:
**Modelling surface aeration rates in shaken microtitre plates using dimensionless groups.***Chem Eng Sci*2005,**60:**2741-2750. 10.1016/j.ces.2004.12.025View ArticleGoogle Scholar - Hermann R, Lehmann M, Büchs J:
**Characterization of gas–liquid mass transfer phenomena in microtiter plates.***Biotechnol Bioeng*2003,**81:**178-186. 10.1002/bit.10456View ArticleGoogle Scholar - Kato Y, Hiraoka S, Tada Y, Koh ST, Lee YS:
**Mixing time and power consumption for a liquid in a vertical cylindrical vessel, shaken in a horizontal circle.***Chem Eng Res Des*1996,**74:**451-455.Google Scholar - Lotter S, Büchs J:
**Utilization of specific power input measurements for optimization of culture conditions in shaking flasks.***Biochem Eng J*2004,**17:**195-203. 10.1016/S1369-703X(03)00178-5View ArticleGoogle Scholar - Peter CP, Lotter S, Maier U, Büchs J:
**Impact of out-of-phase conditions on screening results in shaking flask experiments.***Biochem Eng J*2004,**17:**205-215. 10.1016/S1369-703X(03)00179-7View ArticleGoogle Scholar - Büchs J, Lotter S, Milbradt C:
**Out-of-phase operating conditions, a hitherto unknown phenomenon in shaking bioreactors.***Biochem Eng J*2001,**7:**135-141. 10.1016/S1369-703X(00)00113-3View ArticleGoogle Scholar - Büchs J, Maier U, Milbradt C, Zoels B:
**Power consumption in shaking flasks on rotary shaking machines: II. Nondimensional description of specific power consumption and flow regimes in unbaffled flasks at elevated liquid viscosity.***Biotechnol Bioeng*2000,**68:**594-601. 10.1002/(SICI)1097-0290(20000620)68:6<594::AID-BIT2>3.0.CO;2-UView ArticleGoogle Scholar - Klöckner W, Büchs J:
**Advances in shaking technologies.***Trends Biotechnol*2012,**30:**307-314. 10.1016/j.tibtech.2012.03.001View ArticleGoogle Scholar - Ullisch DA, Müller CA, Maibaum S, Kirchhoff J, Schiermeyer A, Schillberg S, Roberts JL, Treffenfeldt W, Büchs J:
**Comprehensive characterization of two different Nicotiana tabacum cell lines leads to doubled GFP and HA protein production by media optimization.***J Biosci Bioeng*2012,**113:**242-248. 10.1016/j.jbiosc.2011.09.022View ArticleGoogle Scholar - Curtis WR:
**Application of bioreactor design principles to plant micropropagation.***Plant Cell Tiss Org*2005,**81:**255-264. 10.1007/s11240-004-6646-1View ArticleGoogle Scholar - Jamnongwong M, Loubiere K, Dietrich N, Hébrard G:
**Experimental study of oxygen diffusion coefficients in clean water containing salt, glucose or surfactant: Consequences on the liquid-side mass transfer coefficients.***Chem Eng J*2010,**165:**758-768. 10.1016/j.cej.2010.09.040View ArticleGoogle Scholar - Anderlei T, Büchs J:
**Device for sterile online measurement of the oxygen transfer rate in shaking flasks.***Biochem Eng J*2001,**7:**157-162. 10.1016/S1369-703X(00)00116-9View ArticleGoogle Scholar - Weisenberger S, Schumpe A:
**Estimation of gas solubilities in salt solutions at temperatures from 273 K to 363 K.***AlChE J*1996,**42:**298-300. 10.1002/aic.690420130View ArticleGoogle Scholar - Rischbieter E, Schumpe A, Wunder V:
**Gas solubilities in aqueous solutions of organic substances.***J Chem Eng Data*1996,**41:**809-812. 10.1021/je960039cView ArticleGoogle Scholar - Maier U, Losen M, Büchs J:
**Advances in understanding and modeling the gas–liquid mass transfer in shake flasks.***Biochem Eng J*2004,**17:**155-167. 10.1016/S1369-703X(03)00174-8View ArticleGoogle Scholar - Maier B, Dietrich C, Büchs J:
**Correct application of the sulphite oxidation methodology of measuring the volumetric mass transfer coefficient k(L)a under non-pressurized and pressurized conditions.***Food Bioprod Process*2001,**79:**107-113. 10.1205/096030801750286267View ArticleGoogle Scholar - Funke M, Diederichs S, Kensy F, Mueller C, Büchs J:
**The baffled microtiter plate: increased oxygen transfer and improved online monitoring in small scale fermentations.***Biotechnol Bioeng*2009,**103:**1118-1128. 10.1002/bit.22341View ArticleGoogle Scholar - Linek V, Vacek V:
**Chemical engineering use of catalyzed sulfite oxidation kinetics for the determination of mass transfer characteristics of gas–liquid contactors.***Chem Eng Sci*1981,**36:**1747-1768. 10.1016/0009-2509(81)80124-8View ArticleGoogle Scholar - Akita K:
**Diffusivities of gases in aqueous electrolyte solutions.***Ind Eng Chem Fund*1981,**20:**89-94. 10.1021/i100001a017View ArticleGoogle Scholar - Linek V, Vacek V:
**Comments on validity of measuring oxygen diffusion coefficients with polarographic oxygen electrodes.***Biotechnol Bioeng*1988,**31:**1010-1011. 10.1002/bit.260310914View ArticleGoogle Scholar - Anderlei T, Zang W, Papaspyrou M, Büchs J:
**Online respiration activity measurement (OTR, CTR, RQ) in shake flasks.***Biochem Eng J*2004,**17:**187-194. 10.1016/S1369-703X(03)00181-5View ArticleGoogle Scholar - Anderlei T, Mrotzek C, Bartsch S, Amoabediny G, Peter CP, Büchs J:
**New method to determine the mass transfer resistance of sterile closures for shaken Bioreactors.***Biotechnol Bioeng*2007,**98:**999-1007. 10.1002/bit.21490View ArticleGoogle Scholar - Bandyopadhyay B, Humphrey AE, Taguchi H:
**Dynamic measurement of volumetric oxygen transfer coefficient in fermentation systems.***Biotechnol Bioeng*1967,**9:**533-544. 10.1002/bit.260090408View ArticleGoogle Scholar - Van Suijdam JC, Kossen NWF, Joha AC:
**Model for oxygen transfer in a shake flask.***Biotechnol Bioeng*1978,**20:**1695-1709. 10.1002/bit.260201102View ArticleGoogle Scholar - Suresh S, Srivastava VC, Mishra IM:
**Techniques for oxygen transfer measurement in bioreactors: a review.***J Chem Technol Biotechnol*2009,**84:**1091-1103. 10.1002/jctb.2154View ArticleGoogle Scholar - Raven N, Schillberg S, Kirchhoff J, Brändli J, Imseng N, Eibl R:
**Growth of BY-2 suspension cells and plantibody production in single-use bioreactors.**In*Single-Use Technology in Biopharmaceutical Manufacture*. Edited by: Eibl R, Eibl D. Hoboken, NJ, USA: John Wiley & Sons; 2010:251-261.Google Scholar - Murashige T, Skoog F:
**A revised medium for rapid growth and bio assays with tobacco tissue cultures.***Physiol Plantarum*1962,**15:**473-497. 10.1111/j.1399-3054.1962.tb08052.xView ArticleGoogle Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.