Anticancer Activities and QSAR Study of Novel Agents with a Chemical Profile of Benzimidazolyl-Retrochalcone ()
1. Introduction
Cancer is characterized by a wide range of conditions involving abnormally large and unregulated cell proliferation in normal organism tissue. This proliferation can affect any part of the organism [1]. Cancer is currently a major cause of death worldwide. Indeed, it is the second leading cause of death, with 8.8 million deaths in 2015 according to the WHO [1] [2] [3]. Moreover, the growing economic impact of cancer is considerable. In 2010, the total annual cost of the disease was estimated to be approximately $1.16 trillion (US $) [4]. Among the main types of cancer that are fatal are lung, liver and colorectal cancer [1]. Chemotherapy treatment of cancer remains one of the preferred strategies for cancer management. Unfortunately, current cancer drugs are either ineffective for some people or too expensive for many populations.
It is in this perspective of contributing to the accessibility of new high-performance anti-cancer drugs that we have initiated this pharmacochemical and particularly a molecular modelling research [5] [6]. To this end, we have focused on 1,3-diphenylpropenones (chalcones) and the benzimidazole nucleus as pharmacophore and heterocyclic carrier, respectively. Indeed, taking into account their multiple intrinsic biological properties, it was normal to imagine that the combination of these chemical entities according to the concept of hybrid molecule, could generate new molecules likely to present new pharmacological activities, especially interesting anticancer ones.
Indeed, chalcones are an important class of compounds with antimalarial [7], anti-leishmaniac [8], antifungal [9] [10], antibacterial [11] [12] and anticancer [13] [14] [15] activities. Similarly, the benzimidazole-type diaza-heteroaryl is the pharmacophoric vector of numerous compounds with a broad spectrum of pharmacological activity, including anticancer [16], antiviral [17] [18], antibacterial [19] [20], anthelmintic [21] [22] and antifungal [23] activities.
Also, our benzimidazolyl-chalcones, which are actually retrochalcones, have been obtained in one of our previous work [24], by combining 1H-benzimidazo-2-yl-carbaldehyde derivatives with various acetophenones and analogues. These novel retrochalcones were evaluated in vitro for their antitumor activity only on the NCI-H727 lung cancer cell line. In addition, their cytotoxicity was established in normal human skin fibroblast cells [24]. Taking into account our previous results on benzimidazolyl-chalcone series, in particular their anthelmintic and antifungal activities [25] - [30], as well as the first promising anticancer results in this retrochalcone series [24].
We propose in the present work, to extend their screening to six (06) other human cancer cell lines. The other objective of this work was to establish, by a QSAR study, the quantum descriptors according to the DFT calculation method and the B3LYP/6-31+G (d, p) theory and from the cancer line whose proliferation was most inhibited by our compounds. The aim of such a study was to lay the foundations for the development of new anticancer agents with a chemical profile of benzimidazolyl-retrochalcone.
2. Experimental Section
2.1. Chemistry
We have recently published the synthesis and physico-chemical characteristics of the new series of benzimidazolyl-retrochalcone derivatives and their anticancer activity against human lung cancer cell line NCI-H727 [24]. The chemicalstructures of 13 of the benzimidazolyl-retrochalcone derivatives are reported in Table 1 [24].
Table 1. Chemical structures of benzimidazolyl-retrochalcone derivatives.
2.2. Biology
2.2.1. Cell Culture
Skin normal fibroblastic cells were purchased from Lonza (Basel, Switzerland), HuH7, Caco-2, MDA-MB-231, MCF-7, HCT-116, PC3 and NCI-H727 cancer cell lines were obtained from the ECACC (European Collection Authenticated Cell Cultures) collection (Porton, UK).
Cells were grown according to ECACC recommendations in DMEM (Dulbecco’s Modified Eagle’s Medium) for HuH7, MDA-MB-231, MCF-7 and fibroblast, in EMEM (Eagle’s Minimum Essential Medium) for CaCo-2, in McCoy’s for HCT116 and in RPMI (Roswell Park Memorial Institute medium) for PC3, NCI-H727.
All culture media were added with 10% of FBS, 1% of penicillin-streptomycin and 2 mM glutamine and were incubated at 37˚C in a 5% CO2 atmosphere.
2.2.2. Cytotoxic Assay
Chemicals were solubilized in DMSO at a concentration of 10 mM (stock solution) and diluted in culture medium to the desired final concentrations.
The dose effect cytotoxic assays (IC50 determination) were performed by increasing concentrations of each chemical (final well concentrations: 0.1 μM - 0.3 μM - 0.9 μM - 3 μM - 9 μM – 25 μM). Cells were plated in 96 wells (4000 cells/well).
Twenty-four hours after seeding, cells were exposed to chemicals. After 48 h of treatment, cells were washed in PBS and fixed in cooled 90% ethanol/5% acetic acid for 20 minutes. Then, the nuclei were stained with Hoechst 33,342 (B2261 Sigma).
Image acquisition and analysis were performed using a Cellomics Array Scan VTI/HCS Reader (Thermo Scientific). The survival percentages were calculated as the percentage of cell number after compound treatment over cell number after DMSO treatment.
2.3. QSAR and Theoretical Calculations
2.3.1. Materials and Methods of Calculation
The nine (09) compounds of the training set and the four (04) other validation set used in this study, have their concentrations ranging from 0.83 to 8.86 µM. This range of concentrations makes it possible to define a quantitative relationship between the anticancer activity and the theoretical descriptors. Biological data are generally expressed as the opposite of the log 10 base of activity (−log10(C)) in order to obtain higher mathematical values when the structures are biologically very efficient [31] [32].
2.3.2. Level of Calculation
The calculation to establish the relation between the biological activities values of the studied molecules and their molecular structures were obtained by quantum chemistry calculations using Gaussian software 09 [33]. DFT methods are generally known to generate a variety of molecular properties in QSAR studies [34] - [41]. These methods increase predictability, reduce computation time, and involve costs in the conception of new drugs [42] [43]. The B3LYP/6-31+G (d, p) theory level was used to determine the molecular descriptors. The modelling was done using the multilinear regression method implemented in Excel [44] and XLSTAT [45].
2.3.3. Quantum Descriptors
For the development of the QSAR model, fourteen theoretical descriptors have been calculated such as electronic energy (Eelec), HOMO energy (εHOMO), LUMO energy (εHOMO), gap energy (ΔE), ionisation energy (IE), electronic affinity (EA), chemical hardness (η), chemical softness (S), electrophilicity (ω), chemical potential (µpot), electronegativity (χ), dipole moment (µ), polarizability (α) and lipophilicity (logP). Among these descriptors, the combination of four of them permitted to establish a good model. We have the electronic energy (Eelec), the chemical softness (S), the chemical hardness (η) and the lipophilicity (logP). The electronic energy (Eelectr) represents the electronic contribution of all atoms in a given molecule. The chemical softness (S) indicates the capacity of an atom or a molecule to conserve an acquired charge [46] [47]. The higher the chemical softness, the less stable the molecule. The chemical hardness (η) is a measure of stability of the molecule [48]. Lipophilicity is an important measure for the identification of the similarity with a drug [49]. All these descriptors are determined from optimized molecules follow by frequency calculation. The calculation of the partial correlation coefficient between the descriptors studied is less than 0.70 (aij < 0.70); which means that these different descriptors are independent of each other [50].
2.3.4. Estimation of the Predictive Ability of a QSAR Model
The quality of a model is determined according to certain criteria such as the coefficient of determination R2, the standard deviation S, the correlation coefficients of cross-validation
and Fischer coefficient F. The statistical indicators R2, S and F relate to the adjustment of calculated and experimental values. They describe the predictive capacity within the model limits and allow estimating the precision of values calculated on the training set [51] [52]. The cross-validation coefficient gives information on the predictive power of the model. The coefficients R2 express the dispersion of theoretical values around the experimental ones. The quality of a model is better when the points are close to the adjustment line [53]. The adjustment of the points to the line can be evaluated by the coefficient of determination.
(1)
: Experimental value of anticancer activity
: Theoretical value of anticancer activity
: The mean of the experimental values of anticancer activity
The closer the value of R2 to 1, the more the theoretical and experimental values are correlated. In addition, the variance σ2 is determined by the following relationship (2):
(2)
where k is the number of independent variables (descriptors) is, n the number of molecules in the test or learning set and
the degree of freedom. The standard deviation S is another statistical indicator used to evaluate the reliability and precision of a model:
(3)
The Fischer coefficient F is also used to express the level of statistical significance of the model, that is to say the quality of the choice of the descriptors constituting the model.
(4)
The coefficient of determination of cross-validation
, which allows evaluating the accuracy of the prediction on the test set and is calculated by using the following equation:
(5)
The performance of a model, according to Eriksson et al. [54] is characterized by a value of
for a satisfactory model when for the excellent model
is higher than 0.9. A training set of a model will perform well if the acceptance criterion
is respected. Moreover, the predictive power of a model can be obtained from the value of the log(1/C)theo/log(1/C)exp ratio for the test set. The model is acceptable when the value of the ratio of theoretical and experimental activity tends towards unity.
3. Results and Discussion
3.1. Pharmacochemical Section
All retrochalcone derivatives (compound 1 - 13) were evaluated for their anticancer activities in vitro on a cell panel. This is represented by seven (07) human cancer cell lines gathering in five (05) types cancer and skin normal fibroblastic cells. The products Roscovitine and Paclitaxel were used as positive anticancer controls. The results obtained are reported in Table 2. The results of the anticancer activities expressed in IC50 (µM) showed that most retrochalcones exhibit strong antitumor activity in all cell lines. These compounds also showed cytotoxicity to normal fibroblasts with IC50 values ranging from 0.61 to 9.29 µM.
Table 2. In vitro anticancer activities of benzimidazolyl-retrochalcones against cancer cell lines and normal human fibroblasts.
Except compounds 5, 9, 12 and 13 all other benzimidazolyl-retrochalcones showed better anticancer activity compared to the performance of the Roscovitine. However, Paclitaxel remains by far the most effective drug by any cell line. In particular, the conformational blocking of the functional enchainment propenone with phenolic OH in a coumarin-like ring (compounds 13), did not induce excellent activity.
For the retrochalcones that showed promising anticancer activities, we can note that they carry on their benzene ring (Ar) either a hydroxyl group or a methoxyl group in isomeric position 2 or 3 (compounds 2, 3 and 4). When this benzene ring is unsubstituted (compound 1), the activity obtained is moderate against all tumour cell lines with IC50 values ranging from 1.48 to 7.77 µM. Also, the replacement of benzene ring by furanyl (compound 6) and benzimidazolyl (compound 7) heterocycles type showed better activities than retrochalcone 1. However, compound 5 substituted with nitro group, led to a significant decline in anticancer activities compared to compound 1.
The C-5 substitution of benzimidazole with a chlorine, nitro or benzoyl group does not involve an increase in anticancer activities compared to the non-C-5 benzimidazole substituted (compound 1). The compound 2 with hydroxyl group in position 2 in benzene ring, has been particularly effective against prostate, colon (HCT-116) and breast (MDA-MB 231) cancer cell lines with IC50 values of 1.08, 0.83 and 1.13 µM respectively. Biological activities on colon cancer (HCT-116) have been used to establish the QSAR’s model.
3.2. QSAR and Theoretical Calculations Section
3.2.1. Training and Validation Set
The results of training set of nine (09) benzimidazolyl-retrochalcones and the four (04) molecules of the validation set are reported in Table 3. As for Table 4, it presents the values of partial correlation coefficients aij of the descriptors.
The partial correlation aij between the descriptors is less than 0.70. This shows the independence of the descriptors used in the model.
3.2.2. QSAR Model Validation
The positive and the negative sign of the descriptor’s coefficient of the model reflects the proportionality’s effect between anticancer activity evolution and descriptors of the regression equation. A positive sign shows that the descriptor influencing positively the activity. Then, the negative sign indicates the opposite effect.
Table 3. Training set of nine (09) molecules and the four (04) molecules of the validation set associated with the molecular descriptors and anticancer activity (HTC-116).
*Validation set of the model.
Table 4. Values of partial correlation coefficients of the descriptors.
The best QSAR models obtained for anticancer activity against HCT-116 colon cancer is given by the following relationship: log(1/C) = 79.684 + 3.896E−06 * Eelec (kcal∙mol−1) - 20.469 * η(ev) - 64.245 * S(ev)−1 + 0.301 * logP.
The statistical indicators of multilinear regression are shown in Table 5.
The electronic energy (Eelec) and the lipophilicity (logP) have a positive sign in the model, which suggest that increased activity can be achieved by increasing the electronic energy and the lipophilicity.
The chemical hardness (η) and the chemical softness (S) have a negative sign in the model, which suggest that increased activity can be achieved by decreasing the chemical hardness and softness. A good model is qualified by its high value of the coefficient of Fischer F = 170.803, the cross-validation coefficient
which is 0.954. The values of log(1/C)theo/log(1/C)exp ratio of external validation set are shown in Table 6.
All values of log(1/C)theo/log(1/C)exp tend to 1. This indicates the good correlation between the theoretical and experimental anticancer activity of the molecules studied. The regression line between the experimental and theoretical anticancer activity of training set (blue dots) and the test set (red dots) is illustrated in Figure 1.
3.2.3. Analysis of the Contribution of Descriptors in the Model
The relative contribution of the descriptors in the prediction of anticancer activity of the molecules is shown in Figure 2.
The chemical hardness (η) and the chemical softness (S) have a wide contribution. The lipophilicity presents a weak contribution.
Table 5. Statistical indicators of multilinear regression.
Table 6. Values of the relationship between theoretical and experimental anticancer activity of validation set.
Note: The mathematical model used for the theoretical calculations is only valid for the benzimidazolyl-retrochalcone series.
Figure 1. Regression line of the model.
Figure 2. Contribution of different descriptors.
4. Conclusion
In this work, we evaluated a series of benzimidazolyl-retrochalcone derivatives for their in vitro cytotoxic activity against a panel of seven (7) human cancer cell lines including hepatitis cancer (Huh-7), prostate cancer (PC3), non-small cell lung cancer (NCI-H727), colon cancer (CaCo2, HCT-116), breast cancer (MDA-MB 231, MCF-7) and a normal fibroblast cell of human skin. Benzimidazolyl-retrochalcone showed promising anticancer activities with IC50 values ranging from 0.83 to 14.70 µM and also proved to be less toxic than the reference molecule Paclitaxel and Roscovitine. Among these, compounds 2 and 6 were found to be more effective against colon cancer (HCT-116) with IC50 values of 0.83 and 0.86 µM respectively. QSAR study undertaken to establish a relationship between the anticancer activities on the colon cell line (HCT-116) and the chemical profile of the thirteen (13) benzimidazolyl-retrochalcone derivatives, showed a good correlation with the theoretical descriptors calculated by the DFT method. Thus, electronic energy (Eelec), lipophilicity (logP), chemical softness (S) and chemical hardness (η) play an important role in the induction of antitumor activities of retrochalcones. Negative signs of chemical sweetness (S) and chemical hardness (η) increase anticancer activity. Our QSAR model thus shows that it is possible to predict and especially increase the cytotoxicity of the benzimidazolyl-retrochalcones against the colon cancer cell line (HCT-116).
Acknowledgements
The authors express their deep thanks to Campus France and AMRUGE-CI/C2D for financial support. We are very grateful to Dr Rémy Le Guével, responsible for the platform ‘‘ImPACcell’’ of Rennes 1 university, for anticancer assay in human cell lines.