Bayan Al Tbakhia,
Hamdi Nsairatb,
Walhan Alshaerc,
Abdulfattah Al-Kadashd,
Wissam Helala,
Lubna Alrawashdehe,
Anthony Dayf,
Khaleel I. Assafg,
Rola Hassouneha,
Fadwa Odeh*a and
Abeer Al Bawab*ah
aDepartment of Chemistry, School of Science, The University of Jordan, Amman 11942, Jordan. E-mail: f.odeh@ju.edu.jo; Tel: +962-979250409
bPharmacological and Diagnostic Research Centre, Faculty of Pharmacy, Al-Ahliyya Amman University, Amman 19328, Jordan
cCell Therapy Center, The University of Jordan, Amman 11942, Jordan
dSchool of Medicine, The University of Jordan, Amman 11942, Jordan
eDepartment of Chemistry, Faculty of Science, The Hashemite University, P. O. Box 330127, Zarqa 13133, Jordan
fSchool of Science, The University of New South Wales, Canberra, ACT 2600, Australia
gDepartment of Chemistry, Al-Balqa Applied University, Al-Salt 19117, Jordan
hHamdi Mango Center for Scientific Research, The University of Jordan, Amman 11942, Jordan. E-mail: drabeer@ju.edu.jo; Tel: +962-796661601
First published on 8th March 2022
This study aimed to clarify the physico-chemical properties of cucurbit[7]uril (CB[7]) and cinnamaldehyde (Cinn) inclusion complexes (CB[7]–Cinn) and their resulting antitumor activity. CB[7]–Cinn inclusion complexes were prepared by a simple experimental approach and fully characterized for their stoichiometry, formation constant, particle size and morphology. Quantum chemical calculations were performed to elucidate the stable molecular structures of the inclusion complexes and their precursors and to investigate the probable stoichiometry and direction of interaction using three different DFT functionals at the 6-31G(d,p) basis set. The UV-vis spectrophotometric titrations as well as the Job plot, based on 1H NMR spectroscopy, suggested 1:1 and 1:2 stoichiometries of CB[7]:Cinn. The formation constants of the complexes were calculated using Benesi–Hildebrand equations and non-linear fittings. Moreover, the theoretical calculations confirmed the potential formation of 1:1 and 1:2 stoichiometries and clarify the orientation of binding from the Cinn phenyl moiety. The nanoparticles' TEM images showed a crystal-like spherical shape, smooth surface, with a small tendency to agglomerate. CB[7]–Cinn inclusion complexes were analyzed for their antitumor activity against MDA-MB-231 breast cancer and U-87 glioblastoma cell lines. The IC50 values were calculated after 72 hours of incubation with different concentrations of CB[7]–Cinn inclusion complexes and compared to free Cinn and free CB[7]. The IC50 values for free Cinn and CB[7]–Cinn inclusion complexes were 240.17 ± 32.46 μM and 260.47 ± 20.83 μM against U-87 cells and 85.93 ± 3.35 μM and 176.3 ± 7.79 μM against MDA-MB-231 cells, respectively, despite the enhanced aqueous solubility. No significant cytotoxicity was noticed for the free CB[7].
Cinn complexes with cyclodextrins and their derivatives were synthesized and investigated.12–17 However, other macrocyclic compound, such as cucurbit[7]uril (Fig. 1b), may act as drug carriers.
Cucurbit[n]urils (CB[n]), a new family of synthetic macrocyclic hosting agents composed of glycoluril (Fig. 1b) repeated units to form a molecular container.18 The CB[n] family members have common characteristics, such as a hydrophobic cavity and two identical opening hydrophilic entrances. In addition, the varying cavity and portal sizes lead to the formation of inclusion complexes with different molecules through a combination of ion-dipole, hydrogen bonding, and hydrophobic interactions.19 Toxicity studies proved that CB[n]s are biocompatible with low toxicity in the intravenous and oral administration routes.20 The host–guest features of CB[n] have a wide range of applications in nanomachines, self-assembly, drug delivery, the removal of wastes, and in chromatography.19 Moreover, the reversibility of the host–guest complexation process allows full or partial binding of the biological or medicinal activity of the complexed pharmaceutical agents.8
Among CB[n]s, CB[6] and CB[8] are almost insoluble in water, whereas CB[7] has a moderate solubility (2 × 10−2 mol L−1), which is similar to the commonly used complexing host, β-cyclodextrin (with a solubility of 1.6 × 10−2 mol L−1 in water). Moreover, the suitable size of cavity (7.3 Å diameter),21 symmetrical portals (5.4 Å diameter) and the 9.1 Å height of CB[7] make it suitable for binding certain drug molecules, such as Cinn which has a length of 8.3 Å and a height of 4.9 Å.10
Cinn is a planar molecule consisting of one aromatic ring (phenyl) connected with an aliphatic arm of three carbon atoms ending with an aldehyde. The encapsulation of Cinn into host molecules (such as β-cyclodextrin) is usually through the non-polar part of the molecules (the phenyl and aryl),14 which would enhance its water solubility, protect it from oxidation, increase its stability and may-as a potential drug-improve its activity.
In this study, we prepared and characterized CB[7]–Cinn inclusion complex using optical and spectroscopic techniques in addition to electron microscopic imaging. Moreover, a theoretical evaluation of the CB[7]–Cinn inclusion complex stoichiometry and energy measurements was performed. The anti-tumor activity of CB[7]–Cinn inclusion complex against MDA-231 and U87 cancer cell lines was also investigated.
The stoichiometry and the formation constants of CB[7]–Cinn inclusion complex were determined using the Benesi–Hildebrand equations24 based on the data from UV-vis spectrophotometry. Eqn (1) was used for 1:1 complexation and eqn (2) for 1:2 complexation.25–28
(1) |
(2) |
According to the design of Benesi–Hildebrand experiment, Z is the increasing concentration of either the host or the guest, Yo is absorbance (Ao) of Cinn. Y is absorbance (A) with a particular concentration of either the host or the guest, Y′ is absorbance (A′) at the maximum concentration of either the host or the guest. The formation constant (K) is calculated from the Benesi–Hildebrand equation,26 after plotting (1/Y − Yo) against (1/[Z]) or (1/[Z]2):
(3) |
The continuous variation method (Job's method) was also used to determine the complexation stoichiometry24,29 using 1H NMR spectroscopy. The experiments were performed using solutions of equimolar concentrations of the Cinn and CB[7]. The samples were prepared by mixing different volumes of these two solutions so that the total concentration remains constant and the molar fraction of Cinn or CB[7], χ varies in the range 0–1. The chemical shift differences in the presence of CB[7] in respect to the value for the free Cinn, was determined. The variation of the chemical shifts (Δδ) multiplied by Cinn or CB[7] molar ratio, is plotted versus Cinn or CB[7] mole fraction (χ). The value of χ for which the plot presents the maximum deviation gives the stoichiometry of the inclusion complex (χ = 0.5 for 1:1 or 2:2 ratio complexes; χ = 0.33 for 1:2 ratio complexes).30
TEM analysis was accomplished by placing 7 μL of the diluted aqueous suspensions (1% v/v) on 300-mesh formvar copper grids (SPI supplies, USA) and left to dry overnight. The loaded grids were then stained using 2% uranyl acetate dye solution. Excess liquid on the grid was removed with filter paper and the grid was allowed to dry before being observed under Versa 3D (FEI, Netherlands) operated at 30 kV.
The solvent effects were taken into account by using the conductor-like polarizable continuum model (CPCM) method with water as a solvent. Frequency calculations were performed on the optimized geometries at the same level of theory. All computed vibrational transitions have no imaginary frequency implying that each optimized geometry is located at the global minimum point on the potential energy surface. In the second step, an optimization procedure was performed in order to elucidate the stability of the 1:1 complex as a function of the direction of the aldehyde arm toward the macrocycle cavity at constant distances between both compounds. Subsequently, the 1:1 complex potential energy surface was modeled by freezing all geometry parameters, while scanning the single point energy of the complex as a function of the distance between Cinn and the center of the CB[7] macrocycle cavity, with intervals of 0.250 Å. For this step, three different hybrid functionals were used, namely: B3LYP, PBE0, and M06-2X. Finally, a second Cinn molecule was incorporated into the 1:1 complex and a new potential energy curve was computed by freezing all geometry parameters, while scanning the single point energy of the complex as a function of the distance between the second Cinn and the center of the CB[7] cavity, with intervals of 0.250 (Å). All DFT calculations were applied to both sides of Cinn compound, i.e. sides of inclusion in CB[7] host (phenyl side and aldehyde side).
A total of 5 × 103 cells per well were seeded in 96-well plates (TPP, Switzerland) at 37 °C and 5% CO2.33 After 24 hours, cells were treated with different concentrations of CB[7]–Cinn inclusion complexes and their free precursors, Cinn and CB[7] up to 900 μM. The cell growth was analyzed using tetrazolium dye (MTT; 3[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide) assay33,34 after 72 hours. Untreated cells were considered as a negative control and were used as blank to calibrate the spectrophotometer, then the absorbance was read at 570 nm using a using Glomax plate reader (Promega, USA). IC50 values were calculated using GraphPad Prism 6 (GraphPad Software Inc., USA). All experiments were performed in triplicate.
Fig. 2 1H-NMR spectrum of (a) CB[7] and (b) Cinn, in D2O at room temperature, with the protons assigned. |
The formation of the CB[7]–Cinn inclusion complex was clearly evidenced by the chemical shift changes experienced by Cinn protons upon the addition of increasing amounts of CB[7] as shown in Fig. 3. The upfield shift of Cinn protons in the inclusion complex compared to the free Cinn is an indication of a more shielded and restricted microenvironment of the Cinn molecules due to complexation.35 It is clear that the phenyl protons experienced a dramatic change both in shape and chemical shift (Fig. 3) indicating a full inclusion of Cinn in CB[7] from the benzene ring as expected. Both H2 and H3 protons underwent a considerable shift, even though it was less than the shift of the phenyl protons, whereas H1 (the aldehyde proton) witnessed no shift at all.
From the 1H-NMR spectra, it is obvious that the guest (Cinn) enters the CB[7] cavity from the phenyl side. These results are in good agreement with previous studies dealing with aromatic compound having the methylene side chain as reported by Sindelar and coworkers.36 They confirmed, using 1H-NMR titration, that all anthraquinone unit protons are located inside the CB[7] cavity while the methylene protons positioned outside CB[7].36
Fig. 4 Benesi–Hildebrand plots of (a) 1/ΔA versus 1/[CB[7]] and (b) 1/ΔA versus 1/[CB[7]]2 for CB[7]–Cinn inclusion complex using UV-vis spectrophotometric titration. |
The Benesi–Hildebrand plots from UV-vis spectrophotometric spectroscopic data indicate the two possibilities of the formation of complexes with different stoichiometries (1:1 and 1:2 and/or 2:1 of the CB[7]–Cinn inclusion complex). Furthermore, our results are in good consistency with previous studies. For example, Koner et al. reported a binding constant values for the neutral form of albendazole into CB[7] to be 6.5 × 104 M−1 in the UV titrations with CB[7].38
Continuous variation method (Job's method) using 1H-NMR titration data was also employed to determine the complexation stoichiometry. Job's plots are shown in Fig. 5a and b with respect to Cinn and CB[7], respectively, where the mole fraction of either the host or guest (χ) was varied from 0.0 to 1.0 and plotted against χ multiplied by the chemical shift difference of H3 proton of CB[7] (Fig. 6). As shown in Job's plot, two maximum peaks were obtained around χ = 0.5. These two parabola sections may point to the probability of a 1:1 and 1:2 inclusion complex between Cinn and CB[7] as evidenced in Benesi–Hildebrand plots.39,40 This may be attributed to the symmetrical portal sides of CB[7] and the small size of Cinn that led to fast equilibrium of the physical interaction between the complex precursors. The same observation was reported by Romano et al. with Cinn–α-cyclodextrin inclusion complex.41
Fig. 5 Job's plot of CB[7]–Cinn inclusion complex with respect to (a) Cinn and (b) CB[7], where χ is the mole fraction. |
Fig. 6 1H NMR stacked spectra in D2O solvent of Cinn with CB[7] at different mole fractions according to Job's plot method. |
The possibility of the formation of different inclusion complexes is clearly suggested by 1H NMR data. Fig. 6 shows the 1H NMR spectrum of CB[7] in the presence of different concentrations of Cinn. As the Cinn concentration increases, the ‘a’ protons of CB[7] start to broaden with an upfield shift suggesting the presence of different species including empty CB[7], 1:1 and 2:1 Cinn:CB[7] complexes. At a ratio of Cinn:CB[7] of 4:1, a shifted doublet exists suggesting a predominant 2:1 complex. The same trend is observed with ‘c’ protons of CB[7], a singlet at 5.39 ppm (no Cinn present). Upon the addition of Cinn, the signal starts to split in the Cinn:CB[7] ratio from 0.25:1 to 1.5:1, suggesting the presence of different species. However, when the concentration of Cinn is increased to reach 4:1 ratio, a singlet at 5.34 ppm, suggesting a predominant 2:1 complex. These results suggest that an equilibrium is present between free CB[7], a 1:1 complex (I) and 2:1 complex with two possibilities (II and III) as shown in Fig. 7. At a large excess of Cinn, a 2:1 complex is dominant.
Fig. 8 Average size distribution by intensity with different time point for 1:1 inclusion complex nanoparticles of CB[7]–Cinn. |
The charging the particles on the grid in TEM experiment to set and adhere to the sample separate the agglomerated particles into smaller one.42
Fig. 10 A proposed structure for inclusion of Cinn into CB[7] with different orientations (a) aldehyde side and (b) phenyl side (by Gaussian view 09). |
The analysis of the potential energy surfaces corresponding to aldehyde closing-up of the Cinn–CB[7] complex using the PBE0, M06-2X and B3LYP functional with 6-31G(d,p) basis set adapted similar energy behavior (Fig. 11). These energy diagrams showed unstable energy maxima at almost 3 Å before the 7 Å dissociation energy equilibrium distance. This result decreases the probability of formation of an inclusion complex from the aldehyde side of Cinn due to the repulsion interactions between CB[7] glycoluril carbonyl group with aldehyde side of Cinn which is in total agreement with 1H-NMR data that showed no alteration in shape or chemical shift of the aldehyde proton (Fig. 4a).
Clear stable energy minimums were observed at 3 Å before the 7 distance (Å) dissociation energy equilibrium distance obtained from the different DFT functional caculations when considering the phenyl side orientation (Fig. 12). Similar energy curves were produced by the three DFT functional with a small variation in the calculated stabilization energy. These results indicate a strong attraction energy between the Cinn phenyl side and the hydrophobic cavity of CB[7] which is again in total agreement with 1H-NMR data (Fig. 4a).
A comparison of the theoretical calculations of minimum energy, r equilibrium and stabilization energy for inclusion of Cinn into CB[7] with different orientations, phenyl and aldehyde side, are summarize in Table 1. The inclusion complex formation is more favorable from the Cinn phenyl side at 3 Å due to the hydrophobic interaction, where it needs more close distance (0.5–0.75 Å) of Cinn to form inclusion complex from the aldehyde side by the hydrogen bonding of the phenyl group with the glycoluril carbonyl of CB[7].
Method | Minimum energy (a.u) | r equilibrium (Å) | Dissociation energy (a.u) | r equilibrium (Å) | Stabilization energy (Diss-min) (a.u) |
---|---|---|---|---|---|
Phenyl side | |||||
B3LYP | −4635.090656 | 3.00 Å | −4635.086844 | 7.00 Å | 0.003817 |
PBEO | −4630.146811 | 3.00 Å | −4630.141837 | 7.00 Å | 0.004974 |
M06-2X | −4633.449633 | 3.00 Å | −4633.443572 | 7.00 Å | 0.006061 |
Aldehyde side | |||||
B3LYP | −4635.092211 | 0.75 Å | −4635.086678 | 7.00 Å | 0.005533 |
PBE0 | −4630.149158 | 0.75 Å | −4630.141671 | 7.00 Å | 0.007486 |
M06-2X | −4633.456408 | 0.50 Å | −4633.443431 | 7.00 Å | 0.012977 |
The highest stabilization energy was theoretically achieved using M06-2X functionals in both Cinn orientations. This indicates that DFT M06-2X functional is the most suitable for Cinn and CB[7] interaction to build up the potential energy curve of the inclusion complex.
The last computational experiment was performed to investigate the theoretical ability of CB[7] to complex two molecules of Cinn from the phenyl side, a 1:2 stoichiometry, (Fig. 13a) using the DFT, B3LYP functional with 6-31G(d,p) basis set and CPMC/water solvent. The potential energy curve obtained from the CB[7]:Cinn, 1:2 stoichiometry, inclusion complex at the indicated level of theory and basis set showed more stable minima than observed for 1:1 stoichiometry (Fig. 13b). The stabilization energies for 1:1 and 1:2 stoichiometry were 10.02 kJ mol−1 and 11.32 kJ mol−1, respectively, at 3.00 Å distance, at the same level of theory and quantum parameters.
Our study confirmed that free CB[7] has no significant toxicity and CB[7]–Cinn inclusion complex and the free Cinn showed antiproliferative effect against both cell lines. Almost no change in IC50 between free Cinn and CB[7]–Cinn inclusion complex against U-87 cells, while the cytotoxicity was decreased by half with MDA-MB-231 cells. The slight variation in the cytotoxicity of CB[7]–Cinn inclusion complex in the two cell lines related to varied characteristics of different cell lines and may be due to the probability of the two types of stoichiometry that may increase or decrease the concentration of the complexed Cinn. In addition, the CB[7] new complex with hydrophobic Cinn increases the aqueous solubility of the drug and/or drug permeability which could enhance Cinn bioactivity.
The anti-tumor effect of Cinn has been widely investigated against different cancers both in vitro and in vivo.8 Lee et al., 1999 reported the antiproliferative effect of Cinn against 29 cancer cell lines in vitro.43 Several reports described that Cinn exerts the anti-tumor effect through induction apoptosis and cell cycle arrest.44–47 Moreover, the combination of Cinn with other chemotherapies showed a higher anti-tumor effect compared to single drug use.47,48 A further study to investigate the combination of CB[7]–Cinn inclusion complex with other chemotherapies in vitro and in vivo can be of great value.
Our study suggests that CB[7]–Cinn inclusion complex can be an effective as free Cinn against MDA-MB-231 and U-87 cancer cell lines that provide IC50 values 176.3 ± 7.79 μM, 260.47 ± 20.83 μM, respectively. This study suggests that CB[7]-based pharmaceutical formulation is a promising drug delivery system that may improve the stability and bioavailability of hydrophobic drugs in vitro and in vivo. Moreover, the investigation of the synergistic effect between Cinn and other chemotherapies can provide a valuable combination therapy through decreasing the required doses of the cytotoxic drugs and may reduce the unwanted side-effects.
This journal is © The Royal Society of Chemistry 2022 |