Compounds
1–6,
8, and
9 were synthesized as described before [
37], whereas compound
7 was prepared according to a previous study [
38]. All spectral data were in agreement with those reported. Herein, we report the first reductive amination of a PT, gypsogenin, being endowed with a unique carbaldehyde group. The latter was aminated by four different aromatic amines in the presence of sodium triacetoxyborohydride (
Figure 2).
Figure 2. The synthetic route for the preparation of new compounds (10–13).
MTT assay, the most widely used colorimetric assay for in vitro drug screening [
3], was conducted to assess the cytotoxic effects of compounds
1–13 and cisplatin (positive control) on U251, T98G, and U87 human glioblastoma cell lines. As indicated in
Table 1, compound
10 was identified as the most effective anti-glioma agent in this series. This compound exerted antitumor action towards U251, T98G, and U87 cells with IC
50 values of 5.82 µM, 8.19 µM, and 17.04 µM, respectively superior to cisplatin (IC
50 = 7.70 µM, 16.92 µM, and 20.90 µM, respectively). The potency order of cytotoxic effects of the other compounds on U251 cells was determined as compound
4 > compound
13 > compound
1 > compound
9 > compound
6 > compound
7 > compound
2. In accordance with the data recorded for U251 cells, compounds
4,
13, and
1 followed the same order for T98G and U87 cell lines. The IC
50 values of compounds
4,
13, and
1 for U251 cells were detected as 8.06 µM, 9.95 µM, and 13.18 µM, whereas their IC
50 values for T98G cells were found as 9.86 µM, 20.19 µM, and 20.54 µM, respectively. Compounds
4,
13, and
1 exhibited anticancer activity against U87 cells with IC
50 values of 19.54 µM, 21.71 µM, and 22.64 µM, respectively. As depicted in
Figure 3, compounds
4,
10, and
13 showed pronounced antiproliferative effects on U251, T98G, and U87 cells compared to cisplatin at varying concentrations. On the other hand, compounds
3,
5,
8,
11, and
12 displayed no significant anticancer activity against all tested GBM cell lines at 100 µM concentration.
Figure 3. The anticancer effects of compounds 4, 10, 13, and cisplatin at varying concentrations on U251 cells (a), T98G cells (b), U87 cells (c), Jurkat cells (d), and PBMCs (e). All descriptive data were expressed as the mean ± standard deviation (SD). All experiments were repeated three times.
Table 1. The cytotoxic effects of the compounds on U251, T98G, U87, and Jurkat cells, PBMCs.
Compound |
IC50 value (µM) |
SI1 |
U251 cells |
T98G cells |
U87 cells |
Jurkat cells |
PBMCs |
1 |
13.18 ± 3.19 |
20.54 ± 4.34 |
22.64 ± 6.75 |
|
|
|
2 |
24.00 ± 4.98 |
>100 |
>100 |
|
|
|
3 |
>100 |
>100 |
>100 |
|
|
|
4 |
8.06 ± 2.04 |
9.86 ± 2.21 |
19.54 ± 4.52 |
9.97 ± 3.24 |
21.91 ± 5.13 |
2.20 |
5 |
>100 |
>100 |
>100 |
|
|
|
6 |
16.68 ± 3.17 |
64.12 ± 7.36 |
79.70 ± 10.08 |
|
|
|
7 |
17.98 ± 2.23 |
61.11 ± 5.13 |
60.93 ± 8.87 |
|
|
|
8 |
>100 |
>100 |
>100 |
|
|
|
9 |
14.13 ± 3.41 |
56.55 ± 6.08 |
>100 |
|
|
|
10 |
5.82 ± 1.66 |
8.19 ± 2.42 |
17.04 ± 4.92 |
3.56 ± 1.45 |
28.12 ± 5.05 |
7.90 |
11 |
>100 |
>100 |
>100 |
|
|
|
12 |
>100 |
>100 |
>100 |
|
|
|
13 |
9.95 ± 2.04 |
20.19 ± 5.47 |
21.71 ± 6.09 |
12.08 ± 1.64 |
43.15 ± 8.32 |
3.57 |
Cisplatin |
7.70 ± 2.81 |
16.92 ± 3.95 |
20.90 ± 5.16 |
4.87 ± 2.00 |
34.67 ± 7.11 |
7.12 |
In order to determine the selectivity of the mode of anti-glioma action, the most effective anticancer agents (compounds 4, 10, and 13) were further screened for their cytotoxic effects on Jurkat human leukemic T-cells and human peripheral blood mononuclear cells (PBMCs). The selectivity of compound 10 was found as the most promising between Jurkat cells and PBMCs with a selectivity index (SI) value of 7.90. The anticancer effects of compounds 4, 10, and 13 compared to cisplatin at varying concentrations on Jurkat cells and PBMCs also supported this outcome (Figure 3).
As compound 10 was designated as the most potent and selective anti-glioma agent according to MTT results, its further apoptotic and EGFR inhibitory effects were also investigated to provide mechanistic insight. Using the annexin V/ethidium homodimer III staining method, the apoptotic effects of compound 10 on U251 cells were evaluated. This method indicates apoptosis and necrosis based on staining green and red, respectively. The results indicated that compound 10 boosted apoptosis in U251 cells with 10.29% similar to cisplatin (13.83%) (Figure 4).
Figure 4. Alteration of U251 cells following exposure to IC50 concentration of the control (DMSO), compound 10, and cisplatin (a) for 24 h. The percentage of alive (blue), apoptotic (green), necrotic or late apoptotic (both green and red), and necrotic (red) cells (b) was determined by analyzing 100 randomly chosen stained cells in each experiment. Quantification of apoptotic effects of compound 10 and cisplatin (c). Data from three independent experiments were expressed as mean ± standard deviation and p values were determined using the student’s test.
In continuation of our mechanistic research, the inhibitory effects of compound 10 on EGFR were analyzed due to the correlation between diminished EGFR signaling with increased anti-glioma activity. It was observed that compound 10 significantly inhibited EGFR with an IC50 value of 9.43 µM compared to erlotinib (IC50 = 0.06 µM), a first-generation EGFR tyrosine kinase inhibitor (Figure 5). Moreover, Figure 6 also highlighted the significant EGFR activity of compound 10 at 30 µM concentration compared to erlotinib. This outcome also pointed out that newly synthesized compounds displayed higher EGFR inhibition than our previously synthesized compounds. In order to explore the kinase selectivity profiling of compound 10, the inhibition of compound 10 at 30 µM concentration was examined on a large panel of tyrosine kinase enzymes including TK-1 (HER2, HER4, IGF1R, InsR, KDR, PDGFR-α, and PDGFR-β) and TK-2 (ABL1, BRK, BTK, CSK, FYN A, LCK, LYN B, and SRC) compared to erlotinib. Compound 10 showed the most potent inhibitory activity on the InsR followed by KDR, PDGFR-α, LCK, ABL1, HER2, CSK, PDGFR-β, and FYN A. This compound displayed no significant inhibition against HER4, IGF1R, BTK, BRK, LYN B, and SRC. Compound 10 and erlotinib exhibited similar and moderate HER2, ABL1, and FYN A inhibition, whereas compound 10 showed inhibitory effects on InsR, CSK, and LCK stronger than erlotinib. According to the results, compound 10 revealed a different kinase inhibitory profile than erlotinib as depicted in Figure 7. It can be concluded that compound 10 at 30 µM concentration exhibited the most selective inhibition against EGFR (approximately 2-fold stronger inhibition than InsR, the second promising tyrosine kinase target of compound 10) among all the tested tyrosine kinases.
Figure 5. The EGFR kinase inhibition of compound 10 and erlotinib at different concentrations. All descriptive data were expressed as the mean ± SD. All experiments were repeated three times.
Figure 6. The EGFR kinase inhibition of compounds 1–13 and erlotinib at 30 μM concentration. All descriptive data were expressed as the mean ± SD. All experiments were repeated three times.
Figure 7. The inhibition of a panel of tyrosine kinases by compound 10 and erlotinib at 30 µM concentration. All descriptive data were expressed as the mean ± SD. All experiments were repeated three times.
On the basis of its significant in vitro EGFR inhibitory potency, molecular docking studies were also carried out to understand the affinity of compound
10 to the adenosine triphosphate (ATP) binding site of the EGFR, which was acquired from the Protein Data Bank (PDB) server (PDB ID: 4HJO) [
39]. Molecular docking data demonstrated that compound
10 presented strong affinity at ionized state with favorable hydrogen bonding with Cys773, Pro770, and Lys704 by means of carboxylic acid, (phenylethyl)amino group, and hydroxyl substituent, respectively. This strong affinity could also be attributed to its interaction with Cys773 similar to erlotinib. On the other hand, its less in vitro EGFR inhibitory activity compared to that of erlotinib could be explained by the lack of key interaction with Met769 (
Figure 8).
Figure 8. Docking poses of compound 10 and erlotinib (a) and docking interactions of compound 10 (b) and erlotinib (c) in the ATP binding site of EGFR (PDB code: 4HJO). Yellow dashes: hydrogen bonding. Compound 10 and erlotinib were colored in yellow green, and pink, respectively.
As compared to costly and time-consuming absorption, distribution, metabolism, excretion (ADME) experimental procedures [
40], computational models are advantageous approaches to provide access to a set of rapid, yet robust predictive models for physicochemical, and pharmacokinetic properties [
10]. In this direction, we performed in silico predictions of some pharmacokinetic parameters of the compounds by the QikProp, a predictive ADME module within the Maestro suite produced by Schrödinger. As depicted in
Table 2, the brain/blood partition coefficient (QPlogBB) values of compounds
1–13 ranging from −1.175 to −0.029 were found within the specified limits (−3 to 1.2). The central nervous system (CNS) activity values of compounds
1–13 (−2 to 0) were in agreement within the range (−2 to 2). The QPlogPo/w value, which is a crucial parameter for membrane permeability, metabolism, bioavailability, the toxicity of molecules, and a ligand binding to the receptor [
41], was determined as 5.744, 5.453, and 5.745 for compounds
10,
11, and
12, respectively within the specified range (−2 to 6.5). The QPlogPo/w values of other compounds were detected out of limit. Taking into account the importance of hydrogen bonding with pivotal residues in the ATP binding site of EGFR, the number of donor (nHBD) and acceptor (nHBA) sites for hydrogen bonds were calculated. Appropriate nHBD (0 to 4) and nHBA (3.7 to 7.2) values of all compounds within the limits (0 to 6 and 2 to 20, respectively) also supported the outcomes of the molecular docking study. Solvent accessible surface area (SASA) is defined as the accessibility of the residue to the solvent; either it is between lipid or water accessibility and it is also essential to BBB permeability [
42]. The SASA values of all compounds were in an optimal range of the specified values (300 to 1000). Compounds
1–13 violated two parameters of Lipinski’s rule of five (maximum is four) and one parameter of Jorgensen’s rule of three (maximum is three).
Table 2. Predicted ADME properties of compounds 1–13.
Compound |
QPlogBB * (−3 to 1.2) |
CNS * (−2 to 2) |
QPlogPo/w * (−2 to 6.5) |
nHBD * (0 to 6) |
nHBA * (2 to 20) |
SASA * (300–1000) |
Rule of Five ** |
Rule of Three *** |
1 |
−0.595 |
−1 |
5.877 |
3 |
7.1 |
743.017 |
2 |
1 |
2 |
−0.388 |
0 |
8.081 |
1 |
3.7 |
813.580 |
2 |
1 |
3 |
−0.529 |
0 |
7.171 |
1 |
5.7 |
827.146 |
2 |
1 |
4 |
−0.574 |
0 |
6.781 |
2 |
5.4 |
756.322 |
2 |
1 |
5 |
−0.120 |
0 |
7.783 |
1 |
3.7 |
746.935 |
2 |
1 |
6 |
−0.105 |
0 |
7.829 |
1 |
3.7 |
741.337 |
2 |
1 |
7 |
−0.569 |
0 |
6.968 |
0 |
4.7 |
811.227 |
2 |
1 |
8 |
−0.035 |
0 |
8.795 |
0 |
4.7 |
799.503 |
2 |
1 |
9 |
−1.175 |
−2 |
6.275 |
0 |
6.2 |
780.099 |
2 |
1 |
10 |
−0.061 |
−1 |
5.744 |
3 |
5.2 |
863.864 |
2 |
1 |
11 |
−0.038 |
−1 |
5.453 |
3 |
5.2 |
848.674 |
2 |
1 |
12 |
−0.029 |
−1 |
5.745 |
3 |
5.2 |
863.850 |
2 |
1 |
13 |
−1.057 |
−2 |
6.866 |
4 |
7.2 |
890.067 |
2 |
1 |
* QPlogBB: Brain/blood partition coefficient, CNS: Predicted central nervous system activity. QPlogPo/w: Predicted octanol/water partition coefficient. nHBD and nHBA: Estimated number of hydrogen bonds that would be donated and accepted, respectively, by the solute to water molecules in an aqueous solution. Values are averages taken over a number of configurations, so they can be non-integer. SASA: Total solvent accessible surface area in square angstroms using a probe with a 1.4 Å radius. ** Rule of Five: Number of violations of Lipinski’s rule of five. *** Rule of Three: Number of violations of Jorgensen’s rule of three.