Abstract
Background/Aim: Twenty-four 2-azolylchromones were subjected to quantitative structure–activity relationship (QSAR) analysis based on their cytotoxicity and tumor specificity, in order to find their new biological activities. Materials and Methods: Cytotoxicity against two human oral squamous cell carcinoma cell lines and two human normal oral mesenchymal cells was determined by the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide method. Tumor specificity (TS) was evaluated by the ratio of the mean 50% cytotoxic concentration (CC50) against oral cells to that against oral squamous cell carcinoma cell lines. Potency-selectivity expression (PSE) value was calculated by dividing the TS value by CC50 against tumor cells. Apoptosis markers were detected by western blot analysis. Physicochemical, structural and quantum-chemical parameters were calculated based on the conformations optimized by force-field minimization. Results: Three sets of 4H-1-benzopyran-4-ones with indole ring showed much higher TS values than those with pyrrole, pyrazole, imidazole, 1,2,4-triazole, 1,2,3-triazole, indazole and benzimidazole rings. Among those with an indole ring, the compound having a 6-methoxy group, that exhibited the highest cytotoxicity, yielded one to three-order higher PSE values to compared with other groups of compounds. Western blot analysis demonstrated that this compound stimulated the cleavage of caspase-3, suggesting the induction of apoptosis. QSAR analysis demonstrated that TS values were correlated with 3D shape, polarizability, ionization potential and lipophilicity. Conclusion: Chemical modification of the lead compound may be a potential choice for designing a new type of anticancer drug.
4H-1-Benzopyran-4-ones (chromones) are an important class of oxygenated heterocyclic compounds, since the core structure is found ubiquitously in the plant kingdom in notable amounts. They show anti-inflammatory, antitumor and antimicrobial activities, and inhibitory activities towards cyclo-oxygenases, kinases, phosphatases, aromatases, acetylcholinesterases, and monoamine oxidases (1, 2).
Among a total of 133 compounds, (E)-3-[2-(4-hydroxyphenyl)ethenyl]-6-methoxy-4H-1-benzopyran-4-one (classified as 3-styrylchromones), (E)-3-[2-(4-chlorophenyl) ethenyl]-7-methoxy-2H-1-benzopyran (classified as 3-styryl-2H-chromenes) were found to have the highest tumor specificity with the least keratinocyte toxicity (3). This finding urged us to synthesize various compounds derived from chromone.
Recently we synthesized a series of 2-azolylchromone derivatives that showed potent and selective inhibition of monoamine oxidases A and B (4). The publication of data on the biological activity of azolylchroman derivatives, however, has been limited. A series of imidazolylchromanone oximes containing phenoxyethyl ether moiety showed marginal activity against filamentous fungi (Aspergillus fumigatus and Exophiala dermatitidis), while they exhibited potent activity against yeast (Cryptococcus gattii), possibly binding to the homology modeled active site of Cryptococcus CYP51-14α-demethylase (5). 7-Chloro-3-(1H-imidazol-1-yl)chroman-4-one more effectively delayed seizures, as well as decreasing seizure duration, in a model of epilepsy (6).
In continuation of discovering new biological activities of chromone derivatives, a total of 24 2-azolylchromone derivatives were investigated for their cytotoxicity against two human oral squamous cell carcinoma (OSCC) and two human normal oral cells, and then subjected to quantitative structure–activity relationship (QSAR) analysis.
Structure of 24 2-azolylchromones of this study.
Materials and Methods
Materials. The following chemicals and reagents were obtained from the indicated companies: Dulbecco's modified Eagle's medium (DMEM), from GIBCO BRL, Grand Island, NY, USA; fetal bovine serum (FBS), 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), doxorubicin from Sigma-Aldrich Inc., St. Louis, MO, USA; dimethyl sulfoxide (DMSO). Culture plastic dishes and plates (96-well) were purchased from Becton Dickinson (Franklin Lakes, NJ, USA).
Induction of apoptosis by [6c] in HSC-2 human oral squamous cell carcinoma cell line. PARP: Poly (ADP-ribose) polymerase; GAPDH: glyceraldehyde 3-phosphate dehydrogenase; Act D: actinomycin D.
Synthesis of test compounds. 2-(1H-Pyrrol-1-yl)-4H-1-benzopyran-4-one [1a], 7-methoxy-2-(1H-pyrrol-1-yl)-4H-1-benzopyran-4-one [1b], 6-methoxy-2-(1H-pyrrol-1-yl)-4H-1-benzopyran-4-one [1c], 2-(1H-pyrazol-1-yl)-4H-1-benzopyran-4-one [2a], 7-methoxy-2-(1H-pyrazol-1-yl)-4H-1-benzopyran-4-one [2b], 6-methoxy-2-(1H-pyrazol-1-yl)-4H-1-benzopyran-4-one [2c], 2-(1H-imidazol-1-yl)-4H-1-benzopyran-4-one [3a], 2-(1H-imidazol-1-yl)-7-methoxy-4H-1-benzopyran-4-one [3b], 2-(1H-imidazol-1-yl)-6-methoxy-4H-1-benzopyran-4-one [3c], 2-(1H-1,2,4-triazol-1-yl)-4H-1-benzopyran-4-one [4a], 7-methoxy-2-(1H-1,2,4-triazol-1-yl)-4H-1-benzopyran-4-one [4b], 6-methoxy-2-(1H-1,2,4-triazol-1-yl)-4H-1-benzopyran-4-one [4c], 2-(1H-1,2,3-triazol-1-yl)-4H-1-benzopyran-4-one [5a], 7-methoxy-2-(1H-1,2,3-triazol-1-yl)-4H-1-benzopyran-4-one [5b], 6-methoxy-2-(1H-1,2,3-triazol-1-yl)-4H-1-benzopyran-4-one [5c], 2-(1H-indol-1-yl)-4H-1-benzopyran-4-one [6a], 2-(1H-indol-1-yl)-7-methoxy-4H-1-benzopyran-4-one [6b], 2-(1H-indol-1-yl)-6-methoxy-4H-1-benzopyran-4-one [6c], 2-(1H-indazol-1-yl)-4H-1-benzopyran-4-one [7a], 2-(1H-indazol-1-yl)-7-methoxy-4H-1-benzopyran-4-one [7b], 2-(1H-indazol-1-yl)-6-methoxy-4H-1-benzopyran-4-one [7c], 2-(1H-benzimidazol-1-yl)-4H-1-benzopyran-4-one [8a], 2-(1H-benzimidazol-1-yl)-7-methoxy-4H-1-benzopyran-4-one [8b], 2-(1H-benzimidazol-1-yl)-6-methoxy-4H-1-benzopyran-4-one [8c] were synthesized by the conjugated addition reactions of 3-iodochromone derivatives with selected azoles, according to previous methods (4). All compounds were dissolved in DMSO at 40 mM and stored at −20°C before use.
Cell culture. Human normal oral mesenchymal cells (human gingival fibroblast, HGF; human periodontal ligament fibroblast, HPLF) were established from the first premolar tooth extracted from the lower jaw of a 12-year-old girl (7), and cells at 10-18 population-doubling levels were used in this study. Human OSCC cell lines [Ca9-22 (derived from gingival tissue); HSC-2 (derived from tongue)] were purchased from Riken Cell Bank (Tsukuba, Japan). All of these cells were cultured at 37°C in DMEM supplemented with 10% heat-inactivated FBS, 100 units/ml, penicillin G and 100 μg/ml streptomycin sulfate under a humidified 5% CO2 atmosphere.
Cytotoxic activity of 24 2-azolylchromone derivatives against oral malignant and non-malignant cells. Each value represents the mean of triplicate determinations.
Assay for cytotoxic activity. Cells were inoculated at 2.5×103 cells/0.1 ml in a 96-microwell plate. After 48 h, the medium was replaced with 0.1 ml of fresh medium containing different concentrations of single test compounds. Cells were incubated further for 48 h and the relative viable cell number was then determined by the MTT method (8). The relative viable cell number was determined by the absorbance of the cell lysate at 560 nm, using a microplate reader (Infinite F 50 R; TECAN, Kawasaki, Japan). Control cells were treated with the same amounts of DMSO and the cell damage induced by DMSO was subtracted from that induced by test agents. The concentration of compound that reduced the viable cell number by 50% (CC50) was determined from the dose–response curve and the mean value of CC50 for each cell type was calculated from triplicate assays.
Calculation of tumor-selectivity index (TS). TS was calculated using the following equation: TS=mean CC50 against three human normal cells/mean CC50 against for human OSCC cell lines [(D/B) in Table I]. Since both Ca9-22 and HGF cells were derived from the gingival tissue (9), the relative sensitivity of these cells was also compared [(C/A) in Table I]. We did not use human normal oral keratinocytes as controls, since many anticancer drugs showed potent cytotoxicity against normal keratinocytes by inducing apoptosis (8).
Calculation of potency-selectivity expression (PSE). PSE was calculated using the following equation: PSE=TS/CC50 against tumor cells ×100 (10) [that is, (D/B2) ×100 (HGF, HPLF vs. Ca9-22, HSC-2) and (C/A2) ×100 (HGF vs. Ca9-22 in Table I).
Source and explanation of chemical descriptors that correlate with cytotoxicity to tumor cells, normal cells and tumor specificity.
Estimation of CC50 values. Since the CC50 values had a distribution pattern close to a logarithmic normal distribution, we used the pCC50 (i.e., the −log CC50) for the comparison of the cytotoxicity between the compounds. The mean pCC50 values for normal cells and tumor cell lines were defined as N and T, respectively (10).
Calculation of chemical descriptors. The 3D-structure of each chemical structure (drawn by Marvin Sketch ver 16, ChemAxon, Budapest, Hungary, http://www.chemaxon.com) was optimized by CORINA Classic (Molecular Networks GmbH, Germany) and force-field calculations (amber-10: EHT) in Molecular Operating Environment (MOE) version 2015.1001 (Chemical Computing Group Inc., Quebec, Canada). The number of structural descriptors calculated from MOE and Dragon 7.0 (Kode srl., Pisa, Italy) after the elimination of overlapped descriptors were 291 and 2771, respectively. The 13 Dragon descriptors (11) and 1 MOE descriptor (12) listed in Table II were significantly correlated with T, N and T-N.
Western blot analysis. The cells were washed with PBS and processed for western blot analysis, as described previously (8). Antibodies against cleaved caspase-3 (Cell Signaling Technology Inc., Beverly, MP, USA) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH; Trevigen, Gaithersburg, MD, USA) were used as primary antibodies. As secondary antibodies, we used α-rabbit IgG (DAKO, Tokyo, Japan) antibodies which were conjugated with horseradish peroxidase.
Statistical treatment. The relation among cytotoxicity, tumor specificity index and chemical descriptors was investigated using simple regression analyses by JMP Pro version 12.2.0 (SAS Institute Inc., Cary, NC, USA). The significance level was set at p<0.05.
Results
Cytotoxicity. We synthesized a total of 24 2-azolylchromones, classified into the following eight sets of 1H-pyrrol-1-yl [1a-c], 1H-pyrazol-1-yl [2a-c], 1H-imidazol-1-yl [3a-c], 1H-1,2,4-triazol-1-yl [4a-c], 1H-1,2,3-triazol-1-yl [5a-c], 1H-indol-1-yl [6a-c], 1H-indazol-1-yl [7a-c] and 1H-benzimidazol-1-yl [8a-c] compounds. Among these eight groups, [6a-c] showed the highest cytotoxicity against human OSCC cell lines (Ca9-22, HGF) (mean CC50=1.5-6.3 μM), followed by [8a-c] (4.7-74.9 μM), [5a-c] (46.6-50.0 μM), [7a-c] (51.0->214 μM), [4a-c] (105.7->244 μM), [1a-c] (159.7->295 μM), [3a-c] (71.0->400μM) and [2a-c] (>293 μM).
On the other hand, all these groups of compounds generally showed much lower cytotoxicity against human normal oral cells (HGF and HPLF). [6a-c] showed slightly higher cytotoxicity (mean CC50=35.7->229 μM) than other groups: [1a-c] (251.2-288.7 μM), [2a-c] (300.2->400 μM), [3a-c] (>265->400 μM), [4a-c] (>223->254 μM), [5a-c] (>217->252 μM), [7a-c] (>232->400 μM) and [8a-c] (61.3->214 μM).
Tumor specificity. TS values were calculated for the compounds by dividing the average CC50 value towards normal cells into the average CC50 value towards cancer cell lines. When two human OSCC cell lines (Ca9-22, HGF) and two normal oral cells (HGF, HPLF) were used, [6a-c] had much higher TS values (24.1->37.5) than [1a-c] (<0.9-1.8), [2a-c] (><1.0), [3a-c] (><1.0->5.4), [4a-c] (><0.9->2.2), [5a-c] (>4.6->5.0), [7a-c] (><1.1->7.8) and [8a-c] (3.5-13.3). (D/B in Table I)
When Ca9-22 and HGF cells, both derived from gingival tissue, were used, [6a-c] again had higher TS valueS (8.3-40.8) than [1a-c] (1.4-2.6), [2a-c] (><1.0-1.7), [3a-c] (><1.0->6.3), [4a-c] (>2.8->4.9), [5a-c] (3.2-3.6), [7a-c] (><1.0->12.0) and [8a-c] (4.2->17.4) (C/A in Table I).
Determination of correlation coefficient between chemical descriptors and cytotoxicity [50% cytotoxic concentration (CC50)] of 24 2-azolylchromones against tumor cells (defined as T). The mean pCC50, i.e. the −log CC50 values for tumor cell lines, were defined as T.
PSE value. In order to identify the most promising compounds in terms of both good potencies and selective cytotoxicity, the PSE values were calculated. Among [6a-c], [6c] having a 2-(1H-indol-1-yl)-6-methoxy group, exhibited the highest cytotoxicity (CC50=1.5 μM). As expected, [6c] yielded the greatest PSE values:1622.5 and 1346.9, followed by [6b] (384.6 and 627.2) and [6a] (>614.3 and 121.2). On the other hand, PSE values of other groups were two-orders less in most cases (Table I).
Western blot analysis demonstrated that [6c] stimulated the cleavage of caspase-3 (Figure 2), suggesting the induction of apoptosis.
Computational analysis. We next performed the QSAR analysis of 2-azolylchromones in regards to their cytotoxicity against tumor cells and normal cells. Among a total of 3062 descriptors, 14 descriptors described below correlated well with cytotoxicity and TS (Table III).
Cytotoxicity of 2-azolylchromones against human OSCC cell lines was correlated with G3m (3D-shape and molecular weight) (r2=0.612, p<0.0001), G3e (3D-shape and electronegativity) (r2=0.601, p<0.0001), G3v (3D-shape and molecular volume) (r2=0.596, p<0.0001). Gm (symmetry and molecular weight) (r2=0.595, p<0.0001). G3p (3D-shape and polarizability) (r2=0.569, p<0.0001) and G3s (3D-shape and intrinsic state) (r2=0.566, p<0.0001) (Figure 3).
Cytotoxicity of 2-azolylchromones against human normal oral mesenchymal cells was correlated with SpMin8_Bh(s) (topological shape and intrinsic state) (r2=0.536, p<0.0001), Q_RPC-(negative partial charge) (r2=0.485, p=0.0002), G3s (3D-shape and intrinsic state) (r2=0.462, p=0.0003), G3e (3D-shape and electronegativity) (r2=0.453, p=0.0003), G3m (3D-shape and molecular weight) (r2=0.451, p=0.0003) and Gm (symmetry and molecular weight) (r2=0.448, p=0.0003) (Figure 4)
Properties of descriptors that significantly affect the cytotoxicity of compounds towards tumor cells (T), and normal cells (N), and their tumor specificity (T-N).
Determination of correlation coefficient between chemical descriptors and cytotoxicity [50% cytotoxic concentration (CC50)] of 24 2-azolylchromones against normal cells (defined as N). The mean pCC50, i.e. the −log CC50 values for normal cells, were defined as N.
Tumor specificity of 2-azolylchromones was correlated with Kp (3D-shape and polarizability) (r2=0.577, p<0.0001), P1p (3D-shape and polarizability) (r2=0.577, p<0.0001), Mor32i (3D-shape and ionization potential) (r2=0.570, p<0.0001), P2p (3D-shape and polarizability) (r2=0.564, p<0.0001), Mor32u (3D-shape) (r2=0.560, p<0.0001) and CATS2D_02_LL (topological shape and lipophilicity) (r2=0.555, p<0.0001) (Figure 5).
Discussion
The present study demonstrated that the introduction of indole ring to 4H-1-benzopyran-4-ones [6a-c] produced much higher tumor-specific compounds, as compared with those induced with pyrrole [1a-c], pyrazole [2a-c], imidazole [3a-c], 1,2,4-triazole [4a-c], 1,2,3-triazole [5a-c], indazole [7a-c] and benzimidazole [8a-c] rings. Especially, [6c] with methoxy group attached at C-6 position showed the highest PSE value (1622.5 and 1346.9) and induced apoptosis in HSC-2 cell OSCC cell line. This PSE value was about three times higher than that for (E)-3-(4-hydroxystyryl)-6-methoxy-4H-chromen-4-one (PSE=630) that showed high TS and lower keratinocyte specificity (13), and slightly higher than that of (E)-3-(4-chlorostyryl)-7-methoxy-2H-chromene (PSE=1277.5)(14). The cytotoxicity of [6c] against human oral keratinocytes should be investigated.
We also found that three sets of compounds with benzimidazol rings [8a-c] showed some TS (TS=3.5-13.0; PSE=4.7-273.8). It is interesting to state that cytotoxicity and PSE value of [8c] were 4.8-times (22.7/4.7) and 6.6-times (273.8/41.7) higher than those of [8b], although these two molecules have the same chemical formula with the only difference being in where the methoxy group attached. Similarly, cytotoxicity and PSE value of [6c] were 4.2-times (6.3/1.5) and 4.2-times (1,622.5/384.6) higher than those of [6b], both having the same chemical formula with the only difference being in the position of the methoxy group (Table I). This indicates that the chemical structure is not the sole determinant of TS. QSAR analysis demonstrated that TS of 2-azolylchromones was correlated with 3D-shape, polarizability, ionization potential and lipophilicity.
Anticancer drugs in clinical use have been reported to induce many side-effects in patients due to the lack of selectivity of the drugs to cancer tissues. Cytotoxicity of a compound is more valuable when it is selective. Although the chemical structures of all 24 compounds used in this study have already been reported (Table I) (4, 15), the present study is the first that investigated the relative potency of TS of all of these compounds. Chemical modification of [6c] may be a potential choice for designing a new type of anticancer drug.
Determination of correlation coefficient between chemical descriptors and tumor specificity of 24 2-azolylchromones (defined as T-N).
Acknowledgements
This work was partially supported by KAKENHI from the Japan Society for the Promotion of Science (JSPS) (15K08111, 16K11519).
Footnotes
This article is freely accessible online.
Conflicts of Interest
The Authors wish to confirm that there are no known conflicts of interest associated with this publication and there has been no significant financial support for this work that could have influenced its outcome.
- Received November 22, 2017.
- Revision received December 12, 2017.
- Accepted December 13, 2017.
- Copyright© 2018, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved