Next Article in Journal
Continuous-Flow Biocatalytic Process for the Synthesis of the Best Stereoisomers of the Commercial Fragrances Leather Cyclohexanol (4-Isopropylcyclohexanol) and Woody Acetate (4-(Tert-Butyl)Cyclohexyl Acetate)
Next Article in Special Issue
Valorization of Glycerol through the Enzymatic Synthesis of Acylglycerides with High Nutritional Value
Previous Article in Journal
Structural Properties and Catalytic Activity of Binary Poly (vinyl alcohol)/Al2O3 Nanocomposite Film for Synthesis of Thiazoles
Previous Article in Special Issue
Synthesis of Ibuprofen Monoglyceride in Solventless Medium with Novozym®435: Kinetic Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Kinetic Modeling of Dihydroxyacetone Production from Glycerol by Gluconobacter oxydans ATCC 621 Resting Cells: Effect of Fluid Dynamics Conditions

by
Susana de la Morena
,
Mateusz Wojtusik
,
Victoria E. Santos
* and
Felix Garcia-Ochoa
Chemical Engineering and Materials Department, Universidad Complutense de Madrid, 28040 Madrid, Spain
*
Author to whom correspondence should be addressed.
Catalysts 2020, 10(1), 101; https://doi.org/10.3390/catal10010101
Submission received: 14 October 2019 / Revised: 3 January 2020 / Accepted: 7 January 2020 / Published: 10 January 2020
(This article belongs to the Special Issue New Glycerol Upgrading Processes)

Abstract

:
Dihydroxyacetone production from glycerol has been studied. Cultures of Gluconobacter oxydans ATCC 621, a promising microorganism that is able to convert glycerol into dihydroxyacetone, has been employed. In this work, the influence of oxygen transport rate and the fluid dynamic conditions have been studied working with resting cells cultures. Several experiments were carried out at two different scales: 250 mL Erlenmeyer flasks and a 2 L stirred tank bioreactor, varying the agitation speed. Product and substrate concentration were determined employing high-performance liquid chromatography. Additionally, oxygen concentration was measured in the runs carried out in stirred tank reactors. Taking into account the results obtained in these experiments, three different behaviors were observed, depending on the mass transfer and chemical reactions rates. For experiments with low stirring speed (below 200 rpm for shake flasks and 300 rpm for reactors), the oxygen transport rate is the controlling step, while at high stirring speed (over 300 rpm in shake flasks and 560 rpm in the bioreactor), the chemical reaction is controlling the overall process rate. In some runs conducted at medium agitation, a mix control was found. All the kinetic models were able to reproduce experimental data and fulfill thermodynamic and statistical criteria, highlighting the importance of the mass transfer rate upon this system.

1. Introduction

Glycerol has become a waste due to the enormous availability caused by the huge growth of the biodiesel production industry [1], which produces glycerol as a by-product; approximately 10 kg of glycerol are produced per 100 kg of biodiesel [2]. An estimation of OECD-FAO (Organization for Economic Co-operation and Development—Food and Agricultural Organization) assumes that the global production of biodiesel will grow 4.5% annually, reaching 41 Mm3 in 2022 [2]. This makes the development of different sustainable processes that are able to transform the growing amount of glycerol produced into valuable products peremptory.
Glycerol has been employed in many industries (food, pharmaceutical, cosmetic, chemical, etc.), due to its versatile physicochemical properties [3]. Amongst all these applications, one promising route is the bioconversion of glycerol into chemical products and precursors, such as 1,3-dihydroxyacetone (DHA). DHA is extensively used in the cosmetic industry due to its employment as a building block for the production of various fine chemicals such as 1,3-propanediol and citric acid (being an intermediate compound in this processes) [4], with microbial production being the preferred synthesis pathway, in contrast with the chemical processes, which have expensive safety requirements [5].
The most widely employed microorganism to produce DHA from glycerol as raw material is Gluconobacter oxydans, although there are several different works using other microorganisms [5,6,7,8,9,10,11,12,13,14,15,16]. G. oxydans is a Gram-negative aerobic bacterium, whose respiratory metabolism allows it to partially oxidize different organic compounds [17]. The enzyme that transforms glycerol into DHA, glycerol dehydrogenase (corresponding to the classification EC 1.1.1.6), is located in the microorganism membrane [18], and has been employed as a separate catalyst in immobilized configurations to produce DHA [19,20]. This provides G. oxydans with the ability to transform glycerol into DHA rapidly and efficiently. However, as the whole process is intimately associated to the oxygen caption, and this compound is the final donor of electrons in this process, it is crucial to evaluate the influence of oxygen transport from air to the culture media and to the membrane of the microorganism [16].
In the literature, there have been some efforts to establish a valid kinetic model to describe the transformation of glycerol into DHA. In 2003, Hekmat et al. [5] employed a semi-continuous stirred tank bioreactor combined with a column with irregular packing, in which the microorganism was immobilized, in order to minimize operational costs such as cleaning, sterilization, or inoculation. In order to fix the optimal operating conditions, the authors evaluated several kinetic parameters for semi-continuous systems, describing the rate of growth by means of a double substrate limitation Monod model, coupled with a linear product inhibition. Although the experimental approach improved final DHA yield and productivity, the kinetic model was not able to explain some phenomena associated to this process, and the kinetic model employed was purely empirical [5].
Ma et al., in 2009, took a different approach in order to describe DHA production, using a modified Aiba model for cell growth and the Leudeking–Piret growth associate model for the product formation description [21]; in this case, the inhibitory effect of both DHA and glycerol presence on biomass growth and therefore DHA production were considered. Additionally, Ma et al. compared two different G. oxydans strains: the wild type and a genetically modified strain using He–Ne laser irradiation technology [21]. The kinetic model proposed was able to describe correctly the experimental results achieved, both for the natural and for the modified strains.
Other authors also describe a strong substrate and product inhibition on DHA production [12,13,15,22]. Furthermore, Dikshit et al. employed various different kinetic models in order to describe DHA production considering the effect of substrate inhibition [23] and substrate and product combined inhibition [24]. In this work, an immobilized G. oxydans in reticulated polyurethane was used, studying the effect of raw glycerol over the kinetic parameters and both substrate and product inhibitions [23,24].
However, none of the previously described kinetic models consider the oxygen availability. The importance of oxygen transfer rate in aerobic the bioprocess has already been pointed out in many works [16,25,26,27,28]. Recently, several authors have put their effort on describing the influence of oxygen transfer rate on DHA production. Zheng et al. compared several different kinetic parameters for three different values of the oxygen volumetric mass transfer coefficient (kLa), using crude glycerol as the carbon source and employing Gluconobacter frateurii CGMCC 5397 as the biocatalyst [8]. The results shown in this work demonstrate that an increment on mass transfer rate by improving kLa leads to an increase of both biomass growth and DHA production rates [8]. Poljungreed et al. studied the effect of the amount of dissolved oxygen on DHA final concentration and yield in a 7 L stirred tank bioreactor (STBR) using G. frateurii BCC 36199 strain cultured in glycerol and inorganic salts medium. The results shown an improvement in both factors—final concentration and yield—which was directly related with the amount of available oxygen [11].
Regarding the strain employed in this work, previously published results [16] have shown that there is a need to establish a minimum oxygen transfer rate in order to ensure the highest culture activity. Thus, we can avoid the decrease of DHA production rate due to the oxygen transfer rate (OTR) limitation and its control over the overall process rate [16].
The aim of this work is to establish a kinetic model that is able to correlate OTR with DHA production in two different scales: low-scale 250 mL flasks and medium-scale 2 L STBR. In order to achieve this objective, several experiments were carried out both in Erlenmeyer flasks and STBR, varying the agitation speed. Operating conditions were fixed, taking into account the results achieved in previous works [15,16]. Throughout the study, resting cells were employed in operation mode. With this kinetic model, it will be possible to establish a clear dependence between OTR and DHA productivity, highlighting even more the importance of the mass transfer phenomenon in this process. Mass transfer rate has not been considered in previous works studying this bioprocess. Moreover, the improvement of OTR would lead to an increase in the DHA productivity and in the overall feasibility of the process.

2. Results

The experiments realized in this work are shown in Table 1. Based on the results obtained in the aforementioned runs, three different kinetic modeling approaches have been used in order to describe DHA concentration evolution alongside each experiment, considering the different behaviors observed in each case. In all runs, as it was observed in a previous work [16], the glyceric acid (a by-product) concentration was negligible; therefore, this acid was not considered in order to develop the kinetic models in this work.

2.1. Shake Flask Experiments

Figure 1, Figure 2, Figure 3 and Figure 4 show the evolution of both substrate (glycerol) and product (DHA) concentrations for the different experiments carried out at low-scale (shake flasks). A clear linear tendency with time for DHA production as well as for glycerol consumption can be observed in these figures, at least at short times. However, only at the lowest agitation (100 rpm), this linear tendency is maintained throughout the run; while in the rest of the runs, the linear trend is only observed at shorter times. Furthermore, it is important to note that at the lowest agitation speed, both evolutions—glycerol and DHA concentrations—are independent of the biomass concentration employed, which is a clear sign of the existence of a mass transfer limitation phenomenon affecting the DHA production rate. Finally, it is necessary to point out that the glycerol is totally consumed during the experimental time in all of the runs carried out in this work.
On the other hand, as can be seen in Figure 2, Figure 3 and Figure 4, the linear tendency commented does not last during the entirety of the experiments represented in these figures. Although the linear phase is the same for each of the stirring speeds over all the runs, both the slope and duration of this phase vary with agitation speed. Thus, stirring speed—and therefore, mass transfer rate—had an important impact on DHA production rate for the runs shown in Figure 2, Figure 3 and Figure 4.
This effect can also be stated over the specific productivity values (PX), which are shown in Figure 5. In this figure, it can be appreciated that the effect of the stirring speed is more impactful over PX for the amount of biomass employed in each case, which is in line with the observation previously described for the experimental results. In addition, these results are comparable to the work present in the literature. Jackson et al. reached a PX of 3.21 × 10−2 mol·gX−1·h−1 employing an N of 180 rpm and a biomass concentration of 0.66 gX·L−1 using commercial glycerol and G. oxydans resting cells [29], whereas in this work, we have reached a value of 2.25 × 10−2 mol·gX−1·h−1 for similar conditions. However, as can be seen in Figure 5, this productivity can be increased substantially through incrementing stirring rate and, therefore, favoring OTR.

2.2. STBR Experiments

The results of the runs performed at STBR scale are given in Figure 6. As can be seen, both glycerol and DHA concentrations change with the stirring speed employed in each run. For the lowest agitation (300 rpm), both glycerol and DHA concentration maintained a linear time course trend.
However, as stirring increases, the evolutions of both compound concentrations change into hyperbolic type, in which the first part of the time-course evolution for both concentrations resembles a linear evolution with a higher slope than the first experiment at 300 rpm, and it is the same tendency observed in the runs carried out in shake flasks.
For the runs with higher agitation, DHA evolution is the same, which is maintained despite the stirring speed increase. It is important to point out that for all the experiments shown in Figure 6, the biomass concentration employed remains the same (0.5 g·L−1).
As with the experiments carried out in shake flasks, PX was evaluated for each run, and the results are shown in Table 2. As can be seen in this table, the productivity reached for these experiments is lower than that obtained in flask experiments. However, it is clear that the same tendency can be observed in both cases, as the productivity increases with the stirring rate.

3. Discussion

The results observed in both of the experimental devices employed leads to the following considerations, which have to be addressed in order to formulate an adequate kinetic model to describe DHA production:
  • There are three different observed patterns, which depend on the stirring speed employed. This is valid for both shake flasks and STBR.
  • For low agitation, DHA production is strongly related to the stirring speed and, therefore, to the mass transfer rate.
  • For high stirring speeds, biomass concentration exerts a strong influence on DHA evolution. When the amount of biocatalyst is maintained, DHA evolution is not affected by agitation speed. Thus, the production rate of DHA depends on the chemical reaction rate.
  • There is a mixed behavior when the agitation takes a medium value; in this case, the process rate seems to depend on the oxygen availability, or the OTR from gas to the culture media, which is directly related to the stirring speed.
Taking into account these factors, several different kinetic models were proposed to describe each one of the situations described above.

3.1. Kinetic Model for Low Agitation: Mass Transfer Rate Evaluation

OTR can be evaluated by means of the oxygen mass balance [27,28], according to Equation (1):
OTR = k L a · ( C O 2 * C O 2 ) ,
where CO2 is the oxygen concentration in the liquid phase (mol·L−1) and C O 2 * is the saturation oxygen concentration (mol·L−1).
For the experiments carried out at low agitation, it can be assumed that the concentration of dissolved oxygen (DO) will be very low or negligible due to the slower income of oxygen compared to the oxygen consumption by the culture. Therefore, the DHA production rate will be strongly related to OTR. Therefore, the DHA production rate can be described by Equation (2):
dC DHA dt = k L a · C O 2 * ,
where CDHA is DHA concentration (mol·L−1) and kLa is the volumetric mass transfer coefficient (h−1).
As predicted by Equation (2), a linear evolution of CDHA with time is expected. Thus, integrating Equation (2), assuming an initial DHA concentration = 0, Equation (3) is obtained:
C DHA = k L a · C O 2 * · t ,
where t is time (h).
This tendency can be observed in two different ways in the experiments carried out in this work:
  • For low agitation rate, the linear trend of DHA concentration with time adjusts to the tendency shown in Equation (3).
  • For the rest of experiments, only the first stage, at short times, DHA concentration evolves following a linear relationship with time.
Thus, Equation (2) can be used to describe the DHA production rate for the cases in which this linear tendency is observed throughout the whole process, which is evident in runs 1 to 6 (in shake flasks, at 100 rpm stirring speed) and 19 (in STBR, at 300 rpm stirring speed).
In Figure 6, the experimental results of these experiments are represented versus the results achieved using Equation (3) in both cases. The kinetic and statistical parameters values are given in Table 3.
As can be seen in Figure 7, the kinetic model (Equation (3)) is able to reproduce the experimental results of all the runs carried out at low agitation speeds in a suitable way. The results presented in Table 3 show a very small deviation of the kinetic parameter, kLa, and a value of r2 very close to 1. Furthermore, as can be seen in Figure 7a, the biomass concentration does not exert any influence on DHA production rate. Hence, it is possible to say that for this stirring speed, the mass transfer rate is controlling the overall process rate.
On the other hand, it is possible to estimate the value of kLa for the rest of experiments. However, it is necessary to employ only the initial linear stage, since the DHA concentration time evolves with time from a linear to a hyperbolic trend, as has been shown in Figure 2, Figure 3 and Figure 4 and Figure 6.
The results achieved for kLa calculation are presented in Table 4. With the results shown in this table, kLa values were fitted to agitation speed in order to correlate the stirring speed with the mass transfer rate. These correlations are given by Equations (4) and (5).
( k L a ) Flask = 10 4.75 · N 1.2 r 2   =   0.997
( k L a ) STBR = 10 4.90 · N 1.1   r 2   =   0.984
As can be seen in Table 4, all of the calculated kLa values are directly related to the agitation speed. Furthermore, as given by Equations (4) and (5), these relationships are similar in both systems, the effect of stirring speed depending on a function that is almost linear in both cases. While Equation (4) is a reasonable relationship [16,30], Equation (5) is not, which is surely due to the kLa values having been calculated under conditions where mass transfer is not the limiting step; thus, these values are not realistic [16]. Therefore, this method is useful to evaluate the impact of the mass transfer rate over the DHA production rate in a rapid and simple way.

3.2. Kinetic Model for High Agitation: Chemical Reaction Rate Calculation

As the reaction of the transformation of glycerol into DHA is carried out by a membrane-embedded enzyme [18], it can be described by a Michaelis–Menten-type kinetic model. Accordingly, DHA production model is stated by Equation (6):
dC DHA dt = k C · C X · C G K G + C G ,
where kC is the Michaelis–Menten (M-M) kinetic constant (L·gX−1·h−1), CG is glycerol concentration (mol·L−1); and KG is the M-M affinity constant (mol·L−1).
This kinetic model assumes that the enzymes on the membrane have enough oxygen available, which is an essential co-factor for the DHA production reaction. Therefore, this kinetic model only can be applied for high stirring speeds, and will be used only in runs 16 to 18 (agitation of 400 rpm in shake flasks) and runs 20 to 24 (stirring speed higher than 560 rpm in STBR).
The results of the kinetic model fitting to experimental results are shown in Figure 8, and the kinetic and statistical parameters are enclosed in Table 5. As can be appreciated in Figure 8, the model is able to reproduce the experimental results very well, being able to describe the results of different runs carried out at several stirring speeds and employing different biomass concentrations in shake flasks and STBR, respectively. On the other hand, for both scales, the kinetic model satisfies statistical criteria, with very low sum of residuals (SQR), standard error of estimates (SEE), and residual mean square error (RMSE) values, while the values of VE reaches very high values (more than 98%).

3.3. Kinetic Model for Intermediate Cases: Oxygen Concentration Effect

Although the kinetic models formulated for the previously commented fluid dynamic conditions are able to describe experiments conducted at both low and high stirring speeds, in the experimental program, there are experiments performed under intermediate conditions, and the results of these runs cannot be described using the previous kinetic models. In these experiments, the agitation and stirring speeds are neither high nor low. Consequently, in these runs, the OTR and chemical reaction rate will have a similar impact on the DHA production rate. Thus, it is necessary to propose another different kinetic model that is able to reflect the impact of oxygen availability on the DHA production rate.
For this purpose, it was mandatory to evaluate the DO concentration in the culture media during each run. Unfortunately, this was not possible for shake flasks experiments due to a lack of the required equipment. Thus, the next model was obtained only for the results corresponding to STBR experiments. The evolution of DO concentration is given in Figure 9, in which some interesting trends can be observed:
  • In the run with less agitation (300 rpm), DO concentration reduces to 0 shortly after the start of the experiments. This reinforces the idea that in this case, there is a shortage in oxygen availability caused by the low mass transfer rate.
  • In the experiments performed with more than 560 rpm, it can be seen that the DO concentration never reaches 0. Therefore, it is safe to say that the OTR is far greater than the DHA production rate (chemical reaction).
  • There are two runs, corresponding to those carried out at 360 and 420 rpm, in which the DO concentration reaches 0 in the first 3 h. However, in both cases, this oxygen shortage only lasts for a short time period.
  • Therefore, it is necessary to take into account the DO concentration in order to explain the DHA production rate. Thus, Equation (7) was formulated with these premises:
dC DHA dt = k C · C X · C G K G + C G · C O 2 n ,
where CO2 is DO concentration (mol·L−1) and n is an empirical exponent.
The results achieved by fitting Equation (7) to the experimental results of runs 20 and 21 are shown in Figure 10, and the kinetic and statistical parameters obtained are given in Table 6.
As can be seen in Figure 10, the experimental results of both runs are perfectly described by the model shown in Equation (7). This is also reflected by the values of the kinetic and statistical parameters, which can be seen in Table 6. This result highlights even more the importance of the OTR on DHA production, especially in the cases in which the oxygen availability in the culture is not oversaturated. Furthermore, it has been demonstrated that in order to describe and to maximize the DHA production rate, it is critical to choose OTR values that meet the oxygen requirements of the enzyme to perform the glycerol transformation at the maximum reaction rate.

4. Materials and Methods

4.1. Microorganism

The strain employed in this work is Gluconobacter oxydans ATCC 621, which was provided by the Spanish Type Culture Collection (CECT 360). Bacteria were stored in a glycerol solution stock at −80 °C. As the experiments were carried out at resting cell conditions, in which the bacteria cells are able to metabolize glycerol into DHA but can not replicate due to the lack of nitrogen source, it was necessary to grow the cells previously in two successive growth steps, employing the same methodology described in previous works [16,31].

4.2. Experimental Procedure

Two different experimental setups were employed in this work. In the first one, resting cells were transferred to 250 mL Erlenmeyer-type (shake flasks), with a working volume of 50 mL. The medium used in these experiments contains 20 mM of acetate buffer (initial pH of 5.5) and 25 g·L−1 (0.27 mol·L−1) of glycerol (supplied by Panreac AppliChem, Barcelona). The temperature was fixed at 30 °C, and the agitation was varied from 100 to 400 rpm. In the second experimental setup, a 1 L STBR with a working volume of 0.5 L was employed. The medium used in this setup was the same as the one previously described, as well as the temperature. The stirring speed was changed from 360 to 750 rpm, and the aeration was kept at 0.5 L·min-1 of the air flow rate. All of the experiments described in this work have been performed in triplicate, and only the mean value has been presented. The experimental error for all of the data was below 5%.

4.3. Analytical Methods

The concentration of biomass was measured employing spectrophotometry, as the relationship between the biomass concentration and the absorbance at 600 nm is linear. Both the equation employed to relate the biomass concentration to the absorbance and the methodology used to this purpose were stated elsewhere [15,16].
DHA, glycerol, and glyceric acid were determined by HPLC, with an Agilent 1100 series and a Rezex RHM Monosaccharide H+ column. Temperature was maintained at 80 °C, with an 0.5 mL·min−1 of elution rate using 1.2 mM H2SO4 acid water dilution as the mobile phase and 15 µL of injection volume. For DHA concentration determination, a diode array detector (DAD) was employed at 270 nm to measure DHA concentration, with a retention time of 16.8 min, and a refractive index detector (RID) was employed to measure the other compounds, with a retention time of 17.6 and 15.3 min for glycerol and glyceraldehyde, respectively. Commercial glycerol, DHA, and glyceric acid were employed for calibration purposes (supplied by Panreac AppliChem, Barcelona).

4.4. Mathematics and Calculations

Productivity as a specific rate of DHA bioproduction (PX) was evaluated for each experiment employing Equation (8):
P X = C DHA max C X · t max ,
where C DHA max (mol·L−1) is the maximum DHA concentration for each experiment, and tmax (h) is the time in which this concentration is reached.
To calculate kLa in the runs in which the mass transfer rate exerts control over the global process rate, as well as in the applicable zones in which this control is not maintained throughout the whole reaction time, OriginLab software was employed, using a linear fitting algorithm. This software was also used for the correlation studies between kLa and stirring speed.
To fit the kinetic models to experimental data, the Aspen Custom Modeler© software was employed. For the purpose of validation of the different models and in order to quantify the goodness of the fit, several statistical parameters were employed, such as the following referring to DHA concentration: Fischer’s F criterion (F), sum of residuals (SQR), standard error of estimates (SEE), residual mean square error (RMSE), and percentage of explained variance (VE). All of these statistical parameters are defined in Appendix A.

5. Conclusions

In summary, there is an important effect of OTR on DHA production rate employing G. oxydans as the biocatalyst. For low or very low mass transfer rates, the DHA production rate is controlled by the mass transfer rate, limiting the productivity of DHA in each experiment. However, once the maximum mass transfer rate is achieved, an increase in this rate by increasing agitation does not affect the DHA production rate, which is affected by other different factors related to the chemical reaction rate, such as biomass concentration and product/substrate inhibition, as previously reported in the literature [7,14,18,22,23]. These achieved results are in line with the previously published results regarding this system [16].

Author Contributions

Conceptualization: V.E.S. and F.G.-O. Data curation: S.d.l.M. and M.W. Formal analysis: S.d.l.M. and M.W. Funding acquisition: V.E.S. and F.G.-O. Investigation: S.d.l.M. Methodology: S.d.l.M., M.W., V.E.S. and F.G.-O. Project administration: F.G.-O. Resources: V.E.S. and F.G.-O.; Software: S.d.l.M. and M.W.; Supervision: V.E.S. and F.G.-O. Validation: S.d.l.M. and M.W. Visualization: S.d.l.M. and M.W. Writing—original draft: M.W. Writing—review and editing: S.d.l.M., M.W., V.E.S. and F.G.-O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “Ministerio de Economía, Industria y Competitividad”, grant number CTQ2017-84963-C2-1-R and CTQ2010-15460 and the grant FPI BES-2011-044438.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Fischer’s F criterion refers to the adequacy of the model to the observed values of the chosen variable, and it is based on the null hypothesis, and expressed by Equation (A1).
F = i = 1 N ( y i , calc ) 2 P i = 1 N SQR N P ,
where yi, calc is the observed variable value predicted by the model, N is the number of data, P is the parameter number, and SQR is defined by Equation (A2).
SQR = i = 1 N ( y i , exp y i , cal ) 2 ,
where yi,exp is the experimental value for the observed variable and yi,cal is the value predicted by the model.
SEE is calculated by Equation (A3).
SEE = i = 1 N ( y i , exp y i , calc ) 2 N 2 ,
For RMSE calculation, Equation (A4) was employed.
RMSE = SQR N P ,

References

  1. Veluturla, S.; Archna, N.; Subba Rao, D.; Hezil, N.; Indraja, I.S.; Spoorthi, S. Catalytic valorization of raw glycerol derived from biodiesel: A review. Biofuels 2018, 9, 305–314. [Google Scholar] [CrossRef]
  2. Monteiro, M.R.; Kugelmeier, C.L.; Pinheiro, R.S.; Batalha, M.O.; da Silva César, A. Glycerol from biodiesel production: Technological paths for sustainability. Renew. Sustain. Energy Rev. 2018, 88, 109–122. [Google Scholar] [CrossRef]
  3. Nda-Umar, U.; Ramli, I.; Taufiq-Yap, Y.; Muhamad, E. An overview of recent research in the conversion of glycerol into biofuels, fuel additives and other bio-based chemicals. Catalysts 2019, 9, 15. [Google Scholar] [CrossRef] [Green Version]
  4. Mitrea, L.; Trif, M.; Cătoi, A.-F.; Vodnar, D.-C. Utilization of biodiesel derived-glycerol for 1,3-PD and citric acid production. Microb. Cell Factories 2017, 16, 190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Hekmat, D.; Bauer, R.; Fricke, J. Optimization of the microbial synthesis of dihydroxyacetone from glycerol with Gluconobacter oxydans. Bioprocess Biosyst. Eng. 2003, 26, 109–116. [Google Scholar] [CrossRef] [PubMed]
  6. Jain, V.K.; Tear, C.J.Y.; Lim, C.Y. Dihydroxyacetone production in an engineered Escherichia coli through expression of Corynebacterium glutamicum dihydroxyacetone phosphate dephosphorylase. Enzym. Microb. Technol. 2016, 86, 39–44. [Google Scholar] [CrossRef]
  7. Lin, X.; Liu, S.; Xie, G.; Chen, J.; Li, P.; Chen, J. Enhancement of 1,3-dihydroxyacetone production from Gluconobacter oxydans by combined mutagenesis. J. Microbiol. Biotechnol. 2016, 26, 1908–1917. [Google Scholar] [CrossRef] [Green Version]
  8. Zheng, X.J.; Jin, K.Q.; Zhang, L.; Wang, G.; Liu, Y.P. Effects of oxygen transfer coefficient on dihydroxyacetone production from crude glycerol. Braz. J. Microbiol. 2016, 47, 129–135. [Google Scholar] [CrossRef] [Green Version]
  9. Zhou, X.; Zhou, X.; Xu, Y.; Yu, S. Improving the production yield and productivity of 1,3-dihydroxyacetone from glycerol fermentation using Gluconobacter oxydans NL71 in a compressed oxygen supply-sealed and stirred tank reactor (COS-SSTR). Bioprocess Biosyst. Eng. 2016, 39, 1315–1318. [Google Scholar] [CrossRef]
  10. Hu, Z.-C.; Tian, S.-Y.; Ruan, L.-J.; Zheng, Y.-G. Repeated biotransformation of glycerol to 1,3-dihydroxyacetone by immobilized cells of Gluconobacter oxydans with glycerol- and urea-feeding strategy in a bubble column bioreactor. Bioresour. Technol. 2017, 233, 144–149. [Google Scholar] [CrossRef]
  11. Poljungreed, I.; Boonyarattanakalin, S. Dihydroxyacetone production by Gluconobacter frateurii in a minimum medium using fed-batch fermentation. J. Chem. Technol. Biotechnol. 2017, 92, 2635–2641. [Google Scholar] [CrossRef]
  12. Dikshit, P.K.; Moholkar, V.S. Batch and repeated-batch fermentation for 1,3-Dihydroxyacetone production from waste glycerol using free, immobilized and resting Gluconobacter oxydans cells. Waste Biomass Valorizat. 2018, 10, 2455–2465. [Google Scholar] [CrossRef]
  13. Stasiak-Różańska, L.; Berthold-Pluta, A.; Dikshit, P.K. Valorization of waste glycerol to dihydroxyacetone with biocatalysts obtained from Gluconobacter oxydans. Appl. Sci. 2018, 8, 2517. [Google Scholar] [CrossRef] [Green Version]
  14. Sudarshan, B.; Sanjay, K. Optimization of fermentation conditions for production of 1,3-Dihydroxyacetone from glycerol obtained as a byproduct during biodiesel production. Microbiol. Res. J. Int. 2018, 1–9. [Google Scholar] [CrossRef]
  15. de la Morena, S.; Acedos, M.G.; Santos, V.E.; García-Ochoa, F. Dihydroxyacetone production from glycerol using Gluconobacter oxydans: Study of medium composition and operational conditions in shaken flasks. Biotechnol. Prog. 2019, 35, e2803. [Google Scholar] [CrossRef]
  16. de la Morena, S.; Santos, V.E.; García-Ochoa, F. Influence of oxygen transfer and uptake rates on dihydroxyacetone production from glycerol by Gluconobacter oxydans in resting cells operation. Biochem. Eng. J. 2019, 147, 20–28. [Google Scholar] [CrossRef]
  17. Keliang, G.; Dongzhi, W. Asymmetric oxidation by Gluconobacter oxydans. Appl. Microbiol. Biotechnol. 2006, 70, 135–139. [Google Scholar] [CrossRef]
  18. Deppenmeier, U.; Hoffmeister, M.; Prust, C. Biochemistry and biotechnological applications of Gluconobacter strains. Appl. Microbiol. Biotechnol. 2002, 60, 233–242. [Google Scholar]
  19. Rocha-Martin, J.; Acosta, A.; Berenguer, J.; Guisan, J.M.; Lopez-Gallego, F. Selective oxidation of glycerol to 1,3-dihydroxyacetone by covalently immobilized glycerol dehydrogenases with higher stability and lower product inhibition. Bioresour. Technol. 2014, 170, 445–453. [Google Scholar] [CrossRef]
  20. Xu, M.-Q.; Li, F.-L.; Yu, W.-Q.; Li, R.-F.; Zhang, Y.-W. Combined cross-linked enzyme aggregates of glycerol dehydrogenase and NADH oxidase for high efficiency in situ NAD+ regeneration. Int. J. Biol. Macromol. 2019. [Google Scholar] [CrossRef]
  21. Ma, L.; Lu, W.; Xia, Z.; Wen, J. Enhancement of dihydroxyacetone production by a mutant of Gluconobacter oxydans. Biochem. Eng. J. 2010, 49, 61–67. [Google Scholar] [CrossRef]
  22. Stasiak-Różańska, L.; Błażejak, S.; Gientka, I.; Bzducha-Wróbel, A.; Lipińska, E. Utilization of a waste glycerol fraction using immobilized cell extract from Gluconobacter oxydans ATCC 621 and its possible re-use. Electron. J. Biotechnol. 2017, 27, 44–48. [Google Scholar] [CrossRef]
  23. Dikshit, P.K.; Moholkar, V.S. Kinetic analysis of dihydroxyacetone production from crude glycerol by immobilized cells of Gluconobacter oxydans MTCC 904. Bioresour. Technol. 2016, 216, 948–957. [Google Scholar] [CrossRef] [PubMed]
  24. Dikshit, P.K.; Padhi, S.K.; Moholkar, V.S. Process optimization and analysis of product inhibition kinetics of crude glycerol fermentation for 1,3-Dihydroxyacetone production. Bioresour. Technol. 2017, 244, 362–370. [Google Scholar] [CrossRef]
  25. Garcia-Ochoa, F.; Gomez, E. Bioreactor scale-up and oxygen transfer rate in microbial processes: An overview. Biotechnol. Adv. 2009, 27, 153–176. [Google Scholar] [CrossRef]
  26. Rodriguez, A.; Escobar, S.; Gomez, E.; Santos, V.E.; Garcia-Ochoa, F. Behavior of several Pseudomonas putida strains growth under different agitation and oxygen supply conditions. Biotechnol. Prog. 2018, 34, 900–909. [Google Scholar] [CrossRef]
  27. Acedos, M.G.; Hermida, A.; Gomez, E.; Santos, V.E.; Garcia-Ochoa, F. Effects of fluid-dynamic conditions in Shimwellia blattae (p424IbPSO) cultures in stirred tank bioreactors: Hydrodynamic stress and change of metabolic routes by oxygen availability. Biochem. Eng. J. 2019, 149, 107238. [Google Scholar] [CrossRef]
  28. Garcia-Ochoa, F.; Gomez, E.; Santos, V.E.; Merchuk, J.C. Oxygen uptake rate in microbial processes: An overview. Biochem. Eng. J. 2010, 49, 289–307. [Google Scholar] [CrossRef]
  29. Jackson, E.; Ripoll, M.; Betancor, L. Efficient glycerol transformation by resting Gluconobacter cells. Microbiol. Open 2019, 8, e926. [Google Scholar] [CrossRef] [Green Version]
  30. Klöckner, W.; Büchs, J. Advances in shaking technologies. Trends Biotechnol. 2012, 30, 307–314. [Google Scholar] [CrossRef]
  31. de la Torre Pascual, I.; Acedos, M.G.; Ladero, M.; Santos Mazorra, V. On the use of resting L. delbrueckii spp. delbrueckii cells for D-lactic acid production from orange peel wastes hydrolysates. Biochem. Eng. J. 2019, 145, 162–169. [Google Scholar] [CrossRef]
Figure 1. Experimental results. 1,3-dihydroxyacetone (DHA) (●) and glycerol (○) concentrations vs. time for runs 1 to 6 (see Table 1) carried out in shake flasks. All runs conducted at agitation speed of 100 rpm employing different biomass concentrations according to (a) 0.1, (b) 0.2, (c) 0.3, (d) 0.4, (e) 0.5, and (f) 1.0 g·L−1.
Figure 1. Experimental results. 1,3-dihydroxyacetone (DHA) (●) and glycerol (○) concentrations vs. time for runs 1 to 6 (see Table 1) carried out in shake flasks. All runs conducted at agitation speed of 100 rpm employing different biomass concentrations according to (a) 0.1, (b) 0.2, (c) 0.3, (d) 0.4, (e) 0.5, and (f) 1.0 g·L−1.
Catalysts 10 00101 g001
Figure 2. Experimental results. DHA (●) and glycerol (○) concentrations vs. time for runs 7 to 11 (see Table 1) carried out in shake flasks. All runs conducted at agitation speed of 200 rpm employing different biomass concentrations according to (a) 0.2, (b) 0.3, (c) 0.4, (d) 0.5, and (e) 1.0 g·L−1.
Figure 2. Experimental results. DHA (●) and glycerol (○) concentrations vs. time for runs 7 to 11 (see Table 1) carried out in shake flasks. All runs conducted at agitation speed of 200 rpm employing different biomass concentrations according to (a) 0.2, (b) 0.3, (c) 0.4, (d) 0.5, and (e) 1.0 g·L−1.
Catalysts 10 00101 g002
Figure 3. Experimental results. DHA (●) and glycerol (○) concentrations vs. time for runs 12 to 16 (see Table 1) carried out in shake flasks. All runs conducted at agitation speed of 300 rpm employing different biomass concentrations according to (a) 0.3, (b) 0.4, (c) 0.5, and (d) 1.0 g·L−1.
Figure 3. Experimental results. DHA (●) and glycerol (○) concentrations vs. time for runs 12 to 16 (see Table 1) carried out in shake flasks. All runs conducted at agitation speed of 300 rpm employing different biomass concentrations according to (a) 0.3, (b) 0.4, (c) 0.5, and (d) 1.0 g·L−1.
Catalysts 10 00101 g003
Figure 4. Experimental results. DHA (●) and glycerol (○) concentrations vs. time for runs 17 to 19 (see Table 1) carried out in shake flasks. All runs conducted at agitation speed of 300 rpm employing different biomass concentrations according to (a) 0.4, (b) 0.5, and (c) 1.0 g·L−1.
Figure 4. Experimental results. DHA (●) and glycerol (○) concentrations vs. time for runs 17 to 19 (see Table 1) carried out in shake flasks. All runs conducted at agitation speed of 300 rpm employing different biomass concentrations according to (a) 0.4, (b) 0.5, and (c) 1.0 g·L−1.
Catalysts 10 00101 g004
Figure 5. Productivity results corresponding to each biomass concentration employed for shake flask runs: 100 rpm (●), 200 rpm (○), 300 rpm (■), and 400 rpm (□).
Figure 5. Productivity results corresponding to each biomass concentration employed for shake flask runs: 100 rpm (●), 200 rpm (○), 300 rpm (■), and 400 rpm (□).
Catalysts 10 00101 g005
Figure 6. Experimental results. DHA (●) and glycerol (○) concentrations vs. time for runs 20 to 24 (see Table 1 carried out in a stirred tank bioreactor (STBR). All runs conducted with a biomass concentration of 0.5 g·L−1, employing several stirring speeds according to: (a) 300, (b) 360, (c) 420, (d) 560, (e) 650, and (f) 750 rpm.
Figure 6. Experimental results. DHA (●) and glycerol (○) concentrations vs. time for runs 20 to 24 (see Table 1 carried out in a stirred tank bioreactor (STBR). All runs conducted with a biomass concentration of 0.5 g·L−1, employing several stirring speeds according to: (a) 300, (b) 360, (c) 420, (d) 560, (e) 650, and (f) 750 rpm.
Catalysts 10 00101 g006
Figure 7. Reproduction of DHA concentration vs. time for (a) runs 1 to 6 (carried out in shake flasks) and (b) for run 19 (carried out in STBR), the agitation being (a) 100 rpm (flasks) and (b) 300 rpm (STBR). Different biomass concentrations were employed according to the following legends: 0.1 (●), 0.2 (○), 0.3 (■), 0.4 (□), 0.5 (▲), and 1.0 g·L−1 (∆). Points corresponding to experimental results, lines corresponding to kinetic model simulation when the oxygen transport rate (OTR) is controlling the overall process rate. Model Equations (2) and (3).
Figure 7. Reproduction of DHA concentration vs. time for (a) runs 1 to 6 (carried out in shake flasks) and (b) for run 19 (carried out in STBR), the agitation being (a) 100 rpm (flasks) and (b) 300 rpm (STBR). Different biomass concentrations were employed according to the following legends: 0.1 (●), 0.2 (○), 0.3 (■), 0.4 (□), 0.5 (▲), and 1.0 g·L−1 (∆). Points corresponding to experimental results, lines corresponding to kinetic model simulation when the oxygen transport rate (OTR) is controlling the overall process rate. Model Equations (2) and (3).
Catalysts 10 00101 g007
Figure 8. Reproduction of DHA concentration vs. time for runs 16 to 18 (carried out in shake flasks, (upper graphic) and for runs 22 to 24 (carried out in STBR, lower graphic). Runs in shake flasks conducted at 400 rpm of agitation and several biomass concentrations: 0.4 (●), 0.5 (○), and 1.0 g·L−1 (■). Runs carried out in STBR with a biomass concentrations of 0.5 g·L−1 and several stirring speeds: 560 rpm (□), 650 rpm (▲), and 750 rpm (∆). Points correspond to experimental results, and lines correspond to kinetic model simulation when the oxygen transport rate (OTR) is controlling the overall process rate. Model Equation (6).
Figure 8. Reproduction of DHA concentration vs. time for runs 16 to 18 (carried out in shake flasks, (upper graphic) and for runs 22 to 24 (carried out in STBR, lower graphic). Runs in shake flasks conducted at 400 rpm of agitation and several biomass concentrations: 0.4 (●), 0.5 (○), and 1.0 g·L−1 (■). Runs carried out in STBR with a biomass concentrations of 0.5 g·L−1 and several stirring speeds: 560 rpm (□), 650 rpm (▲), and 750 rpm (∆). Points correspond to experimental results, and lines correspond to kinetic model simulation when the oxygen transport rate (OTR) is controlling the overall process rate. Model Equation (6).
Catalysts 10 00101 g008
Figure 9. Experimental variation of dissolved oxygen (DO) concentration vs. time for runs 19 to 24 carried out in STBR employing different stirring speeds: 300 (●), 360 (○) 420 (■), 560 (□), 650 (▲), and 750 rpm (∆).
Figure 9. Experimental variation of dissolved oxygen (DO) concentration vs. time for runs 19 to 24 carried out in STBR employing different stirring speeds: 300 (●), 360 (○) 420 (■), 560 (□), 650 (▲), and 750 rpm (∆).
Catalysts 10 00101 g009
Figure 10. Variation of DHA concentration vs. time for runs 20 and 21, carried out in STBR, with stirring speeds of 360 rpm (●) and 420 rpm (○). Points are experimental results and lines plot the kinetic model given by Equation (7).
Figure 10. Variation of DHA concentration vs. time for runs 20 and 21, carried out in STBR, with stirring speeds of 360 rpm (●) and 420 rpm (○). Points are experimental results and lines plot the kinetic model given by Equation (7).
Catalysts 10 00101 g010
Table 1. Experimental planning: runs carried out under different operating conditions (biomass concentration and agitation) employing both shake flasks (runs 1 to 18) and stirred tank bioreactor (STBR) (runs 19 to 24).
Table 1. Experimental planning: runs carried out under different operating conditions (biomass concentration and agitation) employing both shake flasks (runs 1 to 18) and stirred tank bioreactor (STBR) (runs 19 to 24).
RunCX 1 (g·L−1)N 2 (rpm)Experimental Device
10.1100Shake flasks
20.2
30.3
40.4
50.5
61.0
70.2200
80.3
90.4
100.5
111.0
120.3300
130.4
140.5
151.0
160.4400
170.5
181.0
190.5300STBR (Stirred tank bioreactor)
20360
21420
22560
23650
Table 2. Productivity results achieved in STBR runs.
Table 2. Productivity results achieved in STBR runs.
NPX (mol·gx−1·L−1)·102
3001.65
3602.03
4203.11
5603.34
6503.46
7503.73
Table 3. Kinetic and statistical parameters calculated for low agitation experiments with mass transfer control of the overall process rate. Runs 1 to 6 (shake flasks) and 19 (STBR).
Table 3. Kinetic and statistical parameters calculated for low agitation experiments with mass transfer control of the overall process rate. Runs 1 to 6 (shake flasks) and 19 (STBR).
RunCX (g·L−1) C O 2 *   ( mol · L 1 ) · 10 4 kLa·103·(h−1)r2
10.14.756.92 ± 0.050.997
20.2
30.3
40.4
50.5
61.0
190.54.758.21 ± 0.050.9997
Table 4. Mass transfer coefficient values calculated for runs 7 to 18 (shake flasks) and runs 20 to 24 (STBR); all runs conducted at low agitation.
Table 4. Mass transfer coefficient values calculated for runs 7 to 18 (shake flasks) and runs 20 to 24 (STBR); all runs conducted at low agitation.
RunkLa·103·(h−1)N (rpm)RunkLa·103·(h−1)N (rpm)RunkLa·103·(h−1)N (rpm)
78.2 ± 0.42001314.6 ± 0.73002010.0 ± 0.7360
8142112.5 ± 0.3420
9152217.7 ± 0.5560
101619.9 ± 0.54002319.6 ± 0.6650
11172423.6 ± 0.5750
1214.6 ± 0.730018
Table 5. Kinetic and statistical parameters calculated for high agitation experiments with the reaction rate control of the overall process rate. Runs 16 to 18 (shake flasks) and 22 to 24 (STBR). Model Equation (6). RMSE: residual mean square error, SQR: sum of residuals, SEE: standard error of estimates, VE: percentage of explained variance.
Table 5. Kinetic and statistical parameters calculated for high agitation experiments with the reaction rate control of the overall process rate. Runs 16 to 18 (shake flasks) and 22 to 24 (STBR). Model Equation (6). RMSE: residual mean square error, SQR: sum of residuals, SEE: standard error of estimates, VE: percentage of explained variance.
DeviceParameterValue
Shake flaskskC (L·gx−1·h−1)·101.67 ± 0.71
KG (mol·L−1)·101.94 ± 1.60
SQR 1·1032.07
SEE 2·1039.50
RMSE 3·1024.55
Fcalc/Ftab81.6/2.5
VE 4 (%)98.5
STBRkC (L·gx−1·h−1)·1029.50 ± 2.49
KG (mol·L−1)·101.55 ± 0.82
SQR1·1032.85
SEE 2·1021.00
RMSE 3·1039.74
Fcalc/Ftab3990/2.5
VE 4 (%)98.7
1 Sum of Quadratic Residue; 2 Standard Error of Estimate, 3 Residual Mean Square Error, 4 Variation explained.
Table 6. Kinetic and statistical parameters calculated for medium agitation experiments with a mix control of the mass transfer and the reaction rate. Runs 20 and 21 conducted in STBR. Model Equation (7).
Table 6. Kinetic and statistical parameters calculated for medium agitation experiments with a mix control of the mass transfer and the reaction rate. Runs 20 and 21 conducted in STBR. Model Equation (7).
ParameterValue
kC (L·gx−1·h−1)·105.2 ± 1.3
KG (mol·L−1)·105.6 ± 2.2
n0.14 ± 0.02
SQR·1036.26
SEE·1022.40
RMSE·1022.11
Fcalc/Ftab187/4.3
VE (%)97.2

Share and Cite

MDPI and ACS Style

de la Morena, S.; Wojtusik, M.; Santos, V.E.; Garcia-Ochoa, F. Kinetic Modeling of Dihydroxyacetone Production from Glycerol by Gluconobacter oxydans ATCC 621 Resting Cells: Effect of Fluid Dynamics Conditions. Catalysts 2020, 10, 101. https://doi.org/10.3390/catal10010101

AMA Style

de la Morena S, Wojtusik M, Santos VE, Garcia-Ochoa F. Kinetic Modeling of Dihydroxyacetone Production from Glycerol by Gluconobacter oxydans ATCC 621 Resting Cells: Effect of Fluid Dynamics Conditions. Catalysts. 2020; 10(1):101. https://doi.org/10.3390/catal10010101

Chicago/Turabian Style

de la Morena, Susana, Mateusz Wojtusik, Victoria E. Santos, and Felix Garcia-Ochoa. 2020. "Kinetic Modeling of Dihydroxyacetone Production from Glycerol by Gluconobacter oxydans ATCC 621 Resting Cells: Effect of Fluid Dynamics Conditions" Catalysts 10, no. 1: 101. https://doi.org/10.3390/catal10010101

APA Style

de la Morena, S., Wojtusik, M., Santos, V. E., & Garcia-Ochoa, F. (2020). Kinetic Modeling of Dihydroxyacetone Production from Glycerol by Gluconobacter oxydans ATCC 621 Resting Cells: Effect of Fluid Dynamics Conditions. Catalysts, 10(1), 101. https://doi.org/10.3390/catal10010101

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop