1. Introduction
Almond trees (
Prunus dulcis) are the most cultivated tree nut in the world, with a production of more than 1.5 metric tons of whole nuts in the 2022/2023 season [
1]. The largest producer is the United States, followed by Australia and Spain, with an increasing production trend [
1].
The almond fruit is divided into a greenish-colored casing called the hull, a hard and porous intermediate layer called the shell, a brown skin that surrounds the nut, protecting it and preventing oxidation and contamination, as well as the nut itself. Combined, the hull, shell, and skin represent approximately 85% of the mass of the almond fruit, the hull being the most abundant part [
2,
3,
4]. Finding high-value-added applications for all parts of the almond fruit would provide value in the almond production chain. Additionally, finding useful applications would be in line with the goals of the United Nations (UN) 2030 agenda, more specifically with goals 9 and 12 to build resilient infrastructure, promote inclusive and sustainable industrialization, and foster innovation, as well as to ensure sustainable consumption and production patterns [
5]. Currently, almond components other than the nut are typically used as animal feed or in energy production [
3].
One promising option to add value to the almond production chain is the sustainable exploration of almond by-products as a potential long-term source of bioactive natural products [
3]. The shell has a high xylan content—which can be fractionated into cellulose, pentosans, and lignins—which has been studied for water resistance, demonstrating that the addition of almond shell particles to wood increased the waterproofness of the wooden panels studied [
6]. Lignocellulosic biomass from almond shells has also been suggested as a renewable organic carbon source as its combustion, compared to petroleum derivatives, emits fewer hydrocarbons and monoxides [
4]. Almond hulls are often overlooked as a by-product in almond production. However, they have been studied extensively for their phenolic constituents and antioxidant properties. Sang et al. [
7] isolated several compounds from almond hulls, including protocatechuic acid and catechin, both known for their strong antioxidant capabilities, with their 2,2-diphenyl-1-picrylhydrazyl (DPPH) Radical Scavenging Capacity being above 90%. Barreira et al. evaluated almond hulls for their antioxidant properties, obtaining a DPPH Radical Scavenging Capacity above 80% at a concentration of 0.5 mg/mL [
8]. The development of highly efficient and sustainable processes for identifying and extracting these types of compounds is needed to accelerate efforts in the area.
The U.S. Food and Drug Administration (FDA) recommends that all parameters that may interfere in a process be studied [
9]. To assess multiple parameters in a systematic, simultaneous, and unbiased approach, a Design of Experiments (DoE) protocol can be used. Multivariate optimizations not only allow for detection of synergy between traditional extraction parameters, but they also allow for a greener optimization as they avoid the additional experiments that would be needed to optimize parameters one at a time [
10]. With a focus on green extraction processes, low-energy consuming, fast, and efficient extraction techniques are needed for studying almonds. One example is Microwave-Assisted Extraction (MAE) [
11]. Moreover, the selection of the extraction solvent is important since it has the greatest impact on extraction efficiency and the greenness as such solvents create potentially hazardous waste [
10,
12,
13]. Finding extraction solvents that lead to efficient extractions, while maintaining green processes, is not a trivial task as most extraction solvents are derived from petroleum. One alternative is offered by Natural Deep Eutectic Solvents (NADESs) [
14,
15].
NADESs are eutectic mixtures of two or more natural compounds, generally from the primary metabolism of plants, which can remain liquid at room temperature [
16]. The advantages of NADESs can include their easy preparation, biodegradability, renewability, biocompatibility, and ability to extract secondary metabolites [
12]. In 2018, the company Naturex
® (U. S. Patent No. 0055904, 2018) filed a patent reporting the use of a eutectic solvent composed of natural molecules for the efficient extraction of plant materials, demonstrating the practical and commercial potential of NADESs [
17]. Additionally, Naturex
® has a line of NADES-based extracts for the cosmetic industry known as Eutectys™ [
18]. NADESs can be classified into five categories, depending on the nature of the compounds used in their synthesis: (i) ionic liquid type, composed of an acid and a base; (ii) neutral type, formed by sugars or a combination of sugars and polyalcohols; (iii) neutral-acid type, synthesized from a sugar or polyalcohol mixed with an organic acid; (iv) neutral-basic type, with a sugar or polyalcohol mixed with an organic base; and v) amino acid type, which contains amino acids combined with sugars or organic acids [
19]. This classification reflects the diversity of NADESs available, making it possible to optimize extractions based on the matrix and bioactive compounds desired by selecting different NADESs. Continuous advances in knowledge and understanding of the characteristics and applications of NADESs have driven the exploration of these solvents as promising alternatives for various industrial purposes [
15,
19]. An additional advantage of using NADESs as an extraction solvent is that those made from food components (those that are formulated with natural ingredients and safe for human consumption) can be used for the delivery of bioactive compounds directly without removal of the solvent [
17].
NADESs have great potential as alternative solvents, especially in the removal of bioactive compounds from plants. They improve the degradation capacity, solubility, stability, bioactivity, and bioavailability of these compounds [
20]. Gomez-Urios et al. [
21] evidenced that orange peel extracts prepared with NADESs contained more bioactive compounds, such as catechin and caffeic acid, compared to ethanol, avoiding a high layer between these bioactive compounds and the eutectic solvents. Another study also highlighted the use of NADESs to improve thermal and storage stability in addition to the antimicrobial activity of catechins. The interaction between NADESs and catechins forms hydrogen bonds, explaining the high stability of catechins in these solvents [
22].
One of the disadvantages of NADESs is their non-volatile nature, due to high viscosity and high boiling points, making them difficult to remove after extraction [
23]. Despite this difficulty, studies show great potential to produce ready-to-use extracts without the need to remove the eutectic solvent, increasing the bioavailability of phenolic compounds from plant matrices. [
23]
Herein, this work proposes an innovative green method for the extraction of by-products from almond hulls. The research strategy adopted was based on the principles of sustainable analytical chemistry, using a multivariate experimental optimization (called central composite design, CCD), reducing the overall waste and resources needed for the method development. MAE using NADESs as extracting solvents (NADES–MAE) was evaluated as an improved extraction method for identifying and collecting value-added compounds from almond by-products. The feasibility and effectiveness of this method was thoroughly evaluated, and directly compared to established methods from the literature, both in terms of analytical extraction performance and greenness, using two complementary multiparametric green metrics.
3. Materials and Methods
3.1. Chemicals
Extraction: acetone (AR grade, Synth®, Diadema, SP, Brazil), ethanol (EtOH), and methanol (MeOH), all HPLC grade (LiChrosolv, MerckKGaA, Darmstadt, Germany), ultra-purified water from a Milli-Q system (Millipore, Burlington, MA, USA), glycerol (GLY) ≥ 99.5%, d-glucose (GLU) ≥ 99.5%, sucrose (SAC) ≥ 99.5%, lactic acid (LA) ≥ 85%, d-sorbitol (SOR) ≥ 98%, d-malic acid (MA) ≥ 98%, choline chloride (CL) ≥ 98%, d-fructose (FRU) ≥ 99%, l-proline (PRO) ≥ 99%, l-menthol (MEN) ≥ 99%, thymol (THY) ≥ 99% (all Sigma Aldrich, St. Louis, MO, USA).
Mobile phase HPLC: ethanol (EtOH, HPLC grade) (LiChrosolv, MerckKGaA, Darmstadt, Germany), ultra-purified water from a Milli-Q system and acetic acid (AcOH, HPLC grade) (LiChrosolv, MerckKGaA, Darmstadt, Germany).
3.2. Plant Material
Whole almond nuts including their skins, shells, and hulls were obtained from the producers Yunis Pty (Salisbury North, SA, Australia), CMV farms (Adelaide, SA, Australia), and Taronga Almonds (McLaren Vale, SA, Australia) from summer and autumn harvests in 2019–2021. Almonds were stored at 4 °C in closed containers prior to having the hull and shells separated.
3.3. Sample Preparation
The hulls were separated, cleaned, placed in equal portions, and frozen at a temperature of −20 °C. They were subsequently freeze-dried in a Christ freeze dryer (model Alpha 12LDplus Gefriertrocknungsanlagen, Osterode am Harz, Germany) for a period of 7 days. After lyophilization, the resulting material was pulverized in a knife mill and separated into a powder with particle sizes between 250 and 850 μm.
3.4. Reference Extraction Methods
Three reference methods of extracting secondary metabolites from almond by-products that have multiple citations in the literature were reproduced, with adaptations. Firstly, the method from Pinelo et al. [
24] using methanol and water with an acidic pH as an extracting solvent was used. Secondly, the method from Rubilar et al. [
36] using ethanol as the extraction solvent was used. Lastly, the method from Meshkini [
35] using acetone as the extraction solvent was reproduced. In all methods, dynamic maceration was applied as the extraction technique. Extractions were adapted to microscale and performed using a temperature-controlled Heidolph
® magnetic stirrer (Schwabach, Baviera, Germany) (using 0.1 g of powder, following the solvent volume, temperature, and rpm from each reference).
After dynamic maceration, the extracts were centrifuged at 7200 rpm in a mini-centrifuge (BioPet®, Knoxville, TN, USA), and approximately 1 mL of supernatant was collected and filtered through a 0.45 µm PTFE microfilter (Phenomenex®, Torrance, CA, USA) into 2 mL vials for further analysis by High-Performance Liquid Chromatography (HPLC) coupled to a PhotoDiode Array Detector (HPLC–DAD).
3.5. High-Performance Liquid Chromatography (HPLC–DAD) System
HPLC–DAD analyses were carried out on a Shimadzu system (LC-20AT pump, DGU-20A5R degasser, SIL-20HT sampler, CTA-10AS VP column oven, and interface CBM-20A, Shimadzu, Kyoto, Japan), equipped with a Luna column (2) C18-Phenomenex (250 × 4.6 mm), and a photodiode array detector (DAD, SPD-M20A), using 1% v/v AcOH in water and EtOH as the mobile phases, A and B, respectively. The samples were analyzed in gradient mode, with B varying from 5% to 60% in 60 min at a flow rate of 1 mL∙min−1. The oven temperature was 30 °C. A 50 µL injection volume was used. The detection wavelength was 254 nm. The equilibration was achieved under the initial conditions of the gradient, with 5% B, a flow rate of 1 mL∙min−1, and the column at 30 °C, for 15 min. The wash out was 10 min, with 100% B passing through the column.
3.6. Identification of Compounds Present in the Hydroethanolic Extract of Almond Hulls
The identity of the peaks in the chromatogram was confirmed from UV spectra obtained using the HPLC–DAD system and by liquid chromatography coupled to high-resolution mass spectrometry (LC–HRMS).
LC–HRMS analyses were performed using the conditions described in
Section 3.5 using an Agilent 1200 LC–6520 QTOF HPLC system with an electrospray ionization (ESI) (Agilent Technologies, Santa Clara, CA, USA), source in negative mode. A capillary voltage of 3.5 kV and N
2 gas flow rate (10 L∙min
−1) at 280 °C was used. The nebulizer was set to 45 psi; the fragmentor voltage was 155 V. N
2 was used as the collision gas (30 eV). MassHunter B.07 software (Agilent
®) was used to acquire and process the data acquired.
3.7. Preparation of Natural Deep Eutectic Solvents (NADESs)
The preparation of NADESs was based on the methodology of Gomez et al. [
37] via microwaves, which consists of subjecting the mixture to microwave radiation, varying and optimizing parameters such as irradiation time, temperature, and power. Two or three starting components at specific molar ratios and with water (
Table 8) were placed in borosilicate vials (4 cm i.d., 5 mL), then introduced into Teflon tubes (24.5 cm × 5 cm) prior to being subjected to microwave irradiation (Ethos Easy, Milestone SrL, Milan, Italy) at low power (200 W) in 5 min cycles (3 min of ramp to reach the programmed temperature of 40 or 50 °C followed by 2 min of irradiation with a stable temperature). The cycles were repeated for each NADES until a homogeneous and transparent solution was obtained. The mass percentage of water in the eutectic solvents was 20% (
w/
w), as adding water in quantities exceeding this amount can result in the disruption of hydrogen bonds present. This, in turn, leads to the dissolution of the supramolecular structure of the NADES, as observed in previous studies [
30,
31]. The resulting solutions were used as extraction solvents.
3.8. Conditions Used in NADES Screening
To identify the NADESs with the greatest extractive capacity, all studied NADESs (
Table 8) were used to extract secondary metabolites present in almond by-products, using microwave irradiation in triplicates. For every 0.1 g of almond hull powder, 1 mL of NADES extraction solution was added in a borosilicate vial (4 cm i.d., 5 mL) and introduced into Teflon tubes (24.5 cm × 5 cm). All extractions were carried out under microwave irradiation for 60 min as follows: 10 min of ramp to reach the programmed temperature, followed by 50 min of irradiation with a stable temperature at 50 °C and 500 W. Due to the high viscosity of most NADESs, after extraction, the mixtures were diluted in EtOH:H
2O (7:1,
v:
v) so that they all had the same final concentration, placed in 1.5 mL microtubes (Eppendorf, Hamburg, Germany, Olen, EUA), and centrifuged in a microcentrifuge (Biopet Technologies, Brazil) for 3 min at 7200 rpm. This process was carried out twice to obtain a more homogeneous supernatant. After centrifugation, 500 µL of supernatant was collected, filtered through a 0.45 µm PTFE microfilter (Phenomenex, USA), and analyzed with the HPLC–DAD method described in
Section 3.5.
3.9. Experimental Design
A complete 2
3 factorial design with repetition at the central point was carried out to evaluate the significant variables and perform a curvature test for subsequent optimization. The independent variables studied are presented in
Table 2.
The response used in this study was the sum of the peak areas corresponding to the target compounds (tannins and catechin monomers) obtained from the HPLC–DAD analysis with detection at λ = 254 nm. Furthermore, the response (total peak area of the chromatogram) was evaluated at a 95% significance level. A central composite design (CCD), using the significant variables, was performed to globally optimize the extraction conditions. At this stage, the power was set at 500 W, according to the power used in the NADES screening (
Section 3.8).
Analysis of Variance (ANOVA) calculations, and Response Surface plots (RSM) were performed using Statistica 10 software (Stat-Ease, Inc., Minneapolis, MN, USA), Origin Pro version 8.5, and Excel 2013 (Microsoft, Washington, DC, USA).
3.10. Assessment of Method Greenness
The comparison between the method developed in this work and the reference methods was carried out using two green metrics: Analytical GREEness (AGREE) (v.0.4.2020, Vigo, Pontevedra, Spain) and Green Analytical Procedure Index (GAPI) ComplexGAPI (v.0.2 beta Gdańsk, Poland).
The AGREE metric uses pictograms with colors that vary from green to red and shades that are calculated on a scale of 0 to 1, 1 being a highly green method (represented by dark green coloring) and 0 representing a lack of green chemistry principles (represented by red color) [
38].
The GAPI metric also uses pictograms with a color scale to classify the degree of greenness of each step of an analytical procedure, with three levels of evaluation for each step. GAPI has five pentagrams to evaluate and quantify the environmental impact as low, medium, and high at each stage of the methodology, using the colors green, yellow, and red, respectively. Each field reflects a different aspect of the described analytical procedure and is filled in green if certain requirements are met [
39,
40].