Application of Box-Behnken Design in Optimization of Glucose Production from Oil Palm Empty Fruit Bunch Cellulose

Oil palm empty fruit bunch fiber (OPEFB) is a lignocellulosic waste from palm oil mills. It contains mainly cellulose from which glucose can be derived to serve as raw materials for valuable chemicals such as succinic acid. A three-level Box-Behnken design combined with the canonical and ridge analysis was employed to optimize the process parameters for glucose production from OPEFB cellulose using enzymatic hydrolysis. Organosolv pretreatment was used to extract cellulose from OPEFB using ethanol and water as the solvents. The extracted cellulose was characterized by thermogravimetric analysis, FTIR spectroscopy, and field emission scanning electron microscopy. Hydrolysis parameters including amount of enzyme, amount of cellulose, and reaction time were investigated. The experimental results were fitted with a second-order polynomial equation by a multiple regression analysis and found that more than 97% of the variations could be predicted by the models. Using the ridge analysis, the optimal conditions reaction time found for the production of glucose was 76 hours and 30 min, whereas the optimum amount of enzyme and cellulose was 0.5mL and 0.9 g, respectively. Under these optimal conditions, the corresponding response value predicted for glucose concentration was 169.34 g/L, which was confirmed by validation experiments.


Introduction
Bioconversion of lignocellulosic waste materials to chemicals and fuels are receiving interest as they are low cost, renewable, and widespread in nature [1].Malaysia is well acknowledged for its potential in renewable resources of lignocellulosic materials such as oil palm waste, sugar cane bagasse, and rice straw.At present Malaysia is the largest exporter and producer of palm oil and its production accounts approximately 40-60% of world total palm oil over the 25 years [2,3].In the process of extraction of palm oil from oil palm fruit, a lignocellulosic material, namely, OPEFB, is generated as a main solid residue, with every ton of fresh fruit bunches producing ∼0.22 ton of OPEFB [4].In Malaysia, ∼17 million tons of OPEFB is produced every year.Commonly, this biomass is burnt in incinerators, and it does not only create environmental pollution problems but it also offers limited value to the industry.Several approaches have been developed to utilize OPEFB fibres to produce valuable materials including paper pulp, composite boards, thermoset polymer, and activated carbon [5].
OPEFB consists of 67% of holocellulose (cellulose and hemicellulose) and 24% of lignin [6].Cellulose (-glucan) is a polymer of glucose [7] and can be used as a substrate for the production of a wide variety of compounds by chemical and biochemical processes.One such compound is succinic acid, which is extensively used in food, pharmaceuticals, and manufacturing of biodegradable polymers [8].Production of glucose from OPEFB cellulose increased the yield and concentration of glucose produced compared to using acid or base pretreated OPEFB as reported before.Enzymatic hydrolysis using cellulase from Trichoderma reesei (a mixture of endoglucanases, exoglucanases, and -glucosidase) is widely used for the degradation of cellulose into soluble oligo glucan/oligomers, cellobiose, and glucose [9].In this project, a mixture of two enzymes was used to carry out the reaction.Novozyme Celluclast 1.5 L is a brand for cellulase enzyme that is capable to break the cellulose chains, particularly at middle and chain end sections into cellobiose (the dimer of glucose).Whereas, Novozyme 188 that mainly contains -glucosidase was used to optimize the glucose production since it breaks the cellobiose chains into glucose.
Optimization of parameters in reaction is one of the most important stages in the development of an efficient and economic production of high value products from renewable sources.The traditional "one-factor-at-a-time approach" is time consuming, and moreover the interactions between independent variables are not considered.Response surface methodology (RSM) is an effective optimization tool wherein many factors and their interactions can be identified with fewer experimental trials [10].
RSM has been widely used in various fields including food process operations, new product development, biotechnology-media composition, and bioprocessing such as enzymatic hydrolysis and fermentation [10].The RSM is a collection of mathematical and statistical techniques for designing experiments, building models, evaluating the effects of factors, and searching optimum condition of factors for desired responses.The optimization process of this methodology involves studying the response of statistically designed combinations, estimating the coefficients by fitting it in mathematical model that fits best the experimental conditions, predicting the response of the fitted model, and checking the adequacy of the model.Central composite design (CCD) and Box-Behnken design (BBD) are amongst the most commonly used in various experiments [11].Box-Behnken design allows calculation of the response to be made at intermediate levels which were not experimentally studied.A three-level Box-Behnken design was employed in the present study and the optimal conditions were determined through a minimal experiment number compared with other designs [12].
The objective of this study is to optimize the glucose production from OPEFB cellulose via enzymatic hydrolysis with the help of Box-Behnken design, followed by canonical and ridge analyses.OPEFB was treated by organosolv method for lignin and hemicellulose removal.In this context, reaction time, amount of enzyme, and amount of substrates were optimized.Enzymatic hydrolysis was also carried out on commercially available cellulose for comparison purpose.

Materials and Methods
2.1.Raw Material.OPEFB fibers were collected from local palm oil mill (United Oil Palm Industries SdnBhd, Malaysia), sun dried, and ground to particles with diameter of <1 mm.The OPEFB biomass was then oven dried at 105 ∘ C overnight and was analyzed following the standard method for determination of its main composition [13].

Chemicals and Enzymes.
Two types of standard cellulose from Sigma-Aldrich were used in this study, namely, fiber cellulose and microcrystalline cellulose.Solvent and other chemicals were obtained from R&M Chemicals.Celluclast 1.5 L (cellulase) and Novozyme 188 (cellobiase) were obtained from Novozymes Malaysia SdnBhd.

Autonydrolysis Pretreatment.
Autohydrolysis pre-treatment was conducted for the delignification process of OPEFB.OPEFB fibers (30 grams) were loaded into 4 L stainless steel reactor (98 kPa, 120 ∘ C) and were supplemented with appropriate amount of deionized water.The autohydrolysis was carried out at 120 ∘ C for 1-2 hours.
2.5.Enzymatic Hydrolysis.Enzymatic hydrolysis experiment was carried out in a 100 mL shaking flask at 40 ∘ C and 145 rpm inside an incubator shaker.The cellulase-catalysed hydrolysis of different cellulose substrates (untreated or pretreated) was carried out in a stirred flask.In a typical hydrolysis reaction, 500 mg of cellulose was added to 9 mL acetate buffer (50 mm, pH 4) and incubated for 2 hours (40 ∘ C; 145 rpm).After this preincubation step, hydrolysis was initiated by adding 0.1 mL of 10 mg/mL cellulase (activity: 700 EGU/g).The reaction medium was withdrawn repeatedly at a specific time interval to determine the progress of the reaction which was stopped by incubating the withdrawn sample at 90 ∘ C for 20 min [15].Then, the sample was diluted in ultrapure water and filtered (0.2 m) prior to analysis by high performance liquid chromatography (HPLC).
The same procedure was also performed for the cellulose extracted from OPEFB.The reaction was carried out at 40 ∘ C, for 48-94 hours.The amount of cellulose and enzyme used was varied from 300 mg to 500 mg and from 0.3 mL to 0.7 mL, respectively.

Box-Behnken Design.
Based on the preliminary onefactor-at-a-time experiments, variables such as duration of reaction, amount of samples, and total enzymes were identified to have strong effects on the response.Therefore, these factors were selected as the variables tested in the 15-run experiment of the Box-Behnken design experiment.
Table 1 shows the Box-Behnken design matrix of the experiment of 15 trials.As shown in Table 1, the three factors chosen for this study were designated as  1 ,  2 , and  3 and prescribed into three levels, coded +1, 0, −1 for high, intermediate, and low value, respectively.Three test variables were coded according to the following equation: where   is the coded value of an independent variable;   is the actual value of an independent variable;   is the actual value of an independent variable at centre point; and Δ is the step change value of an independent variable.All experiments were carried out in triplicate, and the averages of glucose concentration were taken as response [12].For predicting the optimal point, a second-order polynomial model was fitted to correlate the relationship between independent variables and response (glucose concentration).The equation for three factors is stated as follows: where  is the predicted response;  0 is model constant; The glucose concentration was determined by using High Performance Liquid Chromatography (HPLC) and glucose meter.The system used for HPLC was Waters HPLC system, and evaporative light scattering detectors (ELSD) was used as detectors.Glucose was determined using a Zorbax NH 2 column and acetonitrile : water 80 : 20 as the mobile phase at a flow rate of 1 mL/min.The retention time of glucose is 13.32 min (rsd ± 0.56%).Quantification was based on calibration curves established using standard glucose purchased from Sigma-Aldrich.

Characterization of OPEFB Cellulose
3.1.1.FTIR Spectroscopy.Figure 1 shows the FTIR spectra of cellulose extracted from OPEFB (OPEFB cellulose), raw OPEFB, and delignified OPEFB.The peaks ascribed to lignin appeared at 1501 and 1512 cm −1 that are due to C=C stretching and C=C aromatic skeletal vibration of lignin, respectively.However, the two peaks of lignin were not observed in the FTIR spectra of the delignified OPEFB (Figure 1,b) and OPEFB cellulose (Figure 1,c).The peak attributed to hemicellulose was observed in the spectra of the raw OPEFB and delignified OPEFB (Figure 1,b) at 1732-1735 cm −1 due to C=O stretching.However, the corresponding peaks were disappeared in the spectrum of OPEFB cellulose (Figure 1,c).This observation indicated that the organosolv treatment was able to remove the lignin as well as hemicellulose from OPEFB fibers.The peak observed at 898 cm −1 is attributed to the presence of -glucoside linkage between glucose units in cellulose [16].All FTIR spectra in this study showed a glucoside linkage peak including the FTIR spectrum for the obtained OPEFB cellulose (Figure 1,c).lignin and hemicellulose were successfully removed.Thermal decomposition at around 310 ∘ C showed that OPEFB cellulose was successfully obtained (Figure 2(a)).

Morphology of OPEFB Fibers. SEM micrographs of raw
OPEFB, delignified OPEFB, and obtained cellulose (OPEFB cellulose) are shown in Figure 3.A rigid appearance can be seen for raw OPEFB material (Figure 3(a)).The figure shows that raw OPEFB fibers exhibit stiff and hard surfaces.Delignification of OPEFB reduced the stiff appearance and some parts of the fibers were split and became more refined.The images also showed that some fibers were broken.Figure 3(c) shows the SEM image of obtained cellulose (OPEFB cellulose) after delignification and organosolv treatment.The appearance is quite different and the length and diameter of fibers were reduced.This is probably due to the decrease in spiral angle around the fibre axis and increase in molecular orientation after pretreatment.A fair amount of randomness is introduced to the orientation of the crystallites due to removal of noncellulosic matters [15] leading to the formation of fibers with smaller diameter and length.

Effect of Pretreatments on Enzymatic Hydrolysis of Cellulose.
Effect of pretreatments on enzymatic hydrolysis of various cellulose substrates is shown in Figure 4. Plots in Figure 4 show that glucose concentration increased with reaction time.Figure 4 also indicated that organosolv treated OPEFB cellulose produced the highest glucose concentration followed by acetic acid treated OPEFB, microcrystalline cellulose, cellulose fiber, and delignified OPEFB.The concentration of glucose produced from the organosolv treated OPEFB cellulose was 36.5 g/L with 48 hours reaction time and used minimum amount of enzyme.The sugar recovery for organosolv treated OPEFB cellulose is 0.73 g/g.Table 2 shows the effect of different parameters on enzymatic hydrolysis of cellulose extracted from OPEFB.It was observed that model run number 8 produced the highest amount of glucose, that is, 167.4 g/L.The lowest glucose concentration was produced by run number 9 with concentration of only 21.6 g/L.The glucose concentration increased with the increasing amount of substrate.The longer reaction time also increased the glucose concentration.The best total enzyme to produce high glucose concentration is 0.2 mL.This means that the use of the same amount of novozyme 188 and celluclast (0.1 mL : 0.1 mL) can improve the activity of the two enzymes as catalyst.appropriate second-order responses surface.The three centre point runs were added to provide as a measure of process stability and inherent variability.The considerable variation in the glucose concentration produced from OPEFB cellulose under different conditions is shown in Table 3.The glucose concentration predicted by the final quadratic model along with corresponding values observed was given too.The agreement between the yield predicted by the model and the experimental data is very strong, with a small difference in glucose concentration.Furthermore, the model capability was checked using the -test and determination coefficient  2 .The analysis of variance (Table 4) showed that this regression model was  (3)

Canonical and Ridge Analysis.
The canonical analysis of the response surface was performed with SAS to determine the shape of the fitted response and the estimated stationary point.According to the model, the predicted response at the The three eigenvalues had different signs, indicating that the stationary point for this model was a saddle point as shown in Table 7.Therefore, the estimated surface did not have a unique optimum, and a ridge analysis was performed to determine the optimum.
The results of ridge analysis (Table 8) indicated that all variables tested, that is, hydrolysis reaction time, total of enzymes, and amount of substrate, were positively related to the response, and the optimal level of them was found as 76 hours and 30 min reaction time, 0.5 mL total enzymes, and 0.9 g of substrate (OPEFB cellulose), respectively, with a predicted glucose concentration of 169.34 g/L.
In order to confirm the predicted values, three further experiments using the optimum parameters found above were conducted and the value of 165.24-172.5 g/L (mean value 168.71 g/L) of glucose concentration was obtained.
The strong relationship between the yield predicted by the final quadratic model and the experimental results (Table 2) and the variance analysis of the second-order polynomial model and the value of lack of fit (Table 3) showed that the accuracy and overall ability of the polynomial model are very good.The analysis of response trends using the model is considered to be acceptable.
Kunamneni and Singh [18] (2005) used a central composite design (CCD) which is one of the models from response surface methodology (RSM) to optimize the effect of four variables (precooking-amylase dose, postcooking-amylase dose, glucoamylase dose, and saccharification temperature) in the optimization of glucose production from maize starch [18].As comparison, they did 28 experiments for 4 factors and finally obtained a coefficient determination of 0.9558, whereas in this project only 3 factors with 15 experiments were used and the coefficient determination value obtained is 0.9765.Although many experiments are completed at the canonical analysis, the ridge analysis is very useful to find out the maximum response which may occur in the experiment when the results of the canonical analysis shows a saddle, stationary point, and no unique optimum in the estimated surface [19].Based on the ridge analysis, the optimal reaction time, total of enzymes, and amount of samples for glucose production were determined in this study with a good predicted hydrolysis reaction rate without further experimental work.An average glucose concentration of 168.71 g/L was obtained in three further validation experiments using the optimum reaction parameters.As compared to those of experimental (167.4g/L) and predicted value (169.34 g/L), there is a good correlation amongst those three results verifying the validity of the response model and the actuality of optimal point.

Conclusions
Enzymatic hydrolysis of OPEFB cellulose was successfully carried out at 40 ∘ C under various amounts of samples, total enzymes, and reaction time (48-94 hours).Box-Behnken design was shown to be useful in investigating those three critical variables.Through the canonical and ridge analyses of the second-order polynomial model, a maximum glucose concentration of 168.71 g/L was obtained under the following conditions: 76 hours and 30 minutes hydrolysis reaction time used 0.5 mL enzymes and 0.9 gram samples/OPEFB cellulose.In summary, within the reaction conditions tested the concentration of glucose produced is proportional to the reaction time, total enzymes, and amount of samples.

Table 1 :
Levels and code of variables chosen for Box-Behnken design.

Table 2 :
Effect of different parameters on enzymatic hydrolysis of OPEFB cellulose.

Table 3 :
Box-Behnken design with experimental and predicted values of glucose concentration.

Table 4 :
Analysis of variance in the regression model for optimisation of glucose production from OPEFB cellulose.

Table 5 :
Regression analysis of a full second-order polynomial model for optimisation of glucose production from OPEFB cellulose.

Table 5 .
The  values were used as a device to check the significance of each coefficient, and the  values showed that  2 1 and  2 2 as well as two cross products ( 1  3 and  2  3 ) were significant. 1 ) − (42.8 3  2 ) .

Table 6 :
Canonical analysis on the coded and actual value.

Table 8 :
Results of ridge analysis.