Food-Safe Process for High Recovery of Flavonoids from Cocoa Beans: Antioxidant and HPLC-DAD-ESI-MS/MS Analysis

Considering the increasing interest in the incorporation of natural antioxidants in enriched foods, this work aimed to establish a food-grade and suitable procedure for the recovery of polyphenols from cocoa beans avoiding the degreasing process. The results showed that ultrasound for 30 min with particle sample size < 0.18 mm changed the microstructure of the cell, thus increasing the diffusion pathway of polyphenols and avoiding the degreasing process. The effect of temperature, pH, and concentration of ethanol and solute on the extraction of polyphenols was evaluated. Through a 24 full factorial design, a maximum recovery of 122.34 ± 2.35 mg GAE/g, 88.87 ± 0.78 mg ECE/g, and 62.57 ± 3.37 mg ECE/g cocoa beans, for total concentration of polyphenols (TP), flavonoids (TF), and flavan-3-ols (TF3), respectively, was obtained. Based on mathematical models, the kinetics of the solid–liquid extraction process indicates a maximum equilibrium time of 45 min. Analysis by HPLC-DAD-ESI-MS/MS showed that our process allowed a high amount of methylxanthines (10.43 mg/g), catechins (7.92 mg/g), and procyanidins (34.0 mg/g) with a degree of polymerization >7, as well as high antioxidant activity determined by Oxygen Radical Absorbance Capacity (1149.85 ± 25.10 µMTrolox eq/g) and radical scavenging activity (DPPH•, 120.60 ± 0.50 µM Trolox eq/g). Overall, the recovery method made possible increases of 59.7% and 12.8% in cocoa polyphenols content and extraction yield, respectively. This study showed an effective, suitable and cost-effective process for the extraction of bioactive compounds from cocoa beans without degreasing.


Introduction
Flavanols are the most abundant substances of polyphenols found in cocoa, with a degree of polymerization ranging from monomers to polymeric proanthocyanidins [1]. Therefore, unfermented cocoa bean is composed of 1.3-3.3% methylxanthines [2] and of about 6% condensed flavan-3-ols [3]; as a result, it is has been listed as the 4th richest dietary source of polyphenols [4]. Cocoa is one of the top crops in Colombia in terms of economic impact, with a national production record of 59,665 tons in 2019 according to National Cacao Producers Federation. There is also a government initiative aiming

Pre-Treatment of Cocoa Beans
Fresh cocoa pods (Trinitary, clone ICS 39) were collected at Villa Santa Monica (San Vicente de Chucurí, Santander, Colombia) and immediately protected from light and transported on ice to CICTA Lab for processing. Beans were removed manually from the pods using a knife, cleaned and separated from the pulp (mucilage) using a mucilage remover (Penagos Ltd.a, Colombia). After that, beans were heat-treated (96 • C for 6.4 min) with an enzyme inhibition solution (70 mM·L-ascorbic acid and l-cysteine, ratio 1:1 v/v) until reaching a polyphenol oxidase inhibition up to 93%. The inactivated beans were frozen at −20 • C and used for further analysis.

Drying and Milling Process
Drying and milling processes were evaluated as follows: (a) Drying: Beans after polyphenol oxidase (PPO) inactivation were used immediately to evaluate the effect of drying technology on the total polyphenol content. To do so, beans were chopped (cross section of 50 × 30 mm 2 ). Then, beans were a) oven dried (FD 23, Binder, Germany) at both 50 • C and 70 • C at atmospheric pressure, and b) freeze-dried (Labconco Corp., Kansas City, MO, USA) at −84 • C processing temperature and 13 Pa constant pressure in the drying chamber to obtain final humidity < 4%. The moisture content was evaluated by AOAC method 931.04 (AOAC, 1990). (b) Milling: As a strategy to avoid the use of a non-GRAS (Generally Recognized as Safe) solvent for removing the fat from cocoa beans, different particle size distributions and ultrasound time were evaluated as a function of total phenol content. Dried beans were milled at max speed at −20 • C in N 2 environment for 30 s during three cycles (Grindomix GM 200, Retsch, Haan, Germany). The milled samples were sieved through steel mesh (W.S. Tyler, Mentor, OH, USA) with a sieve shaker (Gilson, Screen Co., USA) and fractionated in three groups: sieved and retained on 20 to 40-mesh (sample 1); 40 to 80-mesh (sample 2); and 80 to 200-mesh (sample 3). After that, the powdered sample was immersed in 50% aqueous ethanol and ultrasonicated (35 kHz, ice bath at 4 • C, Elma, Ultrasonic LC20H, Germany) for several intervals of time. A defatted cocoa sample (<5 wt.%) was employed as a control sample. To do so, the cocoa bean powder (1.0 g) was three times defatted with n-hexane (10 mL, extraction in an ultrasonic bath at 25 • C for 15 min). The resulting powder was dried overnight at room temperature.
In all experiments, a standard polyphenol extraction step was carried out as follows: 1 g of sample was added to 60 mL of a mixture of 50% ethanol and water (w/w) at 50 • C with stirring at 300 rpm for 30 min using a magnetic stirrer hotplate (IKA C-MAG HS7, Germany) with temperature being monitored with a thermocouple (IKA ETS-D5, Germany). The resulting extract was centrifuged (5000× g, 4 • C, 20 min), then the supernatant was filtered through 0.45 µm hydrophilic PTFE filter (Millipore, Milford, USA), and the total polyphenol content was measured immediately as detailed in Section 2.4.1.

Scanning Electron Microscopy
Scanning electron microscopy (SEM) was used to evaluate the microstructure of (a) raw cocoa beans, (b) cocoa beans with reduced PPO activity, and (c) cocoa beans with reduced PPO and treated by ultrasound. Bean samples cut into longitudinal and transversal sections were mounted on aluminum stubs with double-sided carbon adhesive tape and observed using the XL-30 Environmental SEM (Philips, USA) at 25 kV accelerating voltage with the BSE (backscattered electron) detector. The images were stored in TIFF format at 712 × 484 pixel in grayscale with brightness values between 0 and 255 for each pixel constituting the image.

Solid-Liquid Extraction of Polyphenols
Extraction temperature, pH, solute/solvent ratio, and ethanol/water ratio were evaluated as the major factors that can affect the extraction yield. pH was adjusted with 1 M citric acid. The combination of these factors was modeled through a surface design consisting of 2 4 + four replicates at the central point + triplicates at the start point. Low and high levels for the different factors were, as follows: Temperature (40, 60 • C); pH (3,5); solute/solvent ratio (w/v) (1/60, 1/30); ethanol/water ratio (v/v) (25/75, 75/25).
The resulting extracts were centrifuged at 5000× g, 4 • C for 15 min and the supernatant was collected and filtered using Whatman Nº1 filter paper (Whatman, Inc., USA), and used for further experiments. At maximum extraction conditions, large-scale recovery (up to 10 L) controlled by Bioflow-110 bioreactor (New Brunswick Scientific, Enfield, CT, USA) was used to prove that the scale up of the process is possible.

Total Polyphenol Content by Folin-Ciocalteu
The total polyphenol (TP) content was assayed using the Folin-Ciocalteu reagent according to Singleton et al. [19] with modifications as follows. The reaction was initiated by the addition of 50 µL of the sample with 1.5 mL of 10-fold diluted Folin-Ciocalteu reagent. After 5 min, 1.5 mL of 7.5% (w/v) sodium carbonate was added and vortexed for 10 s. The reaction medium was stored in the dark for 1 h at 25 • C. Absorbance was measured at 765 nm (Genesys 20; Thermo Scientific-Fisher, Waltham, MA, USA) against a blank sample. A calibration curve of gallic acid (ranging from 0.01-0.8 mg/mL, r 2 = 0.99) was prepared, and the results were expressed as mg gallic acid equivalents (GAE) per gram of dried cocoa beans.

Total Flavonoid Assay
Total flavonoid (TF) content was measured according to Zhishen et al. [20]. Sample (500 µL) was added to 5 mL volumetric flask containing 2 mL H 2 O followed by addition of 0.15 mL 5% NaNO 2 . After 5 min, 0.15 mL 10% AlCl 3 was added and 1 min after, 1 mL 1M NaOH was added. The total volume was made up to 5 mL with H 2 O. The reactants were mixed well and stored in the dark for 15 min at 25 • C. Absorbance was measured at 510 nm (Genesys 20; Thermo Scientific-Fisher, Waltham, MA, USA) against a blank sample. Total flavonoid content was expressed as mg (−)-epicatechin equivalents (ECE) per gram of dried cocoa beans.

Total Flavan-3-Ol Assay
Total flavan-3-ol (TF3) content was determined by the vanillin-H 2 SO 4 assay as described by Sun [21]. The reaction consisted of 1 mL sample in methanol with 2.5 mL of 1% vanillin in methanol and 2.5 mL of 9N H 2 SO 4 . The reaction medium was well mixed at 30 • C and allowed to react for 15 min. Absorbance was measured at 500 nm (Genesys 20; Thermo Scientific-Fisher, Waltham, MA, USA) against a blank sample. Total flavan-3-ol content was expressed as mg (−)-epicatechin equivalents (ECE) per gram of dried cocoa beans.

Kinetic of Solid-Liquid (S-L) Extraction of Polyphenols
Once the optimum S-L extraction conditions were selected, the extraction kinetics were evaluated by plotting the concentration of the isolated target analyte versus time. Aliquots (100 µL) of each sample were taken out at various times to measure the total concentration of polyphenol, flavonoid, and flavan-3-ol. The extraction curves were adjusted by several kinetic models previously documented [22][23][24][25][26].
The mass detector was an Agilent 6320 Ion Trap LC/MS (Agilent Technology, Waldbronn, Germany) equipped with an ESI source and ion trap mass analyzer which was controlled by the 6300 series trap control software (Bruker Daltonik GmbH). MS/MS analyses were carried out to obtain the structural information of the separated compounds. Mass spectrometer was operated under positive and negative ESI mode, nebulizer pressure, 40 psi; dry gas, 12 L min −1 ; dry temperature, 350 • C, and mass spectra recorded from 90 to 2200 m/z. MS characterization features were analyzed using extraction ion compound tool and commercial standards. It was also consulted METLIN and HMBD databases for matching exact mass.

Antioxidant Assays
Oxygen radical absorbance capacity (ORAC) assay described by Huang et al. [28] was carried out as follow: 96 well microplates were filled with 50 µL of the daily working fluorescein solution (4 × 10 −3 mM in 75 mM phosphate buffer, pH 7.4), 50 µL of samples at a known concentration, and incubated at 37 • C for 10 min in a microplate reader (Synergy HT Multi-Detection, Biotek Instruments, Inc. Winooski, VT, USA). The reaction was initiated by the addition of 50 µL of 2,2 -Azobis(2-methylpropionamidine) dihydrochloride (221 mM in 75 mM phosphate buffer, pH 7.4) and the fluorescence decay was monitored kinetically for 2 h, using emission and excitation wavelength of 485 nm and 528 nm, respectively. A calibration curve was prepared with 12.5-375 µM Trolox (r 2 = 0.99). Results are expressed as µM Trolox equivalents per g of cocoa beans (dry matter basis) as follows: where k is the sample dilution factor, and AUC is the area below the fluorescence decay curve of the sample, blank, and Trolox, respectively. The area under the curve of normalized data was calculated using GraphPad Prism v. 6.0 (GraphPad Soft. Inc., La Jolla, San Diego, CA, USA). The reducing ability of antioxidants toward the DPPH (2,2-diphenyl-1-picrylhydrazyl) radical was measured according to Brand-Williams et al. [29]. DPPH methanolic solution (684.7 µM) was adjusted with methanol until an absorbance of 0.550 ± 0.01 units at 517 nm (Genesys 20, Thermo Scientific, Waltham, MA, USA) was obtained. Samples (100 µL) were incubated with 1.45 mL of this DPPH solution for 60 min in the dark. A calibration curve was prepared from 1.95-150 µM Trolox (r 2 = 0.99). Results are expressed as µM Trolox equivalent per g of cocoa beans (dry matter basis).

Statistical Analysis
All the experimental measurements were repeated at least three times, and data were expressed as the mean ± standard deviation. Statistical analysis was done using Statistica 7.1 (Stat-Soft Inc., USA). The analysis of variance (ANOVA) and F-test were used to evaluate the influence of the factors and their interactions on the experimental designs. ANOVA one-way and Tukey's multiple range test was also conducted at a 5% level of significance. The kinetic constants in this study were determined from experimental data using non-linear regression employing Quasi-Newton method and least squares as custom loss function (0.0001 for both start value and initial step size). The response surface methodology consisting of full factorial central composite rotatable design with four replicates at the central point was conducted according to a completely randomized model. A second-order polynomial equation was used to fit the experimental data as follows: where Y is the predicted factor, β 0 is the value of the fitted response to the design, β i , β ii , and β ij are the coefficients of linear, quadratic, and cross product terms, respectively. The performance of full factorial central composite design method was measured by r and r 2 . Experimental runs were also randomized to evaluate the concordance of experimental data and predicted values; therefore, the root mean squared error (RMSE) was calculated as follows: where y i andŷ i is the measured value and predicted value by the model, respectively. In addition, n is the number of the set data.

Drying Technology
Beans with reduced enzyme activity were manually chopped, and the total phenol content was evaluated after (a) oven-air drying at 50 • C and 70 • C, and (b) freeze-drying. The initial moisture content in a fresh unfermented bean was 40.1 wt.% and after drying process the final moisture content was between 1.4-3.2 wt.%. Table 1 summarizes the effect of drying process on TP content. Therefore, high temperature air drying led to a great extent of dehydration levels. Visual evidence showed that dried beans were slightly grey/purple color, which is characteristic of dried unfermented cocoa beans with high procyanidin content [12]. However, drying at 50 • C for 20 h resulted in TP content 10% lower than drying at 70 • C for 3 h (p < 0.05), which showed that polyphenols were more sensitive to drying time than temperature. No significant differences in TP amount was observed between freeze drying and oven-dried at 70 • C treatments, which was also reported by Kothe et al. [30], who found that when employing high temperatures and short times (100-160 • C for 30 min), different reactions can occur, such as configuration rearrangement of flavan-3-ols and epimerization, which could led to the formation of new flavan-3-ol products with higher phenol activity. Thereby, oven-air drying at 70 • C was chosen as drying technology for further experiments. Moreover, air drying is a more cost-effective solution (and faster) than freeze-drying.

Impact of Particle Size on Extraction Yield
After the PPO enzymatic inhibition and drying steps, the samples were milled at low temperature for a short time, thus avoiding the fat melting during the reduction of particle size. In an attempt to discriminate which sizes led to better performance on extraction yield, three distinct samples varying on particles size were evaluated as follow: (S 1 ) > 0.42 mm, (S 2 ) 0.42-0.18 mm, and (S 3 ) < 0.18 mm. The results showed that the best treatment was S 3 with a TP recovery of 50 mg GAE/g being 2× and 1.2× fold higher than S 1 and S 2 , respectively. This study reinforces the idea that there is an inverse relationship between the particle size and the extraction yield. This behavior was also observed by Sun et al. [31], who found that extraction yields were 6 fold higher when particle size decreased from 2 to 0.074 mm.

Conditions to Avoid the Degreasing Process
Non-defatted samples and defatted samples with a total fat content of 58.5% and 2.5%, respectively, were used for these experiments. In prior assays, results showed that TP extraction from the defatted samples was 60% higher than in non-defatted samples, which confirms that fat is a barrier for the diffusion of the bioactive compounds from the solid phase to the extractor solvent. Since one of the aims of the present work was to avoid the use of toxic organic solvents, ultrasonic treatment was selected as an alternative clean extraction technology. Thus, non-defatted cocoa samples were placed in a test tube and mixed with 50% ethanol (1/60 w/v ratio solute/solvent) and ultrasonicated (35 kHz, 4 • C). The effect of ultrasonic time (0-60 min) was tested as a function of TP amount. As can be seen in Figure 1, ultrasound increases the polyphenols extraction efficiency up to a 30 min. This is because the cavitation forces improving the polarity of the system; moreover, the bubbles in the liquid/solid extraction can explosively collapse and generate localized pressure causing plant tissue rupture [32], thus increasing the mass transfer rate. Figure 1 also shows that samples ultrasonicated for 30 min had no significant difference with the defatted control (p < 0.05). Thereby, the ultrasonic bath for 30 min was chosen as optimal value for further experiments. 7 Mean value ± standard deviation (n = 5) with similar letters mean are not significantly different by ANOVA One-way Tukey (p > 0.05). *The thermophysical property data for the determination of enthalpy of sublimation of pure compounds (water) was calculated thorough the libraries of the NIST ThermoData Engine (TDE) version 10.

Impact of Particle Size on Extraction Yield
After the PPO enzymatic inhibition and drying steps, the samples were milled at low temperature for a short time, thus avoiding the fat melting during the reduction of particle size. In an attempt to discriminate which sizes led to better performance on extraction yield, three distinct samples varying on particles size were evaluated as follow: (S1) > 0.42 mm, (S2) 0.42-0.18 mm, and (S3) < 0.18 mm. The results showed that the best treatment was S3 with a TP recovery of 50 mg GAE/g being 2× and 1.2× fold higher than S1 and S2, respectively. This study reinforces the idea that there is an inverse relationship between the particle size and the extraction yield. This behavior was also observed by Sun et al. [31], who found that extraction yields were 6 fold higher when particle size decreased from 2 to 0.074 mm.

Conditions to Avoid the Degreasing Process
Non-defatted samples and defatted samples with a total fat content of 58.5% and 2.5%, respectively, were used for these experiments. In prior assays, results showed that TP extraction from the defatted samples was 60% higher than in non-defatted samples, which confirms that fat is a barrier for the diffusion of the bioactive compounds from the solid phase to the extractor solvent. Since one of the aims of the present work was to avoid the use of toxic organic solvents, ultrasonic treatment was selected as an alternative clean extraction technology. Thus, non-defatted cocoa samples were placed in a test tube and mixed with 50% ethanol (1/60 w/v ratio solute/solvent) and ultrasonicated (35 kHz, 4 °C). The effect of ultrasonic time (0-60 min) was tested as a function of TP amount. As can be seen in Figure 1, ultrasound increases the polyphenols extraction efficiency up to a 30 min. This is because the cavitation forces improving the polarity of the system; moreover, the bubbles in the liquid/solid extraction can explosively collapse and generate localized pressure causing plant tissue rupture [32], thus increasing the mass transfer rate. Figure 1 also shows that samples ultrasonicated for 30 min had no significant difference with the defatted control (p < 0.05). Thereby, the ultrasonic bath for 30 min was chosen as optimal value for further experiments.  fractured, as well as appearing more solid and denser due to the cellular contents remaining embedded in the cell (Figure 2A). In general, the cocoa bean is formed of parenchyma cells (containing cocoa butter and proteins), which represent over 80% of the mass [33], and polyphenols and alkaloids, which are placed into the vacuoles [34]. Physical changes were observed in the structure of cocoa starch granules that become larger, smooth and fibrous, quite likely as produced by heat treatment ( Figure 2B). This behavior could be the result of enzyme heat treatment, which opens up the cell and favors diffusion of bioactive compounds from the matrix into the extraction medium, as well as to the mobilization of proteins and polyphenols and redisposition of fat within the cell [34].

Solid-Liquid Extraction of Polyphenols from Cocoa Beans
As mentioned, temperature, solute/solvent ratio, ethanol/water ratio, and pH were evaluated through a surface experimental design. As can be seen in Table 2, the concentration of ethanol had a strong impact on response factors, for instance, TP, TF and TF3 varied from 41.3-107.6 mgGAE/g, 15.2-86.0 mgECE/g, and 21.7-59.8 mgECE/g, respectively, throughout the ethanol concentration range (0-100%). Moreover, TP yield increased with increasing extraction temperature, i.e., extraction at 60 °C (25% ethanol, pH 5, 1/60 solvent/solute ratio) was 1.12 fold higher at 40 °C. Maximum values for TP (107.6 mg GAE/g), TF (86.0 mg ECE/g), and TF3 (59.8 mg ECE/g) were observed at different levels of factors, which indicate that combination of factors plays a key role in the extraction process. The morphological differences are more evident in Figure 2C. Smaller fragments were dispersed within the cell, and the microstructure was more porous, which could be the result of cell disruption and changes in the intercellular spaces. Hence, changes in the cell might increase the permeability and thus diffusion of polyphenols out of the cells. In accordance with our results, Rostagno et al. [35] found that disruption of tissue and cell walls are most efficient when heat treatment together with ultrasound are used, which resulted in a greater penetration of solvent into the sample matrix, increasing the contact surface between solid and liquid phase, and as a result, the solute quickly diffuses from the solid phase to the solvent. The observation of this study confirms that the enzymatic inhibition by heat treatment, milling process, and ultrasounds play a significant role in the change of the internal structure of cocoa beans, thus enhancing a high concentration of polyphenols while avoiding the defatting process.

Solid-Liquid Extraction of Polyphenols from Cocoa Beans
As mentioned, temperature, solute/solvent ratio, ethanol/water ratio, and pH were evaluated through a surface experimental design. As can be seen in Table 2, the concentration of ethanol had a strong impact on response factors, for instance, TP, TF and TF3 varied from 41.3-107.6 mg GAE/g, 15.2-86.0 mg ECE/g, and 21.7-59.8 mg ECE/g, respectively, throughout the ethanol concentration range (0-100%). Moreover, TP yield increased with increasing extraction temperature, i.e., extraction at 60 • C (25% ethanol, pH 5, 1/60 solvent/solute ratio) was 1.12 fold higher at 40 • C. Maximum values for TP (107.6 mg GAE/g), TF (86.0 mg ECE/g), and TF3 (59.8 mg ECE/g) were observed at different levels of factors, which indicate that combination of factors plays a key role in the extraction process.  where TP is total polyphenol TF is total flavonoid and TF3 is total flavan-3-ol. S-L is the solid-liquid extraction. Means (n = 3) within a column (comparison of same bioactive compound family) with different letters were significantly different by ANOVA One-way Tukey (p < 0.05). For more information see the methodology section.

T ( • C) SS (w/v) EW (v/v) pH TP (mgGAE/g) TF (mgECE/g) TF3 (mgECE/g)
Results were fitted to a quadratic model (Equations (S1)-(S3)) and the 3D surface graphs for dependent variables were plotted ( Figure 3). Thus, analysis of variance summarized the observations and differences among the studied factors (Table S1). ANOVA showed that the recovery of bioactive compounds was mainly dependent on the linear effect of temperature (T) and pH, linear and quadratic effect of ethanol/water (EW) ratio, and solute/solvent (SS) ratio. Besides, the interactions of T × EW, and EW × pH for total polyphenol content; SS × T, and SS × EW for total flavonoid concentration; and T × SS, T × EW, SS × pH, and EW × pH for total flavan-3-ol amount were also significant (p < 0.05). ANOVA analysis also showed that the selected quadratic model adequately represented the extraction process. Therefore, the model has a good coefficient of multiple determination of r 2 = 0.957, r adj = 0.941; r 2 = 0.954, r adj = 0.937; and r 2 = 0.905, r adj = 0.872 for total polyphenol, flavonoid and flavan-3-ols content, respectively. Overall, higher accordance regression fit values mean that the model explained most of the variability in the responses.

Effect of Independent Factors on the Recovery of Total Polyphenols and Total Flavonoids
As can be seen in Figure 3, the significant factors evaluated had a parabolic behavior on the extraction of TP, TF, and TF3. For instance, it was confirmed that temperature promoted faster diffusion rate, better mass transfer, reduce viscosity and surface tension ( Figure 3A,C,E), which is in agreement with Cacace and Mazza [36]. Indeed, this factor had a linear effect until it reached a maximum, after which the extraction yield decreased. This phenomenon is explained by the softening of plant tissue at high temperature, while at the same time, it weakens phenol-protein and phenol-polysaccharide interactions [37]. Moreover, results showed that higher diffusion rate is enhanced at lower ratio solute/solvent, which is confirmed by the significance of linear and quadratic effect for both total phenol and total flavonoid content. Also, Fick´s second law of diffusion predicts this phenomenon, that is, a final equilibrium between the concentration of solute in the solid matrix and in the bulk solution after a certain time [37]. Antioxidants 2020, 9,  where TP is total polyphenol TF is total flavonoid and TF3 is total flavan-3-ol. S-L is the solid-liquid extraction. Means (n = 3) within a column (comparison of same bioactive compound family) with different letters were significantly different by ANOVA One-way Tukey (p < 0.05). For more information see the methodology section. Results were fitted to a quadratic model (Equations S1-S3) and the 3D surface graphs for dependent variables were plotted (Figure 3). Thus, analysis of variance summarized the observations and differences among the studied factors (Table S1). ANOVA showed that the recovery of bioactive The influence of ethanol concentration ( Figure 3A,C,E) suggested that the extract contains polyphenols with different polarities, which also confirms the principle "like dissolves like" The lowest concentration of water or ethanol were not efficient for the extraction of flavonoids from cocoa beans. The maximum yield was enhanced at 50% ethanolic aqueous solution. Previous works are consistent with our results, reporting a maximum extraction yield of polyphenols with an ethanol/water ratio of 40-50% [18,38], which is because the concentration of ethanol influences the velocity-transfer kinetic and the thermodynamics of the process [38].
In addition, the interaction of pH with the solvent (Figure 3B,D,F) plays an essential role in controlling adsorption/desorption rate (i.e., changing the surface charge of adsorbent, the degree of ionization and speciation of adsorbate during adsorption) [39], and an impact on the pKa value for the polyphenols. In fact, lower solvent acidity allowed the highest recovery of flavonoids, because acid could hydrolyze the cell wall and thus facilitates the diffusion of the phenolic compounds. However, at higher concentrations of ethanol, a reduction on the extraction yield was observed. Similar results show that the amount of acid in the extraction solvent increases concentration by increasing the stability of the phenolics during extraction and increasing their dissolution [37].
In addition, significant differences (p < 0.05) among control samples and sample at optimum conditions were found; thereby, TP, TF, and TF3 were 59.7, 59.8, and 57.8% higher, respectively, when employing the optimized S-L extraction method (Table 2).

Extraction Kinetics Parameters
Given the maximum conditions for the recovery of polyphenols, different extractions were carried out to evaluate the equilibrium time and kinetic parameters as a function of TP, TF, and TF3 content.
The study of equilibrium time plays an essential role in economizing the energy and cost of the industrial process [37], thus improving the accuracy of the procedure and the quality of the final product. Table 3 summarizes the kinetic parameters, RMSE and the accordance of the model (r) for all the mathematical models selected, which were previously used to model the solid-to-liquid extraction of bioactive compounds. Overall, TP, TF, and TF3 content had the lowest r and highest RMSE for those equations that consider the extraction is occurring in one continuous step, i.e., Equations (1), (2), (4), and (5). Goodness of fit of the model and lower standard deviation of residuals were obtained for the two models that represent the recovery of polyphenols on two different rates (sorption/desorption), for instance, Equation (3) and Equation (6). These equations represent the kinetic process of a liquid/solid system based on the solid capacity. As can be seen in Figure 4, the extraction curve shape had a faster extraction rate followed by a slower extraction rate and asymptotically approached the equilibrium concentration. Generally, the temperature enhances the solubility of polyphenols in the solvent [40], but the rate of extraction yield decreased progressively due to the exposure to both higher temperature and time extraction, which confirms that polyphenols are thermosensitive compounds [41].
Results suggested that the Peleg model (Equation (3)) proved to be most suitable to model the solid-liquid extraction kinetics for the dependent variables (r ≥ 0.98 and RMSE ≤ 0.71). In fact, k 1 and k 2 represent the extraction rate constant and constant of extraction extent, respectively [42]. Table 3 shows that k 1 was similar in all the cases, but k 2 increased by TF3~TF >> TP, which is related to its maximum equilibrium concentration (C t→∞ ). Thereby, the equilibrium extraction time were 45 min, 39 min and 34 min for TP, TF, and TF3 content, respectively. Comparison of these data with those of previous authors shows that our optimized extraction method was 2.7×, 4.4×, 6.0×, and 2.7× faster than the polyphenol recovery of cocoa beans [43], grape seed [42], and mango kernel [44], respectively.
In addition, a large-scale extraction process was also assayed. The results showed that recovery from 0.05 to 5 L did not significantly impact extraction yield (p < 0.05), but the extraction at 10 L was 25.7 ± 0.9% lower for TP content while was 23.0 ± 6.9% and 5.3 ± 2.6% higher for TF and TF3 amount, respectively. These differences could be due to diffusional changes as a result of the scaling up the process.
In general, this work confirmed the importance of studying the impact of extraction parameters on both secondary plant metabolites, thereby increasing the yield of extraction. In terms of going beyond the highlights, the optimized ultrasound-assisted solid-liquid extraction not only allowed a high concentration of both total polyphenols and flavonoids but also employed food-grade solvents, reduced the number of stages (i.e., avoiding the degreasing), the extraction time (<45 min), and therefore energy consumption. Thereby, the extraction process could be suitable for large-scale applications. For example, cocoa polyphenol extract can be used to enrich products in very high demand in the food and cosmetic industries.  where TP is total polyphenol (mgGAE/g) TF is total flavonoid (mgECE/g), and TF3 is total flavan-3-ol (mgECE/g). Antioxidants 2020, 9, x FOR PEER REVIEW 14 of 19 Results suggested that the Peleg model (Equation 3) proved to be most suitable to model the solid-liquid extraction kinetics for the dependent variables (r ≥ 0.98 and RMSE ≤ 0.71). In fact, k1 and k2 represent the extraction rate constant and constant of extraction extent, respectively [42]. Table 3 shows that k1 was similar in all the cases, but k2 increased by TF3 ⁓ TF >> TP, which is related to its maximum equilibrium concentration ( →∞ ). Thereby, the equilibrium extraction time were 45 min, 39 min and 34 min for TP, TF, and TF3 content, respectively. Comparison of these data with those of previous authors shows that our optimized extraction method was 2.7×, 4.4×, 6.0×, and 2.7× faster than the polyphenol recovery of cocoa beans [43], grape seed [42], and mango kernel [44],

Chromatographic and Antioxidant Analysis
Identification of alkaloids, catechins and procyanidins in cocoa, were achieved by HPLC-DAD-ESI-MS/MS. The chromatographic method allowed not only the detection of main alkaloids but also the procyanidins with different degree of polymerization up to 7 ( Figure 5). As can be seen in Table 4, the cocoa extract obtained at maximum conditions contained 65% (w/w) procyanidins followed by methylxanthines (20% w/w) and catechins (15% w/w), which is in agreement with previous authors [45,46]. Theobromine/caffeine ratio was 2.93, which can be used to classify hybrid genotypes [47]. Therefore, it was confirmed that the cocoa sample was of the Trinitarian variety. (−)-Epicatechin (7.30 mg/g) was the main catechin, being 11.8 fold higher than (+)-catechin, which was also confirmed by Romero et al. [48]. The major procyanidin in cocoa polyphenol extract was the Trimer C1 (11.9 mg/g) being 2.9×, 1.3×, 1.7×, 6.6×-fold higher than the dimer, tetramer, pentamer, and hexamer, respectively.
Our results are in agreement with previous works that reported that the main alkaloid and flavanol are the theobromine and (−)-epicatechin, but their concentrations vary considerably depending on the methodology of extraction and cocoa bean variety [49]. Theobromine and caffeine concentrations were in agreement with previous reports [8,15], but epicatechin and catechin amounts were 1.5× and 1.9×, and 2.0× and 2.3× fold higher than those reported by Kothe et al. [30] and Carrillo et al. [15], respectively. In general, procyanidin content was much higher than previous authors, i.e., dimer B2 was 2.0× and 1.5× -fold greater than procyanidin B2 previously reported by Kothe et al. [30] and Tomas-Barberán et al. [43], respectively.
In addition, the identification of oligomeric procyanidins was carried out by mass spectrometry in comparison with commercial standards and published literature. For instance, dimer B2 (B-type, EC-4β→8-EC-4β), trimer C1 (B-type, EC-4β→8-EC-4β→8-EC), tetramer, pentamer, hexamer and heptamer with a molecular ion [M − H] − m/z 577, 865, 1153, 1441, 1729, 2019, respectively were identified. Characterization of larger polymers was not possible due to their low concentration, low ionization, peak signal dispersion, the formation of multiple ions, and limitations of the ion trap MS analyzer. Fragment patterns also suggested that various fragmentation mechanisms are involved in ESI such as quinone methide (QM), retro-Diels-Alder (RDA), as well as heterocyclic ring fission (HRF), which could take place on the extension unit or the terminal unit of the molecule [50]. For instance, the loss of a fragment with m/z 152, 170 and additional loss of water corresponding to RDA fission and loss of 288 Da corresponds to QM cleavage ( Figure 5). 15 [45,46]. Theobromine/caffeine ratio was 2.93, which can be used to classify hybrid genotypes [47]. Therefore, it was confirmed that the cocoa sample was of the Trinitarian variety. (-)-Epicatechin (7.30 mg/g) was the main catechin, being 11.8 fold higher than (+)-catechin, which was also confirmed by Romero et al. [48]. The major procyanidin in cocoa polyphenol extract was the Trimer C1 (11.9 mg/g) being 2.9×, 1.3×, 1.7×, 6.6×-fold higher than the dimer, tetramer, pentamer, and hexamer, respectively. Data expressed as means of triplicate experiments. *Concentration expressed as mg polyphenol per g cocoa beans (dry weight basis).   Among the antioxidant assays, our results showed that the cocoa extract had an ORAC of 1149.85 ± 25.1 µM Trolox eq/g. These values are higher than those reported by Hurst et al. [51] and Carrillo et al. [15], with TP in the range 58.0-61.7 and 45.3-70.0 mg GAE/g, and ORAC in the range 797.0-947.0 and 387.3-618.1 µM Trolox eq/g, respectively. Moreover, DPPH radical scavenging activity was 120.6 ± 0.5 µM Trolox eq/g (equal to 0.72 µM Trolox eq/mg cocoa extract), which was 2.4 fold higher than previously reported by Summa et al. [52]. These differences could be associated with the improved methodology of extraction at optimal conditions, thus enhancing both high composition and concentration of catechins (7.92 mg) and procyanidins (34.0 mg) with a DP ≥ 7.

Conclusions
An effective ultrasound-assisted solid-liquid process for extracting polyphenols from cocoa beans was optimized using an experimental design. The 2 4 factorial design showed that all parameters studied were significant factors in affecting the polyphenolic content as well as enabling us to determine the optimal values for the extraction process. This optimization showed that the best conditions to obtain high polyphenol yield were 50% of ethanol, solid/solvent ratio of 1:120 w/v, pH 6 at 70 • C for a maximum equilibrium time of 45 min. Operating conditions to avoid degreasing and freeze-drying steps were also established, thus leading to a more cost-effective strategy. Overall, the process extraction allowed to increase on 59.7% and 12.8% of cocoa polyphenols amount and extraction yield, respectively. The extract rich in polyphenols could replace the synthetic antioxidants and could be used in the food and cosmetic industries. Our results suggest that ultrasound-assisted solid-liquid extraction is a suitable method for the recovery of cocoa polyphenols and could be used for the scaleup procedure for further research.