Open Access Article
Em Canh Pham
ab,
Song-Thuong Nguyena,
Kim Anh Thi Lec,
Tuoi Thi Hong Doa and
Tuyen Ngoc Truong
*a
aSchool of Pharmacy, University of Medicine and Pharmacy at Ho Chi Minh City, 700000 Ho Chi Minh City, Vietnam. E-mail: truongtuyen@ump.edu.vn
bFaculty of Pharmacy, Hong Bang International University, 700000 Ho Chi Minh City, Vietnam
cUMP Science and Technology Center, University of Medicine and Pharmacy at Ho Chi Minh City, 700000 Ho Chi Minh City, Vietnam
First published on 28th October 2025
Twenty symmetrical chlorophenylamino-s-triazine derivatives were designed and synthesized by reflux (RF) and microwave-assisted (MW) methods. The MW method achieved superior yields (88–95%) in less time (15–30 min) compared to RF (78–86%, 12–24 h), particularly for 3-Cl and 3,4-diCl derivatives with piperidine or diethylamine, due to rapid, uniform heating, enhancing nucleophilic substitution and minimizing side reactions. In particular, compounds 2c (IC50 = 4.14 μM for MCF7, 7.87 μM for C26), 3c (IC50 = 4.98 μM for MCF7, 3.05 μM for C26), and 4c (IC50 = 6.85 μM for MCF7, 1.71 μM for C26) exhibited potent cytotoxic activity with 4c (2,4-diCl, pyrrolidine) surpassing paclitaxel (IC50 = 2.30 μM for C26), and 3c (3,4-diCl, pyrrolidine) rivaling global analogs. Compounds 2f (IC50 = 11.02 μM for MCF7, 4.62 μM for C26) and 3f (IC50 = 5.11 μM for MCF7, 7.10 μM for C26) also showed strong cytotoxicity. QSAR analysis revealed that electron-withdrawing groups (chloro, dichloro) and pyrrolidine enhance C26 potency via improved lipophilicity and π–π stacking, outperforming piperazine and morpholine. Pharmacokinetically, 2c, 3c, and 4c matched Bimiralisib's absorption profiles, surpassing Gefitinib, with 4c showing superior metabolic stability. Compounds 2f and 3c emerged as promising multi-targeted kinase inhibitors, with binding affinities (−7.8 to −9.1 kcal mol−1) closely rivaling Gefitinib, Pazopanib, and Bimiralisib for EGFR, VEGFR2, and PI3K, driven by balanced polar and hydrophobic interactions. Therefore, these findings underscore the potential of 2f and 3c as multi-targeted kinase inhibitors, warranting further mechanistic studies and structural optimization to enhance MCF7 efficacy and reduce toxicity for clinical advancement.
A common approach in the synthesis of symmetrical s-triazine derivatives is the stepwise substitution of 2,4,6-halogenated s-triazines, such as cyanuric chloride, with nucleophiles such as amines, alcohols, or thiols. This method provides precise control over the substitution patterns by exploiting the differential reactivity of chlorine atoms at different temperatures, enabling the sequential introduction of identical substituents to achieve symmetry. In addition, the synthesis from cyanuric chloride stands out due to its operational simplicity, cost-effectiveness, and versatility, especially the possibility of microwave-assisted synthesis to enhance reaction efficiency and yield. Furthermore, this method facilitates the production of high yields of symmetrically substituted s-triazines with excellent reproducibility, making it particularly advantageous for both laboratory-scale experiments and industrial applications.23,29
![]() | ||
| Fig. 2 Rational design of symmetrical chlorophenylamino-s-triazine derivatives from anticancer drugs and the potent compound. | ||
In addition, the strategic incorporation of chlorine (Cl) substituents in the rational drug design of symmetrical s-triazine derivatives plays a key role in enhancing anticancer potency. Compared to other halogens, chlorine strikes an optimal balance between reactivity and stability, avoiding the excessive instability of fluorine or the bulkiness of bromine, making it ideal for iterative drug optimization. The electron-withdrawing nature of chlorine enhances the electrophilicity of the benzene ring, thereby facilitating interactions with nucleophilic sites on biological targets, such as DNA or proteins involved in cancer cell proliferation. This modification influences the lipophilicity of the compound, improving bioavailability and cellular uptake, which are important factors for effective anticancer agents. Furthermore, several studies have demonstrated that chlorine substitution at specific positions on the benzene ring, particularly para- or meta-, optimizes cytotoxic activity against various cancer cell lines, including breast and lung carcinomas.30–32 These findings underscore the strategic importance of chlorine in fine-tuning the physicochemical and pharmacological properties of phenylamino-s-triazine derivatives for anticancer applications.
The strategic design of symmetrical s-triazine derivatives, incorporating one chlorophenylamino group and two saturated cyclic amino substituents, represents a robust approach to developing potent anticancer agents, driven by their tailored chemical and pharmacological properties. The chlorophenylamino moiety, characterized by its electron-withdrawing chlorine atom, enhances the lipophilicity of the triazine core and fine-tunes its electronic profile, enabling precise interactions with hydrophobic regions of target proteins, such as receptor tyrosine kinases. The incorporation of saturated cyclic amino groups, such as piperidine, piperazine, or morpholine, optimizes hydrogen-bonding interactions and steric fit, thereby improving binding affinity and target selectivity. Key therapeutic targets for these derivatives include the PI3K/mTOR signaling pathway, as exemplified by Gedatolisib, where cyclic amines enhance aqueous solubility and receptor engagement, demonstrating efficacy in ovarian and breast cancers.33,34 Additionally, epidermal growth factor receptor (EGFR) and vascular endothelial growth factor receptor 2 (VEGFR2) are promising targets, as the chlorophenylamino group may facilitate specific binding to kinase domains, disrupting oncogenic signaling cascades.35–37 The balanced lipophilicity and solubility conferred by the combination of chlorophenylamino and cyclic amino groups enhance the derivatives' ability to circumvent multidrug resistance mechanisms, such as efflux pumps. Leveraging structure–activity relationship (SAR) analyses and computational modeling, these s-triazine derivatives can be optimized for specific oncogenic pathways, providing a versatile platform for the development of next-generation anticancer therapeutics.
This study aimed to synthesize a series of symmetrical chlorophenylamino-s-triazine derivatives incorporating diverse saturated cyclic amino groups and to assess their anticancer efficacy against MCF7 (human breast cancer) and C26 (colon carcinoma) cell lines. The synthesized compounds were evaluated for their cytotoxic potential, with promising candidates selected for further in silico molecular docking and ADMET (absorption, distribution, metabolism, excretion, and toxicity) analyses. These investigations seek to elucidate the compounds' interactions with key molecular targets and their pharmacokinetic profiles, providing insights into their therapeutic potential and drug-likeness.
![]() | ||
| Scheme 1 Synthesis of chlorophenylamino-s-triazine derivatives (MW – microwave irradiation, RF – reflux, THF – tetrahydrofuran). | ||
The reaction yields of chlorophenylamino-s-triazine derivatives using RF and MW methods showed marked differences between different R and R1 substituents. The MW method consistently outperformed the RF method by 8–13% and achieved higher yields (88–95%) compared to reflux yields (78–86%) across all compounds (2a–4f) (Table 1). The microwave method's efficacy stems from its ability to deliver rapid, uniform heating, which accelerates reaction kinetics, enhances nucleophilic substitution on the s-triazine core, and minimizes side reactions. This results in improved product purity and yield, particularly for compounds with R groups of 3-Cl (e.g., 2a, 95%) and 3,4-diCl (e.g., 3a, 95%) paired with R1 groups of piperidine (Piper) or diethylamine (Dieth). For group R, compounds with 3-Cl (89–95%), 3,4-diCl (88–95%), and 2,4-diCl (88–93%) substitutions showed similar yield trends. However, regarding R1 groups, piperidine (Piper) and diethylamine (Dieth) derivatives (e.g., 2a, 2g, 3a, 3g, 4a) consistently exhibit higher yields in both methods, likely due to their favorable steric and electronic properties, which enhance nucleophilic substitution on the s-triazine core. Conversely, morpholine (Mor) derivatives (e.g., 2d, 3d, 4d) show the lowest yields (78–88%), potentially due to reduced nucleophilicity (Table 1). These findings align closely with those reported by Al-Zaydi et al. (2017), who utilized 4-carboxyaniline as the starting material for tri-substitution reactions with saturated cyclic amines (piperidine and morpholine) in a 1
:
1 mixture of 1,4-dioxane and water, using sodium carbonate (Na2CO3) as the base. Their study achieved yields of 74.5–85.8% over 8–10 h under RF conditions, and 88–93.1% in 10 min using MW irradiation at 400 W.38 In particular, our published research highlights the superior yields of the MW method in synthesizing tri-substituted s-triazine derivatives (R = 4-Cl, 4-F, 4-OCH3, 4-CH3, 4-NO2; R1 = Piper, Mor), achieving excellent yields (91–98%) and rapid reaction times (15–30 min). In contrast, the RF method required significantly longer reaction times (12–24 h) and lower yields (80–88%).23 Therefore, these findings underscore the advantages of the MW method (“green” approach) in terms of efficiency, speed, and environmental protection, offering a scalable, efficient, and environmentally friendly approach for producing potential anticancer agents.
The chemical structures of chlorophenylamino-s-triazine derivatives were suitably elucidated by 1H NMR, 13C NMR, and MS spectroscopy. The 1H and 13C NMR spectra of chlorophenylamino-s-triazine derivatives (2a–2g, 3a–3g, 4a–4f) provide critical insights into their structural characteristics, highlighting the influence of varying R and R1 substituents on the s-triazine core. In the 1H NMR spectra, the –NH– proton, linking the chlorophenyl moiety to the triazine ring, consistently appears as a singlet in the range of 8.05–10.76 ppm, with compounds bearing 3,4-dichlorophenyl (9.03–9.41 ppm) and 3-chlorophenyl (9.13–10.76 ppm) groups showing slight downfield shifts compared to 2,4-dichlorophenyl derivatives (8.00–8.31 ppm). This variation reflected the electron-withdrawing effects of additional chlorine atoms, which modulate the electron density around the –NH– group. Aromatic protons (HAr) on the chlorophenyl ring exhibited characteristic splitting patterns, with 3-chlorophenyl derivatives showing a triplet and doublets at δ 6.90–8.16 ppm, while 3,4-dichlorophenyl and 2,4-dichlorophenyl compounds display simplified patterns due to increased substitution, as seen in the doublet at δ 8.35 ppm for 3a. The R1 cyclic amino groups significantly influence aliphatic signals; for instance, pyrrolidine (2a, 3a, 4a) showed distinct –CH2– triplets at δ 3.46–3.49 ppm, whereas morpholine (2d, 3d, 4d), piperazine (2e, 3e, 4e) and its derivatives (2b, 3b, 4b) exhibited additional oxygen/nitrogen-adjacent –CH2– signals at δ 3.59–3.70 ppm, reflecting its heterocyclic nature. Piperidine derivatives (2f, 3f, 4f) display broader –CH2– signals at δ 1.47–3.70 ppm, indicating greater conformational flexibility. In the 13C NMR spectra, the triazine ring carbons resonated at δ 163.1–165.0 ppm, with minor shifts attributed to the electronic effects of R1 substituents. For example, morpholine derivatives (2d, 3d, 4d) showed a characteristic carbon signal at δ 65.9–66.3 ppm for oxygen-linked –CH2–, absent in piperidine or pyrrolidine analogs. The chlorophenyl carbons appear at δ 117.3–142.6 ppm, with 3,4-dichlorophenyl compounds (e.g., 3a) showing a downfield shift (δ 141.0 ppm) compared to 2,4-dichlorophenyl analogs (e.g., 4a, δ 135.4 ppm), reflecting differences in chlorine positioning. Aliphatic carbons of R1 groups, such as pyrrolidine (δ 24.3–24.8 ppm) and piperidine (δ 24.3–25.8 ppm), are consistent across series, while methyl-substituted piperazines (2b, 3b, 4b) introduced additional signals at δ 42.5–54.5 ppm. Moreover, the mass spectrometry analysis revealed the molecular ion peak (M, m/z) for compounds 2–4, corroborating their proposed structures. Notably, all chlorophenylamino-s-triazine derivatives exhibited physicochemical properties, including molecular weights (MWt < 500), consistent with Lipinski's rule of five (Table 1). These characteristics suggest their potential as promising candidates for further drug development.
| Entry | R group | Code | Physicochemical parametersa | Yield (%) | |||
|---|---|---|---|---|---|---|---|
| R | R1 | RF | MW | ||||
a Calculated using SwissADME, RF – reflux method (/conventional heating), MW – microwave-assisted methods, MWt – molecular weight, nHA – number of hydrogen bond acceptors, nHD – number of hydrogen bond donors, nRot – number rotatable bonds, TPSA – topological polar surface area (Angstroms squared), MR – molar refractivity, log P – log Po/w (XLOGP3), log S – log S (ESOL). |
|||||||
| 1 | 3-Cl | Piper | 2a | MWt: 372.90 | MR: 113.28 | 84 | 95 |
| nHA: 3 | Log P: 5.08 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −5.43 |
||||||
| 2 | 3-Cl | 4-MePiper | 2b | MWt: 400.95 | MR: 122.89 | 82 | 91 |
| nHA: 3 | Log P: 5.95 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −6.13 |
||||||
| 3 | 3-Cl | Pyrro | 2c | MWt: 344.84 | MR: 103.66 | 82 | 94 |
| nHA: 3 | Log P: 4.37 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −4.84 |
||||||
| 4 | 3-Cl | Mor | 2d | MWt: 376.84 | MR: 105.83 | 84 | 92 |
| nHA: 5 | Log P: 2.64 |
||||||
| nHD: 1 | TPSA: 75.64 | ||||||
| nRot: 4 | Log S: −3.92 |
||||||
| 5 | 3-Cl | 4-MePipera | 2e | MWt: 402.92 | MR: 126.90 | 79 | 89 |
| nHA: 5 | Log P: 3.01 |
||||||
| nHD: 1 | TPSA: 63.66 | ||||||
| nRot: 4 | Log S: −4.29 |
||||||
| 6 | 3-Cl | Pipera | 2f | MWt: 374.87 | MR: 117.10 | 82 | 95 |
| nHA: 5 | Log P: 2.08 |
||||||
| nHD: 3 | TPSA: 81.24 | ||||||
| nRot: 4 | Log S: −3.55 |
||||||
| 7 | 3-Cl | Dieth | 2g | MWt: 348.87 | MR: 102.02 | 85 | 94 |
| nHA: 3 | Log P: 4.86 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 8 | Log S: −4.91 |
||||||
| 8 | 3,4-diCl | Piper | 3a | MWt: 407.34 | MR: 118.29 | 82 | 95 |
| nHA: 6 | Log P: 5.71 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −6.03 |
||||||
| 9 | 3,4-diCl | 4-MePiper | 3b | MWt: 435.39 | MR: 127.90 | 79 | 90 |
| nHA: 3 | Log P: 6.58 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −6.73 |
||||||
| 10 | 3,4-diCl | Pyrro | 3c | MWt: 379.29 | MR: 108.67 | 80 | 90 |
| nHA: 3 | Log P: 4.99 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −5.43 |
||||||
| 11 | 3,4-diCl | Mor | 3d | MWt: 411.29 | MR: 110.84 | 78 | 88 |
| nHA: 5 | Log P: 3.27 |
||||||
| nHD: 1 | TPSA: 75.64 | ||||||
| nRot: 4 | Log S: −4.51 |
||||||
| 12 | 3,4-diCl | 4-MePipera | 3e | MWt: 437.37 | MR: 131.91 | 86 | 94 |
| nHA: 5 | Log P: 3.64 |
||||||
| nHD: 1 | TPSA: 63.66 | ||||||
| nRot: 4 | Log S: −4.89 |
||||||
| 13 | 3,4-diCl | Pipera | 3f | MWt: 409.32 | MR: 122.11 | 83 | 92 |
| nHA: 5 | Log P: 2.71 |
||||||
| nHD: 3 | TPSA: 81.24 | ||||||
| nRot: 4 | Log S: −4.15 |
||||||
| 14 | 3,4-diCl | Dieth | 3g | MWt: 383.32 | MR: 107.03 | 83 | 95 |
| nHA: 3 | Log P: 5.48 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 8 | Log S: −5.50 |
||||||
| 15 | 2,4-diCl | Piper | 4a | MWt: 407.34 | MR: 118.29 | 82 | 90 |
| nHA: 3 | Log P: 5.71 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −6.03 |
||||||
| 16 | 2,4-diCl | 4-MePiper | 4b | MWt: 435.39 | MR: 127.90 | 83 | 92 |
| nHA: 3 | Log P: 6.58 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −6.73 |
||||||
| 17 | 2,4-diCl | Pyrro | 4c | MWt: 379.29 | MR: 108.67 | 81 | 90 |
| nHA: 3 | Log P: 4.99 |
||||||
| nHD: 1 | TPSA: 57.18 | ||||||
| nRot: 4 | Log S: −5.43 |
||||||
| 18 | 2,4-diCl | Mor | 4d | MWt: 411.29 | MR: 110.84 | 78 | 88 |
| nHA: 5 | Log P: 3.27 |
||||||
| nHD: 1 | TPSA: 75.64 | ||||||
| nRot: 4 | Log S: −4.51 |
||||||
| 19 | 2,4-diCl | 4-MePipera | 4e | MWt: 437.37 | MR: 131.91 | 82 | 93 |
| nHA: 5 | Log P: 3.64 |
||||||
| nHD: 1 | TPSA: 63.66 | ||||||
| nRot: 4 | Log S: −4.89 |
||||||
| 20 | 2,4-diCl | Pipera | 4f | MWt: 409.32 | MR: 122.11 | 81 | 90 |
| nHA: 5 | Log P: 2.71 |
||||||
| nHD: 3 | TPSA: 81.24 | ||||||
| nRot: 4 | Log S: −4.15 |
||||||
| Entry | Code | Cancer cell line | |
|---|---|---|---|
| MCF7 (IC50, μM) | C26 (IC50, μM) | ||
| a All data are presented as mean ± standard deviation (SD), PTX – paclitaxel, MCF7 – human breast cancer cell line, C26 – colon carcinoma cell line, the values in bold highlight the best compounds with the best IC50 values compared to the positive control, * – statistically significant (p < 0.05) compared to reference drug PTX. | |||
| 1 | 2a | 28.20 ± 1.68 | 23.05 ± 1.25 |
| 2 | 2b | 24.37 ± 0.66 | 20.53 ± 0.65 |
| 3 | 2c | 4.14 ± 1.06 | 7.87 ± 0.96 |
| 4 | 2d | 38.02 ± 3.85 | 21.54 ± 0.81 |
| 5 | 2e | 19.64 ± 1.37 | 47.53 ± 1.03 |
| 6 | 2f | 11.02 ± 0.68 | 4.62 ± 0.65 |
| 7 | 2g | 36.67 ± 1.85 | 58.17 ± 2.19 |
| 8 | 3a | 36.02 ± 2.26 | 29.14 ± 2.63 |
| 9 | 3b | 36.33 ± 0.76 | 21.77 ± 1.10 |
| 10 | 3c | 4.98 ± 0.58 | 3.05 ± 0.61 |
| 11 | 3d | 28.36 ± 1.88 | 30.71 ± 4.25 |
| 12 | 3e | 37.38 ± 3.47 | 24.25 ± 1.34 |
| 13 | 3f | 5.11 ± 0.31 | 7.10 ± 0.48 |
| 14 | 3g | 55.62 ± 5.30 | 37.58 ± 0.88 |
| 15 | 4a | 40.85 ± 0.41 | 46.55 ± 1.42 |
| 16 | 4b | 25.09 ± 1.39 | 48.55 ± 1.30 |
| 17 | 4c | 6.85 ± 0.34 | 1.71 ± 0.88* |
| 18 | 4d | 43.03 ± 0.96 | 39.15 ± 2.96 |
| 19 | 4e | 10.05 ± 0.36 | 16.38 ± 0.63 |
| 20 | 4f | 17.04 ± 1.00 | 17.15 ± 0.83 |
| 21 | PTX | 2.48 ± 0.51 | 2.30 ± 0.27 |
In the series 3 compounds, 3c and 3f stood out as the most potent. Compound 3c exhibited IC50 values of 4.98 ± 0.58 μM (MCF7) and 3.05 ± 0.61 μM (C26), with its C26 activity closely rivaling PTX. Compound 3f showed IC50 values of 5.11 ± 0.31 μM (MCF7) and 7.10 ± 0.48 μM (C26), also indicating strong cytotoxicity. Other compounds, including 3a (36.02 ± 2.26 μM for MCF7, 29.14 ± 2.63 μM for C26), 3b (36.33 ± 0.76 μM for MCF7, 21.77 ± 1.10 μM for C26), 3d (28.36 ± 1.88 μM for MCF7, 30.71 ± 4.25 μM for C26), 3e (37.38 ± 3.47 μM for MCF7, 24.25 ± 1.34 μM for C26), and 3g (55.62 ± 5.30 μM for MCF7, 37.58 ± 0.88 μM for C26), showed moderate to low activity, with IC50 values 10- to 20-fold higher than PTX. The potency of 3c and 3f suggests favorable molecular interactions, possibly involving enhanced binding affinity to apoptotic or mitotic pathways.
Series 4 compounds revealed 4c and 4e as the most effective. Compound 4c was particularly notable, with IC50 values of 6.85 ± 0.34 μM (MCF7) and 1.71 ± 0.88 μM (C26), the latter surpassing PTX's potency against C26. Compound 4e showed IC50 values of 10.05 ± 0.36 μM (MCF7) and 16.38 ± 0.63 μM (C26), indicating moderate activity. Other compounds, such as 4a (40.85 ± 0.41 μM for MCF7, 46.55 ± 1.42 μM for C26), 4b (25.09 ± 1.39 μM for MCF7, 48.55 ± 1.30 μM for C26), 4d (43.03 ± 0.96 μM for MCF7, 39.15 ± 2.96 μM for C26), and 4f (17.04 ± 1.00 μM for MCF7, 17.15 ± 0.83 μM for C26), exhibited lower potency. The exceptional activity of 4c, particularly against C26, suggests a highly optimized structure for targeting colon carcinoma cells, potentially through mechanisms akin to PTX's disruption of microtubule dynamics.
PTX remains the gold standard with IC50 values of 2.48 ± 0.51 μM (MCF7) and 2.30 ± 0.27 μM (C26). Among the tested compounds, 4c stands out as the only compound surpassing PTX's potency against C26, while 3c and 2f closely approach it. Against MCF7, compounds 2c, 3c, and 3f showed promising activity but remain less potent than PTX. The weaker performance of most compounds, particularly 2g, 3g, 4a, and 4d, underscores the challenge of achieving PTX's broad-spectrum efficacy. Structural modifications in 2c, 3c, 3f, and 4c likely enhance their cytotoxic potential, possibly by improving solubility, cellular uptake, or target specificity, but further studies are needed to elucidate their mechanisms and optimize their activity to match or exceed PTX.
In conclusion, compounds 2c (3-Cl, pyrrrolidine), 2f (3-Cl, piperazine), 3c (3,4-diCl, pyrrrolidine), 3f (3,4-diCl, piperazine), and 4c (2,4-diCl, pyrrrolidine) exhibited significant cytotoxic potential, with 4c showing superior activity against C26 compared to PTX (Fig. 3). The results also suggested that the R1 substituents such as pyrrolidine and piperazine in these derivatives may be responsible for their potential cytotoxic activity. These compounds warrant further investigation for structural optimization and mechanism studies to enhance their efficacy, especially against MCF7 cells, where PTX still retains a clear advantage.
![]() | ||
| Fig. 3 The IC50 values (μM) of potential chlorophenylamino-s-triazine derivatives compared to paclitaxel (PTX) (MCF7 – human breast cancer cell line, C26 – colon carcinoma cell line). | ||
:
0.50 μM for HCT-116),39 with the ortho-chlorine likely optimizing steric and electronic interactions for exceptional C26 activity. QSAR analysis reveals that electron-withdrawing groups (e.g., chloro, dichloro) on the phenyl ring enhance cytotoxicity, particularly against C26, as seen in 4-halogeno and 4-nitro substituted triazines.23 Pyrrolidine substitutions consistently improve potency over piperazine, likely due to increased lipophilicity and reduced steric hindrance, while dichlorophenyl groups amplify activity compared to monochlorophenyl, correlating with higher electron-withdrawing capacity. Compared to sulfaguanidine-triazines (IC50: 14.8–33.2 μM for MCF7),21 these compounds exhibited superior potency, though less toxic than doxorubicin (IC50: 0.42 μM),18 suggesting a favorable therapeutic index.
Compared more specifically with our published research, in examining the cytotoxic activities of five chlorophenylamino-s-triazine derivatives (2c, 2f, 3c, 3f, and 4c) against MCF7 and C26 cell lines, these compounds generally exhibited superior potency compared to the three phenylamino-s-triazine analogs (P-2e (4-nitrophenyl and piperidine substituted triazine), P-3a (4-chlorophenyl and morpholine substituted triazine), and P-3e (4-nitrophenyl and morpholine substituted triazine)), as evidenced by their lower IC50 values. For the MCF7 line, the potent chlorophenyl series displayed IC50 ranges from 4.14 ± 1.06 μM (2c) to 11.02 ± 0.68 μM (2f), markedly outperforming the phenyl series, which spans 13.74 ± 1.96 μM (P-3e) to 42.40 ± 4.48 μM (P-3a), suggesting that chlorine substitution on the phenyl ring enhances selectivity and efficacy against this estrogen receptor-positive breast cancer model, potentially through improved lipophilicity or π–π stacking interactions with cellular targets. In contrast, on the C26 line, the chlorophenyl compounds yield IC50 values from 1.71 ± 0.88 μM (4c) to 7.87 ± 0.96 μM (2c), overlapping with but not consistently surpassing the phenyl analogs (1.21 ± 0.47 μM for P-3a to 14.66 ± 1.70 μM for P-3e), indicating a more variable response where certain chlorine positions, such as the 2,4-dichloro motif in 4c, confer exceptional activity possibly via steric hindrance or altered electron withdrawal affecting membrane permeability.
In summary, QSAR analysis of chlorophenylamino-s-triazine derivatives underscores the pivotal role of strategic substitutions in modulating cytotoxic efficacy against MCF7 breast cancer and C26 colon cancer cell lines. Electron-withdrawing chloro groups on the phenyl ring, particularly in di-substituted configurations such as 3,4-dichloro (3c) and 2,4-dichloro (4c), significantly enhance potency against C26. Moreover, pyrrolidine moieties consistently outperform piperazine counterparts and tend to favor lower IC50 on MCF7 compared to piperazine, as seen in the superior activity of 2c and 3c over 2f and 3f, attributable to enhanced molecular flexibility, reduced steric hindrance, and stronger hydrophobic engagements, while piperazine confers selective hydrogen-bonding advantages in colon cancer models. While less potent than doxorubicin, their favorable therapeutic indices position compounds like 4c and 3c as promising scaffolds for lead optimization, warranting deeper mechanistic investigations into tubulin-binding dynamics and in vivo pharmacokinetics to advance anticancer therapeutics.
| Parameter | 2c | 2f | 3c | 3f | 4c | Ged | Bim | |||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
a Ged – Gedatolisib, Bim – Bimiralisib, Caco-2 permeability (optimal: higher than −5.15 log unit), MDCK permeability (low permeability: <2 × 10−6 cm s−1, medium permeability: 2–20 × 10−6 cm s−1, high passive permeability: >20 × 10−6 cm s−1), PAMPA – the experimental data for Peff was logarithmically transformed (log Peff < 2: low-permeability, log Peff > 2.5: high-permeability), Pgp – P-glycoprotein, HIA – human intestinal absorption (−: ≥30%, +: < 30%), F: bioavailability (+: < percent value, −: ≥ percent value), PPB: plasma protein binding (optimal: < 90%), VD: volume distribution (optimal: 0.04–20 L kg−1), BBB: blood–brain barrier penetration, Fu: the fraction unbound in plasms (low: <5%, middle: 5–20%, high: > 20%), CL : Clearance (low: < 5 mLmin−1kg, moderate: 5–15 mL min−1 kg−1, high: > 15 mL min−1 kg−1), T1/2 (ultra-short half-life drugs: 0.5 – < 1 h; short half-life drugs: 1–4 h; intermediate short half-life drugs: 4–8 h; long half-life drugs: >8 h), hERG blockers (IC50 ≤ 10 μM or ≥ 50% inhibition at 10 μM were classified as hERG +, IC50 > 10 μM or <50% inhibition at 10 μM were classified as hERG −–), DILI: drug-induced liver injury, rat oral acute toxicity (0: low-toxicity > 500 mg kg−1, 1: high-toxicity < 500 mg kg−1), FDAMDD – maximum recommended daily dose, BCF – bioconcentration factors, IGC50 – tetrahymena pyriformis 50 percent growth inhibition Concentration, LC50FM – 96 h fathead minnow 50 percent lethal concentration, LC50DM – 48 h daphnia magna 50 percent lethal concentration. The output value is the probability of being inhibitor/substrate/active/positive/high-toxicity/sensitizer/carcinogens/corrosives/irritants (category 1) or non-inhibitor/non-substrate/inactive/negative/low-toxicity/non-sensitizer/non-carcinogens/noncorrosives/nonirritants (category 0). For the classification endpoints, the prediction probability values are transformed into six symbols: 0–0.1(−−−), 0.1–0.3(−−), 0.3–0.5(−), 0.5–0.7(+), 0.7–0.9(++), and 0.9–1.0(+++). Additionally, the corresponding relationships of the three labels are as follows: E − excellent, M – medium, P – poor. |
||||||||||||||
| Druglikeness | ||||||||||||||
| Lipinski | Yes | Yes | Yes | Yes | Yes | No | Yes | |||||||
| Ghose | Yes | No | Yes | No | Yes | No | Yes | |||||||
| Veber | Yes | Yes | Yes | Yes | Yes | Yes | Yes | |||||||
| Egan | Yes | Yes | Yes | Yes | Yes | Yes | Yes | |||||||
| Muegge | Yes | Yes | Yes | Yes | Yes | No | Yes | |||||||
| Pfizer | No | Yes | No | Yes | No | Yes | Yes | |||||||
| Bioavailability score | 0.55 | 0.55 | 0.55 | 0.55 | 0.55 | 0.17 | 0.55 | |||||||
| SAscore | <6 | E | <6 | E | <6 | E | <6 | E | <6 | E | <6 | E | <6 | E |
![]() |
||||||||||||||
| Absorption | ||||||||||||||
| Caco-2 permeability | −4.883 | E | −5.283 | P | −5.105 | E | −5.569 | P | −4.942 | E | −5.322 | P | −4.546 | E |
| MDCK permeability | 0 | E | 0 | E | 0 | E | 0 | E | 0 | E | 0.0 | E | 0 | E |
| PAMPA | −−− | E | −− | E | −−− | E | − | M | −−− | E | −−− | E | −−− | E |
| Pgp-inhibitor | +++ | P | −−− | E | +++ | P | −−− | E | +++ | P | +++ | P | +++ | P |
| Pgp-substrate | −−− | E | +++ | P | −−− | E | +++ | P | −−− | E | +++ | P | −−− | E |
| HIA | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E |
| F20% | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E |
| F30% | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E | −−− | E |
| F50% | +++ | P | ++ | P | +++ | P | +++ | P | −− | E | +++ | P | −−− | E |
![]() |
||||||||||||||
| Distribution | ||||||||||||||
| PPB (%) | 99.60 | P | 78.00 | E | 99.60 | P | 79.50 | E | 99.50 | P | 76.5 | E | 94.20 | P |
| VDss (L kg−1) | 3.546 | E | 6.224 | E | 4.144 | E | 5.6750 | E | 3.3880 | E | 2.6600 | E | 2.576 | E |
| BBB penetration | −−− | E | − | M | − | M | ++ | P | + | M | −−− | E | ++ | P |
| Fu (%) | 0.30 | P | 20.50 | E | 0.30 | P | 16.40 | E | 0.40 | P | 18.0 | E | 6.30 | E |
| OATP1B1 inhibitor | +++ | P | +++ | P | +++ | P | +++ | P | +++ | P | +++ | P | ++ | P |
| OATP1B3 inhibitor | +++ | P | +++ | P | +++ | P | +++ | P | +++ | P | +++ | P | +++ | P |
| BCRP inhibitor | +++ | P | −−− | P | ++ | P | −−− | P | + | M | −−− | P | −−− | P |
| MRP1 inhibitor | ++ | P | ++ | P | − | M | − | M | +++ | P | − | M | +++ | P |
| BSEP inhibitor | +++ | P | + | P | +++ | P | +++ | P | +++ | P | +++ | P | − | P |
![]() |
||||||||||||||
| Metabolism | ||||||||||||||
| CYP1A2 inhibitor | +++ | +++ | +++ | +++ | +++ | +++ | +++ | |||||||
| CYP1A2 substrate | +++ | +++ | +++ | +++ | +++ | −−− | −−− | |||||||
| CYP2C19 inhibitor | −− | −−− | ++ | − | − | −−− | −−− | |||||||
| CYP2C19 substrate | −−− | +++ | −− | +++ | +++ | +++ | +++ | |||||||
| CYP2C9 inhibitor | +++ | + | +++ | +++ | +++ | −−− | ++ | |||||||
| CYP2C9 substrate | −−− | −−− | −−− | −−− | −−− | −−− | +++ | |||||||
| CYP2D6 inhibitor | ++ | −−− | +++ | −− | − | −−− | −−− | |||||||
| CYP2D6 substrate | −−− | −−− | −−− | −−− | −−− | −−− | −−− | |||||||
| CYP3A4 inhibitor | −− | +++ | + | +++ | − | −−− | −−− | |||||||
| CYP3A4 substrate | −−− | +++ | +++ | +++ | +++ | ++ | −−− | |||||||
| CYP2B6 inhibitor | +++ | −−− | +++ | −−− | +++ | +++ | −−− | |||||||
| CYP2B6 substrate | −−− | −−− | −−− | −−− | −−− | −−− | −−− | |||||||
| CYP2C8 inhibitor | +++ | ++ | +++ | −− | −−− | −− | −−− | |||||||
| HLM stability | − | M | −− | P | + | M | −−− | P | ++ | E | −−− | P | −−− | P |
![]() |
||||||||||||||
| Excretion | ||||||||||||||
| CLplasma (mL min−1 kg−1) | 4.52 | E | 4.222 | E | 4.72 | E | 4.293 | E | 5.437 | M | 5.329 | M | 5.41 | M |
| T1/2 | 0.609 | 0.47 | 0.82 | 0.65 | 0.682 | 0.297 | 0.926 | |||||||
![]() |
||||||||||||||
| Toxicity | ||||||||||||||
| hERG blockers | 0.663 | M | 0.953 | P | 0.728 | P | 0.965 | P | 0.639 | M | 0.949 | P | 0.272 | E |
| hERG blockers (10 μm) | 0.846 | P | 0.868 | P | 0.872 | P | 0.890 | P | 0.856 | P | 0.337 | M | 0.359 | M |
| DILI | 0.772 | P | 0.990 | P | 0.845 | P | 0.994 | P | 0.832 | P | 0.998 | P | 0.985 | P |
| AMES toxicity | 0.167 | E | 0.264 | E | 0.145 | E | 0.234 | E | 0.125 | E | 0.545 | M | 0.503 | M |
| Rat oral acute toxicity | 0.202 | E | 0.720 | P | 0.237 | E | 0.756 | P | 0.263 | E | 0.289 | E | 0.428 | M |
| FDAMDD | 0.521 | M | 0.330 | P | 0.541 | M | 0.350 | M | 0.473 | M | 0.457 | M | 0.212 | E |
| Skin sensitization | 0.420 | M | 0.922 | P | 0.462 | M | 0.932 | P | 0.275 | E | 0.635 | M | 0.286 | E |
| Carcinogenicity | 0.652 | M | 0.149 | E | 0.670 | M | 0.162 | E | 0.709 | P | 0.947 | P | 0.955 | P |
| Eye corrosion | 0 | E | 0.088 | E | 0.001 | E | 0.124 | E | 0.001 | E | 0.0 | E | 0 | E |
| Eye irritation | 0.57 | M | 0.743 | P | 0.416 | M | 0.608 | M | 0.343 | M | 0.0 | E | 0.423 | M |
| Respiratory toxicity | 0.484 | M | 0.998 | P | 0.463 | M | 0.998 | P | 0.487 | M | 0.791 | P | 0.733 | P |
| Human hepatotoxicity | 0.782 | P | 0.987 | P | 0.766 | P | 0.986 | P | 0.805 | P | 0.986 | P | 0.967 | P |
| Drug-induced nephrotoxicity | 0.738 | P | 0.999 | P | 0.815 | P | 1 | P | 0.769 | P | 0.997 | P | 0.982 | P |
| Drug-induced neurotoxicity | 0.774 | P | 0.993 | P | 0.844 | P | 0.995 | P | 0.856 | P | 0.991 | P | 0.994 | P |
| Ototoxicity | 0.395 | M | 0.842 | P | 0.479 | M | 0.883 | P | 0.507 | M | 0.864 | P | 0.85 | P |
| Hematotoxicity | 0.238 | E | 0.546 | M | 0.319 | M | 0.642 | M | 0.305 | M | 0.475 | M | 0.323 | M |
| Genotoxicity | 0.925 | P | 0.999 | P | 0.828 | P | 0.999 | P | 0.741 | P | 1.0 | P | 1 | P |
| RPMI-8226 immunotoxicity | 0.073 | E | 0.070 | E | 0.080 | E | 0.076 | E | 0.066 | E | 0.624 | M | 0.36 | M |
| A549 cytotoxicity | 0.617 | M | 0.639 | M | 0.722 | P | 0.738 | P | 0.67 | M | 0.051 | E | 0.074 | E |
| Hek293 cytotoxicity | 0.874 | P | 0.391 | M | 0.9 | P | 0.457 | M | 0.864 | P | 0.749 | P | 0.424 | M |
| BCF | 1.490 | 0.789 | 1.817 | 1.266 | 1.889 | 0.755 | 0.632 | |||||||
| IGC50 | 3.841 | 3.287 | 4.119 | 3.696 | 4.005 | 3.379 | 3.241 | |||||||
| LC50DM | 5.113 | 4.651 | 5.232 | 4.810 | 5.413 | 5.31 | 4.997 | |||||||
| LC50FM | 4.722 | 4.030 | 5.012 | 4.473 | 5.046 | 4.206 | 3.899 | |||||||
P < 5), hydrogen bond donors/acceptors, and topological polar surface area (TPSA < 140 Å2), aligning with characteristics of orally bioavailable drugs. However, these compounds fail the Pfizer rule, suggesting potential toxicity risks due to high lipophilicity or reactive groups, unlike Ged and Bim, which pass the Pfizer rule, indicating lower toxicity potential. Compounds 2f and 3f, with piperazine substitutions, also meet Lipinski, Veber, Egan, and Muegge criteria but fail the Ghose rule, which may limit their pharmacokinetic profiles compared to 2c, 3c, and 4c. All five compounds achieve a bioavailability score of 0.55, matching Bim but surpassing Ged (0.17), suggesting comparable oral absorption potential to Bim, driven by favorable PSA and hydrogen bonding properties. The excellent SAscore across all compounds indicates synthetic feasibility, comparable to Ged and Bim, facilitating potential scale-up for clinical development. Overall, 2c, 3c, and 4c exhibited superior druglikeness compared to Ged, with profiles comparable to Bim, though their Pfizer rule violations warrant further optimization to mitigate potential toxicity while maintaining their promising pharmacokinetic properties.
Papp: −4.883, −5.105, −4.942, respectively), comparable to Bim (−4.546) and superior to Ged (−5.322) and 2f (−5.283), which are rated poor, suggesting enhanced intestinal epithelial transport, likely due to favorable lipophilicity and molecular size. All compounds, including Ged and Bim, demonstrated excellent MDCK permeability, indicating efficient transcellular diffusion across renal epithelial cells. PAMPA results showed excellent permeability for 2c, 2f, 3c, 4c, Ged, and Bim, but 3f is rated medium, possibly due to piperazine-induced polarity reducing membrane penetration. Regarding Pgp interactions, all compounds except 2f and 3f are poor Pgp inhibitors, similar to Ged and Bim, reducing the risk of drug–drug interactions, but 2f and 3f are excellent, suggesting potential inhibition that could enhance bioavailability of co-administered drugs. However, 2f and 3f are poor Pgp substrates, indicating efflux susceptibility, which may reduce their absorption compared to 2c, 3c, 4c, and Bim (excellent substrates), while Ged is a poor substrate, potentially limiting its intestinal absorption. All compounds exhibited excellent HIA, F20%, and F30%, reflecting high intestinal absorption and bioavailability at lower doses, consistent with their favorable physicochemical profiles. However, 2c, 2f, 3c, 3f, and Ged show poor F50% bioavailability, indicating reduced absorption at higher doses, possibly due to saturation of transport mechanisms, whereas 4c and Bim achieved excellent F50%, suggesting robust dose-dependent absorption. Overall, 2c, 3c, and 4c demonstrated superior absorption profiles, closely matching Bim and surpassing Ged, particularly in Caco-2 and F50% metrics, though 2f and 3f are limited by poor substrate profiles and Caco-2 permeability, warranting structural optimization to enhance absorption.| Entry | Compound | EGFR | VEGFR2 | PI3K | ||||||
|---|---|---|---|---|---|---|---|---|---|---|
| BA | Bond type (Å) | AA | BA | Bond type (Å) | AA | BA | Bond type (Å) | AA | ||
| a BA – binding affinity (kcal mol−1), bond type (distance/bond length – Å), AA – amino acid, HB–hydrogen bond (conventional/strong hydrogen bond), CHB – carbon–hydrogen bond, DHB – π-donor hydrogen bond, EGFR – epidermal growth factor receptor, VEGFR2 – vascular endothelial growth factor receptor 2, PI3K – phosphoinositide 3-kinase. | ||||||||||
| 1 | 2c | −7.4 | 1 π-sulfur (3.52), 11 hydrophobic (3.66–5.48) | LEU844, MET790, LEU718, CYS797, LYS745, LEU788, ALA743, VAL726 | −8.8 | 13 hydrophobic (3.72–5.50) | VAL846, VAL914, LEU838, ALA864, LEU1033, LEU887, VAL897, PHE1045, LYS866 | −8.8 | 1 HB (3.58), 10 hydrophobic (3.60–5.50) | ASP964, ILE963, TYR867, ILE879, ALA885, ILE881, VAL882, MET953, PHE961 |
| 2 | 2f | −7.8 | 2 HB (2.27, 3.48), 4 hydrophobic (3.56–4.84) | GLU762, THR854, LEU718, LEU844, ALA743 | −8.9 | 1 HB (2.10), 8 hydrophobic (4.13–5.45) | CYS917, LYS866, VAL912, VAL914, VAL846, CYS1043, VAL897 | −9.0 | 1 HB (2.78), 8 hydrophobic (3.68–5.16) | ASP836, ILE963, TYR867, ALA885, ILE881, VAL882, MET953 |
| 3 | 3c | −7.6 | 8 Hydrophobic (3.64–5.35) | LEU718, LYS745, LEU788, MET790, LEU844, ALA743 | −9.1 | 2 HB (3.53, 3.38), 16 hydrophobic (3.49–5.48) | LEU838, VAL846, VAL914, ARG840, LYS866, VAL912, LEU887, VAL897, LEU1033, ALA864 | −8.7 | 8 Hydrophobic (3.79–5.38) | ILE963, ILE879, VAL882, MET953, TYR867, ILE831 |
| 4 | 3f | −7.8 | 1 HB (2.85), 6 hydrophobic (3.66–5.17) | GLU762, LEU718, LEU792, LEU844, ALA743 | −8.9 | 1 HB (3.54), 11 hydrophobic (3.89–5.47) | ASP1044, VAL846, VAL914, LEU887, VAL897, VAL912, PHE1045, CYS1043 | −8.8 | 2 HB (2.96, 2.91), 1 electrostatic (4.19), 1 DHB (3.06), 5 hydrophobic (3.59–4.34) | SER806, ASP841, ASP964, LYS807, ILE968 |
| 5 | 4c | −7.7 | 12 Hydrophobic (3.63–5.48) | LEU718, ALA743, LYS745, LEU788, MET790, CYS797, LEU844, MET793, VAL726 | −8.6 | 2 HB (3.60, 3.33), 15 hydrophobic (3.76–5.43) | LEU838, VAL846, VAL914, ARG840, LEU1033, LYS866, LEU887, PHE1045, VAL897, CYS1043 | −8.6 | 7 Hydrophobic (3.65–5.44) | ILE963, ILE879, VAL882, MET953, TYR867, ILE831 |
| 6 | Gefitinib | −7.3 | 1 HB (3.54), 1 π-sulfur (3.56), 8 hydrophobic (3.83–5.48) | ASP855, MET790, LEU718, VAL726, ALA743, LEU844, LYS745 | — | — | — | — | — | — |
| 7 | Pazopanib | — | — | — | −10.0 | 1 HB (3.69), 9 hydrophobic (3.64–5.43) | GLU883, LEU838, LEU1033, LYS866, VAL914, ALA864, VAL846, VAL897, CYS1043 | — | — | — |
| 8 | Bimiralisib | — | — | — | — | — | — | −9.1 | 2 HB (2.05, 3.28), 1 electrostatic (4.15), 11 hydrophobic (3.59–5.47) | VAL882, ASP836, ASP964, ILE879, ILE963, ILE881, ALA885, MET953, PRO810, ILE831, LYS833 |
![]() | ||
| Fig. 4 2D and 3D representation of the interaction of potent chlorophenylamino-s-triazine derivatives and reference drug Gefitinib with epidermal growth factor receptor (EGFR) target. | ||
![]() | ||
| Fig. 5 2D and 3D representation of the interaction of potent chlorophenylamino-s-triazine derivatives and reference drug Pazopanib with vascular endothelial growth factor receptor 2 (VEGFR2) target. | ||
![]() | ||
| Fig. 6 2D and 3D representation of the interaction of potent chlorophenylamino-s-triazine derivatives and reference drug Bimiralisib with phosphoinositide 3-kinase (PI3K) target. | ||
The binding affinities of the compounds (2c, 2f, 3c, 3f, and 4c) to EGFR, VEGFR2, and PI3K were evaluated against the reference drugs Gefitinib, Pazopanib, and Bimiralisib, respectively, to assess their inhibitory potential. For EGFR, five compounds exhibited binding affinities ranging from −7.4 to −7.8 kcal mol−1, surpassing Gefitinib's −7.3 kcal mol−1. Compounds 2f and 3f showed the highest affinity at −7.8 kcal mol−1, followed by 4c (−7.7 kcal mol−1), 3c (−7.6 kcal mol−1), and 2c (−7.4 kcal mol−1), indicating a slight but consistent improvement over Gefitinib. For VEGFR2, five compounds showed the binding affinities range from −8.6 to −9.1 kcal mol−1, approaching but not exceeding Pazopanib's −10.0 kcal mol−1. Compound 3c showed the highest affinity (−9.1 kcal mol−1), followed by 2f and 3f (−8.9 kcal mol−1), 2c (−8.8 kcal mol−1), and 4c (−8.6 kcal mol−1), suggesting robust but slightly weaker binding compared to Pazopanib. For PI3K, potent compounds exhibited binding affinities range from −8.6 to −9.0 kcal mol−1, closely matching Bimiralisib's −9.1 kcal mol−1. Compound 2f (−9.0 kcal mol−1) nearly equals Bimiralisib, followed by 2c and 3f (−8.8 kcal mol−1), and 3c and 4c (−8.6 kcal mol−1), indicating high potency across all compounds. These results highlighted 2f and 3c as consistently strong binders across all three targets, closely rivaling the reference drugs.
At the EGFR active site, Gefitinib binds to EGFR with an affinity of −7.3 kcal mol−1, forming one hydrogen bond (HB) with ASP855 (3.54 Å), one π-sulfur interaction with MET790 (3.56 Å), and eight hydrophobic interactions with LEU718, VAL726, ALA743, LEU844, and LYS745 (3.83–5.48 Å). Compound 2c (−7.4 kcal mol−1) mirrored Gefitinib's π-sulfur interaction with MET790 (3.52 Å) and shared hydrophobic interactions with LEU718, ALA743, LEU844, LYS745, and VAL726, supplemented by additional hydrophobic contacts with CYS797 and LEU788, enhancing its binding network. Compound 2f (−7.8 kcal mol−1) formed two HBs with GLU762 (2.27 Å) and THR854 (3.48 Å), which are absent in Gefitinib, and shared hydrophobic interactions with LEU718, LEU844, and ALA743, suggesting a more diverse interaction profile. Compound 3c (−7.6 kcal mol−1) relies on eight hydrophobic interactions, overlapping with Gefitinib at LEU718, LYS745, MET790, LEU844, and ALA743, but lacks polar interactions, potentially limiting its specificity. Compound 3f (−7.8 kcal mol−1) formed one HB with GLU762 (2.85 Å), a unique feature, and shared hydrophobic interactions with LEU718, LEU844, and ALA743. Compound 4c (−7.7 kcal mol−1) formed 12 hydrophobic interactions, including LEU718, ALA743, LYS745, MET790, VAL726, and LEU844, closely resembling Gefitinib's hydrophobic profile but lacking polar interactions. The shared interactions with MET790, LEU718, LEU844, and ALA743 across all compounds indicate a conserved EGFR binding pocket, with 2f and 3f introducing unique polar contacts.
At the VEGFR2 active site, Pazopanib binded to VEGFR2 with an affinity of −10.0 kcal mol−1, forming one HB with GLU883 (3.69 Å) and nine hydrophobic interactions with LEU838, LEU1033, LYS866, VAL914, ALA864, VAL846, VAL897, and CYS1043 (3.64–5.43 Å). Compound 2c (−8.8 kcal mol−1) formed 13 hydrophobic interactions, sharing VAL846, VAL914, LEU838, LEU1033, LYS866, VAL897, and CYS1043 with Pazopanib, but lacks polar interactions, which may explain its lower affinity. Compound 2f (−8.9 kcal mol−1) formed one HB with CYS917 (2.10 Å), absent in Pazopanib, and shared hydrophobic interactions with VAL846, VAL914, LYS866, VAL897, and CYS1043, suggesting a partially conserved binding mode. Compound 3c (−9.1 kcal mol−1) formed two HBs with LEU838 (3.53 Å) and VAL846 (3.38 Å) and 16 hydrophobic interactions, including VAL846, VAL914, LEU838, LYS866, VAL897, LEU1033, and CYS1043, closely resembling Pazopanib's profile while adding polar interactions for enhanced stability. Compound 3f (−8.9 kcal mol−1) formed one HB with ASP1044 (3.54 Å) and 11 hydrophobic interactions, sharing VAL846, VAL914, LEU887, VAL897, and CYS1043 with Pazopanib. Compound 4c (−8.6 kcal mol−1) formed two HBs with LEU838 (3.60 Å) and VAL846 (3.33 Å) and 15 hydrophobic interactions, overlapping with Pazopanib at VAL846, VAL914, LEU838, LEU1033, LYS866, VAL897, and CYS1043. The shared hydrophobic interactions with VAL846, VAL914, LYS866, and VAL897 across all compounds indicate a conserved binding site, with 3c and 4c introducing additional HBs to enhance binding.
At the PI3K active site, Bimiralisib binded to PI3K with an affinity of −9.1 kcal mol−1, forming two HBs with VAL882 (2.05 Å) and ASP836 (3.28 Å), one electrostatic interaction with ASP964 (4.15 Å), and 11 hydrophobic interactions with VAL882, ILE879, ILE963, ILE881, ALA885, MET953, and ILE831 (3.59–5.47 Å). Compound 2c (−8.8 kcal mol−1) formed one HB with ASP964 (3.58 Å) and 10 hydrophobic interactions, sharing VAL882, ILE879, ILE963, ALA885, and MET953 with Bimiralisib, indicating a similar binding mode. Compound 2f (−9.0 kcal mol−1) formed one HB with ASP836 (2.78 Å), mirroring Bimiralisib's interaction, and shared hydrophobic interactions with VAL882, ILE963, ALA885, and MET953, suggesting high similarity. Compound 3c (−8.7 kcal mol−1) showed eight hydrophobic interactions, overlapping with Bimiralisib at VAL882, ILE879, ILE963, MET953, and ILE831, but lacks polar interactions. Compound 3f (−8.8 kcal mol−1) formed two HBs with SER806 (2.96 Å) and ASP841 (2.91 Å), one electrostatic interaction with ASP964 (4.19 Å), one π-donor HB with LYS807 (3.06 Å), and five hydrophobic interactions, sharing ASP964 with Bimiralisib while introducing unique polar contacts. Compound 4c (−8.6 kcal mol−1) formed seven hydrophobic interactions, overlapping with Bimiralisib at VAL882, ILE879, ILE963, MET953, and ILE831, but lacks polar interactions. The shared hydrophobic interactions with VAL882, ILE963, and MET953 across all compounds highlight a conserved PI3K binding pocket, with 2f and 3f closely mimicking Bimiralisib's polar interactions.
In summary, compounds 2f and 3c stand out as the most promising candidates for multi-targeted inhibition of EGFR, VEGFR2, and PI3K, based on their robust binding affinities and interaction profiles, which closely rival or approach those of the reference drugs Gefitinib, Pazopanib, and Bimiralisib. For EGFR, 2f showed the highest affinity (−7.8 kcal mol−1), surpassing Gefitinib (−7.3 kcal mol−1), with unique hydrogen bonds to GLU762 and THR854, alongside hydrophobic interactions with LEU718, LEU844, and ALA743, suggesting enhanced binding specificity. Compound 3c (−7.6 kcal mol−1) also outperformed Gefitinib, sharing key hydrophobic interactions with MET790, LEU718, and LEU844. For VEGFR2, 3c exhibited the highest affinity (−9.1 kcal mol−1), closely approaching Pazopanib (−10.0 kcal mol−1), with two hydrogen bonds (LEU838, VAL846) and an extensive hydrophobic network (VAL846, VAL914, LYS866, VAL897), indicating strong binding stability. Compound 2f (−8.9 kcal mol−1) complemented this with a hydrogen bond to CYS917 and shared hydrophobic interactions, reinforcing its potency. For PI3K, 2f (−9.0 kcal mol−1) nearly matches Bimiralisib (−9.1 kcal mol−1), with a hydrogen bond to ASP836 and hydrophobic interactions with VAL882, ILE963, and MET953, closely mimicking Bimiralisib's profile. Compound 3f (−8.8 kcal mol−1) also showed promise with diverse polar interactions (SER806, ASP841, ASP964, LYS807), though 2f's closer affinity to Bimiralisib gives it an edge. While 2c, 3f, and 4c demonstrated competitive affinities and interactions, 2f and 3c consistently excel across all targets due to their high affinities and balanced polar and hydrophobic interactions, making them prime candidates for further optimization as multi-targeted kinase inhibitors.
Microwave-assisted reactions were performed in a CEM Discover Microwave Synthesizer (USA), equipped with a magnetic stirrer for thorough mixing and an infrared sensor for precise temperature monitoring and control. Melting points were determined using a Sanyo–Gallenkamp apparatus, providing accurate thermal characterization. Nuclear magnetic resonance (NMR) spectra were obtained on a Bruker Avance 600/500 spectrometer, with chemical shifts reported in parts per million (δ, ppm) relative to standard references. High-resolution mass spectrometry was performed on an Agilent Series 1100 LC-MS trap system, ensuring accurate molecular weight determination. Finally, the optical density (A) was measured at 570 nm on a MultiskanTM microplate reader in the anti-cancer activity assay.
:
1 (v/v) ethanol–water mixture, yielding the desired mono-substituted s-triazine derivatives in excellent yields of 95–97%.
:
water mixture (2
:
8, v/v) or by column chromatography on silica gel with a hexane
:
ethyl acetate eluent system. Reaction yields range from 78 to 86%.
:
8, v/v) or column chromatography on silica gel using a hexane–ethyl acetate eluent. Reaction yields range from 88 to 95%.
), 0.89 (6H, d, J = 5.0 Hz, –CH3). 13C NMR (125 MHz, DMSO-d6, δ ppm): 164.3, 164.0, 142.2, 132.7, 129.8, 120.4, 118.6, 117.3, 43.0, 33.6, 30.6, 21.7. LC-MS (m/z) [M + H]+ calcd for C21H30ClN6 401.2215, found 401.2197.
), 1.01 (4H, q, J = 10.5 Hz, –CH2–), 0.90 (6H, d, J = 5.5 Hz, –CH3). 13C NMR (125 MHz, DMSO-d6, δ ppm): 164.3, 163.9, 140.9, 130.5, 130.0, 122.1, 120.2, 119.0, 43.0, 33.5, 30.6, 21.7. LC-MS (m/z) [M + H]+ calcd for C21H29Cl2N6 435.1825, found 435.1827.
), 0.98 (4H, q, J = 8.0 Hz, –CH2–), 0.90 (6H, d, J = 5.0 Hz, –CH3). 13C NMR (125 MHz, DMSO-d6, δ ppm): 164.3, 164.1, 135.5, 128.5, 127.2, 127.1, 126.4, 126.0, 42.9, 33.6, 30.6, 21.7. LC-MS (m/z) [M + H]+ calcd for C21H29Cl2N6 435.1825, found 435.1825.
−
H]– calcd for C17H18Cl2N6O2 409.0952, found 409.0947.
The IC50 values, representing the compound concentration required to inhibit 50% of cell proliferation, were determined by plotting the inhibition percentages against the logarithm of the compound concentrations and fitting the data to a dose–response curve using GraphPad Prism 10.
| Entry | Target | Symbol | PDB ID | Organism | Expression system | Ref drug |
|---|---|---|---|---|---|---|
| 1 | Epidermal growth factor receptor | EGFR | 3UG2 | Homo sapiens | Spodoptera frugiperda | Gefitinib |
| 2 | Vascular endothelial growth factor receptor 2 | VEGFR2 | 3CJG | Homo sapiens | Spodoptera frugiperda | Pazopanib |
| 3 | Phosphoinositide 3-kinase | PI3K | 5OQ4 | Homo sapiens | Spodoptera frugiperda | Bimiralisib |
| Target | Size | Center | ||||
|---|---|---|---|---|---|---|
| x | y | z | x | y | z | |
| EGFR | 25 | 25 | 25 | −0.1842 | 49.3505 | 20.0221 |
| VEGFR2 | 25 | 25 | 25 | 8.1773 | 40.7893 | 7.3664 |
| PI3K | 25 | 25 | 25 | 43.1501 | 15.2623 | 32.1771 |
| This journal is © The Royal Society of Chemistry 2025 |