2.4. Fourier Transform Infrared Spectroscopic (FTIR) Analysis of Naked and Polydopamine Coated CeO2 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].
2.5. 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 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.
2.6. 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].
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].
2.8. 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.
2.9. 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.
2.10. ANOVA for Biodiesel Yield Response
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 (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 B2, C2, and D2 were the most significant quadratic terms, with p-values of 0.0001, 0.0001, and 0.0240, respectively, while A2 and E2 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
2SO
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.
2.11. 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
3OH/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.
2.12. 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
2C 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].
2.13. Major Fatty Acid Methyl Esters of Synthesized Biodiesel
GC-MS analysis for the profiling of major fatty acids methyl esters in product (biodiesel) is presented in (
Table 6). Palmitic acid (1.448% composition), oleic acid (28.181% composition), stearic acid (0.186%), gondoic acid (4.712%), and erucic acid (65.111%) were the major fatty acid methyl esters identified in sample.
Mumtaz et al. reported the fatty acids profile of
E. sativa oil that consists of palmitic acid (2.8%), linoleic acid (10.3%), stearic acid (0.90%), oleic acid (16.3%), linolenic acid (12.56%), and erucic acid (47.7%) [
35]. Chakrabarti and Ahmad also investigated and reported palmitic acid (10.2%), stearic acid (1.6%), oleic acid (22.8%), linoleic acid (6.4%), linolenic acid (11.9%), and erucic acid (40.8%) as major fatty acids in
E. sativa oil [
36].
2.14. 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].