Research Article  Open Access
Statistical Optimization of Process Parameters for LipaseCatalyzed Synthesis of TriethanolamineBased Esterquats Using Response Surface Methodology in 2Liter Bioreactor
Abstract
Lipasecatalyzed production of triethanolaminebased esterquat by esterification of oleic acid (OA) with triethanolamine (TEA) in nhexane was performed in 2 L stirredtank reactor. A set of experiments was designed by central composite design to process modeling and statistically evaluate the findings. Five independent process variables, including enzyme amount, reaction time, reaction temperature, substrates molar ratio of OA to TEA, and agitation speed, were studied under the given conditions designed by Design Expert software. Experimental data were examined for normality test before data processing stage and skewness and kurtosis indices were determined. The mathematical model developed was found to be adequate and statistically accurate to predict the optimum conversion of product. Response surface methodology with central composite design gave the best performance in this study, and the methodology as a whole has been proven to be adequate for the design and optimization of the enzymatic process.
1. Introduction
Triethanolamine (TEA) based esterquat has been the primary ingredient in European fabric softeners and is becoming the global molecule of choice for various industries [1]. Esterquat cationic surfactant was considered as a type of biodegradable material utilized as a textile softening agent. There have been an increasing number of researchers who concern the biodegradable esterquat cationic surfactant since the beginning of the 1990s [2]. In addition, they are highly biodegradable and biocompatible because their ester bonds are easily hydrolyzed [3–5]. Besides biodegradability additional advantages such as excellent softening properties, suitability for various fabrics, and simple preparation procedures have been discovered by the use of esterquat cationic surfactants as textile softening agents [2].
In this work, triethanolamine and oleic acid were chosen as substrates to design an optimal model reaction which will lead to high conversion rate utilizing lipase from Candida antarctica (Novozym 435) as a biocatalyst in the organic solvent system. The investigated reaction conditions included enzyme amount, reaction time, reaction temperature, the molar ratio of substrates, and agitation speed. The major aim of this study was to model the effect of process parameters on the reaction yield. The most important stages in a process were modeling and optimization to improve a system and increase the efficiency of the process without increasing the cost. All process parameters are selected to conduct the optimization by using response surface methodology (RSM) [6–10]. The optimization of process has been reported by artificial neural network (ANN) in our previous study [11]. Subsequently, the simulated result in optimum conditions from the response surface methodology (RSM) and ANN was compared.
Prior to doing the statistical analyses by RSM, experimental data were inspected and explored the nature of variables by several normality tests. It has been observed that only a few researchers have paid attention to the application of right and accurate statistical techniques in order to validate experimental data [12]. Statistical methods are based on various underlying assumptions. One common assumption is that a random variable is normally distributed. In many statistical analyses, normality is often conveniently assumed without any empirical evidence or test. However, normality is critical in many statistical methods. Testing of assumptions usually involves obtaining descriptive statistics on variables [13]. Descriptive statistics provide important information about variables to be analyzed. Mean, median, and mode measure central tendency of a variable. Measures of dispersion include variance, standard deviation, and range. Researchers may draw a histogram, stemandleaf plot, or box plot to see how a variable is distributed [14]. When this assumption is violated, interpretation and inference may not be reliable or valid [15]. The usual processes in the statistical assessment of a data set are as follows: screen the data for outliers or blunders, plot the data to detect asymmetry and tail weight, calculate the indices of sample shape (i.e., skewness and kurtosis), perform tests of normality, and if the data is normal use parametric statistics for further analysis [15]. In order to test the validity of a normal distribution, quantitative tests need to be employed, such as KolmogorovSmirnov, Liliefors, and ShapiroWilks. In this study, normality tests also included the KolmogorovSmirnov (Lilliefors modification) and the ShapiroWilk for checking the normal distribution validity of variables.
2. Methods
2.1. Materials
Novozym 435, Candida antarctica lipase B immobilized on a macroporous acrylic resin (10,000 propyl laurate units per gram), was purchased from Novo Nordisk A/S (Bagsværd, Denmark). The enzyme is a granular product with a particle size of 0.2–0.6 mm. The bulk density of Novozym 435 is 350–450 kg/m^{3}. nHexane obtained from J. T. Baker (USA) was used as the organic solvent. Oleic acid and triethanolamine were purchased from Merck, Germany. All other chemicals used in this study were of analytical reagent grade.
2.2. Experimental Design
The optimization study was carried out in accordance with the experimental design with 5 factors and 5 levels with 50 experimental points. The fractional factorial designs consisted of 32 factorial points, 10 axial points (two axial points on the axis of each design variable at a distance of 1.75 from the design center) and eight center points. The generalized response surface model is shown by (1), and the variables and their levels selected for the study were represented in Table 1: where (conversion %) represents the response variable, is the constant term, represents the coefficients of the linear parameters, represents the variables, represents the coefficients of the quadratic parameter, represents the coefficients of the interaction parameters, and is the residual associated to the experiments.

2.3. Enzymatic Esterification and Analysis of Samples
The reactions were performed in 2000 mL reactor, and specified volumes of hexane were added as solvent. The reactor consisted of a screw cap and a glass flask with a capacity of 2 liters and an inner diameter of 10 cm. A fourbladed impeller (4.5 cm in diameter) was immersed in the reaction mixture a 2 cm height from the bottom of the flask to provide agitation effect. The impeller was connected by a shaft to motor for speed controlling purpose. A baffle was connected to the cap and immersed in the reaction mixture. The reaction temperature was controlled by immersing reactor in a temperaturecontrolled water bath. The reactions were catalyzed by various amounts of Novozym 435 from 1.5 to 8.5% w/w of oleic acid for experimental design at different temperature (51.25–68.75°C) and agitation speed (137.5–662.5 r.p.m.) values. The studied ranges of the substrates were 708 mmol for OA as a constant amount, while concentrations of TEA were varied according to Table 1 for the experimental design. All experiments were carried out in the range of 2–30 h, as shown in Table 1. The basic points for the design were selected from a preliminary study in laboratory scale [16] by using Taguchi design (data not shown).
At the end of the reaction periods, 30 mL aliquot was withdrawn from the system using a syringe. The reaction sample was terminated by dilution with 10 mL of ethanolacetone (50 : 50, v/v). The enzyme particles were then separated by filtration, and the remaining free acid in the reaction mixture was determined by titration of the aliquots of reaction mixture against standard NaOH. The amount of reacted acid was determined from the values obtained for the control (without enzyme) and test samples. The ester formed was expressed as equivalent to conversion of the acid [17]. The ester formation was confirmed by thinlayer chromatography (TLC) using chloroform : methanol (95 : 5) solvent system. Further identification for ester formation was carried out by FTIR (Perkin Elmer, model 1650) and gas chromatography/mass spectroscopy GCMS on a Shimadzu (model GC 17A; model MS QP5050A, Tokyo, Japan) instruments. For purification of the product, after termination of the reaction, the enzyme was filtered and the solvent removed by evaporator under reduced pressure. Product in the remaining mixture was separated via silica gel (Kieselgel 60, Merck, particle size 0.063–0.200 mm) column chromatography (15 cm% 20 mm) using a dichloromethane/methanol (90/10, v/v) mixture as eluent. A sample made up of 1 : 1 (w/w) ratio of silica gel and the free solvent reaction mixture was deposited at the top of the column previously equilibrated with dichloromethane/methanol (90/10, v/v) mixture. Five milliliter fractions were collected and tested using thin layer chromatography to identify the productrich portions. Such fractions were pooled and the solvent evaporated by a rotary evaporator. The purity of product was then checked with TLC and GC before FTIR analysis.
3. Results and Discussion
3.1. Testing Experimental Data for Normality
Normal is used to describe a symmetrical, bellshaped curve, which has the greatest frequency of scores in the middle, with smaller frequencies towards the extremes [18]. Normality can be assessed to some extent by obtaining skewness and kurtosis values.
Table 2 shows descriptive statistics to check the skewness and kurtosis values for five variables at three levels of each of them. For other levels conversions percentage was constant when variables were placed in ±1.75 levels, and they have been omitted.

The results showed that skewness ranged between −0.925 and 0.532 (acceptable range of normality is between −2.0 and +2.0). The values of kurtosis ranged between −0.848 and 1.111 (acceptable range of normality is between −5.0 and +5.0) [19]. As a result, the skewness and kurtosis values indicate almost normal distribution. However, these descriptive statistics do not provide conclusive information about normality, and testing normality needs to use some other statistics tests. SPSS software provides two different statistics for testing normality. The ShapiroWilk and KolmogorovSmirnov tests were used for data distribution analysis. Both tests similarly demonstrated that the data set was normally distributed. As shown in Table 3, the values of ShapiroWilk and KolmogorovSmirnov tests confirm null hypothesis that the variable are normally distributed (). Since the number of observations is less than 2,000, however, ShapiroWilk test will be appropriate to this case.

3.2. Data Processing and Analysis of Variance (ANOVA)
The results at each point based on experimental design for the enzymatic reaction of TEAbased esterquat are presented in Table 4. Evaluation of coefficients of the empirical models and their statistical analyses were carried out using central composite design.

Fitting of the data to various models (linear, 2FI, quadratic, and cubic) and their subsequent analysis of variance showed that TEAbased esterquat synthesis was most suitably described with a quadratic model. The model was modified based on the insignificancy of some model terms. The final reduced model to predict the conversion % of TEAbased esterquat catalyzed by Novozym 435 is shown as follows: where matches product conversion % and , , , , and match to coded values for the enzyme amount (% w/w), reaction time (h), reaction temperature (°C), the molar ratio of substrates (mole), and agitation speed (r.p.m.), respectively. The positive sign in front of the terms indicates a synergistic effect while the negative sign indicates an antagonistic effect. Negative values of coefficient estimates denote negative influence of parameters on the reaction. It was observed that all the linear coefficients from the model gave positive effect except the coefficient estimate for the molar ratio of substrates () in the model of percentage conversion. This may be due to that the percentage of conversion was negatively affected by the presence of the higher ratio of oleic acid as the ratio of oleic acid/triethanolamine. From the equation, the conversion of enzymatic reaction has linear and quadratic effects by the five process variables. The model was found to have coefficient of determination value () of 0.9201, which means that 92.01% of the total variation in the results was attributed to the independent variables investigated. When approaches unity, the better empirical model fits the actual data [20]. Normally, a regression model having an value higher than 0.9 was considered as model having a very high correlation [21]. Hence, the value in this regression model is relatively high, which indicates a good agreement between predicted and experimental conversion of TEAbased esterquat reaction. Figure 1 summarizes correlation between experimental values and predicted values by using the developed model.
(a)
(b)
(c)
Figure 1(a) shows the actual values versus predicted values of the product conversion %, which indicated a good agreement between actual and predicted responses. A residual plot allowed visual assessment of the distance of each observation from the fitted line (Figure 1(b)). The residuals randomly scattered in a constant width band about the zero line. Figure 1(c) shows the histogram of the residuals in allowed visual assessment of the assumption. As observed, the measurement errors in the response variable were normally distributed, and the histogram of the residuals revealed a normal distribution overlay.
Statistical analysis based on ANOVA for the response surface quadratic model is presented in Table 5. The value for the model is less than 0.05, which indicates that it is a significant and desirable model. Besides, the value of indicates that there is only a 0.01% chance that a “model value” this large could occur due to noise in the experiments. The “Lack of Fit value” of 0.43 implies that lack of fit is not significant relative to pure error. Thus, it is possible quantitatively judge if the model represents the observations satisfactorily.

3.3. Mutual Effect of Process Parameters
The terms in (2) show that interactions between variables have significant effect on the conversion% of enzymatic reaction of TEAbased esterquat. Therefore, instead of studying single variable the interactions will be investigated, which is significant and important for a comprehensive optimization study. Figure 2(a) shows the effects of different reaction time and agitation speed on the conversion % of product in threedimensional surface response. Generally, increased reaction time and agitation speed resulted in an increase percentage of conversion until agitation speed reached 523 r.p.m. The response started to decrease after the agitation speed exceeded 523 r.p.m. even at the higher reaction time. However, it was observed that reaction time showed a significant effect to the reaction conversion at the higher agitation speed. Increasing agitation speed had increased the external mass transfer rates between the bulk phase of the reaction mixture and surface of enzyme; moreover, higher reaction time also promoted collision time between enzyme and substrate molecules. As shown in Figure 2(b), the reaction with the enzyme amount of 5.80% w/w led to the maximum percentage of conversion. Response surface plot for interaction between enzyme amount and reaction temperature was generated with reaction time fixed at 16 h, the molar ratio of substrates (OA : TEA) 2 : 1 mole, and agitation speed 400 r.p.m. The percentage conversion of product increased by increase ongoing from 3 to 5.80% w/w and thereafter decreased with further increase to 7% w/w. However, higher temperatures tended to induce enzyme inactivation due to denaturation processes [22, 23]. These results were similar to those in most reviewed papers, namely, that Novozym 435 was optimally used at temperatures between 40°C and 60°C [24, 25]. Figure 2(c) represents the effect of varying amount of enzyme and agitation speed on the synthesis of TEAbased esterquat with constant condition for other independent variables (reaction temperature of 60°C, reaction time of 16 h, and substrate molar ratio of 2 : 1 mole). From Figure 2(c), while the enzyme amount and agitation speed increased, the conversion of esterquat was increased as the agitation speed reached 523 r.p.m. in the enzyme amount of 5.80% w/w. However, the effect of enzyme amount variable was lower than the effect of agitation speed variable. Increase in agitation speed caused the substantial increase in the specific interfacial area between the substrate and the enzyme present in the nonaqueous phase by reducing the droplet size [26, 27]. A negative effect in percentage of conversion was detected with agitation speed greater than 523 r.p.m. This may be due to adverse shear effect caused by impeller at higher agitation speed. Typically, the immobilized enzyme was driven radially from impeller against the wall of the reactor, forcing the breakage, especially at high agitation speed [28]. Finally, Figure 2(d) shows the effect of varying the amounts of enzyme and molar ratio of substrates on the esterification reaction of oleic acid and triethanolamine while reaction time and reaction temperature are fixed at 16 h and 60°C, respectively. It was shown that the maximum conversion of esterquat was obtained when the enzyme amount was 11.6 g and increased with the lower molar ratio of substrates. However, increase in acyl donor showed less significant increase in the esterification conversion, on the other hand, and resulted in slight decrease of percentage conversion at the high amount of enzyme 14 g and the molar ratio of substrates of 3 : 1 mole. This was due to the limiting factor caused by triethanolamine, which was significant at the high amount of oleic acid and hence reduced the percentage of conversion.
(a)
(b)
(c)
(d)
3.4. Optimization by Response Surface Methodology and Model Validation
The next step in the present study was to determine the effects of five independent variables (enzyme amount, reaction time, reaction temperature, molar ratio of substrates, and agitation speed) shown in Table 6, along with the mean predicted values for enzymatic reaction product. For this purpose, the response surface methodology, using a central composite design, was adopted for finding optimal conditions. Experiment was then carried out under the recommended conditions and resulting response was compared to the predicted values. The optimum reaction parameters were enzyme amount of 4.77% w/w, reaction time of 24 h, reaction temperature of 61.9°C, substrates molar ratio (OA : TEA) of 1 : 1 mole (0.708 mole of OA and TEA), and agitation speed of 480 r.p.m. Comparison between RSM and ANN methods was then assessed in optimum conditions point for enzymatic synthesis of TEAbased esterquat at 2000 mL scale. The reaction of experiment gave the reasonable percentage of conversion 63.57%. This result confirmed the validity of the model, and the experimental value was determined to be quite close to the predicted value (65.08%) in comparison with ANN result (61.14%), implying that empirical model derived from RSM experimental design can be used to adequately describe the relationship between the independent variables and response.

4. Conclusion
In the present paper, RSM was used to optimize the enzymatic reaction conditions. A central composite design was applied to optimize the experimental conditions for synthesis of TEAbased esterquat at 2000 mL scale. The normality test was investigated as an initial step in process capability studies for better results and higher accuracy. Considering normality tests, the results indicated that all of the data and distributions were close to expected values under normality. The variables include enzyme amount, reaction time, reaction temperature, substrates molar ratio, and agitation speed. Quadratic mathematical model was suggested for synthesis of TEAbased esterquat. Analysis of variance corroborates the accuracy of the model by using high value (33.60), very low value (<0.0001), nonsignificant lack of fit, and the coefficient of determination (). A conversion percentage of 63.57% was attained, which was good compared to the predicted amount of 65.08%, with the relative standard error percentage (RSE) 2.32%. The comparison of RSM and ANN (QP) indicated that the RSM had less RSE% rather than ANN (QP) method (3.98%). The methodology as a whole has proven that RSM is adequate for the design and optimization of the enzymatic process.
Acknowledgment
The financial assistance provided by Universiti Putra Malaysia under the Research University Grant Scheme (RUGS) is gratefully acknowledged.
References
 F. E. Friedli, R. Keys, C. J. Toney, O. Portwood, D. Whittlinger, and M. Doerr, “Novel new ester quaternaries for improved performance benefits as rinse cycle fabric softeners,” Journal of Surfactants and Detergents, vol. 4, no. 4, pp. 401–405, 2001. View at: Google Scholar
 Z. Miao, J. Yang, L. Wang et al., “Synthesis of biodegradable lauric acid ester quaternary ammonium salt cationic surfactant and its utilization as calico softener,” Materials Letters, vol. 62, no. 19, pp. 3450–3452, 2008. View at: Publisher Site  Google Scholar
 J. Waters, H. H. Kleiser, M. J. How et al., “A new rinse conditioner active with improved environmental properties,” Tenside, Surfactants, Detergents, vol. 28, no. 6, pp. 460–468, 1991. View at: Google Scholar
 R. Puchta, P. Krings, and P. Sandkuehler, “New generation of softeners,” Tenside, Surfactants, Detergents, vol. 30, no. 3, pp. 186–191, 1993. View at: Google Scholar
 M. Levinson, “Rinseadded fabric softener technology at the close of the twentieth century,” Journal of Surfactants and Detergents, vol. 2, no. 2, pp. 223–235, 1999. View at: Google Scholar
 S. KoohiKamali, C. P. Tan, and T. C. Ling, “Optimization of sunflower oil transesterification process using sodium methoxide,” The Scientific World Journal, vol. 2012, Article ID 475027, 8 pages, 2012. View at: Publisher Site  Google Scholar
 Y. Abdollahi, A. Zakaria, A. H. Abdullah et al., “Semiempirical study of orthocresol photo degradation in manganesedoped zinc oxide nanoparticles suspensions,” Chemistry Central Journal, vol. 6, no. 1, p. 88, 2012. View at: Publisher Site  Google Scholar
 S. H. Musa, M. Basri, H. R. F. Masoumi et al., “Formulation optimization of palm kernel oil esters nanoemulsionloaded with chloramphenicol suitable for meningitis treatment,” Colloids and Surfaces B, vol. 112, pp. 113–119, 2013. View at: Google Scholar
 B. Subha and M. Muthukumar, “Optimization of ozonation process for the reduction of excess sludge production from activated sludge process of sago industry wastewater using central composite esign,” The Scientific World Journal, vol. 2012, Article ID 239271, 8 pages, 2012. View at: Publisher Site  Google Scholar
 R. Siala, F. Frikha, S. Mhamdi, M. Nasri, and A. S. Kamoun, “Optimization of acid protease production by aspergillus niger I1 on shrimp peptone using statistical experimental design,” The Scientific World Journal, vol. 2012, Article ID 564932, 11 pages, 2012. View at: Publisher Site  Google Scholar
 H. R. F. Masoumi, A. Kassim, M. Basri, D. K. Abdullah, and M. J. Haron, “Multivariate optimization in the biosynthesis of a triethanolamine (tea)Based esterquat cationic surfactant using an artificial neural network,” Molecules, vol. 16, no. 7, pp. 5538–5549, 2011. View at: Publisher Site  Google Scholar
 D. C. Montgomery, Design and Analysis of Experiments, John Wiley & Sons, New York, NY, USA, 2009.
 J. Pallant, SPSS Survival Manual: A Step by Step Guide to Data Analysis Using SPSS, Allen & Unwin, Crows Nest, Australia, 2005.
 H. M. Park, Univariate Analysis and Normality Test Using SAS, STATA, and SPSS. The University Information Technology Services (UITS) Center for Statistical and Mathematical Computing, Indiana University, Bloomington, Ind, USA, 2008.
 A. R. Henderson, “Testing experimental data for univariate normality,” Clinica Chimica Acta, vol. 366, no. 12, pp. 112–129, 2006. View at: Publisher Site  Google Scholar
 H. R. F. Masoumi, A. Kassim, M. Basri, and D. K. Abdullah, “Determining optimum conditions for lipasecatalyzed synthesis of triethanolamine (TEA)based esterquat cationic surfactant by a Taguchi robust design method,” Molecules, vol. 16, no. 6, pp. 4672–4680, 2011. View at: Publisher Site  Google Scholar
 N. Chaibakhsh, M. B. Abdul Rahman, S. AbdAziz, M. Basri, A. B. Salleh, and R. N. Z. R. A. Rahman, “Optimized lipasecatalyzed synthesis of adipate ester in a solventfree system,” Journal of Industrial Microbiology and Biotechnology, vol. 36, no. 9, pp. 1149–1155, 2009. View at: Publisher Site  Google Scholar
 F. J. Gravetter and L. B. Wallnau, Statistics for the Behavioral Sciences, 5th edition, 2000.
 A. Sher, “Assessing and comparing interaction dynamics, student learning, and satisfaction within Webbased online learning programs,” MERLOT Journal of Online Learning and Teaching, vol. 4, no. 4, pp. 446–458, 2008. View at: Google Scholar
 A. H. Hamzaoui, B. Jamoussi, and A. M'nif, “Lithium recovery from highly concentrated solutions: response surface methodology (RSM) process parameters optimization,” Hydrometallurgy, vol. 90, no. 1, pp. 1–7, 2008. View at: Publisher Site  Google Scholar
 Y. Li, J. Lu, G. Gu, and Z. Mao, “Characterization of the enzymatic degradation of arabinoxylans in grist containing wheat malt using response surface methodology,” Journal of the American Society of Brewing Chemists, vol. 63, no. 4, pp. 171–176, 2005. View at: Publisher Site  Google Scholar
 S. H. Krishna, S. Divakar, S. G. Prapulla, and N. G. Karanth, “Enzymatic synthesis of isoamyl acetate using immobilized lipase from Rhizomucor miehei,” Journal of Biotechnology, vol. 87, no. 3, pp. 193–201, 2001. View at: Publisher Site  Google Scholar
 S. M. Radzi, M. Basri, A. B. Salleh et al., “Large scale production of liquid wax ester by immobilized lipase.,” Journal of Oleo Science, vol. 54, no. 4, pp. 203–209, 2005. View at: Publisher Site  Google Scholar
 S. BourgGarros, N. Razafindramboa, and A. A. Pavia, “Optimization of lipasecatalyzed synthesis of (z)3hexen1yl acetate by direct esterification in hexane and a solventfree medium,” Enzyme and Microbial Technology, vol. 22, no. 4, pp. 240–245, 1998. View at: Publisher Site  Google Scholar
 P. Lozano, T. De Diego, D. Carrié, M. Vaultier, and J. L. Iborra, “Enzymatic ester synthesis in ionic liquids,” Journal of Molecular Catalysis B, vol. 21, no. 12, pp. 9–13, 2003. View at: Publisher Site  Google Scholar
 I. M. Noor, M. Hasan, and K. B. Ramachandran, “Effect of operating variables on the hydrolysis rate of palm oil by lipase,” Process Biochemistry, vol. 39, no. 1, pp. 13–20, 2003. View at: Publisher Site  Google Scholar
 S. AlZuhair, M. Hasan, and K. B. Ramachandran, “Kinetics of the enzymatic hydrolysis of palm oil by lipase,” Process Biochemistry, vol. 38, no. 8, pp. 1155–1163, 2003. View at: Publisher Site  Google Scholar
 P. S. Keng, M. Basri, M. B. A. Rahman, A. B. Salleh, R. N. Z. A. Rahman, and A. Ariff, “Optimization of palmbased wax esters production using statistical experimental designs,” Journal of Oleo Science, vol. 54, no. 10, pp. 519–528, 2005. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2013 Hamid Reza Fard Masoumi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.