Open Access Article
Sherif Ashraf Fahmy
a,
Fortuna Ponteb,
Iten M. Fawzyc,
Emilia Sicilia
b and
Hassan Mohamed El-Said Azzazy
*a
aDepartment of Chemistry, School of Sciences & Engineering, The American University in Cairo, AUC Avenue, PO. Box 74, New Cairo 11835, Egypt. E-mail: sheriffahmy@aucegypt.edu; hazzazy@aucegypt.edu; Fax: +20 2 2795 7565; Tel: +20 2 2615 2559
bDepartment of Chemistry and Chemical Technologies, University of Calabria, Arcavacata di Rende, 87036, Italy. E-mail: fortuna.ponte@unical.it; emilia.sicilia@unical.it
cPharmaceutical Chemistry Department, Faculty of Pharmaceutical Sciences and Pharmaceutical Industries, Future University in Egypt, Cairo 12311, Egypt. E-mail: iten.mamdouh@fue.edu.eg
First published on 14th July 2021
p-Sulfonatocalix[n]arenes have shown excellent potential for accommodating chemotherapeutic drugs through host–guest complexation and enhancing their anticancer activity. Betaine has been reported to exert an anticancer effect at high concentrations. In order to increase its concentration in cancer cells, we have complexed it with p-SC4, which releases its content in an acidic environment typical of cancer tissue. In this work, a host–guest complex of the chemically stable, natural, and safe active methyl donor (betaine) and p-sulfonatocalix[4]arenes (p-SC4) was designed and characterized using 1H NMR, UV, Job's plot analysis, DFT calculations, and molecular modeling for use in cancer therapeutics. The peak amplitude of the prepared host–guest complexes was linearly proportional to the concentration of betaine in the range of 1.0 × 10−5 M−1 to 2.5 × 10−4 M−1. The reaction stoichiometry between p-SC4 and betaine in the formed complex was 1
:
1. The stability constant for the complex is 8.9 × 104 M−1 which corresponds to a complexation free energy of −6.74 kcal mol−1. Complexation between betaine and p-SC4 was found to involve the insertion of the trimethylammonium group of betaine into the p-SC4 cavity, as supported by the experimental data. The complex displayed enhanced cytotoxic activities against breast adenocarcinoma cells (MCF-7) and cervical cancer cells (HeLa) compared to free betaine. In conclusion, the host–guest complexation of betaine with p-SC4 increases its concentration in cancer cells, which warrants further investigation for cancer therapy.
Betaine, trimethylglycine, is a chemically stable, natural, and safe (at a dose up to 15 g per day) active methyl donor that supports the standard DNA methylation patterns.10–12 Betaine is found in many dietary sources such as shrimps, wheat germ, wheat bran, spinach, and sugar beets.12 Many studies reported betaine's ability to prevent and treat cancer, in vitro, via different mechanisms. One study reported the ability of betaine to induce tumor cell apoptosis via altering caspase proteins.13 Betaine was found to affect pro-inflammatory cytokines levels, such as tumor necrosis factor-alpha inside various tumor cells.14 Another study reported betaine's capability to stabilize multiple tumor suppressors and proto-oncogenes' methylation patterns and regulate their expression.15 Researchers reported the cytotoxic activity of different betaine doses on human prostate cancer cells (DU-145).16 However, high doses of betaine are required to exert its anticancer action.17
The use of various delivery vehicles to reformulate natural anticancer agents is a promising approach that may improve their uptake into target tumor cells while simultaneously providing a shield that protects these agents from degradation in the circulation.18 Many studies reported the use of various delivery vehicles to accommodate several anticancer agents.19–24 The host–guest complexation chemistry (macrocycles) has recently attracted much attention as promising drug delivery vehicles.25,26 Many studies reported the possible use of supramolecular structures as hosts for different antineoplastic drugs to enhance their water solubility, chemical stability and improve their selective uptake into tumor cells.26–29 Particularly, calix[n]arenes (n = 4, 6, and 8) signify a crucial class of macrocycles that have been widely utilized in drug delivery. Calixarenes are ideal host molecules formed of phenolic units linked by methylene bridges comprising an upper rim with para-substituent of a phenolic ring, a lower rim with a phenolic hydroxyl group, and a hydrophobic π electron-rich core cavity. This unique structure of calix[n]arenes permits their facile functionalization using sulfonate groups via their phenolic rings' para-position, forming the water-soluble derivative (>0.1 mol L−1), para-sulfonatocalix[n]arenes.26–29 para-Sulfonatocalix[n]arenes are nontoxic to human cells at in vivo doses up to 103 μg kg−1. Furthermore, in aqueous media, para-sulfonatocalix[n]arenes demonstrated an ability to host different guest molecules in their cavities via inclusion complexation.30,31
In this study, the complexation between p-SC4 and betaine in an aqueous solution was studied for potential application in cancer therapeutics. The proposed complex has been characterized by UV, 1H NMR spectroscopy, and computational studies. The stoichiometry and binding constant were computed utilizing Job's plot (continuous variation method). Density functional theory (DFT) calculations were carried out to characterize the structure of the host–guest adduct. Moreover, the anticancer activities of the complex against breast adenocarcinoma (MCF-7) and cervical cancer (HeLa) were studied.
The geometry optimizations were performed using the double-zeta 6-31G(d,p) basis set for all the involved atoms, except for oxygen atoms, for which a diffuse function was included. The solvation effects were computed using the continuum solvation model based on density (SMD),37 as implemented in Gaussian 09. To confirm the minimum character of the intercepted structures, vibrational frequency analysis was carried out at the same level of theory, and frequency calculations, within the harmonic approximation, were used to calculate the Gibbs free energies for the inclusion of betaine (G) into calixarenes (H).
The inclusion Gibbs free energies in solution ΔGsol, in implicit water, was calculated as the difference between the free energy in solvent (GHG complexsol) of the complex and the free energies of the separated species in solvent (GHsol and GGsol):
| ΔGHG complexsol = GHG complexsol − (GHsol + GGsol) |
Counterpoise correction calculations, as formulated by Boys and Bernardi38 were carried out to estimate the Basis Set Superposition Error (BSSE).
The characterization of molecular interactions responsible for the formation of the betaine–p-SC4 adducts was performed using the theory of atoms in molecules (AIM) suggested by Bader,38,39 using the AIMAll program.40
The analysis of the electron density ρ(r) around the interacting species allows identification of the chemical bonding nature of the optimized geometries.
The values of ρ(r) and its derivatives at the CPs (Critical Points; points where the gradient of the density is zero, ∇ρ(r) = 0) provides information about the nature and the strength of the established interactions. To characterize the CPs, it is necessary to explore the second-order derivatives of ρ(r). Only nine derivatives are possible, useful to generate a real and symmetric Hessian matrix of ρ(r). The diagonalization of the matrix, using a unitary transformation, allows finding the corresponding eigenvalues (λn). CPs are marked by the rank ω and the signature σ, respectively. The rank of the diagonalized matrix is given by the number of non-zero eigenvalues (non-zero curvatures of ρ(r) at the CPs). The signature of a CP, σ, is defined as the algebraic sum of the signs of the eigenvalues (signs of curvatures of ρ(r) at the CPs). CP with ω = 3 and σ = +1, corresponding to (3, +1) CP, has two positive curvatures, and ρ is a minimum at CP in the plane defined by their corresponding axes while at CP along the third axis, perpendicular to this plane, ρ is a maximum. The CP (3, −1), with ω = 3 and σ = −1, highlights that electronic charge density is located within a chemical bond between the involved atoms. Such a point is called the bond critical point (BCP).
:
1 molar ratio mixture of p-SC4 and betaine are presented in Fig. S1A–C;† respectively. The 1H NMR findings showed that the protons of betaine had exhibited remarkable upfield shifts upon adding p-SC4 (Table 1 and Fig. S1†) relative to those of the betaine alone. These shifts are attributed to the shielding effect of the aromatic rings of p-SC4. Previous studies reported the direct relationship between the degree of the chemical shift of a guest molecule's protons and their depth of penetration inside the hollow cavity of the host molecule, designating the formation of host–guest inclusion complexes.28,29,41,42 The most significant upfield shift of ∼0.58 ppm was observed for the protons of the methyl groups bonded to the quaternary ammonium head of betaine, Hc (Table 1 and Fig. S1†) showing that they possessed the most significant shielding effect. These indications suggest the possible insertion of betaine inside the p-SC4 cavity through its N-terminal moiety. Moreover, the signals of the aryl and methylene bridge protons (Ha and Hb; respectively) of p-SC4, exhibited some chemical shifts upon betaine addition, which indicates the formation of host–guest inclusion complex (Table 1 and Fig. S1†).27–29
:
1 M ratio mixture in D2O. Proton labels are demonstrated in Fig. S1a
| Proton signals | Individual molecules | 1 : 1 molar ratio |
|---|---|---|
| a Ha and Hb represent the aryl and methylene bridge protons in p-SC4, respectively. Hc and Hd represent the methyl and methylene protons in betaine, respectively. | ||
| Ha | 7.39 | 7.52 |
| Hb | 3.84 | 3.65 |
| Hc | 3.08 | 2.50 |
| Hd | 3.71 | 3.64 |
![]() | ||
| Fig. 1 Absorbance spectra of 0.2 mM betaine, 0.2 mM p-SC4, and a mixture encompassing 0.2 mM of each betaine and p-SC4 all in aqueous media. | ||
The hyperchromic shift observed with increasing betaine concentrations was obtained by dividing the zero-order spectra of the prepared mixtures by the spectrum of p-SC4 where the first derivative of the ratio spectra was computed using a scaling factor of 10 and Δλ = 4 nm.27–29 The values of the peak amplitudes of the first derivative of the ratio spectra for the mixtures were then measured at 264 nm. A plot of the obtained peak amplitudes and the equivalent concentrations of betaine is presented in Fig. 2. The idyllic linear relationships illustrated in Fig. 2 proposed that the hyperchromic shift detected in betaine/p-SC4 is due to the formation of host–guest inclusion complex between betaine and p-SC4.
![]() | ||
| Fig. 2 Plot of peak amplitudes at 264 nm attained from the mixtures containing sequentially increasing concentrations of betaine (0.01–0.25 mM) and a fixed concentration of 0.2 mM p-SC4. | ||
The stoichiometry and the stability constants of the host–guest inclusion complex were examined using Job's plot (the method of continuous variation) utilizing the derivative ratio method as described previously.27–29 From Job's plots, we noticed that the maximum amplitudes were spotted at molar fractions of 0.5, signifying a host–guest inclusion complex stoichiometry of 1
:
1 for betaine/p-SC4. A normalized form of this Job's plot for betaine/p-SC4, where each amplitude value, S, was divided by the corresponding maximum amplitude, Smax, is depicted in Fig. 3. The stability constant of the betaine/p-SC4 inclusion complex computed employing methods detailed previously27–29 was 8.9 × 104 M−1 which corresponds to the complexation free energy of −6.74 kcal mol−1. This value agrees well with the stability constants (0.01 × 103 to 1.7 × 105 M−1) previously reported for host–guest inclusion complexes proposed to be applied in drug delivery. Many of them showed enhanced stability, therapeutic activity, and selective targeting to the desired site of action.27–29,43–54
| Cells | In vitro cytotoxic activity (IC50 in μg mL−1) | ||
|---|---|---|---|
| p-SC4 | Betaine | Betaine/p-SC4 | |
| a Data represent the mean ± standard deviation of triplicate. | |||
| MCF-7 | >300 | >220 | 1.6 ± 1.4 |
| HeLa | >300 | >220 | 3.3 ± 1.1 |
![]() | ||
| Fig. 4 B97-D optimized structures of p-SC4: (a) side view and (b) top view. Upper (d1) and lower (d2) rim diameters of the calixarene are reported. | ||
The inclusion of guest molecules in the host cavity, as was shown recently, can be predicted based on their volumes. For p-SC4, the inner cavity volume results to be 652 Å3, while betaine has a volume of 114 Å3. Based on the calculated volumes, the penetration of betaine inside the p-SC4 receptor, in principle, could take place, even if, in the formation of the host–guest inclusion compounds, several aspects must be taken into consideration as illustrated below.
The potential insertion modes, computationally explored, of the betaine into the host molecule are reported in Scheme 1. In particular, the betaine can assume two different possible orientations with respect to the p-SC4 macrocycle. In the arrangement labelled as 1, the positively charged trimethylammonium group points inwards the host cavity, and the carboxylate group points outwards, while the disposition in which the carboxylate group points towards the inside of the cavity and the trimethylammonium group is outwardly oriented is labelled as 2.
The optimized geometries obtained for the formation of host–guest complexes (or adducts) between betaine and p-SC4 receptor are illustrated in Fig. 5. In the same figure, the most relevant distances together with the highlighted Bond Critical Point (BCPs) and the values of the complexation free energies are reported. The energy of complexation was calculated, including the effect of basis set superposition error (BSSE) and entropy change corrections in the solution.
In the optimization process, any constraints were imposed. In the intercepted structure labelled as 1H-G, the betaine penetrates deeper into the host cavity, maximizing the electrostatic interactions. In particular, in this insertion complex, the quaternary ammonium group is located inside the p-SC4 cavity, while the carboxylate group points towards the upper ring establishing hydrogen bond interactions with the protonated oxygen atoms of the two closest sulfonate groups. The data computationally obtained for this inclusion complex reflects the information coming from 1H NMR analysis. The formation of this betaine/p-SC4 adduct leads to an energy gain compared to separate species, and the value of the calculated complexation free energy is −12.8 kcal mol−1.
In the 2H-G adduct, the betaine guest interacts with the cavity of the receptor through the carboxylate group of betaine, and the trimethylammonium group points outwards. In particular, the carboxylate group of the host interacts with all the sulfonate groups of the guest through intermolecular hydrogen bonds, causing the closure of the upper rim, whose diameter decreases by 0.537 Å. The free energy of complexation for the 2H-G complex calculated to be −3.2 kcal mol−1 is smaller than the 1H-G complex.
With the aim to investigate more in-depth the host–guest interactions, the AIM tool was used.39,40 The electron density (ρ) value and its Laplacian ∇2ρ(r) at the bond critical point (BCP) were estimated by AIM calculation to determine the nature and the strength of the interactions between the species forming both 1H-G and 2H-G complexes. In particular, for van der Waals interactions ρ(r) is 10−3 a.u., while for hydrogen bonding ρ(r) is ∼10−2 a.u. or less. For both interactions, the corresponding Laplacian ∇2ρ(r) is positive. Table 3 reports the results obtained from the AIM analysis.
| Conformation | BCP | ρBCP | ∇2ρBCP |
|---|---|---|---|
| 1H-G | 1 | 0.011 | 0.750 |
| 2 | 0.024 | 0.138 | |
| 2H-G | 1 | 0.054 | 0.138 |
| 2 | 0.054 | 0.148 | |
| 3 | 0.025 | 0.077 | |
| 4 | 0.054 | 0.136 |
In both complexes generated by the betaine guest and the calix p-SC4 host, the electron density values fall in the expected range for H-bonds, namely, 0.0102–0.0642 a.u.,39 and for all the BCPs, the Laplacian sign is positive. The strength of these hydrogen interactions in the identified complexes is established by the highest positive values of the Laplacian.
As reported in Fig. 5, two BCPs, 1BCP1 and 1BCP2, respectively, are predicted for the 1H-G complex. The value of Laplacian is larger for BCP1 than for BCP2, indicating a stronger interaction. In the 2H-G complex, as shown in Fig. 5, four critical points indicated as 2BCP1, 2BCP2, 2BCP3, and 2BCP4 are established. All the BCPs correspond to hydrogen bonds and are defined by values of the Laplacian that are positive and of the same magnitude, except for BCP3 which is described by a smaller Laplacian value indicating a lower interaction. The formation of these H-bonds increases the stability of the adduct with respect to the separate species.
:
1 molar ratio. The stability constant of betaine/p-SC4 was computed to be 8.9 × 104 M−1, which corresponds to the free energy of complexation of −6.74 kcal mol−1. The stability of the prepared complex in an aqueous medium supports its application in cancer drug delivery. Both 1H NMR spectroscopy and DFT calculations indicate that betaine was included within the cavity of p-SC4, and the complexation involves the insertion of the trimethylammonium group of the guest into the host cavity. Data from the SRB assay showed that the IC50 of betaine/p-SC4 was 1.6 ± 1.4 and 3.3 ± 1.1 μg mL−1 against MCF-7 and HeLa cells, respectively, compared to individual betaine, which exerted a negligible decline of cell viability.
In summary, the complexation between p-SC4 and betaine was thoroughly investigated experimentally and computationally. To the best of our knowledge, there are no reports on encapsulation of betaine in nanocarriers or supramolecular structures. A remarkable enhancement of the cytotoxic activities of betaine/p-SC4 was observed in two different cancer cell lines as compared to free betaine. This work presents a novel approach for delivering high betaine concentrations into cancer cells and supports the potential use of the designed complex as a possible effective supplement to cancer therapeutics.
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: 10.1039/d1ra04614d |
| This journal is © The Royal Society of Chemistry 2021 |