Simultaneous Optimization of Extraction Yield, Phenolic Compounds and Antioxidant Activity of Moroccan Propolis Extracts: Improvement of Ultrasound-Assisted Technique Using Response Surface Methodology

: Propolis has given rise to refreshing interest in recent years in the ﬁeld of conventional medicine. Its extraction represents an important process that requires optimal conditions, which strongly affect the yield of extraction, total polyphenols, ﬂavonoid content, and radical scavenging capacity markers. The objective of the present study was to optimize the ultrasound-assisted extraction conditions of Moroccan propolis. The studied responses were the extraction yield, total polyphenols, ﬂavonoid contents (TPC, TFC), and antioxidant activity of the extract evaluated by DPPH-IC 50 and FRAP-EC 50 assays. The response surface methodology (RSM) and speciﬁcally the Box–Behnken design (BBD) were used, taking into account three variables: sonication time (min), solvent/propolis ratio (mL/g), and ethanol concentration (%). After the realization of experiments and data analysis, optimal response values were 15.39%, 192 mg GAE/g of propolis,45.15 mg QEq/g, 29.8 µ g/mL, and 128.3 µ mol Fe 2+ /g for extraction yield, TPC, TFC, DPPH-IC 50 , and FRAP-EC 50 , respectively. Besides, optimal ultrasound extraction conditions were 15 min for sonication time, 30 mL/g for solvent/propolis ratio, and 40% for ethanol concentration. All obtained experimental values were in good agreement with the predicted values, suggesting that using an experimental design in the ultrasound-assisted extraction process and optimization was prudently chosen.


Introduction
Propolis is collected by bees on the buds and bark of certain trees, creating specific modifications in the hive bymixing their salivary secretions with wax, pollen, and collected natural resins [1]. Moreover, bees use propolis to maintain their hives by sealing holes in their honeycombs, by smoothing out the internal walls, and by sheltering the entrance of the hive from intruders [2].
The choice of the extraction solvent type, solvent ratio, and extraction proceduresare key steps in obtaining high-quality propolis with good pharmacological values [12]. Thereby, many traditional and modern methods of extraction of propolis are used, such as

Chemicals and Reagents
All chemicals and reagents used in this study were of analytical grade. The optical density was measured on a spectrophotometer (Spectrophotometer 2005 UV-Vis, Selecta, Barcelona, Spain).

Propolis Preparation
Moroccan propolis samples were directly collected by beekeepers from Apis mellifera hives from the Fez-Meknes region located in the Sidi Hrazem area (33 • 57 49.134 N; 4 • 47 29.758 W), after honey extraction. The collection of propolis was started by scratching hive walls and frames, followed by the removal of debris of wood and bees. The obtained propolis was protected from light and immediately transferred to the laboratory where it was stored at −20 • C for 10 days before being comminuted. Wax was removed by successive multiple extractions with cyclohexane and dichloromethane according to Boissard et al. [20].

Ultrasound-Assisted Extraction (UAE) Conditions and Yield of Extraction
The extractions were carried out in a ultrasonic water bath using a direct sonication method with some modifications [21].The frequency of sonication was 50 kHz, the power was 120 W and the temperature was fixed at 35 • C. Five grams of Moroccan propolis were treated according to the conditions proposed by each experiment of BBD. Then, each mixture was filtered through Whatman paper and dried in a rotavapor at 40 • C to remove the solvent. Finally, the extracts were kept in dark tubes and stored at −20 • C for later use.
The yield (%) of the extraction was evaluated by comparing the dry weight of the extract with the initial weight of propolis using the formula described by Oroian and al [22]: where Wpe is the weight of propolis extract (g) and Wpm is the weight of the initial dried sample.

Total Phenolic Content Determination
Total phenolic concentration in the studied propolis sample was determined spectrophotometrically according to the Folin-Ciocalteu colorimetric method. Then, 50 µL of Processes 2022, 10, 297 3 of 16 each extract solution was mixed with 0.2 mL of the Folin-Ciocalteu reagent and 0.8 mL of 15.9% sodium carbonate. The final volume was made to reach 3 mL with distilled water. The sample was left for 20 min at room temperature, and the absorbance was measured at 760 nm [23]. Calibration was performed using gallic acid as a reference compound.
The results were expressed as mg of gallic acid equivalent per g of crude propolis (mg GAE/g). Triplicate measurements were performed for each extract.

Total Flavonoids Content Determination
The total flavonoid content (TFC) of the propolis sample was determined using aluminum chloride, and the results were expressed as mg quercetin equivalent per g of crude propolis (mg QE/g). Two milliliters (2 mL) of the stock solution were mixed with 3 mL of a 5% aluminum chloride solution. After 30 min of incubation, the absorbance of the reaction mixture was measured at 437 nm against a methanol blank [24]. TFC determination was carried out in triplicate, and a standard curve of quercetin was established.

Diphenyl-1-Picrylhydrazyl (DPPH) Assay
The radical-scavenging activity of selected propolis was determined spectrophotometrically with minor modifications. Crude Moroccan propolis extract was prepared as 1 mg/mL in methanol. Fifty microliters of different concentration samples were added to 5 mL of 0.004% methanol solution of DPPH reagent. After 30 min of incubation in the dark at room temperature, the change in absorbance was measured at 517 nm [25].
The assay was carried out in triplicate, and the percentage of inhibition was calculated using the following formula: where A0 and A1 are the absorbances at 30 min of the control and the sample, respectively.

Ferric-Reducing Antioxidant Power Assay
FRAP working solution was prepared freshly each time: 0.3 M acetate buffer (pH = 3.6), 0.01 M TPTZ (2,4,6-tripyridyl-s-triazine) in 0.04 M of HCl and 0.02 M of FeCl 3 6H 2 O was mixed in 10:1:1 (v/v/v) and kept away from light, and then 0.075 mL of EEP was added to 2.25 mL of FRAP working solution and 0.225 mL of deionized water; the mixture was vortexed and incubated at 37 • C for 30 min. A calibration curve should be prepared with ferrous sulfate (200, 400, 600, 800, and 1000 µM).
The absorbance was measured at 593 nm. FRAP working solution with deionized water instead of a sample was used as a blank. Results were carried out in triplicate and expressed as µmol Fe 2+ /g [26].

HPLC-DAD Analysis
Individual phenolic separation and detection were analyzed using PROMINENCE HPLC: High-Performance Liquid Chromatography system (Shimadzu, Kyoto, Japan) equipped with a DAD detector, a Rheodyne™ Sampler with a loop of 20 µL, a degasser and a quaternary pump. The separation of the polyphenols was carried out on an Agilent Zorbax C18 column with dimension 4.6 mm × 250 mm × 5 um 100 Å, with a flow rate of 1 mL/min of a mobile phase at a ternary gradient of acetonitrile methanol and acidify water. The column temperature was 30 • C, and the sample injection volume was 20 µL. Under the same conditions, syringic acid and tyrosol standards were injected to determine the response factor. A solvent system consisting of 0.2% orthophosphoric acid in water (solvent A) and methanol (solvent B) and acetonitrile (solvent C) was used with the following gradient The gradient was as follows: 0 min-96% A/2% B/2% C; 40 min-50% A/25%B/25% C; 45 min-40% A/30% B/30% C; 60 min-50% B/50% C; 72 min-96% A/2% B/2%. Phenolic compounds were identified based on the retention times of standard materials, and quantification was achieved by the absorbance recorded in the chromatograms relative to external standards, at 280 nm catechin, epicatechin, p-coumaric acid, caffeic acid, pinoresinol, cinnamic acid, apigenin, ferulic acid, and kaempferol. All standards calibration curves showed high degrees of linearity (R 2 > 0.99) [27].

Experimental Design for Optimization
Response surface methodology (RSM) is a multivariate statistical technique that involves complex calculations for the optimization process [28,29]. This approach develops a suitable experimental design that combines all of the independent variables and uses the data input from the experiment to finally produce a set of equations that can give a theoretical value of an output. The outputs are obtained from a well-designed regression analysis that is based on the controlled values of independent variables [30].
In the current study, a three-factors, three-level Box-Behnken design (BBD) was employed with a total of 15 experiments to evaluate the effect of the independent variables, including sonication time(min)/X 1 , solvent/material ratio (mL/g)/X 2 , and ethanol concentration (% v/v)/X 3 on yield of extraction, TPC, TFC and free radical scavenging using DPPH-IC 50 and FRAP-EC 50 tests. The ultrasound-assisted extractions were conducted under experimental conditions, as represented in Table 1.

Mathematical Model
To predict the optimum ultrasound-assisted extraction conditions for studied responses from Moroccan propolis and to express the interaction between dependent and independent factors, the BBD design assumed that interactions with each other are measured using a second-order polynomial model as below [31]. Y = b 0 + b 1 X 1 + b 2 X 2 + b 3 X 3 + b 12 X 1 X 2 + b 13 X 1 X 3 + b 23 X 2 X 3 + b 11 where Y represents the predicted response value; X 1 , X 2 , X 3 are the independent variables; b 0 represents the theoretical mean value of the response when all factors are in the level 0; b 1 , b 2 , b 3 are linear regression coefficients; b 11 , b 22 , b 33 are quadratic regression coefficients; b 12 , b 13 , b 23 are interaction regression coefficients; and ε: is the regression error term.

Statistical Analysis
JMP software V14 and Design EXPERT V.12 were used to analyze the experimental data through the investigation of the influence of the independent variables on the response factors [32]. The software was used to generate response surfaces and contour plots according to the Box-Behnken design, while holding a variable constant in the quadratic polynomial second-order model. The results were expressed as the mean ± standard deviation (SD). The comparison of means was performed by ANOVA and Tuckey's multiple comparison test. The main effect of regression was considered to be statistically significant for p values less than 0.05. In this way, the ratio between the mean square regression (MS R ) and the mean square residual (MS r ) and F-ratio R/r was used to establish if the model was statistically significant. In addition, the ratio between the mean square lack of fit (MS LOF ) and the mean square pure error (MS PE ) and F-ratio LOF/PE was used to evaluate if the model was well adjusted to the observations. High values of F ratio LOF/PE evidenced a lack of fit for the model. The coefficient of determination R 2 was used to judge the adequacy of the postulated model. Finally, the significance of coefficients was analyzed based on their p value using the t test. The coefficient is considered statistically significant if the p value is less than 0.05.

Optimization Tools
Contours plots were used to depict the effects of the three factors on the dependent variables. The red color on the contour plot indicates that the response is increasing, and blue one indicates that it is decreasing; an elliptical contour plot indicates that the interaction between the variables is significant [33]. Next, the simultaneous optimization of the five responses considered in our study was realized by Derringer's desirability function; this approach determines the most desirable response values by setting the optimum conditions [34]. This tool allows us to give the exact optimal adjustment with a percentage between0 and 1. The value 0 is assigned when the factors lead to an unacceptable response, and a value of 1 represents the maximum desired response [35].

Box-Behnken Design Results
To identify the effects of process variables on the ultrasound-assisted extraction (UAE) of total polyphenol and flavonoid content and DPPH-IC 50 and FRAP-EC 50 from Moroccan propolis, the three-level Box-Behnken design was conducted using three variables (sonication time, solvent-to-material ratio, and concentration ratio). The various combinations of coded experimental conditions with their respective experimental observed data (mean from triplicates) are represented in Table 2. The experiments were carried out after randomization, and every response was the average of three replicates. Before proceeding with the design analysis, the recorded means for the experiments of each of the five responses show a statistically significant difference (p value < 0.05). Furthermore, Tuckey's test shows the means, which can be considered statistically identical. The results of the present study indicate that extraction yield, total polyphenol, and flavonoid content range from 6% to 13.1%, from 86.98 to 187.21 mg GAEq/g, and from 9.7 to 38.8 mg QEq/g of crude propolis, respectively, and antioxidant activity tests DPPH-IC 50 and FRAP-EC 50 range respectively from 21 to 39.15 µg/mL and form 100 to 130.21 µmol Fe 2+ /g. This result confirms the influence of the chosen parameters on all studied responses.

Statistical Validation of the Postulated Model
The obtained experimental data were statistically analyzed by analysis of variance (ANOVA), and the significance of the regression coefficients was evaluated by their corresponding p values. The results in Table 3 indicatethat the models were significantly established due to the extremely low probability value (p < 0.05) for all studied responses where p values were 0.001, 0.04, 0.0156, 0.0052, and 0.0370 for extraction yield, TPC, TFC, DPPH-IC 50 , and FRAP-EC 50 , respectively. The calculated F-Ratio(R/r) for every response was higher than the value of F-ratio (0.05; 9.5) at 95%of confidence level, which is equal to 4.77. In addition, all responses show a lack of fit due to their p values, which were more than 0.05, and all their calculated F-ratios (LOF/PE) were lower than the theoretical F value (0.05; 3.2) at 95%of confidence, which was equal to 19.15. The coefficients of determination R 2 were equal to 98%, 94.3%, 93.7%, and 96% for extraction yield TPC, TFC, DPPH-IC 50, and FRAP-EC 50 , respectively.
These values reflect a good relationship between the experimental and predicted values. These results are confirmed by the graphs in Figure 1, demonstrating a linear curve for the actual values in terms of the predicted ones. These values reflect a good relationship between the experimental and predicted values. These results are confirmed by the graphs in Figure 1, demonstrating a linear curve for the actual values in terms of the predicted ones.  Table 4 gives the linear effects of all studied factors, their quadratic and interaction terms, their Student's t statistical values, and the observed probability (pvalue) of each of them. The relationships between the tested parameters and the responses were explained   Table 4 gives the linear effects of all studied factors, their quadratic and interaction terms, their Student's t statistical values, and the observed probability (p value) of each of them. The relationships between the tested parameters and the responses were explained by the following regression equations (2, 3, 4, 5, and 6 for yield of extraction, TPC, TFC, DPPH-IC 50 , and FRAP-EC 50 , respectively). For the response yield of extraction, the constants b 0 , the linear terms b 2 and b 3 ; the binary interaction terms b 12 , b 13 , b 23 ; and the quadratic interaction b 11 were statistically significant coefficients. As a result, the mathematical model is represented by the following equation:

Model Parameter Estimation and Fitted Models
Concerning the TPC response, statistically significant coefficients were the constant b0; the linear terms b 1 , b 3 ; and the quadratic interaction term b 11 . Its mathematical model is described by the equation: Y TPC = 91.14 − 12.52 X 3 + 50.48 X 2 11 + ε As for the TFC response, the significant coefficients were the constant b 0 ; the linear terms b 1 and b 3 ; and the quadratic interaction terms b 11 and b 33 . The following equation represents the model linking TFC to the parameters: Y TFC = 11.39 − 5.58 X 1 − 8.89 X 3 + 11.61 X 11 + 7.08 X 33 + ε Regarding the DPPH-IC 50 response, the constants b 0 ; the linear terms b 2 , b 3 ; and the quadratic terms b 11 and b 33 were the statistically significant coefficients. Therefore, the fitted model is represented by the following equation: Y DPPH-IC50 = 26.43 + 2.70 X 2 + 1.78 X 3 − 2.38 X 11 + 6.61 X 33 + ε Finally, only the constant b 0 and quadratic interaction terms b 11 and b 22 were the statistically significant coefficients for the FRAP-EC 50 response. The following equation represents its mathematical model: Y FRAP-EC50 = 104.20 + 13.88 X 2 11 + 7.72 X 2 22 + ε The effects of the studied parameters, their estimated values, and their observed p values are significantly influenced by the sonication time (X 1 ), the ratio (X 2 ), and the concentration of solvent (X 3 ), which were compatible with the results found in other studies [21,36]. The ethanol-distilled water mixtures give better extraction efficiency of phenolic compounds compared to single-solvent systems [37]. By increasing the proportion of distilled water to ethanol, the polarity of the solvent increases; thus, the solvent system will be able to extract polyphenols with high polarity, those with low polarity, and those with medium polarity [38]. Several studies have shown the importance of the concentration of the solvent compared to other factors [39][40][41]. Moreover, the extraction time, if combined with the right type of solvent, influences extraction yield, TPC, TFC, and antioxidant activity [42].

Optimization of Ultrasound-Assisted Extraction Parameters
To understand the best value that we can obtain, we began with the best value when performing the experiments. Thus, the best recorded values were 13.1%, 187.21 mg GAEq/g, 42.6 mg QEq/g, 21 µg/mL and 130.21 µmol Fe 2+ /g for extraction yield, TPC, TFC, DPPH-IC 50 and FRAP-EC 50 , respectively. Consequently, a setting of the parameters allowing for the obtainment of responses better than or equal to these values were accepted.
For optimizing studied responses, the temperature used during our experiments was fixed at 35 • C and was suitable for extraction, avoiding the loss or degradation of active compounds [43]. After a first analysis of the five-response optimization, a fixation of the extraction time in its lower level (15 min) is strongly recommended for the three first responses (extraction yield, TPC, and TFC). Therefore, this factor will be fixed in its value of 15 min for these responses. However, maximization of time processing is suggested for DPPH-IC 50 and FRAP-EC 50 .

Effect of Operating Conditions on the Propolis Extraction Yield Response
Based on the contours shown in Figure 2a, the highest value of extraction yield observed can reach 15% by minimizing the processing time and solvent concentration and by maximizing the solvent/material ratio. The desirability plot in Figure 3a confirms this finding by showing that it is possible to obtain an extraction yield of 15.35% with desirability of 99.9% by setting the sonication time at 15 min, propolis-ethanol ratio at 30:1 mg/L, and ethanol concentration at a low level of 40%. The increase in extraction yield can be explained by the cavitation effect caused by ultrasonically disrupting the material and allowing the formation of a porous structure [44]. Our optimized propolis extraction yield value (15.35%) was more important than that obtained by shacking the ethanolic extraction [44], and it was similar for polish propolis ultrasound-assisted extracts, which presented an extraction yield ranging between 10.04% and 15.92% [45].

Effect of Operating Conditions on the TPC Response
As shown in the contour plot (Figure 2b), obtaining a TPC that exceeds 190 mgQAEq/g is possible by the same operating conditions described for the extraction yield, namely, a sonication time of 15 min, a solvent/material ratio of 30:1, and an ethanol concentration of 40%. Furthermore, the desirability function (Figure 3b) indicates that a value of 192 mg GAEq/g of TPC was the highest possible value, with 99.7% as a degree of compromise. This high concentration of phenolic compounds due to cavitation, which is easily released from the propolis matrix [46], and the presence of water in the 40% w/w ethanol increased ethanol polarity, which accelerated solvation of the phenolic compounds. In addition, water also caused swelling in the propolis matrix, which allows ethanol to penetrate the cell structure of solid particles more easily [47]. The result of the optimized value of TPC (192 mg GAE/g) was important in comparison with that obtained by ultrasound-assisted extraction from Malaysian propolis, which the predicted value of 162.46 mg GAEq/g [48], which was more important than Lithuanian propolis extract obtained by traditional extraction methods [49]. Besides, it was similar to those obtained by the same extraction method for Chinese propolis, where the TPC range was from 166 to 201.07 mgQEq/g [50].

Effect of Operating Conditions on TFC Response
The maximum level of total flavonoid content that we could obtain exceeded 45 mg QEq/g (Figure 2c). The desirability plot (Figure 3c) confirms this finding by showing that a setting of 15 min, 30 mL/g, and 40% for sonication time, solvent/propolis ratio, and ethanol concentration, respectively, leads to the best value of TFC (45.15 mg QEq, with 99.3% as the degree of compromise). This value was interesting compared to that obtained by maceration from Malaysian propolis, which reported TFC values in the range from 7.68 to 17.23 mgQE/g [51]. Based on the contours shown in Figure 2a, the highest value of extraction yield observed can reach 15% by minimizing the processing time and solvent concentration and by maximizing the solvent/material ratio. The desirability plot in Figure 3a confirms this finding by showing that it is possible to obtain an extraction yield of 15.35% with desirability of 99.9% by setting the sonication time at 15 min, propolis-ethanol ratio at 30:1 mg/L, and ethanol concentration at a low level of 40%. The increase in extraction yield can be explained by the cavitation effect caused by ultrasonically disrupting the material and allowing the formation of a porous structure [44]. Our optimized propolis extraction yield value (15.35%) was more important than that obtained by shacking the ethanolic extraction [44], and it was similar for polish propolis ultrasound-assisted extracts, which presented an extraction yield ranging between 10.04% and 15.92% [45].

Effect of Operating Conditions on DPPH-IC 50 Response
As illustrated in the contour plot (Figure 2d), the optimized DPPH-IC 50 value can reach around 21 µg/mL; this value can be obtained by the following setting: 45 min, 10 mg/L, and 60% for sonication time, solvent/material ratio, and ethanol concentration, respectively. The desirability plot (Figure 3d) shows that we have a 99.7% of chance to reach the DPPH-IC 50 value of 20.27 µg/mL by ensuring the aforementioned operating conditions. This value was more important in comparison with that obtained by conventional extraction from Brazilian propolis, in which DPPH-IC 50 values ranged between 31 and 183 µg/mL [52].On the contrary, it was shown by Yeo

Effect of Operating Conditions on the TPC Response
As shown in the contour plot (Figure 2b), obtaining a TPC that exceeds 190 mgQAEq/g is possible by the same operating conditions described for the extraction yield, namely, a sonication time of 15 min, a solvent/material ratio of 30:1, and an ethanol concentration of 40%. Furthermore, the desirability function (Figure 3b) indicates that a value of 192 mg GAEq/g of TPC was the highest possible value, with 99.7% as a degree of compromise. This high concentration of phenolic compounds due to cavitation, which is easily released from the propolis matrix [46], and the presence of water in the 40% w/w ethanol increased ethanol polarity, which accelerated solvation of the phenolic compounds. In addition, water also caused swelling in the propolis matrix, which allows ethanol to penetrate the cell structure of solid particles more easily [47]. The result of the optimized value of TPC (192 mg GAE/g) was important in comparison with that obtained by ultrasoundassisted extraction from Malaysian propolis, which the predicted value of 162.46 mg GAEq/g [48], which was more important than Lithuanian propolis extract obtained by traditional extraction methods [49]. Besides, it was similar to those obtained by the same extraction method for Chinese propolis, where the TPC range was from 166 to 201.07 mgQEq/g [50].

Effect of Operating Conditions on TFC Response
The maximum level of total flavonoid content that we could obtain exceeded 45 mg QEq/g (Figure 2c). The desirability plot (Figure 3c) confirms this finding by showing that a setting of 15 min, 30 mL/g, and 40% for sonication time, solvent/propolis ratio, and ethanol concentration, respectively, leads to the best value of TFC (45.15 mg QEq, with 99.3% as the degree of compromise). This value was interesting compared to that obtained by maceration from Malaysian propolis, which reported TFC values in the range from 7.68 to 17.23mgQE/g [51].

Effect of Operating Conditions on FRAP-EC 50 Response
As represented in the contour plot (Figure 2e), the optimized FRAP-EC 50 value was around 135 µmol Fe 2+ /g, obtained under the following sitting: 15 min, 30 mL/g, and 40% for sonication time, solvent/material ratio, and ethanol concentration, respectively. Furthermore, 136.26 µmol Fe 2+ /g was the highest value of FRAP-EC 50 , indicated by the desirability function (Figure 3e) with 99.1% as the degree of compromise. This value was similar to that obtained by ultrasound-assisted extraction from Chinese propolis, in which the values ranged from 126 ± 10.60 to 290.34 ± 10.80 µmol Fe 2+ /g of extract.

Simultaneous Multi-Response Optimization
In addition to its ability to give optimal solutions for each of the responses separately, the desirability function also allows for the determination of optimal conditions for all the studied responses simultaneously [54]. Thus, several solutions can be proposed by this function, implying the determination of the priority given to each response in advance. In our case, the optimization of the first three responses yields, TPC, and TFC require the same setting, but the two responses, DPPH-IC 50 and FRAP-EC 50 , require a different one. Among the proposed solutions of simultaneous settings, we have chosen the conditions illustrated in Figure 4. Hence, a sonication time processing of 15 min, a solvent/material ratio of 30:1 mL/g, and a solvent concentration of 58% were chosen as optimal conditions for simultaneous optimization. This setting allows us to obtain 12.32%, 179.6 mg GAEq/g, 32.78 mgQEq/g, 25.36 µg/mL and 123.24 µmol Fe 2+ /g for the extraction yield, TPC, TFC, DPPH-IC 50 , and FRAP-EC 50 , respectively ( Figure 5). Our results confirm that solvent concentration, sonication time, and solvent ratio affect all studied responses as obtained in other studies for the Malaysian propolis [55], Polish propolis [56], and Romanian propolis [57]. Processes 2022, 10, x FOR PEER REVIEW 12 of 17

Phenolic Screening of Optimized Extraction of Moroccan Propolis
To identify and quantify the individual phenolic compounds from Moroccan propolis obtained under optimal BBD conditions, the optimized ultrasound-assisted extract was submitted for HPLC-DAD analysis. As shown in Figure 6 and Table 5, twenty phenolic compounds were identified and quantified; the most abundant phenolic acid was epicatechin (0.1927 ± 0.0037 mg/g) followed by coumaric acid (0.0530 ± 0.0007 mg/g) and rosmarinic acid (0.0429 ± 0.0013 mg/g), and the lowest concentration of phenolics was obtained for vanillin (0.0021 ± 0.0014 mg/g). In a study conducted by Laaroussi and colleagues, the results of the phenolic compounds analysis in seven propolis samples collected from different areas of Morocco revealed the presence of sixteen compounds, including ellagic acid, apigenin, o-coumaric acid, kaempferol, and rosmarinic acid [11]. These results show the special characteristics for Moroccan propolis extract in comparison with those reported for Spanish propolis, which showed a high concentration of pinocembrin (482 mg/g), caffeic acid phenylethyl ester (594 mg/g), and kaempferol (189 mg/g), and Bulgarian propolis contains pinobanksin (0.0147 mg/g) of ethanolic extract and chrysin (0.1204 mg/g) of ethanolic extract [57]. In addition to the extraction operating conditions, this phenolic screening variation can be explained by the diverse flora foraged by the bees [58]. In turn, this variation in the composition of propolis may ultimately affect its biological and pharmacological activities [59,60].

Conclusions
In this study, the optimization of ultrasound-assisted extraction using response surface methodology was performed for the first time on Moroccan propolis. The optimization process has increased the recovery of extracts with high levels of polyphenols and flavonoids and with important antioxidant activity. In addition, using the Box-Behnken design showed clear relationships between all responses and interactions between the independent variables. The results obtained indicate that the sonication time, ethanol concentration, and solvent/material ratio were considered statistically significant and that they can influence the chemical and biological activities of propolis. Therefore, it was important to use an experimental design that increases the process efficiency and reduces costs and the number of experiments. Otherwise, our results can represent the main step in the control of phenolic compounds, such as an alternative to chemical antioxidants mostly used in the food industry and known for their toxicity.

Conclusions
In this study, the optimization of ultrasound-assisted extraction using response surface methodology was performed for the first time on Moroccan propolis. The optimization process has increased the recovery of extracts with high levels of polyphenols and flavonoids and with important antioxidant activity. In addition, using the Box-Behnken design showed clear relationships between all responses and interactions between the independent variables. The results obtained indicate that the sonication time, ethanol concentration, and solvent/material ratio were considered statistically significant and that they can influence the chemical and biological activities of propolis. Therefore, it was important to use an experimental design that increases the process efficiency and reduces costs and the number of experiments. Otherwise, our results can represent the main step in the control of phenolic compounds, such as an alternative to chemical antioxidants mostly used in the food industry and known for their toxicity.