- Open Access
Development of a modified Respiration Activity Monitoring System for accurate and highly resolved measurement of respiration activity in shake flask fermentations
Journal of Biological Engineering volume 6, Article number: 11 (2012)
The Respiration Activity Monitoring System (RAMOS) is an established device to measure on-line the oxygen transfer rate (OTR), thereby, yielding relevant information about metabolic activities of microorganisms and cells during shake flask fermentations. For very fast-growing microbes, however, the RAMOS technique provides too few data points for the OTR. Thus, this current study presents a new model based evaluation method for generating much more data points to enhance the information content and the precision of OTR measurements.
In cultivations with E.coli BL21 pRSET eYFP-IL6, short diauxic and even triauxic metabolic activities were detected with much more detail compared to the conventional evaluation method. The decline of the OTR during the stop phases during oxygen limitations, which occur when the inlet and outlet valves of the RAMOS flask were closed for calibrating the oxygen sensor, were also detected. These declines reflected a reduced oxygen transfer due to the stop phases. In contrast to the conventional calculation method the new method was almost independent from the number of stop phases chosen in the experiments.
This new model based evaluation method unveils new peaks of metabolic activity which otherwise would not have been resolved by the conventional RAMOS evaluation method. The new method yields substantially more OTR data points, thereby, enhancing the information content and the precision of the OTR measurements. Furthermore, oxygen limitations can be detected by a decrease of the OTR during the stop phases.
Shake flasks are widely used in fermentations for biotechnological research and industrial process development[1, 2]. For gaining a better understanding and control of shake flask cultivations, various methods have been recently developed for online monitoring of process parameters.
Monitoring of pH-values in shake flasks has been realized both with standard autoclavable pH-probes that are immersed in the bulk liquid and with fluorescent optodes fixed at the flask wall that allow optical measurement[4–6]. Moreover, dissolved oxygen tension (DOT) can in principle be measured in shake flasks by using either Clark-type electrodes[7–10] or optical sensors based on dynamic quenching of luminescence[11–16]. These non-invasive measurement methods have proven to be more reliable, since they do not alter the hydrodynamics of the culture system due to baffling effects[17, 18].
As shown in numerous studies, almost all metabolic activities of aerobic microorganisms depend on the oxygen consumption of the culture[8, 19–21]. Therefore, online measuring techniques are also useful for determining the gas transfer rates in shake flasks. The company BlueSens GmbH (Herten, Germany) has developed a system that measures the gas transfer rates in the headspace of shake flasks down to nominal flask volumes of 500 mL.
As an alternative method, the Respiration Activity Monitoring System (RAMOS) presented by Anderlei et al.[23, 24], represents another non-invasive measurement technique which allows the online determination of oxygen transfer rate (OTR), carbon dioxide transfer rate (CTR) and respiratory quotient (RQ) in shake flasks down to 100 ml. Since its introduction, RAMOS has been used for multiple applications: e.g., for determining oxygen limitations in shake flask cultivations[18, 25–30], screening of microorganisms[31–33], optimizing media[34–36], investigating secondary substrate limitations[37, 38] and stress phenomena, process development and optimization[2, 40–44] and for monitoring precultures for fermentations in stirred tank reactors[45, 46]. The setup of a RAMOS device is illustrated in Figure1. For determining OTR with RAMOS, a certain measuring cycle is regularly repeated and consists of a rinsing phase and a stop phase[23, 24]. In the rinsing phase air is injected into the RAMOS flask to ensure a headspace gas concentration that equals that in a normal Erlenmeyer flask. During the stop phase the air flow into the system is briefly interrupted, and then the OTR is calculated from the corresponding drop (linear slope) in the oxygen partial pressure as detected by an oxygen sensor. Since these measuring cycles are typically repeated every 30 min, a data density of merely two OTR measurement points per hour is usually achieved.
The objective of this study is to present a new model based evaluation method for generating substantially more data points to enhance the information content and the precision of OTR measurements. First, the conventional way of generating OTR data is recalled. Then, the new evaluation method based on a complete oxygen headspace balance is developed. Microbial cultivations have been devised in order to evaluate the new method.
Results and discussion
Conventional evaluation method of the RAMOS device
In the conventional approach as proposed by Anderlei et al.[23, 24], the OTR in RAMOS is calculated during the stop phase based on equation A.1. The oxygen sensor does not directly provide the oxygen partial pressure pO2 but a voltage signal UO2 that depends on it. Therefore, the expression is introduced to convert this voltage signal into the oxygen partial pressure during the stop phase as described by the right part of equation A.1. Here, VG and VL are the volumes of the gas headspace and the liquid in the flask, respectively, and T is the temperature.
To compensate for the signal drift, the sensors are calibrated by calculating the partial pressure of oxygen at the end of the rinsing phase by applying the steady state gas composition of the headspace volume at 0 h. Figure2 gives an example of the sensor signal during a stop phase. This leads to equation A.2 for calculating the OTR:
Here, pO2,in is the oxygen partial pressure of the inlet flow, and is the volumetric flow into the flask during the rinsing phase, pamb is the ambient pressure and Vm is the molar volume.
A calibration factor (K) can then be calculated by using equation (A.3).
To examine the conventional evaluation method presented by Anderlei et al.[23, 24]E.coli BL21 pRSET eYFP-IL6 was cultivated in Wilms-MOPS-medium as described in “Materials and methods”. Figure3A illustrates the OTRs calculated using equation (A.2). Here, the OTR curve consists of one measuring point every 30 min, because the OTR is only evaluated during the stop phases. This curve shows a typical diauxic growth with the first peak leveling off and forming a horizontal plateau between 9.5 h and 10.5 h after inoculation. At the end of this phase glucose is fully consumed under oxygen limitation[23, 24]. At around 12.5 h after inoculation, a second peak appears indicating the consumption of acetate formed during the first phase of glucose consumption. This typical growth behavior has already been confirmed in other studies[6, 47, 48]. Therefore, the conventional evaluation method of RAMOS is able to measure the respiration activity of shake flask cultures and to detect phenomena as oxygen limitations and diauxic growth behavior, as long as these effects are lasting longer than the typical measuring cycle of 30 min. This condition is valid for many biological systems. However, a higher data density would be desirable when a highly resolved measurement is needed, e.g. for tracking metabolic activities that take place within much shorter periods.
Also illustrated in Figure3A is the calibration factor K as calculated by equation (A.3). This curve also shows two peaks which both occur right after the OTR peaks. These two peaks are totally unexpected, because the calibration factor K is assumed to be steady and independent from the OTR. Consequently, the assumption that the gas headspace volume at the end of the rinsing phase is in a steady state results in an inadequate calibration. Therefore, a calculation method should be developed that increases the data density of the OTR and includes a more accurate calibration.
New model based evaluation method
For increasing the data density of the OTR measurement the data from the rinsing phase should also be utilized for the evaluation. Therefore, the complete oxygen balance of the headspace of a RAMOS flask is derived and resolved for the OTR as follows (equation (B.1)):
and describe the volumetric flows into and out of the flask, respectively, pO2,in the oxygen partial pressure of the inlet flow, K is the calibration factor of the oxygen sensor, UO2 is the oxygen sensor signal, VG and VL are the volumes of the gas headspace and the liquid medium in the flask, respectively and T is the temperature.
The outlet flow is calculated from the slope of the pressure sensor signal of the headspace volume in the adjacent stop phases according to equation (B.2). This change in pressure during the stop phases occurs when the respiratory quotient (RQ) differs from 1. During the rinsing phases, when the valves are opened, this results in a change of the outlet flow. Equation (B.2) also takes into account the water vapor partial pressure that affects the outlet flow.
Here, pamb describes the ambient pressure and pH2O is the water vapor partial pressure in the gas headspace volume.
For determining the slope of the oxygen sensor signal, the sensor data of each rinsing phase are approximated with a smoothing spline function f that minimizes equation (B.3) according to the algorithm of de Boor. This algorithm is implemented into the Curve Fitting Toolbox of MATLAB R2010b (The Math Works, MA, USA). The left part of equation (B.3) consists of a least square method and is responsible for a close approximation of the sensor signal. In contrast, the right part of the equation represents the second derivative of the fit and, therefore, is responsible for the smoothness required for calculating reliable derivations. For adjusting a good trade-off between these two terms the regularization parameter p is used. A regularization parameter of p = 0.999 is suitable to describe the sensor signal of the rinsing phases and is, therefore, used to perform the fitting
p is the regularization parameter for the smoothing spline fitting and f describes the smoothing spline function.
For calculating the OTR in the stop phase, equation (B.1) simplifies to equation (B.4) since no air flow is applied to the flask. Therefore, and are zero and the left part of equation (B.1) vanishes.
The oxygen partial pressure in the stop phase is approximated with a second-order polynomial as described by equation (B.5). The reason for using a second-order polynomial equation is because it also considers a change of OTR during the stop phase and uses just one parameter more than a linear slope. As the stop phases are relatively short a more complicated equation resulting in a higher computational effort does not seem to be reasonable. Since the changes of the OTR are not restricted to the rinsing phase this method can be considered to be more accurate than using a linear slope. The first derivative of this polynomial is inserted into equation (B.4).
Due to the inertia of the oxygen sensor, there is a lag time for the sensor signal which does not allow a direct fit of equation (B.5). This lag time is compensated for by incorporating the first order delay model equation (B.6) into the fitting procedure. Thereby, the calculation of the coefficients a, b and c of the polynomial (B.5) is based on a least-square optimization procedure of these coefficients. Equation (B.6) is used to constantly calculate the corresponding sensor response UO2,delay on equation (B.5) by varying the estimates of a, b and c. As soon as the deviation between the calculated sensor response UO2,delay and the measured signal is small enough, the optimization stops and the course of the oxygen partial pressure can be calculated by using the corresponding estimates of the coefficients for the polynomial (B.5). Figure2 also shows the optimized polynomial of the stop phase according to equation (B.5) and the delayed signal of the stop phase according to the delay model described by equation (B.6). For the oxygen sensors used in this work, a time constant τ of 0.013 h is used (see Materials and methods). This value best describes the dynamic behavior of the sensor.
For solving equations (B.1) and (B.4), the calibration factor K is required. A new calibration method has been developed which also considers the dynamic behavior of the gas headspace volume instead of assuming a steady state. The end of the rinsing phase as well as the start of the subsequent stop phase of a biological experiment is utilized in the new sensor calibration method. This point occurs at 0 h in Figure2.
Figure2 also shows a section of the sensor signal curve during a RAMOS measurement and the mathematical fits including both the rinsing phase and the stop phase. Within the rinsing phase which last up to the beginning of the stop phase at 0 min, the sensor signal is described by a smoothing spline function f according to equation (B.3). During the stop phase, from 0 min to 5 min, the sensor signal strongly decreases. However, due to the lag time at the beginning of this stop phase, the initial slope of the signal curve merely changes slowly. Figure2, thus, shows both the polynomial described by equation (B.5) and the delayed signal described by equation (B.6) after the optimization procedure.
To calculate the calibration factor K, equations (B.1) and (B.4) are used. Whereas equation (B.1) describes the OTR during the rinsing phase and equation (B.4) the OTR during the stop phase, both equations can be equalized at the intersection of the curves of the rinsing phase and the stop phase as illustrated in Figure2. This is justified, because a sudden change in the OTR or K due to a change in the gas flow can be excluded. By resolving equations (B.1) and (B.4) under the assumption of their equality (OTRrinse = OTRstop), equation (B.7) is obtained to calculate K.
Using the same experiment, Figure3B depicts the OTR curve resulting from the new calculation method. Here, an OTR is calculated every 5 min. The OTR data from the stop phase are in good agreement with the OTR data calculated using the method of Anderlei et al.[23, 24]. The OTR data from the rinsing phase mostly concur with that of the stop phase. Only during the phase of oxygen limitation, indicated by the horizontal plateau of the OTR data at ca. 9.5 – 10.5 h[23, 24], the OTR in the stop phase is somewhat lower than that in the rinsing phase.
To understand this phenomenon, Figure4 is considered. The oxygen partial pressure of the gas headspace volume and the dissolved oxygen concentration in the medium of a simulation example are shown. The simulation is based on the model of equations (C.1) to (C.7) as shown in the appendix. This simulation example describes an oxygen limitation between 19 h and 21 h. The dissolved oxygen concentration cO2 was calculated by using the oxygen transfer equation (B.8). A constant mass transfer coefficient kLa was assumed.
Up to the beginning of the oxygen limitation at 19 h both the oxygen partial pressure of the headspace and the dissolved oxygen generally decrease due to increasing respiration activity. However, during the rinsing phases of this oxygen limitation between 19 h and 21 h the oxygen partial pressure of the headspace shows a horizontal plateau of ca. 0.18 bar whereas the dissolved oxygen has very low values which are close to the KO2 value of oxygen. Both the oxygen partial pressure of the headspace and the dissolved oxygen generally increase after the oxygen limitation when the respiration activity decreases.
Up to the oxygen limitation at 19 h, both, the oxygen partial pressure of the gas headspace volume and the dissolved oxygen strongly decrease during the stop phases as can be seen by the spikes in the signal. As both signals decrease the driving force of the oxygen transfer into the medium is not affected. However, during oxygen limitation between 19 h and 21 h the dissolved oxygen in the medium approaches very low values in the range of the KM value for oxygen and cannot decrease further. Nevertheless, the oxygen partial pressure of the headspace decreases, resulting in a lower driving force for the oxygen transfer into the medium and, consequently, in a slightly lower OTR during the stop phases, as depicted for 9.5 h–10.5 h in Figure3B.
In Figure3B also the calibration factor K as calculated with equation (B.7) is shown. In contrast to the steady state calibration using equation (A.3) as shown in Figure3A, the calibration factor K calculated with equation (B.7) has a steady course as expected and does not show any peaks correlated to the OTR. However, the data for K show some fluctuations which are due to the lower accuracy of data fitting at higher sensor dynamics. Therefore, the data for K are smoothed by using a smoothing spline with a regularization parameter of p = 0.1. This value has been found to be a good trade-off between smoothing of K and revealing the effect of sensor drift.
Recombinant E. coli fermentation in mineral medium
To verify the new evaluation method, E.coli BL21 pRSET eYFP-IL6 was cultivated using a RAMOS device with Wilms-MOPS-medium containing glucose and sorbitol as carbon sources. Four shake flasks were examined, each having a different initial sorbitol concentration. Figure5 shows the results. According to Monod, this cultivation leads to a diauxic growth where first glucose is consumed followed by a consumption of sorbitol. In all four cases, the OTR curves show an exponential growth phase during the 8 h after inoculation. Afterwards, the OTR curves form a horizontal plateau for up to 9.75 h after inoculation which indicates an oxygen limitation during this period of time[23, 24]. This oxygen limitation becomes more apparent when considering that the OTR of the stop phases are slightly lower during this period compared to the OTR of the rinsing phases. This difference is attributed to a reduced driving force of the oxygen transfer during the stop phases at oxygen limitation as described above. The oxygen concentration in the medium approaches very low values in the range of the KO2 value for oxygen, whereas the oxygen partial pressure of the headspace is still decreasing. Up to 9.75 h after inoculation when the oxygen limitation ends and the OTR decreases, glucose is completely consumed. Before glucose is depleted, the sorbitol concentrations remain constant and acetate is formed as an overflow metabolite of E. coli fermentation[6, 47, 48]. Furthermore, the optical density increases indicating biomass formation.
After 9.75 h, when the glucose is fully exhausted, sorbitol is consumed as can be seen by a decline in the sorbitol concentration. This sorbitol consumption is also observed in the OTR curve via a peak beginning at 10 h. The higher the initial concentration of sorbitol, the more pronounced the peak of the OTR curve is (from Figure5A to5D). Whereas the second OTR peak reaches a maximum of 0.039 mol/L/h at an initial sorbitol concentration of 0.5 g/L (Figure5A), it attains a maximum of 0.048 mol/L/h at an initial sorbitol concentration of 1 g/L (Figure5B) and a limiting value of 0.063 mol/L/h at an initial sorbitol concentration of 1.5 g/L (Figure5C). At an initial sorbitol concentration of 2 g/L (Figure5D), the OTR peak also reaches limiting values at about 0.062 mol/L/h, but is much wider. Thus, due to the higher data density, the new evaluation method unveils the real size of these peaks.
When sorbitol is depleted, the OTR curve again decreases before rising up to form a third peak which infers that acetate is now being consumed. This third peak has the same size for all four initial sorbitol concentrations, because the acetate is formed during the phase of glucose consumption and not during the sorbitol consumption. When acetate is consumed biomass is not produced anymore. Consequently, in the historical work of Monod, who measured biomass, only diauxic growth was observed and not a triauxic metabolic activity, as observed with the newly proposed calculation method of RAMOS data.
To show the strengths of the new evaluation method, three fermentation examples were considered where the OTR data of both methods are shown. Figure6A illustrates an enlarged section of the OTR curve of Figure5A. This OTR curve indicates an oxygen limitation between 8.5 h and 9.5 h as the OTR data points of the rinsing phases form a horizontal plateau during this time[23, 24]. Subsequently, the OTR drops sharply, indicating a reduced respiration activity until 11 h after inoculation. The OTR curve calculated with the new method generates a measurement point every 5 min, hence, offers a very high data density. Therefore, a second peak between 10 h and 10.5 h could be detected which is due to the consumption of sorbitol as could be shown in Figure5A. This shows that the newly proposed method is very useful in displaying short-term effects of the OTR.
The OTRs calculated with the new evaluation method in the stop phase are slightly higher than those calculated with the method of Anderlei et al.[23, 24]. This deviation is caused by the consideration of the sensor lag time in the new calculation method (equation (B.6)).
Using the new method, it is clearly observed that during the oxygen limitation between 8.5 h and 9.5 h the data points in the stop phase (filled squares) are not in line with the data points from the rinsing phase (open squares), in contrast to the rest of the fermentation. As discussed previously, this difference is due to a reduced driving force of the oxygen transfer during the stop phases of an oxygen limitation. The oxygen concentration in the medium approaches very low values, whereas the partial pressure of oxygen in the gas headspace volume still decreases. This effect does not occur in usual Erlenmeyer flasks, where the diffusive mass transfer through the cotton plug is not interrupted. Consequently, the OTR of the stop phase, gives a slightly different OTR than appearing in a typical Erlenmeyer flask.
Recombinant E. coli fermentation in complex medium
Figure6B illustrates the OTR curves of a RAMOS cultivation of E.coli BL21 pRSET eYFP-IL6 in TB-medium. Instead of the typically chosen value of 25 min, the rinsing phase was set to a duration of 55 min. However, the data density of the OTR curve calculated with the new evaluation method is not affected by the reduced number of stop phases because the data of the rinsing phase are also considered. After the lag phase at 7 hours the OTR shows an exponential growth up to 10.5 h leading directly into a horizontal plateau which can be considered as an oxygen limitation[23, 24]. Afterwards, even two very distinct peaks can be seen between 13 h and 16 h before the OTR drops down to a level of 0.005 mol/L/h at ca. 17 h. This clearly shows that the new method is also independent from the number of stop phases selected in the experiments. This characteristic is advantageous in applications where short-term effects are expected but a high number of stop phases is not desired, e.g. for preventing further oxygen limitations due to frequent interruption of the air flow.
Gluconobacter oxydans fermentation in complex medium
Figure6C illustrates the OTR curves of an RAMOS cultivation of Gluconobacter oxydans 621 H wild type in complex medium with 40 g/L mannitol. A rinsing phase of 25 min was selected. The OTR curve calculated with the method of Anderlei et al.[23, 24] indicates an exponential growth of the culture up to 7 h of cultivation when an OTR of 0.035 mol/L/h is achieved. The next measuring point at 7.5 h also indicates an OTR of 0.035 mol/L/h before the OTR decreases again. Consequently, the time between 7 h and 7.5 h could be interpreted as a short period of oxygen limitation. Even though the OTR curve calculated with the new method basically shows the same course, the time span between 7 h and 7.5 h shows considerably more data points, and more importantly, the OTR curve with the new calculation method does not show an oxygen limitation between the two stop phases.
The newly proposed evaluation method yields substantially more OTR data points than the conventional method by Anderlei et al.[23, 24]. This new evaluation method unveils additional peaks of metabolic activity which otherwise would remain undetected by the former method. Consequently, this new technique is a sophisticated means to generate more detailed information about metabolic activities of any kind of microorganisms and cells during shake flask cultivations. Additionally, possible oxygen limitations can be detected by a decrease of the OTR during the stop phases of the RAMOS measurement.
Materials and methods
E.coli BL21 pRSET eYFP-IL6 was maintained at −80°C in Lysogeny broth (LB) medium with 100 μg/mL ampicillin. Stock solutions contained 200 g/L glycerol. Gluconobacter oxydans 621 H wild type was maintained in its cultivation medium (see below) including 80 g/L mannitol at −80°C. Stock solutions contained 150 g/L glycerol.
LB medium for maintaining E.coli consists of: 5 g/L yeast extract (powdered, Roth, Karlsruhe, Germany), 10 g/L tryptone (pancreatic digest of casein, Roth, Karlsruhe, Germany) and 10 g/L NaCl. Terrific broth (TB) medium was used for cultivating both E.coli precultures and main cultures. The medium consists of: 5 g/L glycerol, 12 g/L tryptone (pancreatic digest of casein, Roth, Karlsruhe, Germany), 24 g/L yeast extract (powdered, Roth, Karlsruhe, Germany), 12.54 g/L K2HPO4 and 2.31 g/L KH2PO4. Additionally, 0.1 g/L ampicillin was added.
Main cultures of E.coli in mineral medium were cultivated in modified Wilms & Reuss synthetic medium (henceforth referred to as Wilms-MOPS medium)[6, 51]. This medium consists of: 20 g/L glucose, 5 g/L (NH4)2SO4, 0.5 g/L NH4Cl, 3 g/L K2HPO4, 2 g/L Na2SO4, 0.5 g/L MgSO4·7H2O, 41.85 g/L (0.2 M) 3-(N-morpholino)-propanesulfonic acid (MOPS), 0.01 g/L thiamine hydrochloride, 1 mL/L trace element solution (0.54 g/L ZnSO4·7H2O, 0.48 g/L CuSO4·5H2O, 0.3 g/L MnSO4·H2O, 0.54 g/L CoCl2·6H2O, 41.76 g/L FeCl3·6H2O, 1.98 g/L CaCl2·2H2O, 33.39 g/L Na2EDTA (Titriplex III)) and concentrations of sorbitol ranging from 0 – 2 g/L. Additionally, 0.1 g/L ampicillin was added.
Gluconobacter oxydans was cultivated in a complex medium consisting of: 5 g/L yeast extract (powdered, Roth, Karlsruhe, Germany), 1 g/L K2HPO4, 1 g/L (NH4)2SO4, 0.5 g/L glycerol, 2.5 g/L MgSO4·7H2O and 40 g/L mannitol. The pH-value was adjusted to 6 with 1 M KOH.
E. coli precultures and main cultures were cultivated in 250 ml shake flasks at a temperature of T = 30°C with a shaking diameter of d0 = 50 mm and a shaking frequency of n = 350 rpm (shaking machine LS-W, Kuehner AG, Birsfelden, Switzerland). For E. coli precultures, 20 ml of TB medium was inoculated with 200 μl of stock solution. Main cultures of E.coli in 10 mL Wilms-MOPS medium were inoculated with preculture broth from the exponential growth phase, resulting in an optical density (OD) of 0.2.
The main culture of E.coli in TB medium was cultivated in 250 ml shake flasks with a filling volume of VL = 10 mL, a shaking diameter of d0 = 50 mm, a shaking frequency of n = 300 rpm, and at a temperature of T = 30°C. This culture was inoculated with 200 μl of stock solution.
The preculture and the main culture of Gluconobacter oxydans were cultivated in 250 ml shake flasks with a filling volume of VL = 10 mL, a shaking diameter of d0 = 50 mm, a shaking frequency of n = 350 rpm, and at a temperature of T = 30°C. The preculture was inoculated with 500 μL of stock solution. The main culture was inoculated with preculture broth from the exponential growth phase resulting in an OD of 0.1.
Respiration Activity Monitoring System (RAMOS)
The respiration activity was measured in modified 250 mL Erlenmeyer flasks using a self-made RAMOS device as introduced by Anderlei et al.[23, 24]. The setup of the RAMOS device is shown in Figure1. The 8 RAMOS flasks in parallel were supplied with air using a mass flow controller (type 5850 TR, Brooks, Hatfield, PA, USA). The air flow through the flasks is controlled by using valves at the inlet and the outlet of each RAMOS flask. For the rinsing phase the RAMOS flasks were flushed for 25 min with air at a flow rate of (low flow), which corresponds to an aeration rate of 1 vvm at a filling volume of VL = 10 mL. The rinsing phase was followed by a 5 min stop phase (tstop) with no air flow through the flasks. To compensate for the drop in oxygen partial pressure during the stop phase, a higher air flow of (high flow) was applied during the first 0.9 min of the rinsing phase. This measuring cycle is regularly repeated. The oxygen partial pressure of the headspace in a RAMOS flask was measured with a MAX250 oxygen sensor from Maxtec (Salt Lake City, Utah, USA). The difference between the headspace pressure and ambient pressure was detected with a 26PCA pressure sensor from Honeywell (Morristown, NJ, USA). The RAMOS flasks ensure the same hydrodynamic conditions and headspace gas concentrations as are found in regular Erlenmeyer flasks with cotton plugs. Commercial versions of the RAMOS device are available from Kuehner AG, Birsfelden, Switzerland and Hitech Zang, Herzogenrath, Germany.
The sensor lag time τ has been determined in step change experiments. Beginning at a steady state of air in the gas headspace volume, oxygen depleted air was flushed into a RAMOS flask under well-defined conditions ( = 0.54 mL/min, pO2,in = 0.197 bar, VL = 10 mL, n = 300 rpm). By calculating the real partial pressure of oxygen in the gas headspace of the flask and then optimizing the sensor lag time τ, the model described by equation (B. 6) was fitted to the oxygen partial pressure signal measured by the sensor. This leads to a sensor lag time of τ = 0.013 h (data not shown).
Parallel shake flask cultivations
Samples were taken from E. coli cultivations in Wilms-MOPS medium in Erlenmeyer flasks in parallel to the RAMOS experiments and cultivated under the same conditions as used in the RAMOS experiments. The OD of the samples at 600 nm was determined with a Thermo Scientific Genesys 20 spectrophotometer (Waltham, MA, USA).
For determining glucose, sorbitol and acetate, the samples were centrifuged for 5 min at 14000 rpm with a Sigma 1–15 Microfuge (Osterode am Harz, Germany) and the supernatants were used for the analysis. Concentrations of glucose, sorbitol and acetate in the respective supernatants were determined using a Dionex HPLC (Dionex, Sunnyvale, USA) with an Organic Acid-Resin 300 × 8 mm (CSChromatographie, Langerwehe, Germany) and a Skodex RI-71 detector. Sulfuric acid in a concentration of 5 mM was used as solvent at a flow rate of 0.6 ml/min and a temperature of 60°C.
OTRs, the calibration factor (K) and simulation examples were calculated with MATLAB R2010b (The Math Works, MA, USA). The mathematical optimization problems were solved using the trust-region-reflective algorithm. Ordinary differential equations where solved using the trapezoidal rule.
Model used for calculating the data of Figure4:
Dissolved oxygen equation:
Oxygen partial pressure equation:
Oxygen transfer rate:
Dissolved oxygen at the gas–liquid interface:
With: Oxygen solubility LO2 = 0.0011 mol/L/bar, mass transfer coefficient kLa = 0.082 1/s, volumetric flow into the flask = 0.67 L/h, volumetric flow out of the flask = 0.7 L/h, maximal growth rate μmax = 0.26 1/h, half velocity constant of the substrate KS = 4 g/L, half velocity constant of oxygen KO2 = 8 ·10-10 mol/L, oxygen partial pressure of inlet pO2,in = 0.2095 bar, yield coefficient of the substrate YXS = 0.5 g/g, liquid volume VL = 10 mL, gas volume VG = 270 mL, gas constant R = 0.0831 L*bar/mol/K, yield coefficient of oxygen YXO2 = 41.6 g/mol, temperature T = 303.15 K, initial conditions: biomass concentration Xo = 0.5 g/L, substrate concentration S0 = 45 g/L, oxygen concentration in the liquid cO2,0 = 0.22 mmol/L, oxygen partial pressure pO2,0 = 0.2006 bar.
a, b, c: Polynomial coefficients; cO2: Dissolved oxygen in the liquid (mmol/L); cO2,0: Initial dissolved oxygen in the liquid (mmol/L); cO2*: Dissolved oxygen at the gas–liquid interface (mmol/L); f: Smoothing spline; kLa: Mass transfer coefficient (1/s); K: Calibration factor (bar/mV); KO2: half velocity constant for oxygen (mol/L); KS: Half velocity constant for the substrate (g/L); LO2: Oxygen solubility (mol/L/bar); OTR: Oxygen transfer rate (mol/L/h); p: Regularization parameter (-); pamb: Ambient pressure (bar); pH2O: Water vapor partial pressure in the headspace (bar); pO2: Oxygen partial pressure in the headspace (bar); pO2,0: Initial oxygen partial pressure in the headspace (bar); pO2,in: Oxygen partial pressure of inlet flow (bar); pO2,real: Steady state oxygen partial pressure (bar); R: Gas constant (bar L/mol/K); RQ: Respiratory quotient (-); S: Substrate concentration (g/L); S0: Initial substrate concentration (g/L); T: Temperature (K); UO2: Oxygen sensor signal (mol/L/h); VG: Gas volume (L); VL: Liquid volume (L); Vm: Molar volume (L/mol);: Volumetric flow into the flask (L/h);: Volumetric flow out of the flask (L/h); X: Biomass concentration (g/L); X0: Initial biomass concentration (g/L); YXO2: Yield coefficient for oxygen (g/mol); YXS: Yield coefficient for the substrate (g/g); μ: Growth rate (1/s); μmax: Maximal growth rate (1/s); t: Lag time of oxygen sensor (h).
Büchs J: Introduction to advantages and problems of shaken cultures. Biochemical Engineering Journal 2001, 7: 91-98. 10.1016/S1369-703X(00)00106-6
Freyer SA, König M, Künkel A: Validating shaking flasks as representative screening systems. Biochemical Engineering Journal 2004, 17: 169-173. 10.1016/S1369-703X(03)00175-X
Weuster-Botz D, Altenbach-Rehm J, Arnold M: Parallel substrate feeding and pH-control in shaking-flasks. Biochemical Engineering Journal 2001, 7: 163-170. 10.1016/S1369-703X(00)00117-0
Kermis HR, Kostov Y, Harms P, Rao G: Dual excitation ratiometric fluorescent pH sensor for noninvasive bioprocess monitoring: Development and application. Biotechnology Progress 2002, 18: 1047-1053. 10.1021/bp0255560
Badugu R, Kostov Y, Rao G, Tolosa L: Development and Application of an Excitation Ratiometric Optical pH Sensor for Bioprocess Monitoring. Biotechnology Progress 2008, 24: 1393-1401. 10.1002/btpr.66
Scheidle M, Klinger J, Büchs J: Combination of on-line pH and oxygen transfer rate measurement in shake flasks by fiber optical technique and Respiration Activity Monitoring System (RAMOS). Sensors 2007, 7: 3472-3480. 10.3390/s7123472
Hirose Y, Sonoda H, Kinoshita K, Okada H: Studies on oxygen transfer in submerged fermentations Part 4 Determination of oxygen transfer rate and respiration rate in shaken cultures using oxygen analyzers. Agricultural and Biological Chemistry 1966, 30: 49-58. 10.1271/bbb1961.30.49
McDaniel LE, Bailey EG: Effect of shaking speed and type of closure on shake flask cultures. Applied Microbiology 1969, 17: 286-290.
Ruottinen M, Bollok M, Kogler M, Neubauer A, Krause M, Hamalainen ER, Myllyharju J, Vasala A, Neubauer P: Improved production of human type II procollagen in the yeast Pichia pastoris in shake flasks by a wireless-controlled fed-batch system. Bmc Biotechnology 2008, 8: 33. 10.1186/1472-6750-8-33
Vasala A, Panula J, Bollok M, Illmann L, Halsig C, Neubauer P: A new wireless system for decentralised measurement of physiological parameters from shake flasks. Microbial Cell Factories 2006, 5: 8. 10.1186/1475-2859-5-8
Klimant I, Wolfbeis OS: Oxygen-sensitive luminescent materials based on silicone-soluble ruthenium diimine complexes. Analytical Chemistry 1995, 67: 3160-3166. 10.1021/ac00114a010
Tolosa L, Kostov Y, Harms P, Rao G: Noninvasive measurement of dissolved oxygen in shake flasks. Biotechnology and Bioengineering 2002, 80: 594-597. 10.1002/bit.10409
Gupta A, Rao G: A study of oxygen transfer in shake flasks using a non-invasive oxygen sensor. Biotechnology and Bioengineering 2003, 84: 351-358. 10.1002/bit.10740
Wittmann C, Kim HM, John G, Heinzle E: Characterization and application of an optical sensor for quantification of dissolved O-2 in shake-flasks. Biotechnology Letters 2003, 25: 377-380. 10.1023/A:1022402212537
Becker J, Klopprogge C, Wittmann C: Metabolic responses to pyruvate kinase deletion in lysine producing Corynebacterium glutamicum. Microbial Cell Factories 2008, 7: 8. 10.1186/1475-2859-7-8
Schneider K, Schuetz V, John GT, Heinzle E: Optical device for parallel online measurement of dissolved oxygen and pH in shake flask cultures. Bioprocess and Biosystems Engineering 2010, 33: 541-547. 10.1007/s00449-009-0367-0
Suresh S, Srivastava VC, Mishra IM: Techniques for oxygen transfer measurement in bioreactors: a review. Journal of Chemical Technology and Biotechnology 2009, 84: 1091-1103. 10.1002/jctb.2154
Hansen S, Kensy F, Kaeser A, Büchs J: Potential errors in conventional DOT measurement techniques in shake flasks and verification using a rotating flexitube optical sensor. BMC Biotechnol 2011, 11: 49-49. 10.1186/1472-6750-11-49
Calik P, Yilgör P, Ayhan P, Demir AS: Oxygen transfer effects on recombinant benzaldehyde lyase production. Chemical Engineering Science 59: 5075-5083.
Garcia-Ochoa F, Castro EG, Santos VE: Oxygen transfer and uptake rates during xanthan gum production. Enzyme and Microbial Technology 2000, 27: 680-690. 10.1016/S0141-0229(00)00272-6
Liu Y-S, Wu J-Y, Ho K-p: Characterization of oxygen transfer conditions and their effects on Phaffia rhodozyma growth and carotenoid production in shake-flask cultures. Microbiology 2009, 155: 3827-3837. 10.1099/mic.0.030213-0
Heyland J, Fu J, Blank L: Correlation between TCA cycle flux and glucose uptake rate during respiro-fermentative growth of Saccharomyces cerevisiae. Biochemical Engineering Journal 2001, 7: 157-162. 10.1016/S1369-703X(00)00116-9
Anderlei T, Büchs J: Device for sterile online measurement of the oxygen transfer rate in shaking flasks. Biochemical Engineering Journal 2001, 7: 157-162. 10.1016/S1369-703X(00)00116-9
Anderlei T, Zang W, Papaspyrou M, Büchs J: Online respiration activity measurement (OTR, CTR, RQ) in shake flasks. Biochemical Engineering Journal 2004, 17: 187-194. 10.1016/S1369-703X(03)00181-5
Peña C, Peter CP, Büchs J, Galindo E: Evolution of the specific power consumption and oxygen transfer rate in alginate-producing cultures of Azotobacter vinelandii conducted in shake flasks. Biochemical Engineering Journal 2007, 36: 73-80. 10.1016/j.bej.2007.02.019
Guez JS, Müller CH, Danze PM, Büchs J, Jacques P: Respiration activity monitoring system (RAMOS), an efficient tool to study the influence of the oxygen transfer rate on the synthesis of lipopeptide by Bacillus subtilis ATCC6633. Journal of Biotechnology 2008, 134: 121-126. 10.1016/j.jbiotec.2008.01.003
Zimmermann H, Anderlei T, Büchs J, Binder M: Oxygen limitation is a pitfall during screening for industrial strains. Applied Microbiology and Biotechnology 2006, 72: 1157-1160. 10.1007/s00253-006-0414-6
Mehmood N, Olmos E, Goergen JL, Blanchard F, Ullisch D, Klöckner W, Büchs J, Delaunay S: Oxygen supply controls the onset of pristinamycins production by Streptomyces pristinaespiralis in shaking flasks. Biotechnology and Bioengineering 2011, 108: 2151-2161. 10.1002/bit.23177
Seletzky JM, Noak U, Fricke J, Welk E, Eberhard W, Knocke C, Büchs J: Scale-up from shake flasks to fermenters in batch and continuous mode with Corynebacterium glutamicum on lactic acid based on oxygen transfer and pH. Biotechnology and Bioengineering 2007, 98: 800-811. 10.1002/bit.21359
Maier U, Losen M, Büchs J: Advances in understanding and modeling the gas–liquid mass transfer in shake flasks. Biochemical Engineering Journal 2004, 17: 155-167. 10.1016/S1369-703X(03)00174-8
Stöckmann C, Losen M, Dahlems U, Knocke C, Gellissen G, Büchs J: Effect of oxygen supply on passaging, stabilising and screening of recombinant Hansenula polymorpha production strains in test tube cultures. FEMS Yeast Research 2003, 4: 195-205. 10.1016/S1567-1356(03)00147-8
Scheidle M, Jeude M, Dittrich B, Denter S, Kensy F, Suckow M, Klee D, Büchs J: High-throughput screening of Hansenula polymorpha clones in the batch compared with the controlled-release fed-batch mode on a small scale. FEMS Yeast Research 2010, 10: 83-92. 10.1111/j.1567-1364.2009.00586.x
Stöckmann C, Scheidle M, Dittrich B, Merckelbach A, Hehmann G, Melmer G, Klee D, Büchs J, Kang HA, Gellissen G: Process development in Hansenula polymorpha and Arxula adeninivorans, a re-assessment. Microbial Cell Factories 2009, 8: 22. 10.1186/1475-2859-8-22
Losen M, Frölich B, Pohl M, Büchs J: Effect of Oxygen Limitation and Medium Composition on Escherichia coli Fermentation in Shake-Flask Cultures. Biotechnology Progress 2004, 20: 1062-1068. 10.1021/bp034282t
Klement T, Dankmeyer L, Hommes R, van Solingen P, Büchs J: Acetate-glycerol cometabolism: Cultivating Schizosaccharomyces pombe on a non-fermentable carbon source in a defined minimal medium. Journal of Bioscience and Bioengineering 2011, 112: 20-25. 10.1016/j.jbiosc.2011.02.014
Kunze M, Huber R, Gutjahr C, Müllner S, Büchs J: Predictive tool for recombinant protein production in E. coli shake flask cultures using an on-line monitoring system. Biotechnology Progress 2011, 28: 103-113.
Kottmeier K, Weber J, Müller C, Bley T, Büchs J: Asymmetric division of Hansenula polymorpha reflected by a drop of light scatter intensity measured in batch microtiter plate cultivations at phosphate limitation. Biotechnology and Bioengineering 2009, 104: 554-561. 10.1002/bit.22410
Huber R, Roth S, Rahmen N, Büchs J: Utilizing high-throughput experimentation to enhance specific productivity of an E.coli T7 expression system by phosphate limitation. Bmc Biotechnology 2011, 11: 22. 10.1186/1472-6750-11-22
Seletzky J, Noack U, Fricke J, Hahn S, Büchs J: Metabolic activity of Corynebacterium glutamicum grown on L-lactic acid under stress. Applied Microbiology and Biotechnology 2006, 72: 1297-1307. 10.1007/s00253-006-0436-0
Huber R, Scheidle M, Dittrich B, Klee D, Büchs J: Equalizing growth in high-throughput small scale cultivations via precultures operated in fed-batch mode. Biotechnology and Bioengineering 2009, 103: 1095-1102. 10.1002/bit.22349
Scheidle M, Dittrich B, Klinger J, Ikeda H, Klee D, Büchs J: Controlling pH in shake flasks using polymer-based controlled-release discs with pre-determined release kinetics. Bmc Biotechnology 2011, 11: 25. 10.1186/1472-6750-11-25
Rechmann H, Friedrich A, Forouzan D, Barth S, Schnabl H, Biselli M, Boehm R: Characterization of photosynthetically active duckweed (Wolffia australiana) in vitro culture by Respiration Activity Monitoring System (RAMOS). Biotechnology Letters 2007, 29: 971-977. 10.1007/s10529-007-9330-9
Jänicke G, Sauter C, Bux R, Haas J: Characterisation of Shake Flasks for Cultivation of Animal Cell Cultures Cell Technology for Cell Products. In Cell Technology for Cell Products. Edited by: Smith R. Netherlands: Springer; 2007:727-731.
Seletzky J, Noack U, Hahn S, Knoll A, Amoabediny G, Büchs J: An experimental comparison of respiration measuring techniques in fermenters and shake flasks: exhaust gas analyzer vs. RAMOS device vs. respirometer. Journal of Industrial Microbiology & Biotechnology 2007, 34: 123-130. 10.1007/s10295-006-0176-2
Knoll A, Bartsch S, Husemann B, Engel P, Schroer K, Ribeiro B, Stöckmann C, Seletzky J, Büchs J: High cell density cultivation of recombinant yeasts and bacteria under non-pressurized and pressurized conditions in stirred tank bioreactors. Journal of Biotechnology 2007, 132: 167-179. 10.1016/j.jbiotec.2007.06.010
Mrotzek C, Anderlei T, Henzler H-J, Büchs J: Mass transfer resistance of sterile plugs in shaking bioreactors. Biochemical Engineering Journal 2001, 7: 107-112. 10.1016/S1369-703X(00)00108-X
Xu B, Jahic M, Enfors S-O: Modeling of Overflow Metabolism in Batch and Fed-Batch Cultures of Escherichia coli. Biotechnology Progress 1999, 15: 81-90. 10.1021/bp9801087
Christensen ML, Eriksen NT: Growth and proton exchange in recombinant Escherichia coli BL21. Enzyme and Microbial Technology 2002, 31: 566-574. 10.1016/S0141-0229(02)00153-9
De Boor C: A Practical guide to splines. New York: Springer; 1978.
Monod J: The Growth of Bacterial Cultures. Annual Review of Microbiology 1949, 3: 371-394. 10.1146/annurev.mi.03.100149.002103
Wilms B, Hauck A, Reuss M, Syldatk C, Mattes R, Siemann M, Altenbuchner J: High-cell-density fermentation for production of L-N-carbamoylase using an expression system based on the Escherichia coli rhaBAD promoter. Biotechnology and Bioengineering 2001, 73: 95-103. 10.1002/bit.1041
We wish to thank Udo Koesfeld, Thomas Heise and René Petri for their essential technical support. Moreover, we are grateful to Mary-Joan Bluemich for the many hours we spent on editing this paper.
The authors declare that they have no competing interests.
SH drafted the manuscript and performed the calculations. SH, IH and BL performed the experiments. Having conceived the study, JB helped to design, coordinate and write the manuscript. All authors read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Hansen, S., Hariskos, I., Luchterhand, B. et al. Development of a modified Respiration Activity Monitoring System for accurate and highly resolved measurement of respiration activity in shake flask fermentations. J Biol Eng 6, 11 (2012). https://doi.org/10.1186/1754-1611-6-11
- Oxygen transfer rate
- Shake flask
- Bioprocess monitoring
- Diauxic growth
- Oxygen limitation