Next Article in Journal
The Therapeutic Potential of Harpagophytum procumbens and Turnera subulata and Advances in Nutraceutical Delivery Systems in Neurodegenerative Diseases
Previous Article in Journal
Effects of Smallanthus sonchifolius Flour on Metabolic Parameters: A Systematic Review
Previous Article in Special Issue
Structural Investigations on 2-Amidobenzimidazole Derivatives as New Inhibitors of Protein Kinase CK1 Delta
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovery and Anticancer Screening of Novel Oxindole-Based Derivative Bearing Pyridyl Group as Potent and Selective Dual FLT3/CDK2 Kinase Inhibitor

1
Department of Medicinal Chemistry, Faculty of Pharmacy, Minia University, Minia 61519, Egypt
2
Biotechnology Institute, Ankara University, Ankara 06135, Turkey
3
Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Erciyes University, Kayseri 38280, Turkey
4
Department of Medicinal Chemistry, Institute of Pharmacy, Martin-Luther-University of Halle-Wittenberg, 06120 Halle, Germany
5
Pharmacology and Experimental Oncology Unit, Department of Cancer Biology, National Cancer Institute, Cairo University, Cairo 11796, Egypt
6
Department of Pharmaceutical Sciences, College of Pharmacy, AlMaarefa University, Riyadh 13713, Saudi Arabia
7
Department of Pharmaceutical Sciences, College of Pharmacy, Princess Nourah bint Abdulrahman University, Riyadh 11671, Saudi Arabia
8
Department of Medicinal Chemistry, Faculty of Pharmacy, Port Said University, Port Said 42511, Egypt
*
Authors to whom correspondence should be addressed.
Pharmaceuticals 2024, 17(5), 659; https://doi.org/10.3390/ph17050659
Submission received: 23 April 2024 / Revised: 14 May 2024 / Accepted: 15 May 2024 / Published: 20 May 2024
(This article belongs to the Special Issue Kinase Inhibitors in Targeted Cancer Therapy)

Abstract

:
Protein kinases regulate cellular activities and make up over 60% of oncoproteins and proto-oncoproteins. Among these kinases, FLT3 is a member of class III receptor tyrosine kinase family which is abundantly expressed in individuals with acute leukemia. Our previous oxindole-based hit has a particular affinity toward FLT3 (IC50 = 2.49 μM) and has demonstrated selectivity towards FLT3 ITD-mutated MV4-11 AML cells, with an IC50 of 4.3 μM. By utilizing the scaffold of the previous hit, sixteen new compounds were synthesized and screened against NCI-60 human cancer cell lines. This leads to the discovery of a potent antiproliferative compound, namely 5l, with an average GI50 value against leukemia and colon cancer subpanels equalling 3.39 and 5.97 µM, respectively. Screening against a specific set of 10 kinases that are associated with carcinogenesis indicates that compound 5l has a potent FLT3 inhibition (IC50 = 36.21 ± 1.07 nM). Remarkably, compound 5l was three times more effective as a CDK2 inhibitor (IC50 = 8.17 ± 0.32 nM) compared to sunitinib (IC50 = 27.90 ± 1.80 nM). Compound 5l was further analyzed by means of docking and molecular dynamics simulation for CDK2 and FLT3 active sites which provided a rational for the observed strong inhibition of kinases. These results suggest a novel structural scaffold candidate that simultaneously inhibits CDK2 and FLT3 and gives encouragement for further development as a potential therapeutic for leukemia and colon cancer.

Graphical Abstract

1. Introduction

Cancer is an extremely diverse disease, encompassing over 200 distinct types, and its origins are quite intricate. Globally, cancer is widely regarded as the most prevalent and concerning disease that impacts human health. It is referred to as a collection of diseases that are characterized by the uncontrolled proliferation and division of cells. Cancer has hallmarks such as exploiting growth signals, evading suppressive mechanisms, resisting cell death, inducing migration and angiogenesis. In addition, emerging phenotypes of cancer continue to be defined especially related to cell metabolism [1,2,3]. Cancer diagnosis and treatment protocols are based on the elements of these features. In such a complex disease, efforts are focused on the development of functional agents such as small molecule combinations, immunochemotherapeutics, dual inhibitors, multi-target peptides and smart drugs [4,5,6]. Thus, pharmacologically active sources, novel designs, repurposing, biological inhibitors, various manipulators and artificial intelligence are used for multifunctional therapeutic formulations and rationalizations [7,8,9,10,11,12].
Protein kinase enzymes play a crucial part in the signal transduction pathways that regulate numerous cellular functions. The majority of cancer types are caused by protein kinases malfunction, with protein kinases being responsible for more than 60% of all oncoproteins and proto-oncoproteins (the main contributors to the cancer condition) [13]. One of these protein kinases is the Fms-like receptor tyrosine kinase 3 (FLT3), which belongs to the class III receptor tyrosine kinase (RTK) group. FLT3 plays a crucial role in regulating hematopoiesis and is commonly found to be overexpressed in the majority of acute leukemia patients. The FLT3 pathway is pivotal for the proliferation and specialization of hematopoietic progenitor cells (HPCs). Multiple studies have demonstrated a strong association between AML and other hematologic malignancies and aberrant FLT3 pathways. Similar to other receptor tyrosine kinases (RTKs), the FLT3 receptor forms dimers upon binding to the FLT3 ligand. This leads to autophosphorylation and the subsequent activation of downstream signaling cascades, including RAS/MEK, PI3K/AKT/mTOR, and JAK/STAT [14]. These pathways have crucial functions in controlling the cell cycle, cell death, and cell specialization [15]. Cyclin-dependent kinases (CDKs) are another type of protein kinases that specifically phosphorylate serine and threonine residues. They play a vital role in regulating the course of the cell cycle and promoting cell proliferation. Blocking the activity of CDKs has been discovered to restrict the uncontrolled cellular proliferation observed in certain forms of malignancies [16]. CDK2 plays a critical role in promoting cancer progression when it is overexpressed. Therefore, suppressing the excessive activity of CDK2 may lead to the reversal of the malignant characteristics of tumor cells. CDK inhibitors, including AT-7519, flavopiridol, AMG-925, and palbociclib, are now undergoing clinical development for the treatment of several types of solid tumors and hematological malignancies. Nevertheless, the majority of CDK inhibitors were terminated during phase II trials due to significant toxicities or restricted clinical efficacy as standalone treatments.
Combination therapies that target multiple signaling pathways can potentially overcome resistance and prolong the short-term response to single FLT3 inhibitors. Although certain combination therapies have exhibited encouraging clinical outcomes, their development may be impeded by unforeseeable drug–drug interactions, pharmacokinetic characteristics, and cumulative toxicity. Multi-target drugs can address the limitations of combination therapies by simultaneously inhibiting multiple pathways. Multiple preclinical multi-target inhibitors display potential in overcoming resistance to FLT3 inhibitors [17].
Over the past few decades, there has been a growing focus on kinase enzyme inhibitors due to their current application in the treatment of cancer. The intricate nature of tumor formation, which encompasses numerous molecular pathways, justifies the use of multi-targeted therapy for cancer. Agents with the ability to suppress various molecular targets, affecting different cell types (such as endothelial cells, pericytes, fibroblasts, and tumor cells), are theoretically more likely to achieve success. Oxindoles, which serve as a fundamental component in the multikinase inhibitor I, sunitinib (Figure 1), have also been employed in the development of several kinase inhibitors. Furthermore, compound II, FN1501, is currently being tested in phase I/II clinical trials (NCT03690154) as a multiple-kinase inhibitor for advanced solid tumors [18]. It has demonstrated strong inhibitory effects against FLT3 (IC50: 0.27 nM), CDK2 (IC50: 2.47 nM), CDK4 (IC50: 0.85 nM), and CDK6 (IC50: 1.96 nM) [19]. Additionally, it has shown remarkable effectiveness in treating leukemia in mouse models with MV4-11 cell xenografts. Our earlier research involved examining a range of synthetic oxindole-based drugs that specifically target FLT3 and had favorable outcomes in treating AML. Compound III, which is based on oxindole, demonstrated remarkable effectiveness against FLT3 kinase, with IC50 values of 2.49 μM. Furthermore, it selectively targeted FLT3 ITD-mutated MV4-11 AML cells, with an IC50 of 4.3 μM after 72 h [20,21]. In addition, we have documented a sequence of substances exhibiting encouraging anticancer properties specifically against the MCF-7 breast cancer cell line. Out of these compounds, IV (Figure 1) demonstrated selectivity towards the MCF-7 breast cancer cell line (IC50 = 17.01 μM) while having no impact on the MCF-12A normal breast cell line [22].
In this current study, we further investigated the anticancer potential of oxindole-based pyridyl derivatives of compound IV. Instead of employing vanillin as a linker, we employed 4-hydroxybenzaldehyde. In addition, different pyridyl groups were attached to this linker to act as an allosteric binding moiety. The objective of our research is to discover novel lead compounds that exhibit promising antiproliferative activity.

2. Results and Discussion

2.1. Chemistry

The preparation methodologies adopted to synthesize the target compounds 5a–p are outlined in Scheme 1. The structures of the final compounds were supported by various spectral and elemental analyses. The final compounds were obtained as a mixture of E and Z isomers, and the spectral data were recorded for the major isomer.
The outline for synthesizing the target compounds 5a–p is illustrated in Scheme 1. The alkylated p-hydroxy benzaldehyde derivatives 3a–d were obtained via an SN2 substitution reaction of various alkyl halides 1a–d with p-hydroxy benzaldehyde’s nucleophilic hydroxyl group. The alkylated p-hydroxy benzaldehyde derivatives 3a–d underwent a base-catalyzed Knoevenagel condensation reaction with the oxindole derivatives 4a–d to afford the compounds 5a–p in a good yield ranging (55–91%). The reported compounds 3a–d were confirmed through measuring the melting point and were compared to the reported result while the newly prepared compounds 5a–p were characterized using various spectroscopic methods. NMR data supported the success of the used route.
In the 1H NMR of all the target compounds, 5a–p showed an almost characteristic singlet signal at δ = 7.6–7.9 ppm for the ethylenic protons; furthermore, the NH protons were noticed as a singlet signal at δ = 9.83–10.74 ppm. An additional singlet signal was observed at δ = 5.1–5.2 ppm for the CH2 protons.
Similarly, the 13C NMR assured the formation of the target scaffolds, showing characteristic ethylenic carbons at δ = 139–140 ppm, the aliphatic-CH2 carbon at δ = 67.20–71.00 ppm and carbonyl carbons of CONH groups at δ = 160.30–169.40 ppm. All the aromatic protons and carbons belonging to the oxindole, phenyl, and benzyl scaffolds were identical with the proposed structures of the compounds.
All the synthesized compounds of 5a–p were obtained as the mixtures of E and Z-diastereomers in different ratios as identified with vinyl -H and 2′ and 6′ Hs chemical shifts (ppm) of E and Z (Table 1). The separation of these diastereomers is not possible, as reported in many previous works that acknowledged that the dynamic interconversion between the two isomers is dependent on solvent, time, and light [23,24,25,26].
1H NMR spectroscopy was employed to validate the assignment of the E/Z-diastereomers; in this analysis, the vinylic protons demonstrated a marginally greater chemical shift in the Z-diastereomer in comparison to the E-diastereomer. The observed variation in chemical shift can be explained by the impact of the carbonyl group situated at position 2 on the oxindole ring. Additionally, it was observed that the E-diastereomers exhibited a more pronounced upfield shift for the 2′ and 6′ ortho-benzylidene protons for the same underlying reason. The 1H NMR spectrum of compound 5l revealed the chemical shifts associated with the principal Z-diastereomer. The signal was 7.80 ppm for vinylic H, 7.44 ppm for the H4 of the oxindole ring, and 8.47 ppm was the signal of H2′,6′. The corresponding signals of the minor E-diastereomer exhibited resonances at 7.63, 7.02, and 7.62 ppm. The ratio of E- to Z-diastereomers for each molecule was ascertained by utilizing the total integration value of the vinylic hydrogen present in both isomers as a reference. In general, it is postulated that a vinylic-H integration of 1 in both E- and Z-diastereomers signifies the complete presence of both E- and Z-diastereomers. Hence, when the upfield-shifted vinylic-H-H integration equals 0.4, it signifies that the proportion of the E-diastereomer is 40%. In a similar fashion, a value of 0.6 for the vinylic-H integration that is further downfield-shifted indicates a 60% Z-diastereomer ratio. NMR data and elemental analysis were employed to validate the purity of the target compounds.

2.2. High Throughput Anticancer Screening for Compounds 5a–p

In a comprehensive test, compounds were assayed at a single dose concentration (10−5 M) in the full panel of NCI-60 cancer cell lines. The detailed anti-tumor activities on the growth inhibition of NCI-60 human cancer cell lines at the single dose of 10 µM are shown in Table 2. The tested compounds showed a diverse but strong antiproliferative effect on the evaluated panel of cell lines and most of the compounds exhibited more than a 50% inhibition of tumor growth at a micromolar concentration. As revealed from the results of 16 tested compounds, six compounds with a wide range in growth percentage displayed strong growth inhibitory activity (−37.79 to 152.57%) at a 10 µM concentration (Table 3).
The 3-pyridyl oxindole hybrids demonstrated superior activity compared to 2-pyridyl and 4-pyridyl oxindole hybrids. The hybridization of 6-Cl oxindole with 3-pyridyl moiety in compound 5l showed powerful activity with cytotoxic effects towards 4 cell lines and with cytostatic effects towards 44 cell lines. Moreover, the replacement of 6-Cl oxindole with 5-fluoro oxindole and the hybridization with 3-pyridyl moiety in compound 5j gives a cytotoxic effect towards only 1 cell line (SK-MEL-5) and a cytostatic effect towards 32 cell lines. The hybridization of 5-Cl oxindole with 2-pyridyl moiety in compound 5g showed powerful activity with cytotoxic effects towards 1 cell line (OVCAR-4), with a growth inhibition of 101.81% and with cytostatic effects towards 17 cell lines. Compound 5g indicated better cytostatic activity than the hybridization of 6-Cl oxindole, 5-F oxindole, and oxindole with 2-pyridyl moiety in the series of compounds 5h, 5f and 5e. Also, the replacement of 3-pyridyl moiety with 4-pyridyl moiety in compounds 5m, 5o, 5p, and 5n demonstrated only cytostatic effects towards 14, 28, 21 and 26 cell lines, respectively. On referring to the total number of sensitive cell lines for tested compounds, it has been found that most of the target compounds exhibited a broad spectrum of anti-tumor activity, covering different cancer subpanels. Amongst these sensitive cell lines, the melanoma cell line (SK-MEL-5) was found to be highly sensitive with a negative growth percentage value (lethal effect) for derivatives 5l, 5j, and 5n. Also, ovarian cancer (OVCAR-4) cell lines were noticed to be the most sensitive cell lines for derivatives 5g and 5h. Moreover, MCF7 (breast cancer) also were proved to be the most responsible cells to compound 5o. Moreover, compounds 5a, 5b, 5c and 5d displayed weak antiproliferative activity. Compound 5l has been shown to have interesting growth inhibitory activity in the preliminary single dose screen, showing effectiveness toward numerous cell lines that belong to different tumor subpanels (Figure 2a). Figure 2b shows the radar chart of 5l activity, indicating that 5l has more overall potent activity (over 90% inhibition) in leukemia, colon cancer and melanoma. Figure 2c expresses the variations resulting from the average of activities on cell lines in each type of cancer. A lower variation in a cancer type indicates consistent inhibition across the cell lines in that panel. While the variations were observed to be over 20% in CNS, melanoma, ovarian and renal cancers, it was observed to be quite low in leukemia, colon and melanoma. Using all this data, it was then evaluated for the advanced 5-dose (10−4–10−8 M) testing mode against the full panel for 5l.

2.3. Dose-Dependent Effects of Compound 5l against NCI 60 Cancer Cell Panel

Compound 5l has shown an interesting growth inhibitory activity in the preliminary single dose screen, showing an effectiveness toward numerous cell lines that belong to different tumor subpanels (Figure 2a–c), and the advanced compound 5l was subjected to a 5-dose (10−4–10−8 M) testing mode against the full panel (Figure 3).
As shown in Table 4, compound 5l showed a high antiproliferative activity at a moderate micromolar concentration on almost all evaluated human cancer types with GI50 values of 2.70–77.10 μM. Interestingly, compound 5l exhibited superior micromolar activity against leukemia (MOLT-4, HL-60, CCRF-CEM, K-562, RPMI-8226), colon (HCT-15, SW-620), melanoma (LOX IMVI, MDA-MB-435, SK-MEL-5), renal cancer (CAKI-1 and UO-31), and breast cancer cell lines (MCF7, T-47D and MDA-MB-468), with a GI50 value range between 2.70 and 4.70 μM. In addition, compound 5l demonstrated a high cytostatic action at a low micromolar concentration against ovarian cancer cell lines (IGROV1 and OVCAR-3) with TGI values of 3.32 and 2.26 μM. Compound 5l showed no cytostatic effect against the left-over cancer cells (TGI > 100 μM). Conversely, compound 5l appeared as a non-lethal agent, which has an LC50 greater than 100 μM for most of the examined cancer cells, excluding melanoma cell lines (LOX IMVI, SK-MEL-5, SK-MEL-28) (LC50 = 96.90, 65.60, and 71.50 μM), respectively, (Table 4).
Regarding the sensitivity towards various tumor cell lines, compound 5l exhibited a potent growth inhibitory activity against the whole NCI panel, with a median GI50 against the full panel (MG-MID) of 12.23 μM, and effective sub-panel median GI50 (MG-MID) values between 3.39 and 29.43 μM. The leukemia sub-panel was the most susceptible cancer type to 5l [GI50 (MG-MID) = 3.39 μM] (Table 5). Moreover, the ratio of the MG-MID of the full-panel to its individual sub-panels offers the selectivity index. Notably, 5l showed the best selectivity index (3.61) for the leukemia cell sub-panel (Table 5).

2.4. Kinase Screening of Compound 5l

The percentage of inhibition of the target compound 5l (10 µM) was tested, versus a panel of ten kinases representing various signaling cascades including PTK2B, JAK1, CDK2, FGFR1, IGF1R, VEDFR-2, PDGFRα, PDGFRβ, FLT3, and SRC kinases, according to Z’-LYTE® technology for all kinases, except where the VEGFR-2, VEGFR-2/KDR ELISA kit is used (Figure 4).
Interestingly, compound 5l exhibited remarkable inhibition percentages on CDK2 and FLT3 kinases with values of 87.71 and 92.59%, respectively. On the contrary, it did not show any promising inhibitory percentages versus the rest. Therefore, as a consequence of the highly promising activities of compound 5l against CDK2 and FLT3, the IC50 value of this compound has to be evaluated versus these two mentioned kinases. Compound 5l was titrated for its potential IC50 value against FLT3 (Figure 5a) and CDK2 (Figure 5b), using sunitinib as a reference. The outcomes were reported as IC50 and are represented in Figure 5c, indicating that compound 5l was presented as being a potent CDK2 inhibitor, showing a single digits nanomolar with an IC50 value (IC50 = 8.17 ± 0.32 nM), compared to sunitinib (IC50 of 27.90 ± 1.80 nM) and FN-1501 (IC50 of 0.79 ± 0.08 nM) as the reference. In contrast, compound 5l possessed moderate FLT-3 inhibitory activity with an IC50 value (IC50 = 36.21 ± 1.07 nM), compared to sunitinib (IC50 of 14.90 ± 0.36 nM) and FN-1501 (IC50 of 1.99 ± 0.12 nM) as the reference.

2.5. Molecular Modeling Studies

In the endeavor to expand the arsenal of targeted treatments for colorectal cancer and acute myeloid leukemia (AML), compound 5l emerges as an indole-2-one derivative with a structural analogy to sunitinib, a potent tyrosine kinase inhibitor. The rationale behind selecting the specific targets for molecular docking studies stems from the direct inhibitory action of sunitinib on the kinase enzymes CDK2 and FLT3 and biological kinase screening. These targets are critically involved in the mechanisms driving cancer cell proliferation and survival, making them prime candidates for therapeutic intervention. To inform and enhance the docking strategy for compound 5l, the crystal structures of the protein CDK2 (PDB ID: 3TI1) [27] were utilized, each of which has sunitinib co-crystallized within. These structures provide a foundational understanding of the key interactions and conformations necessary for effective inhibition. For FLT3, the structure (PDB ID: 6JQR) [28] was selected to explore the binding interactions in silico as in our previous study [20], given the absence of a co-crystallized complex with sunitinib. Through these structural analyses, compound 5l is evaluated to ascertain its potential in replicating or surpassing the inhibitory capabilities of sunitinib, thereby aiding in the design of new inhibitors that could offer enhanced therapeutic benefits for patients suffering from colorectal cancer and AML. To validate the docking protocol for compound 5l, co-crystal ligands were re-docked into the active sites of the CDK2 and FLT3 proteins, resulting in RMSD values of 0.76 and 0.79, respectively. These low RMSD values, which were well below the commonly accepted threshold of 2 Å, indicated a high degree of accuracy in the alignment between the docked conformations and the original co-crystal structures, underscoring the effectiveness of the docking process. The Z isomer of compound 5l exhibited a docking interaction energy of −8.593 kcal/mol with CDK2, whereas its E isomer showed a lesser interaction energy of −6.827 kcal/mol. Furthermore, the Z isomer of 5l demonstrated a docking interaction energy of −9.398 kcal/mol with FLT3, compared to −8.922 kcal/mol for the E isomer. Given that sunitinib, as the reference compound, is also active in the Z configuration, these findings provide a rationale for focusing future molecular docking and dynamics studies on the Z isomer of 5l.
In the molecular docking studies of compound 5l, a comprehensive interaction profile has emerged, delineating the compound’s engagement with several key proteins. These interactions suggest that compound 5l could serve as a potent inhibitor, as evidenced by the specific binding modalities elucidated for each target protein.
Starting with CDK2, as highlighted in Figure 6a, the docking profile is characterized initially by the formation of hydrogen bonds with Glu81 (2.44 Å), Leu83 (1.94 Å), Gln85 (2.25 Å), and Lys89 (2.47 Å). These interactions are complemented by hydrophobic contacts with residues Ile10, Val18, Phe80, Phe82, Leu134, and Asp145. Notably, a π-cation interaction with Lys89 is observed, further solidifying the ligand’s favorable orientation within the binding pocket. Finally, with FLT3, as illustrated in Figure 6b, compound 5l establishes hydrogen bonds with Cys694 (1.98 Å and 2.47 Å) and Asp829 (2.57 Å). This bonding network is supported by hydrophobic interactions with Leu616, Val624, Ala642, Leu818, and Asp829. Additionally, the π-stacking with Phe691 contributes significantly to the binding stability, indicating a highly specific interaction scheme.
In the field of computational drug discovery, molecular dynamics (MD) simulations are invaluable for providing a dynamic perspective on the interactions between drug candidates and their target proteins. Such simulations offer detailed insight into the conformational flexibility, stability, and binding efficacy of compounds at an atomic scale. Within this framework, a 100 ns MD simulation was executed to investigate the binding dynamics of compound 5l with two pivotal kinase targets—FLT3 and CDK2—each playing a crucial role in cancer pathophysiology and thus representing significant therapeutic targets. The primary goal of the simulation was to determine the stability of compound 5l within the active sites of these kinases and to delineate the nature of its binding over the simulated timeframe. In addition, to compare compound 51 with sunitinib, the crystal structure of sunitinib with CDK2 and the protein–ligand complexes obtained by docking sunitinib to FLT3 were also compared using MD simulation. Through meticulous analyses of root-mean-square deviation (RMSD), principal component analysis (PCA), and hydrogen bond dynamics, valuable insights into the binding stability and conformational adaptations of compound 5l were obtained, which are vital for its potential efficacy as a kinase inhibitor.
The RMSD findings revealed that the FLT3-5l complex displayed remarkable stability, with RMSD values persistently under 0.3 nm, denoting a consistently tight binding throughout the simulation (Figure 6d). In contrast, the CDK2-5l complex had slightly higher RMSD values, reaching peaks of up to 0.4 nm, suggesting a stable interaction with a degree of flexibility (Figure 6c). Complementary PCA analysis (Figure 6c,d) reflected these observations, with the FLT3-5l complex demonstrating a densely clustered pattern, indicative of a stable interaction with limited conformational change. The CDK2-5l complex exhibited a more scattered PCA profile, implying a broader range of ligand movements. The hydrogen bond analysis (Figure 6c,d) further substantiated the stability of the FLT3-5l interaction, with a relatively constant hydrogen bond count. The CDK2-5l complex, while maintaining a consistent number of hydrogen bonds, displayed greater variability. Moreover, according to sunitinib’s RMSD profile in the MD simulation (Figure S1), compound 5l was found to be more stable than sunitinib with the target proteins, CDK2 and FLT3.
The stability of compound 5l within the binding pockets of CDK2 and FLT3 is dynamically showcased through animation videos derived from the MD trajectories (Videos S1 and S3). These visualizations add a layer of depth to our understanding by displaying the robustness of the hydrogen bonding interactions in real time. They particularly highlight the persistence of hydrogen bonds with the indole-2-one moiety of compound 5l, while also noting the variable nature of the interactions with the pyridine moiety across the simulation timeframe. Additionally, for the comparison of compound 5l with sunitinib, the animation video obtained from the MD simulation of sunitinib performed with CDK2 is also presented in Videos S2 and S4, respectively.
Further to the dynamic analyses, the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) method was employed to calculate the binding free energies of the complexes, providing an estimate of the binding affinity of compound 5l to the kinases. MM/PBSA calculations were performed on 200 frames extracted between 80 and 100 ns of the simulation, resulting in the plot shown in Figure 6e. This plot offers a quantitative view of the energetic components contributing to the binding affinity, such as van der Waals interactions, electrostatic contributions, and solvation effects, thereby complementing the qualitative insights from the MD simulation. Integrating the MM/PBSA free energy calculations into our analysis further refines our understanding of the interaction between compound 5l and the kinase targets. The MM/PBSA approach is instrumental in quantifying the binding affinities through various energetic contributions, offering a complementary perspective to the dynamic observations from MD simulations. According to MM/PBSA data, the binding free energies (ΔGMM/PBSA) are −22.52 ± 2.20 kcal/mol for CDK2-5l, and −27.09 ± 1.27 kcal/mol for FLT3-5l, indicating a hierarchy of binding affinities, with FLT3-5l exhibiting the strongest interaction. These values reveal that the FLT3-5l complex possesses the most favorable binding affinity with a ΔGMM/PBSA of −27.09 kcal/mol, which correlates with the lower RMSD values and the tighter clustering observed in PCA. The standard deviation of ±1.27 kcal/mol also indicates a consistent binding affinity throughout the sampled frames. The CDK2-5l complex, while displaying a ΔGMM/PBSA of −22.52 kcal/mol, shows the largest standard deviation, reflecting a broader range of interaction energies which may be indicative of the flexibility observed in the MD simulation. Overall, the MM/PBSA calculations align well with the stability trends observed in the MD simulations. The FLT3-5l complex stands out with the highest binding affinity and the most stable interaction, making it a promising candidate for further drug development efforts aimed at FLT3 inhibition. The consistency and favorability of these interactions provide a quantitative foundation for the potential efficacy of compound 5l as a therapeutic agent, confirming the qualitative insights gained from the MD simulation and subsequent analysis.

2.6. Physicochemical, ADME, Pharmacokinetic Characteristics and Drug-Likeness Forecast

The Swiss Institute of Bioinformatics (SIB) hosts the freely accessible SwissADME tool, which serves as an invaluable computational resource for globally assessing the pharmacokinetic profiles and drug-likeness of small molecular scaffolds. This analysis is essential for confirming that newly synthesized compounds exhibit promising physiological and pharmacokinetic properties, thereby streamlining the drug development process in terms of cost, time, and efficacy [29].
Compound 5l, a potent antiproliferative agent within the 3-substituted oxindole series, was characterized by its strong hydrophilicity, adequate gastrointestinal (GI) tract absorption, and competent blood–brain barrier (BBB) penetration. These attributes are quantified by a predicted log Po/w value of 3.10. When evaluated against the Lipinski and Veber rules for drug-likeness, compound 5l demonstrated a profile marked by increased lipophilicity and low solubility. Notably, the compound exhibits a degree of unsaturation outside the optimal range, with a fraction of sp3 hybridized carbons (Fraction Csp3) value of 0.05, which is below the preferred threshold of 0.25, indicating potential challenges for its drug-likeness. Despite this, compound 5l is poorly soluble in water and is classified within the category of hERG I inhibitors, which are known to pose cardiotoxic risks, as well as being associated with AMES toxicity and skin sensitization. In contrast, it does not exhibit hepatotoxic or hERG II channel inhibition, which is often related to adverse drug reactions.
The Brain or Intestinal Estimated permeation method (BOILED-Egg) was employed to graphically project and evaluate compound 5l’s passive GI absorption and BBB penetration [30]. As visualized in Figure 7a, compound 5l is aptly situated within the human intestinal absorption window and demonstrates BBB permeability, while also being identified as a non-substrate for P-glycoprotein (PGP-) efflux transporters, which suggests a reduced likelihood of resistance via an active efflux by cancer cells. Figure 7b further highlights the pharmacokinetic predictions where the lipophilicity (WLOGP) and Total Polar Surface Area (TPSA) are utilized to estimate the compound’s pharmacokinetic behavior. Red points in the bioavailability radar chart for compound 5l are almost entirely encompassed within the ideal property range (pink area), which indicate a favorable oral bioavailability.

3. Materials and Methods

3.1. Chemistry

3.1.1. General

All chemicals and organic solvents utilized to synthesize the desired compounds were of commercial grade, purchased from the Alfa Aesar (Cairo, Egypt), Sigma Aldrich (Darmstadt, Germany), and El-gomhoria (Cairo, Egypt) companies for pharmaceutical use, and were used without further purification. Reaction progress was monitored using pre-coated TLC plates (Kiesslgel 60 F254, Merck, Rahway, NJ, USA) and detecting spots was performed via UV lamp at a wavelength of 254 nm. Melting points were measured on as Stuart Electro-Thermal apparatus and recorded without correction. 1HNMR was carried out on a JOEL JNM-ECA400 (1HNMR: 400 MHz, 13CNMR: 100 MHz) spectrophotometer at Zagazig University, Egypt. Chemical shifts were recorded in parts per million (ppm), compared to the standard TMS using deuterated DMSO (δ: 2.5 ppm) as a solvent, and the coupling constant (J) was measured in hertz (Hz). Multiplicity was designated as: s, singlet; d, doublet; t, triplet; dd, doublet of doublet; and m for multiplet. Elemental analyses were performed on a Perkin Elmer 2400 CHN elemental analyzer at the regional center for mycology and biotechnology, Al-Azhar University, Nasr City, Cairo, Egypt.

3.1.2. General Procedure for the Synthesis of Compounds 3a–d

A mixture of the appropriate alkyl halide (3a–d, 1 mmol), p-hydroxy benzaldehyde 2a (1 mmol, 122 mg), and potassium carbonate (1 mmol, 138 mg) was heated to reflux in acetonitrile (20 mL) for 2 h. The reaction progress was monitored via TLC, then the reaction mixture was concentrated under reduced pressure and the remaining mixture was extracted using ethyl acetate (3 × 30 mL) and washed with water (3 × 30 mL). The organic layer was collected and dried over anhydrous sodium sulfate, the solvent was filtered off and then evaporated to yield compounds 3a–d that were used for further reactions without additional purification [31].
4-(Benzyloxy) benzaldehyde (3a)
White solid; 71% yield (160 mg); mp = 67–69 °C, reported mp: 68–70 °C; Rf = 0.3 (pet.ether:EtOAc, 3:1) [32,33].
4-(Pyridin-2-ylmethoxy) benzaldehyde (3b)
White solid; 53% yield (110 mg); mp = 95–98 °C, reported mp: 93–95 °C; Rf = 0.3 (pet.ether:EtOAc, 3:1) [34].
4-(Pyridin-3-ylmethoxy) benzaldehyde (3c)
Yellow solid; 53% yield (110 mg); mp = 75–78 °C, reported mp: 76–78 °C; Rf = 0.3 (pet.ether:EtOAc, 3:1) [35].
4-(Pyridin-4-ylmethoxy) benzaldehyde (3d)
White solid; 53% yield (110 mg); mp = 99–103 °C, reported mp: 102–104 °C; Rf = 0.3 (pet.ether:EtOAc, 3:1) [36].

3.1.3. General Procedure for the Synthesis of Compounds 5a–p

Indolin-2-ones (4a–d,1 mmol) were added to a solution of appropriate alkylated p-hydroxy benzaldehyde derivatives (3a–d, 1 mmol) in EtOH (10 mL). A catalytic amount of piperidine was added as base, and the reaction mixture was heated at reflux for 5 h. Reaction progress was monitored using TLC. After completion, the reaction mixture was allowed to cool. The formed precipitate was filtered and washed with cold EtOH then recrystallized using absolute ethanol to yield compounds 5a–p [37].
(E/Z)-3-[4-(Benzyloxy) benzylidene] indolin-2-one (5a)
Yellow powder; 89% yield (125 mg); mp: 220–223 °C, Rf = 0.42 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 3:97
Z-diastereomer: 10.59 (s, 1H, NH), 8.49 (d, J = 8.9 Hz, 2H, Ar-H), 7.75 (s, 1H, C=CH), 7.69 (d, J = 7.5 Hz, 1H, Ar-H), 7.50 (d, J = 7.2 Hz, 2H, Ar-H), 7.44 (t, J = 7.3 Hz, 2H, Ar-H), 7.37 (t, J = 7.2 Hz, 1H, Ar-H), 7.20 (t, J = 7.3 Hz, 1H, Ar-H), 7.13 (d, J = 8.9 Hz, 2H, Ar-H), 7.00 (t, J = 7.4 Hz, 1H, Ar-H), 6.83 (d, J = 7.7 Hz, 1H, Ar-H), 5.21 (s, 2H, CH2).
E-diastereomer: 10.59 (s, 1H, NH), 8.47 (d, J = 8.9 Hz, 2H, Ar-H), 7.69 (t, J = 7.2 Hz,1H, Ar-H), 7.67 (d, J = 7.5 Hz, 1H, Ar-H), 7.48 (d, J = 7.2 Hz, 2H, Ar-H), 7.40 (t, J = 7.3 Hz, 2H, Ar-H), 7.34 (s, 1H, C=CH), 7.17 (t, J = 7.3 Hz, 1H, Ar-H), 7.11 (d, J = 8.9 Hz, 2H, Ar-H), 6.96 (t, J = 7.4 Hz, 1H, Ar-H), 6.81 (d, J = 7.7 Hz, 1H, Ar-H), 5.21 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 167.37, 160.25, 140.26, 136.75, 136.68, 134.38, 128.48, 128.25, 127.98, 127.84, 127.17, 125.33, 124.13, 120.89, 119.26, 114.55, 109.19, 69.38.
E-diastereomer: 167.37, 160.25, 140.26, 136.75, 136.68, 134.38, 128.48, 128.25, 127.98, 127.84, 127.17, 125.33, 124.13, 120.89, 119.26, 114.55, 109.19, 69.38.
Anal. Calcd. For C22H17NO2 (327.13): C, 80.71; H, 5.23; N, 4.28, found; C, 80.36; H, 5.42; N, 4.31.
(E/Z)-3-[4-(Benzyloxy) benzylidene]-5-fluoroindolin-2-one (5b)
Yellow powder; 82% yield (91.3 mg); mp: 244–248 °C, Rf = 0.41 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 15:85
Z-diastereomer: 10.59 (s, 1H, NH), 8.50 (d, J = 8.9 Hz, 2H, Ar-H), 7.82 (s, 1H, C=CH), 7.62 (t, J = 9.1, 1H, Ar-H), 7.49–7.47 (m, 2H, Ar-H), 7.44–7.41 (m, 2H, Ar-H), 7.38–7.35 (m, 1H, Ar-H), 7.14 (d, J = 8.9 Hz, 2H, Ar-H), 7.04–7.01 (m, 1H, Ar-H), 6.80 (t, J = 8.5 Hz, 1H, Ar-H), 5.21 (s, 2H, CH2).
E-diastereomer: δ 10.47 (s, 1H, NH), 8.47 (d, J = 8.9 Hz, 2H, Ar-H), 7.82 (m, 1H, Ar-H), 7.59 (t, J = 9.1, 1H, Ar-H), 7.47–7.45 (m, 2H, Ar-H), 7.40–7.38 (m, 2H, Ar-H), 7.33 (s, 1H, C=CH), 7.12 (d, J = 8.9 Hz, 2H, Ar-H), 7.00–6.96 (m, 1H, Ar-H), 6.77 (t, J = 8.5 Hz, 1H, Ar-H), 5.21 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 167.49, 160.65, 160.00, 138.57, 137.76, 136.67, 136.48, 134.76, 131.68, 128.56, 128.07, 127.93, 127.91, 127.00, 126.46, 115.26, 114.71, 69.54.
E-diastereomer: 168.97, 159.17, 158.40, 156.84, 139.02, 134.76, 131.68, 128.56, 127.93, 123.83, 122.22, 115.90, 114.53, 110.73, 110.01, 109.00, 106.89, 69.48.
Anal. Calcd. For C22H16FNO2 (345.12): C, 76.51; H, 4.67; N, 4.06; found: C, 76.23; H, 4.98; N, 4.20.
(E/Z)-3-[4-(Benzyloxy) benzylidene]-5-chloroindolin-2-one (5c)
Yellow powder; 65% yield (87 mg); mp = 214–216 °C, Rf = 0.44 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 35:65
Z-diastereomer: 10.71 (s, 1H, NH), 8.49 (d, J = 8.9 Hz, 2H, Ar-H), 7.89 (s, 1H, C=CH), 7.65–7.54 (m, 1H, Ar-H), 7.49–7.45 (m, 2H, Ar-H), 7.39 (t, J = 7.2 Hz, 2H, Ar-H), 7.35–7.32 (m, 1H, Ar-H), 7.19 (t, J = 8.4 Hz, 2H, Ar-H), 7.14 (d, J = 8.9 Hz, 1H, Ar-H), 6.90–6.88 (t, J = 8.3 Hz, 1H, Ar-H), 5.20 (s, 2H, CH2).
E-diastereomer: 10.71 (s, 1H, NH), 8.62 (d, J = 8.9 Hz, 2H, Ar-H), 7.93–7.85 (m, 1H, Ar-H), 7.75–7.66 (m, 1H, Ar-H), 7.56–7.50 (m, 2H, Ar-H), 7.41 (t, J = 7.2 Hz, 2H, Ar-H), 7.33 (s, 1H, C=CH), 7.24 (t, J = 8.4 Hz, 2H, Ar-H), 7.12 (d, J = 8.9 Hz, 1H, Ar-H), 6.82–6.80 (t, J = 8.3 Hz, 1H, Ar-H), 5.20 (s, 2H, CH2).
  • 13C NMR (101 MHz, DMSO-d6) δ:
Z-diastereomer: 167.14, 160.66, 138.87, 138.77, 136.61, 134.77, 131.62, 128.49, 128.00, 127.89, 127.86, 125.28, 124.87, 119.31, 115.17, 114.66, 110.53, 69.48.
E-diastereomer: 168.56, 159.98, 141.43, 137.90, 131.62, 129.97,128.49, 127.57, 127.32, 127.00, 126.40, 125.28, 124.84, 123.00, 122.81, 121.43, 111.39, 69.43.
Anal. Calcd. For C22H16ClNO2 (361.09): C, 73.03; H, 4.46; N, 3.87. found: C, 72.71; H, 4.17; N, 3.92.
(E/Z)-3-[4-(Benzyloxy) benzylidene]-6-chloroindolin-2-one (5d)
Yellow powder; 80% yield (94 mg); mp: 259–261 °C, Rf = 0.43 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 15:85
Z-diastereomer: 10.74 (s, 1H, NH), 8.47 (d, J = 8.9 Hz, 2H, Ar-H), 7.80 (s, 1H, C=CH), 7.70 (d, J = 8.2 Hz, 1H, Ar-H), 7.49–7.45 (d, J = 7.1 Hz, 2H, Ar-H), 7.41 (t, J = 7.3 Hz, 2H, Ar-H), 7.40–7.37 (t, J = 7.1 Hz, 1H, Ar-H), 7.14 (d, J = 8.9 Hz, 2H, Ar-H), 7.06 (d, J = 1.8 Hz, 1H, Ar-H), 6.83 (d, J = 1.8 Hz, 1H, Ar-H), 5.21 (s, 2H, CH2).
E-diastereomer: 10.74 (s, 1H, NH), 8.45 (d, J = 8.9 Hz, 2H, Ar-H), 7.80 (t, J = 7.1 Hz, 1H, Ar-H), 7.68 (d, J = 8.2 Hz, 1H, Ar-H), 7.44–7.42 (d, J = 7.1 Hz, 2H, Ar-H), 7.36 (t, J = 7.3 Hz, 2H, Ar-H), 7.30 (s, 1H, C=CH), 7.11 (d, J = 8.9 Hz, 2H, Ar-H), 7.02 (d, J = 1.8 Hz, 1H, Ar-H), 6.82 (d, J = 1.8 Hz, 1H, Ar-H), 5.21 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 167.24, 160.51, 141.32, 137.90, 136.62, 134.58, 132.24, 131.66, 128.48, 127.99, 127.85, 127.00, 124.32, 122.91, 120.66, 114.62, 109.13, 69.41.
E-diastereomer: 167.24, 141.32,137.90, 136.62, 134.58, 133.61, 132.24, 131.66, 128.48, 127.99,127.85, 124.32, 122.91, 120.66, 115.15, 114.62, 109.13, 69.41.
Anal. Calcd. For C22H16ClNO2 (361.09): C, 73.03; H, 4.46; N, 3.87. found: C, 72.76; H, 4.17; N, 3.68.
(E/Z)-3-[4-(Pyridin-2-ylmethoxy] benzylidene) indolin-2-one (5e)
Orange powder; 76% yield (90 mg); mp: 221–224 °C, Rf = 0.42 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 31:69
Z-diastereomer: 10.58 (s, 1H, NH), 8.60 (d, 1H, Ar-H), 8.48 (d, J = 8.9 Hz, 1H, Ar-H), 7.90–7.86 (m, 1H, Ar-H), 7.78–7.74 (m, 1H, Ar-H), 7.68 (t, J = 8.5 Hz, 1H, Ar-H), 7.58–7.56 (t, J = 6.4 Hz, 1H, Ar-H), 7.36 (s, 1H, C=CH), 7.26 –7.18 (m, 3H, Ar-H), 6.99–6.97 (t, J = 7.5 Hz, 1H, Ar-H), 6.91–6.88 (m, 1H, Ar-H), 6.82 (d, J = 7.7 Hz, 1H, Ar-H), 5.27 (s, 2H, CH2).
E-diastereomer: 10.56 (s, 1H, NH), 8.59 (d, 1H, Ar-H), 8.46 (d, J = 8.9 Hz, 1H, Ar-H), 7.85–7.81 (m, 1H, Ar-H), 7.73–7.71 (m, 3H, Ar-H), 7.64 (t, J = 8.5 Hz, 1H, Ar-H), 7.55–7.51 (t, J = 6.4 Hz, 1H, Ar-H), 7.35–7.31 (m, 1H, Ar-H), 7.22 ((s, 1H, C=CH), 6.97–6.95 (t, J = 7.5 Hz, 1H, Ar-H), 6.87–6.84 (m, 1H, Ar-H), 6.80 (d, J = 7.7 Hz, 1H, Ar-H), 5.27 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 168.89, 159.41, 156.31, 149.24, 142.70, 137.07, 135.87, 134.39, 131.53, 129.75, 127.04, 125.80, 123.11, 121.85, 121.06, 115.08, 110.05, 109.21, 70.57.
E-diastereomer: 168.89, 160.01, 159.41, 149.19, 142.70, 140.29, 137.05, 136.67, 128.30, 125.31, 124.28, 121.12, 120.91, 119.29, 114.54, 110.05, 109.21, 70.40, 56.04.
Anal. Calcd. For C21H16N2O2 (328.37) C, 76.81; H, 4.91; N, 8.53. found: C, 77.40; H, 5.21; N, 8.14.
(E/Z)-5-fluoro-3-[4-(pyridin-2-ylmethoxy) benzylidene]indolin-2-one (5f)
Yellow powder; 93% yield (170 mg); mp: 226–228 °C, Rf = 0.40 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 33:67
Z-diastereomer: 10.59 (s, 1H, NH), 8.61–8.59 (d, 1H, Ar-H), 7.90–7.86 (m, 2H, Ar-H), 7.73 (d, J = 8.6 Hz, 2H, Ar-H), 7.65 (s, 1H, C=CH), 7.43–7.36 (m, 2H, Ar-H), 7.22 (t, J = 8.9 Hz, 2H, Ar-H), 7.16–7.01 (m, 1H, Ar-H), 6.87–6.82 (t, J = 8.5 Hz, 1H), 5.29 (s, 2H, CH2).
E-diastereomer: 10.58 (s, 1H, NH), 8.49–8.47 (d, 1H, Ar-H), 7.85–7.81 (m, 2H, Ar-H), 7.70 (d, J = 8.6 Hz, 2H, Ar-H), 7.58–7.53 (m, 2H, Ar-H), 7.38 (s, 1H, C=CH), 7.20 (t, J = 8.9 Hz, 2H, Ar-H), 7.00–6.96 (m, 1H, Ar-H), 6.80–6.77 (t, J = 8.5 Hz, 1H), 5.29 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ
Z-diastereomer: 168.86, 167.41, 160.37, 159.72, 156.28, 156.24, 149.20, 138.43, 137.60, 137.07, 134.70, 131.64, 127.14, 126.62, 123.12, 121.87, 115.22, 114.65, 70.49.
E-diastereomer: 158.34, 156.78, 156.00, 139.02, 137.43, 136.49, 134.70, 123.89, 122.06, 120.75, 116.10, 115.87, 114.26, 110.65, 109.93, 108.99, 106.89, 70.49, 52.16.
Anal. Calcd. For C21H15FN2O2 (346.36) C, 72.82; H, 4.37; N, 8.09. found: C, 73.13; H, 4.08; N, 8.20.
(E/Z)-5-Chloro-3-[4-(pyridin-2-ylmethoxy) benzylidene] indolin-2-one (5g)
Yellow powder; 89% yield (150 mg); mp: 219–220 °C, Rf = 0.35 (pet. ether: EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 23:77
Z-diastereomer: 10.71 (s, 1H, NH), 8.61–8.59 (d, 1H, Ar-H), 8.47 (d, 1H, Ar-H), 7.93–7.85 (m, 2H, Ar-H), 7.66 (s, 1H, C=CH), 7.59–7.55 (t, J = 8.8 Hz, 2H, Ar-H), 7.41–7.36 (s, 1H, Ar-H), 7.29–7.24 (t, J = 8.3 Hz, 1H, Ar-H), 7.16 (t, J = 8.9 Hz, 2H, Ar-H), 6.90–6.87 (t, J = 8.3 Hz, 1H, Ar-H), 5.28 (s, 2H, CH2).
E-diastereomer: 10.71 (s, 1H, NH), 8.50–8.47 (t, J = 8.8 Hz, 2H, Ar-H), 7.84–7.80 (m, 2H, Ar-H), 7.69–7.63 (m, 2H, Ar-H), 7.57 (s, 1H, C=CH), 7.35–7.32 (s, 1H, Ar-H), 7.23–7.20 (t, J = 8.3 Hz, 1H, Ar-H), 7.14 (t, J = 8.9 Hz, 2H, Ar-H), 6.82–6.80 (t, J = 8.3 Hz, 1H, Ar-H), 5.28 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 187.44, 176.19, 166.09, 148.62, 142.12, 137.19, 134.07, 132.62, 131.88, 130.74, 129.99, 123.89, 122.72, 118.88, 111.21, 107.14, 106.09, 52.05, 49.10.
E-diastereomer: 187.44, 176.19, 176.19, 166.09, 148.62, 142.12, 137.19, 134.07, 132.62, 131.88, 130.74, 129.99, 123.89, 122.72, 118.88, 111.21, 107.14, 106.09, 52.05, 49.10.
Anal. Calcd. For C21H15ClN2O2 (362.81) C, 69.52; H, 4.17; N, 7.72. found: C, 69.24; H, 4.02; N, 7.58.
(E/Z)-6-Chloro-3-[4-(pyridin-2-ylmethoxy) benzylidene]indolin-2-one (5h)
Yellow powder; 90% yield (213 mg); mp: 260–263 °C, Rf = 0.38 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 47:53
Z-diastereomer: 10.74 (s, 1H, NH), 8.64–8.60 (m, 2H, Ar-H), 8.47 (d, J = 8.9 Hz, 2H, Ar-H), 7.81 (s, 1H, C=CH), 7.72 (t, J = 7.7 Hz, 1H, Ar-H), 7.57 (t, J = 6.3 Hz, 1H, Ar-H), 7.41–7.37 (m, 1H, Ar-H), 7.15 (d, J = 9.0 Hz, 2H, Ar-H), 7.04 (dd, J = 1.9 Hz, 1H, Ar-H), 6.83 (d, J = 1.9 Hz, 1H, Ar-H), 5.28 (s, 2H, CH2).
E-diastereomer: 10.74 (s, 1H, NH), 8.59 (t, J = 6.3 Hz, 1H, Ar-H), 8.45 (d, J = 8.9 Hz, 2H, Ar-H), 7.84–7.79 (m, 2H, Ar-H), 7.71 (t, J = 7.7 Hz, 1H, Ar-H), 7.59 (s, 1H, C=CH), 7.36–7.33 (m, 1H, Ar-H), 7.13 (d, J = 9.0 Hz, 2H, Ar-H), 7.02 (dd, J = 1.9 Hz, 1H, Ar-H), 6.82 (d, J = 1.9 Hz, 1H, Ar-H), 5.28 (s, 2H, CH2).
  • 13C NMR (101 MHz, DMSO-d6) δ:
Z-diastereomer: 167.69, 160.09, 149.66, 144.45, 138.29, 137.55, 137.27, 135.05, 134.11, 132.15, 127.64, 127.23, 125.08, 123.58, 123.52, 122.33, 115.64, 115.09, 70.96.
E-diastereomer: 162.08, 160.09, 149.66, 144.45, 138.29, 137.52, 137.27, 135.05, 132.75, 123.80, 123.52, 121.29, 121.15, 120.53, 115.09, 115.41, 110.44, 109.62, 70.89.
Anal. Calcd. For C21H15ClN2O2 (362.81) C, 69.52; H, 4.17; N, 7.72. found: C, 69.23; H, 4.09; N, 7.53.
(E/Z)-3-[4-(Pyridin-3-ylmethoxy)benzylidene]indolin-2-one (5i)
Yellow powder; 66% yield (57.5 mg); mp: 163–161 °C, Rf = 0.32 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 44:56
Z-diastereomer: 10.59 (s, 1H, NH), 8.58 (m, 3H, Ar-H), 7.92 (d, J = 7.8 Hz, 1H, Ar-H), 7.75–7.71 (m, 3H, Ar-H), 7.58 (s, 1H, C=CH), 7.28–7.18 (m, 3H, Ar-H), 7.02–6.82 (m, 2H, Ar-H), 5.25 (s, 2H, CH2).
E-diastereomer: 10.56 (s, 1H, NH), 8.56 (m, 3H, Ar-H), 7.90 (d, J = 7.8 Hz, 1H, Ar-H), 7.68–7.57 (m, 3H, Ar-H), 7.44–7.41 (m, 3H, Ar-H), 7.13 (s, 1H, C=CH), 6.81–6.73 (m, 2H, Ar-H), 5.25 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 168.89, 159.39, 149.31, 149.21, 142.71, 135.87, 135.85, 134.39, 132.29, 131.51, 129.76, 127.07, 125.85, 123.66, 122.09, 121.13, 121.05, 115.08, 114.56, 110.07, 67.16.
E-diastereomer: 167.38, 165.80, 159.98, 159.39, 140.30, 136.67, 128.30, 127.38, 125.31, 124.31, 120.92, 119.29, 116.77, 109.52, 109.22, 67.10, 64.79, 59.76, 54.90, 48.60, 14.08.
Anal. Calcd. For C21H16N2O2 (328.37) C, 76.81; H, 4.91; N, 8.53. found: C, 76.53; H, 5.21; N, 8.24.
(E/Z)-5-Fluoro-3-[4-(pyridin-3-ylmethoxy)benzylidene]indolin-2-one (5j)
Yellow powder; 75% yield (133 mg); mp: 187–189 °C, Rf = 0.42 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 38:62
Z-diastereomer: 10.59 (s, 1H, NH), 8.98–8.67 (m, 2H, Ar-H), 8.16–7.66 (m, 3H, Ar-H), 7.57 (m, 2H, Ar-H), 7.45 (s,1H, C=CH), 7.13 (m, 3H, Ar-H), 6.85 (d, J = 8.9 Hz, 1H), 5.26 (s, 2H, CH2).
E-diastereomer: 10.59 (s, 1H, NH), 8.66–8.28 (m, 2H, Ar-H), 8.14–7.66 (m, 3H, Ar-H), 7.55 (m, 2H, Ar-H 7.40 (m, 2H, Ar-H), 7.12 (m, 1H, Ar-H), 7.00 (s, 1H, C=CH), 6.83 (d, J = 8.9 Hz, 1H), 5.26 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 168.89, 167.44, 160.34, 159.71, 149.33, 149.23, 139.04, 138.44, 137.62, 136.50, 135.89, 134.72, 131.64, 127.18, 123.96, 123.70, 116.13, 115.22, 114.67, 106.89, 67.21.
E-diastereomer: 156.03, 136.49, 132.25, 126.67, 123.70, 122.17, 115.89,110.77, 109.89, 109.23, 108.98, 67.15, 54.91, 40.15, 39.93, 39.73, 39.51, 39.31, 39.10, 38.89, 14.10.
Anal. Calcd. For C21H15FN2O2 (346.36) C, 72.82; H, 4.37; N, 8.09. found: C, 72.51; H, 4.18; N, 8.24.
(E/Z)-5-Chloro-3-[4-(pyridin-3-ylmethoxy)benzylidene] indolin-2-one (5k)
Yellow powder; 84% yield (120 mg); mp: 199–201 °C, Rf = 0.43 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 41:59
Z-diastereomer: 10.71 (s, 1H, NH), 8.72 (d, J = 9.0 Hz, 1H, Ar-H), 8.64 (d, J = 8.2 Hz, 2H, Ar-H), 8.54 (m, 1H, Ar-H), 7.93–7.89 (m, 2H, Ar-H), 7.67 (s, 1H, C=CH), 7.55 (d, J = 1.9 Hz, 1H, Ar-H), 7.37–7.20 (m, 3H, Ar-H), 6.85 (d, J = 8.3 Hz, 1H, Ar-H), 5.27 (s, 2H, CH2).
E-diastereomer: 10.71 (s, 1H, NH), 8.70 (d, J = 9.0 Hz, 1H, Ar-H), 8.55–8.47 (s, 1H, Ar-H), 7.84–7.81 (d, J = 8.2 Hz, 2H, Ar-H), 7.70–7.63 (m, 2H, Ar-H), 7.67 (m, 3H, Ar-H), 7.54 (d, J = 1.9 Hz, 1H, Ar-H), 7.29 (s, 1H, C=CH), 6.84 (d, J = 8.3 Hz, 1H, Ar-H), 5.24 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 160.38, 149.31, 149.20, 138.89, 138.67, 135.86, 135.83, 134.76, 132.21, 131.60, 129.22, 127.19, 126.62, 125.29, 124.88, 123.63, 123.15, 115.17, 114.66, 110.54, 108.84, 67.12.
E-diastereomer: 161.97, 160.38, 159.72, 149.31, 149.20, 141.44, 138.89, 138.67, 137.81, 131.17, 128.79, 127.61, 121.44, 119.33, 114.22, 111.39, 67.18, 64.26, 59.74, 54.89, 45.93, 20.74.
Anal. Calcd. For C21H15ClN2O2 (362.81) C, 69.52; H, 4.17; N, 7.72. found: C, 69.33; H, 4.02; N, 7.84.
(E/Z)-6-Chloro-3-[4-(pyridin-3-ylmethoxy) benzylidene] indolin-2-one (5l)
Yellow powder; 79% yield (125 mg); mp: 216–218 °C, Rf = 0.38 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 39:61
Z-diastereomer: 10.74 (s, 1H, NH), 8.71 (s, 1H, Ar-H), 8.57–8.55 (t, J = 5.9 Hz, 2H, Ar-H), 7.91 (d, J = 7.7 Hz, 1H, Ar-H), 7.80 (s,1H, C=CH), 7.50 (m, 2H, Ar-H), 7.46 (t, J = 7.6 Hz, 1H, Ar-H), 7.17 (t, J = 8.7 Hz, 2H, Ar-H), 7.08–6.79 (m, 2H, Ar-H), 5.26 (s, 2H, CH2).
E-diastereomer: δ 10.74 (s, 1H, NH), 8.71 (s, 1H, Ar-H), 8.49–8.46 (t, J = 5.9 Hz, 2H, Ar-H), 7.90 (d, J = 7.7 Hz, 1H, Ar-H), 7.63 (s,1H, C=CH), 7.46 (m, 2H, Ar-H), 7.44 (t, J = 7.6 Hz, 1H, Ar-H), 7.20 (t, J = 8.7 Hz, 2H, Ar-H), 7.08–6.92 (m, 2H, Ar-H), 5.26 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 168.83, 159.58, 149.31, 149.20, δ 144.00, 136.78, 135.83, 134.57, 133.66, 132.23,131.65, 126.80, 124.65, 123.65, 120.79, 120.70, 120.71, 115.15, 114.62,67.18, 67.12.
E-diastereomer: 167.23, 167.55, 161.89, 160.23, 160.06, 141.36, 137.79, 132.90, 132.30, 127.21, 124.29, 123.27, 123.08, 120.07, 109.98, 109.16, 64.41, 59.74, 54.89, 46.30, 14.07.
Anal. Calcd. For C21H15ClN2O2 (362.81) C, 69.52; H, 4.17; N, 7.72. found: C, 69.23; H, 4.37; N, 7.64.
(E/Z)-3-[4-(Pyridin-4-ylmethoxy) benzylidene] indolin-2-one (5m)
Yellow powder; 87% yield (160 mg); mp: 238–240 °C, Rf = 0.33 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 37:63
Z-diastereomer: 10.56 (s, 1H, NH), 8.61 (d, J = 4.9 Hz, 2H, Ar-H), 7.73 (d, J = 5.1 Hz, 2H, Ar-H), 7.68–7.59 (m, 3H, Ar-H), 7.58 (s, 1H, C=CH), 7.31–7.20 (m, 3H, Ar-H), 6.88 (t, J = 7.3 Hz, 2H), 5.29 (s, 2H, CH2).
E-diastereomer: 10.56 (s, 1H, NH), 8.59 (d, J = 4.9 Hz, 2H, Ar-H), 7.71 (m, 3H, Ar-H), 7.58–7.55 (m, 3H, Ar-H), 7.47 (d, J = 5.1 Hz, 2H, Ar-H), 7.20 (s, 1H, C=CH), 6.86 (t, J = 7.3 Hz, 2H), 5.29 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 169.34, 159.63, 150.26, 146.28, 143.19, 137.07, 136.28, 134.85, 132.00, 130.25, 127.67, 126.38, 122.57, 122.35, 121.58, 121.52, 121.39, 115.56, 115.03, 110.54, 68.17.
E-diastereomer: 166.27, 160.21, 150.26, 145.87, 136.73,136.02, 134.85, 133.96,131.17, 129.90, 128.81, 127.96, 126.68, 122.35, 120.04, 119.78, 118.42, 112.21, 109.69, 68.17, 51.50.
Anal. Calcd. For C21H16N2O2 (328.37) C, 76.81; H, 4.91; N, 8.53. found C, 77.19; H, 4.72; N, 8.38.
(E/Z)-5-Fluoro-3-[4-(pyridin-4-ylmethoxy)benzylidene] indolin-2-one (5n)
Yellow powder; 71% yield (116 mg); mp: 233–235 °C, Rf = 0.29 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 5:95
Z-diastereomer: 10.59 (s, 1H, NH), 8.61–8.59 (m, 2H, Ar-H), 7.73 (d, J = 8.6 Hz, 2H), 7.58 (s, 1H, C=CH), 7.48 (d, J = 6.0 Hz, 2H, Ar-H), 7.35 (t, J = 9.4 Hz, 1H), 7.21 (d, J = 8.8 Hz, 2H, Ar-H), 7.07 (t, J = 5.8 Hz, 1H, Ar-H), 6.85 (t, J = 4.7 Hz, 1H, Ar-H), 5.29 (s, 2H, CH2).
E-diastereomer: 10.59 (s, 1H, NH), 8.50–8.47 (m, 2H, Ar-H), 7.65 (s, 1H, C=CH), 7.58 (d, J = 8.8 Hz, 2H, Ar-H), 7.46 (d, J = 6.0 Hz, 2H, Ar-H), 7.33 (t, J = 9.4 Hz, 1H), 7.20 (d, J = 8.6 Hz, 2H), 7.05 (t, J = 5.8 Hz, 1H, Ar-H), 6.87 (t, J = 4.7 Hz, 1H, Ar-H), 5.29 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 169.32, 160.55, 159.92, 158.82, 150.25, 146.23, 138.81, 137.99, 136.97, 135.15, 132.10, 127.74, 127.25, 122.36, 116.36, 115.68, 115.12, 110.42, 109.45, 107.36, 68.19.
E-diastereomer:168.61, 162.23, 156.47, 143.28, 139.51, 138.39, 137.99, 135.15, 132.86,127.25, 126.29, 124.51, 123.08, 122.59, 118.52, 116.59, 111.80, 111.21, 110.10, 109.70, 26.21.
Anal. Calcd. For C21H15FN2O2 (346.36) C, 72.82; H, 4.37; N, 8.09. found: C, 72.50; H, 4.58; N, 8.26.
(E/Z)-5-Chloro-3-[4-(pyridin-4-ylmethoxy)benzylidene] indolin-2-one (5o)
Yellow powder; 81% yield (140 mg); mp: 238–240 °C, Rf = 0.30 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 57:43
Z-diastereomer: 10.71 (s, 1H, NH), 8.68–8.59 (m, 2H, Ar-H), 7.88 (s, 1H, C=CH), 7.66 (t, J = 5.3 Hz, 2H, Ar-H), 7.53 (d, J = 1.9 Hz, 1H, Ar-H), 7.46 (t, J = 6.0 Hz, 2H, Ar-H), 7.28 (t, J = 8.3 Hz, 1H, Ar-H), 7.17 (t, J = 8.8 Hz, 2H, Ar-H), 6.89 (t, J = 8.3 Hz, 1H, Ar-H), 5.29 (s, 2H, CH2).
E-diastereomer: 10.71 (s, 1H, NH), 8.51–8.47 (m, 2H, Ar-H), 7.66 (s, 1H, C=CH), 7.64 (t, J = 5.3 Hz, 2H, Ar-H), 7.48 (d, J = 1.9 Hz, 1H, Ar-H), 7.46 (t, J = 6.0 Hz, 2H, Ar-H), 7.19 (t, J = 8.3 Hz, 1H, Ar-H), 7.17 (t, J = 8.8 Hz, 2H, Ar-H), 6.85 (t, J = 8.3 Hz, 1H, Ar-H), 5.29 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 160.61, 150.25, 146.20, 139.38, 139.08, 138.23, 135.23, 132.09, 129.72, 127.78, 127.73, 127.24, 125.79, 123.73, 122.37, 122.35, 119.83, 115.66, 111.03, 68.19.
E-diastereomer: 159.96, 150.25,141.93, 139.08, 135.23, 129.72, 128.13, 127.73, 127.24, 125.54, 125.38, 123.73, 123.24, 122.37, 119.83, 116.45, 115.13, 111.88, 68.19, 68.13, 30.56.
Anal. Calcd. For C21H15ClN2O2 (362.81) C, 69.52; H, 4.17; N, 7.72. found: C, 69.10; H, 4.48; N, 7.57.
(E/Z)-6-Chloro-3-[4-(pyridin-4-ylmethoxy)benzylidene]indolin-2-one (5p)
Yellow powder; 86% yield (154 mg); mp: 233–235 °C, Rf = 0.28 (pet. ether:EtOAc, 1:1). 1H NMR (400 MHz, DMSO-d6) δ: E:Z ratio = 17:82
Z-diastereomer: 10.72 (s, 1H, NH), 8.61–8.59(t, J = 7.0 Hz, 2H, Ar-H), 7.84 (s,1H, Ar-H), 7.67 (d, J = 8.3 Hz, 2H, Ar-H), 7.62 (s, 1H, C=CH), 7.47 (d, J = 5.1 Hz, 2H, Ar-H), 7.15 (t, J = 8.6 Hz, 2H, Ar-H), 7.07–6.84 (m, 2H, Ar-H), 5.29 (s, 2H, CH2).
E-diastereomer: 10.72 (s, 1H, NH), 8.47–8.45 (t, J = 7.0 Hz, 2H, Ar-H), 7.70 (s,1H, Ar-H), 7.63 (d, J = 8.3 Hz, 2H, Ar-H), 7.59 (s, 1H, C=CH), 7.46 (d, J = 5.1 Hz, 2H, Ar-H), 7.12 (t, J = 8.6 Hz, 2H, Ar-H), 6.82 (m, 2H, Ar-H), 5.29 (s, 2H, CH2).
  • 13C NMR (100 MHz, DMSO-d6) δ:
Z-diastereomer: 169.30, 159.83, 150.26, 146.24, 144.48, 138.21, 137.20, 135.05, 134.17, 132.15, 127.40, 125.19, 123.78, 122.35, 121.28, 121.21, 120.52, 115.65, 115.10, 110.47, 68.19.
E-diastereomer: 168.53, 167.70, 162.97, 160.46, 145.54142.91, 141.85, 137.72,136.97, 133.99, 132.81, 127.79, 127.57, 126.62, 124.88, 124.74, 123.65, 122.10,112.84, 109.65, 68.19, 68.12.
Anal. Calcd. For C21H15ClN2O2 (362.81) C, 69.52; H, 4.17; N, 7.72. found: C, 69.83; H, 4.38; N, 7.53.

3.2. Biological Assays

3.2.1. Screening of Antiproliferative Activity by NCI

A total of 16 novel oxindole derivatives were sent to the National Cancer Institute (Germantown, MD, USA) for comprehensive anticancer screening [38]. Anticancer activity was tested against 60 cancer cell lines (NCI-60) after 48h treatments with a single dosage of 10 μM of compounds by using sulforhodamine B assay [39] on the traditional 96-well platform. The NCI-60 panel covers nearly 60 cancer cell models of leukemia, non-small cell lung cancer (NSCLC), colon, central nervous system (CNS), melanoma, ovarian, renal, prostate and breast cancers. The human cancer cell lines of the NCI-60 panel are grown in RPMI 1640 medium supplemented with 5% FBS (fetal bovine serum) and 2 mM L-glutamine. The detailed cell line list, cellular doubling times, and inoculation densities can be found on https://dtp.cancer.gov/discovery_development/nci-60/cell_list.htm (accessed on 10 November 2023). For further analysis of promising compounds, five dose screens were performed against the NCI-60 panel. Drug treatments are performed in 5 doses using 10-fold or ½ log dose dilutions. The relevant protocol can be accessed in full detail at https://dtp.cancer.gov/discovery_development/nci-60/methodology.htm (accessed on 10 November 2023).

3.2.2. Kinase Screening of Compound 5l

The inhibition percentages of the target compound 5l at a concentration of 10 µM were assessed using DMSO as a solvent against ten kinases, namely PTK2B, JAK1, CDK2, FGFR1, IGF1R, VEGFR-2, PDGFRα, PDGFRβ, FLT3, and Src. The Z’-LYTE® technology, which relies on FRET (Fluorescence Resonance Energy Transfer) from Invitrogen/Life Technologies, was employed for this evaluation for all kinases except VEGFR-2. The VEGFR-2 kinase assay was conducted in a 96-well streptavidin-coated plate using a recombinant human VEGFR-2/KDR ELISA kit, following the manufacturer’s instructions. Given the heightened activity of 5l against CDK2 and FLT3, the IC50 values of 5l, sunitinib, and FN-1501 against these two kinases were determined, with sunitinib and FN-1501 serving as the reference drugs. Statistical significance was assessed using one-way analysis of variance (ANOVA), followed by the Dunnett test; data are represented as means ± SD.

3.3. Molecular Modelling

3.3.1. Molecular Docking

Molecular docking commenced with the importation of the target protein structures from the Protein Data Bank (PDB), specifically using PDB IDs: 3TI1 for CDK2, and 6JQR for FLT3. These structures underwent preparation in the Maestro v13.3 interface, which included the removal of water molecules, the addition of hydrogen atoms, and the assignment of appropriate charge states to ensure an accurate physiological representation. Ligand preparation involved optimizing compound 5l’s geometry and confirming its tautomeric and protonation states, crucial for interaction fidelity. A grid representing the active site was created around the known binding pockets, using the co-crystallized ligands as a reference for CDK2 and FLT3. The Glide XP docking [40] protocol then systematically searched for the optimal conformation of compound 5l within these grids, utilizing an empirical scoring function to evaluate each potential binding mode. Validation was performed through the re-docking of the co-crystallized ligands, ensuring that the docking protocol could reliably reproduce known binding poses. Post-docking interaction analyses were facilitated by PLIP [41], which provided detailed insights into the binding interactions such as hydrogen bonds and hydrophobic contacts. The docked conformations were visualized using the PyMOL Molecular Graphics System v2.4, offering a three-dimensional perspective on the spatial arrangement of the interactions within the active sites.

3.3.2. MD Simulations

The MD simulations for the protein–ligand complexes of CDK2-5l and sunitinib, and of FLT3-5l and sunitinib were conducted using Gromacs 2021.2 [42], with the system setup prepared on the CHARMM-GUI [43,44] server using the CHARMM36m [45] force fields. The van der Waals interactions were managed using a cutoff type with a force-switch modifier. Here, a switching distance was set at 1.0 nm where the force-switching function begins, and the cutoff distance was set at 1.2 nm where the van der Waals interaction effectively becomes zero. This setup ensures a smooth transition in the computation of van der Waals forces, minimizing potential artifacts in the simulation. For electrostatic interactions, the Particle Mesh Ewald (PME) method was employed with a cutoff distance of 1.2 nm, ensuring accurate computation of long-range electrostatic forces. Temperature regulation was achieved via the Nose–Hoover thermostat, with separate temperature coupling groups for the solute and solvent at a target of 303.15 K. Pressure was maintained using the Parrinello–Rahman barostat with isotropic coupling, a compressibility setting of 4.5 × 10−5 bar−1, and a reference pressure of 1.0 atm. Constraints on hydrogen bonds were applied using the LINCS algorithm, facilitating a more stable and efficient simulation. The communication between solute and solvent was managed linearly with updates every 100 steps to maintain the center of mass. The ensuing 100 ns MD simulations captured the temporal evolution of the compound’s interaction within the active sites. For post-simulation analysis, the MM/PBSA method [46] was employed to calculate binding free energies, offering insights into the affinity of compound 5l with the kinase targets. MM/PBSA calculations were performed on 200 frames extracted between 80 and 100 ns, providing a comprehensive view of the binding energetics. Graphical representations of RMSD, PCA, and hydrogen bond dynamics were generated using QtGrace v0.2.6. Additionally, animations showcasing the dynamic molecular interactions within the MD simulations were created with PyMOL Molecular Graphics System v2.4.

3.3.3. Prediction of Physicochemical and Pharmacokinetic Properties

To predict the ADME parameters, pharmacokinetic properties, and the drug-likeness of the novel compound 5l, the SwissADME online tool by the SIB was utilized [47]. The molecular structures of the compounds were encoded into SMILES notations and inputted into the server. Toxicity profiles were ascertained using the pkCSM [48] pharmacokinetics website, facilitating a comprehensive assessment of the safety and efficacy potential of compound 5l. This predictive analysis forms the cornerstone of a prudent drug development strategy, allowing for the identification of favorable physicochemical characteristics and a manageable safety profile early in the drug development process.

4. Conclusions

A new series of oxindole-based derivatives linked to various pyridyl groups were synthesized and evaluated for their antiproliferative potential against the NCI 60 cancer cell panel. A subpanel of leukemia was selectively inhibited by compound 5l, which has nanomolar IC50 values for FLT3 (IC50 = 36.21 ± 1.07 nM), and CDK2 (IC50 = 8.17 ± 0.32 nM). A structure–activity–relationship investigation demonstrated that in order to inhibit FLT3 and CDK simultaneously, a 3-pyridyl group should be added to the linker’s 4-hydroxy group in addition to a chloro group at the 6-position of the oxindole ring. Molecular docking and dynamics have shown that compound 5l has a strong binding affinity and is very stable in the FLT3 and CDK2 binding pockets. These results point to the possibility of developing this novel scaffold into a medication for leukemia.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ph17050659/s1, Figure S1: RMSD profiles for compound 5l and sunitinib for each target CDK2 and FLT3 complexes, indicating deviations from initial conformations.; Figures S2–S34: 1H and 13C NMR spectrums of compounds 5a–p; Video S1: MD Simulation animation of 200 snapshots between 0 and 100 ns of CDK2 with compound 5l (PDB ID: 3TI1); Video S2: MD Simulation animation of 200 snapshots between 0 and 100 ns of CDK2 with sunitinib (PDB ID: 3TI1); Video S3: MD Simulation animation of 200 snapshots between 0 and 100 ns of FLT3 with compound 5l (PDB ID: 6JQR); Video S4: MD Simulation animation of 200 snapshots between 0 and 100 ns of FLT3 with sunitinib (PDB ID: 6JQR).

Author Contributions

Conceptualization, O.B. and T.F.S.A.; data curation, A.S., I.C., A.A.A.E.-H., E.M.N.A. and T.F.S.A.; formal analysis, A.S., O.B., I.C., A.A.A.E.-H., R.D. and T.F.S.A.; funding acquisition, E.B.E. and A.A.A.; investigation, A.S., O.B., R.D., E.B.E., A.A.A., E.M.N.A. and T.F.S.A.; methodology, A.S., O.B., A.A.A.E.-H., R.D., A.A., E.B.E., A.A.A., E.M.N.A., O.M.A., W.S. and T.F.S.A.; project administration, O.B., E.B.E., A.A.A. and T.F.S.A.; resources, I.C., A.A.A.E.-H. and T.F.S.A.; software, I.C.; supervision, O.B., A.A., O.M.A., W.S. and T.F.S.A.; validation, O.B., I.C., A.A., O.M.A., W.S. and T.F.S.A.; visualization, O.B. and T.F.S.A.; writing—original draft, A.S., I.C., A.A.A.E.-H. and T.F.S.A.; writing—review and editing, O.B., E.B.E., A.A., W.S. and T.F.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2024R116), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Dataset available on request from the authors.

Acknowledgments

The authors would like to express sincere gratitude to AlMaarefa University, Riyadh, Saudi Arabia, for supporting this research. The authors thank Ankara University-Scientific Research Unity for supplying the Schrödinger software version 2022.2 purchased under grant project number BAP-21B0237004. The molecular dynamics numerical calculations reported in this paper were partially performed at TUBITAK ULAKBIM in TURKEY, High Performance and Grid Computing Center (TRUBA resources). The authors highly thankful to Ramazan Ceylan from Sabanci University, Istanbul, Turkey for the valuable support.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hanahan, D.; Weinberg, R.A. Hallmarks of cancer: The next generation. Cell 2011, 144, 646–674. [Google Scholar] [CrossRef] [PubMed]
  2. Pavlova, N.N.; Thompson, C.B. The emerging hallmarks of cancer metabolism. Cell Metab. 2016, 23, 27–47. [Google Scholar] [CrossRef] [PubMed]
  3. Siegel, R.L.; Miller, K.D.; Wagle, N.S.; Jemal, A. Cancer statistics, 2023. CA Cancer J. Clin. 2023, 73, 17–48. [Google Scholar] [CrossRef]
  4. Fu, R.-G.; Sun, Y.; Sheng, W.-B.; Liao, D.-F. Designing multi-targeted agents: An emerging anticancer drug discovery paradigm. Eur. J. Med. Chem. 2017, 136, 195–211. [Google Scholar] [CrossRef]
  5. Giordano, S.; Petrelli, A. From single-to multi-target drugs in cancer therapy: When aspecificity becomes an advantage. Curr. Med. Chem. 2008, 15, 422–432. [Google Scholar] [CrossRef] [PubMed]
  6. Trajanoska, K.; Bhérer, C.; Taliun, D.; Zhou, S.; Richards, J.B.; Mooser, V. From target discovery to clinical drug development with human genetics. Nature 2023, 620, 737–745. [Google Scholar] [CrossRef]
  7. Aroosa, M.; Malik, J.A.; Ahmed, S.; Bender, O.; Ahemad, N.; Anwar, S. The evidence for repurposing anti-epileptic drugs to target cancer. Mol. Biol. Rep. 2023, 50, 7667–7680. [Google Scholar] [CrossRef]
  8. Atanasov, A.G.; Zotchev, S.B.; Dirsch, V.M.; Supuran, C.T. Natural products in drug discovery: Advances and opportunities. Nat. Rev. Drug Discov. 2021, 20, 200–216. [Google Scholar] [CrossRef]
  9. Mahomoodally, M.F.; Atalay, A.; Picot, M.C.N.; Bender, O.; Celebi, E.; Mollica, A.; Zengin, G. Chemical, biological and molecular modelling analyses to probe into the pharmacological potential of Antidesma madagascariense Lam.: A multifunctional agent for developing novel therapeutic formulations. J. Pharm. Biomed. Anal. 2018, 161, 425–435. [Google Scholar] [CrossRef]
  10. Nicolaou, K.C. Advancing the drug discovery and development process. Angew. Chem. 2014, 126, 9280–9292. [Google Scholar] [CrossRef]
  11. Paul, D.; Sanap, G.; Shenoy, S.; Kalyane, D.; Kalia, K.; Tekade, R.K. Artificial intelligence in drug discovery and development. Drug Discov. Today 2021, 26, 80. [Google Scholar] [CrossRef] [PubMed]
  12. Shen, B. A new golden age of natural products drug discovery. Cell 2015, 163, 1297–1300. [Google Scholar] [CrossRef]
  13. Alotaibi, A.A.; Asiri, H.H.; Rahman, A.F.M.M.; Alanazi, M.M. Novel pyrrolo[2,3-d]pyrimidine derivatives as multi-kinase inhibitors with VEGFR-2 selectivity. J. Saudi Chem. Soc. 2023, 27, 101712. [Google Scholar] [CrossRef]
  14. Wang, Z.; Cai, J.; Cheng, J.; Yang, W.; Zhu, Y.; Li, H.; Lu, T.; Chen, Y.; Lu, S. FLT3 Inhibitors in Acute Myeloid Leukemia: Challenges and Recent Developments in Overcoming Resistance. J. Med. Chem. 2021, 64, 2878–2900. [Google Scholar] [CrossRef] [PubMed]
  15. Ezelarab, H.A.A.; Ali, T.F.S.; Abbas, S.H.; Hassan, H.A.; Beshr, E.A.M. Indole-based FLT3 inhibitors and related scaffolds as potential therapeutic agents for acute myeloid leukemia. BMC Chem. 2023, 17, 73. [Google Scholar] [CrossRef] [PubMed]
  16. Al-Sanea, M.M.; Obaidullah, A.J.; Shaker, M.E.; Chilingaryan, G.; Alanazi, M.M.; Alsaif, N.A.; Alkahtani, H.M.; Alsubaie, S.A.; Abdelgawad, M.A. A New CDK2 Inhibitor with 3-Hydrazonoindolin-2-One Scaffold Endowed with Anti-Breast Cancer Activity: Design, Synthesis, Biological Evaluation, and In Silico Insights. Molecules 2021, 26, 412. [Google Scholar] [CrossRef] [PubMed]
  17. Li, C.; Liu, L.; Liang, L.; Xia, Z.; Li, Z.; Wang, X.; McGee, L.R.; Newhall, K.; Sinclair, A.; Kamb, A.; et al. AMG 925 is a dual FLT3/CDK4 inhibitor with the potential to overcome FLT3 inhibitor resistance in acute myeloid leukemia. Mol. Cancer Ther. 2015, 14, 375–383. [Google Scholar] [CrossRef] [PubMed]
  18. Richardson, G.E.; Al-Rajabi, R.; Uprety, D.; Hamid, A.; Williamson, S.K.; Baranda, J.; Mamdani, H.; Lee, Y.L.; Nitika; Li, L.; et al. A Multicenter, Open-Label, Phase I/II Study of FN-1501 in Patients with Advanced Solid Tumors. Cancers 2023, 15, 2553. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, Y.; Zhi, Y.; Jin, Q.; Lu, S.; Lin, G.; Yuan, H.; Yang, T.; Wang, Z.; Yao, C.; Ling, J.; et al. Discovery of 4-((7H-Pyrrolo[2,3-d]pyrimidin-4-yl)amino)-N-(4-((4-methylpiperazin-1-yl)methyl)phenyl)-1H-pyrazole-3-carboxamide (FN-1501), an FLT3- and CDK-Kinase Inhibitor with Potentially High Efficiency against Acute Myelocytic Leukemia. J. Med. Chem. 2018, 61, 1499–1518. [Google Scholar] [CrossRef]
  20. Bender, O.; Shoman, M.E.; Ali, T.F.S.; Dogan, R.; Celik, I.; Mollica, A.; Hamed, M.I.A.; Aly, O.M.; Alamri, A.; Alanazi, J.; et al. Discovery of oxindole-based FLT3 inhibitors as a promising therapeutic lead for acute myeloid leukemia carrying the oncogenic ITD mutation. Arch. Pharm. 2023, 356, e2200407. [Google Scholar] [CrossRef]
  21. Alanazi, J.; Bender, O.; Dogan, R.; Malik, J.A.; Atalay, A.; Ali, T.F.S.; Beshr, E.A.M.; Shawky, A.M.; Aly, O.M.; Alqahtani, Y.N.H.; et al. Combination of an Oxindole Derivative with (−)-beta-Elemene Alters Cell Death Pathways in FLT3/ITD(+) Acute Myeloid Leukemia Cells. Molecules 2023, 28, 5253. [Google Scholar] [CrossRef] [PubMed]
  22. Bender, O.; Celik, I.; Dogan, R.; Atalay, A.; Shoman, M.E.; Ali, T.F.S.; Beshr, E.A.M.; Mohamed, M.; Alaaeldin, E.; Shawky, A.M.; et al. Vanillin-Based Indolin-2-one Derivative Bearing a Pyridyl Moiety as a Promising Anti-Breast Cancer Agent via Anti-Estrogenic Activity. ACS Omega 2023, 8, 6968–6981. [Google Scholar] [CrossRef] [PubMed]
  23. Chen, X.; Yang, T.; Deivasigamani, A.; Shanmugam, M.K.; Hui, K.M.; Sethi, G.; Go, M.L. N′-Alkylaminosulfonyl Analogues of 6-Fluorobenzylideneindolinones with Desirable Physicochemical Profiles and Potent Growth Inhibitory Activities on Hepatocellular Carcinoma. ChemMedChem 2015, 10, 1548–1558. [Google Scholar] [CrossRef] [PubMed]
  24. Lackey, K.; Cory, M.; Davis, R.; Frye, S.V.; Harris, P.A.; Hunter, R.N.; Jung, D.K.; McDonald, O.B.; McNutt, R.W.; Peel, M.R. The discovery of potent cRaf1 kinase inhibitors. Bioorg. Med. Chem. Lett. 2000, 10, 223–226. [Google Scholar] [CrossRef] [PubMed]
  25. Mansour, H.S.; Abd El-wahab, H.A.; Ali, A.M.; Aboul-Fadl, T. Inversion kinetics of some E/Z 3-(benzylidene)-2-oxo-indoline derivatives and their in silico CDK2 docking studies. RSC Adv. 2021, 11, 7839–7850. [Google Scholar] [CrossRef] [PubMed]
  26. Sun, L.; Tran, N.; Tang, F.; App, H.; Hirth, P.; McMahon, G.; Tang, C. Synthesis and biological evaluations of 3-substituted indolin-2-ones: A novel class of tyrosine kinase inhibitors that exhibit selectivity toward particular receptor tyrosine kinases. J. Med. Chem. 1998, 41, 2588–2603. [Google Scholar] [CrossRef] [PubMed]
  27. Martin, M.P.; Alam, R.; Betzi, S.; Ingles, D.J.; Zhu, J.Y.; Schonbrunn, E. A novel approach to the discovery of small-molecule ligands of CDK2. Chembiochem 2012, 13, 2128–2136. [Google Scholar] [CrossRef] [PubMed]
  28. Kawase, T.; Nakazawa, T.; Eguchi, T.; Tsuzuki, H.; Ueno, Y.; Amano, Y.; Suzuki, T.; Mori, M.; Yoshida, T. Effect of Fms-like tyrosine kinase 3 (FLT3) ligand (FL) on antitumor activity of gilteritinib, a FLT3 inhibitor, in mice xenografted with FL-overexpressing cells. Oncotarget 2019, 10, 6111–6123. [Google Scholar] [CrossRef]
  29. Alafnan, A.; Dogan, R.; Bender, O.; Celik, I.; Mollica, A.; Malik, J.A.; Rengasamy, K.R.R.; Break, M.K.B.; Khojali, W.M.A.; Alharby, T.N. Beta Elemene induces cytotoxic effects in FLT3 ITD-mutated acute myeloid leukemia by modulating apoptosis. Eur. Rev. Med. Pharmacol. Sci. 2023, 27, 3270–3287. [Google Scholar]
  30. Daina, A.; Zoete, V. A BOILED-Egg To Predict Gastrointestinal Absorption and Brain Penetration of Small Molecules. ChemMedChem 2016, 11, 1117–1121. [Google Scholar] [CrossRef]
  31. Ebenezer, O.; Oyetunde-Joshua, F.; Omotoso, O.D.; Shapi, M. Benzimidazole and its derivatives: Recent Advances (2020–2022). Results Chem. 2023, 5, 100925. [Google Scholar] [CrossRef]
  32. Ye, Z.; Lei, Z.; Ye, X.; Zhou, L.; Wang, Y.; Yuan, Z.; Gao, F.; Britton, R. Decatungstate Catalyzed Synthesis of Trifluoromethylthioesters from Aldehydes via a Radical Process. J. Org. Chem. 2022, 87, 765–775. [Google Scholar] [CrossRef]
  33. Sun, Z.; Zaman, M.K.; Khan, S.N.; Cai, Y. Decarboxylative Oxidation of Carboxylic Acids Using Photocatalysis and Copper Catalysis. Synlett 2023, 34, 2029–2033. [Google Scholar] [CrossRef]
  34. Popov, A.B.; Krstulovic, L.; Kostrun, S.; Jelic, D.; Bokulic, A.; Stojkovic, M.R.; Zonjic, I.; Taylor, M.C.; Kelly, J.M.; Bajic, M.; et al. Design, synthesis, antitrypanosomal activity, DNA/RNA binding and in vitro ADME profiling of novel imidazoline-substituted 2-arylbenzimidazoles. Eur. J. Med. Chem. 2020, 207, 112802. [Google Scholar] [CrossRef] [PubMed]
  35. Sundberg, R.J.; Dahlhausen, D.J.; Manikumar, G.; Mavunkel, B.; Biswas, A.; Srinivasan, V.; King, F.; Waid, P. Preparation of 2-aryl and 2-aryloxymethyl imidazo[1,2-a]pyridines and related compounds. J. Heterocycl. Chem. 2009, 25, 129–137. [Google Scholar] [CrossRef]
  36. Taban, I.M.; Zhu, J.; DeLuca, H.F.; Simons, C. Analysis of the binding sites of vitamin D 1alpha-hydroxylase (CYP27B1) and vitamin D 24-hydroxylase (CYP24A1) for the design of selective CYP24A1 inhibitors: Homology modelling, molecular dynamics simulations and identification of key binding requirements. Bioorg. Med. Chem. 2017, 25, 5629–5636. [Google Scholar] [CrossRef] [PubMed]
  37. Amombo, G.M.; Kramer, T.; Lo Monte, F.; Goring, S.; Fach, M.; Smith, S.; Kolb, S.; Schubenel, R.; Baumann, K.; Schmidt, B. Modification of a promiscuous inhibitor shifts the inhibition from gamma-secretase to FLT-3. Bioorg. Med. Chem. Lett. 2012, 22, 7634–7640. [Google Scholar] [CrossRef]
  38. Workman, P. The NCI-60 Human Tumor Cell Line Screen: A Catalyst for Progressive Evolution of Models for Discovery and Development of Cancer Drugs. Cancer Res. 2023, 83, 3170–3173. [Google Scholar] [CrossRef]
  39. Vichai, V.; Kirtikara, K. Sulforhodamine B colorimetric assay for cytotoxicity screening. Nat. Protoc. 2006, 1, 1112–1116. [Google Scholar] [CrossRef]
  40. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra precision glide: Docking and scoring incorporating a model of hydrophobic enclosure for protein-ligand complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef]
  41. Adasme, M.F.; Linnemann, K.L.; Bolz, S.N.; Kaiser, F.; Salentin, S.; Haupt, V.J.; Schroeder, M. PLIP 2021: Expanding the scope of the protein-ligand interaction profiler to DNA and RNA. Nucleic Acids Res. 2021, 49, W530–W534. [Google Scholar] [CrossRef] [PubMed]
  42. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19–25. [Google Scholar] [CrossRef]
  43. Jo, S.; Kim, T.; Iyer, V.G.; Im, W. CHARMM-GUI: A web-based graphical user interface for CHARMM. J. Comput. Chem. 2008, 29, 1859–1865. [Google Scholar] [CrossRef]
  44. Lee, J.; Cheng, X.; Swails, J.M.; Yeom, M.S.; Eastman, P.K.; Lemkul, J.A.; Wei, S.; Buckner, J.; Jeong, J.C.; Qi, Y.; et al. CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field. J. Chem. Theory Comput. 2016, 12, 405–413. [Google Scholar] [CrossRef] [PubMed]
  45. Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B.L.; Grubmuller, H.; MacKerell, A.D., Jr. CHARMM36m: An improved force field for folded and intrinsically disordered proteins. Nat. Methods 2017, 14, 71–73. [Google Scholar] [CrossRef] [PubMed]
  46. Valdes-Tresanco, M.S.; Valdes-Tresanco, M.E.; Valiente, P.A.; Moreno, E. gmx_MMPBSA: A New Tool to Perform End-State Free Energy Calculations with GROMACS. J. Chem. Theory Comput. 2021, 17, 6281–6291. [Google Scholar] [CrossRef] [PubMed]
  47. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef]
  48. Pires, D.E.; Blundell, T.L.; Ascher, D.B. pkCSM: Predicting Small-Molecule Pharmacokinetic and Toxicity Properties Using Graph-Based Signatures. J. Med. Chem. 2015, 58, 4066–4072. [Google Scholar] [CrossRef]
Figure 1. Structure of sunitinib, FN-1501, III, IV, and targeted compounds (5a–p).
Figure 1. Structure of sunitinib, FN-1501, III, IV, and targeted compounds (5a–p).
Pharmaceuticals 17 00659 g001
Scheme 1. Synthetic approach of compounds 5a–p. Reagents and conditions: (a) K2CO3, ACN, reflux, 2 h; (b) Piperidine, EtOH, reflux, 5 h.
Scheme 1. Synthetic approach of compounds 5a–p. Reagents and conditions: (a) K2CO3, ACN, reflux, 2 h; (b) Piperidine, EtOH, reflux, 5 h.
Pharmaceuticals 17 00659 sch001
Figure 2. Dynamics of compound 5l on 60 cancer cell lines from 9 different cancer panels. (a) % growth inhibition exerted by compound 5l at 10 µM concentration over NCI-60 cancer cell lines. Red squares indicate the highest number of consistent inhibitory activity. (b) Radar representation of the inhibitory activity of compound 5l on nine types of cancers. The red line is shaped in relation to increasing activity outward from the radar center. (c) Variations in 5l inhibitory effect among cell lines in each panel of cancer in NCI-60. The turquoise line shows 20% as the threshold.
Figure 2. Dynamics of compound 5l on 60 cancer cell lines from 9 different cancer panels. (a) % growth inhibition exerted by compound 5l at 10 µM concentration over NCI-60 cancer cell lines. Red squares indicate the highest number of consistent inhibitory activity. (b) Radar representation of the inhibitory activity of compound 5l on nine types of cancers. The red line is shaped in relation to increasing activity outward from the radar center. (c) Variations in 5l inhibitory effect among cell lines in each panel of cancer in NCI-60. The turquoise line shows 20% as the threshold.
Pharmaceuticals 17 00659 g002
Figure 3. Five dose (0.01 to 100 μM) assay with compound 5l against NCI-60 cell lines. Cells were treated with five different doses of 5l for 48 h. Charts represent percentage growth between −100 and 100. Each cell line is indicated with different shapes and colors under the corresponding cancer chart. (a) Leukemia; (b) Non-small cell lung cancer; (c) Colon cancer; (d) CNS (central nervous system) cancer; (e) Melanoma; (f) Ovarian cancer; (g) Renal cancer; (h) Prostate cancer; (i) Breast cancer.
Figure 3. Five dose (0.01 to 100 μM) assay with compound 5l against NCI-60 cell lines. Cells were treated with five different doses of 5l for 48 h. Charts represent percentage growth between −100 and 100. Each cell line is indicated with different shapes and colors under the corresponding cancer chart. (a) Leukemia; (b) Non-small cell lung cancer; (c) Colon cancer; (d) CNS (central nervous system) cancer; (e) Melanoma; (f) Ovarian cancer; (g) Renal cancer; (h) Prostate cancer; (i) Breast cancer.
Pharmaceuticals 17 00659 g003
Figure 4. The percentage of inhibition of the tested compound 5l (10 µM) against PTK2B, JAK1, CDK2, FGFR1, IGF1R, VEDFR-2, PDGFRα, PDGFRβ, FLT3, and SRC kinases.
Figure 4. The percentage of inhibition of the tested compound 5l (10 µM) against PTK2B, JAK1, CDK2, FGFR1, IGF1R, VEDFR-2, PDGFRα, PDGFRβ, FLT3, and SRC kinases.
Pharmaceuticals 17 00659 g004
Figure 5. Dose titrations of 5l compared to sunitinib and FN-1501 against (a) FLT3 and (b) CDK2 kinases from 1 to 10,000 nM. (c) The half maximal inhibitory concentration (IC50) of compound 5l against FLT3 and CDK2 kinase activity compared to sunitinib and FN-1501. a IC50 values are of three separate experiments ± SD. * p ≤ 0.05, *** p ≤ 0.001, **** p ≤ 0.0001.
Figure 5. Dose titrations of 5l compared to sunitinib and FN-1501 against (a) FLT3 and (b) CDK2 kinases from 1 to 10,000 nM. (c) The half maximal inhibitory concentration (IC50) of compound 5l against FLT3 and CDK2 kinase activity compared to sunitinib and FN-1501. a IC50 values are of three separate experiments ± SD. * p ≤ 0.05, *** p ≤ 0.001, **** p ≤ 0.0001.
Pharmaceuticals 17 00659 g005
Figure 6. Docking interactions and dynamics simulations of compound 5l with kinase targets related to colorectal cancer and AML. (a) The conformation of compound 5l within CDK2 (PDB ID: 3TI1) active site, showing interaction points. (b) FLT3 (PDB ID: 6JQR) in complex with compound 5l, presenting the compound’s orientation and contact regions. (c,d) Molecular dynamics simulation of compound 5l with kinase targets over 100 ns. RMSD profiles for 5l for each target (green), CDK2 (red line) and FLT3 (blue line) complexes, indicating deviations from initial conformations. Conformational space exploration of 5l with each kinase, represented using PCA. Hydrogen bond count over the simulation period, reflecting interaction consistency for each kinase–ligand complex. (e) MM/PBSA binding free energy components for compound 5l with kinase targets. The bar graph summarizes the van der Waals, electrostatic, and solvation energy contributions to the total binding free energy (ΔTOTAL) of 5l with CDK2 and FLT3.
Figure 6. Docking interactions and dynamics simulations of compound 5l with kinase targets related to colorectal cancer and AML. (a) The conformation of compound 5l within CDK2 (PDB ID: 3TI1) active site, showing interaction points. (b) FLT3 (PDB ID: 6JQR) in complex with compound 5l, presenting the compound’s orientation and contact regions. (c,d) Molecular dynamics simulation of compound 5l with kinase targets over 100 ns. RMSD profiles for 5l for each target (green), CDK2 (red line) and FLT3 (blue line) complexes, indicating deviations from initial conformations. Conformational space exploration of 5l with each kinase, represented using PCA. Hydrogen bond count over the simulation period, reflecting interaction consistency for each kinase–ligand complex. (e) MM/PBSA binding free energy components for compound 5l with kinase targets. The bar graph summarizes the van der Waals, electrostatic, and solvation energy contributions to the total binding free energy (ΔTOTAL) of 5l with CDK2 and FLT3.
Pharmaceuticals 17 00659 g006
Figure 7. ADME estimation of active compound 5l. (a) The BOILED-Egg plot showcases the predicted GI absorption and BBB permeability for 5l, aligning with non-P-glycoprotein substrate characteristics (PGP−), and (b) the bioavailability radar chart, with the pink area depicting the optimal property values and the red lines representing the properties of 5l, which largely fall within the desirable range, indicating good oral bioavailability.
Figure 7. ADME estimation of active compound 5l. (a) The BOILED-Egg plot showcases the predicted GI absorption and BBB permeability for 5l, aligning with non-P-glycoprotein substrate characteristics (PGP−), and (b) the bioavailability radar chart, with the pink area depicting the optimal property values and the red lines representing the properties of 5l, which largely fall within the desirable range, indicating good oral bioavailability.
Pharmaceuticals 17 00659 g007
Table 1. E:Z diastereomeric ratios of the synthesized compounds (5a–p).
Table 1. E:Z diastereomeric ratios of the synthesized compounds (5a–p).
CompoundSubstituentsChemical Shifts (Vinyl H, δ ppm)Chemical Shifts (2′ and 6′-Hs, δ ppm)E:Z Ratio a
RArEZEZ
5aHphenyl7.347.757.367.683:97
5b5-Fphenyl7.337.827.357.6115:85
5c5-Clphenyl7.337.897.357.8135:65
5d6-Clphenyl7.307.807.367.6915:85
5eH2-pyridyl7.227.367.587.7331:69
5f5-F2-pyridyl7.387.657.607.7333:67
5g5-Cl2-pyridyl7.577.667.297.7223:77
5h6-Cl2-pyridyl7.597.817.387.7147:53
5iH3-pyridyl7.137.587.637.7244:56
5j5-F3-pyridyl7.007.457.368.5738:62
5k5-Cl3-pyridyl7.297.677.197.7641:59
5l6-Cl3-pyridyl7.637.807.628.4739:61
5mH4-pyridyl7.207.587.647.2437:63
5n5-F4-pyridyl7.657.587.007.735:95
5o5-Cl4-pyridyl7.667.887.737.8057:43
5p6-Clphenyl7.597.627.047.7317:82
a E:Z ratio was calculated according to the ratio between the relative integration of vinylic H for E and Z isomers in 1H NMR.
Table 2. Growth inhibition (%) induced by 10 µM of compounds 5a–p against NCI-60 cell lines panel (9 different cancer types).
Table 2. Growth inhibition (%) induced by 10 µM of compounds 5a–p against NCI-60 cell lines panel (9 different cancer types).
5a5b5c5d5e5f5g5h5i5j5k5l5m5n5o5p
LeukemiaCCRF-CEM0.52−10.840.21−2.1857.725.7916.8053.7252.4475.0036.3883.9340.5140.6353.8534.36
HL-60(TB)11.517.92−3.22−13.1438.7513.9314.995.3349.8781.6814.8295.5455.2657.7566.5152.44
K-56211.807.19−0.346.9159.2114.6348.9531.3575.9186.1338.7695.0562.6779.3272.4559.07
MOLT-42.73−0.053.280.8662.523.5553.1561.4769.9888.5352.6298.0878.4776.7381.2279.00
RPMI-82261.132.59−3.85−7.3663.9919.3826.1032.4145.3761.2329.6795.5171.4764.5170.5074.88
SR14.63−0.14−0.06−4.9357.72ND45.1042.7872.0472.9966.1184.2073.2985.4488.8872.56
Non-Small Cell Lung CancerA549/ATCC−3.263.90−0.99−4.5929.20−1.0435.6462.5326.1839.3616.1358.8034.5129.7939.9427.89
EKVX0.690.650.85−0.3147.3115.7734.2649.3731.7041.0416.0660.9363.5352.1637.6145.07
HOP-62−4.06−1.203.773.8525.62−7.7154.4418.9237.2445.440.5672.9023.0417.0419.6821.38
HOP-92−7.004.779.75−8.6768.9832.9327.9535.7735.3925.82−1.4056.0640.1936.9216.3049.65
NCI-H2264.459.5311.097.1340.806.6798.7448.6438.1655.5611.8062.9342.4044.4838.3049.48
NCI-H23−2.421.73−1.965.7121.3728.0137.7747.4342.4358.2016.7475.6049.8358.5655.3157.17
NCI-H322M0.442.81−0.070.4929.483.2935.7524.71NDNDNDND20.8527.5640.6422.52
NCI-H460−8.64−0.94−3.66−4.4814.9711.6084.0156.7221.7153.9414.3377.7340.5444.4057.7536.76
NCI-H5222.185.912.100.1924.5668.0833.9531.5849.9089.5319.75103.9632.7081.9886.0345.21
Colon CancerCOLO 205−13.051.92−6.62−9.74−14.73−10.8010.90−3.4726.1346.065.6696.3638.2622.0719.2139.71
HCC-2998−13.07−9.12−3.96−1.791.60−20.25−4.5116.74−2.0823.37−9.2150.8820.6425.7412.3920.13
HCT-1164.244.153.412.3341.1057.8371.5547.8256.5284.6734.6492.9940.5156.1260.5161.51
HCT-150.043.51−4.255.0045.9310.7064.2745.5060.5584.7746.5497.2462.3178.7579.1476.92
HT29−6.642.38−2.95−5.2716.23−0.5132.5319.6632.6684.0918.1090.3311.2536.4439.6937.78
KM12−3.152.42−1.024.9261.2718.8221.2521.6587.4497.3352.8099.7066.2769.8560.6447.83
SW-620−7.33−2.341.080.60ND−2.7952.53ND19.9570.102.6292.7321.9741.1048.3349.50
CNS CancerSF-268−1.745.433.265.5527.3329.8720.6823.17−19.5911.105.1045.3622.7138.8355.7116.72
SF-295−6.781.11−11.14−14.6627.031.9815.5214.2831.2054.814.5163.3158.2751.7044.1732.18
SF-539−3.404.674.125.4835.3447.6726.5057.6836.2063.2820.2359.2624.8841.7443.9594.23
SNB-192.171.26−1.05−3.0939.3012.9549.2619.9937.1262.0622.9256.1843.1360.4673.7828.72
SNB-7514.605.1513.0726.7034.7310.24−12.5663.32−1.11−5.67−10.2418.2835.2426.5721.2414.90
U251−3.490.40−1.54−4.2546.4541.7776.3061.5853.1783.6527.3376.3838.9261.2766.8964.39
MelanomaLOX IMVI−7.76−6.240.147.32ND12.6043.51ND54.9976.9110.0991.0041.3767.2658.7465.24
MALME-3M−7.597.33−12.99−15.3825.407.2461.5318.04NDNDNDND17.4126.7943.0312.53
M145.673.397.343.838.30−5.0910.0720.0037.6267.0612.1777.0823.7032.4933.0516.86
MDA-MB-435−2.335.59−3.790.5910.09−2.7723.1023.8660.9398.2829.29121.0649.6261.6758.7329.23
SK-MEL-2−4.59ND−0.35−11.32−12.23−3.45−11.36−10.2824.1050.380.3666.2327.1928.054.2327.64
SK-MEL-28−7.212.92−8.75−2.97−0.76−8.53−0.48−8.02−2.6528.18−20.4746.27−4.73−9.3914.48−10.92
SK-MEL-51.96−0.591.902.8348.36−5.8510.4923.9859.40112.6511.41152.5890.0293.1162.5664.97
UACC-257−9.452.53−6.58−16.0712.60−4.06−18.19−2.1030.5239.792.4276.6935.7923.4323.6917.22
UACC-626.114.385.694.0940.389.9433.2849.6243.4857.7823.2964.3131.0928.3537.6034.70
Ovarian CancerIGROV1−9.532.92−5.148.4136.846.3545.4026.16NDNDNDND25.3564.6584.7367.88
OVCAR-3NDNDNDND18.168.1742.686.16−9.249.05−10.9361.7821.4939.2467.3653.00
OVCAR-4−7.990.82−7.183.6530.835.38101.8175.8337.8347.2321.2368.4142.6743.4347.3431.80
OVCAR-5−2.24−0.99−2.92−0.555.53−9.4511.7615.63−27.58−11.18−34.5221.73−33.20−27.63−14.36-29.38
OVCAR-82.34−0.84−0.47−0.4129.683.5254.0613.6047.1965.1323.4588.6335.0256.2561.4761.29
NCI/ADR-RES−0.550.011.88−4.2446.7221.7085.8838.9854.9785.5726.04112.4061.8174.1667.6259.29
SK-OV-3−5.780.16−1.70−3.0116.833.7510.804.3913.1513.155.3632.249.3813.6911.4115.54
Renal Cancer786-03.670.43−2.86−4.9717.9716.6374.3018.6617.2524.201.4747.2514.6019.3732.3723.86
A4987.90ND−2.376.9816.4232.5154.92−11.54−17.75−29.11−31.568.902.1712.871.01−1.39
ACHN−4.27−7.30−1.455.1833.6223.7543.049.6137.5451.8416.3855.1324.1070.2770.2334.02
CAKI-12.008.2313.4210.7040.7533.4738.4428.7956.7367.6051.0369.8045.6362.3238.7152.78
RXF 393−16.170.31−5.38−14.2544.7824.20ND38.1824.8637.266.8463.4243.8847.6273.8467.40
SN12C6.160.270.704.0347.674.0455.4824.9027.1446.159.1766.0621.4041.3234.4134.65
TK-10−8.254.82−12.74−53.9812.674.479.89−0.24−15.67−1.92−37.8029.4517.9429.9517.72−0.51
UO-3114.902.5523.9223.8361.1626.5440.4538.05NDNDNDND44.7254.7551.1548.99
Prostate CancerPC-35.85−2.7712.8611.7544.338.7037.9354.6929.7643.6024.3067.1145.9647.3952.4144.98
DU-145−4.380.40−5.61−1.2520.297.7312.1839.5639.8246.090.7558.6133.7138.4449.2829.35
Breast CancerMCF715.767.5716.0038.4250.7837.7571.0863.6657.9591.0455.4690.4571.5991.7792.8171.58
MDA-MB-231−10.012.43−5.30−3.65ND14.9278.27ND35.3642.049.0867.5216.6320.7324.9332.73
HS 578T−6.13−0.32−2.253.27ND21.8885.17ND17.5837.363.6762.6617.9428.8527.0334.99
BT-549−5.74−6.03−2.68−1.2341.2313.1521.2033.8153.2065.7114.5095.71NDNDNDND
T-47D0.853.794.114.3254.9060.7431.0137.9551.0266.7223.2368.9565.1264.3561.0067.13
MDA-MB-4680.751.56−0.961.5261.7328.2249.5760.7851.2458.5112.4985.7191.1261.8244.8565.63
Cells shaded gray for non-determined activity; cells shaded orange for growth inhibition >60%; cells shaded pink for growth inhibition >90%.
Table 3. An overview of 10 µM single dose screening NCI-60 anticancer assay.
Table 3. An overview of 10 µM single dose screening NCI-60 anticancer assay.
CompoundMean Growth PercentRange of Growth Inhibition
on All Cells
The Most Sensitive Cell LinePositive Cytostatic Effect bPositive
Cytotoxic Effect c
Sensitive/Total Cell Lines
5a1.0016.17 to 15.76-0/590/590/59
5b1.6210.84 to 9.53-0/570/570/57
5c0.1012.99 to 23.92-0/590/590/59
5d0.1353.98 to 38.42-0/590/590/59
5e33.3714.73 a to 68.98HOP-9211/560/5611/56
5f13.7420.25 a to 68.87NCI-H5223/590/593/59
5g38.6318.19 a to 101.81OVCAR-417/591/5918/59
5h30.7711.54 a to 75.83OVCAR-411/560/5611/56
5i34.9527.57 a to 87.44KM1217/560/5617/56
5j53.9529.10 a to 112.65SK-MEL-532/561/5633/56
5k14.4937.79 a to 66.11SR5/560/565/56
5l72.818.89 to 152.57SK-MEL-544/564/5648/56
5m38.0433.20 a to 91.11MDA-MB-46814/590/5914/59
5n46.0227.62 a to 93.11SK-MEL-526/590/5926/59
5o47.1914.36 a to 92.81MCF728/590/5928/59
5p41.2729.37 a to 94.23SF-53921/590/5921/59
[a] Negative (−) indicates the cell killed. [b] The ratio between number of cell lines with percent growth from 0 to 50 and total number of cell lines. [c] The ratio between number of cell lines with percent growth of <0 and total number of cell lines.
Table 4. NCI in vitro testing results of GI50, TGI and LC50 of compound 5l.
Table 4. NCI in vitro testing results of GI50, TGI and LC50 of compound 5l.
Subpanel/Tumor Cell LinesCompound
5l
GI50 (µM)TGI (µM)LC50 (µM)
  Leukemia
CCRF-CEM3.72>100>100
HL60(TB)4.08>100>100
K-5623.08>100>100
MOLT-42.99ND>100
RPMI-82263.09>100>100
Non-small cell lung cancer
A549/ATCC15.50>100>100
EKVX6.02>100>100
HOP-6212.9>100>100
HOP-926.5162.70>100
NCI-H2269.0070.00>100
NCI-H235.93>100>100
NCI-H322M23.20>100>100
NCI-H46014.00>100>100
NCI-H52212.3082.20>100
Colon cancer
COLO 2059.62>100>100
HCC-299813.70>100>100
HCT-1166.97>100>100
HCT-153.4625.20>100
HT296.4537.80>100
KM125.30>100>100
SW-6202.76>100>100
CNS cancer
SF-26810.8088.10>100
SF-29510.6039.90>100
SF-53913.60>100>100
SNB-1912.40>100>100
U2517.5539.40>100
Melanoma
LOX IMVI4.7024.8096.60
MELME-3M18.9075.70>100
M1429.00>100>100
MDA-MB-4354.4476.90>100
SK-MEL-29.6765.90>100
SK-MEL-2812.3028.5065.60
SK-MEL-53.4614..2071.50
UACC-2578.0643.80>100
UACC-627.7043.50>100
Ovarian cancer
IGROV110.703.32>100
OVCAR-357.702.26>100
OVCAR-525.00>100>100
OVCAR-87.9645.20>100
NCI/ADR-RES3.4523.60>100
SK-OV-371.80>100>100
Renal cancer
786-026.60>100>100
A49818.0085.20>100
ACHN10.8075.30>100
CAKI-14.4426.90>100
RXF3939.0790.20>100
SN12C9.7466.60>100
TK-1077.10>100>100
UO-312.70>100>100
Prostate cancer
PC-36.37>100>100
DU-14513.2054.90>100
Breast cancer
MCF72.9030.90>100
MDA-MB-231/ATCC15.80>100>100
HS 578T12.8098.10>100
BT-54910.4046.30>100
T-47D3.18>100>100
MDA-MB-4683.7643.00>100
Table 5. Median growth inhibitory concentrations a (GI50, µM) of in vitro subpanel tumor cell lines for compound 5l.
Table 5. Median growth inhibitory concentrations a (GI50, µM) of in vitro subpanel tumor cell lines for compound 5l.
Subpanel Tumor Cell Line5l
MG-MID aSelectivity Index
Leukemia3.393.61
NSCL cancer11.701.05
Colon cancer5.972.05
CNS cancer10.991.11
Melanoma10.911.12
Ovarian cancer29.430.42
Renal cancer19.800.62
Prostate cancer9.781.25
Breast cancer8.141.50
Full panel MG-MID12.23
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Soudi, A.; Bender, O.; Celik, I.; El-Hafeez, A.A.A.; Dogan, R.; Atalay, A.; Elkaeed, E.B.; Alsfouk, A.A.; Abdelhafez, E.M.N.; Aly, O.M.; et al. Discovery and Anticancer Screening of Novel Oxindole-Based Derivative Bearing Pyridyl Group as Potent and Selective Dual FLT3/CDK2 Kinase Inhibitor. Pharmaceuticals 2024, 17, 659. https://doi.org/10.3390/ph17050659

AMA Style

Soudi A, Bender O, Celik I, El-Hafeez AAA, Dogan R, Atalay A, Elkaeed EB, Alsfouk AA, Abdelhafez EMN, Aly OM, et al. Discovery and Anticancer Screening of Novel Oxindole-Based Derivative Bearing Pyridyl Group as Potent and Selective Dual FLT3/CDK2 Kinase Inhibitor. Pharmaceuticals. 2024; 17(5):659. https://doi.org/10.3390/ph17050659

Chicago/Turabian Style

Soudi, Aya, Onur Bender, Ismail Celik, Amer Ali Abd El-Hafeez, Rumeysa Dogan, Arzu Atalay, Eslam B. Elkaeed, Aisha A. Alsfouk, Elshimaa M. N. Abdelhafez, Omar M. Aly, and et al. 2024. "Discovery and Anticancer Screening of Novel Oxindole-Based Derivative Bearing Pyridyl Group as Potent and Selective Dual FLT3/CDK2 Kinase Inhibitor" Pharmaceuticals 17, no. 5: 659. https://doi.org/10.3390/ph17050659

APA Style

Soudi, A., Bender, O., Celik, I., El-Hafeez, A. A. A., Dogan, R., Atalay, A., Elkaeed, E. B., Alsfouk, A. A., Abdelhafez, E. M. N., Aly, O. M., Sippl, W., & Ali, T. F. S. (2024). Discovery and Anticancer Screening of Novel Oxindole-Based Derivative Bearing Pyridyl Group as Potent and Selective Dual FLT3/CDK2 Kinase Inhibitor. Pharmaceuticals, 17(5), 659. https://doi.org/10.3390/ph17050659

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop