Next Article in Journal
Effect of Olive Pomace Extract Application and Packaging Material on the Preservation of Fresh-Cut Royal Gala Apples
Previous Article in Journal
Physicochemical Properties, Antioxidant Capacity, and Bioavailability of Laurus nobilis L. Leaf Polyphenolic Extracts Microencapsulated by Spray Drying
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling and Optimization of Ultrasound-Assisted Extraction of Bioactive Compounds from Allium sativum Leaves Using Response Surface Methodology and Artificial Neural Network Coupled with Genetic Algorithm

1
Department of Food Process Engineering, National Institute of Technology, Rourkela 769008, Odisha, India
2
Department of Food Engineering and Technology, Sant Longowal Institute of Engineering and Technology, Longowal 148106, Punjab, India
*
Authors to whom correspondence should be addressed.
Foods 2023, 12(9), 1925; https://doi.org/10.3390/foods12091925
Submission received: 11 April 2023 / Revised: 20 April 2023 / Accepted: 5 May 2023 / Published: 8 May 2023
(This article belongs to the Section Food Packaging and Preservation)

Abstract

:
This study explains the effect of ultrasound on the extraction of the bioactive compounds from garlic (Allium sativum L.) leaf powder. The experiment was carried out by varying the ultrasound amplitude (30–60%), treatment time (5–15 min), and ethanol concentration (40–60%) required to obtain the maximum extraction yield of total phenol content (TPC), total flavonoid content (TFC), and antioxidant activity. Rotatable central composite design (RCCD) provided experimental parameter combinations in the ultrasound-assisted extraction (UAE) of garlic leaf powder. The values of extraction yield, TPC, TFC, and antioxidant activity for the optimized condition of RSM were obtained at 53% amplitude, 13 min of treatment time, and 50% ethanol concentration. The values of the target compounds predicted at this optimized condition from RSM were 32.2% extraction yield, 9.9 mg GAE/g TPC, 6.8 mg QE/g TFC, and 58% antioxidant activity. The ANN-GA optimized condition for the leaf extracts was obtained at 60% amplitude, 13 min treatment time, and 53% ethanol concentration. The predicted values of optimized condition obtained by ANN-GA were recorded as 32.1738% extraction yield and 9.8661 mg GAE/g, 6.8398 mg QE/g, and 58.5527% for TPC, TFC, and antioxidant activity, respectively. The matured leaves of garlic, if not harvested during its cultivation, often go waste despite being rich in antioxidants and phenolic compounds. With the increased demand for the production of value-added products, the extraction of the bioactive compounds from garlic leaves can resolve waste management and potential health issues without affecting the crop yield through the process for high-end use in value addition.

Graphical Abstract

1. Introduction

The agro-food sector plays a strategic role in the commercial and economic operations of producing and cultivating Allium species. Allium is one of the most produced and consumed vegetative products worldwide, with a world production of 32.4 MMT in the year 2020, with China being the highest producer, followed by India [1].
Recently, the demand for natural compounds in the prevention and cure of various diseases such as asthma, cholera, fever, diabetes, hypertension, and cancer is constantly increasing [2,3]. Besides that, consumers are becoming more health conscious and more aware of the ingredients in food products and are becoming inclined towards natural food additives rather than synthetic ones. This shift is causing an increase in the demand for natural compounds in the industrial processing of food. Allium species is not only a popular vegetable but also has therapeutic values [4]. It has been used in traditional medicinal practices since ancient times, attributed mainly to sulfur-containing compounds [5,6]. Allium sativum L. belongs to the onion family Alliaceae, having volatile compounds, sulfur, antioxidants, and many bioactive compounds [7]. With the plant growth, the older matured leaves, if not harvested, are dried on the plant and go to waste. Harvesting those leaves helps in the protection of plants from lodging and provides economic importance if harvesting occurs at proper intervals during their cultivation. Leaves of A. sativum contain higher levels of bioactive compounds, i.e., polyphenols, ascorbic acid, and dietary fiber, and exhibit more antioxidant activity [8]. Similar results were observed for the aged leaf extracts from garlic for antioxidant activity, followed by the flower and then for the bulb [9]. The presence of antioxidants in food can prevent cell damage by scavenging free radicals, thus preventing various ailments [3,4,10]. Furthermore, more than 20 phenolic compounds are present in garlic, which is a higher number than that in most vegetables [11]. β-resorcylic acid is the major phenolic compound in garlic, accompanied by pyrogallol, gallic acid, rutin, protocatechuic acid, and quercetin [12,13].
The bioactive compounds mainly find application in the development of newer food products, medicine, agriculture, and industry but also are used in the green synthesis of metal nanoparticles [14,15,16,17,18,19]. In the green synthesis of metal nanoparticles, a bioactive compound can serve as a reducing agent by donating electrons and converting metal ions into metal nanoparticles (MNPs) and as a stabilizing agent by preventing the agglomeration of the MNPs [14,15,19]. The bioactive compounds in the extract form are mainly used to synthesize the MNPs [20]. The bioactive compounds present in garlic can be derived in the form of an extract by using different conventional or novel extraction methods. The conventional extraction techniques used have lower efficiency than novel techniques of extraction such as microwave-assisted extraction (MAE) and ultrasound-assisted extraction (UAE) [21,22,23,24]. Unlike high-power ultrasound applied in component extraction, low-power ultrasound is often used as a non-destructive tool in food component analysis and medical imaging [25,26,27]. The extraction using UAE is relatively more suitable than other methods due to its cost-effectiveness and lower instrumental requirements [28]. High-power ultrasound is used to achieve increased extraction efficiency due to the phenomena of cavitation bubbles. These cavitation bubbles compress and decompress with ultrasonic waves and finally collapse, producing a shock wave. These shock waves can enhance the mixing and mechanical effect, causing increased solvent diffusion in the product matrix and subsequently increasing the release of intracellular compounds to the bulk medium [29,30]. To obtain the maximum yield of the intracellular compound and superior extraction, the optimum values of the factors influencing the extract must be determined [31]. Prediction tools such as RSM and ANN-GN are primarily utilized to determine the optimal conditions because of their ability to model complex relationships between variables and predict outcomes based on those relationships [32].
Response surface methodology (RSM) is a mathematical and statistical method used to elucidate the relationship between multiple input variables (also known as factors) and a single output variable (also known as the response). It is often used in engineering and scientific research to optimize the design of a product or process to identify the optimal combination of input variables that will result in the desired output [33,34,35]. The mathematical model that explains the input variables and the response variable is used to make predictions and to identify the optimal combinations of input variables [36,37]. In contrast, an artificial neural network (ANN) is a machine learning method inspired by the structure and function of the human brain. The structure comprises neurons that receive, process, and transmit information [37]. Based on the input data, an ANN is designed to recognize patterns and make predictions. It can be used as a prediction tool even when the data are complex and have a non-linear relationship which requires high levels of prediction accuracy [38,39]. To evaluate the optimized condition from the predicted parameters, a genetic algorithm (GA) is applied.
A GA is an optimization technique that is built on the natural selection and genetics principle. It comprises the application of genetic operators such as selection, crossover, and mutation to create a new generation of candidate solutions [40]. The main aim of a genetic algorithm is to find the optimum solution to the problem by mimicking the natural selection process [41]. An ANN combined with a GA is a hybrid approach to machine learning that combines the capabilities of both the ANN and GA. Combining both allows the incorporation of a combined approach and strengthens prediction and optimization. For the analysis and processing of data, a neural network can be used, whereas for optimizing the parameters, a genetic algorithm can be used [42].
Thus, the aim of this study was to assess and compare the effectiveness of two popular optimization techniques, RSM and ANN-GA, in the modeling and optimization of complex systems for the production of quality garlic (A. sativum) leaf powder extract. Moreover, the effect of different extraction parameters on the properties of the extract was also investigated.

2. Materials and Methods

2.1. Reagents and Standards

To estimate the bioactive compounds in the sample, chemicals such as quercetin dihydrate, Folin–Ciocalteu reagent, gallic acid, DPPH, and other chemicals were procured from Hi-Media Laboratories Pvt. Ltd., Mumbai, India. All the glassware used was from Borosil Glass Works Limited, Mumbai, India.

2.2. Raw Material

The two most matured green leaves from each A. sativum plant leaving four well-developed leaves on the plant of variety Haryana Garlic-17 (HG-17) were harvested in the early morning from a locally grown progressive farmer’s field near the vicinity of SLIET. Longowal Campus, Sangrur (Punjab, India). Before use, the hay-like dried portion at the distal ends of harvested green leaves was removed and then washed in running tap water twice. After apparent surface water drying, the leaves were cut lengthwise into approximately 1.5 cm and dried using a PID-controlled cabinet drier maintained at 50 ± 1 °C. The dried garlic leaves were ground and sieved using a sieve of 100 mesh to obtain garlic leaf powder (Figure 1).

2.3. Ultrasound-Assisted Extraction (UAE)

UAE was carried out using a 20 kHz probe ultrasonicator (Qsonica 700, Qsonica, Newtown, CT, USA) with a probe diameter of 16 mm. A frequency of 20 kHz with a constant power of 700 W was used for the extraction of the target compounds. The temperature of the samples was maintained by keeping them in an ice bath, causing the temperature not to exceed 40 °C. For the experimental runs, the garlic leaf powder was mixed with an ethanol–water mixture and placed in a beaker, maintaining the same distance from the probe for each run. The applied energy for each treatment was recorded as the instrument’s energy reading (J). Calorimetric energy was calculated as per Kikuchi and Uchida [43] for each experiment. The digital panel on the instrument was used to control the ultrasound amplitude and treatment time. After applying the ultrasound, the extract obtained was immediately filtered using the Whatman filter paper 42. The obtained extracts were kept at 4 °C for further analysis.

2.4. One Factor at a Time

The one factor at a time (OFAT) approach was used to identify and select critical factors among the various independent variables affecting the extraction yield. The various experimental variables and their levels selected via OFAT experiments were ultrasonic amplitude (20%, 30%, 40%, 50%, and 60%), ethanol concentration (20%, 40%, 50%, 60%, 80%, and 100%), treatment time (5, 10, 15, and 20 min) and solvent-to-solid ratio (10:1, 20:1, 30:1, 40:1, and 50:1). During this approach only a single factor was varied, keeping other variables constant. The total phenol content of different experiments was evaluated for the selection of independent variables that significantly influenced the extraction efficacy.

2.5. Determination of Proximate Composition

Moisture, protein, fat, ash content, and pH of the garlic leaves were estimated using the standard method [44]. The carbohydrate content of garlic leaves was determined using the difference method.

2.6. Determination of Total Phenol Content

For determining the total phenol content, the Folin–Ciocalteu method was followed [45] with slight modification. Initially, 0.1 mL aliquot of the extract was combined with 2.5 mL Folin–Ciocalteu reagent. The mixture was then incubated for 30 min after adding 2 mL of sodium carbonate (10% w/v). The absorbance of the standard and the reaction mixture was observed using a spectrophotometer (UV-1800 Shimadzu spectrophotometer, Shimadzu, Kyoto, Japan) at 725 nm.

2.7. Determination of Total Flavonoid Content

Total flavonoid content was determined using the aluminum chloride colorimetry assay method [46]. As the standard in this procedure, quercetin was utilized. Initially, 0.5 mL of supernatant was mixed with 4 mL of distilled water. Then, 0.3 mL of 5% sodium nitrite was added to the above mixture. After 5 min, a mixture of 0.3 mL of 10% aluminum chloride was added. At the 6th minute, 2 mL of 1 M sodium hydroxide was added, and 2.4 mL of distilled water was mixed properly. The mixture was kept at room temperature for 15 min, and then absorbance was measured at 510 nm using a spectrophotometer. The calibration curve of quercetin was drawn from 0–1 mg/mL, and the results were expressed as quercetin equivalent.

2.8. DPPH Assay

The antioxidant activity was evaluated using a modified version of the DPPH scavenging assay as proposed in [47]. For this, 0.1 mL of garlic leaf extract was mixed with 3.9 mL of 1 mM DPPH solution. The mixture was maintained in the dark at room temperature for thirty minutes. The absorbance of the solution was then evaluated with a spectrophotometer at 515 nm. For the preparation of the control, instead of the sample, 0.1 mL of distilled water was combined with DPPH solution (3.9 mL) and observed at the same wavelength. Equation (1) was used to calculate the DPPH scavenging activity.
DPPH   scavenging   % = A S A C A C × 100
where, AC and AS are the absorbances of the control and extract, respectively.

2.9. Estimation of Organosulfur and Phenolic Compounds Using HPLC

Organosulfur and phenolic compounds were estimated by using a chromatographic estimation technique. It was performed using Water’s HPLC 2489 equipped with Water’s UV detector (Waters Corporation, Milford, MA, USA). The Luna LC column (150 × 4.6 mm, 5 μ particle size) and C18 column (Phenomix 100 × 4.6 mm, 5 μ particle size) were used for the estimation of organosulfur compounds (OSCs) and phenolic compounds, respectively.
OSCs were detected with a gradient condition of 50% mobile phase A (acetonitrile) and 50% mobile phase B (water) with a constant flow rate of 1.0 mL/min at a detection wavelength of 254 nm. Phenolic compounds were detected at 280 nm wavelength with varying concentrations of solvent A (water/acetonitrile/formic acid = 96/3.9/1; v/v) and solvent B (acetonitrile/formic acid = 99.9/0.1; v/v) with time, at the flow rate of 1.0 mL/min. The gradient program was performed with varying percentage proportions of solvent A and solvent B (at the initial stage, the concentration of solvent A:solvent B was 100:0, at 5 min 95:5, 25 min 85:15, 30 min 80:20, 39 min 75:25, 43 min 55:45, 48 min 5:95, 50 min 5:95, 55 min 80:20, and 60 min 100:0).
All the prepared solvents were ultrasonicated for 7–10 min to remove the entrapped air before the analysis. The extracted samples were filtered through a 0.45 μ nylon filter to eliminate the undesirable contaminants if present. An individual OSC and phenolic standard was first run to estimate the faithful retention time at respective wavelengths. These compounds were identified by comparison with standard compound retention times. The individual OSCs and phenolic compounds were calculated by applying Equations (2) and (3), respectively.
Concentration   of   individual   OSC   p p m = A r e a s a m p l e A r e a S t d × C o n c s t d
Concentration   of   individual   phenolic   compound   p p m = A r e a s a m p l e A r e a S t d × C o n c s t d

2.10. Experimental Design

The input parameters and their ranges were designated based on the method of one factor at a time. The solute-to-solvent ratio of 1:40 was kept constant throughout the experiment. Amplitude (20–60%), time (5–15 min), and ethanol concentration (40–60%) were chosen as independent variables to optimize the extraction yield, total phenol content, total flavonoid content, and antioxidant activity as the dependent parameters. Rotatable central composite design (RCCD) was used to obtain the different compositions of input variables and analyzed using design expert software (version 11.0.5.0, Stat-ease Inc., Minneapolis, MN, USA). According to RCCD, a total of 20 runs consisted of experimental combinations of eight factorial, six axial, and six central points [48] as presented in Table 1 and Table 2.

2.11. Fourier Transform Infrared Spectroscopy Analysis

Fourier transform infrared spectroscopy (FTIR) is an analytical technique that provides information on the functional groups present in the extract, which helps to identify and quantify the chemical components present in the sample. A Bruker Alpha FTIR system (Germany) was used to determine the chemical conformation of the garlic leaf powder extract obtained using UAE. The samples were scanned over a range of 4000–500 cm−1 at a resolution of 4 cm−1 with 64 scans per sample.

2.12. Field Emission Scanning Electron Microscopy (FE-SEM)

Field emission scanning electron microscopy (FE-SEM) is a powerful imaging technique used to visualize the pores or microcracks present on the surface structure of the sample at higher resolution. The structural analysis of the garlic leaf powder with and without ultrasonic treatment was observed using a field emission scanning electron microscope (FEI Novanano SEM 450, Nebraska Center for Materials and Nanoscience, Lincoln, NE, USA). A magnification of 2500× was used to acquire the garlic leaf powder image before and after the application of ultrasound treatment.

2.13. Response Surface Methodology (RSM)

RSM was utilized to model and optimize the process parameters to optimize the dependent parameter values. The response variables were fitted to the second-order polynomial model as given in Equation (4). The experimental data were analyzed using a second-order full polynomial regression model to predict optimal conditions.
Y = β 0 + i = 1 k β i X i + i = 1 k β i i X i 2 + i k 1 j k β i j X i X j
where Y represents the response variable; Xi and Xj are actual values in coded form of the independent variable; k represents the number of considered factors (k = 3); and β 0 ,  β i , β j j , and β i j are the constant and coefficients of linear, quadratic, and interaction terms, respectively.
Analysis of variance (ANOVA) was used to determine the significance of the regression coefficient, and the developed regression models were validated using the conducted statistical analysis. The adequacy of the model was evaluated by the coefficient of multiple determination (R2). The terms in the polynomial were considered statistically different when p < 0.05.

2.14. Artificial Neural Network (ANN) Modeling Coupled with Genetic Algorithm (GA)

The artificial neural network and the genetic algorithm were used for modeling and optimization of parameters by neural net fitting application in machine learning and the optimization tool in maths, statistics, and optimization tool, respectively, of MATLAB software (Version R2018a, The MathWorks, Inc., Natick, MA, USA). The neural network used for the experiment comprised one input, one hidden, and one output layer. Three neurons were present in the input layer, i.e., ultrasound amplitude (X1), treatment time (X2), and ethanol concentration (X3). A single hidden layer was selected for the experiment to prevent the issue of overfitting. The number of neurons in the hidden layer was decided based on the trial and error method, i.e., the heuristic method [49] to obtain the larger R2 and lower MSE values; based on that, 10 neurons were selected for the hidden layer. The training for each response (TPC (Y1), TFC (Y2), Antioxidant (Y3), and Extraction yield (Y4)) was performed separately, and hence the output layer comprised only a single neuron. RCCD was used to build the experimental design for UAE of garlic leaf powder with the above three input parameters. The Levenberg–Marquardt backpropagation algorithm was used for training the algorithm as it is fast and most recommended in the toolbox. The input–output datasets were divided into three subsets, namely training, validation, and testing, consisting of 70%, 15%, and 15%, respectively. The transfer function for the hidden layer was the hyperbolic sigmoid function (tansig), whereas the transfer function for the output layer was linear (purelin). The network was trained continuously until the maximum R2 value and minimum MSE value were attained. After training, bias and the weight values of the generated model were obtained, and Equation (5) was utilized for the prediction of the response values (Yi).
t a n s i g n = 2 1 + e 2 × n 1
p u r e l i n n = n
Y i = p u r e l i n ( W O H × t a n s i g ( U I H × X i + T H ) + T O )
where n is the net input, Yi is the predicted output parameter, and Xi is the input parameter. UIH is the weight between the input layer and hidden layers, and WOH is the weight between the hidden layer and output layer. TH is the bias value of hidden layer neurons, and TO is the bias value of output.
The GA was used to optimize the input variables of the developed ANN model. During the optimization phase, the fitness function (f) was employed to maximize all of the dependent parameters. The migration, crossover function, mutation function, selection function, scaling function, creation function, and population type for the genetic algorithm were chosen as forward migration, scattered, adaptive feasibility, roulette function, rank, feasible population, and double vector. For the best optimization result, a 0.8 value was used for the crossover fraction, and the remaining functions were set to default. The ANN-GA optimization fitness function (f) is as follows:
f = ( Y 1 + Y 2 + Y 3 + Y 4 )
where Y1, Y2, Y3, and Y4 are the ANN predicted extraction yield, TPC, TFC, and antioxidant response, respectively, and the negative sign indicates maximizing of the function f in the GA.

2.15. Statistical Analysis

All experiments were performed at least in triplicates; one-way analysis of variance (ANOVA) was used to examine the data, and SPSS (version 25) statistic software was used for the analysis. Duncan’s multiple range test was used to separate mean values, and significant differences were considered at p < 0.05. Data are presented as mean ± standard deviation (SD). The results were analyzed using Design Expert (version 11.0.5.0, Stat-ease Inc.), and ANOVA was used to determine the model’s significant terms and regression coefficients of linear, interaction, and quadratic terms.
Statistical parameters such as coefficient of multiple determination (R2), normal root mean square (NRMSE), root mean square error (RSME), mean percentage error (MPE), normal mean square error (NMSE), mean square error (MSE), and average absolute deviation (AAD) were used for the comparison of the performance of models built with ANN and RSM. The optimum model for expressing the responses has the highest R2 and lowest NRMSE, RSME, MPE, AAD, and MSE. The formulae of R2, NRMSE, RSME, MPE, AAD, and MSE are presented as follows:
MSE = x p x a 2 n
NMSE = MSE x m
AAD = x p x a n
MPE = 100 n | x p x a x p |
RMSE = x p x a 2 n
NRMSE = RMSE x m
R 2 = 1 x p x a 2 x p x m 2
where xp, xa, xm, and n represent predicted data, experimental data, experimental mean data, and number of experiments, respectively.

3. Results and Discussion

3.1. Proximate Analysis

The chemical parameters of fresh garlic leaves (variety: HG-17) were analyzed for the proximate composition and pH. The moisture content of garlic was found to be 90.04%, which is in the range of other green leafy vegetables [50]. Protein, fat, carbohydrate, ash content, and pH were found to be 2.23%, 0.18%, 5.99%, 1.56%, and 6.08%, respectively (Table 3). Protein content in the garlic leaves was found to be considerably higher and fat content was lower than that in other leafy vegetables, namely Cardiospermum halicacabum, Premna latifolia, Argyreia pomacea, Mollugo pentaphylla, and Pisonia grandis [51].

3.2. One Factor at a Time (OFAT)

The range for the independent variable was determined using the OFAT approach. The range of the independent variable was selected in such a way that it covers the expected operating range and is not wide enough to make the experimental design impractical and complex. Preliminary screening of the independent variable was performed firstly by keeping all the variables constant and varying only the solute-to-solvent ratio. The solid-to-solvent ratio was kept constant at 1:40 as there was no significant effect on the extraction efficiency. After keeping the solid-to-solvent ratio constant, other independent variables were adjusted in such a way as to include the optimal values for each factor, i.e., ultrasound amplitude, treatment time, and ethanol concentration. These identified optimal values of selected critical variables were considered as the respective levels of the central point in the CCD of RSM (Table 1). The optimal values identified in the OFAT approach for ultrasonic amplitude, treatment time, and ethanol concentration were 45%, 10 min, and 50%, respectively.

3.3. RSM Analysis

3.3.1. Model Fitting

The experimental and predicted values of all the dependent variables (extraction yield, TPC, TFC, and antioxidant content) and the independent parameters (ultrasound amplitude, time, and ethanol concentration) along with the applied energy and calorimetric energy are presented in Table 2. Similarly, the equations below represent the second-order full polynomial models with regression coefficients used to predict the dependent variables.
EY = 31.91 + 0.4514X1 + 0.5161X2 + 0.3259X3 − 0.0950X1X2 − 0.2875X1X3 − 0.5175X2X3 − 0.3482X12 − 0.3482X22 − 0.3058X32
TPC = 9.85 + 0.1194X1 + 0.1052X2 + 0.0735X3 + 0.0163X1X2 + 0.0187X1X3 − 0.0662X2X3 − 0.0794X12 − 0.0511X22 − 0.0740X32
TFC = 6.71 + 0.4526X1 + 0.2183X2 + 0.2596X3 − 0.0188X1X2 + 0.0337X1X3 − 0.0188X2X3 − 0.1614X12 − 0.1985X22 − 0.2038X32
Antioxidant activity = 57.84 + 1.01X1 + 0.6629X2 + 1.35X3 + 0.1412X1X2 − 0.3113X1X3 − 0.3438X2X3 − 0.6993X12 − 0.6481X22 − 0.6092X32
where X1, X2, and X3 are independent variables in coded form for ultrasound amplitude, treatment time, and ethanol concentration, respectively.
The regression coefficients, R2, adj. R2, and lack of fit of the obtained second-order polynomial equation for the various responses of garlic leaf extract treated with ultrasound are shown in Table 4. Table 5 shows different parameters such as R2, NRMSE, RSME, MPE, AAD, and MSE using formulae given in Equations (7)–(13).
According to Table 4, the R2 and adj. R2 values of EY, TPC, TFC, and antioxidant activity were 0.96 and 0.93, 0.96 and 0.92, 0.93 and 0.87, and 0.91 and 0.84, respectively, indicating that the generated model is desirable. In addition, the lower values of statistical parameters for the responses (Table 4) verify the model’s desirability. All four responses exhibited a non-significant lack of fit, with values of 0.660 (EY), 0.336 (TPC), 0.051 (TFC), and 0.109 (antioxidant), indicating that the models fit the data well.

3.3.2. Effect of Ultrasound Treatment on Total Phenol Content

Table 2 shows the mean response value of total phenol content (TPC) of garlic leaf extract for 20 different experimental runs. The highest TPC value, 9.9 mg GAE/g, was obtained at an amplitude of 45%, treatment time of 10 min, and ethanol concentration of 50%, whereas the lowest value, 9.38 mg GAE/g, was obtained at an amplitude of 30%, treatment time of 5 min and ethanol concentration of 40%. Ultrasound amplitude (X1) (p < 0.001), treatment time (X2) (p < 0.001), and ethanol concentration (X3) (p < 0.001) exhibited a linear significant positive effect towards extraction of TPC. Similarly, the interaction term X2X3 (p < 0.05) and each quadratic term substantially impacted the TPC yield.
The three-dimensional surface plots shown in Figure 2a–c describe the interaction between amplitude (X1) and treatment time (X2). TPC yield increases as amplitude and treatment time increase; this may be due to increased cavitation and thermal effect at higher ultrasonic amplitude, which induces disruption and softening of the cell walls [52]. However, in the case of interaction between ethanol concentration (X3) and treatment time (X2), total phenol content yield initially increased with an increase in ethanol concentration but subsequently decreased; this could be due to the disintegration of the cell wall as water acts as a swelling agent breaking up the plant cell wall. The ethanol dissolves the polyphenol by disrupting the solvent solid contact area, but a further increase in the ethanol concentration leads to an increase in the dielectric constant of the solvent, which causes a further decrease in the yield of TPC [53]. Shirzad et al. [54] reported a similar trend for extracting TPC from olive leaves using ultrasound.

3.3.3. Effect of Ultrasound Treatment on Total Flavonoid Content (TFC)

The highest TFC value of 6.96 QE/g was obtained at an amplitude of 60%, treatment time of 15 min, and ethanol concentration of 60%, whereas the lowest value, 5.39 QE/g, was obtained at an amplitude of 30%, treatment time of 5 min, and ethanol concentration of 40%. Antioxidant activity was linearly, significantly, and positively affected by ultrasound amplitude (X1) (p < 0.005), treatment time (X2) (p < 0.001), and ethanol concentration (X3) (p < 0.001) and negatively affected by the quadratic terms.
The three-dimensional surface plot represented in Figure 2 shows the interaction between amplitude (X1) and time (X2), which shows an increase in the value of TFC with the increase in the applied amplitude and time (up to about 12 min). The increase in the phenolic compounds can be due to the cavitation effect caused by US disruption of the structure of leaves, which aided in the diffusion of phenolic compounds in the solvent [55]. Further increase in the time caused a reduction in the TFC value. A similar trend was reported by Zimare et al. [56], where ultrasound treatment beyond 10 min decreased the TFC content of Lobelia nicotianifolia leaf extracts. The reduction in the bioactive compounds with treatment time can be attributed to the formation of free radicals during cavitation [57].

3.3.4. Effect of Ultrasound Treatment on Antioxidant Activity

The DPPH (2,2-diphenyl-1-picrylhydrazyl) scavenging assay is a commonly used method for measuring the antioxidant activity of a substance. The mean values of the obtained experimental data and ANOVA analysis are shown in Table 2 and Table 4, respectively. The highest antioxidant value obtained was 58.59% at amplitude, treatment time, and ethanol concentration of 45%, 10 min, and 50%, respectively, whereas the lowest value obtained was 53.12% at amplitude, treatment time, and ethanol concentration of 45%, 10 min, and 33%, respectively. Antioxidant activity was positively affected by ultrasound amplitude (X1) (p < 0.001), treatment time (X2) (p < 0.01), and ethanol concentration (X3) (p < 0.001) and negatively affected by the quadratic terms.
The three-dimensional plot in Figure 2 represents the interaction between amplitude (X1) and time (X2), demonstrating an increase in antioxidant activity with an increase in amplitude and time. The antioxidant values increased up to a certain ultrasonic amplitude of 55%. The surge in the antioxidant activity might be due to cavitation that increases the thermal effect leading to disruption of the structure of plant cells that leads to the release of previously bound antioxidants within the cell [52]. Moreover, increased antioxidant activity can also be attributed to increased polyphenol content in garlic leaf extract due to cavitation generated during ultrasonication [58]. The further decrease in antioxidant activity is due to the detrimental effect of oxidation during prolonged ultrasonication [59]. A similar trend for antioxidant activity in cashew apple bagasse was reported [31].

3.3.5. Effect of Ultrasound Treatment on Extraction Yield

The mean value of extraction yield is presented in Table 2. The highest value of 32.24% was obtained at amplitude, treatment time, and ethanol concentration of 60%, 15 min, and 40%, respectively, whereas the lowest extraction yield of 28.71% was obtained at an amplitude of 30%, treatment time of 5 min, and ethanol concentration of 40%. The ANOVA analysis of the data represents a positive effect of linear terms, i.e., ultrasound amplitude (X1) (p < 0.0001), treatment time (X2) (p < 0.0001), and ethanol concentration (X3) (p < 0.0005), and negative effect on interaction terms such as the amplitude–ethanol concentration (X1X3) (p < 0.01) and treatment time–ethanol concentration (X2X3) (p < 0.0001) and quadratic terms.
Visual analysis of the three-dimensional surface and contour plots show how different parameters affect the extraction yield. The plot showing the interaction between amplitude (X1) and treatment time (X2) with ethanol concentration (X3) shows that there is an increase in the extraction yield with an increase in ethanol concentration to 53%, but further decrease can be observed due to the negative effect of the quadratic term on the ethanol concentration. The decrease in the yield can be due to the fact that at higher ethanol concentrations, dehydration of plant tissue occurs, which results in a decrease in the yield at higher concentrations [60]. A similar trend was reported by Tomšik, Pavlić, Vladić, Ramić, Brindza, and Vidović [29] for bioactive compound extraction using Allium ursinum L.

3.4. Artificial Neural Network Analysis

Model Fitting

The Levenberg–Marquardt (LM) algorithm was used for the prediction of TPC, TFC, antioxidant, and extraction yield. Maximum R2 and minimum MSE values give reliability and applicability to the developed model in prediction adequacy. Accordingly, the best model was selected based on the maximum R2 and minimum MSE values. The final model that was selected for TPC, TFC, and antioxidant activity was obtained at epoch 5, and for the extraction yield, epoch 6. The model had one layer of input (3 neurons), one hidden layer (10 neurons), and one layer of output for each response. The neural network diagram for the experiment is depicted in Figure 3.
The values of MSE for training, testing, and validation of TPC were 0.53554, 0.45408, and 0.38106, and the values of R2 for the same were 0.984, 0.987, and 0.988, respectively. Similarly, for TFC, the values of MSE for training, testing, and validation were 0.33468, 0.66318, and 0.931, and the values of R2 was 0.987, 0.994, and 0.994, respectively. Further, the values of MSE for training, testing, and validation of antioxidant activity were 0.62286, 0.8568, and 0.13710, and the values of R2 for the same were 0.962, 0.996, and 0.884, respectively. Similarly, the MSE values for training, testing, and validation of extraction yield were 0.18343, 0.77744, and 0.12273, and the values of R2 for the same were 0.967, 0.993, and 0.987, respectively. High R2 values confirm the prediction capability of the model. The optimized output data obtained from the model were compared with the target values, and the difference was used to determine the weights for the minimization of the individual sum of square errors (SSE). Figure 4 shows the error histogram and performance of the model.

3.5. Optimization Using Genetic Algorithm (GA)

Using GA, the process parameters (ultrasound amplitude, treatment time, and ethanol concentration) were tuned to maximize the extraction of TPC, TFC, antioxidant activity, and extraction yield. The procedure was repeated until the minimum RSME and MSE between the individual fitness values and mean value were achieved. After the mutation, the optimization cycle was maintained, and if the targeted outcome was not reached, the entire population was employed for reproduction, crossover, and mutation for the following cycle. Figure 5 illustrates the relation between the number of generations and the fitness value.
After the 63rd generation, the figure did not show any variation in the fitness value; a value of −107.413 was obtained as the mean fitness value. The optimization was accomplished under the following conditions: treatment period, 13 min; amplitude, 60%; and ethanol concentration, 53%, with applied energy of 26,150 J and calorimetric energy of 161.37 J. The expected results of EY, TPC, TFC, and antioxidant activity were obtained at this optimum condition as 32.17%, 9.87 mg GAE/g, 6.84 mg QE/g, and 58.55%, respectively.

3.6. Comparison between RSM and ANN Models

The predicted values of RSM and ANN models were statistically compared using statistical measures, such as R2, NRMSE, RSME, MPE, AAD, and MSE.
Statistical parameters of both models are shown in Table 5, from which it can be inferred that both models were entirely predictive. However, when compared, the ANN model’s overall R2 values were more significant than those of RSM. In addition, the ANN model has lower values for NRMSE, RSME, MPE, AAD, and MSE than the RSM model. These findings indicate that the artificial intelligence approach of the ANN model predicts responses far more accurately than the RSM model, the statistical approach in optimization.
The values of extraction yield, TPC, TFC, and antioxidant activity for the optimized condition of RSM were obtained at 53% amplitude, 13 min of treatment time, and 50% ethanol concentration with 22,430 J as applied energy and 163.94 J as calorimetric energy. The values of the target compounds predicted at this optimized condition were 32.2% extraction yield, 9.9 mg GAE/g TPC, 6.8 mg QE/g TFC, and 58% antioxidant activity. Whereas for the ANN-GA optimized condition of 60% amplitude, 13 min treatment time, and ethanol concentration of 53%, having applied energy of 26,150 J and calorimetric energy of 161.37 J, the expected values of EY, TPC, TFC, and antioxidant activity were 32.1738%, 9.8661 mg GAE/g, 6.8398 mg QE/g, and 58.5527%, respectively. The experimental values at these conditions showed that the values obtained at the optimum condition determined by ANN-GA were more significant than those obtained using RSM as ANN predicted more accurately as it had a higher R2 value and lower NRMSE, RSME, MPE, AAD, and MSE values than RSM.

3.7. FE-SEM Analysis

Figure 6 shows the image of the garlic leaf powder at 2500× magnification before and after the application of ultrasound treatment. It was observed that there was a significant change in the garlic leaf powder surface after the ultrasonic treatment application because of the formation of microcracks, pores, and cell wall fragmentation. These changes were caused due to cavitation on the surface-treated leaf extract samples and are depicted in Figure 6b. These microcracks, pores, and cell wall fragmentation assist the extraction of phenolic compounds present in the plant cell’s innermost part by increasing the material’s surface area and mass transfer of the target compounds to the solvent. Similar results were reported by AL-Bukhaiti et al. [61], showing that the application of ultrasound treatment on Cissus rotundifolia increased the phenol diffusion rate from the plant cell into the solvent. Therefore, the benefits of applying ultrasound may result from tissue and cell damage that would speed up the mass transfer of bioactive molecules such as polyphenols.

3.8. FTIR, Polyphenol, and Sulfur Compound Analysis

The infrared spectra of garlic leaf extract obtained at the optimum condition of UAE revealed the presence of functional groups that indicate a polyphenol-rich extract. The spectral region between 3570 and 3200 cm−1 represents stretching vibrations of the hydroxyl group (OH) (Figure 7) [62]. The bands between 3300 and 2500 cm−1 confirmed the presence of alkane (C-H stretch) and carboxylic acid (O-H stretch) [63]. The region of 1711–1511 cm−1 shows C=O stretching [64]. The peak in the region between 1600 and 1500 cm−1 could be due to aromatic nuclei vibration, and that between 1500 and 1400 cm−1 indicates the presence of C–C cm−1 bonds in the phenolic groups [65]. The spectral region between 1044 cm−1 and 1088 cm−1 shows the presence of an alcohol functional group [66].

3.9. Identification and Quantification of Organosulfur Compounds (OSCs)

Table 6 shows the concentrations of different identified organosulfur compounds as observed in the chromatogram of garlic leaf powder extract, which were confirmed as alliin, S-allyl-L-cysteine (SAC), and allicin. The concentrations of alliin, SAC, and allicin were estimated to be 90.207, 4.314, and 219.536 ppm, respectively. Allicin is the major bioactive compound formed on crushing of garlic and is decomposed at the injection port of GC and converted to vinyldithiins [67]; thus, an HPLC detection approach for the organosulfur compounds is preferred.

3.10. Identification and Quantification of Phenolic Compounds

Twelve phenolic compounds were identified in the HPLC chromatogram against the used known phenolic compound standard (Figure 8). The retention times of standard phenolic compounds were 4.225 min, 7.333 min, 13.532 min, 14.201 min, 16.479 min, 23.587 min, 26.421 min, 26.772 min, 30.292 min, 31.905 min, 42.509 min, and 50.179 min against the known concentrations of standard compounds, namely gallic acid, 3,4-dihydroxybenzoic acid, chlorogenic acid, catechin hydrate, syringic acid, p-coumaric acid, rutin, ellagic acid, benzoic acid, hesperidin, quercetin, and β-carotene, respectively. Table 6 shows the concentration of identified phenolic compounds present in the extracted sample. Among the identified phenolic compounds, β-carotene (117.607 ppm) showed the maximum concentration of identified phenolic compounds. Chlorogenic acid, 3,4-dihydroxybenzoic acid, catechin hydrate, rutin, benzoic acid, and gallic acid were found to have concentrations of 36.537 ppm, 34.403 ppm, 26.327 ppm, 21.741 ppm, 15.234 ppm, and 13.591 ppm, respectively. The result obtained showed a similar trend for the phenolic compounds for gallic acid, catechin hydrate, and p-coumaric acid and showed lower concentrations for the phenolic compounds syringic acid and benzoic acid, which were estimated in the aqueous extracts of garlic cultivars [68]. Syringic acid, p-coumaric acid, ellagic acid, hesperidin, and quercetin were found in concentrations of 4.276 ppm, 1.228 ppm, 0.968 ppm, 9.272 ppm, and 6.140 ppm, respectively.

4. Conclusions

Ultrasound-assisted garlic leaf powder bioactive component extraction was optimized using response surface methodology and an artificial neural network coupled with genetic algorithm (ANN–GA) approach. Analysis of variance data for assessing the impact of different independent parameters on the EY, TPC, TFC, and antioxidant content were analyzed, and the developed models were assessed for predicting the responses in accordance with the experimental results. The artificial neural network coupled with genetic algorithm approach was much more efficient in prediction and estimation capabilities even with fewer datasets as compared to RSM, making it a more reliable tool. The optimum values of the independent variables obtained using ANN-GA were 60% ultrasound amplitude, 13 min treatment time, and 53% ethanol concentration, having applied energy of 26,150 J and calorimetric energy of 161.37 J. The expected values for this optimized condition were obtained as 32.1738%, 9.8661 mg GAE/g, 6.8398 mg QE/g, and 58.5527% for EY, TPC, TFC, and antioxidant activity, respectively. The results also showed the potential of garlic leaf powder in the extraction of free radical scavenging characteristics. FTIR and FE-SEM further revealed the positive effect of applied high-power ultrasound in releasing active components, specifically polyphenols and sulfur-containing compounds. This paves the way for the high-end applications of developed garlic leaf powder and extract. Scaling up the developed process thus may have the potential to unlock a multitude of benefits in the exploration of bioactive compounds derived from garlic leaves, preferably at the industrial level.

Author Contributions

Conceptualization, S.S. and S.K.S.; methodology, S.S. and P.P.; software, S.S.; validation, S.S.; formal analysis, S.S.; investigation, S.S.; resources, P.P.; data curation, S.S. and P.P.; writing—original draft, S.S.; writing—review and editing, P.S., K.P. and S.K.S.; supervision, P.S., K.P. and S.K.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All research data has been reported in this manuscript.

Acknowledgments

Special thanks to Niveditha Asaithambi and Abhipriya Patra, PhD research scholars of National Institute of Technology (NIT) Rourkela for their invaluable assistance during the experiment and manuscript preparation. The authors are grateful and would like to thank the Department of Food Process Engineering, National Institute of Technology (NIT) Rourkela, Odisha, for the support provided during the experiment and manuscript preparation. The authors are also thankful to the Centralised Research Facility (CRF) of the National Institute of Technology Rourkela, India, for providing the facility for the use of analytical instruments such as a field emission scanning electron microscope, UV-vis spectrophotometer, and Fourier transform infrared spectrometer.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAO. Food and Agriculture Organization of the United Nations (FAO) 2021. Available online: http://www.fao.org/faostat/en/#data (accessed on 3 December 2022).
  2. Charron, C.S.; Milner, J.A.; Novotny, J.A. Garlic. In Encyclopedia of Food and Health; Academic Press: Cambridge, MA, USA, 2016; pp. 184–190. [Google Scholar]
  3. Singh, R.; Prasad, K.K.; Siddiqui, M.W.; Prasad, K. Medicinal Plants in Preventive and Curative Role for Various Ailments; Apple Academic Press: New York, NY, USA, 2017; pp. 63–100. [Google Scholar]
  4. Kaur, R.; Shekhar, S.; Prasad, K. Secondary Metabolites of Fruits and Vegetables with Antioxidant Potential. Second. Metab. Trends Rev. 2022, 141. [Google Scholar] [CrossRef]
  5. Fukaya, M.; Nakamura, S.; Hayashida, H.; Noguchi, D.; Nakashima, S.; Yoneda, T.; Matsuda, H. Structures of Cyclic Organosulfur Compounds From Garlic (Allium sativum L.) Leaves. Front. Chem. 2020, 8, 282. [Google Scholar] [CrossRef] [PubMed]
  6. Prasad, K.; Kaushik, V.; Prasad, K.K. Garlic: A potential source of herbal medicine. Beverage Food World 2002, 29, 24–26. [Google Scholar]
  7. Prakash, P.; Prasad, K. Comparative elucidation of garlic peeling methods and positioning of quality characteristics using principal component analysis. Acta Sci. Pol. Technol. Aliment. 2023, 22, 119–131, in press. [Google Scholar]
  8. JĘDrszczyk, E.; KopeĆ, A.; Bucki, P.; Ambroszczyk, A.M.; Skowera, B. The Enhancing Effect of Plants Growth Biostimulants in Garlic Cultivation on the Chemical Composition and Level of Bioactive Compounds in the Garlic Leaves, Stems and Bulbs. Not. Bot. Horti Agrobot. Cluj-Napoca 2018, 47, 81–91. [Google Scholar] [CrossRef]
  9. Capasso, A. Antioxidant Action and Therapeutic Efficacy of Allium sativum L. Molecules 2013, 18, 690–700. [Google Scholar] [CrossRef] [PubMed]
  10. Lobo, V.; Patil, A.; Phatak, A.; Chandra, N. Free radicals, antioxidants and functional foods: Impact on human health. Pharmacogn. Rev. 2010, 4, 118. [Google Scholar] [CrossRef]
  11. Liu, J.; Ji, F.; Chen, F.; Guo, W.; Yang, M.; Huang, S.; Zhang, F.; Liu, Y. Determination of garlic phenolic compounds using supercritical fluid extraction coupled to supercritical fluid chromatography/tandem mass spectrometry. J. Pharm. Biomed. Anal. 2018, 159, 513–523. [Google Scholar] [CrossRef] [PubMed]
  12. Nagella, P.; Thiruvengadam, M.; Ahmad, A.; Yoon, J.-Y.; Chung, I.-M. Composition of Polyphenols and Antioxidant Activity of Garlic Bulbs Collected from Different Locations of Korea. Asian J. Chem. 2014, 26, 897–902. [Google Scholar] [CrossRef]
  13. Fernández-López, J.; Botella-Martínez, C.; Navarro-Rodríguez de Vera, C.; Sayas-Barberá, M.E.; Viuda-Martos, M.; Sánchez-Zapata, E.; Pérez-Álvarez, J.A. Vegetable Soups and Creams: Raw Materials, Processing, Health Benefits, and Innovation Trends. Plants 2020, 9, 1769. [Google Scholar] [CrossRef] [PubMed]
  14. Barbhuiya, R.I.; Singha, P.; Asaithambi, N.; Singh, S.K. Ultrasound-assisted rapid biological synthesis and characterization of silver nanoparticles using pomelo peel waste. Food Chem. 2022, 385, 132602. [Google Scholar] [CrossRef]
  15. El-Seedi, H.R.; El-Shabasy, R.M.; Khalifa, S.A.; Saeed, A.; Shah, A.; Shah, R.; Iftikhar, F.J.; Abdel-Daim, M.M.; Omri, A.; Hajrahand, N.H. Metal nanoparticles fabricated by green chemistry using natural extracts: Biosynthesis, mechanisms, and applications. RSC Adv. 2019, 9, 24539–24559. [Google Scholar] [CrossRef] [PubMed]
  16. Jha, A.K.; Prasad, K.; Kumar, V.; Prasad, K. Biosynthesis of silver nanoparticles using eclipta leaf. Biotechnol. Prog. 2009, 25, 1476–1479. [Google Scholar] [CrossRef] [PubMed]
  17. Jha, A.K.; Prasad, K.; Prasad, K. A green low-cost biosynthesis of Sb2O3 nanoparticles. Biochem. Eng. J. 2009, 43, 303–306. [Google Scholar] [CrossRef]
  18. Jha, A.K.; Prasad, K.; Prasad, K. Biosynthesis of Sb2O3 nanoparticles: A low-cost green approach. Biotechnol. J. 2009, 4, 1582–1585. [Google Scholar] [CrossRef]
  19. Jha, A.K.; Prasad, K.; Prasad, K.; Kulkarni, A.R. Plant system: Nature’s nanofactory. Colloids Surf. B Biointerfaces 2009, 73, 219–223. [Google Scholar] [CrossRef] [PubMed]
  20. Marslin, G.; Siram, K.; Maqbool, Q.; Selvakesavan, R.K.; Kruszka, D.; Kachlicki, P.; Franklin, G. Secondary metabolites in the green synthesis of metallic nanoparticles. Materials 2018, 11, 940. [Google Scholar] [CrossRef]
  21. Ramić, M.; Vidović, S.; Zeković, Z.; Vladić, J.; Cvejin, A.; Pavlić, B. Modeling and optimization of ultrasound-assisted extraction of polyphenolic compounds from Aronia melanocarpa by-products from filter-tea factory. Ultrason. Sonochem. 2015, 23, 360–368. [Google Scholar] [CrossRef] [PubMed]
  22. Barbhuiya, R.I.; Singha, P.; Singh, S.K. A comprehensive review on impact of non-thermal processing on the structural changes of food components. Food Res. Int. 2021, 149, 110647. [Google Scholar] [CrossRef] [PubMed]
  23. Asaithambi, N.; Singh, S.K.; Singha, P. Current status of non-thermal processing of probiotic foods: A review. J. Food Eng. 2021, 303, 110567. [Google Scholar] [CrossRef]
  24. Dash, D.R.; Singh, S.K.; Singha, P. Recent advances on the impact of novel non-thermal technologies on structure and functionality of plant proteins: A comprehensive review. Crit. Rev. Food Sci. Nutr. 2022, 1–16. [Google Scholar] [CrossRef]
  25. Prasad, K.; Nath, N.; Prasad, K. Estimation of sugar content in commercially available beverages using ultrasonic velocity measurement. Indian J. Phys. 2000, 74A, 387–389. [Google Scholar]
  26. Prasad, K. Advances in Nondestructive Quality Measurement of Fruits and Vegetables. In Postharvest Biology and Technology of Horticultural Crops: Principles and Practices for Quality Maintenance; Apple Academic Press: New York, NY, USA, 2015; pp. 51–87. [Google Scholar]
  27. Prasad, K. Non-destructive Quality Analysis of Fruits. In Postharvest Quality Assurance of Fruits: Practical Approaches for Developing Countries; Springer International Publishing: Cham, Switzerland, 2015; pp. 239–258. [Google Scholar]
  28. Rostagno, M.A.; Palma, M.; Barroso, C.G. Ultrasound-assisted extraction of soy isoflavones. J. Chromatogr. A 2003, 1012, 119–128. [Google Scholar] [CrossRef] [PubMed]
  29. Tomšik, A.; Pavlić, B.; Vladić, J.; Ramić, M.; Brindza, J.; Vidović, S. Optimization of ultrasound-assisted extraction of bioactive compounds from wild garlic (Allium ursinum L.). Ultrason. Sonochem. 2016, 29, 502–511. [Google Scholar] [CrossRef]
  30. Kaur, R.; Shekhar, S.; Chaudhary, S.; Singh, B.; Prasad, K. Non-thermal Food Preservation Technologies. In Smart and Sustainable Food Technologies; Springer Nature: Basel, Switzerland, 2022; p. 157. [Google Scholar]
  31. Patra, A.; Abdullah, S.; Pradhan, R.C. Optimization of ultrasound-assisted extraction of ascorbic acid, protein and total antioxidants from cashew apple bagasse using artificial neural network-genetic algorithm and response surface methodology. J. Food Process. Preserv. 2022, 46, e16317. [Google Scholar] [CrossRef]
  32. Nelofer, R.; Ramanan, R.N.; Rahman, R.N.Z.R.A.; Basri, M.; Ariff, A.B. Comparison of the estimation capabilities of response surface methodology and artificial neural network for the optimization of recombinant lipase production by E. coli BL21. J. Ind. Microbiol. Biotechnol. 2012, 39, 243–254. [Google Scholar] [CrossRef] [PubMed]
  33. Baş, D.; Boyacı, İ.H. Modeling and optimization I: Usability of response surface methodology. J. Food Eng. 2007, 78, 836–845. [Google Scholar] [CrossRef]
  34. Kukreja, T.R.; Kumar, D.; Prasad, K.; Chauhan, R.C.; Choe, S.; Kundu, P.P. Optimisation of physical and mechanical properties of rubber compounds by response surface methodology––Two component modelling using vegetable oil and carbon black. Eur. Polym. J. 2002, 38, 1417–1422. [Google Scholar] [CrossRef]
  35. Kumar, S.; Baniwal, P.; Nayik, G.A.; Prasad, K.; Khan, K.A.; Ghramh, H.A.; Kumar, H.; Karabagias, I.K. Optimization and Development of Ready to Eat Chocolate Coated Roasted Flaked Rice as Instant Breakfast Food. Foods 2021, 10, 1658. [Google Scholar] [CrossRef]
  36. Bezerra, M.A.; Santelli, R.E.; Oliveira, E.P.; Villar, L.S.; Escaleira, L.A. Response surface methodology (RSM) as a tool for optimization in analytical chemistry. Talanta 2008, 76, 965–977. [Google Scholar] [CrossRef]
  37. Ciric, A.; Krajnc, B.; Heath, D.; Ogrinc, N. Response surface methodology and artificial neural network approach for the optimization of ultrasound-assisted extraction of polyphenols from garlic. Food Chem. Toxicol. 2020, 135, 110976. [Google Scholar] [CrossRef]
  38. Desai, K.M.; Survase, S.A.; Saudagar, P.S.; Lele, S.S.; Singhal, R.S. Comparison of artificial neural network (ANN) and response surface methodology (RSM) in fermentation media optimization: Case study of fermentative production of scleroglucan. Biochem. Eng. J. 2008, 41, 266–273. [Google Scholar] [CrossRef]
  39. Chen, S.; Zhang, H.; Yang, L.; Zhang, S.; Jiang, H. Optimization of Ultrasonic-Assisted Extraction Conditions for Bioactive Components and Antioxidant Activity of Poria cocos (Schw.) Wolf by an RSM-ANN-GA Hybrid Approach. Foods 2023, 12, 619. [Google Scholar] [CrossRef]
  40. Simoncini, D.; Zhang, K.Y.J. Population-Based Sampling and Fragment-Based De Novo Protein Structure Prediction. In Encyclopedia of Bioinformatics and Computational Biology; Elsevier: Amsterdam, The Netherlands, 2019; pp. 774–784. [Google Scholar]
  41. Albadr, M.A.; Tiun, S.; Ayob, M.; Al-Dhief, F. Genetic Algorithm Based on Natural Selection Theory for Optimization Problems. Symmetry 2020, 12, 1758. [Google Scholar] [CrossRef]
  42. Ding, S.; Xu, L.; Su, C.; Zhu, H. Using Genetic Algorithms to Optimize Artificial Neural Networks. J. Converg. Inf. Technol. 2010, 5, 54–62. [Google Scholar] [CrossRef]
  43. Kikuchi, T.; Uchida, T. Calorimetric method for measuring high ultrasonic power using water as a heating material. J. Phys. Conf. Ser. 2011, 279, 012012. [Google Scholar] [CrossRef]
  44. AOAC. Official Methods of Analysis. In Standard Methods for the Examination; Association of Official Analytical Chemists: Washington, DC, USA, 2006. [Google Scholar]
  45. Slinkard, K.; Singleton, V.L. Total phenol analysis: Automation and comparison with manual methods. Am. J. Enol. Vitic. 1977, 28, 49–55. [Google Scholar] [CrossRef]
  46. Marinova, D.; Ribarova, F.; Atanassova, M. Total phenolics and flavonoids in Bulgarian fruits and vegetables. J. Univ. Chem. Technol. Metall. 2005, 40, 255–260. [Google Scholar]
  47. Gorinstein, S.; Vargas, O.J.M.; Jaramillo, N.O.; Salas, I.A.; Ayala, A.L.M.; Arancibia-Avila, P.; Toledo, F.; Katrich, E.; Trakhtenberg, S. The total polyphenols and the antioxidant potentials of some selected cereals and pseudocereals. Eur. Food Res. Technol. 2007, 225, 321–328. [Google Scholar] [CrossRef]
  48. Montgomery, D.C. Design and Analysis of Experiments; John Wiley & Sons: New York, NY, USA, 2017. [Google Scholar]
  49. Li, M.; Huang, X.; Liu, H.; Liu, B.; Wu, Y.; Xiong, A.; Dong, T. Prediction of gas solubility in polymers by back propagation artificial neural network based on self-adaptive particle swarm optimization algorithm and chaos theory. Fluid Phase Equilibria 2013, 356, 11–17. [Google Scholar] [CrossRef]
  50. Gupta, S.; Gowri, B.S.; Lakshmi, A.J.; Prakash, J. Retention of nutrients in green leafy vegetables on dehydration. J. Food Sci. Technol. 2013, 50, 918–925. [Google Scholar] [CrossRef] [PubMed]
  51. Arasaretnam, S.; Kiruthika, A.; Mahendran, T. Nutritional and mineral composition of selected green leafy vegetables. Ceylon J. Sci. 2018, 47, 35. [Google Scholar] [CrossRef]
  52. Samaram, S.; Mirhosseini, H.; Tan, C.P.; Ghazali, H.M.; Bordbar, S.; Serjouie, A. Optimisation of ultrasound-assisted extraction of oil from papaya seed by response surface methodology: Oil recovery, radical scavenging antioxidant activity, and oxidation stability. Food Chem. 2015, 172, 7–17. [Google Scholar] [CrossRef]
  53. Yahya, N.A.; Wahab, R.A.; Xine, T.L.S.; Hamid, M.A. Ultrasound-assisted extraction of polyphenols from pineapple skin. In Proceedings of the 2nd International Conference on Biosciences And Medical Engineering (Icbme2019): Towards Innovative Research and Cross-Disciplinary Collaborations, Bali, Indonesia, 11–12 April 2019. [Google Scholar]
  54. Shirzad, H.; Niknam, V.; Taheri, M.; Ebrahimzadeh, H. Ultrasound-assisted extraction process of phenolic antioxidants from Olive leaves: A nutraceutical study using RSM and LC–ESI–DAD–MS. J. Food Sci. Technol. 2017, 54, 2361–2371. [Google Scholar] [CrossRef] [PubMed]
  55. Gadioli Tarone, A.; Keven Silva, E.; Dias de Freitas Queiroz Barros, H.; Baú Betim Cazarin, C.; Roberto Marostica Junior, M. High-intensity ultrasound-assisted recovery of anthocyanins from jabuticaba by-products using green solvents: Effects of ultrasound intensity and solvent composition on the extraction of phenolic compounds. Food Res. Int. 2021, 140, 110048. [Google Scholar] [CrossRef]
  56. Zimare, S.B.; Mankar, G.D.; Barmukh, R.B. Optimization of ultrasound-assisted extraction of total phenolics and flavonoids from the leaves of Lobelia nicotianifolia and their radical scavenging potential. Curr. Res. Green Sustain. Chem. 2021, 4, 100109. [Google Scholar] [CrossRef]
  57. Rios-Romero, E.A.; Ochoa-Martínez, L.A.; Bello-Pérez, L.A.; Morales-Castro, J.; Quintero-Ramos, A.; Gallegos-Infante, J.A. Effect of ultrasound and steam treatments on bioaccessibility of β-carotene and physicochemical parameters in orange-fleshed sweet potato juice. Heliyon 2021, 7, e06632. [Google Scholar] [CrossRef]
  58. Shen, Y.; Zhu, D.; Xi, P.; Cai, T.; Cao, X.; Liu, H.; Li, J. Effects of temperature-controlled ultrasound treatment on sensory properties, physical characteristics and antioxidant activity of cloudy apple juice. Lwt 2021, 142, 111030. [Google Scholar] [CrossRef]
  59. Um, M.; Han, T.-H.; Lee, J.-W. Ultrasound-assisted extraction and antioxidant activity of phenolic and flavonoid compounds and ascorbic acid from rugosa rose (Rosa rugosa Thunb.) fruit. Food Sci. Biotechnol. 2017, 27, 375–382. [Google Scholar] [CrossRef]
  60. Kumar, K.; Srivastav, S.; Sharanagat, V.S. Ultrasound assisted extraction (UAE) of bioactive compounds from fruit and vegetable processing by-products: A review. Ultrason. Sonochem. 2021, 70, 105325. [Google Scholar] [CrossRef] [PubMed]
  61. AL-Bukhaiti, W.Q.; Noman, A.; Mahdi, A.A.; Ali, A.H.; Abed, S.M.; Rashed, M.M.; Wang, H. Profiling of phenolic compounds and antioxidant activities of Cissus rotundifolia (Forssk.) as influenced by ultrasonic-assisted extraction conditions. J. Food Meas. Charact. 2019, 13, 634–647. [Google Scholar] [CrossRef]
  62. Nandiyanto, A.B.D.; Oktiani, R.; Ragadhita, R. How to read and interpret FTIR spectroscope of organic material. Indones. J. Sci. Technol. 2019, 4, 97–118. [Google Scholar] [CrossRef]
  63. Hari, N.; Nair, V.P. FTIR spectroscopic analysis of leaf extract in hexane in Jasminum azoricum L. Recent Res. Sci. Technol. 2018, 4, 170–172. [Google Scholar]
  64. Alotibi, F.; Rizwana, H. Chemical composition, FTIR studies, morphological alterations, and antifungal activity of leaf extracts of Artemisia sieberi from Saudi Arabia. Int. J. Agric. Biol 2019, 21, 1241–1248. [Google Scholar]
  65. Ghitescu, R.-E.; Volf, I.; Carausu, C.; Bühlmann, A.-M.; Gilca, I.A.; Popa, V.I. Optimization of ultrasound-assisted extraction of polyphenols from spruce wood bark. Ultrason. Sonochem. 2015, 22, 535–541. [Google Scholar] [CrossRef]
  66. Medina-Jaramillo, C.; Gomez-Delgado, E.; López-Córdoba, A. Improvement of the Ultrasound-Assisted Extraction of Polyphenols from Welsh Onion (Allium fistulosum) Leaves Using Response Surface Methodology. Foods 2022, 11, 2425. [Google Scholar] [CrossRef]
  67. Itakura, Y.; Ichikawa, M.; Mori, Y.; Okino, R.; Udayama, M.; Morita, T. How to Distinguish Garlic from the Other Allium Vegetables. J. Nutr. 2001, 131, 963S–967S. [Google Scholar] [CrossRef]
  68. Szychowski, K.; Rybczyńska-Tkaczyk, K.; Gaweł-Bęben, K.; Świeca, M.; Karaś, M.; Jakubczyk, A.; Matysiak, M.; Binduga, U.; Gmiński, J. Characterization of Active Compounds of Different Garlic (Allium sativum L.) Cultivars. Pol. J. Food Nutr. Sci. 2018, 68, 73–81. [Google Scholar] [CrossRef]
Figure 1. Process flow diagram for garlic leaf powder preparation.
Figure 1. Process flow diagram for garlic leaf powder preparation.
Foods 12 01925 g001
Figure 2. The effect of ultrasound amplitude, extraction time, and ethanol concentration on (ac) extraction yield, (df) total phenol content, (gi) total flavonoid content, and (jl) antioxidant activity of garlic leaf powder.
Figure 2. The effect of ultrasound amplitude, extraction time, and ethanol concentration on (ac) extraction yield, (df) total phenol content, (gi) total flavonoid content, and (jl) antioxidant activity of garlic leaf powder.
Foods 12 01925 g002
Figure 3. Optimal neural network architecture for present experiment.
Figure 3. Optimal neural network architecture for present experiment.
Foods 12 01925 g003
Figure 4. Post-training performance and error histogram of (a,b) extraction yield, (c,d) total phenol content, (e,f) total flavonoid content, and (g,h) antioxidant activity of generated artificial neural network model.
Figure 4. Post-training performance and error histogram of (a,b) extraction yield, (c,d) total phenol content, (e,f) total flavonoid content, and (g,h) antioxidant activity of generated artificial neural network model.
Foods 12 01925 g004aFoods 12 01925 g004b
Figure 5. Fitness value variation across generations during optimization with a genetic algorithm.
Figure 5. Fitness value variation across generations during optimization with a genetic algorithm.
Foods 12 01925 g005
Figure 6. FE-SEM images of garlic leaf powder (a) before treatment and (b) after treatment by ultrasonic processing at optimized ANN-GA condition (60% ultrasound amplitude, 13 min treatment time, and 53% ethanol concentration).
Figure 6. FE-SEM images of garlic leaf powder (a) before treatment and (b) after treatment by ultrasonic processing at optimized ANN-GA condition (60% ultrasound amplitude, 13 min treatment time, and 53% ethanol concentration).
Foods 12 01925 g006
Figure 7. FTIR analysis of the UAE extract at optimized ANN−GA condition (60% ultrasound amplitude, 13 min treatment time, and 53% ethanol concentration).
Figure 7. FTIR analysis of the UAE extract at optimized ANN−GA condition (60% ultrasound amplitude, 13 min treatment time, and 53% ethanol concentration).
Foods 12 01925 g007
Figure 8. HPLC chromatogram of identified polyphenols of garlic leaf powder.
Figure 8. HPLC chromatogram of identified polyphenols of garlic leaf powder.
Foods 12 01925 g008
Table 1. The experimental domain of rotatable central composite design.
Table 1. The experimental domain of rotatable central composite design.
Coded Level
Independent VariableUnits−α*−101α*
Ultrasound Amplitude (X1)%19.77 (~20)30456070.23 (~70)
Treatment Time (X2)min1.60 (~2)5101518.40 (~18)
Ethanol Conc. (X3)%33.18 (~33)40506066.82 (~67)
* The (±) α values of the independent parameters are represented as round figures.
Table 2. Rotatable central composite design of three factors with their experimental and predicted responses using RSM and ANN.
Table 2. Rotatable central composite design of three factors with their experimental and predicted responses using RSM and ANN.
RunSpace TypeUltrasound Amplitude (X1)Treatment Time (X2)Ethanol Conc. (X3)Applied Energy (J) *Calorimetric Energy (J) **Yield (%)TPC (mg GAE/g)TFC (mg QE/g)Antioxidant (%)
ExperimentalRSM Pred.ANN Pred.ExperimentalRSM Pred.ANN Pred.ExperimentalRSM Pred.ANN Pred.ExperimentalRSM Pred.ANN Pred.
1Center45105016,744294.3431.6531.9131.909.819.859.856.806.716.6957.0457.8457.85
2Center45105016,744277.7332.1431.9131.909.899.859.856.786.716.6957.7457.8457.85
3Factorial6054011,187302.9830.6530.5130.659.499.489.496.326.096.3255.2354.7055.23
4Center45105016,744837.1931.7431.9131.909.829.859.856.856.716.6957.8757.8457.85
5Factorial30156018,36697.6831.1631.2931.169.639.609.635.976.145.9756.5556.7156.55
6Factorial305606261188.3831.2631.1131.269.579.569.575.725.705.7256.9356.3556.93
7Factorial305407390173.1628.7128.8428.719.389.319.385.395.215.3953.2152.3453.21
8Axial452502846166.5030.4530.4230.459.479.539.475.495.785.4953.8654.8953.86
9Axial70105023,389336.3831.7331.6831.739.829.829.826.887.016.8857.8057.5557.80
10Axial20105010,390134.6630.1130.1730.119.389.429.385.535.495.5353.3954.1653.39
11Factorial30154022,332219.7831.3131.1031.319.629.639.625.685.725.6854.3954.0754.39
12Axial45106713,877158.8531.6231.5931.629.759.769.756.696.576.6958.5858.3858.58
13Axial45185033,414116.8032.1232.1632.129.899.889.896.716.516.7157.6257.1257.62
14Center45105016,744170.0932.2131.9131.909.909.859.856.676.716.6958.5957.8457.85
15Center45105016,744173.1632.0531.9131.909.889.859.856.546.716.6957.7357.8457.85
16Factorial60154030,957109.4132.2432.3932.249.889.869.886.576.526.5756.7956.9956.79
17Axial45103317,820166.5030.4630.5030.469.489.529.485.485.705.4853.1253.8453.12
18Factorial6056010,050173.1631.4231.6231.429.849.809.846.826.716.8257.5257.4757.52
19Factorial60156028,857146.5231.5731.4331.579.889.919.886.967.076.9657.8958.3957.89
20Center45105016,744173.1631.6731.9131.909.809.859.856.636.716.6958.1357.8457.85
Abbreviations: EY, extraction yield; TPC, total phenol content; TFC, total flavonoid content; Antioxidant, antioxidant activity; RSM pred., response surface methodology predicted value; ANN, artificial neural network predicted value; * applied energy obtained from instrument; ** calorimetric energy calculated as per Kikuchi and Uchida [43].
Table 3. Proximate composition of raw garlic leaves.
Table 3. Proximate composition of raw garlic leaves.
Chemical ComponentsAmount (%)
Moisture90.04
Protein2.23
Fat0.18
Ash1.56
Carbohydrate5.99
Table 4. Estimated coefficients of the fitted model and analysis of variance table representing the regression coefficients, R2, adj. R2, and lack of fit of the responses.
Table 4. Estimated coefficients of the fitted model and analysis of variance table representing the regression coefficients, R2, adj. R2, and lack of fit of the responses.
CoefficientYield (%)TPC (mg GAE/g)TFC (mg QE/g)Antioxidant (%)
b0+31.91+9.85+6.71+57.84
b1+0.4514 *+0.1194 *+0.4526 *+1.01 *
b2+0.5161 *+0.1052 *+0.2183 *+0.6629 *
b3+0.3259 *+0.0735 *+0.2596 *+1.35 *
b12−0.0950+0.0163−0.0188+0.1412
b13−0.2875 *+0.0187+0.0337−0.3113
b23−0.5175 *−0.0662 *−0.0188−0.3438
b11−0.3482 *−0.0794 *−0.1614 *−0.6993 *
b22−0.2192 *−0.0511 *−0.1985 *−0.6481 *
b33−0.3058 *−0.0740 *−0.2038 *−0.6092 *
R20.9630.9620.9320.917
Adj. R20.9300.9280.8700.842
Lack of fit0.6600.3360.0510.109
* Significant at p ≤ 0.05. Abbreviations: Yield, extraction yield; TPC, total phenol content; TFC, total flavonoid content. b denotes the regression model; b0 denotes the constant term; b1, b2, and b3 denote the linear terms; b12, b13, and b23 denote the interaction terms; and b11, b22, and b33 denote the quadratic terms. 1, 2, and 3 signify independent parameters, such as ultrasound amplitude, treatment time, and ethanol concentration.
Table 5. List of all statistical parameters (R2, NRMSE, RSME, MPE, AAD, and MSE) of responses.
Table 5. List of all statistical parameters (R2, NRMSE, RSME, MPE, AAD, and MSE) of responses.
EYTPCTFCAntioxidant
RSMANNRSMANNRSMANNRSMANN
RMSE0.16300.12590.03530.02240.14420.06040.52980.2558
NMSE0.00520.00400.00360.00230.02280.00950.00940.0045
MSE0.02660.01580.00120.00050.02080.00360.28070.0655
NRMSE0.00080.00050.00010.00010.00330.00060.00500.0012
MPE0.45000.20990.31600.12182.01210.43880.78600.1799
AAD0.14160.06700.03060.01200.12480.02950.43760.1040
R20.96330.97810.96220.98480.93170.98830.91710.9807
Abbreviations: EY, extraction yield; TPC, total phenol content; TFC, total flavonoid content, Antioxidant, antioxidant activity; RSM, response surface methodology; ANN, artificial neural network; R2, coefficient of multiple determination; NRMSE, normal root mean square; RSME, root mean square error; MPE, mean percentage error; NMSE, normal mean square error; MSE, mean square error; AAD, average absolute deviation.
Table 6. Concentration of organosulfur and phenolic compounds in garlic leaf powder using HPLC at 254 nm and 280 nm.
Table 6. Concentration of organosulfur and phenolic compounds in garlic leaf powder using HPLC at 254 nm and 280 nm.
Peak Number    CompoundsConcentration (ppm)
Organosulfur Alliin90.207
S-Allyl-L-cysteine4.314
Allicin219.536
Phenolic1Gallic acid13.591
23,4-Dihydroxybenzoic acid34.403
3Chlorogenic acid36.537
4Catechin hydrate26.327
5Syringic acid4.276
6p-Coumaric acid1.228
7Rutin21.741
8Ellagic acid0.968
9Benzoic acid15.234
10Hesperidin9.272
11Quercetin6.140
12β-Carotene117.607
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shekhar, S.; Prakash, P.; Singha, P.; Prasad, K.; Singh, S.K. Modeling and Optimization of Ultrasound-Assisted Extraction of Bioactive Compounds from Allium sativum Leaves Using Response Surface Methodology and Artificial Neural Network Coupled with Genetic Algorithm. Foods 2023, 12, 1925. https://doi.org/10.3390/foods12091925

AMA Style

Shekhar S, Prakash P, Singha P, Prasad K, Singh SK. Modeling and Optimization of Ultrasound-Assisted Extraction of Bioactive Compounds from Allium sativum Leaves Using Response Surface Methodology and Artificial Neural Network Coupled with Genetic Algorithm. Foods. 2023; 12(9):1925. https://doi.org/10.3390/foods12091925

Chicago/Turabian Style

Shekhar, Shubhra, Prem Prakash, Poonam Singha, Kamlesh Prasad, and Sushil Kumar Singh. 2023. "Modeling and Optimization of Ultrasound-Assisted Extraction of Bioactive Compounds from Allium sativum Leaves Using Response Surface Methodology and Artificial Neural Network Coupled with Genetic Algorithm" Foods 12, no. 9: 1925. https://doi.org/10.3390/foods12091925

APA Style

Shekhar, S., Prakash, P., Singha, P., Prasad, K., & Singh, S. K. (2023). Modeling and Optimization of Ultrasound-Assisted Extraction of Bioactive Compounds from Allium sativum Leaves Using Response Surface Methodology and Artificial Neural Network Coupled with Genetic Algorithm. Foods, 12(9), 1925. https://doi.org/10.3390/foods12091925

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop