Next Article in Journal
Wikstroemiaganpi Extract Improved Atopic Dermatitis-Like Skin Lesions via Suppression of Interleukin-4 in 2,4-Dinitrochlorobenzene-Induced SKH-1 Hairless Mice
Next Article in Special Issue
Antiangiogenic and Antioxidant In Vitro Properties of Hydroethanolic Extract from açaí (Euterpe oleracea) Dietary Powder Supplement
Previous Article in Journal
A Review on Cs-Based Pb-Free Double Halide Perovskites: From Theoretical and Experimental Studies to Doping and Applications
Previous Article in Special Issue
Combined Antimicrobial Effect of Bio-Waste Olive Leaf Extract and Remote Cold Atmospheric Plasma Effluent
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Statistical Optimization of Flavonoid and Antioxidant Recovery from Macerated Chinese and Malaysian Lotus Root (Nelumbo nucifera) Using Response Surface Methodology

1
Main Campus, Tunku Abdul Rahman University College, Jalan Genting Klang, Kuala Lumpur 53300, Malaysia
2
Bioprocess Division, School of Industrial Technology, Universiti Sains Malaysia, Gelugor 11800, Malaysia
3
Analytical Biochemistry Research Center (ABrC), sains@usm Campus, Universiti Sains Malaysia, University Innovation Incubator Building, Lebuh Bukit Jambul, Bayan Lepas 11900, Malaysia
*
Authors to whom correspondence should be addressed.
Submission received: 15 March 2021 / Revised: 25 March 2021 / Accepted: 26 March 2021 / Published: 1 April 2021

Abstract

:
In this study, the combination of parameters required for optimal extraction of anti-oxidative components from the Chinese lotus (CLR) and Malaysian lotus (MLR) roots were carefully investigated. Box–Behnken design was employed to optimize the pH (X1: 2–3), extraction time (X2: 0.5–1.5 h) and solvent-to-sample ratio (X3: 20–40 mL/g) to obtain a high flavonoid yield with high % DPPHsc free radical scavenging and Ferric-reducing power assay (FRAP). The analysis of variance clearly showed the significant contribution of quadratic model for all responses. The optimal conditions for both Chinese lotus (CLR) and Malaysian lotus (MLR) roots were obtained as: CLR: X1 = 2.5; X2 = 0.5 h; X3 = 40 mL/g; MLR: X1 = 2.4; X2 = 0.5 h; X3 = 40 mL/g. These optimum conditions gave (a) Total flavonoid content (TFC) of 0.599 mg PCE/g sample and 0.549 mg PCE/g sample, respectively; (b) % DPPHsc of 48.36% and 29.11%, respectively; (c) FRAP value of 2.07 mM FeSO4 and 1.89 mM FeSO4, respectively. A close agreement between predicted and experimental values was found. The result obtained succinctly revealed that the Chinese lotus exhibited higher antioxidant and total flavonoid content when compared with the Malaysia lotus root at optimum extraction condition.

1. Introduction

Oxidative diseases such as cancer are the leading cause of deaths, and this accounted for 7.6 million global deaths in the year 2008 alone, as reported by Cancer Research UK [1]. The future projection had forecasted 13.1 million deaths in the year 2030 [1]. According to the study conducted by Sharma et al. [2], low antioxidant status and increased levels of oxidative stress biomarkers were observed in cancer patients. This indicated that the supplementation of natural antioxidants plays an important role in altering the efficacy of cancer chemotherapy. However, synthetic antioxidants that are widely used in the food and pharmaceutical industry were suspected to be carcinogenic [3,4]. Therefore, special attention has been shifted from synthetic antioxidants to the search for natural antioxidants which is advantageous and does not carry any negative effects on human health [5]. Flavonoid is a large class of secondary metabolites produced by higher plants that are challenged with environment stresses. Besides the effects as antioxidants and free radical scavengers, flavonoids also exert other biological effects such as antibacterial, antidiabetic, antiviral and anti-inflammatory actions [6,7].
Lotus plant is widely distributed throughout East Asia, China, Australia and India, not only for holding religious significance but also well known for its nutritional and medicinal value [8]. Various studies were conducted on the isolation of natural flavonoids from lotus (Nelumbo nucifera) parts including seeds, leaves, flowers and rhizomes [9,10]. All the parts of the lotus are reported to possess a significant amount of flavonoids and exhibit convincing antioxidant activities [11]. To date, there are limited studies on the determination of effective conditions in the extraction of flavonoids from lotus root. It is important to know that the lotus species grown in Malaysia was never explored. The preliminary study showed that these flavonoids were highly water-soluble and hardly extracted by organic solvents such as ethanol, methanol and acetone. The main objectives of this study were to compare the flavonoid contents and activities from both Malaysian lotus roots (MLR) and the Chinese lotus roots (CLR) and to access the effect of pH, extraction temperature, extraction time, and buffer-to-sample ratio in the extraction of flavonoids. The extraction conditions were then optimized using Response Surface Methodology (RSM) in order to determine the best combination of parameters for optimal flavonoid recovery.

2. Results and Discussion

2.1. Investigating the Effects of the Parameters on Antioxidant Activities

In this study, the effects of four parameters on the antioxidant capacity of the extracts were tested using a single-factor experiment. Parameters effects of variables such as the pH, extraction time, solvent-to-sample ratio and extraction temperature on the total flavonoid content (TFC), DPPH scavenging activity and ferric-reducing antioxidant power (FRAP) were further succinctly investigated as presented in the preceding sub-sections.

2.2. Effect of pH

The extraction process was conducted by varying the pH condition within the range of 2 and 6. The result obtained showed an increasing trend in the pH values which invariably triggered a corresponding decrease in the TFC yield for both Malaysia lotus (MLR) and China lotus (CLR) roots extracts, as presented in Figure 1. As the pH rises to a threshold of 2, a higher total flavonoid content was recorded for Malaysia lotus (MLR) and China lotus (CLR) roots extracts with values 0.699 mg PCE/g sample and 0.533 mg PCE/g sample, respectively. The results obtained indicated that the acidic condition is more effective in the total recovery of flavonoid content from MLR and CLR. The synthesized flavonoid from the plant matrix is transported to the vacuoles and cell walls as reported by Agati et al. [12]. However, the acidic environment could result in the disruption of the cellulosic cell walls by breaking the chemical bonds between the flavonoid and cell walls [13]. This subsequently led to an eventual exudation of the flavonoid content and their solubilization in the extraction solvent.
Figure 2a,b shows the effect of pH variation on the in vitro DPPH radical scavenging activity and ferric-reducing antioxidant activities of the plant extracts. The result obtained from the China lotus root extract clearly depicted an inverse proportionality. The %DPPHsc and FRAP value significantly dropped from 57.13 to 40.29%DPPHsc and 4.745 to 3.086 mM FeSO4, respectively. Meanwhile, for MLR extract, a similar pattern was observed (61.86% to 38.50%DPPHsc and 4.745 to 3.086 mM FeSO4, respectively). These results suggest that flavonoid extracted from the lotus roots is responsible for the free radical scavenging and ferric-reducing activity. Therefore, the high concentration of flavonoid extracted from the lotus roots is largely responsible for the strong antioxidant activities [14].

2.3. Effect of Extraction Time

In this study, extraction time was varied from 0.25 to 5 h, at constant pH, solvent-to-sample ratio and temperature. Extract from China lotus root gave the highest TFC at 1 h of extraction time (0.185 mg PCE/g sample), while the lowest TFC was obtained at 5 h of extraction time (0.124 mg PCE/g sample). On the other hand, extract from MLR gave the highest TFC at 0.5 h of extraction (0.237 mg PCE/g sample) while the lowest TFC was obtained at 5 h of extraction (0.114 mg PCE/g sample) (Figure 3).
The decrement suggested that the extension of extraction time could cause the released flavonoids to be exposed to the environment at high temperatures for a longer period, thus increasing their chance of decomposition. As discussed previously, antioxidant activities like DPPH radical scavenging and ferric-reducing activity are greatly influenced by the concentration of flavonoid extracted. Therefore, the trends of the two antioxidant activities conducted (Figure 4a,b) were observed to be similar to the trend obtained from TFC studies (Figure 3). From the result, extraction time within 1.5 h was chosen for optimization study. A similar result was obtained by Nepote et al. [15], who found that 10 min was the optimum time to extract phenolic compounds from peanut skin.

2.4. Effect of Solvent-to-Sample Ratio

A range of 10–50 mL/g solvent-to-sample ratio was studied in the extraction of phytochemicals from both Malaysia lotus (MLR) and China lotus (CLR) roots. In general, the increase in the solvent-to-sample ratio encouraged the increase in TFC and antioxidant activities in both Malaysia lotus (MLR) and China lotus (CLR) roots extracts. For China lotus root extract, the highest TFC was obtained with the solvent-to-sample ratio of 40 mL/g (0.393 mg PCE/g sample), followed by 50 mL/g (0.359 mg PCE/g sample) and 30 mL/g (0.207 mg PCE/g sample) (Figure 5).
On the other hand, the highest TFC content obtained from MLR extract was at a solvent-to-sample ratio of 50 mL/g (0.797 mg PCE/g sample), then 40 mL/g (0.699 mg PCE/g sample) and 30 mL/g (0.504 mg PCE/g sample). Moreover, the antioxidant activity of Malaysia lotus (MLR) and China lotus (CLR) roots increase as the solvent-to-sample increases (Figure 6). The result indicated that the increase in solvent volume is directly proportional to the antioxidant activities [4]. The result is consistent with the mass transfer principle, which describes that the concentration gradient between the sample and the solvent is the driving force during mass transfer [16]. Therefore, the higher buffer-to-sample ratio used will create a higher driving force for the extraction process. Thus, a larger volume of solvents will be able to extract more flavonoid compounds [17].

2.5. Effect of Temperature

TFC contents extracted from both Malaysia lotus (MLR) and China lotus (CLR) roots with their DPPH and FRAP assays are presented in Figure 7a–c. There were no significant differences in the range of extraction temperature at 25–60 °C. Thus, the effect of temperature was not studied in optimization studies and was set at room temperature as a constant in the further experiment.

2.6. Comparison of CLR and MLR Extracts

Malaysia lotus (MLR) and China lotus (CLR) roots were compared in this research as these two cultivars are the most common lotus root found in local markets in Malaysia. When comparing the physical appearance of raw lotus root, CLR showed a pinkish-beige and larger (in terms of diameter) root, whereas MLR is relatively slimmer and whitish-beige. However, when the root was cooked, the CLR turned into a purplish-brown colour, whereas MLR remained as beige colour and slightly translucent. The textures of both lotus roots were also different. Therefore, it was believed that these two cultivars of lotus root are different in compositions. According to the results, CLR extract shows higher TFC (0.699 mg PCE/g sample) compared to MLR extract (0.533 mg PCE/g sample) at pH 2. CLR extract also exhibited stronger free radical scavenging activity (21.96 %DPPHsc) compared to MLR extract (13.53 %DPPHsc) at low solvent-to-sample ratio (10 mL/g). At long extraction time (5 h), CLR extract also shows a more stable and higher reducing power (2.18 mM FeSO4) than MLR extract (1.73 mM FeSO4). These results could be due to the climate and botanical aspects of the plantation of different cultivars lotus root. According to the research conducted by Jang et al. [9], the growth of lotus in monsoon climate is slower than the growth in summer and thus affects the synthesis of flavonoids in the plant cell. Also, progressive changes in traits like the production of phytochemicals in the root could happen during seasonal variation especially from the dormant to the active phase. Many other factors could also influence the secondary metabolite profile of the root such as geographic location, age of the plant, season, associated microflora, nutritional status, and environmental stress [17]. Ultimately, the result obtained clearly revealed CLR as a better source of water-soluble flavonoids when compared to MLR. However, the optimization study was subsequently conducted to obtain the highest flavonoid content.

2.7. Optimization Study

The effects of extraction factors, i.e., pH (X1), time (X2), and solvent-to-sample ratio (X3), were included in the optimization study. These parameters (X1: pH 2–3; X2: 0.5–1.5 h; X3: 20–40 mL/g) were selected during the preliminary studies (Section 3.1). The three responses of interest were TFC, %DPPHsc, and FRAP. The results of 17 runs using BBD are shown in Table 1a,b and this includes the design of the experiment, observed responses and predicted values. A close agreement between experimental and predicted values was subsequently attained. From the results, TFC ranged from 0 to 0.68 mg PCE/g CLR extract and 0.1–0.58 mg PCE/g MLR extract were reported. For CLR extract, the highest TFC value (0.68 mg PCE/g sample) was achieved under condition of X1 = pH 2.5, X2 = 1.5 h and X3 = 40 mL/g. Meanwhile, for MLR extract, the highest TFC value was of 0.58 mg PCE/g sample which was obtained under the condition of X1 = pH 2.5, X2 = 0.5 or 1.5 h and X3 = 40 mL/g. As for the antioxidant activities, the ranges of %DPPHsc of CLR and MLR extracts were reported as 7.2–50.5 %DPPHsc and 12.88–30.22 %DPPHsc, respectively. Both highest %DPPHsc (50.5 and 30.22 %DPPHsc) from CLR and MLR extracts, respectively, were found in conditions of X1 = pH 2.5, X2 = 30 min and X3 = 40 mL/g. At the same condition, CLR extract gave the highest FRAP value (2.20 mM FeSO4), but MLR extract gave the highest FRAP value (1.96 mM FeSO4) at condition X1 = pH 2.5, X2 = 60 min and X3 = 30 mL/g. These conditions varied depending on the response required. The optimum extraction condition was therefore investigated in order to obtain desirable TFC, %DPPHsc, and FRAP value.

2.8. Model Fitting

Table 1a,b presents the results of quadratic models fitted to the data. The results of the analysis of variance (ANOVA) indicated that the contribution of the quadratic model was significant for the responses. Equations (1), (3), and (5) represent the fitted quadratic models for TFC, %DPPHsc, and FRAP of CLR and MLR extracts; on the other hand, Equations (2), (4), and (6) are of MLR extracts, respectively. The significance of each coefficient was determined using the F-test and p-value in Table 2a,b. If the absolute F-value becomes greater and the p-value becomes smaller, the corresponding variables would be more significant [4].
It could be observed that the largest effects on the TFC of CLR extract were the linear term of solvent-to-sample ratio ( x 3 ) and quadratic term of pH ( x 1 2 ) . The results shown in Table 2a suggested that the change of pH and solvent-to-sample had significant effects (p < 0.05) on the TFC of CLR extract. Time, on the other hand, did not have any significant (p > 0.05) contribution in TFC of CLR extract. The coefficient of determination (r2) of the predicted models in this response of CLR was 0.9765, indicating a good correlation of the data to the model. However, the p-value for lack of fit was 0.0264 with a lack of good fit for the mathematical models in Equation (1)
T F C =   0.59 + 0.034 x 1 + 0.046 x 2 + 0.091 x 3 0.043 x 1 2 0.055 x 2 2 + 0.019 x 3 2 + 8.96 × 10 3 x 1 x 2 + 0.057 x 1 x 3 + 2.29 × 10 3 x 2 x 3
However, the significant effects on the TFC of MLR extract have linear term for pH ( x 1 ) , linear term for solvent-to-sample ratio ( x 3 ) , quadratic term for pH ( x 1 2 ) , quadratic term for solvent-to-sample ratio ( x 3 2 ) , interaction term for pH and time ( x 1 x 2 ) and interaction term for time and solvent-to-sample ratio ( x 2 x 3 ) (Table 2b). The results showed that all three factors (pH, time and solvent-to-sample ratio) contributed significantly (p < 0.05) in TFC of MLR extract. The coefficient of determination (r2) of the predicted models in the response of MLR was 0.9873. Although the p-value for lack of fit (0.0003) shows not a good fit for mathematic models in Equation (2), verification of an optimized extraction conditions will be carried out.
T F C =   0.49 0.059 x 1 0.018 x 2 + 0.11 x 3 0.21 x 1 2 5.17 × 10 3 x 2 2   0.036 x 3 2 0.036 x 1 x 2 0.024 x 1 x 3 + 0.03 x 2 x 3
For CLR extract, it could be seen that the quadratic term of pH ( x 1 2 ) showed the largest effect on %DPPHsc, followed by the linear term of solvent-to-sample ratio ( x 3 ) and the interaction term of pH and solvent-to-sample ratio ( x 1 x 3 ) . The coefficient of determination (r2) of the predicted models in the response was 0.9937. However, the p-value for the lack of fit was 0.0115, which suggested not a good fit for the mathematical model Equation (3).
% D P P H s c =   45.46 1.00 x 1 0.51 x 2 + 3.26 x 3 31.34 x 1 2 1.53 x 2 2   + 0.48 x 3 2 2.18 x 1 x 2 + 3.73 x 1 x 3 0.016 x 2 x 3
The linear term of pH ( x 1 ) showed the most significant (p < 0.05) effect on %DPPHsc of MLR extract, followed by a linear term of solvent-to-sample ratio ( x 3 ) , a quadratic term of pH ( x 1 2 ) and linear term of time ( x 2 ) . The other quadratic terms ( x 2 2 , x 3 2 ) and all the interaction terms x 1 x 2 ,   x 1 x 3 and x 2 x 3 did not have significant (p > 0.05) effects. The coefficient of determination (r2) of the predicted models in this response was 0.9474, whereas the p-value for lack of fit was 0.1988, which suggests a relatively good fit to the mathematic models in Equation (4).
% D P P H s c =   25.34 4.84 x 1 1.70 x 2 + 4.08 x 3 4.27 x 1 2 + 0.37 x 2 2 1.14 x 3 2 1.40 x 1 x 2 0.076 x 1 x 3   1.76 x 2 x 3
In terms of FRAP assay, the only variable with significant (p < 0.05) effect on FRAP value of CLR extract was the quadratic term of pH ( x 1 2 ) , whereas other terms were not significant (p > 0.05). The coefficient of determination (r2) of the predicted model was 0.9212 and the p-value for lack of fit was 0.4533, which was a good fit for the mathematical model in Equation (5).
F R A P =   1.87 0.011 x 1 0.031 x 2 + 0.11 x 3 0.85 x 1 2 0.034 x 2 2 + 0.053 x 3 2 0.073 x 1 x 2 + 0.14 x 1 x 3 0.031 x 2 x 3
Meanwhile, for MLR extract, the variables that gave the largest effect on the FRAP value was the quadratic term of pH ( x 1 2 ) , followed by a linear term of solvent-to-sample ratio ( x 3 ) and linear term of pH ( x 1 ) . Time, on the other hand, did not contribute any significant (p > 0.05) effect to the response. The coefficient of determination (r2) was 0.8865 and the p-value for lack of fit was 0.9720, which suggested that it was an excellent fit for the mathematic model in Equation (6).
F R A P =   1.71 0.13 x 1 0.022 x 2 + 0.18 x 3 0.34 x 1 2 3.96 × 10 3 x 2 2 0.055 x 3 2 0.027 x 1 x 2 0.015 x 1 x 3 0.012 x 2 x 3

2.9. Response Surface Model and Contour Plot

3D response surfaces shown in Figure 8, Figure 9 and Figure 10 demonstrated the changes in TFC, %DPPHsc and FRAP values as a function of pH and solvent-to-sample ratio, while extraction time was set at 3 different levels (i.e., 0.5, 1 and 1.5 h). The 3D response surfaces for TFC of CLR extract as a function of pH and solvent-to-sample ratio were given in Figure 8a(i–iii). The overall results showed that the extraction time had no significant (p > 0.05) effect on TFC value. Results also showed that a higher solvent-to-sample ratio would give higher TFC yield, where at pH 2.5 and the buffer-to-sample ratio of 40:1, TFC reported the highest for the 0.679 mg PCE/g sample. As mentioned previously, TFC yield increased as pH increased from pH 2 to pH 2.5, ranging from 0 to 0.679 mg PCE/g sample and then decreased to 0.098 mg PCE/g sample as pH increased from pH 2.5 to pH 3. This suggests that pH 2.5 (as described previously) is the optimum pH for the extraction of flavonoids.
For the TFC of MLR extract, the 3D surfaces as a function of pH and extraction time are given in Figure 8b(i–iii). Results showed that at low solvent-to-sample ratio (20:1, v/w), an increase of TFC could be observed by increasing pH up to pH 2.5. Further increment in pH to pH 3.0 led to a decrease in TFC. A similar trend could be observed at a solvent-to-sample ratio at 30:1 and 40:1 (v/w). Extension in the extraction time showed no significance (p > 0.05) at all three levels of the solvent-to-sample ratio. However, the interaction between pH and extraction time showed a significant (p < 0.05) effect in this experiment. At low solvent-to-sample ratio (20:1, v/w), TFC increased with increasing pH up to 2.25 and extraction time (~1 h) and then decreased at higher pH (~pH 2.75) with further extended extraction period. A similar trend was found at a solvent-to-sample ratio of 30:1 (v/w). On the other hand, at a high solvent-to-sample ratio (40:1, v/w), increasing pH and extraction time led to higher TFC (~0.589 PCE/g sample). Longer extraction time showed higher efficiency in the extraction of flavonoids as it provides sufficient time for the solutes to expose to the extraction solvents. Similar results were obtained by Xu et al. [18]. These results suggested a positive effect of increasing extraction time (less than 2 h) on the extraction of flavonoid from lotus roots.
Furthermore, the 3D surfaces for %DPPHsc of CLR extract displayed in Figure 9a(i–iii) showed similar results obtained in TFC. Overall, the time had no significant (p > 0.05) effects, while pH exhibited a major effect on DPPH radical scavenging activity. At fixed extraction time and solvent-to-sample ratio, a maximum %DPPHsc could be observed at region pH 2.5. Lower %DPPHsc were obtained at both the lower and higher pH regions. Interaction between pH and solvent-to-sample ratio contributed a significant (p < 0.05) effect on %DPPHsc. At 0.5 h extraction time, given the increment in solvent-to-sample ratio (from 20 mL/g to 40 mL/g) at low pH (pH 2), only a slight increment in %DPPHsc (~7.74 %DPPHsc) was observed. However, when the pH increased up to 2.5 and above, an apparent increment in %DPPHsc could be observed with the increment of the solvent-to-sample ratio.
For MLR, 3D surfaces of %DPPHsc showed that there was no significant (p > 0.05) effect in the interaction between factors. Significant (p < 0.05) effect could only be observed in pH. At all three levels of extraction time and solvent-to-sample ratio, increment in pH from 2 to 2.5 did not show any significant (p > 0.05) difference in %DPPHsc (~18.81–22.81%), however, further increment in pH up to pH 3, %DPPHsc were significantly (p < 0.05) decreased to ~6.58%. This result was similar to the result obtained from CLR. In general, the results suggested that these flavonoids were favored at ex traction conditions of pH 2.5 which could give a high DPPH radical scavenging activity [19].
Also, the 3D response surfaces for FRAP of CLR extract as a function of pH and solvent-to-sample ratio are given in Figure 10a(i–iii). In general, pH affected the activity significantly (p < 0.05); from all three levels of extraction time, FRAP value reached maximum (~1.63–1.93 mM FeSO4) at pH~2.25, then decreased with further increment of pH. On the other hand, the solvent-to-sample ratio did not increase significantly (p > 0.05) at all three levels of extraction time. The result showed that the interaction between pH and solvent-to-sample ratio did not significantly (p < 0.05) affect the FRAP of CLR. For FRAP of MLR extract, the 3D response surfaces as a function of pH and solvent-to-sample ratio are given in Figure 10b(i–iii). At 0.5 h of extraction time, increment in solvent-to-sample ratio (from 20:1 to 40:1, v/w) significantly (p < 0.05) increased FRAP value from (~1.059–1.681 mM FeSO4). At the same extraction time, the increment of pH up to pH 2.25 significantly (p < 0.05) increased the FRAP value, from ~1.266 to ~1.473 mM FeSO4 at 20 mL/g solvent-to-sample ratio. However, continuous increment from pH ~2.5 led to a significant (p < 0.05) decrement in FRAP (<1.059 mM FeSO4). Similar to the result obtained from CLR, the interaction between pH and solvent-to-sample ratio did not have a significant (p > 0.05) effect on FRAP assay.

2.10. Validation of Predictive Models

Based on the findings, a verification study was performed to evaluate the optimal operating conditions for the extraction with a high yield of antioxidant content (TFC) and high antioxidant activities (%DPPHsc and FRAP) with consideration of efficiency, energy conservation and feasibility of the experiment. An optimal condition to achieve high TFC and high antioxidant activities were determined for both CLR and MLR. Therefore, from the model, the optimal conditions for TFC, %DPPHsc and FRAP of CLR extract were pH 2.5, extraction time of 0.5 h, and the solvent-to-sample ratio of 40 mL/g. On the other hand, the optimal conditions for TFC, %DPPHsc and FRAP of MLR extract were pH 2.4, extraction time of 0.5 h, and the solvent-to-sample ratio of 40 mL/g. These conditions gave TFC, %DPPHsc and FRAP of 0.599 mg PCE/g sample, 48.36 %DPPHsc and 2.07 mM FeSO4 for CLR extract and 0.549 mg PCE/g sample, 29.11 %DPPHsc and 1.89 mM FeSO4 for MLR extract, respectively. Experimental values for TFC, %DPPHsc and FRAP obtained were of the 0.66 mg PCE/g sample, 42.78 %DPPHsc and 1.89 mM FeSO4 for CLR extract and 0.59 mg PCE/g sample, 31.72 %DPPHsc and 1.94 mM FeSO4 for MLR extract, respectively. It could be observed that only small deviations were found between the predicted values and experimental values [20]. Thus, the models could be used to optimize the extraction of flavonoids from both CLR and MLR.

3. Materials and Methods

3.1. Materials

Fresh Chinese and Malaysian lotus root cultivars were purchased from the local market located in Jelutong, Penang, Malaysia. The sample was thoroughly washed with deionized distilled water, subsequently cut into small cubes and lyophilized to remove inherent moisture. The lyophilized samples were then ground into a powder prior to further experiments. All chemicals used in this study were analytical grade and purchased from Sigma-Aldrich (Selangor, Malaysia) and Fluka (Industriestrasse 25, Switzerland).

3.2. Design of Experiments

During the preliminary study, different solvents including methanol, ethanol, acetone, ethyl acetate, hexane, water or enzyme solution (cellulase or hemicellulase, 1% w/v) were used to obtain the most efficient extraction condition. The effect of pH (2–6), temperature (25–60 °C), extraction time (0.25–5 h) and buffer-to-sample ratio (10–50 mL/g) was further investigated using the single-factor experiment. This was achieved by varying one factor and keeping other factors constant using the middle point viz: extraction time (2 h), pH (4), extraction temperature (40 °C), and solvent-to-sample ratio (20 mL/g). The selected factors and levels from the single-factor experiments were subsequently used for the optimization study using the Box–Behnken design.
The parameters of extraction were optimized using response surface methodology (RSM). Box–Behnken design (BBD) was applied while pH (X1: 2–3), extraction time (X2: 0.5–1.5 h) and solvent-to-sample ratio (X3: 20–40 mL/g) were selected as independent variables based on the results of preliminary studies. The experimental design consisted of twelve factorial points and each center point has five replicates (Table 2). Each experiment was carried out in triplicate and the mean values of the observed response. Experiments were conducted randomly to reduce the effects of unexpected variability in the observed responses.
The variables were coded accordingly to the equation:
x   =   X i X o Δ X
where x was the coded value, Xi was the corresponding actual value, Xo was the actual value in the center of the domain, and ΔX was the increment of Xi corresponding to a variation of 1 unit of x.
The mathematical model corresponding to the composite design was:
Y   =   β o   +   i   =   1 3 β i X i + i   =   1 3 β i i X i 2 + i   =   1 2 j   =   1   +   1 3 β i j X i X j + ε
where Y indicated the dependent variables (TFC, %DPPHsc, and FRAP), β o is the model constant, β i , β i i and β i j were the model coefficients which represent linear, quadratic or interaction effects of the variables, whereas ε was the error. This analysis of experimental design data and calculation of predicted responses were carried out using Design-Expert software (version 6.0, Stat-Ease Inc., Minneapolis, MN, USA). In order to verify the validity of the statistical experimental design, additional confirmation experiments were conducted.

3.3. Statistical Analysis

Comparison of means was performed by one-way analysis of variance (ANOVA) followed by Duncan’s test. Statistical analysis (p < 0.05) was performed using Statistical Package for Social Science for Windows, version 20.0 (IBM SPSS, Chicago, IL, USA). The optimal extraction conditions were estimated through regression analysis and three-dimensional response surface plots of the independent variables and each dependent variable.

3.4. Extract Preparation

Briefly, 1 g of the plant sample was added into the extracting solvent at varied solvent-to-sample ratio and pH. The mixture was then homogenized and incubated at different temperatures for selected extraction time and constantly agitated at 250 rpm. After the extraction process had been completed, the mixture was then centrifuged at 4500 rpm for 15 min at 4 °C. The supernatant was then collected and reconstituted to 50 mL of the solvent. The crude extract was stored at −20 °C prior to future analysis. All samples were prepared in triplicate.

3.5. Determination of Total Flavonoid Content (TFC)

Total flavonoid content was determined by colorimetric assay modified from Kim et al. [21]. Distilled water (80 µL) was added to 20 µL of extract. Then, 5% (w/v) of sodium nitrite solution (6 µL) and 10% (w/v) aluminum chloride solution (6 µL) were added to the mixture. Mixtures were incubated at room temperature for 5 min. Then, 40 µL of 1 M NaOH was added to the mixture. The final volume of the reaction mixture was topped up to 200 µL with distilled water and the absorbance measured at 510 nm. All samples were prepared in triplicate. A standard curve was generated with pyrocatechol and results were expressed in mg pyrocatechol equivalent (PCE)/g sample.

3.6. DPPH Radical Scavenging Activity Determination

The scavenging activity on DPPH radical of rambutan peel extracts was determined by the modified methods of Liu et al. [22]. Each extract (5 µL) was added to 150 µL of DPPH solution (0.1 mM). The mixtures were incubated for 30 min in the dark at 30 °C and discolorations were determined at 517 nm. %DPPH scavenged (%DPPHsc) was calculated using Equation (9) below:
D P P H s c = A b s c o n t r o l A b s s a m p l e A b s c o n t r o l   ×   100
where Abscontrol was defined as absorbance of the control at 517 nm and Abssample was defined as the absorbance value of the extracts at 517 nm.

3.7. Ferric-Reducing Power Assay (FRAP)

FRAP assay was conducted by referring to the method of Benzie and Strain [23]. The FRAP reagent was a mixture of 1 mL of 10 mM tripyridyltriazine (TPTZ) solution, 1 mL of 20 mM FeCl3·6H2O and 10 mL of 0.3 M acetate buffer at pH 3.6. FRAP reagent (200 µL) was pre-warmed at 37 °C and then mixed with 2.7 µL of the extract. The reaction mixtures were incubated at 37 °C for 1 h and absorbance at 593 nm was measured. All samples were prepared in triplicate. A standard curve using FeSO4 was obtained.

4. Conclusions

RSM was used to determine the optimum extraction parameters that would give the high yield of TFC and antioxidant activities. According to ANOVA, pH, extraction time and solvent-to-sample ratio showed different effects on the targeted responses depending on the types of lotus root used. Interactions of the extraction parameters were also found. Quadratic models were successfully developed and used in predicting all the responses. The optimal conditions based on the combinations of all responses were obtained (CLR: X1 = 2.5; X2 = 0.5 h; X3 = 40 mL/g; MLR: X1 = 2.4; X2 = 0.5 h; X3 = 40 mL/g). Results showed that the predicted and experimental values were similar. Therefore, it could be concluded that the models obtained could be used to optimize the extraction process of water-soluble flavonoids from both China lotus and Malaysian lotus root.

Author Contributions

S.-J.T.: Conception, writing of the manuscript, C.-K.L.: data analysis and interpretation; C.-Y.G.: supervision and manuscript preparation O.A.O.: manuscript writing and revision. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by an ERGS Grant (203/CAATS/6730103) and Analytical Biochemistry Research Center (ABrC), Universiti Sains Malaysia, Malaysia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Not applicable.

Conflicts of Interest

No potential conflict of interest was reported by the authors.

Sample Availability

Not applicable.

References

  1. Cancer Research UK, 2008. Worldwide Cancer Statistics. Available online: http://www.cancerresearchuk.org/cancer-info/cancerstats/world/ (accessed on 2 March 2021).
  2. Sharma, M.; Rajappa, M.; Kumar, G.; Sharma, A. Oxidant-antioxidant status in Indian patients with carcinoma of posterior one-third of tongue. Cancer Biomark. 2009, 5, 253–260. [Google Scholar] [CrossRef] [PubMed]
  3. Pham, H.N.T.; Van, J.A.S.; Vuong, Q.; Bowyer, M.C.; Scarlett, C.J. Physicochemical Properties, Antioxidant and Cytotoxic Activities of Crude Extracts and Fractions from Phyllanthus amarus. Mol. Biol. Rep. 2019, 46, 3265–3273. [Google Scholar] [CrossRef] [PubMed]
  4. Olalere, O.A.; Gan, C.Y.; Akintomiwa, O.E.; Adeyi, O.; Adeyi, A. Optimization of Microwave-assisted extraction and functional elucidation of bioactive compounds from Cola nitida pod. Phytochem. Analys. 2021, 1–9. [Google Scholar] [CrossRef]
  5. Alara, O.R.; Abdurahman, N.H.; Olalere, O.A. Optimization of microwave-assisted extraction of flavonoids and antioxidants from Vernonia amygdalina leaf using response surface methodology. Food Bioprod. Process. 2018, 107, 36–48. [Google Scholar]
  6. Olalere, O.A.; Gan, C.Y. Microwave-assisted extraction of phenolic compounds from Euphorbia hirta leaf and characterization of its morphology and thermal stability. Sep. Sci. Technol. 2020, 1–13. [Google Scholar] [CrossRef]
  7. Olalere, O.A.; Abdurahman, H.N.; Gan, C.Y. Microwave-enhanced extraction and mass spectrometry fingerprints of polyphenolic constituents in Sesamum indicum leaves. Ind. Crops Prod. 2019, 131, 151–159. [Google Scholar] [CrossRef]
  8. Tseng, H.C.; Tsai, P.M.; Chou, Y.H.; Lee, Y.C.; Lin, H.H.; Chen, H.J. In vitro and in vivo protective effects of flavonoid-enriched lotus seedpod extract on lipopolysaccharide-induced hepatic inflammation. Am. J. Chin. Med. 2019, 47, 153–176. [Google Scholar] [CrossRef] [PubMed]
  9. Lien, D.P.; Lieu, V.T.; Tram, P.B.; Toan, H.T. Influence of Different Drying Techniques on Drying Characteristics and Quality Aspects of Pink Lotus (Nelumbo nucifera) Flowers. J. Adv. Res. Food Sci. Nut. 2020, 3, 16–23. [Google Scholar] [CrossRef]
  10. Chen, G.; Zhu, M.; Guo, M. Research advances in traditional and modern use of Nelumbo nucifera: Phytochemicals, health promoting activities and beyond. Crit. Rev. Food Sci. Nut. 2019, 59, S189–S209. [Google Scholar]
  11. Shen, Y.; Guan, Y.; Song, X.; He, J.; Xie, Z.; Zhang, Y.; Tang, D. Polyphenols extract from lotus seedpod (Nelumbo nucifera Gaertn.): Phenolic compositions, antioxidant, and antiproliferative activities. Food Sci. Nutr. 2019, 7, 3062–3070. [Google Scholar] [PubMed] [Green Version]
  12. Agati, G.; Azzarello, E.; Pollastri, S.; Tattini, M. Flavonoids as antioxidants in plants: Location and functional significance. Plant Sci. 2012, 196, 67–76. [Google Scholar] [CrossRef] [PubMed]
  13. Olalere, O.A.; Abdurahman, N.H.; Bin Mohd Yunus, R.; Alara, O.R.; Gan, C.Y. Synergistic intermittent heating and energy intensification of scale-up parameters in an optimized microwave extraction process. Chem. Eng. Process.-Process Intensif. 2018, 132, 160–168. [Google Scholar]
  14. Akbari, S.; Abdurahman, N.H.; Yunus, R.M.; Alara, O.R.; Olalere, O.A. Extraction, characterization and antioxidant activity of fenugreek (Trigonella-foenum graecum) seed oil. Mater. Sci. Energy Technol. 2019, 2, 349–355. [Google Scholar] [CrossRef]
  15. Nepote, V.; Grosso, N.R.; Guzmán, C.A. Optimization of extraction of phenolic antioxidants from peanut skins. J. Sci. Food Agric. 2005, 85, 33–38. [Google Scholar] [CrossRef]
  16. Al-Farsi, M.A.; Lee, C.Y. Optimization of phenolics and dietary fibre extraction from date seeds. Food Chem. 2008, 108, 977–985. [Google Scholar] [CrossRef] [PubMed]
  17. Sati, P.; Pandey, A.; Rawat, S.; Rani, A. Phytochemicals and antioxidants in leaf extracts of Ginkgo biloba with reference to location, seasonal variation and solvent system. J. Pharmacol. Res. 2013, 7, 804–809. [Google Scholar] [CrossRef]
  18. Xu, Q.; Shen, Y.; Wang, H.; Zhang, N.; Xu, S.; Zhang, L. Application of response surface methodology to optimise extraction of flavonoids from fructus sophorae. Food Chem. 2013, 138, 2122–2129. [Google Scholar] [CrossRef] [PubMed]
  19. Settharaksa, S.; Jongjareonrak, A.; Hmadhlu, P.; Chansuwan, W.; Siripongvutikorn, S. Flavonoid, phenolic contents and antioxidant properties of Thai hot curry paste extract and its ingredients as affected of pH, solvent types and high temperature. Int. Food Res. J. 2012, 19, 1581–1587. [Google Scholar]
  20. Olalere, O.A.; Gan, C.Y. Intensification of microwave energy parameters and main effect analysis of total phenolics recovery from Euphorbia hirta leaf. Food Charact. Meas. 2020, 14, 886–893. [Google Scholar] [CrossRef]
  21. Kim, D.O.; Jeong, S.W.; Lee, C.Y. Antioxidant capacity of phenolic phytochemicals from various cultivars of plums. Food Chem. 2003, 81, 321–326. [Google Scholar]
  22. Liu, L.; Sun, Y.; Laura, T.; Liang, X.; Ye, H.; Zeng, X. Determination of polyphenolic content and antioxidant activity of kudingcha made from Ilex kudingcha CJ Tseng. Food Chem. 2009, 112, 35–41. [Google Scholar] [CrossRef]
  23. Benzie, I.F.; Strain, J.J. Ferric reducing/antioxidant power assay: Direct measure of total antioxidant activity of biological fluids and modified version for simultaneous measurement of total antioxidant power and ascorbic acid concentration. Strain Method Enzymol. 1999, 299, 15–27. [Google Scholar]
Figure 1. Effect of pH on total flavonoid content (TFC).
Figure 1. Effect of pH on total flavonoid content (TFC).
Molecules 26 02014 g001
Figure 2. Effect of pH (a) DPPH scavenging activity and (b) ferric-reducing antioxidant power (FRAP).
Figure 2. Effect of pH (a) DPPH scavenging activity and (b) ferric-reducing antioxidant power (FRAP).
Molecules 26 02014 g002
Figure 3. Effect of extraction time on total flavonoid content (TFC).
Figure 3. Effect of extraction time on total flavonoid content (TFC).
Molecules 26 02014 g003
Figure 4. Effect of extraction time on (a) DPPH scavenging activity and (b) ferric-reducing antioxidant power (FRAP).
Figure 4. Effect of extraction time on (a) DPPH scavenging activity and (b) ferric-reducing antioxidant power (FRAP).
Molecules 26 02014 g004
Figure 5. Effect of solvent-to-sample ratio on total flavonoid content (TFC).
Figure 5. Effect of solvent-to-sample ratio on total flavonoid content (TFC).
Molecules 26 02014 g005
Figure 6. Effect of solvent-to-sample ratio on (a) DPPH scavenging activity and (b) ferric-reducing antioxidant power (FRAP).
Figure 6. Effect of solvent-to-sample ratio on (a) DPPH scavenging activity and (b) ferric-reducing antioxidant power (FRAP).
Molecules 26 02014 g006
Figure 7. Effect of temperature on (a) Total flavonoid content (b) DPPH scavenging activity and (c) ferric-reducing antioxidant power (FRAP).
Figure 7. Effect of temperature on (a) Total flavonoid content (b) DPPH scavenging activity and (c) ferric-reducing antioxidant power (FRAP).
Molecules 26 02014 g007aMolecules 26 02014 g007b
Figure 8. Three-dimensional response surfaces for TFC of (a) China cultivar of lotus root at (i) 0.5, (ii) 1, (iii) 1.5 h extraction time as a function of pH and solvent-to-sample ratio and (b) Malaysia cultivar of lotus root at (i) 20, (ii) 30, and (iii) 40 (mL/g) of solvent-to-sample ratio as a function of pH and extraction time.
Figure 8. Three-dimensional response surfaces for TFC of (a) China cultivar of lotus root at (i) 0.5, (ii) 1, (iii) 1.5 h extraction time as a function of pH and solvent-to-sample ratio and (b) Malaysia cultivar of lotus root at (i) 20, (ii) 30, and (iii) 40 (mL/g) of solvent-to-sample ratio as a function of pH and extraction time.
Molecules 26 02014 g008
Figure 9. Three-dimensional response surfaces for %DPPHsc of (a) China cultivar of lotus root at (i) 0.5, (ii) 1, (iii) 1.5 h extraction time and (b) Malaysia cultivar of lotus root at (i) 0.5, (ii) 1, and (iii) 1.5 h extraction time as a function of pH and solvent-to-sample ratio.
Figure 9. Three-dimensional response surfaces for %DPPHsc of (a) China cultivar of lotus root at (i) 0.5, (ii) 1, (iii) 1.5 h extraction time and (b) Malaysia cultivar of lotus root at (i) 0.5, (ii) 1, and (iii) 1.5 h extraction time as a function of pH and solvent-to-sample ratio.
Molecules 26 02014 g009
Figure 10. Three-dimensional response surfaces for FRAP of (a) China cultivar of lotus root at (i) 0.5, (ii) 1, (iii) 1.5 h extraction time and (b) Malaysia cultivar of lotus root at (i) 0.5, (ii) 1, and (iii) 1.5 h extraction time as a function of pH and solvent-to-sample ratio.
Figure 10. Three-dimensional response surfaces for FRAP of (a) China cultivar of lotus root at (i) 0.5, (ii) 1, (iii) 1.5 h extraction time and (b) Malaysia cultivar of lotus root at (i) 0.5, (ii) 1, and (iii) 1.5 h extraction time as a function of pH and solvent-to-sample ratio.
Molecules 26 02014 g010
Table 1. ANOVA for response surface models: estimated regression model of relationship between response variables (TFC, %DPPHsc and FRAP) and independent variables (X1,X2,X3) from (a) China lotus root and (b) Malaysia lotus root.
Table 1. ANOVA for response surface models: estimated regression model of relationship between response variables (TFC, %DPPHsc and FRAP) and independent variables (X1,X2,X3) from (a) China lotus root and (b) Malaysia lotus root.
(a) China Lotus Root
SourceSum of SquaresDFMean SquareF-Valuep-Value
TFC (mg PCE/g)a
Model0.8990.1032.32<0.0001
Quadratic0.7930.2685.51<0.0001
A0.0110.013.030.1254
B0.0210.025.440.0525
C0.0710.0721.530.0024
A20.7610.76248.23<0.0001
B20.0110.014.140.0814
C20.0010.000.510.5001
AB0.0010.000.100.7558
AC0.0110.014.200.0796
BC0.0010.000.010.9364
Residual0.0270.00
Lack of Fit0.0230.019.670.0264
Total0.9116
%DPPHscb
Model4355.209483.91121.87<0.0001
Quadratic4185.3331395.11351.35<0.0001
A8.0718.072.030.1969
B2.0512.050.520.4957
C85.12185.1221.440.0024
A24134.8614134.861041.35<0.0001
B29.9019.902.490.1583
C20.9610.960.240.6374
AB19.09119.094.810.0644
AC55.53155.5313.990.0073
BC0.0010.000.000.9879
Residual27.7973.97
Lack of Fit25.5938.5315.500.0115
Total4383.0016
FRAP (mM)c
Model3.2690.369.090.0041
Quadratic3.0531.0225.550.0004
A0.0010.000.020.8793
B0.0110.010.190.6725
C0.0910.092.290.1736
A23.0213.0275.99<0.0001
B20.0010.000.120.7379
C20.0110.010.300.6010
AB0.0210.020.540.4865
AC0.0810.081.990.2010
BC0.0010.000.100.7651
Residual0.2870.04
Lack of Fit0.1230.041.080.4533
Total3.5316
(b) Malaysia Lotus Root
SourceSum of SquaresDFMean SquareF-Valuep-Value
TFC (mg PCE/g)d
Model0.3590.03855.79<0.0001
Quadratic0.2130.069100.21<0.0001
A0.0310.02941.830.0003
B0.0010.003.560.1012
C0.1010.097140.58<0.0001
A20.1910.20284.64<0.0001
B20.0010.000.200.6666
C20.0110.007.820.0267
AB0.0110.008.170.0244
AC0.0010.002.940.1301
BC0.0010.004.390.0743
Residual0.0070.00
Lack of Fit0.0030.004.400.0930
Total0.3516
%DPPHsce
Model448.98949.8914.020.0011
Quadratic84.57328.197.920.0119
A187.891187.8952.790.0002
B23.15123.156.510.0381
C133.011133.0137.370.0005
A276.69176.6921.550.0024
B20.5710.570.160.7015
C25.4715.471.540.2549
AB7.8717.872.210.1807
AC0.02210.0220.000.9389
BC12.46112.463.500.1035
Residual24.9173.56
Lack of Fit16.2235.412.490.1988
Total473.8316
FRAP (mM)f
Model0.9290.106.120.0130
Quadratic0.5130.1710.120.0061
A0.1410.148.420.0229
B0.0010.000.220.6559
C0.2610.2615.750.0054
A20.4810.4828.910.0010
B20.0010.000.000.9511
C20.01310.0130.790.4032
AB0.0010.000.180.6831
AC0.0010.000.0730.7942
BC0.0010.000.0540.8230
Residual0.1270.017
Lack of Fit0.0130.000.0710.9723
Total1.0416
Table 2. BBD with the observed responses and predicted values for TFC (mg PCE/g), %DPPHsc and FRAP (mM) for (a) China lotus root and (b) Malaysia lotus root.
Table 2. BBD with the observed responses and predicted values for TFC (mg PCE/g), %DPPHsc and FRAP (mM) for (a) China lotus root and (b) Malaysia lotus root.
RunVariable Levels(a) China Lotus Root(b) Malaysia Lotus Root
Observed (Y1) aPredicted (Y0)Observed (Y1) aPredicted (Y0)
X1 (pH)X2
(time)
X3 (ratio)TFC
(mg PCE/g)
%DPPHscFRAP (mM)TFC (mg PCE/g)%DPPHscFRAP (mM)TFC (mg PCE/g)%DPPHscFRAP (mM)TFC (mg PCE/g)%DPPHscFRAP (mM)
12.50.5400.6650.52.200.6048.182.060.5830.22 1.84 0.5528.58 1.88
22.50.5200.4341.41.890.4241.631.790.3823.54 1.48 0.3923.95 1.49
32.51300.6245.72.100.5945.461.870.5226.41 1.72 0.4925.34 1.71
431400.3418.61.310.3620.581.320.2618.83 1.36 0.2719.09 1.35
531200.107.20.860.076.600.820.1212.88 1.00 0.111.08 1.02
620.530010.20.860.0311.910.960.3126.74 1.52 0.3126.58 1.50
721.5300.1515.61.180.1115.271.050.3727.37 1.49 0.3525.98 1.51
82.51300.5545.21.880.5945.461.870.4825.67 1.54 0.4925.34 1.71
931.5300.2310.60.980.208.890.880.1613.33 1.18 0.1613.49 1.20
102.51.5200.4438.31.640.5040.651.790.2715.38 1.51 0.317.02 1.47
112.51300.5844.41.720.5945.461.870.4925.70 1.76 0.4925.34 1.71
1221200.1418.11.130.1116.061.120.1820.88 1.24 0.1720.62 1.25
132.51300.6145.52.030.5945.461.870.526.15 1.96 0.4925.34 1.71
1421400.1514.61.020.1815.141.060.4127.13 1.67 0.4428.93 1.65
152.51300.5846.51.640.5945.461.870.4822.76 1.58 0.4925.34 1.71
162.51.5400.6847.41.840.6947.141.940.5829.12 1.81 0.5728.71 1.81
1730.5300.0413.90.950.0814.271.090.2518.31 1.32 0.2719.70 1.29
a Mean of triplicate determination.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tan, S.-J.; Lee, C.-K.; Gan, C.-Y.; Olalere, O.A. Statistical Optimization of Flavonoid and Antioxidant Recovery from Macerated Chinese and Malaysian Lotus Root (Nelumbo nucifera) Using Response Surface Methodology. Molecules 2021, 26, 2014. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26072014

AMA Style

Tan S-J, Lee C-K, Gan C-Y, Olalere OA. Statistical Optimization of Flavonoid and Antioxidant Recovery from Macerated Chinese and Malaysian Lotus Root (Nelumbo nucifera) Using Response Surface Methodology. Molecules. 2021; 26(7):2014. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26072014

Chicago/Turabian Style

Tan, Sze-Jack, Chee-Keong Lee, Chee-Yuen Gan, and Olusegun Abayomi Olalere. 2021. "Statistical Optimization of Flavonoid and Antioxidant Recovery from Macerated Chinese and Malaysian Lotus Root (Nelumbo nucifera) Using Response Surface Methodology" Molecules 26, no. 7: 2014. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules26072014

Article Metrics

Back to TopTop