Maximizing Contents of Phytochemicals Obtained from Dried Sour Cherries by Ultrasound-Assisted Extraction

: Sour cherries were ﬁrst dried by vacuum drying and then used as material for obtaining extracts rich in bioactive compounds by ultrasound-assisted extraction (UAE). The ﬁrst step was to apply a factorial design for the preliminary experiments to determine the most inﬂuential UAE factors, and thus the three studied parameters were chosen as the most suitable for the design of the main experiment (temperature, liquid–solid ratio and ethanol concentration). In this part, the contents of total phenols and the total content of monomeric anthocyanins were taken for responses. For the further optimization of UAE, experimental design (face-centered) was applied, and the yield, total phenolics, ﬂavonoid content and content of monomeric anthocyanins and antioxidant activity (DPPH, ABTS and FRAP assays) were analyzed. Temperature (40–80 ◦ C), ethanol concentration (40–80%, w / w ) and liquid–solid ratio (10–30 mL/g) were investigated as independent variables. The obtained experimental results were ﬁtted to a second-order polynomial model and analysis of variance was used to determine the ﬁt of the model and the optimal conditions for investigated responses. High quality extracts with high concentrations of polyphenols and anthocyanins were also obtained, which could be used as food additives.


Introduction
The sour cherry (Prunus cerasus L.) is one of the fruits that originated in areas around the Caspian Sea. Russia, Turkey, Ukraine, USA, Iran and Serbia are the leading producing countries of sour cherries [1]. All good references for the possibility of growing various continental stone fruits Serbia owes to a suitable climate and cultivation possibilities at different altitudes [2]. Sour cherries are considered one of the traditional fruit species, as they are widely grown and highly appreciated in the region [3].
Sour cherries are characterized as an important industrial fruit variety, and they are a valued raw material in the area of food production. Thanks to their high postharvest respiration and susceptible structure, numerous methods must be applied in the fruit processing industry to control and maintain adequate quality of sour cherries [4]. Many different cherry varieties can be processed into different sour cherry products, such as dried or frozen products, jams, juices or jellies [5,6], purees, concentrates and alcoholic beverages [7], and these products can also be a common additional ingredient in the industrial production of chocolate and other sweets [3]. Fruit drying has been successfully applied to reduce their high moisture content and to achieve and maintain microbiological stability of the obtained dried products [5]. Fresh sour cherries can be dried by convective, vacuum and freeze drying. Structural and biochemical changes in the products during drying were significantly reduced in vacuum drying compared to convective drying, with great emphasis on low energy use for economic efficiency [8]. Recently, the results of the comparative analysis of the effects of convective drying, vacuum drying and lyophilization (freeze drying) on the physical, chemical and biological properties of stone fruits in the research of Vakula et al. [9] showed that the highest total phenolic and flavonoid content of all dried stone fruit samples were obtained in the vacuum dried sour cherry samples.
Sour cherries are characterized by excellent sensory and physico-chemical properties as well as high levels of biologically active phytochemicals [1]. The presence of polyphenolic compounds, such as phenolic acids and isoflavonoids, as well as the unique composition of anthocyanins, in sour cherries have beneficial effects on human health and are also associated with their ability to scavenge oxygen free radicals in human cells [6,10,11]. In addition, sour cherries are rich sources of various organic acids, vitamins and minerals [12]. Extraction of the widely present bioactive compounds from the dried plant samples has become a significant area of research due to their priceless bioactive value [13].
Conventional extraction techniques have implied a solvent extraction process using various organic solvents and water. Since the plant raw material contains both extractable and non-extractable bioactive compounds, this technique cannot extract all the valuable antioxidants present. The bound form of these compounds remains unextracted and cannot be analyzed [11]. To this end, several novel extraction techniques, such as ultrasound-assisted extraction (UAE), supercritical fluid extraction (SFE), microwave-assisted extraction (MAE) and pressurized-liquid extraction (PLE) have been successfully used for the extraction of valuable biological components from plant materials with high yield [14]. SFE is a very good alternative to conventional processes because it can prevent chemical changes caused by high temperatures and also by interaction with solvents and/or oxygen. It is also important to consider the right conditions for the extraction, mainly pressure and temperature, e.g., in the case of SFE, where large overall extraction yields can be obtained, but at the same time lower selectivity when higher pressures (up to 500 mbar) are used [15]. One of the crucial facts that enables high efficiency when using ultrasound for extraction is the increased contact between the liquid and solid phases. The facilitated solvent penetration into the sample matrix provides higher extraction yield and quality with reduced process time [16][17][18]. Indeed, the accelerated mass transfer by ultrasonic waves and the resulting cavitation effect during the process are responsible for the release of bioactive compounds from the parenchymal cell wall of the samples [19]. A very important part of the process is the selection of suitable solvents for extraction in accordance with its selectivity and solubility, which will correspond to the targeted extractable compounds. In comparison with methods of extraction that are conventionally used, the application of efficient ultrasound extraction requires less power and solvent consumption, less fossil energy consumption and also excludes the occurrence of wastewater after the end of the extraction process [14]. The importance and advantage of using UAE compared to MAE are presented in the work of Liu et al. [20], where the screening and process optimization of ultrasound-assisted extraction of major antioxidants from sweet tea were studied. On the other hand, a method for the MAE of natural antioxidants from Akebia trifoliate peel was developed in the work of Luo et al. [21].
Nowadays, the emergence of advanced "green" extraction technology increases the competition between large processing companies to be more economical and innovative, as well as to take all necessary precautions to protect the environment as much as possible [22]. The advantages of using UAE can be explained by optimizing important variables that include solvent concentration, solvent-sample ratio, temperature and extraction time [23]. In order to optimize the mentioned variables, response surface methodology (RSM) can be successfully applied. RSM represents a widely known set of statistical and mathematical methods used for the analysis of problems where the influence of multiple input variables on output parameters is observed. The output parameters are referred to as responses of interest or dependent variables. Optimization of the process involves finding an appropriate combination of input factors from given intervals, for which the desired output values are best achieved with the required priority [24]. RSM, as an experimental methodology, has been successfully used for the optimization of ultrasound extraction processes of various raw materials [23,[25][26][27].
The application of RSM on UAE of different classes of bioactives is investigated and described in numerous studies. However, among the available and scientific data sources, there are no works dealing with optimization of the extraction process of phytochemicals from sour cherries. Chemical and biological compounds (total phenolic content (TP), total flavonoid content (TF) and anthocyanins content (TMA)) from sour cherries and their antioxidant activity (measured by DPPH, FRAP and ABTS assay) were investigated. The main goal of this research was optimization of the UAE of valuable bioactive compounds and antioxidants from vacuum dried sour cherries using RSM.

Sample
Sour cherries were purchased from a local market (Novi Sad, Serbia) and after purchase the pits were carefully removed, and the samples were stored in a freezer (−20 • C) until the drying process. Before the extraction process, the samples were dried in a vacuum dryer. Drying was performed in a prototype vacuum dryer designed and installed at the Faculty of Technology Novi Sad, Serbia. This dryer was described in detail in the work of Šumić et al. [28]. Drying conditions were chosen according to previous studies on stone fruit drying presented by Vakula et al. [9], where it was found that vacuum dried sour cherry samples had the highest total phenolic and flavonoid content of all dried stone fruit samples. In accordance with these results, vacuum drying in this study was carried out at 60 • C, 20 mbar and 11 h. The moisture content of the vacuum dried sour cherry sample was 16.53%.

Preliminary Ultrasound-Assisted Extraction (UAE) and Screening of Variables
In all UAE experimental runs, 5.0 g of vacuum dried sample was mixed with extraction solvent (aqueous ethanol) in 250 mL glass flasks. Extractions were performed in a sonication water bath (EUP540A, Euinstruments, France) at a frequency of 40 kHz. The flasks were always positioned in the position of the ultrasonic bath to ensure constant ultrasonic power, according to the procedure described by Zeković et al. [27].
The first step was to screen the independent variables to define the most influential parameters and their domain on already targeted responses (TP and TMA) using a 2 5-1 fractional factorial design. The experimental design consisted of 16 runs using 40 and 60 • C temperature, 20 and 40 min extraction time, 40 and 60% ethanol concentration, 30 and 60 W/L ultrasonic power and 10 and 20 mL/g liquid-solid ratio as numerical input parameters. The following linear model presented in Equation (1) was used to define the influence of parameters of UAE on TP and TMA: where Y represents the response variable, β 0 the intercept, β i the linear regression coefficient, β ij the regression coefficients for cross-product terms and X i and X j variables (independent) that affect the response.

Experimental Design and Statistical Analysis
After the screening, three of the five most influential parameters were selected to be further used in face-centered CCD with RSM. The influence of temperature (40, 60 and 80 • C), extraction time (40, 60 and 80 min) and liquid-solid ratio (10, 20 and 30 mL/g) were used as independent variables. Optimal extraction conditions were determined considering Y and TP, TF, TMA, as well as antioxidant activity parameters obtained by DPPH, ABTS and FRAP assays, while the selection of optimal conditions was based on the desirability function (D). For multiple linear regression analysis, Design-Expert v11 software (Stat-Ease, Minneapolis, MN, USA) was used and the results were fitted to a second-order polynomial model (Equation (2)): where Y represents the response variable, X i and X j are the independent variables affecting the response, and β 0 , β i , β ii , and β ij are the regression coefficients for intercept, linear, quadratic and cross-product terms. Model adequacy was evaluated by the R 2 (coefficient of determination), CV (coefficient of variation) and p values for the model and lack of fit. In order to verify obtained empirical models, validation was performed by using the extracts prepared at optimized UAE conditions.

Yield
The content of yields in the obtained sour cherry extracts was obtained by the process of vacuum vaporization using 10 mL of crude extract. After vaporization, drying was carried out in an oven at 105 • C until a constant mass was obtained. The results are expressed as mass of total extractable solids per 100 g of dry plant material (%; w/w).

Total Phenolic Content
The total phenolic content (TP) of each extract was estimated by the Folin-Ciocalteu method [29] using gallic acid as the standard equivalent. The absorbance of the bluecolored reaction solutions was read at 750 nm wavelength (6300 Spectrophotometer, Jenway, UK). Total phenolic content was calculated as g gallic acid equivalents per 100 g sample (g GAE/100 g) using the prepared gallic acid calibration curve. Each measurement was carried out three times and results were presented as mean values.

Total Flavonoid Content
For determination of total flavonoid content (TF), the colorimetric aluminum chloride assay was used [30]. Absorbance values were observed at 510 nm. Results were expressed as g of catechin equivalents per 100 g of dried sour cherries (g CE/100 g) using the equation obtained from the catechin standard diagram. Each measurement was performed in three replicates and the results were expressed as mean values.

Total Anthocyanin Content
Total monomeric anthocyanins (TMA) content in dried sour cherry extracts was determined according to the pH differential method described by Fuleki and Francis [31]. First, two buffer solutions with pH 1.0 (potassium chloride buffer) and pH 4.5 (sodium acetate buffer) were prepared. Then, 400 µL of each liquid extract was mixed with 3.6 mL of the two buffer solutions separately. The absorbance of the mixture was recorded at wavelengths 510 and 700 nm, using distilled water as a blank. Absorbance (Ad) was calculated according to Equation (3): Total monomeric anthocyanin content in the extracts was calculated according to the following Equation (4): where A d is difference in absorbance, MW is a molecular weight of cyaniding-3-glucoside (449.2 g/mol), DF is the dilution factor of the samples and M a is the molar absorptivity of cyaniding-3-glucoside (26.900 M/cm). Finally, the results obtained were expressed as mg of cyanidin-3-glycoside equivalent (CGE) per 100 g of sample (mg C3G/100 g).
Each measurement was performed in three replicates and the results were expressed as mean values.

DPPH Assay
To determine the antioxidant activity of methanolic sour cherry extracts, a slightly modified DPPH (2,2-diphenyl-1-picrylhydrazyl) free radical scavenging assay was applied, mainly described in the research conducted by Brand-Williams et al. [32]. The methanolic DPPH solution was freshly prepared to perform the analysis. Methanol was used to adjust the final absorbance of the DPPH reagent to 0.70 (±0.02). A total of 0.1 mL of the liquid extract was added to 2.9 mL of DPPH reagent. For the blank sample, an appropriate volume of solvent (methanol) was added instead of 0.1 mL of the sour cherry extract. After incubation for 60 min at room temperature, the absorbance of each sample was observed at 517 nm by using VIS spectrophotometer (6300 Spectrophotometer, Jenway, UK). Trolox was used as an ethanolic standard solution to obtain a calibration curve, and the calculated results have been presented as µM Trolox equivalents per g of dried sour cherries (µM/g). Each measurement was performed in three replicates and the results were expressed as mean values.

FRAP Assay
The potential of ferric reducing antioxidant (FRAP) assay was conducted with the aim of measuring the reducing power of the sour cherry extracts against positive trivalent ferric ion (Fe 3+ ) [33]. The freshly prepared FRAP reagent contained 10 mM TPTZ (2,4,6-tris(2-pyridil)-s triazine) in 40 mM HCl, 20 mM iron (III) chloride (FeCl 3 ) aqueous solution and 300 mM acetate buffer, pH 3.6. The working FRAP reagent just prior to assay was prepared by mixing these three solutions in a 1:1:10 (v/v/v) ratio. Appropriate volumes of sour cherry extracts (0.1 mL) were mixed with 2.9 mL FRAP reagent in glass cuvettes. The absorbance values of the reaction mixture at 593 nm were determined spectrophotometrically (6300 Spectrophotometer, Jenway, UK) after incubation for 10 min in the dark, at 37 • C. Freshly prepared aqueous Fe 2+ (Fe 2 SO 4 ) solution at different concentrations (0-0.23 mM, R 2 = 0.999) was taken to construct the standard calibration curve. The obtained results were presented as µM Fe 2+ equivalents per g of dried sour cherries (µM Fe 2+ /g). Each measurement was performed in three replicates and the results were expressed as mean values.

ABTS Assay
The antioxidant activity assay with 2,2 -azino-bis-(-3-ethylbenzothiazoline-6-sulfonic acid) diammonium salt) (ABTS) was carried out in accordance with modified methodology described by Re et al. [34]. The first step to obtain the ABTS reagent involved reaction of 7 mM aqueous ABTS solution with 2.45 mM aqueous potassium persulfate solution in the dark for 16 h at room temperature. Before the assay, this solution was diluted with acetate buffer (pH 3.6) and the absorbance of the prepared ABTS reagent was adjusted to 0.70 (±0.02). Then, the reaction mixture obtained by adding 2.9 mL of ABTS reagent into 0.1 mL of diluted extracts was incubated for 300 min at room temperature in the dark. After completion of incubation, absorbance values were recorded at a wavelength of 734 nm. Trolox was used as a reference standard to establish the calibration curve. The ability of the sour cherry extracts to scavenge ABTS radical cation was expressed as mM Trolox equivalents per g of sample dry weight (mM TE/g). Each measurement was performed in three replicates and the results were expressed as mean values.

Screening of UAE Factors
The UAE factors' influence during the extraction process has been thoroughly investigated in different fields and for different types of fruit and vegetables, which is described in detail by Kumar et al. [35] and also the obtained results of the different research have been presented in papers by Zou et al. [36], Medina-Torres et al. [37] and Dumitraşcu et al. [38]. The first goal in this research was to reduce the initial experiments resulting from all combinations of the five most important UAE factors (temperature, extraction time, ethanol concentration, ultrasonic power and liquid-solid ratio). This was carried out by a 2 5-1 factorial design that was firstly applied for screening the most influential UAE factors (Table 1). Total phenolic content and total monomeric anthocyanin content were selected as the most appropriate quality indicators and investigated responses. According to Pareto charts, where t values and Bonferroni limits are observed, the influence of the investigated UAE parameters on the contents of total phenolics and monomeric anthocyanins was shown ( Figure 1). According to t values that were observed for total phenolic content, it was found that temperature and liquid-solid ratio were above the Bonferroni limit (17.28), while ethanol concentration and extraction time were slightly below this limit. It was also seen that the ultrasonic power was significantly below the Bonferroni limit but also above t value limit (4.30). All five parameters studied were over the t value limit (Figure 1a). According to this, it was concluded that temperature and liquid-solid ratio had a significant effect on the response, while the concentration of ethanol was on the third place, slightly below the limit. In terms of total monomeric anthocyanin content, temperature, extraction time, ultrasonic power and liquid-solid ratio were above the Bonferroni limit and also t value limit, while the concentration of ethanol was under both limits ( Figure 1b). Based on these results, the three parameters studied, temperature, liquid-solid ratio and ethanol concentration, have been selected as the most important for further experimental design to optimize the UAE process. Linear terms of all three selected parameters showed a positive influence on the total phenolic content.
Accordingly, the following experimental domain of three chosen parameters for RSM study were: temperature (40, 60 and 80 • C), ethanol concentration (40, 60 and 80%) and liquid-solid ratio (10, 20 and 30 mL/g). Similar screening of significant variables was conducted in the research by Chen et al. [39], where ultrasound-assisted extraction of Rubia sylvatica Nakai fruit was performed and the impact of each factor on the yields of total anthocyanins and total phenolics was assessed according to the t value on Pareto charts. liquid-solid ratio had a significant effect on the response, while the concentration of ethanol was on the third place, slightly below the limit. In terms of total monomeric anthocyanin content, temperature, extraction time, ultrasonic power and liquid-solid ratio were above the Bonferroni limit and also t value limit, while the concentration of ethanol was under both limits ( Figure 1b). Based on these results, the three parameters studied, temperature, liquid-solid ratio and ethanol concentration, have been selected as the most important for further experimental design to optimize the UAE process. Linear terms of all three selected parameters showed a positive influence on the total phenolic content. Accordingly, the following experimental domain of three chosen parameters for RSM study were: temperature (40, 60 and 80 °C), ethanol concentration (40, 60 and 80%) and liquidsolid ratio (10, 20 and 30 mL/g). Similar screening of significant variables was conducted in the research by Chen et al. [39], where ultrasound-assisted extraction of Rubia sylvatica Nakai fruit was performed and the impact of each factor on the yields of total anthocyanins and total phenolics was assessed according to the t value on Pareto charts.

Defining the Influence of Process Parameters Extraction Process Modeling
After the previous screening of UAE factors, it was found that three investigated parameters (temperature, liquid-solid ratio and ethanol concentration) showed the greatest effect on total phenolics and monomeric anthocyanin content throughout the ultrasoundassisted extraction of vacuum dried sour cherry samples. Accordingly, these parameters were selected as the most suitable for the further main experimental design. The quality

Defining the Influence of Process Parameters Extraction Process Modeling
After the previous screening of UAE factors, it was found that three investigated parameters (temperature, liquid-solid ratio and ethanol concentration) showed the greatest effect on total phenolics and monomeric anthocyanin content throughout the ultrasoundassisted extraction of vacuum dried sour cherry samples. Accordingly, these parameters were selected as the most suitable for the further main experimental design. The quality indicators, i.e., the investigated responses besides total phenolic and monomeric anthocyanin contents, were also yield content, total flavonoid content and antioxidant activity (DPPH, FRAP and ABTS assays). For the optimization of UAE, experimental design (facecentered) was observed with the main goal of maximization of extraction yield (Y) and the investigated chemical (TP, TF and TMAC) and biological properties (antioxidant activity determined by three assays: DPPH, ABTS and FRAP) of the vacuum dried sour cherry extracts. The parameters of UAE and obtained results of examined responses are shown in Table 2.  (Table 3). These results indicated good fit between the values predicted experimentally and by the model that was applied. Regarding the values of CV obtained for all the investigated responses, presented in Table 3, it was noticed that a good fit of the model applied was obtained for all responses due to relatively low CV (1.68, 3.03, 4.80, 6.34, 10.04, 4.46 and 4.61 for Y, TP, TF, TMAC, DPPH, ABTS and FRAP, respectively), since low CV values indicate the good accordance and reliability of experiments and practical works [40]. According to the results of ANOVA presented in Table 3, a good fit between experimental and predicted values was observed. ANOVA was also used to determine the p values for linear, interaction and quadratic terms (Table 4).

Total Extraction Yield (Y)
The yield content of the extracts analyzed varied from 64.69 to 75.66% ( Table 2). The range of experimental yield of extract of Megalayan cherry fruit by ultrasonically assisted extraction extract obtained in the work of Kashyap et al. [41] was between 29.87 and 43.71%. The lowest yield content was obtained at 80 • C, 80% and 10 mL/g, as was also the case with total phenolic content, and the highest yield content was obtained at 40 • C, 80% and 30 mL/g ( Table 2). These results indicate the higher influence of higher liquid-to-solid ratio on the yield in extracts of vacuum dried sour cherries, compared with temperature and ethanol concentration during extraction. Moreover, with respect to the linear and interaction terms, all three independent variables, temperature, ethanol concentration and liquid-solid ratio, showed a significant influence (p ≤ 0.05) on the yield content during the extraction process. Regarding the p values of the quadratic term, it was noticed that only the quadratic term of the temperature was significant (p ≤ 0.05). According to the results of the study by Kashyap et al. [41], in which the optimization of ultrasound-assisted extraction of polyphenols from Meghalayan cherry fruit was investigated, it was found that the extraction yield of extract from Megalayan cherry fruit was highly influenced by solvent concentration, solvent-to-solid ratio and extraction time, as shown by the regression analysis. The effects of interactions between temperature, extraction time and ethanol concentration were studied in the work of Mrkonjić et al. [42], where it was found that the interaction between extraction time and temperature showed a significant effect, while the interaction between ethanol concentration and temperature showed significance (p ≤ 0.01) on yield content in the Thymus serpyllum extracts. The final equation that can predict the behavior of researched UAE of yield content is shown in Table 4. Plots (response surface) of the contents of yield in investigated vacuum dried sour cherries are presented in Figure 2a, where the influence of three investigated parameters can be seen. centration were studied in the work of Mrkonjić et al. [42], where it was found that the interaction between extraction time and temperature showed a significant effect, while the interaction between ethanol concentration and temperature showed significance (p ≤ 0.01) on yield content in the Thymus serpyllum extracts. The final equation that can predict the behavior of researched UAE of yield content is shown in Table 4. Plots (response surface) of the contents of yield in investigated vacuum dried sour cherries are presented in Figure 2a, where the influence of three investigated parameters can be seen.

Total Phenols Content (TPC)
The total phenolic content in the vacuum dried sour cherry extracts obtained in this research ranged from 1.16 to 1.96 g GAE/100 g ( Table 2). To compare with, similar results for total phenolic content in different types of fruit were obtained from the research of Kashyap et al. [41], where the extraction yield of Megalayan cherry fruit extract ranged from 93.20 to 123.91 mg GAE/g. In the study of UAE of total anthocyanins from Rubia sylvatica Nakai fruit [39], the experimental values of extraction yield of total phenolics varied from 28.52 to 35.23 mg GA/g. In this work, the highest total phenolic content was measured in the extract of vacuum dried sour cherries obtained at 80 °C, 40% and 30 mL/g liquid-solid ratio, while the lowest was obtained in extract that is observed at 80 °C, 80% and 10 mL/g liquid-solid ratio. Table 4 shows the p values of the linear, interaction and

Total Phenols Content (TPC)
The total phenolic content in the vacuum dried sour cherry extracts obtained in this research ranged from 1.16 to 1.96 g GAE/100 g ( Table 2). To compare with, similar results for total phenolic content in different types of fruit were obtained from the research of Kashyap et al. [41], where the extraction yield of Megalayan cherry fruit extract ranged from 93.20 to 123.91 mg GAE/g. In the study of UAE of total anthocyanins from Rubia sylvatica Nakai fruit [39], the experimental values of extraction yield of total phenolics varied from 28.52 to 35.23 mg GA/g. In this work, the highest total phenolic content was measured in the extract of vacuum dried sour cherries obtained at 80 • C, 40% and 30 mL/g liquid-solid ratio, while the lowest was obtained in extract that is observed at 80 • C, 80% and 10 mL/g liquid-solid ratio. Table 4 shows the p values of the linear, interaction and quadratic terms for the regression coefficients from the second-order polynomial model and predictive model equation with neglected, insignificant coefficients. Based on these results obtained for the total phenolic content, it was found that positive effects on the total phenolic content were observed for all the parameters studied (temperature, ethanol concentration and liquid-solid ratio). Based on the p values of linear term, it was also found that all three investigated parameters had a significant effect (p ≤ 0.05) on the total phenolic content during UAE of vacuum dried sour cherry samples. Moreover, corresponding to p values of the interaction term, it has been noticed that interactions between temperature and ethanol concentration as well as between ethanol concentration and liquid-solid ratio were significant (p ≤ 0.05). Moreover, interaction between temperature and liquid-solid ratio, the quadratic term of liquid-solid ratio as an investigated response, was significant (p ≤ 0.05). Table 4 also shows the predictive model equation from the second-order polynomial model for the content of TP obtained by UAE. According to the results showed by Kashyaop et al. [41], solvent-to-solid ratio, solvent concentration and amplitude significantly influenced the total phenolic content of the extract from Meghalayan cherries, as shown by the regression analysis. In the work of Chen et al. [39], investigating the UAE of Rubia sylvatica Nakai fruit, after screening significant variables using the Plackett-Burman design, the ethanol concentration, extraction temperature and pH were further optimized, and it was concluded that the coefficients of the linear term of these three variables were significant for the yields of total phenolic content. The linear effects of ultrasonic power, solvent-to-sample ratio and ethanol concentration, as well as the quadratic effects of ultrasonic power and ethanol concentration, showed significant effects on total phenolic content in the study, in which the simultaneous optimization of UAE for phenolic compound content and antioxidant activity of Lycium ruthenicum Murr. fruit was investigated [43]. In the work of Mrkonjić et al. [42], it was noticed that there was a very small difference between the total phenolic values when using 30% and 60% ethanol at the same extraction time and temperature, and regarding the obtained p values, the significance of the linear term of ethanol concentration on the total phenolic content could be neglected. Two other linear terms, temperature and extraction time, exhibited high significant effect on total phenolic content (p ≤ 0.01). Response surface plots were presented to better visualize the multiple influences of all independent variables (Figure 2b). A significant effect (p ≤ 0.05) of temperature, ethanol concentration and liquid-solid ratio was also observed in Figure 2b. A significant interaction between temperature and liquid-solid ratio as well as between ethanol concentration and liquid-solid ratio (p ≤ 0.05) was also observed in Figure 2b. In this figure, the higher TPC on higher values of liquid-solid ratio can also be noticed, which can be explained by the fact that higher liquid-solid ratio leads to an increase of concentration gradient between solid (plant material) and solvent (liquid) phases. This higher concentration difference of the solute increases the dissolution and diffusivity of the solute in the solvent, increasing the extraction process [35].

Total Flavonoids Content (TFC)
The total flavonoid content in the vacuum dried sour cherry extracts ranged from 0.60 and 0.87 g CE/100 g, presented in Table 2. The highest total flavonoid content was obtained in the extract observed at 60 • C, 40% and 20 mL/g, while the lowest content was noticed in the extract observed at 40 • C, 80%, and 10 mL/g. Chen et al. [43] investigated the UAE of Lycium ruthenicum Murr. fruit and they noticed that the range of experimentally obtained total flavonoid values in the obtained extracts were between 15.20 and 26.66 mg RE/g. Just as in the case of total phenolic content, according to the corresponding p values of all terms, the positive influences on total flavonoid content in vacuum dried sour cherry extracts were found.
Two independent variables, ethanol concentration and liquid-solid ratio, had a significant influence (p ≤ 0.05) on total flavonoid content during the extraction process, while temperature had no significant effect on flavonoid content during extraction. According to the p values obtained for the interaction term, it was found that the interaction between ethanol concentration and liquid-solid ratio was significant (p ≤ 0.05). At the same time, the interactions between temperature and ethanol concentration as well as between temperature and liquid-solid ratio were not significant. In contrast to the results for total phenolic content, according to the quadratic term it was found that the influences of none of the three independent variables were significant.
The linear effects of solvent-to-sample ratio, ethanol concentration and the quadratic effects of ultrasonic power exhibited a significant (p ≤ 0.05) positive effect on the total flavonoid content in the obtained extracts, while the quadratic effects of ethanol concentration showed a highly significant (p ≤ 0.01) negative effect on the total flavonoid content in the extracts in the work of Chen et al. [43]. The UAE of polyphenols from garlic was investigated by Ciric et al. [44] and it was concluded that the interactions between extraction time and methanol concentration, extraction temperature and methanol concentration and extraction temperature and solvent-solid ratio showed significant (p ≤ 0.01) effects on total flavonoid content. Table 4 shows the predictive model equation from the second-order polynomial model for the UAE of total flavonoid content. Response surface plots are shown in Figure 2c, where the influences of three investigated parameters on the total flavonoids content have been presented. In this figure can be noticed the similar influence of the liquid-solid ratio, as it was in the case of TPC.

Total Monomeric Anthocyanins Content (TMAC)
The content of total monomeric anthocyanins ranged from 9.80 to 30.17 mg C3G/100 g ( Table 2). The total monomeric anthocyanins highest content was noticed in the extract obtained at the highest temperature studied, 80 • C, 40% and 10 mL/g, while the lowest content was observed in the vacuum dried sour cherry extract obtained at 60 • C, 60% and 30 10 mL/g. In the study by Chen et al. [39], the results for extraction yields of total anthocyanins ranged from 17.32 to 22.13 mg Cy3G/g from Rubia sylvatica fruit. The storage stability of blood fruit anthocyanins, contained in conventional extract obtained with methanol and in an ultrasound-assisted extract obtained with methanol in different buffer solutions, was studied over 90 days. It was found that the degradation of anthocyanins was lower in acidic buffers than in near neutral and alkaline buffers [45]. According to the p values of the linear term presented in Table 4, it was observed that only the liquid-solid ratio, as one of the three parameters studied, had a significant influence (p ≤ 0.05) on the total monomeric anthocyanins content of the vacuum dried sour cherry samples.
From the same table it was also noticed that in both the linear and quadratic terms, none of the investigated parameters had an insignificant effect on the total content of monomeric anthocyanin in vacuum dried sour cherry samples, except for the quadratic term of liquid-solid ratio, which showed a significant influence (p ≤ 0.05) on the total content of monomeric anthocyanin in the extracts. Table 4 also shows the equation of the predictive model from second-order polynomial model for the UAE of total monomeric anthocyanin content. In the investigation by Chen et al. [39] on the UAE of total anthocyanins from Rubia sylvatica Nakai fruits, it was found that for the yield of total anthocyanins, the pH showed negative effects, while ultrasonic power, ethanol concentration, liquid-solid ratio, extraction temperature and extraction time showed positive effects. Response surface plots of the levels of total monomeric anthocyanins of the studied vacuum dried sour cherries are shown in Figure 2d, where the influence of three investigated parameters can be seen. In the case of TMAC, it was noticed that the liquid-solid ratio influenced differently from TPC and TFC, in term that the lower the values of liquid-solid ratio that were obtained, the higher the values of TMAC.

Antioxidant Activity (DPPH, FRAP and ABTS Assay)
The lowest antioxidant activity determined by DPPH assay in analyzed extracts was obtained in the extract observed at 40 • C, 80% and 10 mL/g ( Table 2). This was also the extract in which the lowest content of total flavonoids was noticed. The highest antioxidant activity obtained by DPPH assay was obtained in the extract in which the highest content of total phenolics was also obtained, at 80 • C, 40% and 30 mL/g. From the p values shown in Table 4, it was concluded that for the linear term, all three independent variables, temperature, ethanol concentration and liquid-solid ratio, had a significant influence (p ≤ 0.05) on the antioxidant activity obtained (DPPH assay) during extraction. Regarding the p values of the interaction term, it was not found to be a significant interaction for antioxidant activity obtained by DPPH assay, while none of the investigated variables had a significant influence in the quadratic term. The equation that predicted the behavior of the studied extracts, with respect to antioxidant activity obtained by DPPH assay, is shown in Table 4.
The antioxidant activity determined by FRAP assay ranged from 30.19 to 51.33 µM Fe 2+ /g ( Table 2). It was also observed that the lowest and the highest antioxidant activity obtained by FRAP assay were obtained in the same extracts in which the highest and the lowest antioxidant activity obtained by DPPH assay were obtained (40 • C, 80% and 10 mL/g and 80 • C, 40% and 30 mL/g, respectively). With regard to the p values of the linear term presented in Table 4, it was observed that all three investigated variables had a significant influence (p ≤ 0.05) on the antioxidant activity observed by FRAP assay of the vacuum dried sour cherry samples. It was concluded that in the case of interaction term, only the interaction between ethanol concentration and liquid-solid ratio was significant (p ≤ 0.05), while in the quadratic term, the influence of temperature was not significant. Liquid-solid ratio had a significant influence (p ≤ 0.05) on antioxidant activity observed by the FRAP assay of the vacuum dried sour cherry samples. Table 4 shows the equation of the predictive model from the second-order polynomial model for the UAE of antioxidant activity obtained by FRAP assay.
The antioxidant activity determined by ABTS assay ranged from 76.93 to 125.22 µM/g is presented in Table 2. The highest content was obtained in the extract observed at 80 • C, 80% and 10 mL/g. This was also the extract in which the highest content of total monomeric anthocyanins was observed. The lowest antioxidant activity detected by ABTS assay was obtained in the extract observed, at 80 • C, 80% and 30 mL/g. According to the results presented in Table 4, it was found that for the linear term, only one independent variable, ethanol concentration, had a significant influence (p ≤ 0.05) on the antioxidant activity obtained by ABTS assay during extraction. Regarding the p values of the interaction term, it was concluded that only the interaction between temperature and ethanol concentration had significant influence on antioxidant activity obtained by ABTS (p ≤ 0.05).
With regard to the p values of the quadratic term, it was found that the extraction time had a significant influence (p ≤ 0.05) on antioxidant activity obtained with ABTS. The equation that predicted the behavior of the studied UAE of antioxidant activity obtained by ABTS assay is shown in Table 4. Response surface plots of antioxidant activity (ABTS assay) in the investigated vacuum dried sour cherries are given in Figure 2c, where the influence of three investigated parameters can be seen. According to the results obtained by Mrkonjić et al. [42], the antioxidant activity of T. serpyllum extracts, obtained by DPPH, FRAP and ABTS assays, varied between 0.2431 and 0.2914 mMTE/g, 0.6726 and 0.935 7mM Fe 2+ /g and 0.4362 and 0.6482 mM TE/g, respectively. Moreover, according to the results of Chen et al. [43], FRAP values in the obtained extracts were found to be primarily affected by the solvent-to-sample ratio, and FRAP was also found to be higher at a high level of ultrasonic power and solvent-to-sample ratio. In the same work, it was noticed that with respect to ABTS assay, ABTS is also largely determined by solvent-to-sample ratio.
The response surface plots are shown in Figure 3a-c, where the correlation between investigated parameters and analyzed antioxidant activity (DPPH, FRAP and ABTS assays) of vacuum dried sour cherry extracts can be seen. In terms of all three investigated assays, in Figure 3a-c can be seen the influence of temperature, where the higher values of antioxidant activity were obtained at higher temperatures, which can be explained with the multiple influences of temperature on the process of mass-transfer, such as improved diffusion, improvement of solvent characteristics and degradation of the plant matrix [23].
sults of Chen et al. [43], FRAP values in the obtained extracts were found to be primarily affected by the solvent-to-sample ratio, and FRAP was also found to be higher at a high level of ultrasonic power and solvent-to-sample ratio. In the same work, it was noticed that with respect to ABTS assay, ABTS is also largely determined by solvent-to-sample ratio.
The response surface plots are shown in Figure 3a-c, where the correlation between investigated parameters and analyzed antioxidant activity (DPPH, FRAP and ABTS assays) of vacuum dried sour cherry extracts can be seen. In terms of all three investigated assays, in Figure 3a-c can be seen the influence of temperature, where the higher values of antioxidant activity were obtained at higher temperatures, which can be explained with the multiple influences of temperature on the process of mass-transfer, such as improved diffusion, improvement of solvent characteristics and degradation of the plant matrix [23].

Process Optimization and Experimental Verification
Ultrasound-assisted extraction parameters of vacuum dried sour cherry were optimized in order to maximize yields of the investigated chemical and biological properties. The optimized extraction conditions obtained in this research were the following: a maximum applied temperature of 80 °C, a minimum applied ethanol concentration of 40% and an 18.88 mL/g liquid-solid ratio, which presents approximately half of applied ratio values. The predicted values of the investigated responses, according to the optimized extraction conditions, were: 74.64% yield; 1.59 g GAE/100 g total phenolic content; 0.85 g CE/100 g total flavonoid content; 17.99 mg C3G/100 g total monomeric anthocyanins content; 54.90 µ M/g antioxidant activity (DPPH assay); 44.72 µ M Fe 2+ /g antioxidant activity (FRAP assay); 130.29 µ M/g antioxidant activity (ABTS assay). Determination of the predicted values and optimal conditions was obtained based on the function of desirability

Process Optimization and Experimental Verification
Ultrasound-assisted extraction parameters of vacuum dried sour cherry were optimized in order to maximize yields of the investigated chemical and biological properties. The optimized extraction conditions obtained in this research were the following: a maximum applied temperature of 80 • C, a minimum applied ethanol concentration of 40% and an 18.88 mL/g liquid-solid ratio, which presents approximately half of applied ratio values. The predicted values of the investigated responses, according to the optimized extraction conditions, were: 74.64% yield; 1.59 g GAE/100 g total phenolic content; 0.85 g CE/100 g total flavonoid content; 17.99 mg C3G/100 g total monomeric anthocyanins content; 54.90 µM/g antioxidant activity (DPPH assay); 44.72 µM Fe 2+ /g antioxidant activity (FRAP assay); 130.29 µM/g antioxidant activity (ABTS assay). Determination of the predicted values and optimal conditions was obtained based on the function of desirability that is in the range from zero to one. The desirability obtained in this research was 0.679. Furthermore, UAE was performed at the obtained optimized conditions of temperature, ethanol concentration and liquid-solid ratio, with the goal to validate the mathematical model predicted in this research. The results of experimentally observed responses were in good accordance with predicted values (74.72%, 1.56 g GAE/100 g, 0.86 g CE/100 g, 17.11 mg C3G/100 g, 53.84 µM/g (DPPH), 43.31 µM Fe 2+ /g and 131.87 µM/g (ABTS), respectively). This confirmed the successful optimization of the recovery of dried sour cherry polyphenolic antioxidants by UAE.

Conclusions
After the preliminary screening of the studied UAE factors (temperature, extraction time, ethanol concentration, ultrasonic power and liquid-solid ratio), it was found that three investigated parameters (temperature, liquid-solid ratio and ethanol concentration) showed the highest influence on total phenolic and monomeric anthocyanin content throughout the UAE of previously vacuum dried sour cherry samples. Therefore, mentioned parameters were selected in the main experimental design, while the analysis of responses was extended on yield content, total flavonoid content and antioxidant activity (DPPH, FRAP and ABTS assays). Face-centered experimental design was applied in order to optimize the UAE of vacuum dried sour cherry samples.
According to the obtained results, it was also found that significantly high values of the coefficients of multiple determinations (R 2 ) were obtained for all investigated responses, between 0.8391 (DPPH) and 0.9877 (TMAC), indicating a good fit between experimentally obtained values and values predicted by the model that was applied in this research. According to the CV values, a similar conclusion was drawn as the CV values varied from 1.68% (Y) to 10.04% (DPPH). In agreement, ANOVA results were also obtained, where a good fit between experimental and predicted values was observed.
Finally, it can be concluded that multi-response optimization of UAE can be successfully applied to obtain maximum bioactivity in vacuum dried sour cherry extracts. These results provide a concept for the possible application of the extracts obtained and their use in various forms in both scientific and industrial fields. The extracted compounds can be used directly in most industries, for example in the spirits industry, or as food additives, as in essential oils and high value-added compounds.