Synthesis of Lipase-Immobilized CeO 2 Nanorods as Heterogeneous Nano-Biocatalyst for Optimized Biodiesel Production from Eruca sativa Seed Oil

: Biodiesel has emerged as one of the most attractive alternative energy sources to meet the growing needs of energy. Many approaches have been adopted for biodiesel synthesis. In the present work, biodiesel was produced from non-edible Eruca sativa oil using nano-biocatalyst-catalysed transesteriﬁcation. Nano-biocatalyst (CeO 2 @PDA@ A. terreus Lipase) was developed via the immobilization of lipase on polydopamine coated ceria nanorods, and CeO 2 nanorods were developed via a hydrothermal process. The mean diameter of nanorods were measured to be 50–60 nm, while their mean length was 150–200 nm. Lipase activity before and after immobilization was measured to be 18.32 and 16.90 U / mg / min, respectively. The immobilized lipase depicted high stability at high temperature and pH. CeO 2 @PDA@ A. terreus Lipase-catalysed transesteriﬁcation resulted in 89.3% yield of the product. Process optimization through response surface methodology was also executed, and it was depicted that the optimum / maximum E. sativa oil-based biodiesel yield was procured at conditions of 10% CeO 2 @PDA@ A. terreus Lipase, 6:1 methanol / oil ratio, 0.6% water content, 35 ◦ C reaction temperature, and 30 h reaction time. The fuel compatibility of synthesized biodiesel was conﬁrmed via the estimation of fuel properties that were in agreement with the ASTM D standard. The nanorods and dopamine-modiﬁed nanorods were characterized by FTIR spectroscopy, SEM, and energy dispersive X-ray (EDX), while conversion of E. sativa oil to biodiesel was conﬁrmed by GC / MS and FTIR spectroscopy. Conclusively, it was revealed that CeO 2 @PDA@ A. terreus Lipase has potential to be employed as an emphatic nano-biocatalyst.


Energy Dispersive X-Ray Spectroscopy (EDX)
The constituent elements of the nanorods were determined using energy dispersive X-ray (EDX) spectroscopy. The plot between energy of X-rays vs. X-rays counts is presented in (Figure 2). Each peak in the plot corresponds to a specific element. EDX plot confirms the presence of cerium and oxygen in the product. Absence of extra peaks may be attributed to the purity of synthesized nanorods, while peak for carbon is due to the carbon plates of instrument. Atomic% and weight% of elements obtained by EDX are shown in Table 1. Ratio of atomic% between cerium and oxygen affirms the formation of CeO2 nanorods.

X-Ray Diffraction (XRD)
XRD pattern (Figure 3) was used to characterize the CeO2 nano-powder, and instrument was operated within range of 20-80°. The peaks appeared in the XRD pattern with 2θ values at 27.

Energy Dispersive X-Ray Spectroscopy (EDX)
The constituent elements of the nanorods were determined using energy dispersive X-ray (EDX) spectroscopy. The plot between energy of X-rays vs. X-rays counts is presented in (Figure 2). Each peak in the plot corresponds to a specific element. EDX plot confirms the presence of cerium and oxygen in the product. Absence of extra peaks may be attributed to the purity of synthesized nanorods, while peak for carbon is due to the carbon plates of instrument. Atomic% and weight% of elements obtained by EDX are shown in Table 1. Ratio of atomic% between cerium and oxygen affirms the formation of CeO 2 nanorods.

Energy Dispersive X-Ray Spectroscopy (EDX)
The constituent elements of the nanorods were determined using energy dispersive X-ray (EDX) spectroscopy. The plot between energy of X-rays vs. X-rays counts is presented in (Figure 2). Each peak in the plot corresponds to a specific element. EDX plot confirms the presence of cerium and oxygen in the product. Absence of extra peaks may be attributed to the purity of synthesized nanorods, while peak for carbon is due to the carbon plates of instrument. Atomic% and weight% of elements obtained by EDX are shown in Table 1. Ratio of atomic% between cerium and oxygen affirms the formation of CeO2 nanorods.

X-Ray Diffraction (XRD)
XRD pattern (Figure 3) was used to characterize the CeO2 nano-powder, and instrument was operated within range of 20-80°. The peaks appeared in the XRD pattern with 2θ values at 27.

X-Ray Diffraction (XRD)
XRD pattern (Figure 3) was used to characterize the CeO 2 nano-powder, and instrument was operated within range of 20-80 • . The peaks appeared in the XRD pattern with 2θ values at 27

Fourier Transform Infrared Spectroscopic (FTIR) Analysis of Naked and Polydopamine Coated CeO2 Nanorods
Both the polydopamine-coated and uncoated CeO2 nanorods were scanned for FTIR spectroscopy in the range of 600-4000 cm −1 (Figure 4). FTIR spectrum (blue colored) reveals the functional groups and the chemical bonds that are present in the synthesized, uncoated CeO2 nanorods, whereas FTIR (black colored) spectrum confirms the functionalization of CeO2 nanorods with polydopamine layer. The broad band at 3402.47 cm −1 corresponds to (the stretching vibration of O-H bond) OH-groups. The peak around 1569.67 cm −1 is ascribed to the bending vibration of N-H group. The absorption peak around 1496.92 cm −1 indicates the -CH2 vibrations. The intense band at 691.17 cm −1 corresponds to the Ce-O stretching vibrations [11]; another additional peak at 1312.27 cm −1 corresponds to the C-O bonds in dopamine molecule that verifies the dopamine addition on CeO2 nanorods [12].

Lipase Activity Assay
The activity titer of free lipase (which was produced from Aspergillus terreus AH-F2) was found to be 18.32 U/mg/min, while the activity titer of CeO2@PDA@A. terreus Lipase was found to be 16.90 U/mg/min. It was found that, by using polydopamine, high immobilization efficiency was achieved by using CeO2 nanorods because polydopamine and CeO2 nanorods formed the complex, which was found to be efficient regarding immobilization of lipase. The wider surface was available during the reaction, which gave higher conversion rate in a short period of time by using the lipase-immobilized nanorods for biodiesel production [13]. To the best of our knowledge, the synthesized nano-

Fourier Transform Infrared Spectroscopic (FTIR) Analysis of Naked and Polydopamine Coated CeO 2 Nanorods
Both the polydopamine-coated and uncoated CeO 2 nanorods were scanned for FTIR spectroscopy in the range of 600-4000 cm −1 (Figure 4). FTIR spectrum (blue colored) reveals the functional groups and the chemical bonds that are present in the synthesized, uncoated CeO 2 nanorods, whereas FTIR (black colored) spectrum confirms the functionalization of CeO 2 nanorods with polydopamine layer. The broad band at 3402.47 cm −1 corresponds to (the stretching vibration of O-H bond) OH-groups. The peak around 1569.67 cm −1 is ascribed to the bending vibration of N-H group. The absorption peak around 1496.92 cm −1 indicates the -CH 2 vibrations. The intense band at 691.17 cm −1 corresponds to the Ce-O stretching vibrations [11]; another additional peak at 1312.27 cm −1 corresponds to the C-O bonds in dopamine molecule that verifies the dopamine addition on CeO 2 nanorods [12].

Fourier Transform Infrared Spectroscopic (FTIR) Analysis of Naked and Polydopamine Coated CeO2 Nanorods
Both the polydopamine-coated and uncoated CeO2 nanorods were scanned for FTIR spectroscopy in the range of 600-4000 cm −1 (Figure 4). FTIR spectrum (blue colored) reveals the functional groups and the chemical bonds that are present in the synthesized, uncoated CeO2 nanorods, whereas FTIR (black colored) spectrum confirms the functionalization of CeO2 nanorods with polydopamine layer. The broad band at 3402.47 cm −1 corresponds to (the stretching vibration of O-H bond) OH-groups. The peak around 1569.67 cm −1 is ascribed to the bending vibration of N-H group. The absorption peak around 1496.92 cm −1 indicates the -CH2 vibrations. The intense band at 691.17 cm −1 corresponds to the Ce-O stretching vibrations [11]; another additional peak at 1312.27 cm −1 corresponds to the C-O bonds in dopamine molecule that verifies the dopamine addition on CeO2 nanorods [12].

Lipase Activity Assay
The activity titer of free lipase (which was produced from Aspergillus terreus AH-F2) was found to be 18.32 U/mg/min, while the activity titer of CeO2@PDA@A. terreus Lipase was found to be 16.90 U/mg/min. It was found that, by using polydopamine, high immobilization efficiency was achieved by using CeO2 nanorods because polydopamine and CeO2 nanorods formed the complex, which was found to be efficient regarding immobilization of lipase. The wider surface was available during the reaction, which gave higher conversion rate in a short period of time by using the lipase-immobilized nanorods for biodiesel production [13]. To the best of our knowledge, the synthesized nano-

Lipase Activity Assay
The activity titer of free lipase (which was produced from Aspergillus terreus AH-F2) was found to be 18.32 U/mg/min, while the activity titer of CeO 2 @PDA@A. terreus Lipase was found to be 16.90 U/mg/min. It was found that, by using polydopamine, high immobilization efficiency was achieved by using CeO 2 nanorods because polydopamine and CeO 2 nanorods formed the complex, which was found to be efficient regarding immobilization of lipase. The wider surface was available during the reaction, which gave higher conversion rate in a short period of time by using the lipase-immobilized nanorods for Catalysts 2020, 10, 231 5 of 16 biodiesel production [13]. To the best of our knowledge, the synthesized nano-biocatalyst with CeO 2 nanorods as immobilizing support for Lipase (Aspergillus terreus AH-F2) may be novel nano-biocatalyst for biodiesel synthesis.

Effect of pH and Temperature on Activity of Nano-Biocatalyst
The impact of pH on the activity of free lipase and synthesized nano-biocatalyst is presented in (Figure 5a). The effect of pH was studied within pH range of 5 to 10, and it was found that maximum lipase activity, i.e., 18.32 U/mg/min, was observed at pH 7.0 in case of free lipase, while at pH 8, CeO 2 @PDA@A. terreus Lipase showed maximum activity (i.e., 16.90 U/mg/min). The relationship plot between activities of Lipase (Free and immobilized) and pH depicted that CeO 2 @PDA@A. terreus Lipase could tolerate high pH values. The comparison of free and immobilized lipase revealed that immobilizing A. terreus Lipase on CeO 2 @PDA increased the flexibility of lipase to a wide pH range, compared to free lipase. Our results were comparable to the studies of Baharfar and Mahajer [13].

Effect of pH and Temperature on Activity of Nano-Biocatalyst
The impact of pH on the activity of free lipase and synthesized nano-biocatalyst is presented in (Figure 5a). The effect of pH was studied within pH range of 5 to 10, and it was found that maximum lipase activity, i.e., 18.32 U/mg/min, was observed at pH 7.0 in case of free lipase, while at pH 8, CeO2@PDA@A. terreus Lipase showed maximum activity (i.e., 16.90 U/mg/min). The relationship plot between activities of Lipase (Free and immobilized) and pH depicted that CeO2@PDA@A. terreus Lipase could tolerate high pH values. The comparison of free and immobilized lipase revealed that immobilizing A. terreus Lipase on CeO2@PDA increased the flexibility of lipase to a wide pH range, compared to free lipase. Our results were comparable to the studies of Baharfar and Mahajer [13].
The impact of temperature (ranging from 25 to 50 °C) on the activity of free lipase (A. terreus Lipase) and CeO2@PDA@A. terreus Lipase is presented in Figure 5b. It was depicted that highest free lipase activity was revealed at 30 °C, while the CeO2@PDA@A. terreus Lipase showed highest/maximum activity at 40 °C, so it showed that CeO2@PDA@A. terreus Lipase was tolerant to high temperature conditions and was stable at a wider temperature range. The increased tolerability of lipase may have been due to the formation of covalent bonds during immobilization. Comparable observations have been reported by Baharfar & Mahajer [13] and Dumri & Hung [14].

Physico-Chemical Properties of E. sativa Seed Oil
Physico-chemical properties of E. sativa oil are described below in (Table 2). Density is measure of mass per unit value of a sample (E. sativa oil) [15]. Peroxide value is the measure of amount of peroxides in the sample, i.e., oil or fats [16]. The ratio of density of sample to the water density is described by specific gravity [17]. Acid value is the amount of potassium hydroxide (mg) required for the deactivation of the fatty acids present in the 1.0 g of sample (oil) [18]. Iodine value is the measurement of unsaturation in the given oil; if the iodine value is high, then degree of unsaturation is high. The milligram of KOH required to saponify each gram of oil under specific set of conditions is referred as saponification value [19]. The impact of temperature (ranging from 25 to 50 • C) on the activity of free lipase (A. terreus Lipase) and CeO 2 @PDA@A. terreus Lipase is presented in Figure 5b. It was depicted that highest free lipase activity was revealed at 30 • C, while the CeO 2 @PDA@A. terreus Lipase showed highest/maximum activity at 40 • C, so it showed that CeO 2 @PDA@A. terreus Lipase was tolerant to high temperature conditions and was stable at a wider temperature range. The increased tolerability of lipase may have been due to the formation of covalent bonds during immobilization. Comparable observations have been reported by Baharfar & Mahajer [13] and Dumri & Hung [14].

Physico-Chemical Properties of E. sativa Seed Oil
Physico-chemical properties of E. sativa oil are described below in (Table 2). Density is measure of mass per unit value of a sample (E. sativa oil) [15]. Peroxide value is the measure of amount of peroxides in the sample, i.e., oil or fats [16]. The ratio of density of sample to the water density is described by specific gravity [17]. Acid value is the amount of potassium hydroxide (mg) required for the deactivation of the fatty acids present in the 1.0 g of sample (oil) [18]. Iodine value is the measurement of unsaturation in the given oil; if the iodine value is high, then degree of unsaturation is high. The milligram of KOH required to saponify each gram of oil under specific set of conditions is referred as saponification value [19].

Biodiesel Yield (%) and Optimum Reaction Conditions
Optimum reaction conditions for CeO 2 @PDA@A. terreus Lipase-catalyzed transesterification of E. sativa are given below in (Table 3), with maximum biodiesel yield of 89.3%. Optimum reaction conditions revealed were reaction time (30 h), temperature (35 • C), alcohol/oil ratio (6:1), water content (0.6%), and CeO 2 @PDA@A. terreus Lipase concentration (10%). The same enzyme (A. terreus Lipase) has also been immobilized on magnetite nano-support in another study performed by our research group, and comparable reaction conditions were revealed for optimal biodiesel production (92%) with slight variation in reaction temperature (37 • C). When methanol is added in oil, the viscosity of oil decreases and reaction rate increases, which is why methanol is significant for the reaction rate, but if an excessive amount of methanol is added in oil then emulsification of the glycerol may take place, deactivating the lipase and minimizing the biodiesel yield [20]. The other factor is catalyst concentration that influences the reaction rate. In the current study, 10% CeO 2 @PDA@A. terreus Lipase concentration was observed as optimum concentration for the synthesis of E. sativa seed oil-based biodiesel. The maximum yield of biodiesel (89.3%) may be increased by taking higher concentration of nano-biocatalyst or by using different lipase source for immobilization on functionalized CeO 2 . Another factor is reaction temperature, when the reaction is carried out at a more optimum temperature than reaction rate and biodiesel yield is maximum. However, high temperature denatured the lipase and reduced the biodiesel yield. Water content is also significant for the reaction rate, as it affects the activity and stability of lipase and protects it from deactivation via short chain alcohols [21].
Previously different researchers have employed response surface methodology (RSM) for biodiesel production using various feedstock oils and catalysts [22,23]. Aghababaie et al. has described bio-catalytic biodiesel synthesis from the crude E. sativa oil and obtained the highest FAME yield at 3:1 methanol to oil ratio, 5 mg lipase, 40% water content, and 21 • C temperature. The results that have been reported in literature are comparable to the present work [24][25][26][27], but still some variations are present that might be due to the varying fatty acid profiles of feedstock, different activities of the catalysts, and ranges of reaction parameters selected for optimization process.

RSM Model Fitting
Central composite design was applied to optimize biodiesel production procedure. Among the different models (viz; linear, quadratic, cubic polynomial, and two-factor interaction 2FI), quadratic model was revealed to be the most significant model for the responses with p-value < 0.0001 for CeO 2 @PDA@A. terreus Lipase catalyzed transesterification of sample oil. Fitness of quadratic model was also confirmed by the lack of fitness test, which was insignificant for the model with p-value 0.0728 in addition to higher R 2 and adjusted R 2 values (Table 4). Fitness of the model was also affirmed by the normality and predicted vs. actual values graph ( Figure 6). The linear distribution of data along the straight line in normality plot of the model depicts the fitness of quadratic models, while the small difference b/w predicted, and actual values of biodiesel yield further advocates the significance of quadratic models.
Catalysts 2020, 10, x FOR PEER REVIEW 7 of 16 was also confirmed by the lack of fitness test, which was insignificant for the model with p-value 0.0728 in addition to higher R 2 and adjusted R 2 values (Table 4). Fitness of the model was also affirmed by the normality and predicted vs. actual values graph ( Figure 6). The linear distribution of data along the straight line in normality plot of the model depicts the fitness of quadratic models, while the small difference b/w predicted, and actual values of biodiesel yield further advocates the significance of quadratic models.  Table 5 represents the ANOVA for the selected quadratic model describing the impact of various reaction parameters (as linear terms, 1st order interaction and quadratic terms) on response (% biodiesel yield) catalyzed by CeO2@PDA@A. terreus Lipase.

ANOVA for Biodiesel Yield Response
A (methanol/oil ratio), B (CeO2@PDA@A. terreus Lipase concentration), and D (reaction time) were significant linear terms, while C and E were non-significant (˃0.05) with p-value 0.1344 and 0.6091, respectively. The AB, AC, AE, and BC were significant interaction terms, having p-values of 0.0069, 0.0283, 0.0243, and 0.0013, respectively (<0.05). As for the quadratic terms, quadratic terms B 2 , C 2, and D 2 were the most significant quadratic terms, with p-values of 0.0001, 0.0001, and 0.0240, respectively, while A 2 and E 2 were non-significant quadratic terms that had p-values > 0.05.
Comparable results have been obtained in previous research. Mehmood et al. (2018) described the ANOVA for the significant quadratic model for biodiesel production from E. sativa oil catalysed by H2SO4 and reported that A (amount of catalyst) and D (methanol/oil ratio) impact significantly on the response (biodiesel yield), while B (temperature) and C (time) were non-significant terms. Regarding the 1st order interaction terms, BC and CD were significant, whereas AB, AC, AD, and BD were non-significant. Among the quadratic terms, D 2 and B 2 were significant [28].   Table 5 represents the ANOVA for the selected quadratic model describing the impact of various reaction parameters (as linear terms, 1st order interaction and quadratic terms) on response (% biodiesel yield) catalyzed by CeO 2 @PDA@A. terreus Lipase.
A (methanol/oil ratio), B (CeO 2 @PDA@A. terreus Lipase concentration), and D (reaction time) were significant linear terms, while C and E were non-significant (>0.05) with p-value 0.1344 and 0.6091, respectively. The AB, AC, AE, and BC were significant interaction terms, having p-values of 0.0069, 0.0283, 0.0243, and 0.0013, respectively (<0.05). As for the quadratic terms, quadratic terms B 2 , C 2, and D 2 were the most significant quadratic terms, with p-values of 0.0001, 0.0001, and 0.0240, respectively, while A 2 and E 2 were non-significant quadratic terms that had p-values > 0.05.
Comparable results have been obtained in previous research. Mehmood et al. (2018) described the ANOVA for the significant quadratic model for biodiesel production from E. sativa oil catalysed by H 2 SO 4 and reported that A (amount of catalyst) and D (methanol/oil ratio) impact significantly on the response (biodiesel yield), while B (temperature) and C (time) were non-significant terms. Regarding the 1st order interaction terms, BC and CD were significant, whereas AB, AC, AD, and BD were non-significant. Among the quadratic terms, D 2 and B 2 were significant [28]. Mumtaz et al. also reported the ANOVA for bio-catalytic synthesis of palm oil-based biodiesel. According to their findings, liner terms that showed a significant effect on response (% biodiesel) were A (amount of bio-catalyst), B (reaction time), and D (methanol/oil ratio). Among first-order interacting terms, AD showed significant influence on biodiesel yield. Among quadratic terms, B 2 and D 2 were depicted to be significant [29]. In another previous study, Chang et al. reported the analysis of variance for experimental data obtained by Lipase catalyzed biodiesel production from Novozym 435. The significant linear terms were X 2 (reaction temperature), X 3 (enzyme concentration), X 4 (substrate molar ratio), and X 5 (water content), while the non-significant term was X 1 (reaction time) [30].
Razack and Duraiarasan also reported ANOVA for the biodiesel production from waste cooking oil catalyzed by encapsulated mixed enzyme in which the only significant linear factor was C (reaction temperature); while A (enzyme), B (molar ratio) and D(time) had p-value > 0.05. Significant quadratic factors were revealed to be A 2 , B 2 , and C 2 . AC was the only significant first order interaction term [31].
The results reported previously are comparable with the results obtained in present work with some variation. The significance of methanol-to-oil ratio and enzyme concentration is obvious, as they directly influence the reaction rate. The enzymatic transesterification is a slow process so the significance of reaction time for enzymatic transesterification cannot be denied, as reported by other researchers as well [32]. Water is also considered as an imperative factor; however, water concentration (as linear term) has not been proved significant in the present work.

3D Surface Graphs for % Biodiesel Yield
Response surface plots of interaction terms depicting significant influence on response (biodiesel yield) are presented in (Figure 7). 3D graph describing the cumulative impact of CeO 2 @PDA@A. terreus Lipase and CH 3 OH/oil ratio indicated that the % biodiesel increases with an increase in nano-biocatalyst level and methanol/oil ratio, and optimal biodiesel yield is obtained when nano-biocatalyst level and methanol/oil ratio were 10% and 6:1, respectively. Deviation from these values, however, results in lower biodiesel yield. 3D graph between temperature and methanol/oil ratio showing their impact on response is shown in (Figure 7b), which shows that the highest response was observed at 6:1 methanol/oil ratio and 35 • C reaction temperature. However, high temperature can denature the active sites of the enzyme. Response surface plot for methanol-to-oil ratio and water content is presented in Figure 7c, which predicts that the combination of these two variables also affect the response. Figure 7d shows the combined effect of temperature and nano-biocatalyst concentration on the response, which indicates that the yield rises with rise in reaction temperature and nano-biocatalyst amount until it reaches an optimum point, i.e., 35 • C and 10%, respectively, while beyond this point a decline in % biodiesel yield is observed.
Catalysts 2020, 10, x FOR PEER REVIEW 9 of 16 terreus Lipase and CH3OH/oil ratio indicated that the % biodiesel increases with an increase in nanobiocatalyst level and methanol/oil ratio, and optimal biodiesel yield is obtained when nanobiocatalyst level and methanol/oil ratio were 10% and 6:1, respectively. Deviation from these values, however, results in lower biodiesel yield. 3D graph between temperature and methanol/oil ratio showing their impact on response is shown in (Figure 7b), which shows that the highest response was observed at 6:1 methanol/oil ratio and 35 °C reaction temperature. However, high temperature can denature the active sites of the enzyme. Response surface plot for methanol-to-oil ratio and water content is presented in Figure 7c, which predicts that the combination of these two variables also affect the response. Figure 7d shows the combined effect of temperature and nano-biocatalyst concentration on the response, which indicates that the yield rises with rise in reaction temperature and nano-biocatalyst amount until it reaches an optimum point, i.e., 35 °C and 10%, respectively, while beyond this point a decline in % biodiesel yield is observed. Figure 7. 3D graphical plots describing the impact of significant first-order interactions between (a) enzyme concentration × methanol to oil ratio (b) reaction temperature × methanol to oil ratio (c) water content × methanol to oil ratio (d) reaction temperature × enzyme concentration on biodiesel yield nano-biocatalyst reactions.

FTIR Spectroscopic Analysis of Eruca Sativa Oil and Biodiesel
The transesterification of E. sativa oil and biodiesel was monitored by FTIR spectroscopy ( Figure  8). The IR absorption band at 1438.0935 cm −1 corresponds to the -CH3 asymmetric bending that was present in E. sativa oil-based biodiesel but absent in E. sativa oil FTIR spectra. Similarly, absorption band at 1196.7202 cm −1 ascribed to O-CH2 stretching which was absent in E. sativa oil but present in E. sativa biodiesel spectra. The bands at 1161.1 cm −1 , 1459.3 cm −1, and 1099.6 cm −1, corresponding to C-O stretching, C=C stretching, and OCH2C asymmetric stretching, respectively, were absent in E. sativa oil biodiesel but present in E. sativa oil FTIR spectra. The bands at 1700-1800 cm −1 and 2800-3000 cm −1 (vibrational frequency band of C=O stretching and CH2 symmetric stretching, respectively) were present in both E. sativa oil and E. sativa oil biodiesel spectra. Tariq et al. reported the conversion of Figure 7. 3D graphical plots describing the impact of significant first-order interactions between (a) enzyme concentration × methanol to oil ratio (b) reaction temperature × methanol to oil ratio (c) water content × methanol to oil ratio (d) reaction temperature × enzyme concentration on biodiesel yield nano-biocatalyst reactions.

FTIR Spectroscopic Analysis of Eruca Sativa Oil and Biodiesel
The transesterification of E. sativa oil and biodiesel was monitored by FTIR spectroscopy (Figure 8). The IR absorption band at 1438.0935 cm −1 corresponds to the -CH 3 asymmetric bending that was present in E. sativa oil-based biodiesel but absent in E. sativa oil FTIR spectra. Similarly, absorption band at 1196.7202 cm −1 ascribed to O-CH 2 stretching which was absent in E. sativa oil but present in E. sativa biodiesel spectra. The bands at 1161.1 cm −1 , 1459.3 cm −1, and 1099.6 cm −1, corresponding to C-O stretching, C=C stretching, and OCH 2 C asymmetric stretching, respectively, were absent in E. sativa oil biodiesel but present in E. sativa oil FTIR spectra. The bands at 1700-1800 cm −1 and 2800-3000 cm −1 (vibrational frequency band of C=O stretching and CH 2 symmetric stretching, respectively) were present in both E. sativa oil and E. sativa oil biodiesel spectra. Tariq et al. reported the conversion of E. sativa oil into biodiesel by disappearance of peaks in E. sativa oil at 1465 and 1095 cm −1 and formation of new peaks at 1435 and 1195 cm −1 in the biodiesel [33]. Similar results have been reported in other research [34].
Catalysts 2020, 10, x FOR PEER REVIEW 10 of 16 E. sativa oil into biodiesel by disappearance of peaks in E. sativa oil at 1465 and 1095 cm −1 and formation of new peaks at 1435 and 1195 cm −1 in the biodiesel [33]. Similar results have been reported in other research [34].

Recovery and Reusability of Nano-Biocatalyst
Centrifugation was used to recover the nano-biocatalyst, which was then examined for lipase activity. It was revealed that after first use, there was no change in the CeO2@PDA@A. terreus Lipase activity, i.e., 16.90 U/mg/min (Table 7). So, recovered CeO2@PDA@A. terreus Lipase was reused multiple times for the production of biodiesel and after each reuse, the activity of lipase was assayed, which showed that a considerable decrease in the activity of lipase was observed after five uses, and after seven uses the activity was decreased to 3.9 U/mg/min. The immobilization of lipase is important for the economy of the biodiesel production process, as was observed/revealed by the results. The nano-biocatalyst can be used up to five times without significant decrease in activity, while free enzymes are readily deactivated at high temperatures, with little variations in pH and in the presence of short-chain alcohols, especially methanol [20]. So nano-immobilized lipases can appreciably help in cost reduction of biodiesel production process. After five uses, the biodiesel

Recovery and Reusability of Nano-Biocatalyst
Centrifugation was used to recover the nano-biocatalyst, which was then examined for lipase activity. It was revealed that after first use, there was no change in the CeO 2 @PDA@A. terreus Lipase activity, i.e., 16.90 U/mg/min (Table 7). So, recovered CeO 2 @PDA@A. terreus Lipase was reused multiple times for the production of biodiesel and after each reuse, the activity of lipase was assayed, which showed that a considerable decrease in the activity of lipase was observed after five uses, and after seven uses the activity was decreased to 3.9 U/mg/min. The immobilization of lipase is important for the economy of the biodiesel production process, as was observed/revealed by the results. The nano-biocatalyst can be used up to five times without significant decrease in activity, while free enzymes are readily deactivated at high temperatures, with little variations in pH and in the presence of short-chain alcohols, especially methanol [20]. So nano-immobilized lipases can appreciably help in cost reduction of biodiesel production process. After five uses, the biodiesel synthesis rate was reduced, which might have been due to the contact of nano-biocatalyst to organic compounds present in the reaction mixture during biodiesel production or recurring contact with heat. Similar studies have also been reported by Dumri and Hung [14].

Fuel Characteristics
Biodiesel fuel properties/characteristics were evaluated based on standard methods. The estimated fuel properties values of E. sativa-based biodiesel are presented in (Table 8).

Materials and Methods
Research/analytical grade chemicals were used in the whole experimental work. CeNO 3 .6H 2 O, NaOH, ethanol, lipase, methanol, acetone, polydopamine, phosphate buffer, tris-HCl buffer, and phenolphthalein were obtained from Sigma Chemical Co. (St. Louis, MO, USA). Eruca sativa seeds were obtained from Directorate of Land Reclamation Agriculture Department, Lahore, Pakistan.

Preparation of CeO 2 Nanorods
CeO 2 nanorods were prepared by using the hydrothermal method, during which 7.2 g of NaOH was dissolved in 20.0 mL distilled H 2 O along with stirring in a beaker. In another beaker, 2.17 g of CeNO 3 ·6H 2 O was dissolved in 10 mL distilled H 2 O. Then, the aqueous solution of CeNO 3 ·6H 2 O was added dropwise in the above solution of NaOH with continuous stirring. After 10 min, white precipitates were obtained. Then, the mixture was filled in 60 mL Teflon lined autoclave and kept in oven for 24 h at 120 • C. Light-yellow-colored product was acquired that was washed with distilled H 2 O and ethanol. For calcination, the dried product was also kept in furnace for 2 h at 500 • C [37].

Coating of Dopamine on CeO 2 Nanorods
For dopamine coating on CeO 2 nanorods, 0.3 g of prepared CeO 2 nanorods was dispersed in 20 mL of distilled H 2 O in a flask, then 20 mL of 20 mM tris buffer having pH 8.5 was added in the flask. 0.1 g of polydopamine was added in the above mixture and stirred for 1 h. The obtained suspension was separated by centrifugation and washed with the tris-HCl buffer to remove the unwanted polydopamine. Addition of polydopamine on nanorods was due to the polymerization of dopamine in the basic conditions [38].

Characterization of CeO 2 Nanorods
CeO 2 nanorods were characterized by means of the SEM, XRD, FT-IR, and EDX spectroscopy. XRD pattern for CeO 2 nanorods was obtained by using the X'pertpro (PANalyatical) with radiations of Cu k alpha of wavelength 1.54 Å in the scan range of 20-80 • , 2θ with scan step size of 0.02. XRD gave the information about the particle size, crystal phase, and dimensions of CeO 2 nanorods. The particle size and surface morphology of CeO 2 nanorods were confirmed by using the scanning electron microscope (JSM5910, JEOL, Tokyo, Japan) with 30 kV energy, (300,000×) maximum magnification, and 2.3 nm resolving power. The formation of CeO 2 nanorods and their modification with polydopamine were performed by using the Cary 630 Agilent FT-IR spectroscopy. EDX gives information about the purity, composition, and elemental analysis of the CeO 2 nanorods. EDX spectroscopic analysis was performed using EDX (JSM5910) (INCA200/Oxford instruments, High Wycombe, UK).

Immobilization of Lipase on Modified Nanorods
0.4 g of lipase was added in the 40 mL of phosphate buffer. The mixture of polydopamine-coated nanorods was slowly added in the lipase mixture with continuous stirring for 3 h at 4 • C. The resulting nano-biocatalyst (CeO 2 @PDA@A. terreus Lipase) was washed with the phosphate buffer several times to remove the un-reacted lipase and dried in desiccator at low temperature [39]. Figure 9 shows the lipase immobilization process. The immobilization efficiency was revealed by investigating the protein content in solution (before and after immobilization) using Bradford's method.

Characterization of CeO2 Nanorods
CeO2 nanorods were characterized by means of the SEM, XRD, FT-IR, and EDX spectroscopy. XRD pattern for CeO2 nanorods was obtained by using the X'pertpro (PANalyatical) with radiations of Cu k alpha of wavelength 1.54 Å in the scan range of 20-80°, 2θ with scan step size of 0.02. XRD gave the information about the particle size, crystal phase, and dimensions of CeO2 nanorods. The particle size and surface morphology of CeO2 nanorods were confirmed by using the scanning electron microscope (JSM5910, JEOL, Tokyo, Japan) with 30 kV energy, (300,000×) maximum magnification, and 2.3 nm resolving power. The formation of CeO2 nanorods and their modification with polydopamine were performed by using the Cary 630 Agilent FT-IR spectroscopy. EDX gives information about the purity, composition, and elemental analysis of the CeO2 nanorods. EDX spectroscopic analysis was performed using EDX (JSM5910) (INCA200/Oxford instruments, High Wycombe, UK).

Immobilization of Lipase on Modified Nanorods
0.4 g of lipase was added in the 40 mL of phosphate buffer. The mixture of polydopamine-coated nanorods was slowly added in the lipase mixture with continuous stirring for 3 h at 4 °C. The resulting nano-biocatalyst (CeO2@PDA@A. terreus Lipase) was washed with the phosphate buffer several times to remove the un-reacted lipase and dried in desiccator at low temperature [39]. Figure  9 shows the lipase immobilization process. The immobilization efficiency was revealed by investigating the protein content in solution (before and after immobilization) using Bradford's method. Figure 9. Schematic presentation of lipase immobilization process. Figure 9. Schematic presentation of lipase immobilization process.

Lipase Activity Assay in Free and Immobilized Form
Titrimetric method was performed for free and immobilized lipase activity assay [40]. The assay mixtures consisted of specific amounts of nano-biocatalysts, along with 10 mL of homogeneous mixture of olive oil in gum acacia, 5 mL of phosphate buffer (pH 7), and 2.0 mL of CaCl 2 (0.6%); then, the mixture was incubated at 37 • C for 1 h. Afterward, the reaction was stopped using 20 mL of ethanol:acetone (1:1) followed by the titration with 0.1 N NaOH solution. The lipase activity assay was executed by the titration of fatty acids produced from olive oil after reaction with the lipase.
One unit of lipase activity was defined as "the amount of enzyme which released one micro mole (µmol) of fatty acid per min. under specified assay conditions".
Lipase units were determined as follows: where ∆V = V 2 − V 1 ; V 1 = vol. of NaOH used against control flask; V 2 = vol. of NaOH used against experimental flask; N = normality of NaOH; M(sample) = mass of enzyme extract; and 60 = time of incubation (min) for bacterial lipase.

Effect of pH and Temperature on Activity of Free and Immobilized Enzyme
The effects of temperature and pH were studied on the activity of free lipase (A. terreus Lipase) and immobilized lipase (CeO 2 @PDA@A. terreus Lipase). The impact of temperature in the range of 25 to 50 • C and the pH in the range of 5 to 10 (using phosphate buffer) were investigated. Reactions were performed in triplicates.

Collection of Feedstock, Extraction of Oil, and Its Quality Assessment
The Eruca sativa seeds were subjected to solvent extraction based on Soxhlet technique for the extraction of oil using hexane as solvent. The extracted E. sativa seed oil was analyzed for its initial quality check by evaluating its iodine value, peroxide value, viscosity, acid value, specific gravity, and saponification value by means of AOCS standard protocols [35].

Central Composite Response Surface Methodology (CCRSM) Experimental Design
CCRSM was used for optimized biodiesel production from CeO 2 @PDA@A. terreus Lipase catalyzed transesterification of E. sativa oil. Five independent variables-A (methanol/oil ratio), B (biocatalyst amount), C (reaction temperature), D (reaction period), and E (water content)-were optimized within the ranges 3:1-9:1, 1-10%, 20 • C-50 • C, 12-48 h, and 0.2-1%, respectively. Fifty reactions were performed as per CCRD experimental design. In a typical biodiesel production reaction, a conical flask was used containing the reaction mixture oil, methanol, bio-catalyst, and distilled H 2 O. The reaction conditions were set according to the RSM experimental design. When the reactions were completed, glycerol was separated from biodiesel. Biodiesel was washed with warm water and using the rotary evaporator the residual methanol was recovered under the reduced-pressure conditions. The model used for response surface methodology is given below: 3.9. Recovery and Recycling of CeO 2 @PDA@A. Terreus Lipase After completion of transesterification reactions, the nano-biocatalyst (CeO 2 @PDA@A. terreus Lipase) was recovered by centrifugation of glycerol and methyl esters. The recovered nano-biocatalyst was washed with water, dried in air, and reused for transesterification to produce biodiesel [41].

Characterization of Biodiesel
The resulted biodiesel was characterized by FTIR spectroscopic and gas chromatographic techniques. Fuel properties were also determined using standard ASTM methods [42]. Fourier transform infrared spectroscopic analysis was performed by using the Cary 630 Agilent FTIR spectrometer for the E. sativa oil and biodiesel in the range of 400-4000 cm −1 .
Gas chromatography-mass spectrometry (GC-MS) analysis was carried out for the estimation of fatty acids methyl ester in the synthesized biodiesel. For this purpose, GCMS QP 2010 instrument dB 5 column of diameter 0.15 mm was used. Sample size and split ratio were 1 µL and 1:100, respectively. Helium was the carrier gas, and 1.2 mL/min flow rate was selected for sample elution. Oven temperature was set between 150 to 250 • C, at rate of 4 • C per minute. MS mass scanning range was 30-550 m/z. NIST MS library of GCMS was used to identify the alkyl esters. Biodiesel produced was further analyzed according to ASTM D methods to estimate the fuel properties to check their suitability as fuel.
Gas chromatography was performed using polar BPX_70 capillary column 30 m × 0.25 mm and FID detector, to determine biodiesel yield/ester content. Helium at 1.5 mL/min flow rate was used as carrier gas. The column oven temperature was kept at 100 • C and then increased to 260 • C, at rate of 10 • C/min. 1 g sample was taken in hexane having methyl heptadecanoate as internal standard, and 1 µL of this mixture was injected in column. FAME (%) yield was determined by using the following formula: where A ME denotes the sum of peak areas of all FAMAs. A is the peak area of methyl heptadecanoate, C is the concentration of the internal standard, V is the volume of the internal standard, and M is mass of biodiesel.

Conclusions
In the current study, CeO 2 @PDA@A. terreus Lipase was utilized for conversion of E. sativa oil into methyl esters and the optimization of process through response surface methodology was carried out. 89.3% FAME yield was obtained by carrying out reactions under optimum conditions of 10%CeO 2 @PDA@A. terreus Lipase, 6:1 methanol/oil ratio, and 0.6% water at 35 • C for 30 h reaction time. The fuel properties of synthesized product (biodiesel) were compatible with ASTM standards. CeO 2 @PDA@A. terreus Lipase was ascertained to exhibit efficient catalytic characteristics. Hence, CeO 2 @PDA@A. terreus Lipase may further be explored for the transesterification of other feedstock to attain a maximum yield of environmentally friendly fuel.