Improved Sugar Recovery of Alkaline Pre-Treated Pineapple Leaf Fibres via Enzymatic Hydrolysis and Its Enzymatic Kinetics

: The bioﬁbre complex structure in pineapple leaf ﬁbres (PALFs) can be disrupted using the alkaline pre-treatment method with sodium hydroxide (NaOH). Nonetheless, the pre-treated structures have an impact on the bioconversion of PALFs into sugar. Thus, the impact necessitates precise reaction conditions, which are required for the enzymatic hydrolysis of PALFs. In this study, the Box–Behnken design (BBD) was utilised to achieve maximum sugar yield from PALFs. The optimised enzymatic hydrolysis conditions were 229 µ g/mL of cellulase, a reaction temperature of 45 ◦ C, and a hydrolysis time of 69 h. Resultantly, an 84% increase in sugar yield was observed, from 17.26 mg/mL to 108.74 mg/mL after the optimisation process. Moreover, a signiﬁcant inﬂuence ( p = 0.0009) of the enzymatic hydrolysis on hydrolysis time and temperature was observed. Additionally, the kinetic study analysis of NaOH pre-treatment and the enzymatic hydrolysis process were investigated. The kinetics of enzyme concentrations on total reducing sugar (TRS) production using an insoluble substrate were investigated based on modiﬁed Michaelis–Menten and Chrastil models. The result of the models is in good agreement with the experimental data, as the PALF hydrolysis process was reasonably well predicted. This study provides valuable information for predicting the sugar yield of the PALF bioconversion.


Introduction
From 2016 to 2020, Malaysia produced 1505.3 metric tonnes of pineapples, which makes pineapple leaf fibres (PALFs) among the country's most abundant agricultural waste in the growing pineapple industry [1]. Nevertheless, the abundance of pineapples produces residues, which results in landfilling and environmental issues. Thus, PALF was identified as a promising biomass source candidate for bioproduct conversion in the renewable energy sector. Typically, the bioconversion process of PALF to reduced sugar involves two main steps, which are the pre-treatment and hydrolysis processes. During the alkaline or acid pre-treatment process, the lignocellulosic biomass structure is disrupted, and the polysaccharides are converted to reducing sugar through the hydrolysis process. Although the acid hydrolysis process can be utilised for the pre-treatment step, the enzymatic hydrolysis process was recognised as a more promising approach in multiple biomass conversion processes into sugar [2]. The enzymatic hydrolysis can be completed with no corrosion issues under moderate conditions. Moreover, the process requires less energy, which offers a lower processing cost than acid hydrolysis [3]. Additionally, several studies have reported that enzymatic hydrolysis improves the bioconversion of lignocellulosic material into fermentable sugar by combining the pre-treatment process [4][5][6][7].
Generally, the enzymatic hydrolysis process is difficult to carry out due to the complex structure of lignocellulosic biomass. Therefore, several hydrolytic enzymes are usually required in performing the enzymatic hydrolysis [8,9]. To improve the sugar yield, the optimisation of the enzymatic hydrolysis can be achieved by using a structured procedure, such as the Box-Behnken Design (BBD) as part of the response surface methodology (RSM). The BBD is known as a higher-order regression model that does not contain a fractional factorial embedded within the design. Furthermore, the BBD only requires three levels of each factor. The BBD defines a mid-level between a low and high level of the factors, which reduces an extreme axial point. Moreover, the BBD is preferred over other designs due to the effectiveness of the experimental model as it is interpreted through the number of coefficients divided by the run numbers [10]. In terms of its application, the use of BBD has been applied in rice straw and wheat straw hydrolysis processes to estimate the optimum conditions for enzymatic reactions [6,7].
An important parameter in determining the optimal amount of reducing sugar in the hydrolysis reaction is enzyme loading. The amount of reducing sugar produced can be enhanced by applying a relatively high enzyme concentration. Despite that, the continuous increase in enzyme dosage does not always significantly increase sugar production. Qi et al. reported that when cellulase concentrations of 10 and 30 FPU/g substrates were used, the reducing sugar yield amounts were 72.65% and 90.94%, respectively [7]. Nonetheless, the continuous increase of the cellulase concentration to 50 FPU/g exhibits a plateau trend in the amount of sugar yield (90.07%). Hence, the enzyme loading process in the enzymatic hydrolysis reaction should be optimised to obtain the maximum yield of reducing sugar.
Another parameter that influences enzymatic activity during the hydrolysis process is temperature. One reported work reveals that the optimal temperature of cellulase to hydrolyse the substrate in an enzymatic hydrolysis process was between 40 and 50 • C [8]. According to several studies, the optimum temperature between 43 and 45 • C produces the highest reducing sugar yield in the enzymatic hydrolysis of lignocellulosic biomass [4,11]. Interestingly, in another study, the amount of the reducing sugar continued to increase when the temperature rose to 60 • C [10]. Thermal inactivation caused by an increase in temperature for a certain period can eventually influence the efficiency of enzymatic hydrolysis [12]. The temperatures ranging from 40 to 60 • C were identified as contributing to the high reducing sugar production. However, the production of reducing sugar also depends on the source of the enzymes, whether the enzyme is mesophilic or thermophilic.
Hydrolysis time also has a considerable effect on the enzymatic hydrolysis process. Particularly, the hydrolysis time allows substrates to be completely converted into their products. The two main aspects that influence the hydrolysis rate are enzymatic activity and substrate accessibility. In previous studies, the increment in hydrolysis yield was observed for wheat straw and narrow-leaf cattail at 24 h and 72 h, respectively. On the contrary, when the hydrolysis time was extended to 96 h, the yield amount decreased [7,13]. Other similar findings were also investigated in the relation to the effect of hydrolysis time. For example, the conversion of potato peel residues and sweet sorghum bagasse was at the optimum level after 48 h of reaction [4,14]. Meanwhile, the hydrolysis reactions of Cistus ladanifer and Cytisus striatus reported the highest reducing sugar yield after 72 h of reaction [10].
A high conversion rate is an important aspect that should be highlighted during the conversion of polysaccharides into monomeric sugar. Hence, the understanding of hydrolysis kinetics is essential in determining the mechanism of the bioconversion process. The bioconversion of polysaccharides may face a difficult issue if the hydrolysis reaction occurs heterogeneously between soluble enzymes and insoluble substrates. Factors that cause the modelling of the enzymatic hydrolysis process to be difficult include the activity of enzymes, enzyme-substrate interaction, and substrate complexity [15]. Thus, an improvised classical Michaelis-Menten (MM) model can be utilised. The MM model determines the initial reaction velocity (v 0 ) as a function of initial enzyme loading (E 0 ) rather than substrate concentration [16]. The enzymatic kinetic property of the heterogenous system has been studied by several researchers using the modified MM model in the bioconversion process of cellulose from corn stover, sugarcane, and wheat straw [17][18][19]. Moreover, the modified MM model should be applied if the enzymatic hydrolysis is a diffusion-limited type process as it depends on the heterogenous rate-limiting reaction of the substrateenzyme complex [18]. Additionally, the Chrastil model is a diffusion-limited type of process. The model is used to determine additional structural and kinetic characteristics of a heterogeneous reaction system. Particularly, the Chrastil kinetic model investigates the dependence of enzymatic reactions on the heterogenous rate-limiting factors such as adsorption and diffusion processes [15].
The optimisation of reducing sugar production from the hydrolysis of PALF hydrolysate is crucial because it has potential as a prebiotic for improving probiotic growth [20] and as a precursor for biofuel production via fermentation of pineapple leaf lignocellulosic hydrolysate [21]. In this study, the optimum operating conditions were determined for the bioconversion of the pre-treated PALFs into reducing sugar. The influence of enzyme concentration, temperature, and hydrolysis time on the reducing sugar yield was studied using BBD. The optimised conditions that were developed for the bioconversion of PALFs greatly improved the TRS yield. Then, the pre-treated PALF hydrolysis reaction kinetics were evaluated using different initial enzyme concentrations. The improvised models of classical MM and the Chrastil model proposed for the diffusion-limited enzymatic reaction were employed and evaluated to analyse the experimental data.

Sample Preparation Process
Pineapple leaves were obtained from the Pineapple Cannery of Malaysia Sdn. Bhd., Malaysia. The pineapple leaves were hard-pressed using a two-roll mill machine to eliminate the water content of around 90%. The PALFs were physically separated from the semi-dried leaves. Next, the PALFs were thoroughly rinsed with water. They were dried in an oven for 24 h at 105 • C until reaching a steady weight value. The PALF was cut into smaller sizes and ground using a mortar and pestle. Then, the sample was stored in a dry place for further analysis. The composition of the PALF was investigated using the technique proposed by Nashiruddin et al. [22]. According to the composition analysis, the contents of cellulose, hemicellulose, and lignin in the PALF were 72.76 ± 0.01, 17.15 ± 0.04, and 4.76 ± 0.02% (w/w). From these components, the monosaccharides present in the PALF biomass were glucose from the cellulose part, xylose, arabinose, and D-glucuronic acid from the hemicellulose part. The enzyme used for the enzymatic hydrolysis process was fungal cellulase from the Aspergillus species (C2605), which was purchased from Sigma-Aldrich (Merck, Darmstadt, Germany). All chemicals used in this study were of pure and analytical grade.

Alkaline Pre-Treatment Process
The PALFs were pre-treated using an alkaline pre-treatment method, according to our previous report [22]. The 1% (w/v) PALF samples were soaked in a 2.43% (w/v) sodium hydroxide (NaOH) solution in a 250-mL beaker. The samples were treated for 57.15 min at 87 • C with shaking at 100 rpm. Subsequently, they were neutralised and thoroughly rinsed with distilled water. The pH solution of the samples was measured using a pH metre (Mettler Toledo, Columbus, OH, USA) to confirm it reached a neutral pH of 7. The samples were then air-dried in a 60 • C oven until they reached their constant weight and were stored properly.

Optimisation of Enzymatic Hydrolysis Operating Parameters
The 10% (w/v) pre-treated PALFs were dissolved in a 0.05 M sodium acetate buffer (pH 4.8) containing various enzyme loading per g of PALF (15 to 25 µg/g). An incubator shaker was used to incubate the enzymatic mixture for 24 to 72 h. The process was carried out at 155 rpm using different temperatures (40 to 60 • C). Then, 5 mL were taken from each sample at 24-h intervals to determine the TRS content. The BBD from Design Expert Software 6.0.4 (Stat-Ease Inc., Minneapolis, MN, USA) was used to establish the optimum enzymatic hydrolysis conditions. The operating ranges for enzyme loading per g of PALF (15 to 25 µg/g), hydrolysis times (24 to 72 h), and temperatures (40 to 60 • C) were determined for the enzymatic hydrolysis optimisation process. Subsequently, the parameters were studied at low (−1), centre (0), and high (1) levels. By equating the reducing sugar yield as an output with independent variables (enzyme loading, temperature, and hydrolysis time), a statistical model can be obtained as follows Equation (1): where Y i is the dependent variable, β o , β i , β ii , and β ij are the regression coefficients, x i is the independent variable i, and x j is the independent variable j.
The experiments were performed for 17 runs, as shown in Table 1. The relationship between the coded values and actual values is described according to Equation (2): where x i is the independent variable of the coded value, X i is the independent variable actual of the value, X 0 i is the independent variable of the actual value at the centre point, and ∆X i is the step change value. The analysis of variance (ANOVA) and lack of fit (LOF) presented in Table 2 were used to estimate the statistical parameters. The regression analysis of the data obtained for studying TRS yield optimisation resulted in the following second-order polynomial model (Equation (3)):

Kinetics Study of Enzymatic Reactions
The enzymatic hydrolysis reaction was conducted using various concentrations of the enzyme, as reported by Carillo et al. [17]. The pre-treated PALF were hydrolysed in a 0.05 M sodium acetate buffer (pH 4.8) with different enzyme loadings per g of PALF (0.5, 15.0, 20.0, 22.9 (from enzymatic hydrolysis optimisation), 25.0, and 30.0 µg/g). They were incubated in an incubator shaker at 155 rpm for 72 h at 45 • C. The experiments were run in triplicate, and the hydrolysis rates were calculated from the product formation curve.

Quantification of Total Reducing Sugar (TRS)
The response variable from the optimisation of enzymatic hydrolysis was the TRS concentration. After PALF hydrolysis using cellulase, 3 mL of the sample reaction was taken out at different time intervals. The enzyme mixture was boiled for 10 min to deactivate the enzyme activity. The TRS amount produced from the enzymatic hydrolysis was determined using the 3,5-dinitrosalicylic acid (DNS) procedure [23,24] with glucose as the standard. The specific sugar detected from the hydrolysis of PALF in this study is glucose. An equal amount (1 mL) of sample reaction and DNS reagent was mixed and boiled for 5 min at boiling water (>100 • C). The absorption of the reaction mixture was determined at 575 nm using a spectrometer.

Kinetics Modelling: Estimation of the Parameters
A modified Michaelis-Menten model was applied in this study as we used PALF as a heterogenous substrate and cellulase as a soluble enzyme. The substrate turnover rate at maximum enzyme saturation (V emax ) and the half-saturation constant (K e ) were determined by performing a linear regression analysis on the experimental data using Equation (4): where v 0 is the initial hydrolysis velocity, V emax is the maximum velocity at the initial concentration of adsorption sites on the substrate, and K e is the half-saturation constant for that initial concentration.
The TRS amount at 60 min was chosen as an estimation value. The estimation value was carried out as v o is undefined when t→0 (division by zero) [17]. Non-linear regression analysis was used to define the model variables, which were derived from experimental production curves based on Equation (5): where P and P ∞ are the products at time t and at equilibrium, respectively. In Fick's law, k is the rate constant that changes with the diffusion coefficient, E 0 is the initial enzyme concentration, and n is a structural diffusion resistance constant. The parameter n represents the characteristics of the reaction order.
If the system has low diffusion resistance, then n tends to approach unity. On the other hand, when the system is constrained by diffusion resistance, n tends to be small (n = 0.5 to 0.6) [20].

Optimisation of the Enzymatic Hydrolysis Process
The parameters that influence reducing sugar production were analysed using BBD. The data obtained for enzymatic hydrolysis are summarised in Table 3. The enzyme loading was adjusted between 15.0 and 25.0 µg/g. The temperature was varied between 40 and 60 • C, while the hydrolysis time was varied from 24 to 72 h. The analysis of variance (ANOVA) was utilised to examine the adequacy of the developed model (Table 2).

ANOVA and Statistical Data Analysis
As shown in Table 3, the ANOVA analysis was performed to determine the influence of the variables and interactions among all the independent variables. The results of the ANOVA analysis consist of 17 cycles obtained from Table 2. The F value calculated for the model was 27.82, with a prob > F value of 0.0009. The results indicate that the model terms are statistically significant. Subsequently, the accuracy of the model was checked using the coefficient of determination, R 2 . The R 2 value calculated for the reducing sugar yield was 0.98, which represents 98% of the variability in the response. Therefore, the high value of R 2 reveals that the experimental and projected data are highly correlated. Additionally, the expected versus actual data for TRS was depicted in Next, the lack-of-fit (LOF) test was estimated to be insignificant to ensure that the developed model satisfactorily fits the data. As shown in Table 2, the results of the LOF test for all the responses were insignificant. The prob > F value obtained was 0.6057 (p > 0.05), which reveals the significance of the correlation between all the variables and responses. Table 2 presents the p values of each variable, which suggests that all three variables (enzyme loading (X 1 ), temperature (X 2 ), and time (X 3 )) produce a significant (p < 0.05) influence on TRS yield amount. Furthermore, the interaction of X 2 2 , X 1 X 2 , X 2 X 3 , and X 1 X 2 2 significantly affects the amount of TRS yield. From Equation (3), the X 1 X 2 , X 2 X 3 , X 1 2 X 2 , and X 1 X 2 2 interactions have negative values. Hence, the negative values suggest that a maximum TRS yield was observed at the centre point. Moreover, if one of the variables increases while the others decrease above the centre point, the reducing sugar yield amount will not be affected. The negative values observed in the regression coefficient for the square terms (X 1 2 , X 2 2 , X 3 2 ) indicate that the optimum response at the centre level decreases the TRS yield amount. The occurrence of this observation could be due to the increase or decrease of parameter levels from the centre points. Figure 1. variables (enzyme loading (X1), temperature (X2), and time (X3)) produce a significa < 0.05) influence on TRS yield amount. Furthermore, the interaction of X2 2 , X1 X2, X and X1X2 2 significantly affects the amount of TRS yield. From Equation (3), the X1X2, X2X3, X1 2 X2, and X1X2 2 interactions have negative va Hence, the negative values suggest that a maximum TRS yield was observed at the ce point. Moreover, if one of the variables increases while the others decrease above the tre point, the reducing sugar yield amount will not be affected. The negative values served in the regression coefficient for the square terms (X1 2 , X2 2 , X3 2 ) indicate tha optimum response at the centre level decreases the TRS yield amount. The occurren this observation could be due to the increase or decrease of parameter levels from centre points.

Effect of Variables on TRS Yield
The study on the correlations between response and independent variables was ualised using interactive graphs and 3D response plots. The graphs were plotted base the relationships between two independent variables, while another parameter rema constant under optimum conditions. Figures 2 and 3 illustrate the interactions of enzyme loading, hydrolysis time, temperature on the TRS yield. In Figure 2C, the interaction of the two parameters, hyd ysis time and temperature, at the centre point of 20.0 µg/g enzyme loading was depic On the contrary, Figure 2A,B does not exhibit any interaction of the enzyme loading the temperature and time of 48 h and 50 °C at the centre point, respectively. Meanw Figure 3A reveals a higher TRS yield amount based on the higher enzyme activity (u 22.5 µg/g), the temperature at around 45 °C, and the hydrolysis time of 60 h. In Figur and Figure 3C, they indicate that the prolonged hydrolysis time of 72 h had less o effect on the TRS yield amount in the enzymatic hydrolysis process. The effect on the yield amount was noticed to be smaller when the temperature was increased to 40 °C the enzyme loading was lower (<15.0 µg/g). From the interactive level of the ANOVA (Table 2), the interactions of X1X2 and X2X3 were significant. Hence, the results imply the interactions between (i) enzyme loading and temperature and (ii) temperature time may have a substantial effect on the amount of TRS yield. For the interaction of X it was observed to be insignificant. Nonetheless, this factor cannot be fully elimin from the model hierarchy.

Effect of Variables on TRS Yield
The study on the correlations between response and independent variables was visualised using interactive graphs and 3D response plots. The graphs were plotted based on the relationships between two independent variables, while another parameter remained constant under optimum conditions. Figures 2 and 3 illustrate the interactions of enzyme loading, hydrolysis time, and temperature on the TRS yield. In Figure 2C, the interaction of the two parameters, hydrolysis time and temperature, at the centre point of 20.0 µg/g enzyme loading was depicted. On the contrary, Figure 2A,B does not exhibit any interaction of the enzyme loading with the temperature and time of 48 h and 50 • C at the centre point, respectively. Meanwhile, Figure 3A reveals a higher TRS yield amount based on the higher enzyme activity (up to 22.5 µg/g), the temperature at around 45 • C, and the hydrolysis time of 60 h. In Figures 3B and 3C, they indicate that the prolonged hydrolysis time of 72 h had less of an effect on the TRS yield amount in the enzymatic hydrolysis process. The effect on the TRS yield amount was noticed to be smaller when the temperature was increased to 40 • C and the enzyme loading was lower (<15.0 µg/g). From the interactive level of the ANOVA test (Table 2), the interactions of X 1 X 2 and X 2 X 3 were significant. Hence, the results imply that the interactions between (i) enzyme loading and temperature and (ii) temperature and time may have a substantial effect on the amount of TRS yield. For the interaction of X 1 X 3 , it was observed to be insignificant. Nonetheless, this factor cannot be fully eliminated from the model hierarchy.
In Figure 3A, increased enzyme loading was observed to greatly improve cellulose conversion at the optimum temperature (>45 • C). On the other hand, Figure 3C illustrates that the cellulose conversion ratio improves with the extended hydrolysis time at the optimum temperature of more than 45 • C. Additionally, the figure reveals that the TRS yield was considerably affected by the interaction between temperature and time. The influence factors on the TRS yield were proven through the interaction between these two variables in Figure 2C and the p-value (0.0009) in the ANOVA test. In Figure 3A, increased enzyme loading was observed to greatly improve cellulose conversion at the optimum temperature (>45 °C). On the other hand, Figure 3C illustrates that the cellulose conversion ratio improves with the extended hydrolysis time at the optimum temperature of more than 45 °C. Additionally, the figure reveals that the TRS yield was considerably affected by the interaction between temperature and time. The influence factors on the TRS yield were proven through the interaction between these two variables in Figure 2C and the p-value (0.0009) in the ANOVA test.

Validation of Optimisation Conditions
The expected hydrolysis conditions for the highest TRS yield amount (105.57 mg/mL) were measured, in which enzyme loading: X1 = 22.9 µg/g, temperature: X2 = 45 °C, and time: X3 = 69.06 h. Thus, the models employed in this study established the ideal operating

Validation of Optimisation Conditions
The expected hydrolysis conditions for the highest TRS yield amount (105.57 mg/mL) were measured, in which enzyme loading: X 1 = 22.9 µg/g, temperature: X 2 = 45 • C, and time: X 3 = 69.06 h. Thus, the models employed in this study established the ideal operating conditions for the enzymatic hydrolysis process. Resultantly, the model produces a high output of TRS yield. Nonetheless, an additional experimental procedure was performed to verify the expected TRS yield amount at the optimum hydrolysis conditions. The actual experimental value obtained for TRS was 108.74 mg/mL. The result was observed to be near the estimated value (105.57 mg/mL) and has a predicted error of ±2.91 mg/mL. Therefore, the validation data carried out in this study was significant. Subsequently, the validation data was further used in the kinetic analysis of enzymatic hydrolysis.

Kinetic Study of the Enzymatic Hydrolysis
The experimental findings in Figure 4 reveal the influence of the NaOH pre-treatment process and enzyme concentration on TRS yield during the enzymatic hydrolysis process. The effect of the NaOH pre-treatment process on the substrate surface was highlighted for all enzyme concentrations in all-time curves using the Chrastil model. It was demonstrated that the TRS yield increases with time when the enzyme concentration increases from 0.5 to 22.9 µg/g. Additionally, the hydrolysis processes of untreated and pre-treated PALF by NaOH was compared. At the same cellulase concentration (22.9 µg/g), the pre-treated PALF observed a 93% higher TRS yield amount than the untreated PALF. successfully achieved, which was indicated by the non-linear regression graph. Moreover, the reaction order characteristics were measured using the parameter n. The changes in the structural diffusion resistance coefficient were observed for n at various enzyme concentrations. Hence, if n tends to approach unity (0.9 to 1.0), then the diffusion resistance becomes negligible, and the reaction is a first-order reaction. Nonetheless, when n is small (0.5 to 0.6), the system was constrained by diffusion resistance. Additionally, the increasing value of rate constant k was observed from 1.30 × 10 −3 mg/mL.h (untreated PALF) to 5.77 × 10 −3 mg/mL.h (pre-treated with NaOH).   The Chrastil model used in Equation (4) was implemented to obtain the product formation results. The diffusion properties of the enzyme system were determined by analysing the curves and the parameters n (structural diffusion resistance constant) and k (rate constant). The results in Table 4 demonstrate excellent agreement of the experimental values for each enzyme concentration. In Figure 4, the experimental result of R 2 > 0.90 was successfully achieved, which was indicated by the non-linear regression graph. Moreover, the reaction order characteristics were measured using the parameter n. The changes in the structural diffusion resistance coefficient were observed for n at various enzyme concentrations. Hence, if n tends to approach unity (0.9 to 1.0), then the diffusion resistance becomes negligible, and the reaction is a first-order reaction. Nonetheless, when n is small (0.5 to 0.6), the system was constrained by diffusion resistance. Additionally, the increasing value of rate constant k was observed from 1.30 × 10 −3 mg/mL.h (untreated PALF) to 5.77 × 10 −3 mg/mL.h (pre-treated with NaOH).  Table 4 investigates the diffusion-limited system by summarising the kinetic variables, which were determined by the regression analysis in Equation (4). Figure 5 illustrates the experimental results of the initial hydrolysis rate and enzyme concentration. The kinetic parameters of the modified MM model were obtained using Equation (3), which gave a maximum rate of V emax = 93.45 mg/mL and a half-saturation constant of K e = 191.05 mg/mL. According to the obtained kinetic parameters, the correlations between the initial enzyme concentration and the hydrolysis rate can be constructed as continuous graphs. Additionally, the modified MM model enables a reasonable PALF kinetic hydrolysis behaviour estimation even when the diffusion and structural features were neglected. Moreover, using the modified MM model, the kinetic analysis of enzyme concentration and TRS production depicts a good correlation (R 2 = 0.99). Nevertheless, the sole use of the modified MM model was unable to determine whether the enzymatic system is limited by the diffusion process. enzyme concentration and the hydrolysis rate can be constructed as continuous graphs. Additionally, the modified MM model enables a reasonable PALF kinetic hydrolysis behaviour estimation even when the diffusion and structural features were neglected. Moreover, using the modified MM model, the kinetic analysis of enzyme concentration and TRS production depicts a good correlation (R 2 = 0.99). Nevertheless, the sole use of the modified MM model was unable to determine whether the enzymatic system is limited by the diffusion process.

Discussion
Generally, the lignocellulosic biomass structure, which includes PALF, is primarily constituted of hemicellulose, cellulose, and lignin. In the production of reducing sugar, the disruption of the lignocellulosic cell wall was a crucial factor in optimising the yield amount of reducing sugar from cellulose. Moreover, the enzymatic approach applied after the pre-treatment process for lignin removal can be developed to become a potential green

Discussion
Generally, the lignocellulosic biomass structure, which includes PALF, is primarily constituted of hemicellulose, cellulose, and lignin. In the production of reducing sugar, the disruption of the lignocellulosic cell wall was a crucial factor in optimising the yield amount of reducing sugar from cellulose. Moreover, the enzymatic approach applied after the pre-treatment process for lignin removal can be developed to become a potential green technology. Hence, the enzymatic approach offers the elimination of excessively strong acids or alkalis during the hydrolysis process. Furthermore, this study was conducted to improve the enzymatic hydrolysis of the PALF cellulose. From the ANOVA analysis, the most significant parameter observed in TRS yield amount was hydrolysis time (X 3 ), followed by temperature (X 2 ), and enzyme loading (X 1 ). Nonetheless, when compared to the other parameters, the hydrolysis time of the enzymatic hydrolysis process can have a greater effect on the TRS yield amount. Similar results were reported on these parameters, which have a substantial effect on the reducing sugar yield of potato peel residues and Ayous cellulose, respectively [4,24]. Nevertheless, the enzymatic hydrolysis process of alkaline pre-treated corn stover powder reported by Zhang and Wu was also affected by biomass loading [25]. The study demonstrated a significant (p < 0.05) result for TRS yield, while temperature revealed an insignificant (p > 0.05) effect on TRS yield. Therefore, it was concluded that the results obtained may vary depending on the type of lignocellulosic biomass used and the choice of pre-treatment process.
Based on the observation of the optimisation analysis study, the combined effect of a 72-h reaction time and 20.0 µg/g enzyme loading at 50 • C resulted in the highest TRS yield of 105.06 mg/mL from PALFs. After validation of the optimisation process, the TRS yield depicts a tremendous increase of up to 84% from the initial TRS value. The increment of the TRS yield amount reveals that the optimisation process using BBD was successful in improving the enzymatic hydrolysis process of PALF. Furthermore, the dominant effect of hydrolysis time observed in this study was probably due to the longer time required for cellulase to reach the cellulose molecules within the pre-treated PALF structure. The dominant effect of hydrolysis time was assumed when the mass transfer hindrance was high. Hence, the TRS yield amount observed in this study increased when the hydrolysis time was increased from 24 h to 72 h. Moreover, Liu et al. reported that a specific duration was required since enzymatic hydrolysis is a heterogeneous reaction [26]. Additionally, the enzymatic hydrolysis process usually involves cellulase adsorption, the formation of substrate-cellulase complexes, and saccharide dissolution.
For heterogeneous hydrolysis reactions, the Chrastil model was implemented to determine whether the enzyme system was limited by the diffusion process. Carillo et al. described that when the n value changes, alterations occur in the substrate-enzyme system structure and diffusion system [17]. In this study, the pre-treated PALFs displayed a higher n value (0.41) when compared to the untreated PALFs (0.35). The higher n value indicates that the pre-treated PALFs have less resistance to the diffusion process. Therefore, it can be concluded that the NaOH pre-treatment process enhances the movement of molecules in pre-treated PALF. The enhancement movement of molecules was due to the swelling effect and hydrolysis of hemicellulose and lignin. Furthermore, the increase in enzyme concentration produces a decrease in the n value from 0.79 to 0.31. The decrease in the n value suggests that the diffusion resistance increases, which was possibly caused by molecular steric factors induced by enzyme molecules. Lee and Fan interpreted this phenomenon by assuming that the initially adsorbed enzymes produce a single layer on the surface of the fibre [27]. Additionally, the excessive enzyme molecules were absorbed, which built an additional layer. Likewise, the adsorbed enzymes on the first layer could directly contribute to the hydrolysis process by restricting enzyme diffusion through the structure of the substrate. Hence, the adsorbed enzymes on the first layer increased the diffusion resistance.
The reduction in k values at higher enzyme concentrations was due to the critical role of hydrolysis inhibition mechanisms in product formation [28]. The reduced accessibility towards cellulose was due to the increase in enzyme concentration, which increased limiting diffusion, thus the catalytic activity of the enzyme decreased. In this work, the pre-treatment process of PALF using NaOH enhances the catalytic activity of the hydrolysis process. The improved catalytic activity results in the formation of a more accessible substrate for the hydrolysis process. Additionally, the enzymes were made to be more favourable to the substrate. Nevertheless, the k value obtained in this study decreased (4.98 × 10 −4 mg/mL.h). The decrease in k value was caused by the increase in enzyme concentration due to the limiting substrate saturation, specifically when 30.0 µg/g enzyme was used.

Conclusions
The enzymatic hydrolysis conditions of the pre-treated PALF were optimised using BBD to maximise the TRS yield amount. The optimisation process was carried out by manipulating the enzyme loading, temperature, and hydrolysis time. Hence, the highest TRS yield amount achieved was 108.74 mg/mL through control of the enzyme loading, temperature, and hydrolysis time. Additionally, the optimisation of the enzymatic hydrolysis conditions greatly improved the TRS yield from its initial value before optimisation. Furthermore, the kinetic evaluation was investigated using modified MM and Chrastil models for a limited diffusion system of enzymatic hydrolysis. The kinetic results using the modified MM and Chrastil models demonstrated good agreement with the experimental data. In the heterogeneous enzymatic hydrolysis of PALF, the interaction of cellulase and substrate was found to be highly dependent on the substrate's accessibility to cellulase molecules. Therefore, the optimisation and kinetics data obtained in this study provide essential information for the yield amount prediction in the bioconversion of PALF to sugar. Since PALF is a potential source for bioproduct formation, the obtained PALF hydrolysate could be used as a promising carbon source for probiotic growth and biofuel production.