Application of Multivariate Optimization for Phenolic Compounds and Antioxidants Extraction from Moroccan Cannabis sativa Waste

A statistical simplex centroid design methodology was applied to determine the eﬀects of diﬀerent solvents and their mixtures on the yield, total polyphenol content, 2 ′ 2-dipheny-l-picrylhydrazyl (DPPH) radical scavenging activity, and ferric reducing antioxidant power (FRAP) of extracts from the waste of Cannabis sativa . The diﬀerent extractor solvents (ethanol, methanol, water, and hexane) and their binary and ternary combinations were evaluated. The experimental results and their response surface models showed that the highest TPC yield values occur with the binary interaction between water and ethanol around the proportion of (ethanol, 70%; water, 30%). The desirability function showed that the optimal conditions were for TPC extraction ternary mixtures which consisted of 75% ethanol, 12.5% methanol, and 12.5% water. Ternary mixtures including water and binary mixture (ethanol 50% to 75%) yielded extracts with the best DPPH antioxidant activity, whereas pure methanol was the best solvent for extracting molecules with FRAP antioxidant capacity. The desirability function including all responses showed that the optimal solvent mixture consisted of 25% ethanol and 75% methanol.


Introduction
Cannabis sativa L. is a dioecious plant belonging to the Cannabaceae family. Cannabinoids, flavones, and terpenes are the main phytochemicals found in this plant [1]. is plant was used as a medicine before the Christian era in Asia, principally in India, and after it was introduced to Western medicine in the midst of the 19th century. Recently, the interest in C. sativa L. has drastically increased, due to psychoactive and nonpsychoactive compounds, namely, Δ9-tetrahydrocannabinol (Δ9-THC) and cannabidiol (CBD) [2]. is interest had increased in the 1990s through the description of cannabinoid receptors and the identification of an endogenous cannabinoid system in the brain [3]. In the last years, Cannabis sativa L. has been the center of the attention of the scientific community and definitely became one of the most studied plants. In Morocco, cannabis cultivation in the northern Rif region is currently undergoing its most significant evolution since the hashish industry emerged in the 1960s and dramatically developed in the 1980s [4]. To separate the golden beige powder called "Hashish" resulting from the solidification of the droplets of resin exuding from the tops of the female plants, these plants are first dried, and the resin and trichomes become dusty and more brittle; thus the resin can be separated from the plant material by beat or shake over a thin stretched nylon veil which acts as a sieve [5,6]. After this process, the plant residue (whole plant fragments) is separated from seeds and discarded.
Phenolic compounds are associated with a high number of biological activities including antioxidant capacity, which may help to protect the cells against the oxidative damage caused by free radicals [7]. Several novel extraction techniques have been developed in an attempt to obtain a more efficient extraction of target compounds by reducing both extraction time and used solvent [8]. e efficient recovery of bioactive polyphenolic phytochemicals, which occur widely in several agri-food residues that may serve as both food additives and bioactive substances in cosmetics and pharmaceuticals, is one of the higher value options [9]. ese secondary metabolites, which usually occur in low concentrations, are medicinally useful. e yields of these metabolites depend on the solvents and methods of extraction [10]. e quantity and quality of plant extracts are depending on the extraction protocol, which includes several factors such as the type of solvent, temperature, pH, the number of extraction steps, liquid-to-solid ratio, and the particle size of the solute contributing to the efficacy of the extraction process [11]. Solvent extraction is the most frequently used technique for the isolation of plant antioxidant compounds, due to the presence of different antioxidant compounds of varied chemical characteristics and polarities. Polar solvents like methanol and ethanol have been extensively employed to extract antioxidant compounds from various plants and plant-based foods [12].
Response surface methodology (RSM) is a wellestablished tool for the optimization of analytical methods, which is widely applied for the analysis of foods and herbal medicine [13]. e mixture design is a class of response surface experiments whose aim is to develop better or innovative formulations providing optimal requests and to create general conceptions about responses and interactions between independent factors allowing the modelization of the studied interaction [9]. To our knowledge, its use to study the effect of solvent combinations on cannabis waste has not been previously reported. e goal of this experiment is to perform the lowest number of experiments while gaining the maximum amount of data for the development of an efficient and reproducible model with the desired properties, to obtain clear information on the proportions of individual solvents including ethanol, methanol, water, and hexane in a mixture characterized by its capacity to yield extracts with the highest antioxidant activity and containing the greatest amounts of TPC from Cannabis sativa plant residue.

Materials and Methods
2.1. Extraction. Cannabis sativa plant residue without resin was collected from farmers and was ground into fine powder. e extracts were prepared in triplicate by adding 1 mL of the solvent (pure solvent and mixture) to 50 mg of plant residue powder and subjected to sonication for 30 minutes in an ultrasonic bath at ambient temperature, so the extracts were recovered after centrifugation for 10 minutes at 10,000 rpm and stored in the dark at 4°C.

Total Phenolic Compounds (TPC)
. Total phenolic content (TPC) was determined by spectrophotometry, using the colorimetric method, based on the Folin-Ciocalteu reagent [14] with modifications. 50 μL of extract or gallic acid standard solution was mixed with 450 μL of Folin-Ciocalteu reagent solution diluted 10 times. Subsequently, the mixture was shaken on a vortex mixer and left to incubate for 5 minutes at room temperature, and then 450 μL of a Na 2 CO 3 solution (75 g·L −1 ) was added and the mixture was shaken again. After incubation for 2 hours at room temperature their absorption was measured at 760 nm in a UV/visible Jenway 6505 scanning spectrophotometer. e gallic acid calibration curve (standard curve equation, y � 2.8388x + 0.0556, R 2 � 0.9994) was prepared with a concentration ranging from 0.062 to 1 mg·mL −1 in ethanol.
e results are expressed in gallic acid equivalent (GAE) g −1 dry plant mg. All analyses were done in triplicate.

DPPH Free Radical Scavenging Activity.
e free radical scavenging activity of all extracts was performed following the procedure described by Brand-Williams et al. [15]. Briefly, 25 µL of each extract (with different dilutions) was added to 1 mL of an ethanol solution of DPPH (60 µM). e absorption measurements were made at 515 nm, after the incubation time of 60 minutes at room temperature. e absorption of a blank sample containing the same amount of ethanol and DPPH solution served as a negative control. e percentage inhibition of activity trapping free radicals in each extract was calculated as follows: Triple measurements were made for each sample and the results were expressed in terms of average ± standard error.

Ferric Reducing Antioxidant Power (FRAP).
e reductive capacity of liquid extracts was determined by the reducing power test described by Oyaizu [16]. is method is based on the transformation of ions, Fe 2+ into Fe 3+ by polyphenol antioxidants. 0.2 mL of the extract was mixed with 1 mL of 0.2 M phosphate buffer and 1 mL 1% potassium ferricyanide in glass tubes. e prepared reaction mixture was incubated for 20 minutes at 50°C in a water bath. After incubation, 1 mL of 10% trichloroacetic acid solution was added to the reaction mixture. Tubes were centrifuged at 3000 rpm for 10 minutes after 2 mL of the supernatant was mixed with 2 mL of bidistilled water and 0.4 mL of 0.1% ferric chloride solution. e absorbance was measured at a wavelength of 700 nm. e reductive capacity of the extracts is expressed as ascorbic acid equivalent (Acs. E) g −1 dry plant. All analyses were done in triplicate.

Experimental Design and Optimization.
Mixture models were revealed to be useful tools in describing the behavior of solvent mixtures to extract bioactive compounds from plants [17]. In the present study, a simplex centroid design with a response surface methodology was employed to determine the effect of interactions among extraction solvents on TPC and antioxidants activities to optimize the extraction conditions. Mixture design experiments were designed and analyzed using the free version StatSoft, Inc. (2011), STA-TISTICA (data analysis software system), version 10. A total of seven combinations are used. Results were expressed as mean ± SD (RSD (%)) and were analyzed using ANOVA. Significant differences were determined by Tukey's test, with p < 0.05 as the significance criterion. e following 2 Journal of Chemistry polynomial equation of function Xi was fitted for each factor assessed at each experimental point: where Y is the predicted response and β1, β2, β3, β12, β13, and β23 are constant coefficients for each linear and nonlinear interaction term. Principal component analysis (PCA) was carried out to understand the influence of the extraction solvents on the antioxidant compounds and identify the relationships between the variables.

Results and Discussion
3.1. Solvent Screening. Solvent selection is one of the most important steps in extraction of phenolics and other bioactive compounds from fruits, vegetables, byproducts, and their recovery from plant-based materials [3]. Generally, desired compound extraction efficiency is influenced by multiple parameters, such as temperature, time, and solvent polarity, and their effects may either be independent or interactive [18]. e results of total phenolic contents and the antioxidant activity of the seven pure solvent extracts are shown in (Figure 1). It is obvious to note that the phenolic compounds and antioxidant capacity of C. sativa extracts were highly and significantly influenced (p < 0.05) by the type of extraction solvent. e same conclusions have been reported by Liyanapathirana and Shahidi [19]. Estimated values of total phenolic content in different extracts yielded using different solvents ranged from 1.90 ± 00 to 19.07 ± 00 mg GAE/g DP, pointing out that ethanol extract was done ten times more TPC than ethyl acetate. Other works corroborating our findings reported TPC results for C. sativa ranging from 5.85 to 17.05 mg GAE/g dw [20]. However, higher TPC in inflorescence samples from industrial hemp was quantified in the range of 10.51 to 52.58 mg GAE/g [21].
Ethanol extract exhibited the highest phenolic content followed by methanol extract. Wile, water, and hexane come, respectively, in third and fourth place, whereas dichloromethane was the less efficient solvent for extraction of TPC from Cannabis sativa waste. According to [22], the extract obtained by 100% ethanol from L. aromatica, showed the highest phenolic content and exhibited the highest total antioxidant activity, reducing power, and DPPH radical scavenging activity. e seven solvents were also screened using two antioxidant methods. According to the total antioxidant activity ( Figure 1), methanol extracts demonstrated the best antioxidant activity (73.47 mg Asc.E/g DP), followed by ethanol. Conversely, dichloromethane extracts exhibited the lowest antioxidant activity. As shown in Figure 1, the DPPH free radical scavenging activity revealed that ethanol and methanol were the best solvents for extracting molecules endowed with antioxidant properties followed by hexane and then water. According to this screening test, methanol, ethanol, water, and hexane were the best solvents for extracting phenolic compounds and other molecules with antioxidant activity, which allowed us to select them to perform the solvent mixture design.

Mixture
Design. RSM is considered the best statistical approach to assess the influence of different experimental factors and their linear and quadratic interactions on analytical responses and to model and optimize the extraction conditions of phenolics that display important biological activities [3]. is approach gathers the maximum amount of information in the minimum number of analyses. e first mixture design was perfumed using hydroalcoholic solvents by mixing water, ethanol, and methanol, whereas the second was done using three organic solvents, ethanol, methanol, and hexane. e results of the effects of the solvent composition on the yield of TPC and antioxidant activities are shown in Table 1.

Analysis of Variance (ANOVA).
To improve the extraction conditions allowing the best antioxidant compounds recovery, the simplex centroid mixture design was used to construct a contour plot by regression model analysis. Experimental data were fitted to the linear, quadratic, and special cubic models. e statistical significance of the model and equation terms was analyzed based on p value (Prob > F). e goodness of fit of the model to justify its robustness was evaluated by the coefficient of determination (R 2 ) and the adjusted correlation coefficient (Adj-R 2 ), whereas the model's significance was checked using the Ftest. e analysis of variance (ANOVA) results based at 95% confidence interval, including regression model terms, R 2 , Ftest, and probability values are depicted in Table 2. In all mixtures, the linear model explained the variance up to good levels, and expanding from linear to a quadratic and special cubic model improved the fit for the regression analysis. Concerning the hydroalcoholic mixtures, the coefficient of determination (R 2 ) value of the quadratic model was (R 2 � 1) and the adjusted R 2 was 0.92 for both TPC and DPPH assays. In the case of FRAP assay, R 2 and adjusted R 2 were equal to 0.96 and 0.95, respectively, which indicated that the model adequately represented the real relationship between the parameters chosen. Furthermore, those two coefficients were equal to 1 (R 2 � 1) for the three responses in the special cubic model. Hence, this special cubic model was chosen since it could explain 100% of the variability for TPC, DPPH, and FRAP assays, as well as presenting small p values (p value <0.000001) and high F values (Table 3) indicating that this model was very good in predicting the behavior of the mixtures.
Regarding the mixtures where only organic solvents (ethanol, methanol, and hexane) were used, ANOVA indicated that the quadratic model explained 91%, 100%, and 99% of the total variation of the results, for TPC, DPPH, and FRAP, respectively, but this model was not significant for total phenolic contents (p value �0.08), whereas the special cubic model explained the result's variability at 100% level for all responses, besides presenting small p values (p value <0.000001) and high F values (    In general, a positive sign for the coefficient in the fitted model indicates the ability of the variable to increase the response, whereas the negative sign indicates the ability of a variable to decrease the response [17]. e equation model showed that the phenolic content and antioxidants activities (DPPH and FRAP) were positively and linearly influenced by ethanol, methanol, and water, respectively, noting that, in linear trials, water resulted in the lowest coefficient, thus, the smallest TPC amounts and the lowest antioxidant activities. Among the binary interactions, water and ethanol interactions resulted in the highest positive effect in all responses, followed by water-methanol interactions for TPC and DPPH, whereas, in the case of FRAP assay, this coefficient was negative in the binary interactions of methanol-ethanol and water-methanol indicating an antagonistic effect between the two solvents. In addition, the ternary mixtures influenced positively the extraction of TPC and antioxidant compounds; this effect was much higher in the FRAP assay. e regression equations obtained for the special cubic model in organic mixtures correlating the three variables and the analytical responses are as follows: Regarding the linear interaction, in contrast to hexane, ethanol exhibited the highest positive effect on TPC extraction and DPPH free radical scavenging activity, while methanol demonstrated the highest positive influence in the FRAP assay. e strongest positive effect of binary interaction was found in ethanol-hexane mixtures for TPC and DPPH and methanol-hexane mixtures for FRAP assay, whereas the ternary mixtures had high and negative influence and displayed an antagonistic effect among the components of the mixture for TPC extraction and FRAP assay of the yielded extracts. ese results clearly indicate that the extraction is affected by the polarity of the solvents.

Total Phenolic Compounds.
Among the major classes of plant chemicals, namely, terpenoids, phenolics, and alkaloids phenolic compounds including phenolic acids (hydroxybenzoic and hydroxycinnamic acids), polyphenols (hydrolyzable and condensed tannins), and flavonoids are the most important for dietary applications and the most extensively researched [22]. C. sativa inflorescence could be considered as a potential novel source of polyphenols intended for nutraceutical formulations [21]. e recovery of phenols depends on the solvent system used since the solvent mixture extracts higher TPC amounts than did pure solvent [23]. us, the development of a selective method for the extraction of phenolic compounds is required. e threedimensional (3D) interaction contour plots are given in Figure 2(a) as a function of the ethanol, methanol, and water interactions. ese three-dimensional response surface plots can exhibit how the initial factors affected different responses more clearly. e interpretation of contour plots in Figure 2(a) demonstrated the total polyphenols extraction using the combination of ethanol and water resulted in the highest antioxidant activities, which can confirm a synergistic effect between these two solvents, in which the interaction between them demonstrates a better result than their isolated actions. Adding water to organic solvents improves extraction rate, but too high water content brought an increased concomitant extraction of other compounds and then lower phenols concentrations [23]. We also noticed that joining water to methanol enhanced its ability to extract phenolic compounds. e highest TPC yield value on the contour graph is seen to occur with the binary interaction between water and ethanol around the proportion of ethanol 70%, water 30%, whereas the best solvent stated by other authors, for extraction of phenolic compounds from hemp, was 50% ethanol [20]. According to Wong et al. [24], increasing ethanol concentration from 20% to 80% was favorable to increase total phenolic content recovery from palm kernel byproducts from 7.8 to 11.5 mg/g GAE.
According to the desirability analysis (results shown in Table 4), the best solvent combination for optimal extraction of the phenolic compound is the ternary mixture which consists of 75% ethanol, 12.5% methanol, and 12.5% water. Total phenolic compound contour plot, extracted using organic solvent mixtures including ethanol, methanol, and hexane, is shown in Figure 3(a). As shown in Figure 3(a), ethanol extracted the highest amounts of TPC, while hexane was the worst solvent for TPC recovery, which rapidly increases when the proportion of the ethanol increases in the  Journal of Chemistry binary mixture. e same observation can be noticed also when adding methanol to hexane, but with less effectiveness. us, adding ethanol and methanol to hexane increases its capacity to extract TPC. Desirability analysis (Table 4) revealed the equivalent binary interaction to be the optimal solvent mixture in the absence of water.

DPPH Free Radical Scavenging Activity.
Generally, plants bioactive compounds with higher polar character can be easily extracted with water; however highly hydroxylated phenolic compounds, such as catechins, are more soluble in alcohols such as ethanol and methanol [25]. e DPPH assay is employed to test the ability of compounds to act as free radical scavengers and is frequently used to evaluate the antioxidant capacity of foods. e recovery of antioxidants from plant residues is interesting from a technological point of view as valuable components of nutraceuticals in food and pharmaceutical preparations or the cosmetics industry [26].
To investigate the effects of solvents and their interactive effects on antioxidant activity measured by DPPH assay of the extracts, the three-dimensional plots were depicted in Figures 2(b) and 3(b). Concerning the hydroalcoholic mixtures, as evident in Figure 2(b), the DPPH activity values of water extract increased when the methanol or ethanol concentration increases, indicating that adding these organic solvents to water rises its capacity to extract phenolic compounds, although it starts to decrease when ethanol or methanol concentration in the solvent mixture exceeds 75%. e solvent mixtures yielding extracts with the highest DPPH free radical scavenging activity occur between the equivalent ternary mixture, the equivalent binary mixture water-ethanol, and the binary mixture (water 25%, ethanol 75%). Moreover, the optimum solvent mixture given by the desirability analysis (Table 4) was a ternary mixture that consisted of 50% ethanol, 16.667% methanol, and 33.33% water.
Similarly, in the organic mixtures (Figure 3(b)) adding methanol or ethanol to hexane increases its capacity to extract molecules with DPPH free radical scavenging activity. However, in this case the quadratic interaction between methanol and hexane was more effective. e best solvent mixture for yielding extract with optimal DPPH antioxidant activity was the equivalent binary mixture of methanol and ethanol.  e same sightings have been observed when using hexane instead of water. FRAP desirability analysis demonstrated that pure methanol was the optimal solvent for extracting bioactive compounds with FRAP antioxidant activity. Regarding FRAP activity, methanol proportion was one of the most important variables affecting extraction from ymelaea hirsuta L [27]. According to the planned mixture design performed to extract total anthocyanins, phenolic compounds, and antioxidants from x´kijit peels, methanol was the solvent producing extracts with the highest extraction yields of total phenolic compounds and other bioactive compounds, exhibiting the highest antioxidant capacity [28].

Pareto Chart Analysis.
To examine the relative importance of the main effects and their interactions using statistical significance (p < 0.05), a standardized Pareto chart was employed.
According to Figures 4(a) and 4(b), the result showed the effect of the independent variables (hydroalcoholic mixtures) and their interactions on TPC recovery and DPPH antioxidant activity. We notice for the two responses that all parameter was statistically significant with 95% confidence and had a positive effect on the two responses and was ranked as follows: ethanol > methanol > water > ethanol * water > methanol * water > ethanol * methanol, whereas in the case of FRAP assay (Figure 4(c)), methanol had the greatest significant positive effect followed by ethanol. e binary interaction methanol * water had a low significant negative effect, meaning that adding water to methanol slightly lowered the recovery on molecules having a FRAP antioxidant activity.
Concerning organic solvent (figures not shown), the Pareto analysis revealed that all parameter was statistically significant with 95% confidence. In the linear interactions, mostly ethanol exhibited the highest positive influence, followed by methanol; likewise, the quadratic interactions influenced positively the TPC recovery, but with smaller impact, in contrast to ternary interaction which had a negative effect. In the DPPH assay, methanol and ethanol greatly influenced extractions of antioxidant molecules, followed by hexane with lesser effect. e interaction of ethanol-methanol showed a significant negative influence. Methanol was the pure solvent that highly influenced FRAP antioxidant molecules, followed by ethanol.   Figure 5) show that 94.84% of the variability in the data is accounted by the first two principal components PC1 (Factor 1) and PC2 (Factor 2) which was high enough.   Furthermore, the first three PCs explained 100% of the variance in the data, which was high enough to represent all the variables. In this analysis, if two vectors subtend a small angle to each other, this means that the two variables they represent are strongly correlated. e score plot for PC1 versus PC2 clearly shows a strong positive correlation between TPC and DPPH was observed, indicating that fractions with the highest potency in scavenging DPPH radical presented the highest TPC. e two parameters were in turn positively correlated with ethanol concentration in the solvent extractor and negatively correlated to hexane. On the other hand, FRAP shows a positive correlation with methanol concentration in extraction solvent and was negatively correlated with water content, indicating that those ferric reducing molecules have more affinity to methanol than water and ethanol.

Desirability Analysis.
Optimization using the desirability function was performed in order to maximize the total phenolic compounds and also to maximize the antioxidant capacity (DPPH and FRAP) in the extracts. Results are depicted in Figure 6 for hydroalcoholic mixtures. e final result for the simultaneous optimization including all responses, using the desirability function, suggested that the solvent mixture consisting of 25% ethanol and 75% methanol was the most adequate solvent to achieve the best solution for this combination of variables. However, pure methanol was the solvent when mixtures were performed with the three organic solvents.

Conclusions
is is the first report on the optimization of phenolic antioxidants from Cannabis sativa residue, in which we revealed that this plant's waste is a useful source for bioactive compounds. e methods can be utilized for further isolation of active fraction/compounds from the cannabis waste. e highest TPC yield values occur with the binary interaction between water and ethanol around the proportion of ethanol, 70%, water, 30%, whereas the best DPPH scavenging activity is obtained with extraction using ternary mixtures including water and binary mixture (ethanol 50% to 75%). Meanwhile, pure methanol was the best solvent for extracting molecules with FRAP antioxidant capacity.

Data Availability
e table data used to support the findings of this study are included within the article.

Conflicts of Interest
e author declares that he has no conflicts of interest.