Document Type : Research Paper
Authors
1 Department of Food Technology Research, National Nutrition and Food Technology Research Institute, Faculty of Nutrition Sciences and Food Technology, Shahid Beheshti University of Medical Sciences, Tehran, Iran
2 Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Tehran University of Medical Sciences, Tehran, Iran
Abstract
Keywords
INTRODUCTION
Geraniol (GR) is a monoterpene alcohol with pharmacological properties that is derived from the essential oils of aromatic plants such as Kansho-shochu, sweet basil, hopped beer and green tea [1–3]. The US Food and Drug Administration (FDA) has certified geraniol as a food additive for flavoring foods such as beverages, candies and ice creams [4]. Promising bioactive substances such as essential oils (EO) are used as preservatives in foods or medicines due to the various functional activities of these compounds, including antioxidant, antimicrobial and anticancer effects [5, 6]. Biological activities of the EOs against foodborne pathogens such as Staphylococcus aureus, Bacillus cereus, Salmonella enteritidis, Escherichia coli, Candida albicans and Aspergillus niger have extensively been investigated [2, 7]. Encapsulation simultaneously increases antimicrobial potencies of the aromatic oily components of EOs such as geraniol and their certainly controlled/sustained releases, facilitating their close interactions with microorganisms. Nanobiomaterials used in designation of formulations include various advantages and important functional properties, which are achieved by increased surface areas and binding of the bioactive molecules, enhanced permeabilities and retention effects [8–10]. When encapsulated in β-cyclodextrin (β-CD), EOs and their components present high stabilities against degradative actions of oxygen and temperature as well as promising various uses in pharmaceutical and food industries. Moreover, better retention and increased protection of EOs in colloidal matrices depend on the type of the technique [11, 12]. In the last few decades, various formulation approaches have been introduced for the development of β-cyclodextrin with various EOs compounds [11–14]. For poorly soluble compounds such as geraniol, β-CD is a useful material for complexation because of its various advantages [15, 16]. Several studies have reported that geraniol includes antimicrobial activities against foodborne pathogens. A summary of these studies is provided in Table 1. These studies have majorly included pure forms of GR instead of GR/CD molecular-inclusion complexes [17–28].
As previously reported by Menezes et al., and Hadian et al., GR can be loaded in β-cyclodextrin using precipitation in nanoparticle (NP) complex forms [29, 30]. Use of systemic experimental designs with a mathematical optimization approach such as RSM helps determine acceptable formulations of flavor compounds [31]. Advances of using experimental design methods include decreases in number of experiments, identification of interactions between multiple factors and detection of optimal responses within the experimental regions by decreases in variability [32–34]. Development of encapsulated holy basil essential oil in gelatin using RSM with a maximum encapsulation efficiency has been reported by Sutaphanit et al. [35]. Alves et al., have successfully optimized encapsulation of the cocoa waste extracts using chitosan, maltodextrin and RSM as a nutrient source for aquatic animal feeds [33]. Nowadays, use of natural preservatives in food industries is a major concern of the global public health. Thus, a novel plan for the development of effective carriers of bioactive compounds with various advantages can be the use of nanotechnology for the inclusion complex formation with β-CDs. Based on the previous reports by the authors [30], simultaneous use of β-CD complexation based nanocarriers for the development of natural antimicrobial delivery systems has been investigated in this study. The objective of the current study was to investigate encapsulation efficiency and loading of GR in β-CDs molecular-inclusion complexes using RSM. Furthermore, in vitro antimicrobial activities of GR inclusion complexes against S. aureus, B. cereus, S. enteritidis, E. coli, C. albicans and A. niger were studied. The present study is the first study to report antimicrobial activities of the GR inclusion complexes in nanosize formulations.
MATERIALS AND METHODS
Materials
Geraniol (98% purity) and β-cyclodextrin (≥ 97% purity), menthol and tetramethylsilane (TMS) were purchased from Sigma-Aldrich, Steinheim, Germany. Anhydrous sodium hydroxide, hydrochloric acid and dimethyl sulfoxide-d6 (DMSO-d6, at least 99.5% purity) were supplied by Merck, Darmstadt, Germany. Demineralized water (pH 7.6 ±0.2) was used in the study. Other solvents included HPLC grades. Staphylococcus aureus PTCC 1112, Bacillus cereus PTCC 1015, Salmonella enteritidis PTCC 1709, Escherichia coli PTCC 1330, Candida albicans PTCC 5027 and Aspergillus niger PTCC 5012 were provided by the Iranian Research Organization for Science and Technology (IROST), Tehran, Iran.
Preparation and optimization of the GR-β-CD inclusion complexes
The GR-β-CD NPs were prepared using precipitation method with minor modifications according to Rachmawati et al. [36]. Briefly, 500 mg of the β-CD were dissolved in 20 mL of ethanol and distilled water (DW) (20:80 v/v) for 30 min at 45 °C. Then, GR was added to the hydro-alcoholic solution of β-CD with continuous agitation. The mixture was stirred at 300 rpm for 1 h at 37 ºC and set for 4 h at room temperature. After sonication, solution was centrifuged at 3000 rpm for 15 min to precipitate free GR. The supernatant was filtered using 0.2-μm filters. Then, sample was transferred to a vacuum chamber at 40 °C for 6 h. In this study, RSM and rotatable central composite design (RCCD) with two independent variables were used to optimize formulations of the GR-β-CD complexes [37]. Using multivariate quadratic equation, main and reciprocal effects of the independent variables on each dependent variable were separately assessed. The experimental design and data analysis were carried out using Design-Expert Software v.8.0.5 (Stat-Ease, Minneapolis, USA). All experiments were carried out in triplicate. Independent variables included β-CD (X1) and GR (X2) and dependent variables included loading (Y1) and encapsulation (Y2) efficiencies. Design of the statistical test was carried out to determine optimal conditions for the achievement of the highest percentage of dependent variables. Five coded levels of the two variables were incorporated into the design, resulting in 13 experiments. Based on the number of independent variables and their levels, number of the test samples included 13 treatments with five central and eight non-centered points. Equations 1 and 2 show effects of the variables on the percentage of loading and GR encapsulation, respectively. In analysis of variance (ANOVA), the second power of independent variable B2 (GR2) was significant at 0.01. The regression model was as follows:
Y1 = +32.76 + (-0.69)X1 + (- 0.15)X2 + 4.8X1X2 + 4.48X12 + (- 0.01)X22 (1)
For encapsulation efficiency, the quadratic model was significant at 0.01 (< 0.0002). In ANOVA, the second power of independent variable B2 (GR2) was significant at 0.01. The regression model was as follows:
Y2 = +476.39 + (-9. 93)X1+ (- 7.62)X2 + 0.14X1X2 + 0.056X12 + (- 0.19)X22 (2)
Coefficient of determination (R2) and adjusted and predicted coefficients of determination (R2-adj and R2-pred) were used to compare various models and select an appropriate model for the response. Based on the preliminary tests, ranges of the changes were assessed at 5–23.11% by weight in GR and 76.89–95% by weight in β-CD.
Characterization of GR-β-CD complexes
Assessment of the geraniol EE and loading
Total quantity of the GR in GR-β-CD inclusion complex powders and surface oil extractions was measured using gas chromatographer (Agilent, Waldbrornn, Germany) Hewlett-Packard Model 6890 equipped with an FID detector and stationary phase with HP-5 Column (30 m length × 0.25 mm ID × 0.25 µm film thickness), according to Bhandari et al. [38, 39]. Concentration of the GR was quantified using standard calibration curves. Calibration curves were plotted for each compound using six various concentrations of GR and the internal standard in hexane in the range of 0.5–10 μg mL-1. Loading and EE of the GR-β-CD complexes for GR were respectively calculated using Equations 1 and 2 according to Tao et al. [40]:
EE (%) = A (mg) / B (mg) × 100 (1)
Loading (%) = A (mg) / C (mg) ×100 (2)
Where, A, B and C were quantities of the trapped GR (differences between the total oil and surface oil) in the complexes, primary feeding quantities of the GR and number of the produced particles, respectively. Calibration curves were plotted for each compound using six various concentrations of geraniol and the internal standard in hexane in the range of 0.5–10 μg mL-1.
Particle size distribution, Zeta-potential and polydispersity index assessment
Particle size and polydispersity index (PDI) of the optimized GR-β-CD inclusion complexes were assessed using photon-correlation spectroscopy (PCS). Suspension was serially diluted in MiliQ water and analyzed using Malvern ZetaSizer Nano Series (Malvern, UK) [36].
Infrared spectroscopy
Geraniol, β-CD, physical mixture and nanoinclusion complexes of the GR were studied using infrared (IR) spectroscopy. Geraniol was scanned in IR spectra of 4,000–400 cm-1 (Spectrum One, Perkin Elmer, Waltham, MA, USA). Spectra included an average of 32 scans at a resolution of 4 cm-1. The IR spectral analysis was carried out using KBr pellet press method and the IR spectra were recorded using Nicolet Magna-IR System 550 equipped with Nicolet OMNIC Software [41].
Scanning electron microscopy
Hitachi SU 3500 Scanning Electron Microscopy (SEM) (Hitachi, Japan) instrument was used for the surface morphology. Samples were attached to an aluminum stub using double-sided carbon tape and were electrically conductive by coating with thin layers of gold (33 Aº thickness) prior to study with SEM.
Differential scanning calorimetry
Thermal analysis of β-CD, geraniol and their physical mixture and inclusion complex was carried out using Mettler Toledo DSC System (DSC, 823E, Mettler Toledo, Switzerland). Mettler Stare Software v.9.x was used for data acquisition and indium for calibrating the instrument. Sample (2 mg) was sealed in the aluminum pans and heated at the rate of 10 °C min-1 from -20 to 300 °C under 8-kPa nitrogen atmosphere [36]. Physical mixtures of the GR and β-CD were prepared using simple mechanical mixing. The mole ratio of GR to HP-β-CD was 1:2.
Antimicrobial activity assay
Bacteria were cultured on soybean-casein digest agar and incubated at 30–35 oC for 20–24 h. A loopful of bacteria was suspended in sterile normal saline and homogenized using vortex. Density of the microbial suspension was adjusted to 0.5 McFarland turbidity standard. The 0.5 McFarland standard included an absorbance reading of 0.08–0.1 at 625 nm [42]. Furthermore, C. albicans was cultured on Sabouraud dextrose agar (SDA) and incubated at 35 oC for 24 h. Then, a suspension of 1-5 × 106 cells mL-1 was prepared using 0.5 McFarland turbidity standard (light transmittance 76% at 530 nm). The A. niger was cultured on SDA and incubated at 20–25 oC until spore formation (72 h). The spore suspension was prepared by washing surface of the cultured slant using normal saline containing 0.05% Tween 80 and sterile glass beads. Suspension was transferred to a sterile tube. After 3-5 min, rough particles were precipitated and the supernatant was transferred to another sterile tube. Turbidity of the spore suspension was adjusted to transmittance of 80-82% at 530 nm (equal to 4 × 105 to 5 × 106 spore mL-1) [43]. Minimum inhibitory concentrations (MICs) of the GR and GR-β-CD against the bacteria were assessed using agar dilution method [44].
Statistical analysis
In this study, all experiments were carried out in triplicate and mean ±SD (standard deviation) was reported. The ANOVA, regression coefficient calculation and performance stepwise procedure were used to simplify the models and generating of three-dimensional surface plots using Design-Expert Software v.7.1.1 (Stat-Ease, Minneapolis, USA). Significance of the equation parameters for each response was reported using F-value at a probability (P) of 0.05. In addition, adequacy of the models was calculated using model analysis, lack-of-fit test and coefficient of determination (R2) analysis.
RESULTS AND DISCUSSION
Characterization of GR-β-CD inclusion complexes
To overcome problems such as time and cost linked to the formulation, effective and reliable RSM was used to assess consequences of the independent variables [β-CD (X1) and GR (X2)] on the responses [loading (Y1) and encapsulation efficiency (Y2)] of the NP formulation. Effects of the proportions of GR and β-CD on the encapsulation efficiency and loading of the guest molecules (GR) were the major objectives of the first part of this study. Correlation coefficient of the model for the loading and encapsulation efficiency included 0.83 and 0.88, respectively. The P-value for the loading and encapsulation proportions shows appropriately and efficiency of the computational model; hence, it can be used as a completely successful model for the calculation and estimation of the proportions of loading and encapsulating. Using quadratic equation of multivariate, main and reciprocal effects of the independent variables were tested separately for each dependent variable. Significance of each term was reported using F-value and P-value; where, small P-values and large F-values illustrate further significant effects on the respective independent variables (Table 2). For loading and EE proportions, the quadratic model was significant at 0.05 (< 0.0001).
The EE% of GR-β-CD formulations was calculated using various GR to β-CD ratios. Total oil and surface oil were detected using GC/FID. Based on the results of early stages and available resources, effects of the two independent variables, β-CD and GR were considered as 80–95 and 5–20% w/w, respectively (Table 1). Three-dimensional response surface plots of the β-CD on encapsulation of GR are shown in Fig. 1a. As stated, effects of the independent variables (GR and β-CD) on loading and EE proportions showed a descending trend to the central point and then changed.
As seen in Fig. 1, GR included obvious effects on loading and encapsulation efficiencies of the NP formulation. Based on the model, values of loading and encapsulation efficiency respectively included 10.45 and 87.25% under optimal conditions. To investigate effects of independent variables on dependent ones, surface response plots of the quadratic polynomial models were created by varying two of the independent variables within the experimental ranges (Fig. 1b). As Figs. 1a and 2b illustrate, loading and encapsulation efficiencies increased by increasing the weight proportion of β-CD to GR; therefore, interaction of these two variables was significant. Previously, several attempts have been made to complexation of constituent of several EOs with cyclodextrin [12, 14, 29, 30, 41]. The highest responses were achieved for F3 formulation at 95 and 12.50% for β-CD and GR, respectively. Furthermore, the lowest loading (28.5%) and encapsulation (2.2%) efficiencies were achieved at 1.89:87.5 GR-β-CD ratio for F4 formulation. GR, as a primary component of the rose, palmarosa and citronella oils is able to increase transdermal penetration due to its lipophilic nature [14, 29]. Zhang et al. reported the formation of inclusion compounds in terms of baicalin-ß-CD ratio, time and temperature to achieve the maximum encapsulation efficiency using RSM [45]. A similar report was published by Menezes et al., demonstrating that the EE of GR encapsulated in β-CD was nearly 60% at the presence of surfactants [29]. It has been concluded that differences in proportion of EE may be due to the method of formulation and the solvent used for extraction. Based on the recent studies, optimized nano GR-β-CD inclusion complexes using RSM have shown a higher encapsulation efficiency, compared to the previous studies [30]. The current results are partly similar to the results by other researchers. As shown in Figs. 1c and 1d, the mean particle size and zeta-potential of optimized nanoinclusion complexes (F3 formulation) was in the range of 117 nm ±1 and -17.7 mV ±2.3, respectively. Furthermore, NPs were uniform and monodispersed and PDI was in the range of 0.14 ±0.03.
Scanning electron microscopy and infrared spectroscopy
Fig. 2 demonstrates images of the pure β-CD, GR-β-CD physical mixture and F3 formulation. Pure β-CD showed crystalline particles with various sizes with no characteristic shapes (Fig. 2a). As illustrated in Fig. 2c, particle shape and morphology of F3 formula exhibited various forms, compared to that of physical mixture of GR-β-CD (Fig. 2b). Pellets of the GR-β-CD included a uniform porous surface (Fig. 2c).
The IR spectra for GR (Fig. 2d) showed absorption bands at 3354 (O-H stretching vibration), 2968, 2919 and 2729 (C-H stretching vibration), 1668 and 1444 (C=C stretching vibration), 1235 (C-C stretching vibration) and 1118, 1098 and 1001 (C-O stretching vibration) cm-1. The IR spectra for β-CD (Fig. 2e) included absorption bands at 3582 (O-H stretching vibration), 2914 (C-H stretching vibration), 1652 (H-O-H stretching vibration), 1159 (C-O stretching vibration) and 1091 (C-O-C stretching vibration) cm-1. As shown in Fig. 2f, changes were detected in shape, position and intensity of the absorption bands of pure β-CD, the physical mixture and the inclusion complexes. The physical mixture (Fig. 2f) included a combination of the individual patterns of β-CD and GR (F3 formulation). In IR spectrum the peak at 3000 cm-1 corresponding to the C-H group stretching [12, 29]. Although obvious changes were seen in spectra of the inclusion complex of GR and β-CD (Fig. 2g), absorption peaks shifted to lower frequencies at 3389, 2924, 1645 and 1416 cm-1. These changes could demonstrate formation of the intramolecular hydrogen bonds in GR-β-CD inclusion complex. When the functional groups of the drug were included within the apolar cavity of β-CD in the complex an apparent decrease of the absorption intensities of the corresponding bands is appeared [32]. This results are similar to other findings [12, 28]. In solid state complexes, the only host-guest interactions show are O H· · ·O hydrogen bonds between the hydroxyl groups of geraniol and cyclodextrins [16]. As previously reported by Chadha et al., the present IR results have demonstrated that the inclusion complex prepared using coprecipitation method changed overall intensity of the IR absorption patterns of GR complex with β-CD and presence of interactions between GR and β-CD (46). The SEM micrographs have verified interactions between the matrix and the guest in the inclusion complexes. A similar phenomenon has been reported in literatures [12, 13, 29, 30].
Differential scanning calorimetry
In fact, DSCs are used for the investigation of stability and physical and chemical interactions between β-cyclodextrin and guest molecules. If the guest molecules are well placed in polymer matrices, no exothermic or endothermic peak should be observed in guest-polymer or bioactive-polymer thermograms [9, 45]. The DSC thermograms of samples are presented in Fig. 3. The guest substance might change by melting, evaporation, decomposition, oxidation or polymorphic transition. It has generally been reported that disappearance of endothermic peaks, appearance of new peaks and peak broadening or shifting to various temperatures indicate transition in crystal-lattice, melting, boiling and sublimation points of the inclusion complexes formation. In this study, DSC thermograms of geraniol (Fig. 3a) showed endothermic melting peaks at 248 and 155 °C, possibly due to the oxidation or elimination of water [9, 47].
Thermograms of β-CD and geraniol physical mixture (Fig. 3b) have shown two endothermic peaks at nearly 130 and 240 °C, which were close to those of pure β-CD and geraniol. These may indicate less or no interactions between the drug and β-CD in the physical mixture. In contrast, DSC curves for F1 formulation (Fig. 3c) have shown no endothermic peaks, which suggest that the guest was incorporated into the NPs and verify formation of a host-guest inclusion complex. As reported by the previous studies, elimination of the free-geraniol and β-cyclodextrin peaks demonstrates formation of a complex between these two compounds [12, 41, 47].
Antimicrobial activity of the GR-β-CD inclusion complexes
In this study, antimicrobial properties of F3 formulation against two Gram-positive bacteria of S. aureus and B. cereus and two Gram-negative bacteria of S. enteritidis and E. coli were investigated. As illustrated in Fig. 4, diameters of the GR inhibition zones were significantly larger than those of the DMSO inhibition zones, indicating that antibacterial effects of GR against Gram-positive and Gram-negative bacteria were significantly greater than those of DMSO and water.
In all experiments, diameters of the inhibition zones of gentamicin showed gentamicin as an effective antibiotic against the highlighted bacteria. This has verified that the current method for the assessment antimicrobial activity was carried out correctly. Moreover, MIC of the encapsulated GR was significantly larger than that of its free form, representing a stronger antibacterial effect of the GR-β-CD inclusion complexes (P < 0.05) (Table 3). Comparisons between the antibacterial activities of free and complexed GR and the gentamicin showed significant differences (P < 0.001). This study showed that the MIC of GR for Gram-positive bacteria was mostly lower than that for Gram-negative bacteria. The MIC varied from 3 × 103 to 5 × 103 μg mL-1 for free GR, while it varied from 2 × 103 to × 103 μg mL-1 for the GR-β-CD inclusion complex.
Considering MIC values and diameters of the inhibition zones (Fig. 4 and Table 3), B. cereus and S. entritidis were the most sensitive and the most resistant bacteria to GR either in complex or free form, respectively. These results are likely linked to the fact that further resistance of Gram-negative bacteria to GR can be attributed to the presence of their outer membranes, preventing penetration of lipophilic EOs into the bacterial cells [48, 49]. The β-CD is a relatively large molecule with a number of hydrogen donors and acceptors in its chemical structure. When poorly soluble components such as essential oils are complexed with β-CD, aqueous solubility and efficacy of the poorly soluble components are improved [44, 45]. Results from the current study are similar to results from Hill et al. study, who showed that the antimicrobial efficiency of GR in β-cyclodextrin complexes was higher than that in its free forms since it increased water solubility of the lipophilic compounds and hence helped its easier penetration into the bacterial cells [13]. In another study, it was indicated that antibacterial activity of the monoterpenoid alcohols such as linalool, nerol, citronellol and GR was stronger than antifungal activity of these compounds [50]. Wang et al. reported that β-CD inclusion complexes improved antimicrobial activity of the aromatic oily liquids extracted from EOs such as eugenol through increasing its accessibility to susceptible sites such as membrane and cytoplasm of the microorganisms by enhancing EO solubility [12]. However, advantages of the β-CD encapsulation are not limited to increasing solubility of the EO components. Improvement of the antimicrobial mode of action, decrease of the minimal concentration needed for antimicrobial inhibition and protection of EOs against degradation are other important advantages [12, 41]. The greater susceptibility of Gram-positive bacteria to EOs has been reported severally in literatures [49, 50]. Results of the current study have suggested that the MIC of free GR and βCD-GR nanoinclusion complexes against the Gram-positive and Gram-negative bacteria is much higher than 100 μg mL-1 indicating its weak antimicrobial activity; as previously reported by other studies [48, 51].
In the present study, antifungal activities of GR and GR-β-CD inclusion complexes against C. albicans and A. niger were assessed using microdilution method at concentrations of 1000–6000 μg mL-1. Free GR and GR-β-CD complexes showed antifungal activities at completely different concentrations (P < 0.05). The MIC values of free GR and GR-β-CD complexes against C. albicans were seen at concentrations greater than 4000 and 3000 μg mL-1, respectively (Table 4). A similar result was observed for A. niger at concentrations greater than 4000 and 3000 μg mL-1, respectively. These MIC values were significantly greater than those of amphotericin B. As Table 3 demonstrates, diameters of the antifungal inhibition zones against C. albicans and A. niger induced by GR-β-CD complexes were greater than those induced by its free form (P < 0.05).
This may support inhibitory effects of nanocomplex and free GR against C. albicans growth as 3 × 103 and 4 × 103 μg mL-1, respectively. These values are much higher, compared to the value of 100 μg mL-1. This finding was previously reported by Ayala-Zavala et al. for the inhibitory effects of EO inclusion complexes against fungal growth [52]. Singh et al. showed that disruption of membrane and lose of its integrity were the possible action mechanisms of GR against C. albicans cells using micrographic images [53]. According to Zanetti et al. who studied the antibacterial effects of pure geraniol and cinnamic acid, GR showed a greater antimicrobial activity against a broad spectrum of pathogenic bacteria, similar to that demonstrated in the present study. They suggested this function as a principal mechanism; by which, geraniol and the nerol included double bonds in their molecules. Therefore, they showed a greater antimicrobial activity against the highlighted bacteria, compared to that the other essential oils with single bonds did [54]. Other studies support relationships of the chemical structure of EOs and associated compounds such as position of the OH group in structure of geraniol, possibly showing roles in their antimicrobial activities. Nanoinclusion complexes of geraniol included a higher antimicrobial activity than that the pure or free form of geraniol did. Nanoinclusion complex containing bioactive compounds is an efficient methodology approach to increase physical stability of the active substances, retain them from interactions with other ingredients and increase their bioactivity [47, 54–56]. In general, it seems that biomolecules in nanoparticle structures, as appropriate carrier systems, have improved their delivery characteristics such as antimicrobial activity, compared to single forms of the compounds.
CONCLUSION
In the present study, precipitation method were used for the fabrication of GR-β-CD inclusion complexes with a size distribution of 117 nm ±1. Physical characteristics of GR-β-CD inclusion complexes were verified using SEM, IR spectroscopy and DSC analysis. Response surface design has been shown to be effective for optimization of factors involved in developing of GR-β-CD inclusion complexes with high loading and encapsulation efficiencies. Furthermore, the optimized formulation has demonstrated significant antimicrobial activities against two Gram-positive and two Gram-negative bacteria, compared to free geraniol with a similar concentration. Based on the optimized GR-β-CD nanoinclusion complexes (F3 formulation), the maximum and the minimum inhibition zones of the GR-β-CD complexes have been recorded for B. cereus and E. coli, respectively. Based on these results and biological activities of the developed formulation, further in vivo studies on the use of this compound in food industries are necessary.
ACKNOWLEDGMENT
This study was carried out as a part of a research project supported by the National Nutrition and Food Technology Research Institute, Shahid Beheshti University of Medical Sciences.
CONFLICT OF INTEREST
The authors declare that there is no conflict of interests regarding the publication of this manuscript.