Next Article in Journal
Structural and Functional Changes of Reconstituted High-Density Lipoprotein (HDL) by Incorporation of α-synuclein: A Potent Antioxidant and Anti-Glycation Activity of α-synuclein and apoA-I in HDL at High Molar Ratio of α-synuclein
Next Article in Special Issue
A Dispersion Corrected DFT Investigation of the Inclusion Complexation of Dexamethasone with β-Cyclodextrin and Molecular Docking Study of Its Potential Activity against COVID-19
Previous Article in Journal
Conjugated Metabolites of Hydroxytyrosol and Tyrosol Contribute to the Maintenance of Nitric Oxide Balance in Human Aortic Endothelial Cells at Physiologically Relevant Concentrations
Previous Article in Special Issue
NHC Polymeric Particles Obtained by Self-Assembly and Click Approach of Calix[4]Arene Amphiphiles as Support for Catalytically Active Pd Nanoclusters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

DFT-D4 Insight into the Inclusion of Amphetamine and Methamphetamine in Cucurbit[7]uril: Energetic, Structural and Biosensing Properties

by
Abdelkarim Litim
1,
Youghourta Belhocine
2,
Tahar Benlecheb
1,*,
Monira Galal Ghoniem
3,
Zoubir Kabouche
1,
Fatima Adam Mohamed Ali
3,*,
Babiker Yagoub Abdulkhair
4,5,
Mahamadou Seydou
6 and
Seyfeddine Rahali
7,*
1
Laboratory of Sensors, Instrumentations and Process (LCIP), University of Abbes Laghrour, Khenchela 40000, Algeria
2
Department of Petrochemical and Process Engineering, Faculty of Technology, 20 August 1955 University of Skikda, El Hadaik Road, Skikda 21000, Algeria
3
Department of Chemistry, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
4
Chemistry Department, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
5
Chemistry Department, Sudan University of Science and Technology (SUST), Khartoum 13311, Sudan
6
Université de Paris, CNRS, ITODYS, F-75013 Paris, France
7
Department of Chemistry, College of Science and Arts, Qassim University, Ar Rass, Saudi Arabia
*
Authors to whom correspondence should be addressed.
Molecules 2021, 26(24), 7479; https://doi.org/10.3390/molecules26247479
Submission received: 13 November 2021 / Revised: 4 December 2021 / Accepted: 8 December 2021 / Published: 10 December 2021
(This article belongs to the Special Issue Synthesis and Molecular Recognition of Macrocyclic Compounds)

Abstract

:
The host–guest interactions of cucurbit[7]uril (CB[7]) as host and amphetamine (AMP), methamphetamine (MET) and their enantiomeric forms (S-form and R-form) as guests were computationally investigated using density functional theory calculations with the recent D4 atomic-charge dependent dispersion corrections. The analysis of energetic, structural and electronic properties with the aid of frontier molecular orbital analysis, charge decomposition analysis (CDA), extended charge decomposition analysis (ECDA) and independent gradient model (IGM) approach allowed to characterize the host–guest interactions in the studied systems. Energetic results indicate the formation of stable non-covalent complexes where R-AMP@CB[7] and S-AMP@CB[7] are more stable thermodynamically than R-MET@CB[7] and S-MET@CB[7] in gas phase while the reverse is true in water solvent. Based on structural analysis, a recognition mechanism is proposed, which suggests that the synergistic effect of van der Waals forces, ion–dipole interactions, intermolecular charge transfer interactions and intermolecular hydrogen bonding is responsible for the stabilization of the complexes. The geometries of the complexes obtained theoretically are in good agreement with the X-ray experimental structures and indicate that the phenyl ring of amphetamine and methamphetamine is deeply buried into the cavity of CB[7] through hydrophobic interactions while the ammonium group remains outside the cavity to establish hydrogen bonds with the portal oxygen atoms of CB[7].

1. Introduction

In 1887, the Romanian chemist, Lazăr Edeleanu discovered synthetic psychostimulants which are known as amphetamines [1,2,3]. This class of drugs is used as pharmacotherapy for obesity, narcolepsy and attention-deficit/hyperactivity disorder (ADHD) [4]. Despite their therapeutic benefits, amphetamines have been associated with several adult health issues and an increasing concern about the growing number of babies who are subjected to these drugs due to maternal use during the pregnancy period [5]. These stimulant drugs work by stimulating the central nervous system [6,7], they increase symptoms such as alertness and focus and decrease others such as tiredness and appetite [8].
As reported by the World Health Organization, amphetamine-type stimulants (ATS) belong to a class of drugs that include amphetamine (AMP), methamphetamine (MET) and some stimulants which are structurally related [9].
Methamphetamine is commonly referred to as “meth” or “crystal meth”. It is a very strong stimulant of the central nervous system (psychostimulant) which is derived from phenethylamine and amphetamine group of substances. Similarly to amphetamine, methamphetamine effects result in an increased level of wakefulness, attention and energy; however, taking high doses will expose users to higher levels of euphoria, sexual desire and self-esteem [10,11]. Methamphetamine is the most used illegal stimulant after Cannabis [12]. It is also the most produced among ATS around the world [13].
Immunochromatography and enzyme-linked immunosorbent assay (ELISA) are the most-used techniques to identify ATS using oral fluid [14]. Classical ELISA assays are laboratory tests that are based on the rivalry between the drug present in the sample and the same drug tagged with the enzyme (added to the system), for points of attachment of antibodies immobilized in the wells of a plate. Urine toxicology assessments deliver sensitive and accurate measures of various types of drugs; however, certain drugs such as ATS are confined or constrained to a limited timeline because of the fast metabolizing. The identification window can also be dependent on elements such as purity, the usage frequency, the person’s physical state (e.g., hydration degree), and the way the drug was ingested (e.g., oral, intranasal, intravenous injections) [15].
Macrocyclic compounds are important structures and key building blocks in supramolecular chemistry [16,17,18]. Among them, cucurbiturils are a relatively novel family of barrel-shaped macrocyclic hosts that can encapsulate guest molecules in their unique large cavities, which are edged by carbonyl oxygen atoms [19,20,21]. Cucurbiturils can be synthesized in a variety of different morphology sizes, with either 5, 6, 7, 8 or 10 glycoluril subunits [21,22,23,24]. In recent years, cucurbit[n]urils (CB[n]) were used as drug delivery vehicles for neutral molecules and cations due to their high affinity to form host–guest inclusion complexes [25,26,27]. Among the numerous cucurbiturils of different sizes, CB[7] which comprises seven glycoluril units is highly explored in the literature as a successful host molecule [28,29] with decreased toxicity and increased chemical and thermal stability [30,31].
In their study about the development of new point-of-use, portable and wireless drug sensor, Jang et al. showed the high sensitivity of the host molecule CB[7] (combined with organic field-effect transistors) toward ATS, thus providing a viable methodology for the fabrication of selective and sensitive drug sensors [32]. The potential of functionalized CB[7] was reported in some studies as a drug biosensor for the recognition of amphetamine and methamphetamine stimulants [32,33]. It is possible to better detail and quantify this type of interaction using quantum chemistry based on DFT study [34,35,36,37]. Indeed, atomistic calculations are necessary in providing answers on the nature and effect of inter-molecular interactions [38,39,40,41,42].
Herein, we present a theoretical study on the structural, electronic, biosensing and energetic properties of the inclusion of amphetamine and methamphetamine inside the cavity of CB[7]. The present work aims to investigate the stability of the formed complexes and to examine the nature of the intermolecular interactions involved in the complexation process using the new DFT-D4 dispersion correction model.

2. Computational Procedures

The initial molecular structure of CB[7] was taken from the crystal structure determined by X-ray diffraction methods (CCDC 1542694) [32]. Both S and R enantiomers of amphetamine (AMP) and methamphetamine (MET) ((S)-amphetamine CID: 5826; (R)-amphetamine CID: 32893; (S)-methamphetamine CID: 10836 (R)-methamphetamine CID: 36604) were investigated in this study and their chemical 3D structures were retrieved from PubChem database server [43]. Then, protonated forms of AMP and MET were considered and built with Avogadro program [44] assuming that their amino groups are protonated as observed in the solid-state crystal structures of 1:1 inclusion complexes of AMP and MET with CB[7]. The geometries of all initial structures were subjected to prior geometry optimization with ORCA program (version 4.2.0) [45,46] in gas phase using BLYP functional [47,48] in conjunction with the atom-pairwise dispersion correction based on tight-binding partial charges (D4) [49,50] and def2-SVP basis set [51,52] without any symmetry constraints. A geometrical counterpoise correction (gCP) to the def2-SVP for the intra- and inter-molecular basis set superposition error was used [53]. Then, the computational process followed to calculate the most stable configurations (configurations of minimum energy) corresponding to the inclusion complexes of AMP and MET with CB[7] was based on the method developed by Liu and Guo [54]. In addition to Liu and Guo method, other techniques have been developed recently for the theoretical modelling of host–guest binding affinity [55,56,57,58,59,60]. Liu and Guo proposed the semi-empirical methods such as AM1 and PM3, which are fast and could be applied for several positions along the z-axis. We have extended the use of this method to examine accurately the non-covalent interactions (electrostatic interactions, van der Waals interactions, hydrogen bonds…) by applying exclusively DFT approach associated with the most recent dispersion correction functionals (BLYP-D4).
The center of AMP, MET and CB[7] was set as the center of the coordination system, then AMP and MET guests, as well as their enantiomers, were translated along the Z-axis from −10 to +10 Å with 2 Å step as shown in Figure 1 using Gabedit program [61]. Each generated structure in this process was fully optimized at BLYP-D4/def2-SVP-gCP level of theory. Frequency computations were performed on AMP, MET, CB[7] and the most stable inclusion complexes to ensure that the optimized structures are minima on the potential energy surface. The most stable configuration for each complex (R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-MET@CB[7]) was further reoptimized at BLYP-D4/def2-TZVP-gCP in gas phase and aqueous phase. The water solvent effects were considered by applying the conductor-like polarizable continuum model (CPCM) [62,63]. The following Equation (1) was used for the calculation of complexation energies:
ΔEComplexation = EComplex − (Efree guest + Efree CB[7])
where ∆Ecomplexation represents the complexation energy, Ecomplex, Efree guest and Efree CB[7] are respectively the optimized energies of the complex, the free amphetamine or methamphetamine and the free cucurbit[7]uril.
The aqueous phase geometries were used to characterize the nature of the interactions occurring between the guest molecules (AMP and MET) and CB[7] by performing an independent gradient model (IGM) analysis based on Hirshfeld partition of electron density [64]. The non-covalent interactions (NCI) [65,66] maps were generated with Multiwfn software [67] and visualized using the Visual Molecular Dynamics (VMD) program [68]. To quantify the amount of charge transfer between CB[7] and AMP and MET, the charge decomposition analysis (CDA) [69,70] and extended charge decomposition analysis (ECDA) [71] were performed using the Multiwfn program [67].

3. Results and Discussions

3.1. Searching for the Most Stable Complexes

The complexation energies calculated in gas phase as a function of Z coordinate are reported in Table 1. The optimized structures show negative energies of formation, indicating therefore a favorable thermodynamic process. The most stable configurations for R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-MET@CB[7] with respective energies of −337.43, −346.79, −335.39 and −336.52 kJ/mol are located respectively at Z = 0 Å, Z = −6 Å, Z = −6 Å and Z = −4 Å (Table 1).
The S-enantiomer complexes (S-AMP@CB[7] and S-MET@CB[7]) are respectively slightly more stable by −9.36 and −1.13 kJ/mol in the gas phase than their corresponding R-enantiomers (R-AMP@CB[7] and R-MET@CB[7]).
For higher accuracy, we reoptimized each of the four most stable complexes at BLYP-D4/def2-TZVP-gCP level of theory in gas and aqueous phases. The calculated complexation energies are reported in Table 2.
In gas phase, we note that as the basis set improves from def2-SVP to def2-TZVP, the computed complexation energies of R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-MET@CB[7] are −338.85, −349.58, −334.01 and −334.53 kJ/mol, respectively. The improvement of basis set quality does not lead to significant variation of the complexation energies (Table 2).
In aqueous phase, a significant decrease in the complexation energy values was noted (Table 2); however, the results show a similar trend with those obtained in gas phase calculations. Indeed, S-AMP@CB[7] and S-MET@CB[7] are found to be 3.18 and 1.31 kJ/mol more stable than the complexes R-AMP@CB[7] and R-MET@CB[7]. The small difference in their complexation energies suggests the formation of a racemic mixture of AMP and MET enantiomers [72], as observed in the X-ray crystal structures of the host–guest complexes where both forms of R and S enantiomers were found. Based on the frequency calculations, we determined ΔGComplexation at temperature 298K for each complex and compared it with the experimental values. The experimental ΔGComplexation values of AMP@CB[7] and MET@CB[7] are −34.7 and −33.8 kJ/mol, which are very close to the computed values (−45.5 and −35.3 kJ/mol, respectively). Therefore, our calculations results are in good agreement with the experimental data.
It is also clear that complexes of R-AMP@CB[7] and S-AMP@CB[7] are more stable than their analogs R-MET@CB[7] and S-MET@CB[7] in gas phase while in water solvent the stability is reversed, and the latter complexes become more stable. The water solvent significantly decreased the complexation energy, most likely due to the large dipole moments of the complexes in water. Indeed, the calculated dipole moments of R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-MET@CB[7] in water are 18.03, 14.93, 17.32 and 17.88 D, respectively, and are larger than their corresponding values of 9.23, 8.73, 9.26 and 9.19 D in gas phase.
Further, single-point energy calculations were performed on BLYP-D4/def2-TZVP optimized geometries at B3LYP-D4/def2-TZVP level of theory [73], both functionals show the same behavior for the evaluation of the complexation energies as shown in Table S1 (Supplementary Data).
Using a simple thermodynamic study, it is possible to determine which guest interacts most strongly with CB[7]. The idea will be based on the computation of standard free energy change ΔΔG*Sol for the following guest exchange reaction in solution (2):
MET @ CB [ 7 ] ( aq ) + AMP   ( aq ) AMP @ CB [ 7 ] ( aq ) + MET   ( aq )
As shown in Scheme 1, ΔΔG*Sol can be determined from its components by introducing a thermodynamic cycle using Equation (3):
Δ Δ G sol * = Δ Δ G g 0 + Δ G sol ( AMP @ CB [ 7 ] ) * + Δ G sol ( MET ) * Δ G sol ( MET @ CB [ 7 ] ) * Δ G sol ( AMP ) *
where ΔΔG0g is the change of free energy in the gas-phase reaction; ΔG*sol(AMP@CB[7]), ΔG*sol(MET), ΔG*sol(MET@CB[7]) and ΔG*sol(AMP) are, respectively, the solvation free energies of AMP@CB[7], (MET), MET@CB[7] and (AMP) species in water.
MET and AMP species exhibit several isomers, their Gibbs free energies have been evaluated using a Boltzmann distribution according to the following relation (4):
G [ A ] 0 = R T ln i [ A ] e G i 0 / R T
where [A] emphasizes calculation over the population of the two enantiomers of each guest. Gas phase free energies of all species were estimated at T = 298.15 K and p = 1 atm through the computation of their energy and geometry at BLYP-D4/def2-SVP-gCP level of theory. Gibbs free energies of aqueous solvation were computed for all species using the conductor-like polarizable continuum model (CPCM) using BLYP-D4/def2-TZVP-gCP. Table 3 reports calculated values of all thermodynamic parameters used in the current study. From this Table, we have noticed that the standard free energy change ΔΔG*sol is 12.7 kJ/mol. The positive value of ΔΔG*sol indicates that the exchange reaction moves in the direction of the formation of MET@CB[7] complex, therefore, CB[7] forms complex with methamphetamine more favorably than with amphetamine.

3.2. Geometries of the Most Stable Complexes of AMP and MET with CB[7]

The geometries of fully optimized structures of all inclusion complexes calculated at BLYP-D4-def2-TZVP-gCP level in water solvent (Supplementary Table S2) are visualized using Jmol viewer [74] and illustrated in Figure 2.
DFT simulations show that AMP and MET are encapsulated inside CB[7]. Indeed, the benzyl moiety of AMP and MET is deeply inserted into CB[7] due to hydrophobic interactions between the CB[7] cavity and the phenyl ring, whereas the ammonium group of AMP and the methylammonium group of MET are found outside the cavity, at the top of CB[7] leaning to the closest carbonyl groups at the portals, in agreement with the experimental structures. The calculated shortest N…O bond lengths for AMP@CB[7] and MET@CB[7] are both equal to 2.85 Å and are close to the experimental bond values of 2.76 and 2.78 Å, respectively, involving effective ion–dipole interactions between CB[7] and the ammonium groups.

3.3. Analysis of Non-Covalent Interactions

The characterization of non-covalent interactions (NCI) is of great importance in molecular recognition systems. We have employed NCI analysis to understand the nature of intermolecular interactions occurring in AMP@CB[7] and MET@CB[7] complexes through the isosurface map of the independent gradient model (IGM) approach. IGM Plots for all compounds were generated with an isovalue of 0.005 a.u. and shown in Figure 3, where blue and green colors represent respectively strong attractive interactions and weak van der Waals interactions.
As observed in Figure 3, the IGM isosurfaces indicate the presence of weak van der Waals interactions (green colored fields) and hydrogen bonds (blue disks) between the hydrogen atoms of the ammonium groups of AMP and MET and the lone pairs on the oxygen atoms of CB[7]. The calculated hydrogen bonds (N—H···O) of the studied complexes (Supplementary Figure S1) are within the range 1.87–2.72 Å. Both weak van der Waals interactions and intermolecular hydrogen bonding [34] stabilize the formation of AMP@CB[7] and MET@CB[7] complexes.

3.4. Biosensing Properties

The real-time application of drug biosensing based on organic field-effect transistors sensors is a promising platform for detecting illegal drugs such as amphetamine and methamphetamine. Jang et al. showed that OFET based sensors functionalized with CB[7] derivatives for the binding of amphetamine and methamphetamine exhibit highly sensitive and selective sensing properties in water, urine and physiological buffers [32]. Analysis of frontier molecular orbitals (the HOMO and LUMO) is important for the computational design of drug sensors.
According to the energies of HOMO and LUMO levels reported in Table 4, a decrease of the HOMO–LUMO gaps is observed upon the formation of the inclusion complexes, due to the lowering of the LUMO levels, whereas the energies of the HOMO levels remain almost unchanged, suggesting a charge transfer between CB[7] and AMP and MET, therefore, the variation of electric conductivity can be influenced owing to the charges transfer between AMP and MET and CB[7]. These results indicate that miniaturized electrical biosensors devices can be used to detect AMP and MET in urine or biological fluid samples.
Frontier molecular orbital (HOMO and LUMO) plots for the studied complexes are shown in Figure 4. The analysis of HOMO and LUMO plots reveals that the HOMOs are almost exclusively spread over the CB[7] while the LUMOs are mainly localized on AMP and MET and almost entirely on their aromatic ring moieties.

3.5. Charge Decomposition and Extended Charge Decomposition Analysis

The CDA and ECDA analysis methods can provide insight into the charge transfer (donation and back donation) between fragments in non-covalent complexes. The calculated results indicate that the CDA and ECDA analysis of R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-AMP@CB[7] follow the same trend as shown in Table 4. An electron donation (d) from CB[7] to guest (AMP or MET) in the range [0.21–0.23] was obtained from CDA analysis, while the back-donation (b) occurring from the guests (AMP or MET) to the CB[7] is in the range [0.042−0.071]. The ECDA analysis confirms the net charge transfer obtained by AMP or MET from CB[7] with an approximate value of 0.2 (Table 5), indicating that intermolecular charge transfer is involved in stabilizing the formed complexes.

4. Conclusions

DFT calculations including D4 dispersion correction were carried out to investigate the inclusion process between CB[7] and amphetamine-type stimulant drugs. The conclusions are as followed:
According to the computed complexation energies, the inclusion process of amphetamine and methamphetamine drugs in CB[7] is favorable in gas phase and in water solution.
In gas phase, R-AMP@CB[7] and S-AMP@CB[7] are thermodynamically more favored than R-MET@CB[7] and S-MET@CB[7], while the latter are more favored than the former in water solvent.
The presence of a mixture of R and S enantiomers was evidenced by DFT-D4 calculations, with a small energy difference in favor of the S enantiomers.
X-ray data support the theoretically predicted inclusion conformations of amphetamine and methamphetamine in the CB[7] with the phenyl ring moiety deeply inserted inside the CB[7] cavity and the ammonium group outside the cavity.
Structural analysis indicated that van der Waals forces, intermolecular hydrogen bonds, charge transfer interactions and ion-dipole interactions between CB[7] and amphetamine and methamphetamine play a major role in the process of complex stabilization.
Frontier molecular orbital analysis of all complexes revealed that the LUMO is shifted to lower energy and stabilized upon complexation, the resulting decrease in the HOMO-LUMO gap leads to a variation of the conductivity, making them suitable for biosensing applications.
These theoretical results highlight the high affinity of cucurbiturils for cationic guests through cooperative non-covalent interactions and their potential as stabilizing hosts and a good sensing candidates for the detection of drugs in real samples for practical applications.

Supplementary Materials

The following are available online. The intermolecular hydrogen bonds between the hydrogen atoms of the ammonium groups of AMP and MET and the oxygen atoms of CB[7] are shown in Supplementary Figure S1. Table S1 represents the B3LYP-D4/ def2-TZVP single-point calculations (in italic) on BLYP-D4/def2-TZVP (in bold) optimized structure. The atomic coordinates of optimized structures of S-AMP@CB[7], R-AMP@CB[7], S-MET@CB[7] and R-MET@CB[7] obtained at BLYP-D4-def2-TZVP-gCP level of theory in water solvent are listed in Supplementary Table S2.

Author Contributions

Conceptualization, Y.B. and T.B.; methodology, A.L.; software, Y.B. and S.R.; validation, T.B., M.S. and M.G.G.; formal analysis, Y.B., S.R.; investigation, F.A.M.A.; resources, B.Y.A.; data curation, Y.B.; writing—original draft preparation, A.L.; writing—review and editing, Y.B., Z.K. and M.S.; visualization, S.R.; supervision, Y.B.; project administration, F.A.M.A.; funding acquisition, Y.B. and M.G.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Scientific Research at Imam Mohammad Ibn Saud Islamic University, grant number RG-21-09-75.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

Authors extend their appreciation to the Deanship of Scientific Research at Imam Mohammad Ibn Saud Islamic University for funding this work through Research Group no RG-21-09-75.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

No samples are available Samples of the compounds are not available from the authors.

References

  1. Heal, D.J.; Smith, S.L.; Gosden, J.; Nutt, D.J. Amphetamine, past and present–a pharmacological and clinical perspective. J. Psychopharmacol. 2013, 27, 479–496. [Google Scholar] [CrossRef] [Green Version]
  2. Sitte, H.H.; Freissmuth, M. Amphetamines, new psychoactive drugs and the monoamine transporter cycle. Trends Pharmacol. Sci. 2015, 36, 41–50. [Google Scholar] [CrossRef] [Green Version]
  3. Ruan, X.; Xing, L.; Peng, J.; Li, S.; Song, Y.; Sun, Q. A simplified fabric phase sorptive extraction method for the determination of amphetamine drugs in water samples using liquid chromatography-mass spectrometry. RSC Adv. 2020, 10, 10854–10866. [Google Scholar] [CrossRef]
  4. Faraone, S.V. The pharmacology of amphetamine and methylphenidate: Relevance to the neurobiology of attention-deficit/hyperactivity disorder and other psychiatric comorbidities. Neurosci. Biobehav. Rev. 2018, 87, 255–270. [Google Scholar] [CrossRef] [PubMed]
  5. Oei, J.L.; Kingsbury, A.; Dhawan, A.; Burns, L.; Feller, J.M.; Clews, S.; Falconer, J.; Abdel-Latif, M. Amphetamines, the pregnant woman and her children: A review. J. Perinatol. 2012, 32, 737–747. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bodzon-Kulakowska, A.; Antolak, A.; Drabik, A.; Marszalek-Grabska, M.; Kotlińska, J.; Suder, P. Brain lipidomic changes after morphine, cocaine and amphetamine administration—DESI—MS imaging study. Biochim. Biophys. Acta (BBA)-Mol. Cell Biol. Lipids 2017, 1862, 686–691. [Google Scholar] [CrossRef]
  7. Berman, S.; O’Neill, J.; Fears, S.; Bartzokis, G.; London, E.D. Abuse of amphetamines and structural abnormalities in brain. Ann. N. Y. Acad. Sci. 2008, 1141, 195. [Google Scholar] [CrossRef] [Green Version]
  8. Avois, L.; Robinson, N.; Saudan, C.; Baume, N.; Mangin, P.; Saugy, M. Central nervous system stimulants and sport practice. Br. J. Sports Med. 2006, 40, i16–i20. [Google Scholar] [CrossRef] [Green Version]
  9. World Health Organization (WHO). Management of substance abuse—amphetamine-type stimulants. Available online: https://www.who.int/home/cms-decommissioning (accessed on 30 March 2021).
  10. Cruickshank, C.C.; Dyer, K.R. A review of the clinical pharmacology of methamphetamine. Addiction 2009, 104, 1085–1099. [Google Scholar] [CrossRef]
  11. Courtney, K.E.; Ray, L.A. Methamphetamine: An update on epidemiology, pharmacology, clinical phenomenology, and treatment literature. Drug Alcohol Depend. 2014, 143, 11–21. [Google Scholar] [CrossRef] [Green Version]
  12. Smart, R.G.; Ogborne, A.C. Drug use and drinking among students in 36 countries. Addict. Behav. 2000, 25, 455–460. [Google Scholar] [CrossRef]
  13. Chomchai, C.; Chomchai, S. Global patterns of methamphetamine use. Curr. Opin. Psychiatry 2015, 28, 269–274. [Google Scholar] [CrossRef] [PubMed]
  14. Nieddu, M.; Burrai, L.; Baralla, E.; Pasciu, V.; Varoni, M.V.; Briguglio, I.; Demontis, M.P.; Boatto, G. ELISA detection of 30 new amphetamine designer drugs in whole blood, urine and oral fluid using Neogen®“amphetamine” and “methamphetamine/MDMA” kits. J. Anal. Toxicol. 2016, 40, 492–497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Harrison, L. The validity of self-reported drug use in survey research: An overview and critique of research methods. NIDA Res Monogr 1997, 167, 17–36. [Google Scholar]
  16. Yu, J.; Qi, D.; Li, J. Design, synthesis and applications of responsive macrocycles. Commun. Chem. 2020, 3, 1–14. [Google Scholar] [CrossRef]
  17. Sakagami, R.; Saito, Y.; Mori, R.; Satake, M.; Okayasu, M.; Kikkawa, S.; Hikawa, H.; Azumaya, I. Cylindrical macrocyclic compounds synthesized by connecting two bowl-shaped calix[3]aramide moieties: Structures and chiroptical properties. RSC Adv. 2020, 10, 34549–34555. [Google Scholar] [CrossRef]
  18. Hussain, A.; Yousuf, S.; Mukherjee, D. Importance and synthesis of benzannulated medium-sized and macrocyclic rings (BMRs). RSC Adv. 2014, 4, 43241–43257. [Google Scholar] [CrossRef]
  19. Jie, K.; Zhou, Y.; Yao, Y.; Huang, F. Macrocyclic amphiphiles. Chem. Soc. Rev. 2015, 44, 3568–3587. [Google Scholar] [CrossRef] [PubMed]
  20. Lee, J.W.; Samal, S.; Selvapalam, N.; Kim, H.-J.; Kim, K. Cucurbituril homologues and derivatives: New opportunities in supramolecular chemistry. Acc. Chem. Res. 2003, 36, 621–630. [Google Scholar] [CrossRef]
  21. Masson, E.; Ling, X.; Joseph, R.; Kyeremeh-Mensah, L.; Lu, X. Cucurbituril chemistry: A tale of supramolecular success. Rsc Adv. 2012, 2, 1213–1247. [Google Scholar] [CrossRef]
  22. Kim, J.; Jung, I.-S.; Kim, S.-Y.; Lee, E.; Kang, J.-K.; Sakamoto, S.; Yamaguchi, K.; Kim, K. New cucurbituril homologues: Syntheses, isolation, characterization, and X-ray crystal structures of cucurbit[n]uril (n = 5, 7, and 8). J. Am. Chem. Soc. 2000, 122, 540–541. [Google Scholar] [CrossRef]
  23. Kircheva, N.; Dobrev, S.; Dasheva, L.; Koleva, I.; Nikolova, V.; Angelova, S.; Dudev, T. Complexation of biologically essential (mono-and divalent) metal cations to cucurbiturils: A DFT/SMD evaluation of the key factors governing the host–guest recognition. RSC Adv. 2020, 10, 28139–28147. [Google Scholar] [CrossRef]
  24. Liu, S.; Zavalij, P.Y.; Isaacs, L. Cucurbit[10]uril. J. Am. Chem. Soc. 2005, 127, 16798–16799. [Google Scholar] [CrossRef] [PubMed]
  25. Assaf, K.I.; Nau, W.M. Cucurbiturils: From synthesis to high-affinity binding and catalysis. Chem. Soc. Rev. 2015, 44, 394–418. [Google Scholar] [CrossRef] [Green Version]
  26. Honda, Y.; Hanaya, T.; Sueishi, Y. Inclusion complexation abilities of cucurbit[6]uril for various aromatic amines in the presence of alkali metal cations. J. Incl. Phenom. Macrocycl. Chem. 2017, 88, 253–257. [Google Scholar] [CrossRef]
  27. Zhang, S.; Grimm, L.; Miskolczy, Z.; Biczók, L.; Biedermann, F.; Nau, W.M. Binding affinities of cucurbit[n]urils with cations. Chem. Commun. 2019, 55, 14131–14134. [Google Scholar] [CrossRef] [PubMed]
  28. Jeon, Y.J.; Kim, S.-Y.; Ko, Y.H.; Sakamoto, S.; Yamaguchi, K.; Kim, K. Novel molecular drug carrier: Encapsulation of oxaliplatin in cucurbit[7]uril and its effects on stability and reactivity of the drug. Org. Biomol. Chem. 2005, 3, 2122–2125. [Google Scholar] [CrossRef] [Green Version]
  29. Saleh, N.I.; Khaleel, A.; Al-Dmour, H.; Al-Hindawi, B.; Yakushenko, E. Host–guest complexes of cucurbit[7]uril with albendazole in solid state. J. Therm. Anal. Calorim. 2013, 111, 385–392. [Google Scholar] [CrossRef]
  30. Uzunova, V.D.; Cullinane, C.; Brix, K.; Nau, W.M.; Day, A.I. Toxicity of cucurbit[7]uril and cucurbit[8]uril: An exploratory in vitro and in vivo study. Org. Biomol. Chem. 2010, 8, 2037–2042. [Google Scholar] [CrossRef]
  31. Zhang, P.; Qin, S.; Qi, M.; Fu, R. Cucurbit[n]urils as a new class of stationary phases for gas chromatographic separations. J. Chromatogr. A 2014, 1334, 139–148. [Google Scholar] [CrossRef]
  32. Jang, Y.; Jang, M.; Kim, H.; Lee, S.J.; Jin, E.; Koo, J.Y.; Hwang, I.-C.; Kim, Y.; Ko, Y.H.; Hwang, I. Point-of-use detection of amphetamine-type stimulants with host-molecule-functionalized organic transistors. Chem 2017, 3, 641–651. [Google Scholar] [CrossRef] [Green Version]
  33. Du, X.; Hao, H.; Qin, A.; Tang, B.Z. Highly sensitive chemosensor for detection of methamphetamine by the combination of AIE luminogen and cucurbit[7]uril. Dye. Pigment. 2020, 180, 108413. [Google Scholar] [CrossRef]
  34. Assaba, I.M.; Rahali, S.; Belhocine, Y.; Allal, H. Inclusion complexation of chloroquine with α and β-cyclodextrin: Theoretical insights from the new B97-3c composite method. J. Mol. Struct. 2021, 1227, 129696. [Google Scholar] [CrossRef]
  35. Allal, H.; Belhocine, Y.; Rahali, S.; Damous, M.; Ammouchi, N. Structural, electronic, and energetic investigations of acrolein adsorption on B 36 borophene nanosheet: A dispersion-corrected DFT insight. J. Mol. Modeling 2020, 26, 1–12. [Google Scholar] [CrossRef] [PubMed]
  36. Bouhadiba, A.; Rahali, S.; Belhocine, Y.; Allal, H.; Nouar, L.; Rahim, M. Structural and energetic investigation on the host/guest inclusion process of benzyl isothiocyanate into β-cyclodextrin using dispersion-corrected DFT calculations. Carbohydr. Res. 2020, 491, 107980. [Google Scholar] [CrossRef] [PubMed]
  37. Rahali, S.; Belhocine, Y.; Seydou, M.; Maurel, F.; Tangour, B. Multiscale study of the structure and hydrogen storage capacity of an aluminum metal-organic framework. Int. J. Hydrog. Energy 2017, 42, 15271–15282. [Google Scholar] [CrossRef]
  38. Liu, Y.; Zhao, J.; Li, F.; Chen, Z. Appropriate description of intermolecular interactions in the methane hydrates: An assessment of DFT methods. J. Comput. Chem. 2013, 34, 121–131. [Google Scholar] [CrossRef]
  39. Fang, D.; Piquemal, J.-P.; Liu, S.; Cisneros, G.A. DFT-steric-based energy decomposition analysis of intermolecular interactions. Theor. Chem. Acc. 2014, 133, 1484. [Google Scholar] [CrossRef]
  40. Rahali, S.; Belhocine, Y.; Touzeau, J.; Tangour, B.; Maurel, F.; Seydou, M. Balance between physical and chemical interactions of second-row diatomic molecules with graphene sheet. Superlattices Microstruct. 2017, 102, 45–55. [Google Scholar] [CrossRef]
  41. Rahali, S.; Seydou, M.; Belhocine, Y.; Maurel, F.; Tangour, B. First-principles investigation of hydrogen storage on lead (II)-based metal-organic framework. Int. J. Hydrogen Energy 2016, 41, 2711–2719. [Google Scholar] [CrossRef]
  42. Imane, D.; Leila, N.; Fatiha, M.; Abdelkrim, G.; Mouna, C.; Ismahan, L.; Abdelazize, B.; Brahim, H. Investigation of intermolecular interactions in inclusion complexes of pyroquilon with cucurbit[n]urils (n = 7, 8) using DFT-D3 correction dispersion. J. Mol. Liq. 2020, 309, 113233. [Google Scholar] [CrossRef]
  43. PubChem. Available online: https://www.pubchem.ncbi.nlm.nih.gov/ (accessed on 3 June 2021).
  44. Hanwell, M.D.; Curtis, D.E.; Lonie, D.C.; Vandermeersch, T.; Zurek, E.; Hutchison, G.R. Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. J. Cheminformatics 2012, 4, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Neese, F. The ORCA program system. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 73–78. [Google Scholar] [CrossRef]
  46. Neese, F. Software update: The ORCA program system, version 4.0. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2018, 8, e1327. [Google Scholar] [CrossRef]
  47. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098. [Google Scholar] [CrossRef] [PubMed]
  48. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 1988, 37, 785. [Google Scholar] [CrossRef] [Green Version]
  49. Caldeweyher, E.; Bannwarth, C.; Grimme, S. Extension of the D3 dispersion coefficient model. J. Chem. Phys. 2017, 147, 034112. [Google Scholar] [CrossRef]
  50. Caldeweyher, E.; Mewes, J.-M.; Ehlert, S.; Grimme, S. Extension and evaluation of the D4 London-dispersion model for periodic systems. Phys. Chem. Chem. Phys. 2020, 22, 8499–8512. [Google Scholar] [CrossRef]
  51. Weigend, F.; Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. Phys. Chem. Chem. Phys. 2005, 7, 3297–3305. [Google Scholar] [CrossRef]
  52. Weigend, F. Accurate Coulomb-fitting basis sets for H to Rn. Phys. Chem. Chem. Phys. 2006, 8, 1057–1065. [Google Scholar] [CrossRef]
  53. Kruse, H.; Grimme, S. A geometrical correction for the inter-and intra-molecular basis set superposition error in Hartree-Fock and density functional theory calculations for large systems. J. Chem. Phys. 2012, 136, 04B613. [Google Scholar] [CrossRef] [Green Version]
  54. Liu, L.; Guo, Q.-X. Use of quantum chemical methods to study cyclodextrin chemistry. J. Incl. Phenom. Macrocycl. Chem. 2004, 50, 95–103. [Google Scholar] [CrossRef]
  55. Sun, Z.; Huai, Z.; He, Q.; Liu, Z. A General Picture of Cucurbit[8]uril Host–Guest Binding. J. Chem. Inf. Modeling 2021. [Google Scholar] [CrossRef] [PubMed]
  56. Gonzalez, D.; Macaya, L.; Vöhringer-Martinez, E. Molecular Environment-Specific Atomic Charges Improve Binding Affinity Predictions of SAMPL5 Host–Guest Systems. J. Chem. Inf. Modeling 2021, 61, 4462–4474. [Google Scholar] [CrossRef]
  57. Huai, Z.; Yang, H.; Li, X.; Sun, Z. SAMPL7 TrimerTrip host–guest binding affinities from extensive alchemical and end-point free energy calculations. J. Comput.-Aided Mol. Des. 2021, 35, 117–129. [Google Scholar] [CrossRef]
  58. Yin, J.; Henriksen, N.M.; Slochower, D.R.; Shirts, M.R.; Chiu, M.W.; Mobley, D.L.; Gilson, M.K. Overview of the SAMPL5 Host-Guest Challenge: Are We Doing Better? J. Comput.-Aided Mol. Des. 2017, 31, 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Rizzi, A.; Murkli, S.; McNeill, J.N.; Yao, W.; Sullivan, M.; Gilson, M.K.; Chiu, M.W.; Isaacs, L.; Gibb, B.C.; Mobley, D.L.; et al. Overview of the SAMPL6 host– guest binding affinity prediction challenge. J. Comput.-Aided Mol. Des. 2018, 32, 937–963. [Google Scholar] [CrossRef]
  60. Sun, Z. SAMPL7 TrimerTrip Host-Guest Binding Poses and Binding Affinities from Spherical-Coordinates-Biased Simulations. J. Comput.-Aided Mol. Des. 2021, 35, 105–115. [Google Scholar] [CrossRef] [PubMed]
  61. Allouche, A.R. Gabedit—A graphical user interface for computational chemistry softwares. J. Comput. Chem. 2011, 32, 174–182. [Google Scholar] [CrossRef] [PubMed]
  62. Cossi, M.; Rega, N.; Scalmani, G.; Barone, V. Energies, structures, and electronic properties of molecules in solution with the C-PCM solvation model. J. Comput. Chem. 2003, 24, 669–681. [Google Scholar] [CrossRef]
  63. Takano, Y.; Houk, K. Benchmarking the conductor-like polarizable continuum model (CPCM) for aqueous solvation free energies of neutral and ionic organic molecules. J. Chem. Theory Comput. 2005, 1, 70–77. [Google Scholar] [CrossRef] [PubMed]
  64. Lefebvre, C.; Khartabil, H.; Boisson, J.-C.; Contreras-García, J.; Piquemal, J.-P.; Hénon, E. The independent gradient model: A new approach for probing strong and weak interactions in molecules from wave function calculations. ChemPhysChem 2018, 19, 724–735. [Google Scholar] [CrossRef] [PubMed]
  65. Contreras-García, J.; Boto, R.A.; Izquierdo-Ruiz, F.; Reva, I.; Woller, T.; Alonso, M. A benchmark for the non-covalent interaction (NCI) index or… is it really all in the geometry? Theor. Chem. Acc. 2016, 135, 242. [Google Scholar] [CrossRef] [Green Version]
  66. Al-Hamdani, Y.S.; Tkatchenko, A. Understanding non-covalent interactions in larger molecular complexes from first principles. J. Chem. Phys. 2019, 150, 010901. [Google Scholar] [CrossRef] [Green Version]
  67. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 2012, 33, 580–592. [Google Scholar] [CrossRef]
  68. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  69. Dapprich, S.; Frenking, G. Investigation of donor-acceptor interactions: A charge decomposition analysis using fragment molecular orbitals. J. Phys. Chem. 1995, 99, 9352–9362. [Google Scholar] [CrossRef]
  70. Xiao, M.; Lu, T. Generalized charge decomposition analysis (GCDA) method. J. Adv. Phys. Chem 2015, 4, 111–124. [Google Scholar] [CrossRef]
  71. Gorelsky, S.I.; Ghosh, S.; Solomon, E.I. Mechanism of N2O reduction by the μ4-S tetranuclear CuZ cluster of nitrous oxide reductase. J. Am. Chem. Soc. 2006, 128, 278–290. [Google Scholar] [CrossRef]
  72. Guo, X.; Wang, Z.; Zuo, L.; Zhou, Z.; Guo, X.; Sun, T. Quantitative prediction of enantioseparation using β-cyclodextrin derivatives as chiral selectors in capillary electrophoresis. Analyst 2014, 139, 6511–6519. [Google Scholar] [CrossRef]
  73. Rezac, J. Non-Covalent Interactions Atlas Benchmark Data Sets 2: Hydrogen Bonding in an Extended Chemical Space. J. Chem. Theory Comput. 2020, 16, 6305–6316. [Google Scholar] [CrossRef] [PubMed]
  74. Jmol: An Open-Source Java Viewer for Chemical Structures in 3D. Available online: http://www.jmol.org/ (accessed on 6 July 2021).
Figure 1. Docking orientations of S-AMP@CB[7] from −10 Å to +10 Å. Color code: carbon, grey; nitrogen, blue; oxygen, red; hydrogen, white.
Figure 1. Docking orientations of S-AMP@CB[7] from −10 Å to +10 Å. Color code: carbon, grey; nitrogen, blue; oxygen, red; hydrogen, white.
Molecules 26 07479 g001
Scheme 1. Thermodynamic cycle used to calculate the free energy of reaction 1 from its components.
Scheme 1. Thermodynamic cycle used to calculate the free energy of reaction 1 from its components.
Molecules 26 07479 sch001
Figure 2. Structures of the most stable conformations optimized at BLYP-D4-def2-TZVP-gCP level in aqueous phase.
Figure 2. Structures of the most stable conformations optimized at BLYP-D4-def2-TZVP-gCP level in aqueous phase.
Molecules 26 07479 g002
Figure 3. IGM isosurfaces plots (isovalue = 0.005 a.u.) of the studied complexes.
Figure 3. IGM isosurfaces plots (isovalue = 0.005 a.u.) of the studied complexes.
Molecules 26 07479 g003
Figure 4. HOMO and LUMO plots of R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-AMP@CB[7].
Figure 4. HOMO and LUMO plots of R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-AMP@CB[7].
Molecules 26 07479 g004
Table 1. Calculated complexation energies (in kJ/mol) of CB[7] with R-AMP, S-AMP, R-MET and S-MET at BLYP-D4/def2-SVP-gCP level.
Table 1. Calculated complexation energies (in kJ/mol) of CB[7] with R-AMP, S-AMP, R-MET and S-MET at BLYP-D4/def2-SVP-gCP level.
Docking ConfigurationsR-AMP@CB[7]S-AMP@CB[7]R-MET@CB[7]S-MET@CB[7]
−10 Å−301.08−323.76−334.75−284.33
−8 Å−299.69−295.34−269.57−331.66
−6 Å−324.09−346.79−335.39−320.60
−4 Å−324.33−346.60−334.80−336.52
−2 Å−337.34−337.14−334.73−317.27
0 Å−337.43−337.02−334.71−317.46
+2 Å−320.76−337.54−334.64−322.59
+4 Å−300.51−300.90−282.05−276.98
+6 Å−300.58−300.91−277.80−291.80
+8 Å−300.60−289.05−278.07−287.32
+10 Å−300.71−295.21−278.39−286.65
Table 2. Calculated complexation energies and standard Gibbs free energy (in kJ/mol) of CB[7] with AMP and MET at BLYP-D4/def2-TZVP-gCP level.
Table 2. Calculated complexation energies and standard Gibbs free energy (in kJ/mol) of CB[7] with AMP and MET at BLYP-D4/def2-TZVP-gCP level.
Complex ΔEComplexation (kJ/mol)ΔGComplexation (kJ/mol) (a)
Gas Phase Aqueous PhaseExperimentalCalculated
R-AMP@CB[7]−338.85−73.23−34.7−45.5
S-AMP@CB[7]−349.58−76.41
R-MET@CB[7]−334.01−84.01−33.8−35.3
S-MET@CB[7]−334.53−85.32
(a) ΔGComplexation values are calculated in aqueous phase at 298 k.
Table 3. Thermodynamic parameters used for the computation of standard free energy change.
Table 3. Thermodynamic parameters used for the computation of standard free energy change.
SpeciesΔG0gas (a.u.)ΔG*sol (a.u.)ΔΔG0gas (kJ/mol)ΔΔG*sol (kJ/mol)
AMP−405.16178−405.26800−10.212.7
MET−444.38724−444.48477
AMP@CB[7]−4612.20227−4612.38314
MET@CB[7]−4651.42384−4651.604737
Table 4. Calculated HOMO and LUMO energy levels (eV) and HOMO–LUMO energy gaps (eV) of CB[7], R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-MET@CB[7].
Table 4. Calculated HOMO and LUMO energy levels (eV) and HOMO–LUMO energy gaps (eV) of CB[7], R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-MET@CB[7].
System EHOMO (eV)ELUMO (eV)∆ELUMO-HOMO (eV)
CB[7]−5.62−0.525.10
R-AMP@CB[7]−5.68−0.954.73
S-AMP@CB[7]−5.68−1.134.55
R-MET@CB[7]−5.67−0.954.72
S-MET@CB[7]−5.68−0.984.70
Table 5. Charge decomposition (CDA) analysis and extend charge decomposition analysis (ECDA) of CB[7], R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-MET@CB[7].
Table 5. Charge decomposition (CDA) analysis and extend charge decomposition analysis (ECDA) of CB[7], R-AMP@CB[7], S-AMP@CB[7], R-MET@CB[7] and S-MET@CB[7].
CDAECDA
Systemdbd-brNet Electrons Obtained by AMP or MET
R-AMP@CB[7]0.2320.0710.161−0.3200.196
S-AMP@CB[7]0.2090.0420.167−0.3060.206
R-MET@CB[7]0.2320.0690.163−0.3230.199
S-MET@CB[7]0.2310.0650.166−0.3350.203
d: number of electrons donated from CB[7] to guest (AMP or MET); b: number of electrons back-donated from the guest (AMP or MET) to CB[7]; r: number of electrons involved in repulsive polarization.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Litim, A.; Belhocine, Y.; Benlecheb, T.; Ghoniem, M.G.; Kabouche, Z.; Ali, F.A.M.; Abdulkhair, B.Y.; Seydou, M.; Rahali, S. DFT-D4 Insight into the Inclusion of Amphetamine and Methamphetamine in Cucurbit[7]uril: Energetic, Structural and Biosensing Properties. Molecules 2021, 26, 7479. https://doi.org/10.3390/molecules26247479

AMA Style

Litim A, Belhocine Y, Benlecheb T, Ghoniem MG, Kabouche Z, Ali FAM, Abdulkhair BY, Seydou M, Rahali S. DFT-D4 Insight into the Inclusion of Amphetamine and Methamphetamine in Cucurbit[7]uril: Energetic, Structural and Biosensing Properties. Molecules. 2021; 26(24):7479. https://doi.org/10.3390/molecules26247479

Chicago/Turabian Style

Litim, Abdelkarim, Youghourta Belhocine, Tahar Benlecheb, Monira Galal Ghoniem, Zoubir Kabouche, Fatima Adam Mohamed Ali, Babiker Yagoub Abdulkhair, Mahamadou Seydou, and Seyfeddine Rahali. 2021. "DFT-D4 Insight into the Inclusion of Amphetamine and Methamphetamine in Cucurbit[7]uril: Energetic, Structural and Biosensing Properties" Molecules 26, no. 24: 7479. https://doi.org/10.3390/molecules26247479

APA Style

Litim, A., Belhocine, Y., Benlecheb, T., Ghoniem, M. G., Kabouche, Z., Ali, F. A. M., Abdulkhair, B. Y., Seydou, M., & Rahali, S. (2021). DFT-D4 Insight into the Inclusion of Amphetamine and Methamphetamine in Cucurbit[7]uril: Energetic, Structural and Biosensing Properties. Molecules, 26(24), 7479. https://doi.org/10.3390/molecules26247479

Article Metrics

Back to TopTop