Propagation characteristics of M. caribbica isolate MJTm3
To achieve a maximal yeast concentration for effective conversion of the substrate to ethanol, the yeast cells have been propagated through five phases that each contain varying concentrations of molasses (viz., 8, 10, 12, 14, and 16 oBrix). Overall, the findings showed that the yeast cell number increased exponentially up to 12 oBrix, but started to decline sharply after 14 oBrix. The results revealed that at 8 oBrix, 10 oBrix, and 12 oBrix, the yeast cells population increased by 34, 32, and 8%, respectively after 24 h propagation time (Fig. 3).
The pH value showed a decreasing trend when the propagation period was extended to 24 h. Specifically, the pH demonstrated increasing dynamics from 5.5 to 3.5, 3.6, and 3.5, at 8 oBrix, 10 oBrix, and 12 oBrix, respectively after a 24 h propagation time (Fig. 3). Furthermore, the pH value of the propagated broth increased to 3.2 as propagation time exceeded 24 h, making it more acidic. As a result, the number of viable yeast cells started to decrease. Although the molasses propagated media (MPM) has an acidic pH, the yeast's cellular morphology is maintained (Fig. 1a and b). This is consistent with our earlier findings, which stated that M. caribbica MJTm3 was found to be tolerant to acidic pH in YPD broth (Hawaz et al., 2022).
Regarding the reduction of molasses concentration in function to cell viability and incubation period, the Brix was decreased to 4.25 oBrix, 4.5 oBrix, 5.5 oBrix, and 6.25 oBrix from 8 oBrix, 10 oBrix, 12 oBrix, and 14 oBrix, respectively, during a 24 h propagation period. Consequently, viable yeast cell values of a density of 4.26 x 108 cells/mL, 8.75 x 108 cells/mL, 9.52 x 108 cells/mL, and 2.1 x 107 cells/mL were produced (Fig. 3). The propagated yeast density under different molasses concentration can be seen in Fig. 2. Even though changes in pH, and cell density were observed until the end of the propagation time, the molasses concentration (Brix content) stabilized after 20 h propagation time.
Model Validation And Optimization Of Fermentation Parameters
The aim of this study was to maximize the yield of bioethanol produced from sugarcane molasses by optimizing the conditions of the fermentation parameters essential for efficiently converting all of the available fermentable sugar to ethanol. The complete design matrix of the independent variables in actual values corresponding with predicted responses of the bioethanol yield is presented in the supplementary data. A second-order quadratic model equation was generated based on the experimental data and CCD and this indicated linear, interaction and quadratic effect of variables on bioethanol yield as (+ ve) or (-ve) under Eq. 1. The quadratic equations had high regression coefficients, and the lack of fit values was insignificant (p > 0.05), demonstrating that the quadratic models fit the data well.
Where Y is the bioethanol yield (%) and a positive sign indicates a synergetic effect, whereas a negative sign indicates an antagonistic effect.
The statistical significance of the quadratic model was determined using analysis of variance (ANOVA). The statistical significance was controlled by F-test and P-values and the model was found to be highly significant at a 95% confidence level (Eq. 5), with an F-value of 65.48 and a very low probability (p < 0.0001) (Table 2). This indicated that there is less than 0.01% chance that this error is caused by noise. The significance of the current model agrees with results of Hamouda et al. (2015) who developed a model for optimization ethanol production parameters with statistically significance at a 95% confidence level, an F-value of 29.1, and a very low probability of p < 0.0001 with < 0.01% chance that of which error caused by noise.
The model fitting reliability was evaluated using R2 and adjusted-R2 and found to be 0.9682 and 0.9534, respectively (Table 2). These values indicated that approximately 96.82% of the variability in the response obtained is explained by the model ensuring perfect fit of the model to the observed data. This is in line with results from Hamouda et al. (2015), the R2 value 0.953, which demonstrated the high significance of the model with 95.3% of confidence for bioethanol production from sugarcane molasses. Cavalaglio et al. (2016) have developed a model fitted with R2 equal to 0.970 with 97% significance for optimization of ethanol production from a cellulosic substrate. Flayeh (2017) has designed a model that fitted with a slightly lower confidence of 90.25% compared to the present reported values.
The regression model's suitability between the experimental and predicted data of the response parameter suggests a reasonable correlation over the tested experimental ranges (Table 2). According to the current analysis, a ratio of 38.42 adequate signals (Eq. 4) was achieved for the CCD consisting of 86 trials, for which the overall experimental bioethanol yield ranged between 8.41 and 46.59% and corresponding predicted value of 8.28 to 39.98%, respectively. An adequate signal-to-noise ratio of greater than four (R > 4) and an adequate accuracy value of 38.42 were found using ANOVA measurement. The current model can therefore be reliable and employed to navigate the design space. This is in agreement with Hamouda et al. (2015), who obtained an experimental ethanol yield of 8.20–41.4% with the corresponding predicted values ranging from 9.26–39.1%, respectively with an adequate signal of 17.1. On the other hand, the current developed model indicated that the lack of fit (F = 0.0525) was found significant relative to the pure error.
The actual and predicted values, as well as the normal plot of experimental design residuals, which are shown in Figs. 5 (a and b), further supported the aforementioned ANOVA analysis. A plot of the predicted and experimental values of the bioethanol yield is shown in Fig. 5 (a). The plot showed a strong correlation (R2 = 0.9682) between the experimental and predicted data, demonstrating that the model accurately predicted the bioethanol production within the experimental range under consideration. This demonstrates that the experimental outcomes were largely consistent. The residuals for the bioethanol production are normally distributed on a normal plot, as shown in Fig. 5 (b), with results extremely closely spaced to a straight line with no substantial departure.
In the current study, the temperature (A), pH (B), inoculum size (C), molasses concentration (D), mixing rate (E), and incubation period (F) were selected as key factors to maximize the bioethanol yield (%) using CCD. Results revealed that all linear and interaction factors, except for temperature (A), inoculum size (C), pH and molasses concentration (BD), inoculum size and molasses concentration (CD), inoculum size and mixing rate (CE), and mixing rate and incubation period (EF), were significant at the 95% confidence level. The independent factors, such as linear (B, D, E, F), interactive (AB, AC, AD, AE, AF, BC, BE, BF, CF, DE, and DF), and quadratic (A2) were found significant model terms (p < 0.05), whereas all quadratic variables were showed insignificant (p > 0.05), except for temperature (A2) (p < 0.0001) (Table 2). This is in agreement with Kamal et al. (2021) who reported that linear pH, interactive inoculum size, and pH were significant model terms for efficient bioethanol production from sugarcane molasses.
Table 2
Analysis of variance for response surface of quadratic model for the CCD experiments
Source
|
SS
|
Df
|
MS
|
F-value
|
P-value
|
|
Model
|
4731.99
|
27
|
175.26
|
65.48
|
< 0.0001
|
Significant
|
A
|
13.39
|
1
|
13.39
|
5.00
|
0.0292
|
|
B
|
30.96
|
1
|
30.96
|
11.57
|
0.0012
|
|
C
|
15.29
|
1
|
15.29
|
5.71
|
0.0201
|
|
D
|
802.05
|
1
|
802.05
|
299.65
|
< 0.0001
|
|
E
|
1101.06
|
1
|
1101.06
|
411.36
|
< 0.0001
|
|
F
|
422.51
|
1
|
422.51
|
157.85
|
< 0.0001
|
|
AB
|
230.10
|
1
|
230.10
|
85.97
|
< 0.0001
|
|
AC
|
866.95
|
1
|
866.95
|
323.90
|
< 0.0001
|
|
AD
|
32.85
|
1
|
32.85
|
12.27
|
0.0009
|
|
AE
|
49.76
|
1
|
49.76
|
18.59
|
< 0.0001
|
|
AF
|
81.03
|
1
|
81.03
|
30.27
|
< 0.0001
|
|
BC
|
178.25
|
1
|
178.25
|
66.60
|
< 0.0001
|
|
BD
|
9.86
|
1
|
9.86
|
3.68
|
0.0599
|
|
BE
|
59.31
|
1
|
59.31
|
22.16
|
< 0.0001
|
|
BF
|
79.90
|
1
|
79.90
|
29.85
|
< 0.0001
|
|
CD
|
2.51
|
1
|
2.51
|
0.9383
|
0.3367
|
|
CE
|
10.70
|
1
|
10.70
|
4.00
|
0.0503
|
|
CF
|
328.57
|
1
|
328.57
|
122.76
|
< 0.0001
|
|
DE
|
150.18
|
1
|
150.18
|
56.11
|
< 0.0001
|
|
DF
|
29.31
|
1
|
29.31
|
10.95
|
0.0016
|
|
EF
|
13.12
|
1
|
13.12
|
4.90
|
0.0308
|
|
A²
|
190.47
|
1
|
190.47
|
71.16
|
< 0.0001
|
|
B²
|
2.79
|
1
|
2.79
|
1.04
|
0.3117
|
|
C²
|
11.11
|
1
|
11.11
|
4.15
|
0.0462
|
|
D²
|
16.45
|
1
|
16.45
|
6.15
|
0.0161
|
|
E²
|
12.50
|
1
|
12.50
|
4.67
|
0.0348
|
|
F²
|
21.44
|
1
|
21.44
|
8.01
|
0.0064
|
|
Residual
|
155.24
|
58
|
2.68
|
|
|
|
Lack of fit
|
145.56
|
49
|
2.97
|
2.76
|
0.0525
|
Not significant
|
Pure error
|
9.68
|
9
|
1.08
|
|
|
|
Corrected total
|
4887.23
|
85
|
|
|
|
|
R-Squared (R2) = 0.9682; Adjusted R2 = 0.9534; Adeq Precision = 38.42
|
Note: SS: sum of squares; DF: degree of freedom; MS: mean square. |
The Pareto chart was plotted to highlight the most significant independent variables and their main and interaction effects on bioethanol production (Fig. 4). The effect of each independent parameter on bioethanol production was confirmed by the coefficient of the quadratic equation (Eq. 1.). As a result, linearly; mixing rate (rpm) and incubation period (h) were considered significant variables (p < 0.0001). The mixing rate revealed the highest positive impact on the bioethanol yield (%), followed by the incubation period. This suggests that increasing the mixing rate and incubation time beyond their preset values will increase the bioethanol yield. Other fermentation factors, i.e., inoculum size, initial pH, and temperature, showed a slightly positive impact on the bioethanol production (p = 0.0201, 0.0012, and 0.0292, respectively). In the present study, molasses concentration exerted an adverse negative impact (p < 0.0001) on the ethanol production process, while its quadratic effect demonstrated the highest positive impact on ethanol fermentation (p = 0.0161). The quadratic effects of the incubation period, initial molasses concentration, mixing rate, and inoculum size showed a positive impact on the fermentation process in declining order (p = 0.0064, 0.0161, 0.0348, and 0.0462, respectively). On the other hand, the incubation temperature showed the highest negative quadratic effect (p < 0.0001) on the bioethanol yield.
The positive interactive effect of all the parameters within the studied range of the experiment ranked in declining order: inoculum size and mixing rate (p < 0.0001) > temperature and initial pH (p < 0.0001) > initial pH and inoculum size (p < 0.0001) > initial pH and incubation period (p < 0.0001) > temperature and mixing rate (p < 0.0001) > temperature and initial molasses concentration (p = 0.0009) > inoculum size and mixing rate (p = 0.0503) > initial pH and initial molasses concentration (p = 0.599) > inoculum size and initial molasses concentration (p = 0.3367). The negative interactive effect of the investigated parameters ranked in the following decreasing order: temperature and inoculum size (p < 0.0001) > initial molasses concentration and mixing rate (p < 0.0001) > temperature and incubation period (p < 0.0001) > initial pH and mixing rate (p < 0.0001) > initial molasses concentration and incubation period (p = 0.0016) > mixing rate and incubation period (p = 0.0308).
Surface Plot Analysis
An extensive experimental trial was conducted over the considered range of 86 experimental runs to determine the optimum operating value for the factors that will maximize the bioethanol yield from sugarcane molasses by M. caribbica MJTm3. To elucidate the optimum condition of each factor for a maximum bioethanol yield (%) production, thee-dimensional response surfaces (3D) was plotted based on the predicted second-degree polynomial equation (Fig. 6).
According to the response surface plots for bioethanol yield as a function of temperature and molasses concentration in Fig. 6 (a), lowering the temperature from 35-30.25oC and the molasses concentration from 35-25.13 oBrix resulted in an increase of the bioethanol yield. At 34.90 oBrix and 25.05oC a lower bioethanol production of 17.40% was observed, but a higher bioethanol yield of 28.70% was shown at 25.06 oBrix of the molasses and 29.78oC. There is a negative correlation between molasses concentration and bioethanol yield, meaning that a decrease in molasses concentration will result in higher bioethanol yield. The yield began to decrease at 25.09 oBrix and 29.96oC, which is because high temperatures and high substrate concentration have a negative effect on the fermentation capacity of the yeast cells (Cazetta et al., 2007). Low bioethanol yield signified that the yeast was sensitive to the inhibitory compound present in the fermentation medium. Moreover, the enzymes that regulate the fermentation process are sensitive to high temperatures which can result in the denaturation of their tertiary structure (Phisalaphong et al., 2006). According to Liu and Shen (2008), the optimal operating temperature for free fermenting yeast cells was near 30°C. It has been shown that increasing the initial molasses concentration to 25 oBrix resulted in a maximum bioethanol yield of 25% (Hamouda et al., 2015). In agreement with our finding, Morimura et al. (1997) reported the highest ethanol yield and ethanol concentration at an optimum of 25 oBrix molasses concentration by a yeast strain K211 of Saccharomyces species.
A positive interacting effect between temperature and mixing rate was seen for the bioethanol yield Fig. 6 (b). The bioethanol yield increased with an increase of the mixing rate from 111-149.77 rpm and temperatures ranging from 25-30oC. Results showed that at 110 rpm and 25°C, a bioethanol production of 18.20% was obtained. At a mixing rate of 149.77 rpm and a temperature of 30oC, a maximum yield of 29.26% of bioethanol was observed. When the mixing speed and temperature reached 30.2oC and 149.83 rpm, respectively, the bioethanol production started to decline. It is obvious that as the agitation rate increases, the diffusion of the necessary nutrient from the fermentation broth to the yeast cells is increased. However, at the same time, this also increased the release of ethanol from the cells to the fermentation broth. The ideal agitation rate for ethanol fermentation by yeast cells was reported as 150–200 rpm (Zabed et al., 2014). In the current situation, the fermentation broth was efficiently mixed and distributed when the agitation rate increased ≈ 150 rpm.
Figure 6 (c) demonstrates the negative interactive effect of molasses concentration and mixing rate for bioethanol production from molasses. The bioethanol yield increased with an increase in mixing rate from 110.19-149.67 rpm and a decrease of the molasses concentration from 35-25.18 oBrix. This finding revealed that the maximum bioethanol yield of 34.84% was demonstrated at 25.19 oBrix and 149.67 rpm, which is because the yeast cells are in contact with vital nutrients, like sugars, which resulted in an effective bioethanol yield (Kopsahelis et al., 2007). Moreover, at 25.25 oBrix and 147.32 rpm, the bioethanol output started to decrease.
Figure 6 (d) shows the negative interactive effect of temperature and incubation period on bioethanol production. Results showed that the yield increased with a decrease of the incubation temperature from 35-29.86oC, and a slight decrease in the incubation period from 72-71.78 h. At 29.86oC and a 71.78 h incubation period, a maximum yield of 27.85% ethanol was produced. However, at 25.05°C and 48.19 h, the lowest bioethanol yield of 18.04% was measured. An earlier study using Saccharomyces species indicated that fermentation at 31°C for 72 h produced 26.4% of bioethanol yield (Flayeh, 2017). According to Zabed et al. (2014), a longer fermentation period is required to recover a high ethanol yield with the highest productivity using a batch fermentation system. This suggests that the use of a short fermentation time and a short incubation temperature cause inefficient ethanol fermentation due to inadequate growth of microorganisms (Zabed et al., 2014).
Figure 6 (e) illustrates the interaction effect of molasses concentration and incubation time on the bioethanol yield. The plot showed that ethanol production increased as the incubation period was extended from 48.16 to 71.99 h, while the concentration of molasses decreased from 35 to 25.10 oBrix. At 48.16 h of incubation time and 25.11 oBrix, a lower ethanol yield of 26.71% was obtained. A maximum bioethanol yield of 32.98% was recorded at 25.10 oBrix and 71.99 h of molasses concentration and incubation period, respectively. When the fermentation flask was overloaded with the substrate, a continuous fermentation rate caused the yeast cells to experience osmotic shock that has an inhibitory effect on the yeast (Azhar et al., 2017; Cavalaglio et al., 2016). The results revealed that the production of bioethanol started to decline immediately after 72 h and 26 oBrix of incubation time and molasses concentration, respectively.
A positive interactive effect of mixing rate and incubation time on the production of bioethanol from sugarcane molasses is demonstrated in Fig. 6 (f). The production of bioethanol improved by increasing the mixing rate from 110.43 to 149.94 rpm and the incubation time from 48.15 to 71.96 h. A lower bioethanol yield of 19.60% was shown at 48.15 h and 110.43 rpm. However, a maximum bioethanol production of 32.45% was produced at 149.94 rpm and 71.96 h of mixing rate and incubation time, respectively. The plot also demonstrated that the bioethanol yield started to decrease at 147.92 rpm and 70.89 h of mixing rate and incubation time, respectively. Furthermore, when the incubation period reached 72 h, the production of bioethanol started to decline. This might be due to ethanol oxidation, organic acid formation, and substrate depletion in the fermentation broth that could potentially accelerate the deactivation of enzymes and thereby lowering both ethanol yield and yeast cell viability (Kopsahelis et al., 2007). Low ethanol yield may also occur due to the formation of secondary byproducts that limit ethanol productivity (Ergun & Mutlu, 2000; Hamouda et al., 2015).
Fermentation Under Optimum Conditions
Ethanol fermentation parameters, such as pH, inoculum size, molasses concentration, temperature, mixing rate, and incubation period, were optimized in the above described experiments and applied to evaluate the reliability of the model equation using batch fermentation. The maximum predicted bioethanol concentration of 49 g− 1 L, bioethanol yield of 78.6% was obtained under the predicted optimal conditions of pH 5.5, 20% inoculum size, 25 oBrix initial molasses concentration, 30oC temperature, 72 h incubation period, and 159 rpm with the desirability of 1.0.
Fermentation was conducted using Erlenmeyer flask with a working volume of 2 L to validate the predicted optimal conditions in the actual experiment using the parameters specified above. Results showed that, during the fermentation process, the pH value slightly declined from 5.5 to 5.27 at 24 h. However, after 48 h of the incubation period, the pH returned back to the optimal condition (i.e., 5.5) without adjustment and was kept constant until completion of the ethanol fermentation process, which might be due to the production of enzymes and other chemicals required for adaptation to the new environment to facilitate the overall ethanol fermentation process. On the other hand, molasses might exhibit a buffering effect attributed to its acid composition (weak acids and amino acids) and phosphates that would regulate the pH values to 3–5 and 6–7, respectively (Cazetta et al., 2007).
In the present study, the key process variables were experimentally supported to produce a maximum ethanol concentration of 56 g.L− 1 with a bioethanol yield of 86% and a percent theoretical yield of 95% from the 25 obrix molasses concentration within 72 h Fig. 7 (d). Hamouda et al. (2015) obtained a bioethanol concentration and bioethanol yield of 32.32 g.L− 1 and 44%, respectively, from sugarcane molasses using Pichia veronae strain HC-22 during 60 h fermentation time. The bioethanol yield in our study was found higher due to the differences in fermentation period, inoculum size, fermenting microorganisms, and pretreatment methods. Over all, the current study suggests that the process optimization using RSM was applicable to optimize the bioethanol fermentation reliably from molasses by M. caribbica MJTm3.