Efficient extraction of oleoresin from Ferula gummosa roots by natural deep eutectic solvent and its structure and chemical characterizations

Deep eutectic solvents in the extraction of plant metabolites have found many advantages, such as low toxicity, biodegradability, low cost and ease of preparation over the conventional methods. This work aims to compare natural deep eutectic solvents in extraction and optimization of oleoresin from Ferula gummosa and determining its chemical and structure properties. Box–Behnken design was applied to optimize the extraction of oleoresin from Ferula gummosa using eutectic solvents. The variables of extraction were extraction time, temperature, and ratio of eutectic solvents. Six mixtures of eutectic solvents including choline chloride/urea, acetic acid, lactic acid, formic acid, formamide and glycerol at ratios of 2:1 and 3:1 were evaluated. The highest yields were obtained for choline chloride/formic acid, choline chloride/formamide. The quadratic regression equation was set up as a predictive model with an R2 value of 0.85. The optimum condition was 6 h, 40 °C, and ratio 12.5% (w/v). No significant difference was found between the predicted and experimental yield. The main components of the oleoresin were β-pinene (40.27%), cylcofenchen (11.93%) and α-pinene (7.53%) as characterized by gas chromatography-mass spectrometry. The chemical structure study by spectroscopy showed that no solvents remained in the oleoresin. Therefore, F. gummosa oleoresin can be explored as a novel promising natural pharmaceutical ingredient extracted with eutectic solvents.

Ferula gummosa Boiss. is a native indigenous plant that grows mainly in the north and southwest of Iran and is known as "Barijeh".Ferula essential oil is an effective ingredient in treatment of diarrhea and has spasmolytic activity 1 .Due to its antispasmodic, expectorant, anticonvulsant, anti-catarrhal, anti-nociceptive, and antimicrobial properties, it is also used in the adhesive, paint, cosmetic, perfume, and pharmaceutical industries 2,3 .It has also been used as flavoring agent or emulsifier in foods and beverages 4,5 .The chemical compositions of all essential oils of F. gummosa have shown that it contains mainly components of α-and β-pinene [6][7][8][9] .Several studies of the volatile components of F. gummosa oleoresin have shown that the major fraction consists of monoterpene hydrocarbons (~ 73-88%) 2,10 .
Environmentally friendly extraction technology from natural products is essential and is a high priority in chemistry due to the environmental sustainability and reduced negative impact on humans 11,12 .Natural deep eutectic solvents (NADES) are green solvents that can be used in many chemical processes, such as extraction 13,14 .NADES has many advantages over conventional solvents, such as safety, economy, inflammability, biodegradability, thermal and chemical stability, environmental friendliness, and food grade, and no further procedures are required to isolate and purify the compound 15,16 .DES is a mixture of two or more solvents that, in a given composition, exhibit a strong reduction in melting point and become liquid at ambient temperature 17 .NADES components are primary metabolites such as sugars (glucose, sucrose, and fructose), organic acids (lactic, malic, and citric acids), urea, and quaternary ammonium salts (choline chloride or betaine) [18][19][20][21] .NADES has been applied in the extraction of widely used metabolites, polar and nonpolar compounds, from plants such as green coffee beans 12,22,23 , Pseudowintera colorata 24 , Ginkgo biloba 25 , mangosteen pericarp 26 , cannabis 27 , Olea europaea 28 , Citrus grandis L. Osbeck 29 , Herba Artemisiae Scopariae 30 , Cynara cardunculus L. 31 , Carthamus tinctorius 32 , Tartary www.nature.com/scientificreports/buckwheat hull 33 , fig leaves 34 , Sophora japonica 35 , and Cajanus cajan 36 .Furthermore, the deacidification of palm and soybean oil to preserve antioxidants have been accomplished by NADES 37,38 .
The selection of the right solvents is critical to the success of the extraction process 23 .Choline chloride as a hydrogen bond acceptor (HBA) has received considerable attention due to its innocuousness, biocompatibility, and biodegradability.Furthermore, various hydrogen bond donors (HBD) such as urea, acetic acid, lactic acid, formic acid, and glycerol have also been used together with HBD to develop a eutectic solvent 39 .In this regard, the mixture of choline chloride with other solvents has been shown to improve the extraction yield compared to other conventional extractions 18,36 .Moreover, the mixture of choline chloride-lactic acid has shown higher extraction efficiency in the extraction of flavonoids from Radix Scutellariae 40 and choline chloride-alcohol for the extraction of phenolic acids 41 .
Despite numerous research works on DES extraction of different bioactive compounds such as anthocyanins, polyphenols, flavonoids, and catechins 42 , to the best of our knowledge, no research work has been conducted on the extraction of oleoresin from F. gummosa using eutectic solvents.Therefore, the aim of the present work is to apply tailored DES to extract the oleoresin from F. gummosa, compare it with the conventional method, and characterize the chemical and structure of F. gummosa oleoresin by GC and FTIR spectroscopy.

Materials
Ferula gummosa was collected in September 2021 at an elevation of approximately 2340 m above sea level in Farouj, Iran.A voucher specimen was sent to the herbarium of the Research Institute of Forests and Rangelands (TARI), Tehran, Iran.According to the previous work of Jalali 2,6,9,10 , the soil around the lower part of the herb stem was removed and the surface near the root was scratched to collect the exudate in a stainless steel container for one week.The oleoresin was obtained from approximately 100 healthy, 4-6 year old plants.The exudates were mixed and stored in a hermetic plastic container in a refrigerator (4 °C).The eutectic solvents, including analytical grade choline chloride (CL) (99% purity) as HBA and urea (purity 99-100%), lactic acid (purity > 85%), acetic acid (purity > 85%), formic acid (purity > 85%), and glycerol (purity > 90%), as HBDs, were obtained from Merck (Darmstadt, Germany).

DES synthesis
DES were prepared following previous work with some modifications 43 .The solutions of CL-urea, CL-acetic acid, CL-lactic acid, and CL-glycerol were prepared in a 2:1 ratio, while CL-formic acid and CL-formamide were obtained in a 3:1 ratio.The molar ratio of DES was chosen based on our preliminary experiments.Choline chloride (HBA) and other solvents such as HBD were mixed separately in an optimized molar ratio in a reactor and heated at 90 °C for 4 h via a magnet stirrer at 400 rpm until a transparent liquid was obtained.The DES was stored in a glass with a screw cap vial in the dark condition.The DES samples were examined intermittently over a period of several weeks for the appearance of the crystals.

Conventional extraction
The oleoresin was extracted from the exudates by hydrodistillation in a Clevenger-type apparatus for 3 h and recovered with diethyl ether, dried over anhydrous sodium sulphate, and the solvent carefully removed with a rotary evaporator.The oleoresin was stored at 4 °C until further analysis.
Extraction of F. gummosa oleoresin using DES All 6 DES mixtures were used separately to extract F. gummosa oleoresin.In a typical experiment, 1000 mg of F. gummosa exudates were taken separately in 5 g DES in a beaker.The mixture was then heated at 30-50 °C for various periods of time (3, 6, and 9 h).The mixtures were centrifuged at 3000 rpm for 5 min and the mass obtained at the bottom of the centrifuge tube was separated, washed several times with IPA and dried under vacuum.In a control experiment, 500 mg of F. gummosa exudates were mixed with 10 ml of ethanol, keeping the temperature and time the same as in the previous reactions.The supernatant was precipitated in IPA (1:3 v/v).The precipitated oleoresin dried in vacuum.

Experimental design
The effects of three factors, time (3-9 h), temperature (30-50 °C), and DES percentage (5-20% w/v), on the extraction of F. gummosa oleoresin were analyzed using a single-factor experimental design performed with Design-Expert version 13 (Statease Inc., Minneapolis, MN, USA).In each experiment, a single factor was changed while the other factors remained constant.A three-level Box-Behnken design (BBD) was further optimized.Time (A), temperature (B), and DES content (C) were independent variables, while oleoresin extraction yield was the dependent variable.The values of the variables are given in Table 1.

GC-MS analysis
Gas chromatography-mass spectrometry analysis (GC-MS) was performed using a gas chromatograph (7890B, Agilent, Santa Clara, CA, USA) connected to a mass detector (5977A, Agilent technologies, USA).The gas chromatograph was equipped with a HP 5-ms capillary column (phenylmethylsiloxane, 30 m length, 0.25 mm inner diameter, and 0.25 µm film thickness, Agilent technologies).The injector temperature was 270 °C and the oven temperature was programmed from 60 (0 min) to 200 °C at a rate of 5 °C /min.Helium was selected as the carrier gas, while the flow rate was set to 1 mL/min with an injection volume of 1 µL.However, the mass spectrometer was set to ionization mode with a voltage of 70 eV.The interfacial temperature was adjusted to 280 °C, and the mass range was between 35 and 500 m/z.The oleoresin components were identified based on a comparison of their retention indices (C 7 to C 20 n-alkanes) and their mass spectral fragmentation patterns, which were used to calculate the Kovats Indices from the gas chromatographic analysis.The Kovats indices were calculated using the Kovates equations 44,45 .

FTIR spectroscopy
FTIR spectroscopy of oleoresin was performed using a Perkin-Elmer FTIR spectrometer (Spectrum GX, USA).The measurement was performed according to the previous work with some modifications 46,47 .Briefly, the samples were on the KBr disc to scan the spectral range from 400 to 4000 cm −1 and 50 scans were measured with a resolution of 1 cm −1 .

Statistical analysis
Data were expressed as mean ± standard deviation.Experimental design and regression analysis were done with Design-Expert software version 13 (Statease Inc., Minneapolis, MN, USA), which was used for response surface method.Statistical analysis was conducted using Origin 8.0 and statistical significance was set at p < 0.05.The fitness of the proposed model was calculated by evaluating the coefficient of determination (R 2 ), lack of fit, and F-value based on analysis of variance (ANOVA).

Selection of DES solvents
As mentioned earlier, DES solvents have been widely used to extract antioxidants, polyphenols, and many polar compounds from plant materials.In the present work, 6 types of DES were used for the extraction of F. gummosa oleoresin (Table 2).The results showed that the oleoresin yield was highest when the solvents CL/formic acid and CL/formamide were used.This might be due to the stronger hydrogen bonding of DES solvents developed by choline chloride and formic acid or formamide, which can form a strong intermolecular bond with oleoresin in Ferula gummosa 48 .However, the yield increased when the molar ratio of the solvents CL/formic acid and CL/ formamide was changed from 2:1 to 3:1.These results may be attributed to acid hydrolysis of the bonds between cell wall components due to the acidifying properties of formic acid in DES (i.e., a higher amount of formic acid was present in DES), which released the oleoresin into the extraction medium and increased the extraction yield, consistent with previous results 49 .Increasing the molar ratio of choline chloride (data not shown) was also decreased the extraction rate of oleoresin, which is due to the fact that the increasing ratio of CL increases the pH of the solvent and thus affects the oleoresin yield.Therefore, the highest extraction amount of oleoresin was achieved with the ratios shown in Table 2. Consequently, the solvent DES, consisting of CL/formic acid and CL/ formamide, was selected for further research.www.nature.com/scientificreports/

Single factor experiments
The effects of time, temperature, and the ratio of DES solvent to F. gummosa on the extraction yield of F. gummosa are shown in Fig. 1.As can be seen, the extraction yield increased with time and temperature; however, higher temperatures and longer periods were not used due to the changes in the physical properties of DES.The results showed that the oleoresin yield increased significantly (P < 0.05) when the extraction time was increased from 3.0 to 9.0 h, and the maximum yield (18.67%) was obtained at 9 h (Fig. 1a).This is due to the time (3 h) required to completely release the oleoresin into the extraction medium.During the process, the solution of the DES penetrated the F. gummosa matrices, hydrolyzed the bonds between the solubilized cell walls, and then diffused out of the cell wall.However, when the extraction time was increased up to 10 h, no further increase in yield was observed.When the temperature was increased from 30 to 50 °C, the yield increased significantly (P < 0.05)  www.nature.com/scientificreports/from 9.2 to 17.4% (Fig. 1b).This increasing trend in yield with increasing temperature was related to the increase in flowability of oleoresin.The growth in temperature decreased the viscosity reduction the diffusion of DES, which assisted to destroy the chemical bond between oleoresin and cell walls, thus improving the dissolution of oleoresin.During extraction, heat enhanced the rupture of bonds between cell walls.The result is similar to previous findings in the extraction of polysaccharides from black mulberry, where the yield increased when the extraction temperature was raised from 60 to 90 °C50 .As shown in Fig. 1c, the extraction yield decreased when the ratio of DES solvent to F. gummosa exceeded 20%, while at 10%, the oleoresin yield reached a maximum of 18%.Based on the results of the one single-factor experiment, the ranges of response surface methodology (RSM) factors were applied for further evaluation as follows: DES ratio of solvent to F. gummosa (5-20%), extraction temperature (30-50 °C), and extraction time (3-9 h).

Response surface optimization
The response surface method was generally applied in optimization of the extraction parameters.Thus, the optimal conditions for extraction of oleoresin were evaluated using RSM.The experimental matrix design with a total of 34 runs (17 runs for formic acid and formamide) between the experimental and predicted values.The results showed good agreement between experimental and predicted values.As shown in Table 3, multiple regression models were fitted to the data using Design-Expert software and the quadratic polynomial regression equation corresponding to the coding independent variable of extraction yield (Y) of the samples was obtained: Y = 18.0 The analysis of variance of the quadratic response surface model is presented in Table 4.The mathematical model P < 0.05 revealed the significance of the model.The primary terms B, B 2 , and C 2 had a significant effect on Table 3. BBD with actual and predicted values for the yield of the F. gummosa Oleoresin.www.nature.com/scientificreports/ the quadratic terms; the misfit had no significant effect (P < 0.05).Model R 2 = 0.8542, R 2 Adj = 0.6562, indicating that 85.42% of the change in response value came from the selected variables.This equation fitted well with the actual situation, and could be used to analyze the test results instead of the real test 51 .The coefficient of variation (CV) represents the ratio of the standard deviation to the mean, which shows the extent of variability in data.The model's coefficient of variation was 5.67%, indicating that the experiment had good accuracy and high reliability.
According to the F value, the order of influence of each factor on the extraction rate of oleoresin was: extraction temperature > solvent ratio > extraction time.
The interaction between two independent variables could be directly observed from the 3D response surface map and the 2D response contour map.When the contour line was saddle-shaped or oval, the interaction between two factors was significant, whereas when the contour line was round, the interaction was not significant.The higher the slope of response surface, the more significant the interaction was.The effects of extraction time (A), extraction temperature (B), and solvent content (C) on the extraction rate of oleoresin from F. gummosa are shown in Fig. 2. The slope of the response surface of B was the highest, indicating that B had the greatest effect on the extraction rate.The contour lines of AB and AC were elliptical, indicating significant interactions between extraction time and extraction temperature as well as reagent content.The contour line of BC was nearly a circular circle, indicating little interaction between extraction temperature and reagent content.All these results were in agreement with the results of ANOVA in Table 4.
When the fitting model was checked (Fig. 3), the predicted value and the simulated value were almost on a straight line, indicating that the model was in good agreement with the actual situation.Considering the extraction efficiency, energy saving effect and feasibility of the experiment, the optimal extraction parameters were proposed as follows: extraction time 6 h, temperature 40 °C, solvent content 12.5% (w/v), and the predicted yield was 17.24%.Independent t-test analysis showed that there was no significant difference between the experimental value (17.24%) and the predicted value (p > 0.05).Therefore, the model could be used to optimize the extraction process of oleoresin from F. gummosa in a certain range, and the optimal extraction parameters were determined.

Ferula gummosa oleoresin characterization by FTIR
The infrared spectrum of the oleoresin obtained under optimal extraction conditions is illustrated in Fig. 5. FTIR was performed to determine structural information about the functional groups and the orientation of the groups in the oleoresin of F. gummosa.A broad band in the region of 3400 cm −1 indicates stretching vibration of the hydroxyl group (O-H).A similar result was reported for the FTIR spectrum of pomegranate peel pectin 54 .
The absorption peak of the O-H stretching vibrations was found at 3449 cm −1 .The small band was found in the 2806 cm −1 region represents the C-H absorption, which includes CH, CH 2 and CH 3 stretching and bending vibrations.It was reported that the peaks at 2800 cm −1 and 1320 cm −1 are associated with the CH stretching vibrations of CH 2 groups 55 .Moreover, a small band at around 2637 cm −1 for F. gummosa oleoresin indicates the existence of an aliphatic C-H bond.This C-H absorption and aliphatic C-H bond could be due to the methyl ester groups (OCH 3 ).In addition, two peaks were indicated between 1678 and 1460 cm −1 , characteristics of    54 .In addition, strong absorption peaks between 1149 and 1018 cm −1 indicated the potential presence of pyranose ring.The absorption bands in the range of 1300 cm −1 and 1450 cm −1 for both high and low methyl groups were assigned to CH 3 corresponding to asymmetric stretching vibrations, while the bands around 1020 cm −1 were allocated to the saccharide structure (C-O-C).IR spectrum for DES (Fig. 5) was also determined in order to check any DES residues in the extracted oleoresin.The result showed that the IR spectra of DES were not found in the IR spectra of oleoresin.Therefore, it was assumed that no DES residue was detectable in the oleoresin.

Conclusion
Summary of the work was provided in a separate figure (Fig. 6).The use of DES proved to be an effective extraction medium for the extraction of oleoresin from F. gummosa.Optimal conditions were determined: to be a solvent ratio 12.5% (w/v), an extraction temperature of 40 °C, extraction time of 6 h, and a molar ratio of the components of DES 3:1.Among the different DES solvents, choline chloride/formic acid and choline chloride/ formamide in the ratio of 3:1 showed the highest yield of in the extraction of oleoresin.The optimum condition gave a yield of oleoresin of 17.24%.The GC-MS analysis also exhibited its potential as a functional oleoresin that could be used as a food ingredient.The oleoresin was mainly composed of β-pinene (40.27%), cylcofenchen  www.nature.com/scientificreports/(11.93%) and α-pinene (7.53%), with the main group being monoterpene hydrocarbons (62.70%), including α-thujene, α-pinene and β-pinene potential.These findings demonstrated that oleoresin extracted from F. gummosa could be explored as novel promising natural pharmaceutical ingredient that can be extracted with eutectic solvents.

Natural deep eutectic solvent preparation
Extraction tim e, h

Figure 3 .
Figure 3. Actual value and predicted value fitting model diagram.
www.nature.com/scientificreports/an ester carbonyl group (C=O) and an unesterified carboxylate ion (COO-), respectively56 .The prominent bands were also detected in the range between 1300 and 1000 cm −1 attributed to C-O stretching in C-O-C and C-O-H.These could be expected owing to the typical presence of pyranose ring in oleoresin, which linked with one another via C-O-C glycosidic linkages and side group C-O-H linked bonds.As mentioned by Yang et al. (2018), the significant absorption band located at 1677 cm −1 was assigned to C-O stretching vibration of ester carbonyl, and the absorption band at around 1677 cm −1 was represented by C-O stretching vibration of non-methylated carboxyl in pomegranate peel pectin

aFigure 6 .
Figure 6.Summary of the research work on natural deep eutectic solvent (NADES) in oleoresin extraction from Ferula gummosa roots.

Table 1 .
Coded and actual levels of independent variables used in BBD design.

Table 4 .
Response surface quadratic model variance analysis.

Table 5 .
Chemical composition of F. Gummosa oleoresin by GC-MS analysis.