- Open Access
Correlation between mass transfer coefficient kLa and relevant operating parameters in cylindrical disposable shaken bioreactors on a bench-to-pilot scale
Journal of Biological Engineering volume 7, Article number: 28 (2013)
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 kLa and relevant operating parameters in such bioreactor systems, the aim of this current study is to experimentally determine a universal kLa correlation.
A Respiration Activity Monitoring System (RAMOS) was used to measure kLa values in cylindrical disposable shaken bioreactors and Buckingham’s π-Theorem was applied to define a dimensionless equation for kLa. In this way, a scale- and volume-independent kLa 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.
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 kLa correlation that has been defined and validated for the cited bioreactor system on a bench-to-pilot scale.
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 . Using disposable equipment can help reduce investment costs and increase flexibility and process safety  ultimately leading to enhanced competitiveness of the products. Consequently, disposable bioreactors are increasingly being used in biotechnological production processes within the past ten years .
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 . 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 . By contrast, stirred stainless steel bioreactors have been extensively characterized and optimized over the past 60 years . 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 . Surface aerated reactors without complex built-in components fulfil the requirement for a cost-efficient reactor design . In these systems, oxygen transfer and power input are either introduced by a wave, rocking or shaking motion of the bioreactor . 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  and scale-up performance [15, 17–20] of disposable orbitally shaken reactors. In several studies the volumetric mass transfer coefficient kLa 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.  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 kLa 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 kLa 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 according to Buckingham’s π -Theorem was used to define a scale- and volume-independent kLa correlation [24, 25]. According to the rules of Buckingham’s theory, the resulting correlation is restricted to systems with geometric similarity. Influencing variables for kLa in cylindrical orbitally shaken bioreactors and their corresponding units are the volumetric mass transfer coefficient for oxygen , the reactor diameter d [m], the shaking diameter d0 [m], the shaking frequency , the liquid volume VL [m3], the diffusion coefficient for oxygen , the kinematic viscosity and the gravitational acceleration . 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  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 kLa correlations for bubble aerated bioreactors. These aspects are not relevant for the utilized surface aerated bioreactors that are operated without bubble aeration.
The numbers are independent from one another and the set of numbers can be transformed into a different set of dimensionless numbers, with no influence on the resulting exponents of the included variables, as long as the rules of Buckingham’s π -Theorem are preserved . In this set of dimensionless numbers the Froude number is used to express the ratio of centrifugal to gravitational forces during shaking. The relative filling volume is considered with the Volume number and the Geometric number is used to describe the ratio of shaking diameter to reactor diameter. The height and shape of the free liquid surface during shaking is influenced by the ratio of gravitational to viscous forces, expressed by the Galilei number. The ratio of momentum transfer to diffusive mass transfer in the gas/liquid surface is described by the Schmidt number. In bioreactors without bubble aeration, the influence of the surface tension on the mass transfer coefficient decreases with increasing reactor size as already described by Doig et al.  and Hermann et al.  for microtiter plates. Consequently, the influence of the surface tension was not considered in orbitally shaken vessels with volumes from 2 L to 200 L. The occurrence of gas bubbles was avoided by using surface aeration only and the smooth surface of the inner reactor wall led to a well-defined and homogeneous liquid distribution during shaking. The following power law function was used to describe the correlation among the dimensionless numbers:
Each exponent in Eq. 1 has to be determined separately by empirically varying the corresponding dimensionless number.
Influence of the critical circulation frequency
As described by several authors, a critical circulation frequency Nc 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 kLa values that were measured at shaking frequencies above the critical circulation frequency were considered for determining the kLa correlation. The following equation was used to calculate the critical circulation frequency :
According to Eq. 2, the critical circulating frequency Nc is a function of the inner reactor diameter d, the filling volume VL 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 (d0) 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 kLa correlation
Values for kLa 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 OTRmax measurements. Values for kLa in cylindrical reactors of different size are presented in Figure 1. A decrease in kLa with increasing filling volume was observed. However, similar kLa 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 kLa values in different scales for the same relative filling volume can be explained by an increasing value for the coefficient kL with increasing reactor size because of a higher power input . A precise characterization of the mass transfer coefficient using dimensional analysis is presented in the following section.
The influence of the Froude number on the kLa 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.
The influence of the Volume number in Eq. 1 was determined by varying the filling volume, which only affects the Volume number. Thus, the resulting average exponent can be directly specified from the measuring values presented in Figure 3. An average exponent of β = -1.20 was determined for the influence of the Volume number. The range of applicability for the exponent β is restricted to values of the Volume number of between 0.18 and 0.65.
The influence of the Geometric number on the kLa number was investigated by altering the shaking diameter d0. 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 kLa number had to be considered with the exponent α = 1.06. This calculation was conducted by multiplying the kLa 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 (do) 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.
Vessels of different sizes were used to vary the reactor diameter and, in this way, the Galilei number. The influence of the Volume number and the Geometric number were considered with their respective exponents and multiplied with the kLa 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∙1010 and 1.88∙1012.
The variation of the Schmidt number in Eq. 1 was realized by varying the diffusion coefficient for oxygen . Sodium sulfite (Na2SO3) 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 as stated in Table 1. A variation of the Na2SO3 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 kLa number is also notably smaller.
With the determined exponents, the complete kLa correlation can be written as:
The constant coefficient of 1.06∙10-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 kLa correlation:
The same exponents for the Froude number and the Geometric number in Eq. 3 led to an elimination of the shaking diameter d0 in Eq. 4. Equations 3 and 4 are only applicable for shaking frequencies that are higher than the critical frequency NC according to Eq. 2 and are restricted to the following variation ranges for the Froude number: 0.013 < < 0.097, the Volume number: 0.18 < < 0.65, the Geometric number: 0.06 < < 0.42, the Galilei number: 6.2∙1010 < < 1.88∙1012 and the Schmidt number: 417 < < 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 . For instance, Henzler and Schedel  proposed the following kLa 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.
Experimentally determined kLa values in this study were compared with calculated values using Eq. 4. A comparison between measured and calculated kLa 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 NC 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 kLa correlation can be used to calculate suitable shaking parameters for a sufficient oxygen supply by calculating the maximum oxygen transfer capacity (OTRmax) according to Eq. 10. To avoid an oxygen limitation during scale-up it is essential to select shaking parameters where the OTRmax is higher than the oxygen uptake of the culture.
The newly developed kLa 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 . Thus, a maximum oxygen transfer capacity of OTRmax = 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 OTRmax is a function of the kLa, the dissolved oxygen concentration at equilibrium conditions () and the oxygen partial pressure in the gas phase (). 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  as: = 7.92 mg/L : 32 g/mol = 0.248 mmol/L. For the minimum required kLa value it follows:
Consequently, during scale up from shake flasks to orbitally shaken bioreactors a minimum kLa 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 m2/s, a diffusion coefficient for oxygen of DO2 = 1.72 ∙ 10−9 m2/s estimated according to Jamnongwong et al.  for Murashige and Skoog medium and the gravitational acceleration of g = 9.81 m/s2. The required shaking frequency for a 10 L reactor system (d = 0.24 m) with a filling volume of VL = 3 L results from Eq. (4):
A minimum shaking frequency of n = 156 rpm is required according to Eq. (7) to avoid an oxygen limitation at the aforementioned conditions. Figure 8 shows a comparison of the OTR signals of BY-2 cells cultivated in a 10 L orbitally shaken bioreactor and a 250 mL RAMOS shake flask.
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 , was not observed.
The minimum required shaking frequency for the cultivation of BY-2 cells in a 20 L reactor system (d = 0.286) with a filling volume of VL = 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 kLa value according to Eq. (7) of kLa =17 1/h. This value is significantly lower than the required kLa 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 kLa correlation can be used to determine suitable cultivation conditions for a sufficient oxygen transfer at different scales.
A universally applicable equation was defined for calculating the mass transfer coefficient kLa 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 kLa equation can be applied to determine the OTRmax 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 kLa 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
Table 2 presents an overview of the different geometries and material properties of the utilized disposable reactor systems and the respective shaking platform used in this study.
Model for the gas/liquid oxygen transfer
The oxygen transfer rate (OTR) from the gas phase to the liquid phase is commonly described as a function of the volumetric mass transfer coefficient (kLa) and the concentration difference between oxygen at the saturated gas/liquid interface and oxygen in the bulk phase (CL) as:
The oxygen concentration is a function of the oxygen solubility and the oxygen partial pressure in the gas phase With Eq. 8, it follows:
The oxygen solubility depends on the concentrations of salts and organic compounds in the cultivation media and can be estimated for salt solutions according to Weisenberger and Schumpe  or for aqueous solutions of organic compounds according to Rischbieter et al. . The influence of media components on the oxygen solubility, as discussed by Curtis for plant cell cultivation media , depends on the specific media composition. The maximum oxygen transfer capacity (OTRmax) is reached at a dissolved oxygen concentration in the liquid bulk phase close to zero (CL ≈ 0 mol/L) [28, 41]. For the OTRmax follows from Eq. 9:
Equation 10 can be rewritten as:
Values for the mass transfer coefficient kLa were calculated from the measured OTRmax signal according to Eq. 11.
The sulfite reaction system
A sodium sulfite (Na2SO3) reaction can be used to reduce the dissolved oxygen in the bulk phase (CL), 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 :
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 kLa measurements . In this range, the sulfite reaction is able to reduce the dissolved oxygen concentration (CL) to values close to zero. An increase in the oxygen transfer rate leads to a slight increase in CL as described by Maier et al. . Thus, CL needs to be considered during measurements with high oxygen transfer rates, as they are usually reached in bubble aerated stirred tank reactors . 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 (CL ≈ 0 mol/L) is, in this case, applicable for the determination of kLa 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 in the liquid phase. Both solutions were buffered with 12 mmol Na2HPO4/NaH2PO4 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  and oxygen diffusion coefficients according to Akita . A ratio of = 0.886 was determined for the 0.5 mol/L sulfite solution and a ratio of = 0.788 was found for the 1 mol/L sulfite solution in agreement with values calculated by Linek and Vacek  using the same model proposed by Akita . 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)
The Respiration Activity Monitoring System (RAMOS), initially developed for shake flasks, was adapted to the different disposable bioreactor systems with volumes ranging from 2 L to 50 L to measure the OTRmax 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  and Anderlei et al. . 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 . Measurement values for OTRmax were used to determined kLa 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.  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 OTRmax 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 kLa measurements with water
A defined and constant headspace volume is required for OTRmax measurements with RAMOS . 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. , was used for kLa 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 (kLa) to ensure a sufficient mixing of the liquid phase. According to Tissot et al. , 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 kLa were calculated from the recorded DOT signal over time as described by Van Suijdam et al.  and recently summarized by Suresh et al. .
Cultivation of Nicotiana tabacum cv. BY-2 plant cell suspension cultures
For evaluation of the kLa 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. . BY-2 plant cell suspension cultures were cultivated at 26°C in the dark using Murashige and Skoog media  with minimal organics (# M6899, Sigma Aldrich, Saint Louis, MO, USA) supplemented with 30 g/L sucrose, 0.2 g/L KH2PO4, 0.6 mg/L thiamine-HCl, 0.2 mg/L 2.4-Dichlorophenoxyacetic acid (2.4-D) and, where stated, additional 100 mmol/L KNO3. 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  and Anderlei and Büchs 2001 . 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 . 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).
Oosterhuis NMG, Van den Berg HJ: How multipurpose is a disposable bioreactor? Biopharm International 2011, 24: 51-56.
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-9
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.
Kybal J, Vlcek V: Simple device for stationary cultivation of microorganisms. Biotechnol Bioeng 1976, 18: 1713-1718. 10.1002/bit.260181206
Kybal J, Sikyta B: A device for cultivation of plant and animal cells. Biotechnol Lett 1985, 7: 467-470. 10.1007/BF01199860
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.201100229
Eibl R, Eibl D: Disposable bioreactors in cell culture-based upstream processing. Bioprocess Int 2009, 7: 18-23.
Rao G, Moreira A, Brorson K: Disposable bioprocessing: the future has arrived. Biotechnol Bioeng 2009, 102: 348-356. 10.1002/bit.22192
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-6
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.
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.011
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-3
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.
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.017
Raval KN, Liu C, Büchs J: Large-scale disposable shaking bioreactors: a promising choice. Bioprocess Int 2006, 4: 46. 48,50
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.005
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-X
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.003
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.819
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/bp070219i
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.010
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.004
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.001
Grimvall G: Characteristic quantities and dimensional analysis. Sci Model Simul 2008, 15: 21-39. 10.1007/s10820-008-9102-4
Buckingham E: On physically similar systems; illustrations of the use of dimensional equations. Phys Rev 1914, 4: 345-376. 10.1103/PhysRev.4.345
Henzler HJ, Schedel M: Suitability of the shaking flask for oxygen supply to microbial cultures. Bioprocess Eng 1991, 7: 123-131. 10.1007/BF00369423
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.025
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.10456
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.
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-5
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-7
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-3
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-U
Klöckner W, Büchs J: Advances in shaking technologies. Trends Biotechnol 2012, 30: 307-314. 10.1016/j.tibtech.2012.03.001
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.022
Curtis WR: Application of bioreactor design principles to plant micropropagation. Plant Cell Tiss Org 2005, 81: 255-264. 10.1007/s11240-004-6646-1
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.040
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-9
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.690420130
Rischbieter E, Schumpe A, Wunder V: Gas solubilities in aqueous solutions of organic substances. J Chem Eng Data 1996, 41: 809-812. 10.1021/je960039c
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-8
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/096030801750286267
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.22341
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-8
Akita K: Diffusivities of gases in aqueous electrolyte solutions. Ind Eng Chem Fund 1981, 20: 89-94. 10.1021/i100001a017
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.260310914
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-5
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.21490
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.260090408
Van Suijdam JC, Kossen NWF, Joha AC: Model for oxygen transfer in a shake flask. Biotechnol Bioeng 1978, 20: 1695-1709. 10.1002/bit.260201102
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.2154
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.
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.x
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).
The authors declare that they have no competing interests.
WK conducted experiments, made the data evaluation and drafted the manuscript. RG and TA provided additional data and helped to conduct the experiments. NR provided plant cell suspension cultures, helped to conduct the experiments and reviewed the manuscript. SS coordinated the experiments and reviewed the manuscript. CL calculated oxygen diffusion coefficients and solubilities of the applied solutions. JB helped to design and coordinate the study and reviewed the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Klöckner, W., Gacem, R., Anderlei, T. et al. Correlation between mass transfer coefficient kLa and relevant operating parameters in cylindrical disposable shaken bioreactors on a bench-to-pilot scale. J Biol Eng 7, 28 (2013). https://doi.org/10.1186/1754-1611-7-28
- Dimensional analysis
- Maximum oxygen transfer capacity
- Orbitally shaken
- Plant cell suspension cultures
- Volumetric mass transfer coefficient