Next Article in Journal
Anti-Aging Effects of Anthocyanin Extracts of Sambucus canadensis Caused by Targeting Mitochondrial-Induced Oxidative Stress
Next Article in Special Issue
Thermal Titration Molecular Dynamics (TTMD): Not Your Usual Post-Docking Refinement
Previous Article in Journal
Analysis of the Toll and Spaetzle Genes Involved in Toll Pathway-Dependent Antimicrobial Gene Induction in the Red Flour Beetle, Tribolium castaneum (Coleoptera; Tenebrionidae)
Previous Article in Special Issue
Ligand-Based Drug Design of Novel Antimicrobials against Staphylococcus aureus by Targeting Bacterial Transcription
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

School of Bio-Chemical Engineering and Technology, Sirindhorn International Institute of Technology, Thammasat University, Pathum Thani 12120, Thailand
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(2), 1525; https://doi.org/10.3390/ijms24021525
Submission received: 14 December 2022 / Revised: 8 January 2023 / Accepted: 9 January 2023 / Published: 12 January 2023
(This article belongs to the Special Issue New Avenues in Molecular Docking for Drug Design 2022)

Abstract

:
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.

1. 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 pro-apoptotic 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 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 host-guest 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.
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.

2. Results and Discussions

2.1. 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 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.
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).
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.

2.2. 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.
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 β-CDTop, while its aliphatic chain is located in β-CDBottom, 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 β-CDTop and its aliphatic chain in β-CDBottom (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 β-CDBottom, and the hydroxy-substituted benzene ring located at the wider rim of the same β-CDBottom, while its aliphatic chain settled in β-CDTop, 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.
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 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 (Ks = 300 M−1) is greater than that of 2:1 β-CDs/CBD (Kss = 0.833 M−1) [22].

3. Materials and Methods

3.1. 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.

3.2. 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].

3.3. 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):
Δ E = E complex opt ( E host opt + E guest opt )
where E complex opt , E host opt , and E guest opt 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):
E DEF ( component ) = E ( component ) sp opt     E ( component ) opt
where E ( component ) sp opt 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.

4. 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.

Author Contributions

Writing—original draft preparation: N.T. and L.L. Writing—review and editing: P.T. and L.L. Research fund support: L.L. All authors have read and agreed to the published version of the manuscript.

Funding

This study was supported by the Thammasat University Research Fund, Contract No. TUFT 42/2565.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

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.

References

  1. Mlost, J.; Bryk, M.; Starowicz, K. Cannabidiol for Pain Treatment: Focus on Pharmacology and Mechanism of Action. Int. J. Mol. Sci. 2020, 21, 8870. [Google Scholar] [CrossRef] [PubMed]
  2. Atalay, S.; Jarocka-Karpowicz, I.; Skrzydlewska, E. Antioxidative and Anti-Inflammatory Properties of Cannabidiol. Antioxidants 2020, 9, 21. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Henshaw, F.R.; Dewsbury, L.S.; Lim, C.K.; Steiner, G.Z. The Effects of Cannabinoids on Pro- and Anti-Inflammatory Cytokines: A Systematic Review of In Vivo Studies. Cannabis Cannabinoid Res. 2021, 6, 177–195. [Google Scholar] [CrossRef]
  4. Seltzer, E.S.; Watters, A.K.; MacKenzie, D., Jr.; Granat, L.M.; Zhang, D. Cannabidiol (CBD) as a Promising Anti-Cancer Drug. Cancers 2020, 12, 3203. [Google Scholar] [CrossRef]
  5. Elinav, E.; Nowarski, R.; Thaiss, C.A.; Hu, B.; Jin, C.; Flavell, R.A. Inflammation-Induced Cancer: Crosstalk between Tumours, Immune Cells and Microorganisms. Nat. Rev. Cancer 2013, 13, 759–771. [Google Scholar] [CrossRef]
  6. Porta, C.; Larghi, P.; Rimoldi, M.; Totaro, M.G.; Allavena, P.; Mantovani, A.; Sica, A. Cellular and Molecular Pathways Linking Inflammation and Cancer. Immunobiology 2009, 214, 761–777. [Google Scholar] [CrossRef]
  7. Mantovani, A.; Allavena, P.; Sica, A.; Balkwill, F. Cancer-related inflammation. Nature 2008, 454, 436–444. [Google Scholar] [CrossRef]
  8. Chen, J.W.; Borgelt, L.M.; Blackmer, A.B. Cannabidiol: A New Hope for Patients with Dravet or Lennox-Gastaut Syndromes. Ann. Pharmacother. 2019, 53, 603–611. [Google Scholar] [CrossRef] [PubMed]
  9. Koubeissi, M. Anticonvulsant Effects of Cannabidiol in Dravet Syndrome. Epilepsy Curr. 2017, 17, 281–282. [Google Scholar] [CrossRef]
  10. Pertwee, R.G. Cannabinoid pharmacology: The first 66 years. Br. J. Pharmacol. 2006, 147 (Suppl. S1), S163–S171. [Google Scholar] [CrossRef]
  11. Gaoni, Y.; Mechoulam, R. Isolation, Structure, and Partial Synthesis of an Active Constituent of Hashish. J. Am. Chem. Soc. 1964, 86, 1646–1647. [Google Scholar] [CrossRef]
  12. Almogi-Hazan, O.; Or, R. Cannabis, the Endocannabinoid System and Immunity-the Journey from the Bedside to the Bench and Back. Int. J. Mol. Sci. 2020, 21, 4448. [Google Scholar] [CrossRef]
  13. Freitas, H.R.; Ferreira, G.D.C.; Trevenzoli, I.H.; Oliveira, K.J.; de Melo Reis, R.A. Fatty Acids, Antioxidants and Physical Activity in Brain Aging. Nutrients 2017, 9, 1263. [Google Scholar] [CrossRef] [Green Version]
  14. Freitas, H.R.; Isaac, A.R.; Malcher-Lopes, R.; Diaz, B.L.; Trevenzoli, I.H.; De Melo Reis, R.A. Polyunsaturated Fatty Acids and Endocannabinoids in Health and Disease. Nutr. Neurosci. 2018, 21, 695–714. [Google Scholar] [CrossRef]
  15. Aizpurua-Olaizola, O.; Elezgarai, I.; Rico-Barrio, I.; Zarandona, I.; Etxebarria, N.; Usobiaga, A. Targeting the Endocannabinoid System: Future Therapeutic Strategies. Drug. Discov. Today 2017, 22, 105–110. [Google Scholar] [CrossRef]
  16. Donvito, G.; Nass, S.R.; Wilkerson, J.L.; Curry, Z.A.; Schurman, L.D.; Kinsey, S.G.; Lichtman, A.H. The Endogenous Cannabinoid System: A Budding Source of Targets for Treating Inflammatory and Neuropathic Pain. Neuropsychopharmacology 2018, 43, 52–79. [Google Scholar] [CrossRef] [Green Version]
  17. Chen, C.; Pan, Z. Cannabidiol and Terpenes from Hemp—Ingredients for Future Foods and Processing Technologies. J. Food Agric. Soc. 2021, 1, 113–127. [Google Scholar] [CrossRef]
  18. Koch, N.; Jennotte, O.; Gasparrini, Y.; Vandenbroucke, F.; Lechanteur, A.; Evrard, B. Cannabidiol Aqueous Solubility Enhancement: Comparison of Three Amorphous Formulations Strategies Using Different Type of Polymers. Int. J. Pharm. 2020, 589, 119812. [Google Scholar] [CrossRef]
  19. Chen, P.X.; Rogers, M.A. Opportunities and Challenges in Developing Orally Administered Cannabis Edibles. Curr. Opin. Food Sci. 2019, 28, 7–13. [Google Scholar] [CrossRef]
  20. Bragança, V.A.; França, T.G.; de Jesus, A.C.; Palheta, I.C.; Melo, F.P.; Neves, P.A.; Lima, A.B.; Borges, R.S. Impact of Conformational and Solubility Properties on Psycho-Activity of Cannabidiol (CBD) and Tetrahydrocannabinol (THC). Chem. Data Collect. 2020, 26, 100345. [Google Scholar] [CrossRef]
  21. Garrett, E.R.; Hunt, C.A. Physiochemical Properties, Solubility, and Protein Binding of Δ9-Tetrahydrocannabinol. J. Pharm. Sci. 1974, 63, 1056–1064. [Google Scholar] [CrossRef]
  22. Hatziagapiou, K.; Bethanis, K.; Koniari, E.; Christoforides, E.; Nikola, O.; Andreou, A.; Mantzou, A.; Chrousos, G.P.; Kanaka-Gantenbein, C.; Lambrou, G.I. Biophysical Studies and In Vitro Effects of Tumor Cell Lines of Cannabidiol and Its Cyclodextrin Inclusion Complexes. Pharmaceutics 2022, 14, 706. [Google Scholar] [CrossRef]
  23. Lv, P.; Zhang, D.; Guo, M.; Liu, J.; Chen, X.; Guo, R.; Xu, Y.; Zhang, Q.; Liu, Y.; Guo, H.; et al. Structural Analysis and Cytotoxicity of Host-Guest Inclusion Complexes of Cannabidiol with Three Native Cyclodextrins. J. Drug Deliv. Sci. Technol. 2019, 51, 337–344. [Google Scholar] [CrossRef]
  24. Li, H.; Zhao, Q.-S.; Chang, S.-L.; Chang, T.-R.; Tan, M.-H.; Zhao, B. Development of Cannabidiol Full-Spectrum Oil/2,6-Di-O-Methyl-β-Cyclodextrin Inclusion Complex with Enhanced Water Solubility, Bioactivity, and Thermal Stability. J. Mol. Liq. 2022, 347, 118318. [Google Scholar] [CrossRef]
  25. Mannila, J.; Järvinen, T.; Järvinen, K.; Tarvainen, M.; Jarho, P. Effects of RM-β-CD on Sublingual Bioavailability of Δ9-Tetrahydrocannabinol in Rabbits. Eur. J. Pharm. Sci. 2005, 26, 71–77. [Google Scholar] [CrossRef]
  26. Shoyama, Y.; Morimoto, S.; Nishioka, I. Cannabis. XV. Cannabis. XV. Preparation and Stability of Δ9-Tetrahydrocannabinol-β-Cyclodextrin Inclusion Complex. J. Nat. Prod. 1983, 46, 633–637. [Google Scholar] [CrossRef]
  27. Loftsson, T.; Jarho, P.; Másson, M.; Järvinen, T. Cyclodextrins in Drug Delivery. Expert Opin Drug Deliv 2005, 2, 335–351. [Google Scholar] [CrossRef]
  28. Petov, G.M.; Matrosov, E.I.; Kabachnik, M.I. Determination of dipole moments and stability constants for molecular complexex. Russ. Chem. Bull. 1978, 27, 2020–2024. [Google Scholar] [CrossRef]
  29. Steiner, T.; Koellner, G. Crystalline. β-Cyclodextrin Hydrate at Various Humidities: Fast, Continuous, and Reversible Dehydration Studied by X-ray Diffraction. J. Am. Chem. Soc. 1994, 116, 5122–5128. [Google Scholar] [CrossRef]
  30. Mayr, T.; Grassl, T.; Korber, N.; Christoffel, V.; Bodensteiner, M. Cannabidiol Revisited. IUCrData. 2017, 2, x170276. [Google Scholar] [CrossRef]
  31. Gul, W.; Carvalho, P.; Berberich, D.; Avery, M.A.; ElSohly, M.A. (6aR,10aR)-6,6,9-Trimethyl-3-pentyl6a,7,8,10a-tetrahydro-6H-benzo[c]chromen-1-yl 4-methylbenzenesulfonate. Acta Cryst. 2008, E64, o1686. [Google Scholar] [CrossRef] [Green Version]
  32. Makedonopoulou, S.; Mavridis, I.M. The Dimeric Complex of Cyclomaltoheptaose with 1,14-Tetradecanedioic Acid. Comparison with Related Complexe. Carbohydr. Res. 2002, 335, 213–220. [Google Scholar] [CrossRef]
  33. Grishina, A.; Stanchev, S.; Kumprecht, L.; Budesinsky, M.; Pojarova, M.; Dusek, M.; Rumlova, M.; Krizova, I.; Rulisek, L.; Kraus, T. β-Cyclodextrin Duplexes that are Connected through Two Disulfide Bonds: Potent Hosts for the Complexation of Organic Molecules. Chem. Eur. J. 2012, 18, 12292. [Google Scholar] [CrossRef]
  34. BIOVIA, Dassault Systèmes, Discovery Studio Visualizer, San Diego: Dassault Systèmes, 2020. Available online: https://discover.3ds.com/discovery-studio-visualizer-download (accessed on 12 May 2022).
  35. Morris, G.M.; Huey, R.; Lindstrom, W.; Sanner, M.F.; Belew, R.K.; Goodsell, D.S.; Olson, A.J. Autodock4 and AutoDockTools4: Automated Docking with Selective Receptor Flexiblity. J. Comput. Chem. 2009, 16, 2485–2791. [Google Scholar] [CrossRef] [Green Version]
  36. Morris, G.M.; Goodsell, D.S.; Halliday, R.S.; Huey, R.; Hart, W.E.; Belew, R.K.; Olson, A.J. Automated Docking Using a Lamarckian Genetic Algorithm and an Empirical Binding Free Energy Function. J. Comput. Chem. 1998, 19, 1639–1662. [Google Scholar] [CrossRef]
  37. Sinlikhitkul, N.; Toochinda, P.; Lawtrakul, L.; Kuropakornpong, P.; Itharat, A. Encapsulation of Plumbagin Using Cyclodextrins to Enhance Plumbagin Stability: Computational Simulation, Preparation, Characterization, and Application. J. Incl. Phenom. Macrocycl. Chem. 2019, 93, 229–243. [Google Scholar] [CrossRef]
  38. Jitapunkul, K.; Toochinda, P.; Lawtrakul, L. Molecular Dynamic Simulation Analysis on the Inclusion Complexation of Plumbagin with β-Cyclodextrin Derivatives in Aqueous Solution. Molecules 2021, 11, 6784. [Google Scholar] [CrossRef] [PubMed]
  39. Chulurks, S.; Jitapunkul, K.; Katanyutanon, S.; Toochinda, P.; Lawtrakul, L. Stability Enhancement and Skin Permeation Application of Nicotine by Forming Inclusion Complex with β-Cyclodextrin and Methyl-β-Cyclodextrin. Sci. Pharm. 2021, 9, 43. [Google Scholar] [CrossRef]
  40. Srihakulung, O.; Maezono, R.; Toochinda, P.; Kongprawechnon, W.; Intarapanich, A.; Lawtrakul, L. Host-Guest Interactions of Plumbagin with β-Cyclodextrin, Dimethyl-β-Cyclodextrin and Hydroxypropyl-β-Cyclodextrin: Semi-Empirical Quantum Mechanical PM6 and PM7 Methods. Sci. Pharm. 2018, 86, 20. [Google Scholar] [CrossRef] [Green Version]
  41. Zhao, Y.; Truhlar, D.G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition Elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Account 2008, 120, 215–241. [Google Scholar] [CrossRef]
  42. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Petersson, G.A.; Nakatsuji, H.; et al. Gaussian 16, Revision C.01; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  43. Cova, T.F.; Cruz, S.M.; Valente, A.J.; Abreu, P.; Marques, J.M.; Pais, A.A. Aggregation of Cyclodextrins: Fundamental Issues and Applications. In Cyclodextrin—A Versatile Ingredient [Internet]; Arora, P., Dhingra, N., Eds.; IntechOpen: London, UK, 2018. [Google Scholar] [CrossRef]
Figure 1. Schematic representation of the chemical structures and the dimensions of the energy minimized molecular conformation: (a) β-CD (C42H70O35) and its position of the wider rim (Head, H) and the narrow rim (Tail, T); (b) CBD (C21H30O2); (c) THC (C21H30O2).
Figure 1. Schematic representation of the chemical structures and the dimensions of the energy minimized molecular conformation: (a) β-CD (C42H70O35) and its position of the wider rim (Head, H) and the narrow rim (Tail, T); (b) CBD (C21H30O2); (c) THC (C21H30O2).
Ijms 24 01525 g001
Figure 2. Schematic representations of the conformations of the 1:1 inclusion complex: (a) β-CD/CBD Conf. I (b) β-CD/CBD Conf. II; (c) β-CD/THC Conf. II.
Figure 2. Schematic representations of the conformations of the 1:1 inclusion complex: (a) β-CD/CBD Conf. I (b) β-CD/CBD Conf. II; (c) β-CD/THC Conf. II.
Ijms 24 01525 g002
Figure 3. M062X/6-31g(d,p) energy-optimized structures of 1:1 inclusion complex: (a) β-CD/CBD Conf. I; (b) β-CD/CBD Conf. II; (c) β-CD/THC Conf. II.
Figure 3. M062X/6-31g(d,p) energy-optimized structures of 1:1 inclusion complex: (a) β-CD/CBD Conf. I; (b) β-CD/CBD Conf. II; (c) β-CD/THC Conf. II.
Ijms 24 01525 g003
Figure 4. Schematic representation of the pairwise initial arrangement of β-CDs: (a) Head-to-Head (HH) configuration; (b) Tail-to-Tail (TT) configuration; (c) Tail-to-Head (TH) configuration.
Figure 4. Schematic representation of the pairwise initial arrangement of β-CDs: (a) Head-to-Head (HH) configuration; (b) Tail-to-Tail (TT) configuration; (c) Tail-to-Head (TH) configuration.
Ijms 24 01525 g004
Figure 5. Schematic representation of five conformations of the 2:1 complex: (a) HH/guest Conf. III; (b) TT/guest Conf. IV; (c) TH/guest Conf. V; (d) TH/THC Conf. VI; (e) TH/THC Conf. VII.
Figure 5. Schematic representation of five conformations of the 2:1 complex: (a) HH/guest Conf. III; (b) TT/guest Conf. IV; (c) TH/guest Conf. V; (d) TH/THC Conf. VI; (e) TH/THC Conf. VII.
Ijms 24 01525 g005
Figure 6. M062X/6-31g(d,p) energy-optimized structure of the 2:1 complex: (a) HH/CBD Conf. III; (b) TT/CBD Conf. IV; (c) TH/CBD Conf. V; (d) HH/THC Conf. III; (e) TT/THC Conf. IV; (f) TH/THC Conf. V; (g) TH/THC Conf. VI; (h) TH/THC Conf. VII.
Figure 6. M062X/6-31g(d,p) energy-optimized structure of the 2:1 complex: (a) HH/CBD Conf. III; (b) TT/CBD Conf. IV; (c) TH/CBD Conf. V; (d) HH/THC Conf. III; (e) TT/THC Conf. IV; (f) TH/THC Conf. V; (g) TH/THC Conf. VI; (h) TH/THC Conf. VII.
Ijms 24 01525 g006
Table 1. The lowest and the average values of free energy of binding (ΔG) of β-CD/CBD and β-CD/THC inclusion complexes (1:1 host-guest ratio), and the frequency of conformations in a cluster obtained from molecular docking calculations at 298.15 K.
Table 1. The lowest and the average values of free energy of binding (ΔG) of β-CD/CBD and β-CD/THC inclusion complexes (1:1 host-guest ratio), and the frequency of conformations in a cluster obtained from molecular docking calculations at 298.15 K.
Host/GuestClusterConformationFrequency (%)ΔG (kcal/mol)
LowestAverage
β-CD/CBD1I41−6.87−6.57
2I46−6.81−6.62
3I12−6.73−6.49
4II1−6.07−6.07
β-CD/THC1II44−6.77−6.73
2II26−6.72−6.64
3II25−6.71−6.68
4II4−6.41−6.35
Table 2. Distance of hydrogen bonds between host (β-CD) and guest (CBD and THC), obtained from M062X/6-31g(d,p) energy-optimized structures.
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 (Å)
β-CD/CBD Conf. IO4(β-CD)n ··· H(O1H-CBD)1.84
O2(β-CD)n+1 ··· H(O3H-CBD)1.99
O4(β-CD)n+2 ··· H(O3H-CBD)2.50
β-CD/CBD Conf. IIO4(β-CD)n ··· H(O1H-CBD)1.90
O4(β-CD)n+3 ··· H(O3H-CBD)1.92
β-CD/THC Conf. IIO4(β-CD)n ··· H(O1H-THC)2.04
O2(β-CD)n+1 ··· H(O1H-THC)2.24
O(THC) ··· H(O3H-β-CD)n-22.85
Table 3. Complexation, deformation, HOMO and LUMO energies and dipole moment (µ) of 1:1 inclusion complex of β-CD/CBD and β-CD/THC calculated by M062X/6-31g(d,p) method.
Table 3. Complexation, deformation, HOMO and LUMO energies and dipole moment (µ) of 1:1 inclusion complex of β-CD/CBD and β-CD/THC calculated by M062X/6-31g(d,p) method.
β-CDCBDTHCβ-CD/CBD
Conf. I
β-CD/CBD
Conf. II
β-CD/THC
Conf. II
µ (debye)3.53321.75151.06174.76513.57663.1873
µsp (guest) (debye) 3.07443.32161.0512
µsp (host) (debye) 3.28813.45263.5170
HOMO (eV)−8.7310−6.9408−6.9302−7.2456−7.1514−7.1133
LUMO (eV)2.25501.26401.25830.96250.80410.8041
Δ|HOMO-LUMO| (eV)10.98608.20488.18848.20807.95557.9174
Energy (hartree) a−4273.8373−968.3414−968.3654−5242.2266−5242.2302−5242.2583
ΔE (kcal/mol) −29.9964−32.3017−34.9201
E DEF ( guest ) (kcal/mol) 12.031911.68362.9349
E DEF ( host ) (kcal/mol) 8.69294.71263.4940
a 1 hartree = 627.5095 kcal/mol.
Table 4. The lowest and the average values of free energy of binding (ΔG) of dimeric β-CDs/CBD and dimeric β-CDs/THC inclusion complexes (2:1 host-guest ratio), and the frequency of conformations in a cluster obtained from molecular docking calculations at 298.15 K.
Table 4. The lowest and the average values of free energy of binding (ΔG) of dimeric β-CDs/CBD and dimeric β-CDs/THC inclusion complexes (2:1 host-guest ratio), and the frequency of conformations in a cluster obtained from molecular docking calculations at 298.15 K.
Host/GuestClusterConformationFrequency (%)ΔG (kcal/mol)
LowestAverage
HH/CBD1III55−8.97−8.63
2III16−8.69−8.48
3III14−8.66−8.50
4III10−8.64−8.44
5III5−8.61−8.42
TT/CBD1IV38−8.98−8.72
2IV61−8.52−8.12
3IV1−7.82−7.82
TH/CBD1V91−7.74−7.50
2V5−7.65−7.40
3V4−7.57−7.38
HH/THC1III51−9.61−9.51
2III43−9.49−9.41
3III6−9.12−9.04
TT/THC1IV59−8.77−8.70
2IV37−8.67−8.58
3IV4−8.23−8.18
TH/THC1V41−8.32−8.12
2VI16−8.30−8.20
3VII4−8.22−8.15
4V33−8.11−8.02
5V6−8.09−8.03
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.
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 (Å)
HH/CBD Conf. IIIO3(β-CD-Top)n ··· H(O1H-CBD)1.99
O4(β-CD-Top)n ··· H(O1H-CBD)3.10
TT/CBD Conf. IVO4(β-CD-Top)n ··· H(O3H-CBD)1.86
TH/CBD Conf. VO4(β -CD-Top)n ··· H(O1H-CBD)1.82
HH/THC Conf. IIIO4(β -CD-Bottom)n ··· H(O1H-THC)2.78
TT/THC Conf. IVO4(β -CD-Top)n ··· H(O1H-THC)2.08
TH/THC Conf. VO4(β -CD-Top)n ··· H(O1H-THC)1.92
TH/THC Conf. VIO4(β -CD-Bottom)n ··· H(O1H-THC)2.42
TH/THC Conf. VIIO4(β -CD-Bottom)n ··· H(O1H-THC)2.22
Table 6. Complexation, deformation, HOMO and LUMO energies, and dipole moment (µ) of 2:1 complexes of β-CD/CBD and β-CD/THC calculated by M062X/6-31g(d,p) method.
Table 6. Complexation, deformation, HOMO and LUMO energies, and dipole moment (µ) of 2:1 complexes of β-CD/CBD and β-CD/THC calculated by M062X/6-31g(d,p) method.
HH/CBD
Conf. III
TT/CBD
Conf. IV
TH/CBD
Conf. V
HH/THC
Conf. III
TT/THC
Conf. IV
TH/THC
Conf. V
TH/THC
Conf. VI
TH/THC
Conf. VII
µ (debye)4.94029.06123.81151.85048.41621.79252.24405.8098
µsp (guest) (debye)2.45240.91962.60021.05261.03901.02103.06582.0765
µsp (host) (debye)2.66269.36172.75403.01199.49032.54703.87163.4883
HOMO (eV)−7.3291−7.1884−7.4964−7.4328−7.1522−7.3013−7.6812−7.6439
LUMO (eV)0.67050.83460.49910.44840.90940.63180.23920.2963
Δ|HOMO-LUMO| (eV)7.99968.02307.99557.88128.06167.93327.92047.9403
Energy (hartree) a−9516.2118−9516.1492−9516.2214−9516.2239−9516.1777−9516.2558−9516.2402−9516.2290
ΔE (kcal/mol)−122.7734−83.5332−128.7802−115.3329−86.3642−135.3628−125.5592−118.542
E DEF ( guest ) (kcal/mol)10.01638.23488.94267.40274.73772.28928.836612.7880
E DEF ( host ) (kcal/mol)−70.2202−37.4730−95.2610−72.5865−38.3873−96.9119−88.7794−87.8827
a 1 hartree = 627.5095 kcal/mol.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Triamchaisri, N.; Toochinda, P.; Lawtrakul, L. 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. Int. J. Mol. Sci. 2023, 24, 1525. https://doi.org/10.3390/ijms24021525

AMA Style

Triamchaisri N, Toochinda P, Lawtrakul L. 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. International Journal of Molecular Sciences. 2023; 24(2):1525. https://doi.org/10.3390/ijms24021525

Chicago/Turabian Style

Triamchaisri, Nat, Pisanu Toochinda, and Luckhana Lawtrakul. 2023. "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" International Journal of Molecular Sciences 24, no. 2: 1525. https://doi.org/10.3390/ijms24021525

APA Style

Triamchaisri, N., Toochinda, P., & Lawtrakul, L. (2023). 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. International Journal of Molecular Sciences, 24(2), 1525. https://doi.org/10.3390/ijms24021525

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop