Structural Investigation of Beta-Cyclodextrin Complexes with Cannabidiol and Delta-9-Tetrahydrocannabinol in 1:1 and 2:1 Host-Guest Stoichiometry: Molecular Docking and Density Functional Calculations

The complexation of β-cyclodextrin (β-CD) with cannabidiol (CBD) and Δ9-tetrahydrocannabinol (THC) was investigated using molecular docking and M062X/6-31G(d,p) calculations. The calculations suggested two possible complex formations of 1:1 and 2:1 host-guest molecular ratio of β-CD with CBD and THC. The preferred orientation of all complexes in this study exhibited the hydrogen bonding between hydroxy-substituted benzene ring of CBD and THC with the β-CD’s secondary hydroxy groups at the wide rim. The calculated complexation energies indicate that formation of the 2:1 complexes (−83.53 to −135.36 kcal/mol) was more energetically favorable and chemically stable than the 1:1 complexes (−30.00 to −34.92 kcal/mol). However, the deformation energies of the host and the guest components in the 2:1 complexes (37.47–96.91 kcal/mol) are much higher than those in the 1:1 complexes (3.49–8.69 kcal/mol), which means the formation processes of the 2:1 complexes are more difficult due to the rigidity of the dimeric β-CDs. Therefore, the inclusion complexes of β-CD with CBD and THC are more likely to be in 1:1 host-guest ratio than in 2:1 molecular ratio. The results of this study supported the experimental results that the complexation constant of 1:1 β-CD/CBD (Ks = 300 M−1) is greater than that of 2:1 β-CDs/CBD (Kss = 0.833 M−1). Altogether, this study introduced the fitting parameters that could indicate the stability of the molecular fits in complex formation of each stoichiometry host-guest ratio, which are important for the assessment of the inclusion mechanisms as well as the relationships of reactants and products in chemical reactions.


Introduction
Cannabidiol, often known as CBD, is a phytocannabinoid that is derived from the Cannabis plant. It does not produce any euphoric effects, but it does have pain-relieving [1], anti-inflammatory [2,3], anti-cancer [4], and cancer-preventing properties [5][6][7]. When CBD is taken orally, it immediately provides anti-proliferative, anti-angiogenic, and proapoptotic effects through a variety of mechanisms. Additionally, there is some evidence that CBD can encourage the uptake of cytotoxic chemicals by cancer cells. CBD is a cannabinoid that can be used as a treatment for children who suffer from refractory epilepsy as a result of Lennox-Gastaut syndrome or Dravet syndrome [8,9]. ∆ 9 -tetrahydrocannabinol, also known as THC, is the principal component of the marijuana plant that is responsible for its euphoric and psychoactive effects. In 1964, Raphael Mechoulam, an Israeli chemist born in Bulgaria, made the initial discovery of THC and isolated it [10,11]. It was later discovered that when THC is smoked, it enters the circulation, travels to the brain, and attaches itself to naturally occurring cannabinoid receptors that are a part of the endocannabinoid system [12]. These receptors are found in the cerebral cortex, cerebellum, and basal ganglia [13,14]. These 2 of 12 are the areas of the brain that are in charge of thinking, remembering, feeling pleasure, coordinating movement, and other functions [15,16]. However, the low water solubility of CBD (0.0126 mg/mL) [17][18][19][20] and THC (0.0028 mg/mL) [21] becomes an obstacle to their applications in pharmaceutical fields.
In recent years, there have been several studies indicating that the water solubility and stability of these cannabinoids can be increased through the inclusion complex with β-cyclodextrins (β-CD), about 2.1 mg/mL [22][23][24][25][26]. These inclusion complexes also exhibit greater in vitro effects against human cancer lines Hep G2 and A549 [22,23]. β-CD is a macrocyclic ring consisting of seven glucose subunits joined by α-1, 4 glycosidic bonds. It has the shape of a truncated cone with a hydrophilic outer surface and lipophilic cavity [27]. All the secondary hydroxyl groups (corresponding to the C2 and C3 carbon atoms of the glucose units) are on the same edges of the cavity with wider rim (Head), whereas the primary hydroxyls are on the other end of the cavity with narrower rim (Tail). The cavity size is a major determinant as to which cyclodextrin (CD) should be used in complexation with hydrophobic molecules. The inner diameter of the wide rim of β-CD is 8.74 Å with a depth of 6.47 Å, as illustrated in Figure 1a. Several studies reported that CBD and THC can form complexes with β-CD in 1:1 and 2:1 host-guest stoichiometric ratios. The stoichiometry of the complex is given by the number of host and guest molecules contained in the inclusion complex. Lv et al. [22] reported that β-CD complexes with CBD in both of 1:1 and 2:1 host-guest stoichiometry. Li et al. [24] revealed that CBD formed inclusion complex with β-CD in 1:1 stoichiometric ratio. Mannila et.al. [25] achieved hostguest inclusion complex of randomly methylated β-CD with CBD and THC in 1:1 and 2:1 molecular ratio. Shoyama et al. [26] found a 2:1 stoichiometry in β-CD/THC system.

One-to-One Inclusion Complexes of β-CD with CBD and THC
The only difference between CBD and THC is a chemical bond (Figure 1b,c), which results in different structures, and chemical and pharmacological properties. Molecular docking is used to calculate the possibility of binding between β-CD either with CBD or THC molecules by fixing the host structure and allowing the guest to be flexible in the host's cavity. The calculations revealed two possible orientations in the 1:1 host-guest stoichiometry case of β-CD/CBD and only one orientation for β-CD/THC, as shown in Table 1 and Figure 2.  As the previous studies on the formation of β-CD complexes with CBD and THC implies the simultaneous presence of 1:1 and 2:1 complex, it is necessary to introduce fitting parameters to determine the reliability of the complex formation. Although many articles study the inclusion complexes of CBD and THC with CDs, as mentioned above, currently there are no articles investigating the molecular structures of β-CD/CBD and β-CD/THC in 1:1 and 2:1 host-guest stoichiometry, which are important for the assessment of the inclusion mechanisms as well as the relationships of reactants and products in chemical reactions. Therefore, in this work, molecular docking calculation was used to establish the most possible modes and calculate the binding energy of 1:1 and 2:1 complex formation of β-CD with CBD and THC. Furthermore, the energies of complexation, energy gap between HOMO (the highest occupied molecular orbital) and LUMO (the lowest unoccupied molecular orbital), and the deformation energies of the inclusion complex were assessed by density functional theory (DFT) calculations, to identify the fitting parameters that can indicate the stability of the molecular fits in complex formation of each stoichiometry host-guest ratio.

One-to-One Inclusion Complexes of β-CD with CBD and THC
The only difference between CBD and THC is a chemical bond (Figure 1b,c), which results in different structures, and chemical and pharmacological properties. Molecular docking is used to calculate the possibility of binding between β-CD either with CBD or THC molecules by fixing the host structure and allowing the guest to be flexible in the host's cavity. The calculations revealed two possible orientations in the 1:1 host-guest stoichiometry case of β-CD/CBD and only one orientation for β-CD/THC, as shown in Table 1 and Figure 2.   In term of configurations, in conformation I (Conf. I), the limonene group of CBD is located inside the cavity closer to the narrower side of β-CD (Tail), while its aliphatic chain is elongated outside the cavity closer to the wider rim of β-CD (Head), as displayed in Figure 2a. In conformation II (Conf. II), the benzenediol group of CBD stayed inside the cavity near the wider rim, while its limonene group stayed outside the cavity closer to the wider rim, and its aliphatic chain elongated outside the cavity closer to the narrow rim of β-CD, as presented in Figure 2b. Molecular docking calculations indicated that β-CD/CBD in Conf. I (99% frequency; the lowest ΔG = −6.87 kcal/mol) was more favorable than Conf. II (1% frequency; ΔG = −6.07 kcal/mol). The calculation on β-CD/THC shows that the preferred configuration is Conf. II (100% frequency) with the lowest ΔG = −6.77 kcal/mol, as reported in Table 1. The complex's orientation, which the aliphatic chain of guest molecule located closer to the narrow rim, would be referred to as Conf. II. Note that, in the β-CD/THC Conf. II, the aliphatic chain of THC was puckered inside the cavity, as depicted in Figure 2c, while in the β-CD/CBD Conf. II, the aliphatic chain of CBD is elongated outside the cavity of β-CD (Figure 2b).
Because the rigidity of β-CD molecule in molecular docking simulations was not realistic, the flexibility of the complex geometries and the intermolecular interactions of the In term of configurations, in conformation I (Conf. I), the limonene group of CBD is located inside the cavity closer to the narrower side of β-CD (Tail), while its aliphatic chain is elongated outside the cavity closer to the wider rim of β-CD (Head), as displayed in Figure 2a. In conformation II (Conf. II), the benzenediol group of CBD stayed inside the cavity near the wider rim, while its limonene group stayed outside the cavity closer to the wider rim, and its aliphatic chain elongated outside the cavity closer to the narrow rim of β-CD, as presented in Figure 2b. Molecular docking calculations indicated that β-CD/CBD in Conf. I (99% frequency; the lowest ∆G = −6.87 kcal/mol) was more favorable than Conf. II (1% frequency; ∆G = −6.07 kcal/mol). The calculation on β-CD/THC shows that the preferred configuration is Conf. II (100% frequency) with the lowest ∆G = −6.77 kcal/mol, as reported in Table 1. The complex's orientation, which the aliphatic chain of guest molecule located closer to the narrow rim, would be referred to as Conf. II. Note that, in the β-CD/THC Conf. II, the aliphatic chain of THC was puckered inside the cavity, as depicted in Figure 2c, while in the β-CD/CBD Conf. II, the aliphatic chain of CBD is elongated outside the cavity of β-CD (Figure 2b).
Because the rigidity of β-CD molecule in molecular docking simulations was not realistic, the flexibility of the complex geometries and the intermolecular interactions of the 1:1 β-CD/CBD and β-CD/THC systems were further examined by M062X/6-31g(d,p) density functional calculations. The energy-optimized structures of the β-CD/CBD and β-CD/THC inclusion complexes are depicted in Figure 3. The M062X/6-31G(d,p) energy-optimized conformations of the inclusion complex systems demonstrated that the arrangements of guest molecules were related to the starting docking geometries, as presented in Figure 3. This indicated a good performance and reliability of the molecular docking techniques for this investigation.
wider rim, and its aliphatic chain elongated outside the cavity closer to the narrow rim of β-CD, as presented in Figure 2b. Molecular docking calculations indicated that β-CD/CBD in Conf. I (99% frequency; the lowest ΔG = −6.87 kcal/mol) was more favorable than Conf. II (1% frequency; ΔG = −6.07 kcal/mol). The calculation on β-CD/THC shows that the preferred configuration is Conf. II (100% frequency) with the lowest ΔG = −6.77 kcal/mol, as reported in Table 1. The complex's orientation, which the aliphatic chain of guest molecule located closer to the narrow rim, would be referred to as Conf. II. Note that, in the β-CD/THC Conf. II, the aliphatic chain of THC was puckered inside the cavity, as depicted in Figure 2c, while in the β-CD/CBD Conf. II, the aliphatic chain of CBD is elongated outside the cavity of β-CD (Figure 2b).
Because the rigidity of β-CD molecule in molecular docking simulations was not realistic, the flexibility of the complex geometries and the intermolecular interactions of the 1:1 β-CD/CBD and β-CD/THC systems were further examined by M062X/6-31g(d,p) density functional calculations. The energy-optimized structures of the β-CD/CBD and β-CD/THC inclusion complexes are depicted in Figure 3. The M062X/6-31G(d,p) energyoptimized conformations of the inclusion complex systems demonstrated that the arrangements of guest molecules were related to the starting docking geometries, as presented in Figure 3. This indicated a good performance and reliability of the molecular docking techniques for this investigation.   Table 2 presents the distance of the intermolecular hydrogen bonds (H-bonds), which are found in M062X/6-31G(d,p) energy-optimized structures. In the studied complex systems, three types of H-bonds were identified. The first one, which is found in all systems, is between an ether-like anomeric oxygen atom in the n glucose unit of β-CD and a hydrogen atom of CBD's hydroxyl group (O4 (β-CD)n ··· H (OH-guest) ). The second one, which is found in β-CD/CBD Conf. I and β-CD/THC Conf. II, is from an oxygen atom of the secondary hydroxyl group at O2 of β-CD and the hydrogen atom of guest's hydroxy group (O2 (β-CD)n+1 ··· H (OH-guest) ). The third one, which is found only in β-CD/THC Conf. II, occurs between the oxygen atom in benzo[c]chromene group of THC and the hydrogen atom of the secondary hydroxyl at O3 of β-CD (O (THC) ··· H (O3H-β-CD)n-2 ). Table 2. Distance of hydrogen bonds between host (β-CD) and guest (CBD and THC), obtained from M062X/6-31g(d,p) energy-optimized structures.

Inclusion Complex
Distance (Å) The results obtained from M062X/6-31G(d,p) calculations for the 1:1 complexes are shown in Table 3, where µ sp (component) is the dipole moment of the single-point energy of the guest or the host taken from the optimized complex. Dipole moment results from non-uniform distributions of positive and negative charges on the atoms inside the system. The dipole moment values of β-CD/CBD and β-CD/THC systems are closer to the dipole moment of free β-CD component than those of free CBD or THC, which means that the guest molecules are likely encapsulated inside the hydrophobic cavity of β-CD. The large |HOMO-LUMO| gap indicated that all inclusion complexes are as chemically stable as the free host and guest molecules. Complexation energy (∆E) values are negative for all inclusion complexes which suggests that their formations are energetically favorable. Deformation energy (E DEF ) is the energy required to distort each monomer into the structure it adopts within the complex. The deformation energy of CBD molecule is higher than that of β-CD in both Conf. I and Conf. II, which indicates that the flexibility of the CBD structure plays an important role in stabilizing the whole system upon complexation. In the β-CD/THC complex, both components exhibit low value of deformation energies, which means that the complex formation is very favorable. This information is supported by small changes of µ sp (component) values compared to the dipole moments of isolated guest and host molecules. According to the values of ∆E and E DEF , the stability of 1:1 inclusion complexes is in the following order: β-CD/THC Conf. II > β-CD/CBD Conf. II > β-CD/CBD Conf. I.

Two-to-One Inclusion Complexes of β-CDs with CBD and THC
In this section, 2:1 host-guest stoichiometry of β-CDs with CBD and THC inclusion complexes, wherein two β-CD molecules are fetched closer together, is investigated with three dimerized form of β-CDs: facing secondary portals (Head-to-Head, HH), facing two primary portals (Tail-to-Tail, TT), and facing primary and secondary portals (Tail-to-Head, TH), as shown in Figure 4a-c, respectively.
Molecular docking calculation results in Table 4 indicate only one possible arrangement of CBD molecule inside HH, TT, and TH dimeric cavities, which hereafter would be referred to as Conformation III (Conf. III), Conformation IV (Conf. IV), and Conformation V (Conf. V), respectively. Figure 5 illustrates the orientations of CBD and THC inside the β-CDs dimer. The docking simulations of HH/THC and TT/THC complex systems also found only one preferred conformation, while in TH/THC system, three conformations were found, including Conf. V (80% frequency), Conformation VI (Conf. VI, 16% frequency), and Conformation VII (Conf. VII, 4% frequency), as presented in Figure 5. The binding energy (∆G) of the 1:1 complexes is in the range of −6.07 to −6.87 kcal/mol (Table 1), while ∆G of the 2:1 complexes is in the range of −7.57 to −9.61 kcal/mol (Table 4), which demonstrates that the 2:1 host-guest ratio complexation is more energy favorable than the 1:1 host-guest ratio.

Two-to-One Inclusion Complexes of β-CDs with CBD and THC
In this section, 2:1 host-guest stoichiometry of β-CDs with CBD and THC inclusion complexes, wherein two β-CD molecules are fetched closer together, is investigated with three dimerized form of β-CDs: facing secondary portals (Head-to-Head, HH), facing two primary portals (Tail-to-Tail, TT), and facing primary and secondary portals (Tail-to-Head, TH), as shown in Figure 4a-c, respectively.   Table 4 indicate only one possible arrangement of CBD molecule inside HH, TT, and TH dimeric cavities, which hereafter would be referred to as Conformation III (Conf. III), Conformation IV (Conf. IV), and Conformation V (Conf. V), respectively. Figure 5 illustrates the orientations of CBD and THC inside the β-CDs dimer. The docking simulations of HH/THC and TT/THC complex systems also found only one preferred conformation, while in TH/THC system, three conformations were found, including Conf. V (80% frequency), Conformation VI (Conf. VI, 16% frequency), and Conformation VII (Conf. VII, 4% frequency), as presented in Figure 5. The binding energy (ΔG) of the 1:1 complexes is in the range of −6.07 to −6.87 kcal/mol (Table  1), while ΔG of the 2:1 complexes is in the range of −7.57 to −9.61 kcal/mol (Table 4), which demonstrates that the 2:1 host-guest ratio complexation is more energy favorable than the 1:1 host-guest ratio.    Starting from the molecular docking geometry of each complex, the geometry was fully optimized without any constraint using M062X/6-31G(d,p) calculation. Figure 6 presents the energy-optimized conformations of 2:1 host-guest stoichiometry complex systems. The intermolecular H-bond distances between the host and the guest molecules, which are related to the complex stability, are reported in Table 5. Two types of H-bonds were found. The first type, which occurs in all complex systems, is between an ether-like anomeric oxygen atom in the n glucose unit of β-CD and a hydrogen atom of guest's hydroxyl group (O4(β-CD)n ··· H(OH-guest)). The second type, found only in HH/CBD Conf. III, is between the oxygen atom at O3 of β-CD and the hydrogen atom of the CBD's hydroxyl group (O3(β-CD)n ··· H(O1H-CBD)). Due to its size, the hydroxy-substituted benzene ring of CBD and THC is always located near the wider side of a β-CD molecule. In terms of docking conformations, only one molecular arrangement was found in each of HH/CBD, TT/CBD, and TH/CBD complexes. In all of them, the cyclohexene group of CBD was located in β-CD Top , while its aliphatic chain is located in β-CD Bottom , as shown in Figure 5a-c. Due to its length, the benzenediol group of CBD was always found by the side of the wider rim of a β-CD molecule, when CBD is incorporated into the host cavity. Only one conformation of HH/THC and TT/THC inclusion complexes were obtained from molecular docking simulations, with similar molecular arrangement as β-CDs/CBD. However, in TH/THC system, the docking results obtain three different conformations, as depicted in Figure 5c-e. The first is TH/THC Conf. V, where THC was entrapped inside the dimeric cavity with its cyclohexene group located in the β-CD Top and its aliphatic chain in β-CD Bottom (Figure 5c), as previously observed in TH/CBD systems. On the other hand, in TH/THC Conf. VI, the cyclohexene group of THC was immersed deep into the narrow side of β-CD Bottom , and the hydroxy-substituted benzene ring located at the wider rim of the same β-CD Bottom , while its aliphatic chain settled in β-CD Top , as shown in Figure 5d. The third conformation, TH/THC Conf. VII, is analogous to 1:1 β-CD/THC Conf. II, where the aliphatic chain of THC molecule puckered inside the cavity with the hydroxy-substituted benzene ring locating at the wider rim of the same β-CD, as shown in Figure 5e.

Molecular docking calculation results in
Starting from the molecular docking geometry of each complex, the geometry was fully optimized without any constraint using M062X/6-31G(d,p) calculation. Figure 6 presents the energy-optimized conformations of 2:1 host-guest stoichiometry complex systems. The intermolecular H-bond distances between the host and the guest molecules, which are related to the complex stability, are reported in Table 5. Two types of H-bonds were found. The first type, which occurs in all complex systems, is between an ether-like anomeric oxygen atom in the n glucose unit of β-CD and a hydrogen atom of guest's hydroxyl group (O4 (β-CD)n ··· H (OH-guest) ). The second type, found only in HH/CBD Conf. III, is between the oxygen atom at O3 of β-CD and the hydrogen atom of the CBD's hydroxyl group (O3 (β-CD)n ··· H (O1H-CBD) ). Due to its size, the hydroxy-substituted benzene ring of CBD and THC is always located near the wider side of a β-CD molecule.   Table 6 presents the energy results from M062X/6-31G(d,p) calculations, expressed as complexation, HOMO and LUMO, and deformation energies. The formation of 2:1 host-guest ratio complexes is energetically favorable and chemically stable, as indicated by the negative complexation energy values (ΔE) and the large |HOMO-LUMO| gaps. The deformation of β-CDs is much stronger than the deformation of guest molecule. This is likely due to the rigidity of the dimeric β-CDs. The highest deformation energies for TH configurations (−88.7794 to −96.9119 kcal/mol, negative sign referring to a change in energy) means that complex formation should be difficult. TH and HH configurations exhibit more deformability upon complexation as suggested by higher deformation energy compared to TT configuration. Although the inclusion complexes with CBD and THC in  Table 5. Distances of hydrogen bonds between the host (dimeric β-CDs) and the guest molecules (CBD and THC) obtained from the M062X/6-31g(d,p) energy-optimized structures.

Inclusion Complex
Distance (Å)  Table 6 presents the energy results from M062X/6-31G(d,p) calculations, expressed as complexation, HOMO and LUMO, and deformation energies. The formation of 2:1 host-guest ratio complexes is energetically favorable and chemically stable, as indicated by the negative complexation energy values (∆E) and the large |HOMO-LUMO| gaps.
The deformation of β-CDs is much stronger than the deformation of guest molecule. This is likely due to the rigidity of the dimeric β-CDs. The highest deformation energies for TH configurations (−88.7794 to −96.9119 kcal/mol, negative sign referring to a change in energy) means that complex formation should be difficult. TH and HH configurations exhibit more deformability upon complexation as suggested by higher deformation energy compared to TT configuration. Although the inclusion complexes with CBD and THC in TT configuration exhibit low deformation energy (−37.4730 and −38.3873 kcal/mol, respectively), they have unbalanced charge distribution as indicated by the high dipole moment (9.0612 debye and 8.4162 debye for TT/CBD and TT/THC, respectively), which ultimately lower the stability of compounds [28]. That is the reason why TT/guest systems have lower complexation energy than HH/guest and TH/guest systems. According to the values of ∆E, E DEF , and µ the orders of 2:1 inclusion complex stability is TH/guest Conf. V > HH/guest Conf. III > TT/guest Conf. IV. Finally, as suggested by molecular docking and density functional calculations, both of 1:1 and 2:1 host-guest stoichiometry of β-CD with CBD and THC are feasible. The calculated complexation energies indicate that the formation of the 2:1 host-guest ratio complex is more energetically favorable and chemically stable than the 1:1 host-guest ratio complex. However, the deformation energies of the host and the guest components in the 2:1 complex are much higher than those in the 1:1 complex, which means the formation process of 2:1 is more difficult. This explains why the experimental complexation constant of 1:1 β-CD/CBD (K s = 300 M −1 ) is greater than that of 2:1 β-CDs/CBD (K ss = 0.833 M −1 ) [22].

Molecular Structure Construction
The starting geometries of β-CD, CBD, and THC monomers were taken from the Cambridge Crystallographic Data Centre (CCDC) with deposition number 1,107,192 [29], 1,533,487 [30], and 702,456 [31], respectively. The pairwise initial arrangement of β-CDs, with facing secondary portals (Head-to-Head, HH) and facing primary portals (Tail-to-Tail, TT), were extracted from CCDC 166671 [32] and CCDC 864,041 [33], respectively. Modification of each x-ray structure for atoms elimination and atoms addition was performed with Discovery Studio 2020 Client (DSC) [34]. Furthermore, DSC was also used to duplicate and rotate β-CD molecule to construct the dimeric β-CDs with facing primary and secondary portals (Tail-to-Head, TH), of which the X-ray structure was not available.

Molecular Docking Calculation
The inclusion complex between β-CDs with CBD and THC molecules were constructed by molecular docking calculations using AutoDockTools 1.5.7 and AutoDock 4.2 [35] software packages with the Lamarckian genetic algorithm [36]. The first step was the preparation of the host (β-CD) and the guests (CBD and THC) coordinate files using AutoDockTools. Non-polar hydrogens were deleted, and their charges were merged with carbon atoms. Atom types, hydrogen bond donors and acceptors, aliphatic, and aromatic carbon atoms, as well as rotatable bonds of the guest molecules, were also defined, while the hosts were kept fixed. AutoGrid was used to calculate the grid maps, one for each atom type present in the guest being docked-typically carbon, oxygen, and hydrogen atoms. The systems were investigated in a volume divided into many small grid boxes with a grid spacing of 0.375 Å, where the grid center was set at the center of the host molecules. The investigated systems were in the size of 19.5Å × 19.5Å × 19.5Å and 24.75Å × 24.75Å × 24.75Å for 1:1 and 2:1 host-guest stoichiometry complex systems, respectively. Finally, a hundred docking calculations were performed on each host-guest complex. At the end of docking simulations, the findings were grouped to identify related conformations, or clusters, based on the all-atom root mean square deviation within 2 Å. The binding energy (∆G) and clustering information along with the coordinates for each docked conformation were recorded into the result files. The docked conformations with the lowest energy were selected for full geometry optimization, according to the previous studies [37][38][39][40].

Complexation Energy Calculation
Full geometry optimization was performed on the monomers as well as the 1:1 and 2:1 host-guest stoichiometry inclusion complex structures with the Minnesota density functional [41] M062X/6-31g(d,p) using the Gaussian 16 molecular modeling package [42]. The interaction between the host and the guest in the optimized geometries can be quantified by the complexation energy (∆E) in Equation (1): where E opt complex , E opt host , and E opt guest are the optimized energy of the inclusion complex, the isolated host molecule (β-CD), and the isolated guest molecule (CBD and THC), respectively. The lower complexation energy value, the more stable the inclusion complex. For the 2:1 host-guest stoichiometry inclusion complexes, the values of E host were double the optimized energy of isolated β-CD. According to the DFT calculations in this study, an energy-optimized dimeric β-CDs structures without the guest was not obtained, as the β-CD molecules were separated from each other during the optimization. This behavior is supported by the information that aggregation of cyclodextrins can be promoted by the presence of guest compounds [43].
The deformation energy (E DEF ) for each host and guest component throughout the formation of the inclusion complex was defined as in Equation (2): where E(component) opt sp is the single-point energy of the component taken from the optimized complex and E(component) opt is the energy of the optimized geometry of each free component.

Conclusions
In summary, we investigated the molecular structures of β-CD/CBD and β-CD/THC complexes in 1:1 and 2:1 host-guest stoichiometry. The calculated complexation energy and the HOMO and LUMO gaps of the inclusion complexes were used to indicate the stability of the molecular fits of the complexes' formation with each stoichiometry host-guest ratio, which demonstrate that the 2:1 host-guest ratio complexation is more energy favorable than the 1:1 host-guest ratio. According to the values of ∆E, E DEF , and µ, the stability of 2:1 inclusion complex is in the following order: TH/guest Conf. V > HH/guest Conf. III > TT/guest Conf. IV. However, the deformation energies of dimeric β-CDs in 2:1 ratio were significantly higher than those of β-CD molecule in 1:1 ratio, which means the formation process with dimeric β-CDs is more difficult. Therefore, the inclusion complex of β-CD with CBD and THC is more likely to be in a 1:1 host-guest ratio than in a 2:1 molecular ratio.

Data Availability Statement:
The data presented in the study are included in the article. Further inquiries can be directed to the corresponding author.

Conflicts of Interest:
The authors declare no conflict of interest.