Open accessResearch articleFirst published online 2021-5
Synthesis,in vitro anticancer activities,and quantum chemical investigations on 1,3- bis -(2-methyl-2-propenyl)benzimidazolium chloride and its Ag(I) complex
1,3-Bis-(2-methyl-2-propenyl)benzimidazolium chloride and its Ag(I) complex are synthesized and the structures are elucidated using spectroscopies techniques. The molecular and crystal structures of the benzimidazolium salt are confirmed by X-ray crystallography. The molecular geometries of the benzimidazolium and its Ag(I) salt are analyzed using the B3LYP functional with the 6–311+G(d,p)/LANL2DZ basis set. The observed Fourier transform infrared and nuclear magnetic resonance isotropic shifts are compared with the calculated values. Besides, the quantum chemical identifiers, significant intramolecular interactions, and molecular electrostatic potential plots are used to show the tendency/site of the chemical reactivity behavior. The three-dimensional Hirshfeld surfaces and the associated two-dimensional fingerprint plots are applied to obtain an insight into the behavior of the interactions in the crystal. Both compounds are tested for their in vitro anticancer activities against DU-145 and MCF-7 cancer cells and L-929 non-cancer cells using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assay.
N-heterocyclic carbenes (NHCs) are aromatic organic compounds that contain at least three carbon and two nitrogen atoms. NHCs were first reported by Öfele, Wanzlick, and Schönherr.1,2 In general, NHCs are not easy to isolate as single carbene monomers due to their decomposition in the presence of water or their dimerization. All approaches to obtain a stable N-heterocyclic carbene had failed, until Arduengo3 reported the first isolated NHC in 1991. Recently, the use of N-heterocyclic carbenes in coordination chemistry has attracted considerable attention due to their large steric demand and excellent electronic properties, namely, high σ-basicity and low π-acidity,4–8 which foster increased activity in catalytic systems.9–16 Benzimidazole-based N-heterocyclic carbenes are stable systems and are the subject of significant interest because of pharmacological activities such as antitumor,17 antibacterial,18 and antifungal.19 Most NHCs are prepared by deprotonation of azolium precursors with a simple metal salt or a strong base.20,21 Also, NHCs can be obtained from imidazolidine sources.22
Halide anions have been the subject of increasing research in investigations of both environmental and supramolecular chemistry, since halides are among the most common anions in natural environments. Among halides, chloride plays a vital role in all body fluids because it is responsible for maintaining acid/base balance, transmitting nerve impulses, and so on.23,24 In addition, water is an essential molecule in all aspects of human life, so these water−chloride interactions are of fundamental importance to be able to understand solution phenomena, materials chemistry, drug design, catalysis, and atmospheric research.25,26 In recent years, monochloride hydrates [Cl(H2O)n]– have been extensively investigated using experimental and theoretical techniques, and they have addressed issues such as stability and structure. In addition, the Pt(0) complexes of benzimidazolylidene carbene ligands have been developed and investigated in terms of their ability to catalyze alkene hydrosilylation reactions.27 Recently, the catalytic activity of Pd(II)28 and Co(II)29 complexes of structurally related benzimidazole-based NHCs have been reported. In a study related to different types of Ag(I)–L complexes, the antiproliferative activity of such metal complexes against cancer cells was investigated and the results showed that they had selectivity toward the human breast cancer cells.30
The aim of this paper is to report the preparation and characterization of benzimidazolium salt 1 and its Ag(I)–NHC complex 2. Characterization was performed using 1H nuclear magnetic resonance (NMR), 13C NMR, and Fourier transform infrared (FTIR) spectroscopies for both compounds. Single-crystal X-ray diffraction was used to elucidate the molecular and crystal structure of the benzimidazolium salt. The X-ray results revealed that the asymmetric unit hosts a cation, a chloride anion, and a water molecule, which are bonded to each other through an intricate hydrogen-bonding network. It was reported that the chloride hydrate structures ([Cl2(H2O)2]2–) in the crystal structure, due to the formation of the hydrogen-bonded cluster, have caused a supermolecular structure.31 Crystallographic studies are useful to investigate the supramolecular assembly of such complexes. Quantum chemical studies were conducted on two compounds to estimate the electronic, spectroscopic, and biological reactivity behavior. With this aim, the geometric structures of the two compounds were predicted, then the assigned NMR shifts and IR assignments of these compounds were compared with the corresponding experimental values. A frontier molecular orbital (FMO) investigation was applied, and molecular electrostatic potential (MEP) diagrams were drawn to estimate/evaluate the chemical reactive behavior/site of two compounds. The in vitro anticancer activity of the benzimidazolium salt 1 and the Ag(I) complex 2 was investigated against DU-145 human prostate cancer cells, MCF-7 human breast cancer cells, and mouse L-929 non-cancer adipose cells from the mouse for 24, 48, and 72 h using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay.32
Results and discussion
Molecular structures of the NHC ligands and Ag(I)–NHC complexes
NHC ligand 1 was obtained by binding two 2-methyl-2-propenyl units to benzimidazole. The synthesis reaction steps are shown in Scheme 1.
Synthesis of NHC ligand 1.
The synthesized NHC ligand 1 and Ag2O were mixed in dichloromethane with the exclusion of the light to give Ag(I)–NHC complex 2 in excellent yield. The synthetic route is shown in Scheme 2.
Synthesis of Ag(I)–NHC complex 2.
The optimized parameters of the compounds were verified by frequency calculations at the B3LYP/6–311+G(d,p)/LANL2DZ are summarized in Table 1. These values were then compared with the experimental values of compound 1 determined by X-ray crystallography. The optimized geometries and atom labeling of each compound are shown in Figure 1.
The selected geometric parameters of compounds 1 and 2.
The experimental parameters are based on the single-crystal X-ray structure of compound 1.
The optimized structures of 1 and 2 at the B3LYP/6–311+G(d,p)/LANL2DZ level in CHCl3.
First, it should be noted that there are small differences between the experimental and computational data because the experimental results were recorded for the solid-state structure of the molecule. From Table 1, the N11–C13 bond length was recorded as 1.321 Å and calculated (for 1 and 2) as 1.338 and 1.357 Å. The N11–C16 single bond lengths characterized by single bonding for 1 and 2 were computed as 1.476 and 1.470 Å, with the corresponding experimental bond length being 1.47 Å. In the literature, the C–N bond length and N–C–N bond angle for benzimidazolium chloride hydrate have been observed at 1.4 Å and 111.1°33 and calculated for the structurally related benzimidazolium chloride at 1.4 Å and 110.8° at the B3LYP/6–311++G(d,p) level of theory.34 Furthermore, the C13–N11–C16 bond angles for 1 and 2 are calculated as 125.3° and 124.8°, which are ~1° different from the observed angle of 126.0°.
Similarly, the C13–N12–C25 bond angle observed at 126.4° for 1 was predicted at 125.9° for 1 and at 124.8° for 2. However, the N11–C16–C17 and N12–C25–C26 bond angles observed at 112.0° and 111.5° were estimated as 113.4° and 112.6° for 1 and as 114.00° for 2. Here, it should be noted that the calculated angles deviate from the experimental angle value by approximately 2°. However, the C13–N12–C25-C26 dihedral angles for the two compounds were calculated as −98.7° and −110.7°, being almost planar to the ring of each compound, while this angle deviated by 3.7° and 15.7° from the observed value of −95° for 1. Furthermore, the dihedral angles of N11–C16–C17-C20 for compounds 1 and 2 were calculated as −55.7° and −54.4° with a deviation of −13.9° and −15.2° from the recorded value for 1 of −69.6°. Recently, the C3–C2–C1–N11 and C6–C1–C2–C3 dihedral angles of a structurally related compound were calculated as 179.0°.34 In this study, the same dihedral angle was recorded as 177.4° and calculated as 177.9° for compound 1 and as 179.1° for compound 2, respectively. Here, it is worth mentioning that the simulated data are in good agreement with the counterparts observed in the experiment, even though there are some differences between them.
Vibrational analysis
The infrared spectra of the NHC 1 and Ag(I)–NHC 2 are presented in Figures S3(a) and (b) in supplemental material; a comparison of the experimental and the scaled theoretical vibrational frequencies can be found in Table S1 (see the supplemental material).
In the literature, the stretching vibrations of the C–H bonds have been reported in the region of 3000–3100 cm−1.35–37 In this study, the vibrations of the C–H aromatic stretching (νCH RA) modes of 1 and 2 occurred at 3113 and 3104 cm−1, whereas the counterparts of these modes for 1 and 2 were simulated in the range of 3113–3088 cm−1 and at 3112 cm−1, respectively. Besides, the νasC–H (νasC24–H2 and νasC33–H2) vibrational modes of the carbene parts of 1 and 2 were assigned as 3116 and 3118 cm−1, and 3112 cm−1 as a pure mode, whereas the symmetric stretching modes of these groups (νC24–H2 and νC33–H2) were predicted as 3037–3028 cm−1(1) and as 3033–3032 cm−1 (2). The symmetric stretching modes of the methyl groups (νC20–H3 and νC29–H3) for the compounds were recorded as 2937 cm−1 (1) and as 2914 cm−1 (2), and assigned in the range of 2926 cm−1 (85%)–2916 cm−1 (87%) for 1 and as 2927 cm−1 (87%) for 2. From Table S1, the appearing peaks for 1 at 3052 and 3032 cm−1 were assigned as the νC24–H2 and νC33–H2 modes calculated at 3037 cm−1 (98%) and 3028 cm−1 (99%).
The aromatic ring C–C stretching (νCC RA) mode for 1 was assigned as 1682 cm−1 (63%) and 1053 cm−1 (72%) as a pure mode, while this mode has a contribution to the bending modes in the spectral region of 1671–427 cm−1. For example, the observed peaks for 1 at 1552, 1186, and 800 cm−1were assigned by potential energy distribution (PED) analysis at 1549, 1176, and 801 cm−1 as combined modes with ipb (in-plane bending) HCC RA. However, the νCC RA mode for 2 was assigned at 1676 cm−1 (59%) and 1057 cm−1 (56%) as a pure mode, calculated at 1672, 1541, 1534, and 1168 cm−1 as a mixed with ipb HCC RA and at 1400 and 818 cm−1 as contaminated with the νNC mode. It is challenging to predict the NC stretching modes because these modes are generally coupled with C–C stretching modes and with the other stretching modes of the unsaturated ring. In our previous study, the νNC stretching mode of a benzimidazolium ligand was observed at 1553–1464 cm−1 and was predicted to be at 1586–1458 cm−1.38 Here, the observed peaks at 1616 cm−1 with a strong IR intensity for 1 was assigned at 1619 cm−1 (58%). The νNC modes for 2 were observed at 1426, 1396, 1377, 1212, 1190, and 800 cm−1 and assigned at 1414, 1400, 1379, 1211, 1186, and 813 cm−1, in addition to the assigned modes at 1534, 1446, and 1441 cm−1. Here, the electropositive metal center caused the negative shifts because it attracts electron density.39 In Table S1, it can be seen that the assigned modes at 1446 and 1186 cm−1 for 2 are predicted as a pure νNC mode, and the remaining assigned modes are contaminated with the other vibrational modes.
In this study, the methyl group symmetric bending (umbrella, νCH) for compounds 1 and 2 was observed at 1441 and 1439 cm−1, respectively. Moreover, the corresponding mode for two compounds was assigned at 1443 cm−1 (72%)–1435 cm−1 (80%) for 1, and at 1439 cm−1 (60%)–1438 cm−1 (74%) for 2. In the literature, the νCH (umbrella) mode has been observed at 1424 cm−1,38 and it is predicted at the B3LYP level as 1368 and 1414 cm−1.38 The scissoring modes (for σC16–H2 and σC25–H2) of the methylene group for compounds 1 and 2 were computed in the ranges of 1527–1477 and 1444–1349 cm−1, respectively, and were observed at 1489 and 1478 cm−1 for 1 and at 1478 cm−1 for 2. Moreover, the twisting modes for the same bonds were determined at 1549–1263 cm−1 for 1 and at 1441–1231 cm−1 for 2, the wagging modes in the range of 1415–1368 cm−1 for 1 and 1441–1231 cm−1 for 2, and the rocking modes were in the range of 984–969 cm−1 for 1 and 1268–700 cm−1 for 2. Thus, not only methyl and methylene groups but also the other vibrational modes were in good agreement with both the experimental values and previous reports.
NMR spectral analysis
The NMR shifts for compounds 1 and 2 were obtained using the gauge-independent atomic orbital (GIAO)40,41 approaches to (δ) tetramethylsilane (TMS) [δiso = σTMS – σiso]. The observed NMR spectra of both compounds are given in Figures S4(a)–(d). The chemical shifts were simulated in CHCl3 and recorded in CDCl3 and can be found in Table S2. The correlation equations obtained from regression analysis are given in the supplementary material (Table S2). Accordingly, it can be said that the observed and simulated chemical shifts of both compounds are very comparable with each other; the regression coefficients with regard to the 13C isotropic shifts for compounds 1 and 2 were calculated as R2 = 0.9942 and R2 = 0.9939, respectively, whereas 1H isotropic shifts for these compounds were predicted as R2 = 0.8786 and R2 = 0.9835, respectively, in CHCl3 (Table S2(b)).
Accordingly, the characteristic sharp singlet for the acidic NCHN proton of 1 was observed at 11.75 ppm in the 1H NMR spectrum. The absence of this peak was also proof of the formation of Ag(I)–NHC 2. The 13C NMR spectrum of NHC ligand 1 exhibited the characteristic signal of the NCHN carbon at 137.4 ppm. As reported in the literature,42–44 this peak was not observed for the Ag–carbene because of the fluxional behavior of the NHC complexes.45–47 Furthermore, the chemical shifts of the acidic NCHN proton of the benzimidazolium ligand derivative has been reported at 11.0448 and 11.25 ppm.49 In a recent study, it was reported that the carbon atom chemical shifts for benzimidazole complexes were reported in the range of 142.6–149.5 ppm.48 As expected, the unsaturated ring 13C shifts (atom nos. 1–6 and13) for compounds 1 and 2 occurred in the ranges of 116.1–137.4 and 116.1–137.5 ppm, whereas they were simulated in the ranges of 118.1–140.4 and 118.4–142.5 ppm in CHCl3. Furthermore, the aromatic 1H shifts for compounds 1 and 2 occurred in the NMR spectra between 7.64–11.75 and 7.39–7.71 ppm, respectively, whereas these shifts were simulated at 7.95–9.40 and 7.64–7.81 ppm (CHCl3), respectively. Here, it is worth mentioning that the methylene (C16 and C25) and methyl groups (C20 and C29) were responsible for the 13C NMR spectral peaks at 53.5 and 19.7 ppm for 1 and at 55.5 and 19.7 ppm for 2, respectively. Also, the methylene group protons of these compounds produced the peaks at 5.32 ppm (H27 and H28) for 1 and at 4.99–5.31 ppm (H18 and H19) for 2.
Natural bond orbital (NBO) analysis
The natural bond orbital analysis and second-order perturbative theory results defined by Weinhold et al.50,51 are increasingly used to elucidate the possible intramolecular interactions for a specific organic or inorganic molecular system via the prediction of the electronic parameters such as the stabilization energy (E(2)), the donor/acceptor orbital occupancy (qi/qj), the donor and acceptor orbital energies (εi and εj), and off-diagonal Fock matrix element (Fij). In this study, NBO analysis was applied to elucidate and compare the significant intramolecular interactions of the studied ligand and complex molecules (Table S3).
From Table S3, the possible numbers of intramolecular interactions for 1 are higher than for 2. Accordingly, the most robust resonance interaction between the donor–acceptor orbital was determined as LP(1) N11 → π* N12–C13 for 1 (E(2) = 88.76 kcal/mol) and π C2–N12 → π* N11–C13 for 2 (E(2) = 37.84 kcal/mol), respectively, which means that a significant electron density will be present on the antibonding π* N12–C13 orbital for compound 1 in comparison with the same bond for compound 2. In other words, the electron delocalization for compound 1 is higher than that of compound 2, which causes the N12–C13 bond length for compound 1 to be smaller in comparison to the same bond for compound 2. Here, it is worth remembering that the bond lengths of compounds 1 and 2 were calculated as 1.332 and 1.357 Å, respectively. Moreover, the E(2) value of the resonance interaction occurring in the benzene ring of 1 was calculated in the range of 15.40 and 22.54 kcal/mol; the lowest energy interaction is predicted as π C1–C2 → π* N12–C13 and the highest energy interaction is estimated as π C5–C6 → π* C1–C2. However, the π C1–C2 → π* C3–C4 interaction contributing to the stabilization energy for 2 does not have as much energy as the other interactions occurring in this compound. From Table S3, the highest contribution to the stabilization energy for the benzene ring of 2 has been estimated as then π C5–C6 → π* C3–C4 resonance interaction with an energy of 19.27 kcal/mol and a remarkable orbital occupancy. It can also be seen from Table S3 that there is an anomeric interaction (n → σ*) for 1 of 12.02 kcal/mol (EDj = 0.03647e) for the intramolecular charge transfer from a quite polarizable donor orbital LP(4) Cl to the antibonding orbital σ* O38–H40. Moreover, the polarity of this orbital has been reduced by the presence of the Ag atom of 2, and the hybridization of the LP(4) Cl orbital has been calculated to be sp3.52. An important interaction occurring in 2 is that of LP(5) Ag → π* N11–C13 with a stabilization energy of 5.06 kcal/mol (EDi = 1.97270).
Frontier molecular orbital analysis
FMO investigations have been widely used to provide information on the chemical stability and reactivity tendency of molecular systems via evaluation of the quantum chemical tensors.52–55Table 2 shows the results of the FMO analysis for compounds 1 and 2.
The calculated quantum chemical parameters for 1 and 2.
Gas
CHCl3
1
2
1
2
HOMO (-I)
−0.191
−0.228
−0.234
−0.255
LUMO (-A)
−0.086
−0.067
−0.076
−0.056
ΔE (energy gap)
2.837
4.400
4.293
5.422
Χ
−3.767
−4.013
−4.216
−4.224
Η
1.419
2.200
2.147
2.711
Ω
5.002
3.660
4.139
3.291
ΔNmax
2.656
1.824
1.964
1.558
HOMO: highest occupied molecular orbital; LUMO: lowest unoccupied molecular orbital; I: ionization energy; A: electron affinity; Χ: electronic chemical potential; H: global hardness; Ω: electrophilicity index; ΔNmax: charge transfer capability.
The HOMO and LUMO energies are in au; the other parameters are given in eV.
Accordingly, the ionization energy of the compounds changed as 2 (0.228 au) > 1 (0.191 au) in the gas phase and 2 (0.255 au) > 1 (0.234 au) in CHCl3. The electron affinity changed as 1 (0.086 au) > 2 (0.067 au) in the gas and 1 (0.076 au) > 2 (0.056 au) in CHCl3. Besides, the energy gap value of each compound rose with an increase in the polarity of the stimulation media. As far as the energy gap values in both the vacuum and the polar environment are concerned, it can be observed that intermolecular interactions for 2 are more likely than in 1 because the ΔEgap for the compounds is calculated as 2 (5.422 eV) > 1 (4.293 eV) in CHCl3 and as 2 (4.400 eV) > 1 (2.837 eV) in the gas phase. However, it can be said that compound 2 is harder than compound 1 for both phases, as expected. As seen from the MEP diagrams and the net charge analysis, compound 2 is more acidic than compound 1 because of the electropositive Ag atom. It is well known that the electron acceptors are Lewis acids and that electron donors are the Lewis bases.56 In this study, the electrophilicity index values of the compounds in the CHCl3 (condensed) phase also support the relative acidic characterization. In Table 2, compound 1 (5.002 eV) has a more electrophilic character than compound 1 (3.660 eV) in the gas phase as well as in CHCl3 phase. Furthermore, compound 1 (2.656 eV) has more capability of charge transfer than compound 2 (1.824 eV) in the gas phase. It is well known57–60 that the ΔEgap has been commonly used to provide information on the kinetic stability and reactivity of related molecular systems. Thus, the nucleophilicity of compound 2 is greater than compound 1, which can be the reason for the anticancer activity of compound 2.
The MEP plots also imply the electrophilic and nucleophilic sites.61 The positive potential is specified with a blue color which demonstrates the electrophilic attack center, whereas the negative potential is visualized by red, indicating the nucleophilic attack center. In Figure 2, the MEP plot for 1 shows that the negative potential is mostly due to the chloride ion and the water molecule and that the positive potential regions are mainly over the two N atoms belonging to the heterocyclic part of the compound. Similarly, the red color for 2 seems to be mostly over the Cl atom, but it should also be noted that the medium sized-electrostatic potential with the orange color, implying the electrophilic center, is concentrated on the aromatic part of 2 and the blue color, indicating the nucleophilic center, is over the remaining part of compound 2. Also, Figure 2 shows the essential net charges obtained from the natural bond orbital calculations; the full atomic charges for both compounds are given in Table S4. Here, the natural atomic charge calculations revealed that the positive charge was located on the electropositive Ag atom as expected; the net atomic charges for the Cl atoms of compound 2 were calculated as +0.517e and −0.761e. However, the net charge for the Cl atom of compound 1 was predicted as −0.910e. The net charges of the N atoms in compounds 1 and 2 were calculated as −0.367e and −0.426e, respectively. Furthermore, the charge distribution over the imidazoline moiety of compound 2 changed due to the presence of the electropositive Ag atom, and the N atom charges for compound 2 are more negative than those of compound 1.
HOMO and LUMO (isoval: 0.02), MEP (isoval: 0.0004) plots, and the NBO charges for selected atoms of compounds 1 and 2 in CHCl3.
Hirshfeld surface analysis
Hirshfeld surfaces provide a useful perspective to represent the intermolecular interactions in the crystal structure of a specific molecular system. In this context, the two-dimensional (2D) histograms, known as fingerprint plots, are used to explain the types of intermolecular interactions and show the relative region of the related interactions.62,63
Figure 3 shows the Hirshfeld surfaces for 1 generated using Crystal Explorer (version 17.5), which uses the crystallographic information file (CIF) file as the input.64 Then, Hirshfeld surfaces are visualized over the dnorm ranges from −0.4522 to 1.3003 Å to give details of the interactions. Accordingly, the most effective interactions, which occurred between the oxygen (O) and chloride (Cl) atoms, are marked as bright red areas. Furthermore, Figure 4 (the 2D fingerprint plots) illustrate the possible contributions of the intermolecular interactions to the Hirshfeld surfaces. Long spikes are characteristic of hydrogen bonds, which represent the H···Cl/Cl···H type interactions. The left spike (near the de axis) represents the interaction of the molecule as a donor, while the other represents as an acceptor. In addition to this interaction, there are two pairs of little wings of the H···C/C···H type, which represent the C–H···pi type interactions. These C–H···pi type interactions are not strong enough to warrant extensive discussion herein.
dnorm mapped on the Hirshfeld surface for visualizing the intermolecular contacts of compound 1.
The 2D fingerprint plots of the title compound 1 showing all interactions with the percentage contribution to the total Hirshfeld surface area. The parameter di is the closest internal distance from a given point on the Hirshfeld surface, while de is the closest external contact.
In vitro anticancer activities
In this study, three different cell lines were used to determine the anticancer activities of ligand 1 and complex 2. DU-145 cells are androgen insensitive with high metastatic potential and do not express prostate-specific antigen (PSA). MCF-7 breast cancer cells are ER(+). L-929 non-cancer cells were used to determine the toxicities of ligand 1 and complex 2 in normal cells. The evaluation of the effect of ligand 1 and complex 2 on cancer cell viability was applied using the MTT assay at 24, 48, and 72 h.32Figure 5(a)–(f) shows the dose and time-dependent anticancer activities of ligand 1 and complex 2 toward the cancer cells and non-tumorigenic cells. The IC50 (concentration of the test compound to achieve 50% of cell death) values for the compounds are listed in Table 3 for all the cell lines.
The dose and time-dependence of the in vitro anticancer activities of ligand 1 and complex 2 against DU-145 and MCF-7 cancer cells and L-929 non-cancer. Control cells were treated with DMSO. Data are representative of the mean of three separate experiments performed in triplicate and are reported as SEM. (a), (c), (e) in vitro anticancer activities of ligand I against DU-145, MCF-7 cancer cells and L-929 non-cancer cells, respectively. (b), (d), (f) in vitro anticancer activities of complex 2 against DU-145, MCF-7 cancer cells and L-929 non-cancer cells, respectively.
In vitro anticancer activities of ligand 1 and complex 2 against DU-145, MCF-7, and L-929 cells.
MTT: 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide; SEM: standard error mean.
Each IC50 value represents the mean ± SEM of three independent experiments (nine replicates).
Cell viability after treatment for 24, 48, and 72 h was determined by MTT staining as described in section “Experimental” (µM).
Cancer cells.
Non-cancer cells.
Although the IC50 values of ligand 1 against all the tested cell lines were >20 µM at all the time points, complex 2 showed different IC50 values depending on the time and cell line type. Complex 2 did not have IC50 values against DU-145 and L-929 cells at 24-h exposure, at least at doses equal to 20 µM. These results showed that complex 2 had lower IC50 values compared to ligand 1. As lower IC50 values indicate higher activity, complex 2 showed higher anticancer activity than ligand 1. The lowest IC50 values were determined against MCF-7 breast cancer cells with complex 2 (13.6, 4.83, and 2.18 µM after 24, 48, and 72 h exposure, respectively), compared to the DU-145 and L-929 cells. Therefore, MCF-7 breast cancer cells were most susceptible to complex 2. The selectivity index (SI) values of complex 2 were calculated as >1.47, 2.6, and 2.97 against MCF-7 breast cancer cells and >1.2, and 1.15 against DU-145 prostate cancer cells at 24, 48, and 72 h, respectively. The SI values indicated that complex 2 was more cytotoxic and selective against MCF-7 cells than DU-145 cells. More importantly, complex 2 had higher IC50 values for L-929 normal cells and was more cytotoxic toward DU-145 and MCF-7 cancer cells. Therefore, complex 2 demonstrated selectivity between both prostate and breast cancer cells as well as between healthy and cancer cells. Although ligand 1 had no anticancer activity in any cell lines, complex 2 showed dose and time-dependent anticancer activity against all cell lines. More importantly, complex 2 displayed lower anticancer activities on L-929 non-cancer cells than cancer cells.
Conclusion
This paper has focused on the synthesis of benzimidazolium salt and its Ag(I)–NHC complex, and their structural determination, theoretical studies, the comparison of in vitro anticancer activities of both compounds, and the Hirshfeld surface analysis of the benzimidazolium salt. Both compounds were characterized by spectroscopic techniques in addition to characterizing compound 1 by single-crystal X-ray diffraction. NBO analysis revealed that the main contribution to the stabilization energy lowering was due to the resonance interactions for both compounds; the greatest contribution to the lowering of the stabilization energy was calculated as LP(1) N11 → π* N12–C13 (E(2) = 88.76 kcal/mol) for 1 and π C2–N12 → π* N11-C13 (E(2) = 37.84 kcal/mol) for 2, respectively. Moreover, it can be concluded that compound 2 is more chemically reactive than compound 1 as it has a higher ΔEgap value: the ΔEgap values for the compounds were calculated as 2 (5.422 eV) > 1 (4.293 eV) in CHCl3 and 2 (4.400 eV) > 1 (2.837 eV) in a vacuum. The in vitro anticancer activities of the compounds were also investigated against DU-145 prostate cancer cells, MCF-7 breast cancer cells, and L-929 non-cancer cells using the MTT cell viability assay for 24, 48, and 72 h. The results showed that complex 2 demonstrated time and dose-dependent anticancer activities toward DU-145 and MCF-7 cancer cells. The results also indicated that ligand 1 had lower anticancer activity against cancer cells compared to the Ag(I) complex 2. As the lower IC50 values indicate higher anticancer activity, complex 2 displays the highest anticancer activity against MCF-7 cell lines at all time points studied. More importantly, complex 2 displayed lower anticancer activity toward L-929 non-cancer cells than cancer cells.
Experimental
Materials and measurements
All experiments were performed under argon in flame-dried glassware using standard Schlenk techniques. All reagents were purchased from Sigma-Aldrich Co. (Dorset, UK). The solvents used were purified by distillation over appropriate drying agents and were transferred under argon. Melting points were determined using an Electrothermal 9100 melting point detection apparatus in capillary tubes, and the melting points are uncorrected. FTIR spectra were recorded in the range of 400–4000 cm−1 on a Perkin Elmer Spectrum 100 FTIR. 1H NMR and 13C NMR spectra were obtained using a Bruker As 400 Mercury spectrometer operating at 400 MHz (1H) and 100 MHz (13C) in CDCl3 with TMS as the internal reference. 1H NMR peaks are labeled as singlet (s) and multiplet (m). Chemical shifts and coupling constants are reported in ppm and Hz, respectively. All the measurements were taken at room temperature using freshly prepared solutions.
Synthesis
NHC ligand 1 and its Ag(I)–NHC complex 2 were prepared under an argon gas atmosphere according to the literature.45
Benzimidazole (10 mmol) was added to a solution of NaH (10 mmol) in dry tetrahydrofuran (THF) (30 mL) and the mixture was stirred for 1 h at room temperature. 2-Methyl-2-propenylchloride (10.1 mmol) was added dropwise and the obtained solution was heated for 24 h at 60 °C. Next, the solvent was removed under vacuum. Dichloromethane (50 mL) was added to the solid. The mixture was filtered, and the obtained clear solution was concentrated under vacuum. The remaining solution was distilled to give 1-(2-methyl-2-propenyl)benzimidazole. 1-(2-Methyl-2-propenyl)benzimidazole (1 mmol) and 2-methyl-2-propenylchloride (1 mmol) were stirred in dimethylformamide (DMF) (5 mL) for 24 h at 80 °C, and the product precipitated. The solution was filtered and the solid was rinsed with diethyl ether and dried under vacuum. The crude product was recrystallized from dichloromethane/diethyl ether. Yield: 84%, m.p. 168–169 °C. IR: ν(CN): 1552 cm−1. 1H NMR (400 MHz, CDCl3): δ = 1.90 (s, 6H, NCH2C(CH3)CH2), 4.99 (s, 2H, NCH2C(CH3)CH2), 5.14 (s, 2H, NCH2C(CH3)CH2), 5.32 (s, 4H, NCH2C(CH3)CH2), 7.64–7.65 and 7.75–7.76 (m, 4H, NC6H4N), 11.75 (s, 1H, NCHN).13C NMR (100 MHz, CDCl3): δ = 19.7 (NCH2C(CH3)CH2), 53.5 (NCH2C(CH3)CH2), 113.8 (NCH2C(CH3)CH2), 137.4 (NCHN), 116.1, 127.3, and 131.5 (ArC) and 144.3 (NCH2C(CH3)CH2).
A solution of Ag2O (0.5 mmol) and of NHC (1 mmol) ligand 1 in CH2Cl2 (25 mL) was stirred at room temperature for 24 h in the dark. Next, the mixture was filtered through Celite. The clear filtrate was evaporated under vacuum to afford the crude product, which was then recrystallized from dichloromethane/diethyl ether. Yield: 79%, m.p. 151–153 °C, IR: ν(CN): 1378 cm−1. 1H NMR (400 MHz, CDCl3): δ = 1.73 (s, 3H, NCH2C(CH3)CH2), 1.80 (s, 3H, NCH2C(CH3)CH2), 4.83 (s, 1H, NCH2C(CH3)CH2), 4.98 (s, 1H, NCH2C(CH3)CH2), 4.99 (s, 2H, NCH2C(CH3)CH2), 5.03 (s, 1H, NCH2C(CH3)CH2), 5.14 (s, 1H, NCH2C(CH3)CH2), 5.31 (s, 2H, NCH2C(CH3)CH2), 7.37–7.40, 7.45–7.48, 7.61–7.64, and 7.69–7.73 (m, 4H, NC6H4N). 13C NMR (100 MHz, CDCl3): δ = 19.7 and 20.0 (NCH2C(CH3)CH2), 53.6 and 55.5 (NCH2C(CH3)CH2), 112.1 (NCH2C(CH3)CH2), 113.8 (NCH2C(CH3)CH2), 139.1 (NCH2C(CH3)CH2), 144.3 (NCH2C(CH3)CH2), 114.5, 116.1, 124.3, 127.2, 133.9, and 137.5 (ArC), Ag–C (carbene) not observed.
X-ray crystallography and refinement of 1
Single-crystal X-ray diffraction data of the NHC ligand were collected at room temperature on a Rigaku–Oxford Xcalibur diffractometer with an electro-optical system (EOS)-charge-coupled device (CCD) detector using graphite-monochromated MoKα radiation (λ = 0.71073 Å) with CrysAlisPro software.65 Data reduction and analytical absorption correction were performed using the CrysAlisPro program.66Utilizing OLEX2,67 the structure was solved using the intrinsic phasing method with SHELXT and refined by full-matrix least squares on F2 in SHELXL.68,69 Anisotropic thermal parameters were applied to all non-hydrogen atoms. Hydrogen atoms were placed using standard geometric models and with their thermal parameters riding on those of their parent atoms (C–H = 0.93–0.96–0.97 Å). The details of the crystal data and structure refinement of the title compound are given in Table 4.
Crystallographic data and structure refinement parameters for 1.
Crystallographic data as .cif file for the structure reported in this paper have been deposited at the Cambridge Crystallographic Data Center with CCDC 1822823.
Figure 6 displays the overlay of the X-ray and calculated molecules in the asymmetric unit, which reveals that their structures are virtually equivalent. The root-mean-square deviation (RMSD) value was found to be 0.356 Å with inversion and 0.384 Å without inversion. The crystal structure of the title molecule is consolidated by hydrogen bonds involving the chloride anion, propenyl moieties of the cation, and the water molecule, which are linked through an intricate H-bonding network consisting of two Ow–Hw∙∙∙Cl and three C–H∙∙∙Cl interactions (Table 5). Accordingly, the chloride anion behaves as an H-bond acceptor in the crystal structure, resulting in the formation of the one-dimensional supramolecular array (see Figure S1). Besides, the Cl− ions were hydrated by forming an anionic hydrogen-bonded cluster [Cl2(H2O)2]2− via the intermolecular O1w–H1wA∙∙∙Cl1i and O1w–H1wB∙∙∙Cl1ii hydrogen bonds. These chloride–water tetrameric clusters have a rectangular-like geometry and generate a cyclic graph-set motif in Etter’s graph notation.70 They were also held by two cations through the C8–H8B∙∙∙Cl1ii and C12–H12A∙∙∙Cl1i interactions to form a -type hydrogen-bonded motif, as shown in Figure S2. There is another tetrameric cluster formed by the C12–H12A∙∙∙Cl1i and C12–H12B∙∙∙Cl1iii hydrogen bonds, which reveals the graph-set notation. All these bifurcated hydrogen bonds are responsible for the stabilization and packing in the supramolecular architecture of the crystal structure.
Atom-by-atom superimposition of the calculated structure (red) on the X-ray structure (green) for 1.
Hydrogen-bonding interactions (Å, °) for compound 1.
D–H∙∙∙A
D–H
H∙∙∙A
D∙∙∙A
D–H∙∙∙A
O1w–H1wA∙∙∙Cl1i
0.85
2.35
3.182(3)
168
O1w–H1wB∙∙∙Cl1ii
0.85
2.58
3.187(3)
129
C8–H8B∙∙∙Cl1ii
0.97
2.75
3.668(3)
159
C12–H12A∙∙∙Cl1i
0.97
2.78
3.687(3)
156
C12–H12B∙∙∙Cl1iii
0.97
2.72
3.681(3)
173
Symmetry codes: (i) −1 + x, y, z; (ii) 1 − x, 2 − y, 1 − z; and (iii) 1 − x, 1 − y, 1 − z.
Density-functional theory studies
The molecular geometries of compounds 1 and 2 were optimized using B3LYP71,72 functional and the 6–311+G(d,p)/LANL2DZ basis set using Gaussian 09W software73 and verified by the non-negative in frequency. The optimized structures of the compound were used for further calculations in both a vacuum and in CHCl3. The vibrational modes of both compounds have been scaled74 with factor of 0.9688 for high frequencies and 1.0189 for low frequencies and assigned by PED analysis using the vibrational energy distribution analysis (VEDA)75 program. NMR chemical shifts of both compounds were obtained using the GIAO36,37 approach by subtracting the shielding constants of TMS. NBO50,51 and FMO52–55 analyses were conducted to investigate the intramolecular interactions and biological activity tendency of both compounds. FMO amplitudes and MEP plots provided information on the possible reactive regions of the compounds, and this was visualized using GaussView 6.0.16.76 All quantum chemical calculations in the polar environment were performed using the polarized continuum model (PCM).77,78
Cell cultures
The human cancer cell lines DU-145 (HTB-81, human prostate carcinoma), and MCF-7 (HTB-22, human breast adenocarcinoma), were purchased from American Type Culture Collection (ATCC) (Manassas, VA, USA). L-929 (non-cancer cells adipose from mouse) were purchased from European Collection of Animal Cell Culture (ECACC, Salisbury, UK). Eagle’s minimum essential medium (EMEM, 30-2003), fetal bovine serum (FBS, 30-2020), and penicillin and streptomycin (30-2300) were purchased from ATCC. Dulbecco’s modified Eagle’s medium (DMEM, D6429) and trypsin-ethylenediaminetetraacetic acid (EDTA) solution (T-3924) were purchased from Sigma-Aldrich (Sigma-Aldrich Chemie GmbH, Steinheim, Germany).
MTT assay
The MTT method was used to determine the anticancer activity of ligand 1 and complex 2.32 The DU-145 cell line was cultured in EMEM; MCF-7 and L-929 cells were cultured in the DMEM supplemented with 10% FBS and 1% penicillin/streptomycin solution. All cells were cultured in an incubator with 5% CO2 at 37 °C. Cells were passaged when the confluence of the cells reached 80% or higher. Cells were seeded in a 96-well plate at a density of 1 × 105 cells/well and allowed to adhere for 24 h at 37 °C in a CO2 incubator. 1 μL of different concentrations (1–20 µM) of the compounds were added to each well, and the cells were treated for 24, 48, and 72 h. Control and negative control wells were treated with culture medium and sterile DMSO. At the end of the time points, 10 μL of the MTT solution was added to each well and allowed to incubate at 37 °C for a further 2 h. After complete removal of the media, DMSO was added to each well to dissolve the formazan; then, the dye plates were incubated for 15 min at room temperature. Optical density was measured at 570 nm in enzyme-linked immunosorbent assay (ELISA) (Biotek, Epoch, USA). Data represent the average values of three independent measurements with standard error means (±SEM).
Statistical analysis
All experiments were carried out in triplicate, and the results are expressed as means ± SEM. Data were analyzed using one-way analysis of variance and differences were considered statistically significant at *p < 0.05, **p < 0.005, ***p < 0.0005, and ****p < 0.0001. The IC50 values were determined using statistical software, GraphPad Prism7 (GraphPad Software, San Diego, CA, USA).
Supplemental Material
sj-pdf-1-chl-10.1177_1747519820950219 – Supplemental material for Synthesis, in vitro anticancer activities, and quantum chemical investigations on 1,3-bis-(2-methyl-2-propenyl)benzimidazolium chloride and its Ag(I) complex
Supplemental material, sj-pdf-1-chl-10.1177_1747519820950219 for Synthesis, in vitro anticancer activities, and quantum chemical investigations on 1,3-bis-(2-methyl-2-propenyl)benzimidazolium chloride and its Ag(I) complex by Goncagül Serdaroğlu, Serap Şahin-Bölükbaşı, Duygu Barut-Celepci, Resul Sevinçek, Neslihan Şahin, Nevin Gürbüz and İsmail Özdemir in Journal of Chemical Research
Footnotes
Acknowledgements
Dokuz Eylül University is acknowledged for the use of the Rigaku–Oxford Xcalibur EoS Diffractometer (purchased under University Research Grant No. 2010.KB.FEN.13). Also, the authors thank Sivas Cumhuriyet University Scientific Research Projects Department (Project No. CUBAP: EĞT-086) and İnönü University Scientific Research Projects Department (Project No. FOA-2018-1135) for their financial support. All calculations were carried out at TUBITAK ULAKBIM, High Performance and Grid Computing Center (TR-Grid e-Infrastructure). In vitro anticancer activity studies were performed at Advanced Technology Research and Application Center (CUTAM), Sivas Cumhuriyet University.
Author contributions
All authors contributed to the study conception and design. Synthesis, 1H NMR, 13C NMR, and FTIR spectra were performed by Neslihan Şahin, Nevin Gürbüz, and İsmail Özdemir. Crystallographic studies were performed by Duygu Barut-Celepci and Resul Sevinçek. Density-functional theory studies were performed by Goncagül Serdaroğlu. Cell culture studies were performed by Serap Şahin-Bölükbaşı. The first draft of the paper was written by Serap Şahin-Bölükbaşı, Goncagül Serdaroğlu, Neslihan Şahin, and Duygu Barut-Celepci. All authors commented on the paper and read and approved the final paper.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this paper.
Ethical approval
All ethical guidelines have adhered.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this paper.
ORCID iD
Goncagül Serdaroğlu
Supplemental material
Supplemental material for this paper is available online.
References
1.
WanzlickHWSchönherrHJ.Angew Chem Int Ed1968; 7: 141.
2.
ÖfeleK.J Organomet Chem1968; 12: 42.
3.
ArduengoAJIIIHarlowRLKlineMJ.J Am Chem Soc1991; 113: 361.
ForesmanJBKeithTAWibergKB, et al. J Phys Chem1996; 100: 16098.
78.
TomasiJMennuciBCammiR.Chem Rev2005; 105: 2999.
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.