Comprehensive Characterization of Linalool-HP-β-Cyclodextrin Inclusion Complexes

The objective of the present study is to obtain linalool- cyclodextrin (CDs) solid complexes for possible applications in the food industry. For this purpose, a detailed study of linalool complexation was carried out at different pH values, to optimize the type of CDs and reaction medium that support the highest quantity of encapsulated linalool. Once demonstrated the ability of hydroxypropyl-β-cyclodextrin (HP-β-CDs), to form inclusion complexes with linalool (KC = 921 ± 21 L mol−1) and given their greater complexation efficacy (6.788) at neutral pH, HP-β-CDs were selected to produce solid inclusion complexes by using two different energy sources, ultrasounds and microwave irradiation, subsequently spraying the solutions obtained in the Spray Dryer. To provide scientific solidity to the experimental results, the complexes obtained were characterized by using different instrumental techniques in order to confirm the inclusion of linalool in the HP-β-CDs hydrophobic cavity. The linalool solid complexes obtained were characterized by using 1H nuclear magnetic resonance (1H-NMR) and 2D nuclear magnetic resonance (ROSEY), differential scanning calorimetry, thermogravimetry and Fourier transform infrared spectrometry. Moreover, the structure of the complex obtained were also characterized by molecular modeling.


Introduction
Nowadays, the demand for natural compounds in the food industry has increased due to the need to obtain more biodegradable preservatives [1,2]. Essential oils (EOs) have been known since the Middle Ages due to their antiseptic, therapeutic properties and their intense aroma, arousing great interest in food preservation. The difference in their organoleptic profiles is conditioned by their chemical composition and richness, characteristics that determine their antioxidant and antimicrobial properties [3], which have aroused the interest of large companies related to human and animal nutrition, pharmacy and cosmetics [4]. The antimicrobial action of the EOs is due to their ability to penetrate through the bacterial membranes into the cell, causing the inhibition of their vital functions [5][6][7][8]. EOs are very complex natural mixtures characterized by two or three main components in fairly high concentrations (20-60%) [9]. As observed in Table 1, the KC values obtained increased when the pH approaches neutrality (7.0) for all CDs studied (Table 1). This result agrees with the fact that the solubility (S0) of linalool varies with the medium pH, being more soluble at acidic or alkaline pH (S0 = 9.34 mmol L −1 ) than at neutral pH 7.0 (S0 = 7.37 mmol L −1 ) ( Table 1). It must be taken into account that linalool has a hydroxyl group at the C3 chiral carbon, which implies that the molecule behaves as a weak acid (pKa = 18.46), so its dissociation degree will be conditioned by the pH of the medium, and therefore, its solubility. The main factors limiting complex formation are the relationship between the size of internal cavity of the CDs and the host molecule, and the polarity of the host molecule. So, the protonation of linalool and ultimately its solubility is decisive in the complexation process and the stability of the complexes formed. Table 1. Recovery aqueous solubility (S0), complexation constant (Kc), complexation efficiency (CE) and molar ratio of linalool for α-, β-and HP-β-CDs at different pH values. Standard deviation of triplicate diagrams. As observed in Table 1, the K C values obtained increased when the pH approaches neutrality (7.0) for all CDs studied (Table 1). This result agrees with the fact that the solubility (S 0 ) of linalool varies with the medium pH, being more soluble at acidic or alkaline pH (S 0 = 9.34 mmol L −1 ) than at neutral pH 7.0 (S 0 = 7.37 mmol L −1 ) ( Table 1). It must be taken into account that linalool has a hydroxyl group at the C3 chiral carbon, which implies that the molecule behaves as a weak acid (pKa = 18.46), so its dissociation degree will be conditioned by the pH of the medium, and therefore, its solubility. The main factors limiting complex formation are the relationship between the size of internal cavity of the CDs and the host molecule, and the polarity of the host molecule. So, the protonation of linalool and ultimately its solubility is decisive in the complexation process and the stability of the complexes formed. As can be seen in the results presented in Table 1, independently of the CDs used, the highest K C value was obtained for the complexes formed at pH 7.0 (148 ± 12 L mol −1 for α-, 241 ± 26 L mol −1 for βand 921 ± 21 L mol −1 for HP-β-CDs), which differs significantly from the Kc values obtained at either acidic or alkaline pHs. These results were similar than those obtained for thymol complexation, which showed the highest Kc values at neutral pH as well as linalool [9].
Despite this improvement at neutral pH, in general, the K C values were low, probably conditioned by the structure and instability of linalool, since it is isomerized easily giving rise to primary alcohols (geraniol and nerol), even forming cycles, not favoring its inclusion in the hydrophobic CDs cavity. It is important to note that linalool is characterized by its chemical instability, since it is capable of oxidizing, reducing and polymerizing easily when interacts with CDs in acidic and alkaline conditions. It can suffer intramolecular transformations that can modify its structure and, consequently, modify its ability to be included in the hydrophobic cavity of CDs.
When compared K C values obtained for native and modified β-CDs, the results were 241 ± 26 L mol −1 for βand 921 ± 21 L mol −1 for HP-β-CDs at neutral pH ( Table 1). The differences observed were related with the change of H atoms by hydroxypropyl groups in the case of HP-β-CDs favoring the entrapment of linalool. The K C value obtained for HP-β-CDs (4-fold than that obtained for β-CDs), indicated the formation of a more stable complex due to the intensity of hydrophobic and Van der Waals interactions involved. In addition, the presence of hydroxypropyl groups is enough to break the hydrogen bridges close to the CDs cavity, making the entrance of linalool molecules more accessible to the apolar cavity. These results agree with those obtained for resveratrol [30], kaempferol and p-coumaric acid that present cyclic structures, or citral with a non-cyclic structure [31,32].
In summary, the stability of the soluble complexes formed between linalool and CDs depends both on the medium pH and the CDs type, being the optimal conditions pH 7.0 and HP-β-CDs.
In order to complete the solubility study, not only the stability of the complexes formed (K C ) was analyzed, but also the complexation efficiency (CE) between linalool and CDs was calculated according to Equation (2). For 1:1 stoichiometry complexes the CE is independent of the aqueous solubility of the encapsulated compound. The comparison of CE values gives more information than comparing K C values when study different types of CDs for the complexation of the same compound. As can be seen in Table 1, HP-β-CDs showed the highest CE values at all pHs studied (0.41 for pH 3.5, 0.443 for pH 5.5, 1.836 for pH 6.5, 6.788 for pH 7.0 and 5.215 for pH 8.5) ( Table 1). In the case of pHs 6.5, 7.0 and 8.5, the values obtained above 1 indicated that the quantity of CDs forming complexes with linalool where higher than free CDs in solution. Similar results were previously described for eugenol with HP-β-CDs (2.81), ferulic acid with HP-β-CDs and RAMEB (1.43 and 1.79 respectively), and coumaric acid with HP-β-CDs, RAMEB and CRYSMEB (2.30, 2.57, 1.89 respectively) [33]. The rest of CDs studied, αand β-CDs, presented CE values below 1 at most of studied pHs, indicating that most did not increase linalool solubility (Table 1).
Moreover, CE was also used to calculate the linalool:CDs molar ratio (MR) by using Equation (3), which can be correlated to the expected increase in formulation bulk. As in the case of CE, the highest value for MR, was also obtained for HP-β-CDs at pHs 6.5, 7.0 and 8.5 (1:1 in all cases) ( Table 1). This 1:1 value indicated that almost all CDs in solution were forming soluble complexes with linalool. These results agree with those obtained for eugenol and coumaric acid, in which the highest EC values coincide with an MR of 1:1 [33]. In summary, taking into account K C , CE and MR (linalool:CDs) values, HP-β-CDs at pH 7.0 are the best conditions to increase the aqueous solubility of linalool, forming the most stable and soluble complexes.

Obtaining Solid Inclusion Complexes Linalool-HP-β-CDs
In order to facilitate the handling and preservation, the formation of linalool complexes in solid state was also studied. As described previously, linalool complexes could be made by using different energy sources to reach equilibrium: microwave irradiation (MWI) or ultrasound (US). Figure 2 shows the solubility profiles of linalool in the presence of HP-β-CDs a pH 7.0, by reaching the equilibrium using both methods. As can be seen in Figure 2A, the solutions treatment with MWI (24 h MWI or 48 h MWI) to reach the equilibrium, much less increased the linalool solubility than the application of US ( Figure 2B). solubility of the encapsulated compound. The comparison of CE values gives more information than comparing KC values when study different types of CDs for the complexation of the same compound. As can be seen in Table 1, HP-β-CDs showed the highest CE values at all pHs studied (0.41 for pH 3.5, 0.443 for pH 5.5, 1.836 for pH 6.5, 6.788 for pH 7.0 and 5.215 for pH 8.5) ( Table 1). In the case of pHs 6.5, 7.0 and 8.5, the values obtained above 1 indicated that the quantity of CDs forming complexes with linalool where higher than free CDs in solution. Similar results were previously described for eugenol with HP-β-CDs (2.81), ferulic acid with HP-β-CDs and RAMEB (1.43 and 1.79 respectively), and coumaric acid with HP-β-CDs, RAMEB and CRYSMEB (2.30, 2.57, 1.89 respectively) [33]. The rest of CDs studied, α-and β-CDs, presented CE values below 1 at most of studied pHs, indicating that most did not increase linalool solubility (Table 1).
Moreover, CE was also used to calculate the linalool:CDs molar ratio (MR) by using equation 3, which can be correlated to the expected increase in formulation bulk. As in the case of CE, the highest value for MR, was also obtained for HP-β-CDs at pHs 6.5, 7.0 and 8.5 (1:1 in all cases) ( Table 1). This 1:1 value indicated that almost all CDs in solution were forming soluble complexes with linalool. These results agree with those obtained for eugenol and coumaric acid, in which the highest EC values coincide with an MR of 1:1 [33]. In summary, taking into account KC, CE and MR (linalool:CDs) values, HP-β-CDs at pH 7.0 are the best conditions to increase the aqueous solubility of linalool, forming the most stable and soluble complexes.

Obtaining Solid Inclusion Complexes Linalool-HP-β-CDs
In order to facilitate the handling and preservation, the formation of linalool complexes in solid state was also studied. As described previously, linalool complexes could be made by using different energy sources to reach equilibrium: microwave irradiation (MWI) or ultrasound (US). Figure 2 shows the solubility profiles of linalool in the presence of HP-β-CDs a pH 7.0, by reaching the equilibrium using both methods. As can be seen in Figure 2A, the solutions treatment with MWI (24 h MWI or 48 h MWI) to reach the equilibrium, much less increased the linalool solubility than the application of US ( Figure 2B). When the increase in solubility of linalool was compared through the preparation of complexes by using MWI or US, as an energy source, at a volume of 100 mL, it was found to be 18.2 and 13.1fold, respectively, thus achieving a higher increase in solubility of linalool by MWI. This result is due to the fact that the S0 in the case of MWI (0.75 mmol L −1 ) is lower than that obtained by ultrasound (7.37 mmol L −1 ). In the case of MWI the phenomenon of molecular heating may not be enough to When the increase in solubility of linalool was compared through the preparation of complexes by using MWI or US, as an energy source, at a volume of 100 mL, it was found to be 18.2 and 13.1-fold, respectively, thus achieving a higher increase in solubility of linalool by MWI. This result is due to the fact that the S 0 in the case of MWI (0.75 mmol L −1 ) is lower than that obtained by ultrasound (7.37 mmol L −1 ). In the case of MWI the phenomenon of molecular heating may not be enough to cause the linalool oil micro drops disaggregation. However, in the case of US, the molecular movement favors the drops disaggregation.
Therefore, the contact surface is reduced and consequently its aqueous solubility, as shown in Figure 2A. In any case, in order to obtain an aqueous solution with high concentration of linalool, the use of US should be selected as appropriate method. As can be seen in Figure 2, by using 100 mM Molecules 2020, 25, 5069 6 of 18 of HP-β-CDs, the linalool concentration solved by using MWI irradiation was 13.61 mM, whereas in the case of US it was 96.72 mM (7-fold higher).
Soluble complexes prepared with different concentrations of HP-β-CDs were subjected to an atomization process (spray drying) to obtain solid complexes [34,35]. The shape of the particles obtained after the atomization process can be influenced by different factors such as the solid content in the solution, the inlet temperature, the spraying speed and the air nozzle of the equipment, which affect to the evaporation rate [36]. The wrinkled shape of the particles obtained for the linalool-HP-β-CDs solid complexes (Figure 3), was due to the fact that the evaporation rate on the surface of the droplets was faster than the rate at which the dissolved components were diffuse towards the center, to form a crust around the droplets [37].
Molecules 2020, 24, x FOR PEER REVIEW 6 of 19 cause the linalool oil micro drops disaggregation. However, in the case of US, the molecular movement favors the drops disaggregation. Therefore, the contact surface is reduced and consequently its aqueous solubility, as shown in Figure 2A. In any case, in order to obtain an aqueous solution with high concentration of linalool, the use of US should be selected as appropriate method. As can be seen in Figure 2, by using 100 mM of HP-β-CDs, the linalool concentration solved by using MWI irradiation was 13.61 mM, whereas in the case of US it was 96.72 mM (7-fold higher).
Soluble complexes prepared with different concentrations of HP-β-CDs were subjected to an atomization process (spray drying) to obtain solid complexes [34,35]. The shape of the particles obtained after the atomization process can be influenced by different factors such as the solid content in the solution, the inlet temperature, the spraying speed and the air nozzle of the equipment, which affect to the evaporation rate [36]. The wrinkled shape of the particles obtained for the linalool-HPβ-CDs solid complexes (Figure 3), was due to the fact that the evaporation rate on the surface of the droplets was faster than the rate at which the dissolved components were diffuse towards the center, to form a crust around the droplets [37]. As shown in Figures 3B and C, the external surfaces showed continuous walls without cracks, which significantly influence the loss of volatile compounds. This fact has also been described for oregano essential oil [38], bay leaf infusions [39] and thymol [9].
Subsequently, the performance of the drying process was calculated using equations (4) and (5), preparing the complexes using different energy sources to achieve balance: MWI or US, without taking into account the dust that remains inside the equipment, since it is considered a negligible amount. As can be seen in Figure 4A, the yield is higher than 56% in all cases, so it can be considered such a good performance for both methods. However, it is evident that by increasing the concentration of HP-β-CDs a higher yield is obtained for both MWI and US, reaching the maximum yield at 100 mM for the US (80%). If compare the yield for the complexes obtained by using MWI and US at 100 mM HP-β-CDs, values of 71.9 and 80.5% were obtained, respectively ( Figure 4A). Similar results were obtained for thymol complexes obtained by using the same energy sources as linalool [9]. Furthermore, it must be taken into account that the volume introduced into the equipment were small (a total volume of 100 mL), therefore losses can always occur, since on a laboratory scale it is difficult to achieve yields greater than 80% [40]. As shown in Figure 3B,C, the external surfaces showed continuous walls without cracks, which significantly influence the loss of volatile compounds. This fact has also been described for oregano essential oil [38], bay leaf infusions [39] and thymol [9].
Subsequently, the performance of the drying process was calculated using Equations (4) and (5), preparing the complexes using different energy sources to achieve balance: MWI or US, without taking into account the dust that remains inside the equipment, since it is considered a negligible amount. As can be seen in Figure 4A, the yield is higher than 56% in all cases, so it can be considered such a good performance for both methods. However, it is evident that by increasing the concentration of HP-β-CDs a higher yield is obtained for both MWI and US, reaching the maximum yield at 100 mM for the US (80%). If compare the yield for the complexes obtained by using MWI and US at 100 mM HP-β-CDs, values of 71.9 and 80.5% were obtained, respectively ( Figure 4A). Similar results were obtained for thymol complexes obtained by using the same energy sources as linalool [9]. Furthermore, it must be taken into account that the volume introduced into the equipment were small (a total volume of 100 mL), therefore losses can always occur, since on a laboratory scale it is difficult to achieve yields greater than 80% [40]. Therefore, the determination of the losses of linalool due to the drying treatment was analyzed and the concentration of linalool in the solid complexes after spraying was evaluated ( Figure 4B.). These results were used to calculate the encapsulation efficiency of the complexes obtained by using both energy sources, based on the theoretical amount of linalool in the initial solution and its relationship with the amount of linalool present in the solid complexes at the end of the process. The encapsulation efficiency ranged from 42.6 to 60.8 g/100 g for the MWI and from 61.3 to 93.7 g/100 g for the US ( Figure 4B).
These results can be justified by the fact that CDs form both inclusion and non-inclusion complexes, and the complexes form soluble aggregates capable of also solubilizing linalool molecules forming structures similar to micelles. This justifies the high concentration of linalool in the hydrophobic cavity when the complexes were obtained by using US, while in the case of MWI more molecules remained outside the cavity forming complexes outside of the hydrophobic cavity, favoring the linalool losses in the atomization process, thus decreasing the efficiency [41]. Therefore, although by means of MWI a greater increase in the solubility of linalool was obtained, it was observed that the complexes obtained by US contained higher concentration of linalool, since a better performance and encapsulation efficiency is obtained.

Stability of Solid Complexes
The stability of the solid complexes was also studied, evaluating the linalool content in the powder with the storage time at different temperatures. As shown in Figure 5 the powder obtained by the US maintained a greater amount of linalool retained in the samples stored at 8 °C. This was probably because at lower temperatures, there is also a lower humidity relative (RH) and therefore, fewer water molecules to compete with linalool for the hydrophobic cavity of the CDs. In fact, it has been previously described that HP-β-CDs have a higher affinity for water molecules as RH increases due to the amount of hydroxyl groups in the CDs surface, which enables the release of the complexed compound [42]. Therefore, the determination of the losses of linalool due to the drying treatment was analyzed and the concentration of linalool in the solid complexes after spraying was evaluated ( Figure 4B).
These results were used to calculate the encapsulation efficiency of the complexes obtained by using both energy sources, based on the theoretical amount of linalool in the initial solution and its relationship with the amount of linalool present in the solid complexes at the end of the process. The encapsulation efficiency ranged from 42.6 to 60.8 g/100 g for the MWI and from 61.3 to 93.7 g/100 g for the US ( Figure 4B).
These results can be justified by the fact that CDs form both inclusion and non-inclusion complexes, and the complexes form soluble aggregates capable of also solubilizing linalool molecules forming structures similar to micelles. This justifies the high concentration of linalool in the hydrophobic cavity when the complexes were obtained by using US, while in the case of MWI more molecules remained outside the cavity forming complexes outside of the hydrophobic cavity, favoring the linalool losses in the atomization process, thus decreasing the efficiency [41]. Therefore, although by means of MWI a greater increase in the solubility of linalool was obtained, it was observed that the complexes obtained by US contained higher concentration of linalool, since a better performance and encapsulation efficiency is obtained.

Stability of Solid Complexes
The stability of the solid complexes was also studied, evaluating the linalool content in the powder with the storage time at different temperatures. As shown in Figure 5 the powder obtained by the US maintained a greater amount of linalool retained in the samples stored at 8 • C. This was probably because at lower temperatures, there is also a lower humidity relative (RH) and therefore, fewer water molecules to compete with linalool for the hydrophobic cavity of the CDs. In fact, it has been previously described that HP-β-CDs have a higher affinity for water molecules as RH increases due to the amount of hydroxyl groups in the CDs surface, which enables the release of the complexed compound [42].

1 H and 2D NMR Spectroscopy
To confirm the formation of the linalool-HP-β-CDs inclusion complexes, the interaction of both molecules was studied using the NMR technique, since it offers evidence of the inclusion of the host molecule inside the CDs [43]. Thus, the results obtained after subjecting the complexes to 1 H-NMR confirmed that the stoichiometry was 1:1. Table 2 shows the chemical shift values of linalool and HP-β-CDs, both in their free state and complexed in methanol-d4 solution, as well as the differences between the signals of the free and complexed molecules. The change that occurs in the chemical shifts of the free molecule with respect to the complex, verify the formation of the inclusion complex.
Once the chemical shifts were obtained, a study was performed using two-dimensional (2D) NMR spectroscopy, since it provides important information on the spatial proximity between the host (HP-β-CDs) and linalool, through observation of cross correlations. This technique presents the so-called Overhauser effect (NOE) that originates from the interaction between two protons close in space, being able to observe this interaction in the nuclear Overhauser effect spectroscopy (NOESY) or rotating-frame nuclear Overhauser effect spectroscopy (ROESY) spectra, considering that the presence of crossed peaks (Figure 6), between the protons of two species in the NOE spectrum, confirms the existence of a spatial contact at 0.4 nm [9].

1 H and 2D NMR Spectroscopy
To confirm the formation of the linalool-HP-β-CDs inclusion complexes, the interaction of both molecules was studied using the NMR technique, since it offers evidence of the inclusion of the host molecule inside the CDs [43]. Thus, the results obtained after subjecting the complexes to 1 H-NMR confirmed that the stoichiometry was 1:1. Table 2 shows the chemical shift values of linalool and HP-β-CDs, both in their free state and complexed in methanol-d 4 solution, as well as the differences between the signals of the free and complexed molecules. The change that occurs in the chemical shifts of the free molecule with respect to the complex, verify the formation of the inclusion complex. Once the chemical shifts were obtained, a study was performed using two-dimensional (2D) NMR spectroscopy, since it provides important information on the spatial proximity between the host (HP-β-CDs) and linalool, through observation of cross correlations. This technique presents the so-called Overhauser effect (NOE) that originates from the interaction between two protons close in space, being able to observe this interaction in the nuclear Overhauser effect spectroscopy (NOESY) or rotating-frame nuclear Overhauser effect spectroscopy (ROESY) spectra, considering that the presence of crossed peaks (Figure 6), between the protons of two species in the NOE spectrum, confirms the existence of a spatial contact at 0.4 nm [9].  In order to obtain more conformational information, the 2D ROESY of the inclusion complex was obtained (Figure 6), showing an appreciable correlation (spatial proximity) of the -OH proton of linalool with the H-5 protons of HP-β-CDs. Similar interactions appeared between L4, L8 and L11 protons of linalool with the H1, H2, H3 and H4 protons of HP-β-CDs. These results confirm that linalool was included in the pocket of HP-β-CDs. In order to obtain more conformational information, the 2D ROESY of the inclusion complex was obtained (Figure 6), showing an appreciable correlation (spatial proximity) of the -OH proton of linalool with the H-5 protons of HP-β-CDs. Similar interactions appeared between L4, L8 and L11 protons of linalool with the H1, H2, H3 and H4 protons of HP-β-CDs. These results confirm that linalool was included in the pocket of HP-β-CDs.

Molecular Docking
In order to understand how HP-β-CDs interact with linalool once complexed, molecular coupling simulations (docking) were performed (Figure 7). This technique showed the possible three-dimensional interaction between the HP-β-CDs atoms with linalool. It was confirmed that the hydrogen of the hydroxyl group of linalool (oxygen in red and hydrogen in gray), interacts with the H5 hydrogen atoms of HP-β-CDs. In fact, the simulation showed how the oxygen from the hydroxyl group of the CDs (red arm), interacts with linalool by forming hydrogen bonds (dashed yellow line), with the H of the -OH group (Figure 7), interacting the L4, L8 and L11 hydrogens of linalool with the hydrogen atoms H1, H2, H3 and H4 of HP-β-CDs. The simulation results agree with the data obtained by 1 H-NMR and 2D-ROESY spectroscopy ( Figure 6). Furthermore, as shown in the spherical conformation of the molecules (Figure 7), specific molecular sites in linalool are strongly attached to the hydrophobic sites in HP-β-CDs.
coupling simulations (docking) were performed (Figure 7). This technique showed the possible threedimensional interaction between the HP-β-CDs atoms with linalool. It was confirmed that the hydrogen of the hydroxyl group of linalool (oxygen in red and hydrogen in gray), interacts with the H5 hydrogen atoms of HP-β-CDs. In fact, the simulation showed how the oxygen from the hydroxyl group of the CDs (red arm), interacts with linalool by forming hydrogen bonds (dashed yellow line), with the H of the -OH group (Figure 7), interacting the L4, L8 and L11 hydrogens of linalool with the hydrogen atoms H1, H2, H3 and H4 of HP-β-CDs. The simulation results agree with the data obtained by 1 H-NMR and 2D-ROESY spectroscopy ( Figure 6). Furthermore, as shown in the spherical conformation of the molecules (Figure 7), specific molecular sites in linalool are strongly attached to the hydrophobic sites in HP-β-CDs.

Differential Scanning Calorimetry (DSC) and Thermogravimetric Analysis (TG)
DSC and TG are commonly used to determine the physical properties of a substance. Through this study, the DSC thermograms of the complexes and the isolated reactants were obtained and shown in Figures 8 and 9. Linalool DSC curve presented two endothermic band, one close to 100 °C and other near 210 °C, which corresponds to the oxidation and transition from a liquid to a vapour phase, respectively (Figure 8a).

Differential Scanning Calorimetry (DSC) and Thermogravimetric Analysis (TG)
DSC and TG are commonly used to determine the physical properties of a substance. Through this study, the DSC thermograms of the complexes and the isolated reactants were obtained and shown in Figures 8 and 9. Linalool DSC curve presented two endothermic band, one close to 100 • C and other near 210 • C, which corresponds to the oxidation and transition from a liquid to a vapour phase, respectively (Figure 8a). For HP-β-CDs, due to its amorphous nature, a broad endothermic peak was observed at approximately 70 °C (Figure 8c) associated with the loss of water molecules. Furthermore, a small variation appears at 210 °C due to the transformation of the molecule [44,45], and a small decrease in the curve around 300 °C, probably due to the breakdown of HP-β-CDs. However, when linalool was complexed inside HP-β-CDs, an evident reduction of these signals occurred, suggesting a process of

Fourier Transform Infrared Spectroscopy (FTIR)
FTIR is a technique commonly used to confirm the formation of inclusion complexes [49]. For this reason, it was used to demonstrate the complexes formation between linalool and HP-β-CDs. Complexes and isolated molecules infrared spectra were shown in Figure 10. For HP-β-CDs, due to its amorphous nature, a broad endothermic peak was observed at approximately 70 • C (Figure 8c) associated with the loss of water molecules. Furthermore, a small variation appears at 210 • C due to the transformation of the molecule [44,45], and a small decrease in the curve around 300 • C, probably due to the breakdown of HP-β-CDs. However, when linalool was complexed inside HP-β-CDs, an evident reduction of these signals occurred, suggesting a process of water exclusion during complex formation (Figure 8b).
The DSC curve of the inclusion complex formed by linalool-HP-β-CDs (Figure 8b) did not show the characteristic endothermic peaks of linalool (Figure 8a), indicating that the monoterpene was protected of heat treatment due to the formation of the inclusion complex with HP-β-CDs [46]. Similar results to those were described previously for linalool with β-CDs [47].
Once the results obtained by DSC were justified, the derivatives of the thermogravimetric analysis were interpreted. As can be seen in Figure 9a, the TG curve of HP-β-CDs showed a weight loss of 5%. However, the curve of the complex formed by linalool and HP-β-CDs showed a weight loss of 2.58%, because in the complex there was lower water molecules inside of the CDs cavity, since they have been displaced by linalool. In addition, a full CDs (linalool and HP-β-CDs) would be denser than an empty CDs (HP-β-CDs). Similar results were described by Fernandes et al. (2004), for linalool complexes with β-CDs [48].

Fourier Transform Infrared Spectroscopy (FTIR)
FTIR is a technique commonly used to confirm the formation of inclusion complexes [49]. For this reason, it was used to demonstrate the complexes formation between linalool and HP-β-CDs. Complexes and isolated molecules infrared spectra were shown in Figure 10.

Fourier Transform Infrared Spectroscopy (FTIR)
FTIR is a technique commonly used to confirm the formation of inclusion complexes [49]. For this reason, it was used to demonstrate the complexes formation between linalool and HP-β-CDs. Complexes and isolated molecules infrared spectra were shown in Figure 10. Respect to the linalool IR spectrum (Figure 10), stretching bands associated with the -OH bonds appeared at λ = 3396 cm −1 , out-of-plane flexion bands of -OH at λ = 1457 cm −1 , and two bands of symmetric and asymmetric stretching of the methyl to λ = 2968 and λ = 2923 cm −1 , respectively. A tension band of the -CO group appears at λ = 917 cm −1 . On the other hand, the tension bands of the allyl group (C=C) are reflected in a peak that appears at λ = 1644 cm −1 . With respect to methyl (-CH 3 ) substituents appear at λ = 1375 cm −1 . As shown in Figure 10, the spectra of HP-β-CDs and the complexes formed by linalool and HP-β-CDs were practically the same, except for some bands. The first of these is the -OH group tension band that was at 3342 cm −1 and moved to 3351 cm −1 , this change occurred due to the intermolecular hydrogen bridge formed between linalool and HP-β-CDs. Another band that was modified is the one that corresponds to the -CO tension vibration zone, which appeared at 1147 cm −1 and moved at 1149 cm −1 in the complex, indicating that various interactions occur between the CDs and linalool. The characteristic band of the hydroxypropyl group of CDs did not undergo any change due to complexation, indicating that this group does not interact with linalool, confirming that the molecule is inside HP-β-CDs.

Solubility Studies
The study of linalool solubility in the presence of CDs was carried out according to the method described by Higuchi and Connors in 1965 with slight modifications [50]. CDs were dissolved in buffered solutions at different pHs: 3.5, 5.5 (100 mmol L −1 sodium acetate buffer), 6.5, 7.0 (100 mmol L −1 sodium phosphate buffer) and 8.5 (100 mmol L −1 sodium borate buffer). The αand HP-β-CDs were prepared up to 100 mmol L −1 . However, β-CDs were prepared up to 15 mmol L −1 , since this is their limit of aqueous solubility.
A saturating concentration of linalool was added to each CDs solution and they were kept in an ultrasound bath for 60 min in the dark, at 25 • C, until equilibrium was reached. Subsequently, the solutions were filtered through a 0.45 µm nylon membrane filter (Chromafil, Macherey-Nagel, Düren, Germany) to remove excess linalool.
Once filtered, the solutions were diluted with ethanol until reach ethanol 80% (solution:ethanol, 20:80, v:v) to quantify by gas chromatography-mass spectrometry (GC-MS) the linalool concentration in each solution. Phase solubility diagrams were made in triplicate.

Complexation Process by Using Microwave Irradiation (MWI)
The encapsulation of linalool with HP-β-CDs by using MWI was carried out following the method described by Hernández-Sánchez et al. in 2017 [32]. Aqueous solutions of HP-β-CDs (100 mL, 0-100 mmol L −1 ) were irradiated with microwave (LG Grill Wavedom, LG Electronics Las Rozas, Spain), at 700 W for 30 s, at 10 s intervals, until reaching a temperature of 70 • C in the solution. Later on, linalool was added to each aqueous solution, in a molar ratio of 1:1 and then, were irradiated again for 30 s at 10 s intervals, until reaching again 70 • C. Subsequently, the samples were shaken and kept for 12 h in dark sealed vials at 25 • C, before dividing them into two groups. After that, the first group, named 24 h MWI, was filtered using 0.45 µm nylon membrane filter (Chromafil, Macherey-Nagel). The second group, named 48 h MWI, were irradiated again, until reaching 70 • C, and maintain 12 h in the dark before they were filtered through a 0.45 µm nylon filter to remove linalool excess. Once filtered, both groups of samples were diluted with ethanol (20:80, v:v) until the GC-MS quantification of the linalool. All measurements were made in triplicate.

Quantification of Linalool by GC-MS Analysis
The quantification of linalool was carried out by using a Shimadzu GC-QP 2010 (Kyoto, Japan) gas chromatographer combined with a mass spectrometer. Helium was used as carrier gas, at a flow rate of 0.5 mL min −1 . A ω-WAX 250 fused silica supelco column (30 m × 0.25 mm × 0.25 µm thickness), was used. The temperature conditions were as follows: initial temperature at 70 • C, raised to 160 • C at 4 • C min −1 , raised to 280 • C at 30 • C min −1 , and maintained finally at 280 • C for 6 min. Injector temperature was 250 • C and injector mode was Split 1:20.
The peak area of each sample was used for linalool quantification (mmol L −1 ), by interpolating in the calibration curve obtained using a standard of linalool, defined by equation: Area = 2,223,720 − 2.50 × 10 4 [linalool (mmol L −1 )] and (R 2 = 0.9989) for linalool concentration from 0.0 to 0.5 mmol L −1 .

Complexation Constant (K C ) and Complexation Efficiency (CE) Calculation
Kc between linalool and CDs was calculated from the slope of the phase solubility profile and the solubility of linalool in aqueous solution (S 0 ) by using Equation (1) [50]: CE is the ratio between dissolved complex and free CDs concentration. It is independent of S 0 , and was calculated by using Equation (2) [51]: The molar ratio (MR) linalool:CDs, was calculated using CE values with Equation (3) [51].

Preparation of Linalool-HP-β-Cyclodextrin Inclusion Complexes by Spray Dryer
The equipment used was Buchi B-290 (Flawil, Switzerland). The variables used by the device were the following: inlet and outlet air temperature of 170 and 68 • C, 35 m 3 h −1 for inlet air flow, 5 mL min −1 for pump flow and 360 L h −1 for compressed air flow. The recovered powder was stored in an airtight glass container at 25 and 4 • C.
The drying process yield was calculated using Equation (4): Drying process yield = dehydrated complexes obtained (g) total solids in solution (kg) The encapsulation yield was calculated using Equation (5): Linanool yield = total linalool in dehydrated complexes (g) total linaool in disolved complexes (kg) To quantify the linalool present in the solid complexes, they were solved in water (complex:water, 1:1, w:v). Solutions were then filtered using 0.45 µm nylon membrane filter (Chromafil, Macherey-Nagel). Subsequently, the dissolved filtered complexes were diluted with ethanol (20:80, v:v). The quantification of the linalool content in each solution was carried out by GC-MS in triplicate. Finally, the stability of the dehydrated HP-β-CDs-linalool complexes at 4 • C and 25 • C was measured for 17 months. All samples were analyzed in triplicate.

Field Emission Scanning Electron Microscope (FESEM) Images
The solid complexes were examined by field emission scanning electron microscopy (FESEM) using a MERLIN™ VP COMPACT system (Carl Zeiss Microscopy SL, Kelsterbach, Germany). The microscopy images were taken using a SE2 detector under an accelerating voltage of 1 kV.

1 H and 2D NMR Spectroscopy
1 H-NMR spectra of linalool, CDs, and the inclusion complexes (dissolved in D 2 O) were recorded on an Avance 600 MHz spectrometer (Bruker, Karlsruhe, Germany) at 25 • C. Chemical shifts given in parts per million (ppm), are relative to a tetramethyl silane internal standard (δ = 0.0), and NMR data were processed with MestReNova software (6.0.2-5475 version). Two-dimensional rotational frame nuclear Overhauser effect spectroscopy (2D ROESY) spectra using the standard Bruker pulse program roesygpph were acquired at 32 scans, an acquisition time of 0.150 s and a pulse delay of 2.3 s.

Molecular Docking
In this study the molecular structures for linalool and CDs were manually constructed using AutoDockTools together with the structural information derived from experimental data [52]. The β-CDs structure was extracted from the protein structure of the Protein Data Bank (PDB) with the code 3CGT. The structure of the HP-β-CDs model was constructed by adding hydroxypropyl groups to the β-CDs model. Molecular coupling calculations were carried out using predetermined parameters [53]. The hydroxypropyl groups of HP-β-CDs were considered explicitly flexible during coupling simulations. Graphical representations of the coupling results were prepared using PyMOL (Molecular Graphics System, version 1.3, Schrödinger, LLC).

Differential Scanning Calorimetry (DSC) and Thermogravimetric Analysis (TG)
For thermal analyses, 4-5 mg of linalool, HP-β-CDs, or complexes were weighed, to the nearest 0.1 mg in aluminum capsules. They were analyzed by differential scanning calorimetry (DSC) using a Mettler DSC Q100 analyzer (TA Instruments, Cerdanyola del Valles, Spain), and by thermogravimetry with a Hi Res TAG 2950 thermogravimetric analysis equipment (TA Instruments). The scanning temperature range was from 25 • C to 300 • C, with a scanning temperature variation of 10 • C min −1 , using nitrogen as the carrier gas. The thermal stability of the respective components was shown using first derivative graphs (DTG) of weight (%) versus temperature ( • C).

Fourier Transform Infrared Spectroscopy (FTIR)
To determine the structural changes produced by the linalool complexation process with HP-β-CDs, the FTIR spectra were analyzed. The device used was a Varian FTIR 670 spectrophotometer (Agilent Tech, Amstelveen, The Netherlands) in combination with an accessory to analyze the attenuated total reflectance (ATR) with a wave number resolution of 0.10 cm −1 in the range of 250-4000 cm −1 . A minimum of 32 scans with a resolution of 4 cm −1 were averaged over the above ranges.

Conclusions
This study demonstrated that linalool could be complexed with different types of native and modified CDs, increasing its aqueous by means the formation of 1:1 complexes.
The inclusion complexes of linalool and CDs depend on the pH of the reaction medium, obtaining better Kc values of at a neutral pH. The highest Kc values were obtained with HP-β-CDs, followed by β-CDs, and α-CDs. Attending the results obtained of Kc values, encapsulation efficiency and stability of solid complexes, the US method is the most adequate for preparing linalool solid complexes. The results obtained by NMR and 2D-ROESY, in conjunction with molecular docking report the most probable conformation of linalool after complexation with HP-β-CDs. In addition, the analyses carried out by TG, DSC, and FTIR on the free and complexed molecule, confirm that the complex formation between linalool and HP-β-CDs was effective.