High-Performance Extraction Process of Anthocyanins from Jussara (Euterpe edulis) Using Deep Eutectic Solvents

New strategies for obtaining target bioactive compounds and natural pigments with the use of “green solvents” are consistently being developed, and deep eutectic solvents are (DES) a great alternative. This work established the significant variables and models for anthocyanin extraction, using DES and experimental design, of Euterpe edulis Mart. (jussara) fruit pulp, an endangered palm tree from the Brazilian Atlantic Forest. From a screening of seven initially tested DES, choline chloride/xylitol-based solvents had the best results with up to 42% increase in the total anthocyanin yield compared to methanolic extraction. Antioxidant assays also revealed a maximum antioxidant capacity of 198.93 mmol Trolox/100 g dry weight basis. The DES extract showed slower degradation to heat at 60◦ and 90 ◦C (2.5 times) and indoor constant light source (1.9 times) than methanolic extracts. The optimal extract also revealed slight inhibition of S. enterica and S. aureus growth in the


Introduction
Volatile Organic Solvents (VOS) are estimated to account for 60% of all industrial emissions and 30% of all volatile organic compounds emitted worldwide [1,2]. Although commonly used to extract organic compounds, such solvents cause environmental and social damage due to their limitations on reuse, toxicity, flammability, volatility, and presence of residues in the final product [3][4][5]. These risks have led to restrictive legislation such as Directive 2010/75/EU proposed by the European Parliament, which aims to limit the emission of certain VOS by industries.
Seeking to work around this problem, new technologies aiming at safer and more sustainable alternatives for VOS use have emerged, among them the so-called deep eutectic solvents (DES). These liquids have more desirable solvent characteristics such as negligible volatility, non-flammability, chemical, thermal, and electrochemical stability, and the ability to be synthesized from non-petroleum derived organic and/or inorganic ionic species [6]. Considering their unique properties and attractive price, DES have become a promising topic for academic research and industry for both separation processes and extraction of plant bioactive compounds [7,8].
Anthocyanins are a class of natural pigments with underexplored use in the food and cosmetics industry, which commonly uses artificial colorants. DES have been successfully tested for anthocyanin extraction, contributing to the "clean label" of the product, a strategy used by industry to attract increasingly demanding consumers seeking environmentally responsible consumption and healthy eating choices [7,9]. Potential health risks involving

Experimental Design: Response Surface Methodology
We evaluated the seven synthetized DES potentials for extracting anthocyanins by ultrasound-assisted extraction (Eco-sonics Ultronique, Campinas, Brazil), with fixed potency of 400 W for all experiments. An initial screening to choose 2 of the seven initial DES was conducted comparing anthocyanin content in the DES extracts and control (MeOH) with 30% co-solvent (acidified water), 1:15 sample:solvent ratio, and 3 extraction repetitions, during 5 min in the ultrasonic probe. Two of the seven initial DES were chosen (DES1 and DES2), based on anthocyanin yield, then two experimental designs were conducted: first, a fractional factorial experimental design 2 5-1 (FFED) to evaluate possible non-significant variables, and second, a central composite rotatable design (CCRD) with remaining significant variables. The FFED 2 5−1 with three central points and three levels for each independent variable was adopted, giving a total of 19 trials (Tables 2 and S1A). The independent variables evaluated for each extraction were the following: (i) eutectic solvent composition; (ii) solid-liquid ratio (R(S/L)), meaning the mass (g) of fruit pulp per mass of solution (g), (iii) number of extraction repetitions, (iv) time of extraction, and (v) co-solvent content (%) were chosen. The co-solvent (water) was acidified with HCl 3% to ensure the protonated state of anthocyanins [29]. Methanol with the same co-solvent proportion and parameters was used as a control.

Experimental Design: Response Surface Methodology
We evaluated the seven synthetized DES potentials for extracting anthocyanins by ultrasound-assisted extraction (Eco-sonics Ultronique, Campinas, Brazil), with fixed potency of 400 W for all experiments. An initial screening to choose 2 of the seven initial DES was conducted comparing anthocyanin content in the DES extracts and control (MeOH) with 30% co-solvent (acidified water), 1:15 sample:solvent ratio, and 3 extraction repetitions, during 5 min in the ultrasonic probe. Two of the seven initial DES were chosen (DES1 and DES2), based on anthocyanin yield, then two experimental designs were conducted: first, a fractional factorial experimental design 2 5-1 (FFED) to evaluate possible non-significant variables, and second, a central composite rotatable design (CCRD) with remaining significant variables. The FFED 2 5−1 with three central points and three levels for each independent variable was adopted, giving a total of 19 trials (Tables 2 and S1A). The independent variables evaluated for each extraction were the following: (i) eutectic solvent composition; (ii) solid-liquid ratio (R(S/L)), meaning the mass (g) of fruit pulp per mass of solution (g), (iii) number of extraction repetitions, (iv) time of extraction, and (v) co-solvent content (%) were chosen. The co-solvent (water) was acidified with HCl 3% to ensure the protonated state of anthocyanins [29]. Methanol with the same co-solvent proportion and parameters was used as a control.

Experimental Design: Response Surface Methodology
We evaluated the seven synthetized DES potentials for extracting anthocyanins by ultrasound-assisted extraction (Eco-sonics Ultronique, Campinas, Brazil), with fixed potency of 400 W for all experiments. An initial screening to choose 2 of the seven initial DES was conducted comparing anthocyanin content in the DES extracts and control (MeOH) with 30% co-solvent (acidified water), 1:15 sample:solvent ratio, and 3 extraction repetitions, during 5 min in the ultrasonic probe. Two of the seven initial DES were chosen (DES1 and DES2), based on anthocyanin yield, then two experimental designs were conducted: first, a fractional factorial experimental design 2 5-1 (FFED) to evaluate possible non-significant variables, and second, a central composite rotatable design (CCRD) with remaining significant variables. The FFED 2 5−1 with three central points and three levels for each independent variable was adopted, giving a total of 19 trials (Tables 2 and S1A). The independent variables evaluated for each extraction were the following: (i) eutectic solvent composition; (ii) solid-liquid ratio (R(S/L)), meaning the mass (g) of fruit pulp per mass of solution (g), (iii) number of extraction repetitions, (iv) time of extraction, and (v) co-solvent content (%) were chosen. The co-solvent (water) was acidified with HCl 3% to ensure the protonated state of anthocyanins [29]. Methanol with the same co-solvent proportion and parameters was used as a control.

Experimental Design: Response Surface Methodology
We evaluated the seven synthetized DES potentials for extracting anthocyanins by ultrasound-assisted extraction (Eco-sonics Ultronique, Campinas, Brazil), with fixed potency of 400 W for all experiments. An initial screening to choose 2 of the seven initial DES was conducted comparing anthocyanin content in the DES extracts and control (MeOH) with 30% co-solvent (acidified water), 1:15 sample:solvent ratio, and 3 extraction repetitions, during 5 min in the ultrasonic probe. Two of the seven initial DES were chosen (DES1 and DES2), based on anthocyanin yield, then two experimental designs were conducted: first, a fractional factorial experimental design 2 5-1 (FFED) to evaluate possible non-significant variables, and second, a central composite rotatable design (CCRD) with remaining significant variables. The FFED 2 5−1 with three central points and three levels for each independent variable was adopted, giving a total of 19 trials (Tables 2 and S1A). The independent variables evaluated for each extraction were the following: (i) eutectic solvent composition; (ii) solid-liquid ratio (R(S/L)), meaning the mass (g) of fruit pulp per mass of solution (g), (iii) number of extraction repetitions, (iv) time of extraction, and (v) co-solvent content (%) were chosen. The co-solvent (water) was acidified with HCl 3% to ensure the protonated state of anthocyanins [29]. Methanol with the same co-solvent proportion and parameters was used as a control.

Experimental Design: Response Surface Methodology
We evaluated the seven synthetized DES potentials for extracting anthocyanins by ultrasound-assisted extraction (Eco-sonics Ultronique, Campinas, Brazil), with fixed potency of 400 W for all experiments. An initial screening to choose 2 of the seven initial DES was conducted comparing anthocyanin content in the DES extracts and control (MeOH) with 30% co-solvent (acidified water), 1:15 sample:solvent ratio, and 3 extraction repetitions, during 5 min in the ultrasonic probe. Two of the seven initial DES were chosen (DES1 and DES2), based on anthocyanin yield, then two experimental designs were conducted: first, a fractional factorial experimental design 2 5-1 (FFED) to evaluate possible non-significant variables, and second, a central composite rotatable design (CCRD) with remaining significant variables. The FFED 2 5−1 with three central points and three levels for each independent variable was adopted, giving a total of 19 trials (Tables 2 and S1A). The independent variables evaluated for each extraction were the following: (i) eutectic solvent composition; (ii) solid-liquid ratio (R(S/L)), meaning the mass (g) of fruit pulp per mass of solution (g), (iii) number of extraction repetitions, (iv) time of extraction, and (v) co-solvent content (%) were chosen. The co-solvent (water) was acidified with HCl 3% to ensure the protonated state of anthocyanins [29]. Methanol with the same co-solvent proportion and parameters was used as a control.

Experimental Design: Response Surface Methodology
We evaluated the seven synthetized DES potentials for extracting anthocyanins by ultrasound-assisted extraction (Eco-sonics Ultronique, Campinas, Brazil), with fixed potency of 400 W for all experiments. An initial screening to choose 2 of the seven initial DES was conducted comparing anthocyanin content in the DES extracts and control (MeOH) with 30% co-solvent (acidified water), 1:15 sample:solvent ratio, and 3 extraction repetitions, during 5 min in the ultrasonic probe. Two of the seven initial DES were chosen (DES1 and DES2), based on anthocyanin yield, then two experimental designs were conducted: first, a fractional factorial experimental design 2 5-1 (FFED) to evaluate possible non-significant variables, and second, a central composite rotatable design (CCRD) with remaining significant variables. The FFED 2 5−1 with three central points and three levels for each independent variable was adopted, giving a total of 19 trials (Tables 2 and S1A). The independent variables evaluated for each extraction were the following: (i) eutectic solvent composition; (ii) solid-liquid ratio (R(S/L)), meaning the mass (g) of fruit pulp per mass of solution (g), (iii) number of extraction repetitions, (iv) time of extraction, and (v) co-solvent content (%) were chosen. The co-solvent (water) was acidified with HCl 3% to ensure the protonated state of anthocyanins [29]. Methanol with the same co-solvent proportion and parameters was used as a control.

Experimental Design: Response Surface Methodology
We evaluated the seven synthetized DES potentials for extracting anthocyanins by ultrasound-assisted extraction (Eco-sonics Ultronique, Campinas, Brazil), with fixed potency of 400 W for all experiments. An initial screening to choose 2 of the seven initial DES was conducted comparing anthocyanin content in the DES extracts and control (MeOH) with 30% co-solvent (acidified water), 1:15 sample:solvent ratio, and 3 extraction repetitions, during 5 min in the ultrasonic probe. Two of the seven initial DES were chosen (DES1 and DES2), based on anthocyanin yield, then two experimental designs were conducted: first, a fractional factorial experimental design 2 5-1 (FFED) to evaluate possible non-significant variables, and second, a central composite rotatable design (CCRD) with remaining significant variables. The FFED 2 5−1 with three central points and three levels for each independent variable was adopted, giving a total of 19 trials (Table 2 and S1A). The independent variables evaluated for each extraction were the following: (i) eutectic solvent composition; (ii) solid-liquid ratio (R (S/L) ), meaning the mass (g) of fruit pulp per mass of solution (g), (iii) number of extraction repetitions, (iv) time of extraction, and (v) co-solvent content (%) were chosen. The co-solvent (water) was acidified with HCl 3% to ensure the protonated state of anthocyanins [29]. Methanol with the same co-solvent proportion and parameters was used as a control. After the extractions, the extracts were filtered with filter paper (Unifil qualitative 80 g/m 2 ) and under a vacuum, then concentrated using a rotary evaporator (<37 • C), removing the water content of the co-solvent for quantification. The results were evaluated using ANOVA followed by Dunnett's test considering a confidence level of 95%, comparing the DES results to the control (MeOH extractions) in terms of total anthocyanin level (TAL). The response was expressed in µg anthocyanins /mL. An estimate of the main effect was obtained by evaluating the difference in the total anthocyanin yield caused by a change from the low (−1) to the high (+1) levels of the corresponding variable. The fractional factorial experimental design's focus was to discover the significant variables and their main interactions with the total yield of extracted anthocyanins. Therefore, after analyzing the data, the variables that had a significative effect on the total yield of the process were fixed and the others removed (including the removal of DES2 from further trials) for an additional assay CCRD 2 2 plus axial and central composites, with three replicates at the central point, giving a total of 11 additional trials. Table 3 shows the real values (coded values in Table S1B) and levels used at CCRD 2 2 . After analyzing the CCRD results, the best conditions for the anthocyanin extraction with the best DES was determined, and the model was validated in triplicate to be used in further trials. Table 2. Independent variable values of fractional factorial experimental design (2 5-1 ) with three central points for DES1 and DES2 were used as independent variables as well as R (S/L) , Repetitions, time, and R (CoSolvt%) .  The main response evaluated for optimization of the process was based on the total yield of cyanidin 3-glucoside (C3G), cyanidin 3-rutinoside (C3R), and total anthocyanin level (TAL), considering only these two major anthocyanins, attempting to obtain models for those responses as well. Statistica 14.0 software was used to analyze the results and plot the response surfaces.

Identification and Quantification of Anthocyanins
For all the samples, anthocyanin separation was carried out on a 250 × 4.6 i.d. mm, the latter proportion being maintained for a further 15 min, at a flow rate of 0.9 mL/min and a column temperature set at 28 • C. The chromatograms were processed at 280 and 520 nm, and the spectra were obtained between 250 and 800 nm [19,30]. The anthocyanins were identified based on the combined information provided by the elution order in the reversed phase column, co-chromatography with standards, UV-VIS, and mass spectra compared to the literature data. Anthocyanins were quantified by HPLC as cyanidin 3-glucoside and cyanidin 3-rutinoside, using an external calibration curve for both with a minimum of six concentration levels. All analyses were performed in triplicate.

Antioxidant Activity
In order to compare the antioxidant activity of extracts obtained by DES in CCRD trials, ABTS [2,2 -azinobis-(3-ethylbenzothiazoline-6-sulfonate)] assay was performed as described by [31]. Antioxidant activity (AA) of the extracts was expressed in µmol Trolox/g equivalents. Results were then plotted in CCRD again for the same effects, using antioxidant activity as the dependent variable and analyzed by ANOVA with a 90% confidence level.
The oxygen radical absorbance capacity (ORAC) method was also performed as previously described by [32], which consists of measuring the decrease in the fluorescence of a protein as a result of the loss of its conformation when it suffers oxidative damage caused by a source of peroxyl radicals. The measurements were taken in triplicate. The ORAC values, expressed as µM Trolox equivalents (µM TE), were calculated by applying the following Equation (1): where C Trolox is the concentration (µM) of Trolox (20 µM), k is the sample dilution factor, and AUC is the area below the fluorescence decay curve of the sample, blank, and Trolox, respectively. Results were then plotted in CCRD again for the same effects, using antioxidant activity as the dependent variable and analyzed by ANOVA with a 90% confidence level.

Anthocyanin Thermal Stability
The determination of the anthocyanin's thermal stability was performed for both best DES and methanolic extracts. The stability at 60 • C and 90 • C was performed as described by De Rosso and Mercadante [33], with modifications. After the extraction using the validated condition from the CCRD, the best DES and methanolic extracts were prepared and pH adjusted to 3.5 for all samples to keep anthocyanins stable (protonated). A food preservative (10 mg potassium sorbate) was added to both extracts to avoid spoilage. Then, each solution absorbance was adjusted to 0.8 through dilution, measured at the maximum absorption wavelength in the visible region for anthocyanins (520 nm) to standardize the assays, ensuring that Lambert-Beer law is considered and the citrate-phosphate buffer was utilized since small variations of pH affects anthocyanin concentration. The solutions were allowed to rest for 3 h to attain equilibrium among the different forms of anthocyanin. The solutions were then distributed in 10 mL centrifuge tubes with a total sample volume of 3 mL.
All tubes were placed in a water bath and heated at 60 • C or 90 • C, removing triplicates of each extract (DES or MeOH) at random times. The samples were immediately cooled in an ice bath for 10 min after removal. After this time, they were placed at room temperature for 50 min before reading the absorbance in a spectrophotometer (520 nm). When the absorbance pointed to half the original value (0.4), the assay was finished and results were plotted, and the emerging curve behavior was analyzed.

Anthocyanin Photostability
The anthocyanin's photostability was performed for both DES and methanolic extracts in the presence and absence of light. The extracts were prepared as described in Section 2.7. The solutions were then distributed in 15 mL transparent tubes with a total sample volume of 10 mL.
The tubes were placed both in constant light and constant dark conditions, and samples of each triplicate of each extract (DES or MeOH) were removed at random times over the course of 25 days until absorbance reached half of its initial value (0.4). Samples had their absorbance were read in a spectrophotometer (520 nm) at these intervals. Results were plotted, and the emerging curve behavior was analyzed.

Disc-Diffusion Antimicrobial Sensitivity Test
The growth method was performed as follows: (1) three colonies, well isolated, of the same morphological type of Salmonella enterica subsp. (ATCC 13076) and Staphylococcus aureus (ATTCC 19095) were selected from the agar plate with Müeller-Hilton culture medium. Microorganisms were transferred to a tube containing 4-5 mL of culture medium.
(2) The culture was incubated in broth at 35 • C until reaching or exceeding the standard McFarland 0.5 solution's turbidity measured in a spectrophotometer at 600 nm. (3) The turbidity of the growing culture was adjusted with a sterile saline solution in order to obtain optical turbidity comparable to that of the 0.5 McFarland standard solution. A sterile cotton swab was then dipped into the adjusted suspension up to 15 min after changing the inoculum suspension's turbidity. The dry surface of the Müeller-Hinton agar plate was inoculated by rubbing the swab across the agar's sterile surface ensuring uniform distribution of the inoculum. A set of antimicrobial discs was placed on the surface of a seeded agar plate, applying samples of the best DES extract, SF (saline), DES (deep eutectic solvent only), and positive control (vancomycin for S. aureus and meropenem for S. enterica), as well as samples of the optimized extract. The plates were inverted and placed in an oven at 35 • C for up to 15 min after applying the discs. After 16-18 h of incubation, each plate was examined. The diameters of the total inhibition halos were measured by the diameter (mm) of the disc.

Screening Process
A chromatographic profile of the methanolic jussara extract from DES screening analyzed the two major anthocyanins: C3G, with a retention time of 13.39 min, and C3R, with a retention time of 13.99 min ( Figure S1 of Supporting Information). Quantification of DES and methanolic extract's peak areas to C3G and C3R, respectively, are shown in Figure 1. A confidence level of 95% was used for all analyses.  Figure S1 of Supporting Information). Quantification of DES and methanolic extract's peak areas to C3G and C3R, respectively, are shown in Figure 1. A confidence level of 95% was used for all analyses.
The results of C3G content in DES extracts were not significant (compared to MeOH; Mean: 114.79 ± 3.56 mg/100 g dry matter, ANOVA followed by Dunnett's test) for any of the extracts tested. The results of C3R in the extracts, on the other hand, were significant for Ch-Leu (785.16 ±99.2 mg/100 g dry matter, lower) and Lp-But (1092.41 ± 88.33 mg/100 g dry matter, lower), Ch-Xyl (1343.92 ±172.42 mg/100 g dry matter, higher), and Lp-Gly (1224.55 ± 80.20 mg/100 g dry matter, higher) in comparison to MeOH extract (935.48 ± 22.28 mg/100 g dry matter). C3R is the main anthocyanin of the pulp, and its yield significantly influenced the total anthocyanin level (TAL) of extracts when compared to methanolic extract. After the sum of both anthocyanins, only the TAL of Ch-Xyl (1490.65 ± 198.27 mg/100 g dry matter) and Lp-Gly (1359.06 ± 92.92 mg/100 g dry matter) showed significant differences compared to the MeOH extract (1050.28 ± 29.32 mg/100 g dry matter). As our samples had 88.7% humidity, the TAL content in the methanolic extract can be expressed as The results of C3G content in DES extracts were not significant (compared to MeOH; Mean: 114.79 ± 3.56 mg/100 g dry matter, ANOVA followed by Dunnett's test) for any of the extracts tested. The results of C3R in the extracts, on the other hand, were significant for Ch-Leu (785.16 ±99.2 mg/100 g dry matter, lower) and Lp-But (1092.41 ± 88.33 mg/100 g dry matter, lower), Ch-Xyl (1343.92 ±172.42 mg/100 g dry matter, higher), and Lp-Gly (1224.55 ± 80.20 mg/100 g dry matter, higher) in comparison to MeOH extract (935.48 ± 22.28 mg/100 g dry matter). C3R is the main anthocyanin of the pulp, and its yield significantly influenced the total anthocyanin level (TAL) of extracts when compared to methanolic extract.
After the sum of both anthocyanins, only the TAL of Ch-Xyl (1490.65 ± 198.27 mg/100 g dry matter) and Lp-Gly (1359.06 ± 92.92 mg/100 g dry matter) showed significant differences compared to the MeOH extract (1050.28 ± 29.32 mg/100 g dry matter). As our samples had 88.7% humidity, the TAL content in the methanolic extract can be expressed as 118.68 ± 3.3 mg/100 g fresh matter, being compatible with the MeOH extracts from the literature. One study [21] obtained 14.84 to 409.85 mg/100 g fresh matter with acidified methanolic extract, using samples gathered in different months. More recently [34], a change was observed in the total monomeric anthocyanins during E. edulis ripening. That study found that the anthocyanin content was very positively correlated with the ripening process, with a maximum of 634.26 mg TAL/100 g fresh matter. A progressive increase in the concentration of anthocyanins is believed to occur during ripening, and changes in the color of the fruits are easily observable [35]. Increased anthocyanin levels during ripening have also been previously noted for açai (Euterpe oleracea) [35,36].
The medium absolute values of Ch-Xyl anthocyanin extracts were 42% higher than methanolic extract and Lp-Gly was 29% higher. These results clearly lead us to choose Ch-Xyl and Lp-Gly as the best green solvents for further experimental optimization of the extraction process (DES1 and DES2 of FFED, respectively).

Fractional Factorial Experimental Design (FFED)
After choosing the DES with a higher capacity to extract total anthocyanins from E. edulis lyophilized pulp, the optimization of different process conditions was conducted through a fractional factorial experimental design with three central points 2 5-1 , following the standardized effects described in the methodology section. The extraction process in different experimental conditions of the 2 5-1 and the yields are shown in Table 4. Table 4. Experimental values for C3G, C3R, and TAL obtained from fractional experimental design (2 5-1 ) with three central points for Ch-Xyl and Lp-Gly used as independent variables as well as R (S/L) , extraction repetitions, time, and R (CoSolvt%) .  The experimental data shows that total anthocyanins extracted in all factorial design trials ranged from 829.29 to 1394.97 mg/100 g dry matter, indicating that the selected variables and levels could change the process efficiency up to 60%. As mentioned before, the proportion of DES is a qualitative variable and the levels studied were DES1 (−1) and DES2 (+1), and for the central point (0), we prepared a mixture of both DES. The negative effect of the proportion of DES, obtained from the FFED means that when we use DES2 instead of DES1, i.e., when we change the level from −1 to +1, the extraction of anthocyanins is less efficient. The fact that the DES2 variable has a negative effect implies that DES1 (Ch-Xyl) is more efficient in extracting the target bioactive compounds. Co-solvent usage had a positive influence, moving towards level +1 (more co-solvent in the total solvent used) and repetition (more repetitions). The solid/liquid ratio and time had no significant effects in this factional design. Considering these results, only two variables were carried out for the next step: co-solvent and repetitions.
Although the DES choice has had a significant effect, this variable could be fixed with Ch-Xyl for the CCRD (2 2 ) due to the negative mean effect of the Lp-Gly (DES2) solvent as discussed above. Chlorine chloride and xylitol were present in other tailor-made DES [37,38] for phenolic compounds and anthocyanins [39], where they are described as showing a high affinity for phenolic structure. Analysis obtained with the use of nuclear magnetic resonance spectroscopy by [40] shows that a Ch-Xyl-Water mixture shows great affinity for phenolic compounds compared to many other combinations of DES, including the Ch-Gly tested here. Additionally, Ch-Xyl and water exhibit a more cohesive structure when dissolving quercetin (a phenolic compound) than a citric acid/cholinechloride/water mixture [41]. However, in analyzing extraction of anthocyanins from mulberries, Ch-Gly gave a slightly better anthocyanin extraction than Ch-Xyl and a better extraction than mixtures of chloride-mannitol, chloride-fructose and choline chlorideglucose [39], indicating that food matrix and the type of anthocyanins contained can be important factors when tailoring DES for extraction of these compounds.

Central Composite Rotatable Design (CCRD)
As the literature indicates that using more than 50% of the water in the mixture can cause problems in DES stability [41,42], the co-solvent effect levels in CCRD were defined as 16% to 49.7%. For repetitions, previous experiments observed that beyond the fourth repetition, the solvent was almost colorless after extraction. A hypothesis was formulated then, assuming that after a hypothetical optimum of repetitions (four), more repetitions could generate TAL loss by excess manipulation. To verify this, a quick assay was done, with fixed variables (R(S/L) 1:15; Time: 2 min, co-solvent: 40%) with different repetitions (three, five, and seven times). Results are shown in Figure S2 of the Supporting Information, confirming a significant loss of efficiency in seven repetitions, compared to five repetitions (ANOVA).
After set repetition effect levels (2, 5 and 7 times), Table 5 shows the total yield obtained for the extraction process in different experimental conditions of the 2 2 with central and axial points, for C3G, C3R, and TAL content from the E. edulis used as a response to the construction of the predictive model. Table 5. Experimental values for C3G, C3R, and TAL extracted by CH-Xyl in the CCRD assays, and relative deviation (%) among experimental and predictive values of TAL. The main effects and interactions were estimated for TAL are described in Equation (2).

Assay C3G (µg/mL) C3R (µg/mL) TAL (µg/mL) Relative Deviation (%)
where X 1 is co-solvent participation (%) and X 2 is the number of repetitions. ANOVA was used to evaluate the data's adequacy to verify the possibility of obtaining the model considering the different responses; additionally, the R 2 value provided a confidence measure of the model. In this study, an R 2 value of 0.6969 was obtained. The F value 94.9 (F value tabulated: 4.103) for TAL exceeded the 95% confidence level. The TAL response attended the criteria to generate the model and, thus, can be considered predictive; therefore, the response surface was plotted in Figure 2A.  As observed in the model, repetitions largely influenced the TAL yield of extracts. On the other hand, next to maximum levels, the efficiency is lost; this is attributed to excess manipulation of extracts. The co-solvent had a mild influence on the model, as observed in the Pareto chart of effects. The negative influence when the water participation (%) rises followed other studies, in which it is observed that it can disrupt the DES structure and thus reduce the DES effect in the extraction [25,43]. The best results were obtained around the central point level, with four repetitions and 30% co-solvent participation. The validation assays, performed in triplicate, were reproduced applying the best conditions. The model was effectively validated (mean of TAL was 2350.61 ± 53.05 µg/mL) with low relative deviation (3.10%), evidencing its high confidence in predicting data. These results for co-solvent participation are compatible with others in the literature for DES extraction for phenolic compounds that report best yields at 25-30% water usage [28,44]. Choline cations are highly capable of coordinating the surrounding environment: strong hydrogen bondingmediated correlations between the hydroxyl group and water or chloride are believed to be formed. In addition, it is possible that the ammonium group drives the formation of a solvating environment, with water, chloride, and hydroxyl moieties approaching it, between the methyl groups [45]. The total anthocyanin content extracted by these assays varied from 1215.51 to 2318.48 µg/mL (assays 7 and 11, respectively), corresponding to 117.301 100 g −1 fresh weight and 223.734 100 g −1 fresh weight, respectively. Compared to the literature, a careful approach is necessary since the total anthocyanins contained in fruits can vary greatly due to edaphoclimatic variables and maturation stages [46]. The results obtained by [21] also using the experimental design ranged from 25.54 mg 100 g −1 fresh weight (methanol) up to 418.52 mg 100 g −1 fresh weight (Methanol/HCl 1.5 M). The study used an ultrasonic bath and the best condition with a 1:50 solvent ratio, which is considerably higher than that used in our study (1:15). In a study using the same solvent ratio [47] and 70% ethanolic solvent in ultrasonic bath, with 360 W potency, similar results were found: 190.92 ± 13.48 mg 100 g −1 fresh weight to 282.64 ± 11.07 mg 100 g −1 fresh weight. When authors tested a higher potency in the bath (900 W), a reduction in the anthocyanin content was observed.

Antioxidant Activity
After quantifying anthocyanin content in Ch-Xyl extracts for CCRD, the antioxidant activity of the samples was measured by ABTS and ORAC assays. Another CCRD assay was performed, using ABTS and ORAC assays as a response instead of TAL for the very same effects. Results are also shown in Figure 2B,C. Coded values can be found in Table  S2A of Supporting Information.
In the ABTS assay, ANOVA resulted in an R 2 value of 0.88 and F value of 20.5, for a 95% confidence level; therefore, it is able to generate a model with a predictive equation (Equation (3)).
Antioxidant activity ( mmol Trolox 100 /g) = 189.73 − 36.6 (X 1 ) 2 + 46.5( where X 1 is co-solvent participation (%), and X 2 is the number of repetitions. The graphical model and the Pareto chart of effects are shown in Figure 2B. For the ORAC assay (Table S2B of Supporting Information), ANOVA resulted in an R 2 value of 0.81 for a 95% confidence level and is therefore able to generate a model with a predictive Equation (4): Antioxidant activity ( mmol Trolox 100 /g) = 159.32 − 31.3 (X 1 ) 2 + 38.0(X 2 ) − 27.9 (X 2 ) 2 (4) where X 1 is co-solvent participation (%), and X 2 is the number of repetitions. The precision and exactness of the model were validated by comparing the experimental results with the theoretical data predicted by Equations (2)-(4). As observed, response surfaces of the same extracts show slightly different performances for anthocyanin content and antioxidant activity. It is known that jussara pulp contains multiple hydrophilic compounds such as phenolic compounds [19,21,48,49], and even anthocyanin degradation products, which possess antioxidant properties [50]. Such compounds can explain the dissonant behavior and different maximum values of antioxidant activity when reacting to ABTS and ORAC assays. The results had values between 69.92 and 198.93 for ABTS and between 63.09 and 196.436 mmol Trolox/100 g dwb (dry weight base) for ORAC. Values between 67.7 and 154.4 mmol Trolox/100 g dwb using TEAC (Trolox equivalent antioxidant capacity) and ORAC assays were found for jussara pulp [51]. While analyzing the ripening of the fruit, values between 108.8 (immature) and 207.1 (complete maturation) mmol Trolox/100 g dwb using the ORAC assay were found [52].Considering the evaluated dependent variables in CCRD (TAL and antioxidant activity by the methods ABTS and ORAC), the maximization of anthocyanin content and antioxidant activity are around the central point (applying around four repetitions and 30% co-solvent participation). This way, as presented for the TAL, the validation assays for antioxidant activity were also performed in triplicate for ABTS and ORAC. The result mean of the antioxidant activity was 198.82 ± 10.97 and 155.71 ± 21.77 (mmol Trolox 100 g −1 dwb), representing a relative deviation from the predicted values of 5.06% and 10.51% for ABTS and ORAC, respectively. This confirms that the models were validated with low relative deviation, evidencing their statistical significance and predictiveness.

Determination of Anthocyanins Thermal Stability
The changes in anthocyanin concentration from jussara extracts caused by the temperature increase with the time was expected, following a biphasic behavior that was best fitted by a biexponential equation (Equation (5)), considering both evaluated temperatures (60 and 90 • C). Figure 3 shows the degradation kinetics of total anthocyanins extracted with MeOH and with Ch-Xyl, and the calculated kinetic parameters obtained by fitting the data from Equation (5) can be found in Table S3 of Supporting Information. The biexponential decay behavior, verified for the initial concentrations of anthocyanins, showed the same fast and slow lifetimes at each temperature, thus confirming the occurrence of a shared overall degradation mechanism.
yt and y∞ are anthocyanin absorbance values at real time and infinite time, respectively. A1 and A2 are the pre-exponential factors, whereas γ1 and γ2 are the observed rate constants for the fast and slow decays, respectively.
The degradation rate constants (Kd) at 60 °C and 90 °C from anthocyanins present in The biexponential decay behavior, verified for the initial concentrations of anthocyanins, showed the same fast and slow lifetimes at each temperature, thus confirming the occurrence of a shared overall degradation mechanism.
y t and y ∞ are anthocyanin absorbance values at real time and infinite time, respectively. A 1 and A 2 are the pre-exponential factors, whereas γ 1 and γ 2 are the observed rate constants for the fast and slow decays, respectively.
The degradation rate constants (K d ) at 60 • C and 90 • C from anthocyanins present in Ch-Xyl extract are shown in Table 6. The K d values were estimated from the slope of the semi-natural logarithmic plot of residual activity vs. time, assuming classical Arrhenius linear behavior and using the equation between both temperatures. Similarly, the half-life (t 1/2 ) values were calculated using the degradation rate constants, also presented in Table 4. A comparison of the t 1/2 values of both extracts at different temperatures showed that the results of anthocyanins extracted with Ch-Xyl are more stable than those obtained with MeOH, with half-life 2.5 times greater, in both temperatures.
The activation energy (E a ) for the Ch-Xyl extract (115.35 kJ/mol) was higher compared to the MeOH extracts (53.94 kJ/mol), confirming that the Ch-Xyl utilization provides more thermally stable anthocyanins.
The literature shows that in several studies, the thermal degradation of bioactive compounds, such as enzymes and bioactive peptides, as well as other bioproducts, the kinetics followed a linear adjusted model. Meanwhile, the number of articles evaluating anthocyanins' stability using DES is limited, with different DES combinations and food matrices. Anthocyanins extracted from Catharanthus roseus with lactic acid-glucose (5:1) showed better stability in this DES than in acidified ethanol, about three times at 60 • C [9], although the curve showed a first-order model. Previous works from our research group showed that the thermal stability of carotenoids from tomatoes extracted using an ultrasound-assisted approach and 1-Butyl-3-methyl-imidazolium-chloride [BMIM][Cl] followed a first-order model [53]. However, they also showed the biexponential behavior in two studies regarding carotenoid extraction using distinct ionic liquids of orange peel [54] as well as from pupunha (Bactris gasipaes) fruit [55]. The higher stabilization ability of Ch-Xyl for anthocyanins may be correlated with the interactions between cyanidin and these molecules. The two components of Ch-Xyl may have intermolecular interactions, mainly hydrogen bonding with the carboxyl and hydroxyl groups of anthocyanins, as was observed for quercetin in Ch-Xyl [40]. This interaction is expected to decrease the movement of solute molecules, reducing its contact time with oxygen at the interface of anthocyanins and air, and consequently reducing oxidative degradation, which is the major degradation mechanism. In a study with anthocyanins extracted from black carrots using choline chloride/citric acid DES [56], the samples also exhibited improved stability.

Determination of Anthocyanins Photostability
Light is reported to have two different effects on anthocyanins: it accelerates the rate of the thermal reactions and the formation of new molecules, including chalcones [57]. Anthocyanins from jussara pulp changed with the time both in light and dark conditions and followed a first order curve and were best fitted by a linear equation (Equation (6)), considering both situations. Figure 4 and Table 7 present the degradation kinetics of total anthocyanins extracted with MeOH and with Ch-Xyl.  When comparing t1/2 values of Ch-Xyl and MeOH, the latter perished significantly (p < 0.05, ANOVA) faster-about 1.9 times in dark and 1.5 in light-which indicates that Ch-Xyl extract could have a better storage time than MeOH extract, especially in dark conditions. Anthocyanins from blueberry extract exposed to natural and UV light had the same pattern of degradation after 25 days [58], while they retained about 60% of their initial anthocyanin count in dark conditions and lower than 20% in light conditions. This study obtained 78% (MeOH) and 85% (Ch-Xyl) retention of initial absorbance after 25 days in dark and 43% (MeOH) and 63% (Ch-Xyl) retention of absorbance after 25 days, for jussara extracts exposed to light conditions. Degradation of anthocyanins of purple sweet potato analyzed for 15 days at 25 °C retained more than 80% anthocyanin content both under light and dark conditions [59]. The same was true for anthocyanins of Hibiscus sabdariffa for 10 days [60], and similar to results in this study, except for MeOH extract under light exposure, which retained only 63% of its initial absorbance after 15 days.  From a semi-natural logarithmic plot of the absorbance vs. time, the activation rate constants (K d ) were calculated, and the half-lives were estimated using Equation (6). The halflife (t 1/2 ) is defined as the time taken for the absorbance to reach 50% of the initial measurement.
When comparing t 1/2 values of Ch-Xyl and MeOH, the latter perished significantly (p < 0.05, ANOVA) faster-about 1.9 times in dark and 1.5 in light-which indicates that Ch-Xyl extract could have a better storage time than MeOH extract, especially in dark conditions. Anthocyanins from blueberry extract exposed to natural and UV light had the same pattern of degradation after 25 days [58], while they retained about 60% of their initial anthocyanin count in dark conditions and lower than 20% in light conditions. This study obtained 78% (MeOH) and 85% (Ch-Xyl) retention of initial absorbance after 25 days in dark and 43% (MeOH) and 63% (Ch-Xyl) retention of absorbance after 25 days, for jussara extracts exposed to light conditions. Degradation of anthocyanins of purple sweet potato analyzed for 15 days at 25 • C retained more than 80% anthocyanin content both under light and dark conditions [59]. The same was true for anthocyanins of Hibiscus sabdariffa for 10 days [60], and similar to results in this study, except for MeOH extract under light exposure, which retained only 63% of its initial absorbance after 15 days.
The strong hydrogen bond formed between Ch-Xyl and anthocyanins can once again explain the slower reaction rate of these extracts. This is pointed to as a possible strategy to deal with anthocyanins' relative instability compared to artificial colorants [39]. Evaluation of the thermal and light stability of anthocyanins from different sources is essential in order to determine their properties and hence apply those pigments, especially the most stable ones, for scaling-up processes and industrial applications. The bioaccessibility of anthocyanins during digestion, usually reported as very low [61], should be further evaluated when used in conjunction with Ch-Xyl solvents, evaluating their stabilization power in this environment to better understand the structure-stability relationships of the anthocyanin/extraction solvent.

Antimicrobial Activity
In order to evaluate the antimicrobial activity of extracts, a standard disk diffusion test for microbial sensitivity was conducted, and their inhibitory halo was measured ( Figure 5 and The strong hydrogen bond formed between Ch-Xyl and anthocyanins can once again explain the slower reaction rate of these extracts. This is pointed to as a possible strategy to deal with anthocyanins' relative instability compared to artificial colorants [39]. Evaluation of the thermal and light stability of anthocyanins from different sources is essential in order to determine their properties and hence apply those pigments, especially the most stable ones, for scaling-up processes and industrial applications. The bioaccessibility of anthocyanins during digestion, usually reported as very low [61], should be further evaluated when used in conjunction with Ch-Xyl solvents, evaluating their stabilization power in this environment to better understand the structure-stability relationships of the anthocyanin/extraction solvent.

Antimicrobial Activity
In order to evaluate the antimicrobial activity of extracts, a standard disk diffusion test for microbial sensitivity was conducted, and their inhibitory halo was measured (Figure 5 and Table S4 of Supporting Material). The bacterial pathogens, Salmonella enterica and Staphylococcus aureus, were tested, and the extract displayed moderate inhibition for both (8.6 mm and 7.6 mm halo respectively). S. enterica is one of the most frequent bacterial food-borne pathogens in humans. Salmonella infections range from gastrointestinal infections that are accompanied by inflammation of intestinal epithelia, diarrhea, and vomiting, to typhoid fever, a life-threatening systemic infection [62]. S. aureus is indeed found in the skin, hair, and nostrils of warm-blooded animals. This bacterium is a significant cause of nosocomial infections, as well as community-acquired diseases. The extract resistance to this microorganism is a good indicator for stability and avoidance of microbial growth during storage. The results for S. aureus are according to [18], where Euterpe edulis crude extract displayed slightly positive inhibition against this microorganism. This can be attributed to the antimicrobial activity of anthocyanins and other polyphenols found in the extract, which is usually an The bacterial pathogens, Salmonella enterica and Staphylococcus aureus, were tested, and the extract displayed moderate inhibition for both (8.6 mm and 7.6 mm halo respectively). S. enterica is one of the most frequent bacterial food-borne pathogens in humans. Salmonella infections range from gastrointestinal infections that are accompanied by inflammation of intestinal epithelia, diarrhea, and vomiting, to typhoid fever, a life-threatening systemic infection [62]. S. aureus is indeed found in the skin, hair, and nostrils of warmblooded animals. This bacterium is a significant cause of nosocomial infections, as well as community-acquired diseases. The extract resistance to this microorganism is a good indicator for stability and avoidance of microbial growth during storage. The results for S. aureus are according to [18], where Euterpe edulis crude extract displayed slightly positive inhibition against this microorganism. This can be attributed to the antimicrobial activity of anthocyanins and other polyphenols found in the extract, which is usually an antimicrobial agent in plants, and shows that DES extraction did not remove significant antimicrobial components from the crude extract [63]. Coupled with the relevant information of the stability of the extracts, further studies are necessary to determine the minimum inhibitory concentration (MIC) of the extract and to determine the ideal storage concentrations.
Together, these experiments build an additional step in the knowledge of how efficient and stable anthocyanin-rich extracts obtained with DES can be while keeping its antioxidant and antimicrobial properties. Further experiments should be conducted in order to investigate practical applications in a number of cosmetical and food products regarding the stability and interactions in these systems. Reducing the use of VOS is a worldwide effort, and general industry faces additional appeals for "green label" procedures by its consumers. Alternative green solvents fit this necessity well, and natural colorants obtained from jussara are a promising technology that may provide another sustainable economic option for the tropical forest remnants.

Conflicts of Interest:
The authors declare no conflict of interest.