Ting-Ting Huang1
Xin Wang1
Shao-Jia Qiang1
Zhen-Nan Zhao1
Zhuo-Xun Wu2
Charles R. Ashby Jr.2
Jia-Zhong Li1*
Zhe-Sheng Chen2*- 1School of Pharmacy, Lanzhou University, Lanzhou, China
- 2College of Pharmacy and Health Sciences, St. John’s University, Queens, NY, United States
Chronic myelogenous leukemia (CML) typically results from a reciprocal translocation between chromosomes 9 and 22 to produce the bcr-abl oncogene that when translated, yields the p210 BCR-ABL protein in more than 90% of all CML patients. This protein has constitutive tyrosine kinase activity that activates numerous downstream pathways that ultimately produces uncontrolled myeloid proliferation. Although the use of the BCR-ABL tyrosine kinase inhibitors (TKIs), such as imatinib, nilotinib, dasatinib, bosutinib, and ponatinib have increased the overall survival of CML patients, their use is limited by drug resistance and severe adverse effects. Therefore, there is the need to develop novel compounds that can overcome these problems that limit the use of these drugs. Therefore, in this study, we sought to find novel compounds using Hypogen and Hiphip pharmacophore models based on the structures of clinically approved BCR-ABL TKIs. We also used optimal pharmacophore models such as three-dimensional queries to screen the ZINC database to search for potential BCR-ABL inhibitors. The hit compounds were further screened using Lipinski’s rule of five, ADMET and molecular docking, and the efficacy of the hit compounds was evaluated. Our in vitro results indicated that compound ZINC21710815 significantly inhibited the proliferation of K562, BaF3/WT, and BaF3/T315I leukemia cells by inducing cell cycle arrest. The compound ZINC21710815 decreased the expression of p-BCR-ABL, STAT5, and Crkl and produced apoptosis and autophagy. Our results suggest that ZINC21710815 may be a potential BCR-ABL inhibitor that should undergo in vivo evaluation.
Introduction
Chronic myelogenous leukemia (CML) is a hematopoietic neoplastic disease that primarily results from the reciprocal translocation between chromosomes 9 and 22 (t9;22) (q34; q11) in a hematopoietic stem cell (HSC), resulting in the formation of the bcr-abl oncogene (Deininger et al., 2005; Ren, 2005; Kimura, 2006). Typically, the most common BCR-ABL proteins produced in CML patients, based on the loci of the chromosomal breaks, are p210, p190, and p230 (Winter et al., 1999; Arana-Trejo et al., 2002). The p210 protein occurs was observed in more than 90% of CML patients and has been hypothesized to be the initial step in the pathogenesis of CML (Arana-Trejo et al., 2002). The p210 protein is a constitutively active, non-receptor tyrosine kinase that phosphorylates its Tyr177 residue, which subsequently activates a number of downstream signaling pathways that causes HSC cells undergo significant abnormal proliferation and differentiation compared to normal cells, producing the pathogenic changes reported in CML patients (Deininger et al., 2000; Gora-Tybor and Robak, 2008; Cao et al., 2015; Wang et al., 2016).
Imatinib (Gleevec), which selectively targets the ATP binding site of the BCR-ABL protein, is a first generation BCR-ABL tyrosine kinase inhibitor (TKI) that was approved in 2001 by the FDA for treating CML (Gontarewicz et al., 2008; Wu et al., 2014). Initially, it was reported that imatinib significantly increased the survival time of CML patients (Machova Polakova et al., 2015; Deng et al., 2020). However, subsequent clinical trials indicated that patients with advanced-stage CML experience relapse due to development of resistance caused by point mutations within the BCR-ABL domain (notably T315I), amplification of the BCR-ABL protein, overexpression of BCR-ABL protein or the ABC transporter, P-glycoprotein (Coutre et al., 2000; von Bubnoff et al., 2002; Apperley, 2007; Tanaka and Kimura, 2008).
The second generation TKIs, dasatinib, nilotinib, and bosutinib, have been shown to have efficacy in treating CML patients with certain BCR-ABL mutations, but not those with the T315I mutation (O’Hare et al., 2005; Abbas and Hsyu, 2016). The mutations of BCR-ABL are a recognized cause of resistance to second-generation, especially T315I (Baccarani et al., 2009; Khorashad et al., 2013). Ponatinib, a third-generation TKI, is also efficacious in a number of mutations in the BCR-ABL protein, including T315I (Okabe et al., 2014). Compound mutations of BCR-ABL (T315I/F359V or Y253H/T315I) are usually resistant against ponatinib therapy (Zabriskie et al., 2014). However, ponatinib has been reported to significant adverse cardiovascular effects, notably serious venous thromboembolic events (Bu et al., 2014). Therefore, there is still the need to develop highly efficacious and tolerable BCR-ABL inhibitors for the treatment of CML.
Hypogen pharmacophore model were built based on known inhibitors, the structures and activity data of known active compounds were analyzed to generate quantitative 3D pharmacophore models with common pharmacodynamic action models of the compounds, which can be used to predict the efficacy of unknown compounds (Rampogu et al., 2018; Musoev et al., 2019). The Hiphop pharmacophore models were qualitative pharmacophore models generated according to the common characteristics of known inhibitors without considering the activity data of compounds (Zheng et al., 2009). Pharmacophore models can provide guidances for drug design and virtual screening.
In this study, the quantitative Hypogen pharmacophore model and the qualitative Hiphop model were built, these models were then used to screen the ZINC database to identify potential BCR-ABL inhibitors. The identified hit compounds were further evaluated using Lipinski’s rule of five, ADMET and CDOCKER docking. The molecules obtained from the Hypogen and Hiphop pharmacophore models mapping were selected to determine their anti-proliferative efficacy in specific leukemia cell lines. Subsequently, the skeleton of the most active compound was chosen to identify additional derivatives from the database to perform structural optimization. Finally, our results led us to evaluate the efficacy of the compound ZINC217108151 to: (1) inhibit the proliferation of certain leukemia cell lines; (2) inhibit the phosphorylation of the p210 BCR-ABL tyrosine; and (3) affect the activity of the downstream proteins, STAT5 and Crkl, which are activated by BCR-ABL.
Materials and Methods
Virtual Screening
The best Hypogen and Hiphop pharmacophore models (the method of building pharmacophore shown in Supplementary Material) were used as 3D structures to search the ZINC clean drug–like database (1,469,373 molecules) to obtain the compounds that matched the chemical features (Bhadauriya et al., 2013; Kumar et al., 2015). The workflow of each screening step and the number of compounds is shown in detail in Figure 1. Compounds that satisfied the fit value of Hypogen and Hiphop were retained. All hit compounds were evaluated using Lipinski’s rules of five and absorption, distribution, metabolism, excretion, and toxicity (ADMET) screening process (Rampogu et al., 2018). Lipinski’s rule of five was predicted by using calculate molecular properties in small molecules protocol in DS2.5. Lipinski’s rule of five consists of molecular weight ≤500 Daltons, rotatable bonds number ≤10, hydrogen bond acceptors ≤10, and hydrogen bond donors ≤5 (Rampogu et al., 2017). ADMET was predicted by using ADMET descriptors in calculate molecular properties protocol in DS2.5. The cut-off value of the solubility, absorption, and the BBB were 3, 0, and 3, respectively (Rampogu et al., 2018). Finally, the database was further screened using CDOCKER. Docking analysis was used to determine whether small molecules can enter the active site of the protein (Rampogu et al., 2020). The three-dimensional structure of ABL tyrosine kinase receptor (PDB: 1IEP) was obtained from Protein Data Bank (Protein Data Bank, 1971; Machova Polakova et al., 2015). All ligands were minimized by Chemistry at Harvard Macromolecular Mechanics (CHARMm) force field and the protein was prepared by removing all of the water molecules, adding hydrogen atoms and filling the loop area (Al-Balas et al., 2013). The active site was defined as 12.9 Å around the co-crystal ligand. The prepared protein and ligands were imported to perform CDOCKER docking. Higher docking energy and interaction energy represent the favorable binding between the protein and the ligands (Rampogu et al., 2018).
Figure 1. The flow chart for the virtual screening process.
Reagents and Cell Culture
After virtual screening, the most optimal compounds were purchased from TargetMol company (Boston, MA, United States) and dissolved in DMSO to make a stock solution (1 × 10–1 mol/L) and stored at −20°C for future use.
K562, BaF3/WT, and BaF3/T315I leukemia cells lines were grown in RPMI-1640 medium supplemented with 10% fetal bovine serum (Sijiqing, China) and 1% penicillin and streptomycin. Normal human embryonic liver diploid cells (CCC-HEL-1) were grown in Dulbecco’s modified Eagle’s medium (DMEM, HyClone) supplemented with 10% fetal bovine serum (GIBCO, Australia) and 1% penicillin and streptomycin, all cells were maintained in humidified air at 37°C with 5% CO2 (Gupta et al., 2020).
Antibodies C-Abl and p-C-Abl were purchased from Cell Signaling Technology (MA, United States), p-Crkl was purchased from Affinity Biosciences (OH, United States), p-STAT5A and β-actin were purchased from Bioss (MA, United States), LC3 was purchased from Sigma-Aldrich (United States), Beclin1 was purchased from Abcam (Cambridge, United Kingdom), Caspase3 was purchased from Cell Signaling Technology (MA, United States). HRP-conjugated Affinipure Goat Anti-Rabbit antibody was purchased from Proteintech (IL, United States).
3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium Assay
K562, BaF3/WT, and BaF3/T315I leukemia cells were seeded in 96-well plates at a density of 5 × 103 cells/well and 50 μL of various drug concentrations were added. CCC-HEL-1 cells were grown in 96-well plates at the same density and exposed to the same drug concentrations. The plates were incubated at 37°C for 72 h. Cell proliferation was evaluated using MTT [3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium] assay. Fifteen μL of the MTT solution was added to each well of 96-well and plates were continually incubated at 37°C for 4 h and three liquid solutions. SDS, HCl, and isobutanol were added to each well and the plates were incubated at 37°C overnight. The optical density (OD) of each sample was measured at a wavelength of 570 nm and the 50% inhibitory concentration (IC50) was determined by from the concentration–response curve.
Cell Cycle and Cell Apoptosis Analysis
K562 cells were plated in 6-well plates with 1 × 106 cells/well, BaF3/WT and BaF3/T315I cells were plated in 6-well plates with 5 × 105 cells/well. After all of the cells were incubated with 0 (vehicle), 0.1, 1, or 10 μM concentration of the test compounds for 24 h, they were harvested and washed with PBS and 1 ml of 70% cold ethanol was added to the cells (Kim et al., 2015; Jiang et al., 2016). The ethanol-fixed cells were stored at −20°C overnight. Subsequently, the cells were washed with PBS and centrifuged at 1,800 rpm for 8 min to obtain the pellet. The pellet was resuspended in 100 μL of RNase A, incubated at 37°C for 30 min, and mixed with 400 μL of propidium iodide (PI) and incubated at 4°C in the dark for 30 min. DNA content was determined using flow cytometry, at an emission wavelength of 488 nm and analyzed using ModFit LT software.
All cells were plated in 6-well plates with 5 × 105 cells/well and incubated with 0 (vehicle), 0.1, 1, or 10 μM of the test compounds for 24 h. The cells were harvested, washed twice with cold PBS, resuspended with binding buffer, stained using Annexin V- fluorescein isothiocyante (FITC) and PI and mixed with PBS to obtain a final volume of 500 μL. Cellular apoptosis was determined using flow cytometry method within 1 h of staining and dot plots were set up to detect FITC and PI. Cells stained with FITC and without PI (FITC+, PI-) were categorized as being in early apoptosis, and cells stained with both dyes (FITC+, PI+) were categorized as being in late apoptosis or necrotic (Peng et al., 2001).
Western Blot Analysis
Cells were plated in 6-well plates with 5 × 105 cells/well and incubated with 0 (vehicle), 0.1, or 1 μM of the test compounds for 48 h. The harvested cells were washed with PBS, and lysed in a RIPA lysis buffer and collected by centrifugation at 14,000 × g for 15 min at 4°C. The cell extracts were separated by SDS-polyacrylamide gel electrophoresis (PAGE) and transferred to PVDF membranes. The membranes were blocked with 5% non-fat milk at room temperature for 75 min. The membranes were incubated with C-Abl, p-C-Abl, p-Crkl, p-STAT5A, β-actin, LC3, Beclin1 or caspase-3 antibodies overnight at 4°C, then washed three times with TBST and incubated for 1 h with the HRP-conjugated Affinipure Goat Anti-Rabbit antibody at room temperature. The protein expressions were quantified by ImageJ software.
Results and Discussion
Hypogen Pharmacophore Model
Construction of Hypogen Model
The training set samples for the Hypogen modeling included 21 diverse BCR-ABL inhibitors obtained from different publications (Buchdunger et al., 1996; Nagar et al., 2002; Gumireddy et al., 2005; Barnes et al., 2007; Katsoulas et al., 2008; Azam et al., 2010; Ren et al., 2013; Bu et al., 2014; Pan et al., 2015; Zabriskie et al., 2015; Halbach et al., 2016; Salerno et al., 2017; Rossari et al., 2018) with diverse structures and broad efficacy values (IC50) that ranged from 0.0002 μM to 9.8 μM. The chemicals used to construct the Hypogen pharmacophore models are shown in Supplementary Figure 1. The top 10 pharmacophore models containing HBA, AR and H chemical features are shown in Supplementary Table 1. Among all of these models, Hypogen1, composed of HBA and 4H chemical features, had the lowest total cost (91.20), highest cost difference (69.50), lowest RMS (0.50), best correlation coefficient (0.99), and highest fit value (11.89). Cost analysis was applied to determine the statistical significance of the Hypogen models (Ye et al., 2010). The cost difference of Hypogen1 (69.4973), the difference between the null cost and total cost indicated more than 90 instances of statistical significance of achieving an acceptable pharmacophore (Kavitha et al., 2015; Kumar et al., 2015). The low RMS and high correlation coefficient suggested that the model is stable and internally predictive (Kandakatla and Ramakrishnan, 2014). The fit value was indicative of the overall fitness of the training set compounds on a pharmacophore model during pharmacophore generation (Dube et al., 2012). Therefore, the Hypogen1 was selected based on the highest correlation coefficient, lowest total cost and RMS value (Aparoy et al., 2010), 3D spatial relationship and distance constraints of Hypogen1 are shown in Figure 2A. The Figures 2B,C show the most active and inactive compounds in the training set aligned with Hypogen1.
Figure 2. (A) 3D spatial relationship and geometric parameters of Hypogen1. (B) Hypogen1 mapping with the most active compound. (C) Hypogen 1 mapping with the least active compound.
If the compounds were divided into three categories based on their efficacies as being highly active (IC50 ≤ 0.1 μM, +++), moderately active (0.1 μM < IC50 < 3 μM, ++), and least active (IC50 ≥ 3 μM, +), this categorization could be used to rapidly estimate the predictive accuracy of the pharmacophore (Huang et al., 2012). Table 1 shows that the efficacies of the training set compounds could be accurately predicted.
Table 1. Actual and estimated efficacies of the training set molecules based on the Hypogen1 model.
Validation of Hypogen1
The test set, Fischer randomization and decoy set were used to verify the Hypogen pharmacophore model.
Test set prediction
Thirty eight compounds, with different scaffold and efficacy ranges, were used as a test set (Supplementary Figure 2) to assess whether the pharmacophore model can accurately predict the efficacy of the compound (Huang et al., 2012). The experimental and predicted efficacies of the test set are shown in Table 2. A significant correlation coefficient of 0.847 was obtained for the test set compounds using regression analysis. The error factors of all compounds were less than 5%. Figure 3 shows the scatter plot of the experimental and the predicted efficacies, indicating that the samples are clustered near the diagonal line, y = x. These results indicated that the Hypogen1 model has excellent predictive validity.
Table 2. Actual and estimated efficacy of the test set molecules based on the Hypogen1 model.
Figure 3. Correlation between the experimented and predicted efficacies of the training set and the test set based on the Hypogen1 model.
Fischer randomization test
The Fischer randomization test is a cross validation method that can be used to evaluate the statistical correlation between structures and activities in training set compounds (Pal et al., 2019). To achieve a confidence level of 95% for Fischer’s randomization, 19 random hypotheses were generated (Sakkiah et al., 2012). Figure 4 clearly shows that the Hypogen1 pharmacophore model was not generated by chance as the total cost of these randomly generated hypotheses is much higher than the initial hypothesis.
Figure 4. The difference in cost between the initial hypothesis and 19 random hypotheses. The 95% confidence level was selected to validate the Fischer randomization test results (Hypogen1).
Decoy set
A decoy set was used to validate the efficiency of Hypogen1 by computing various parameters such as false positive, false negative, enrichment factor (EF), goodness of fit score (GF), total number of hit molecules (Ht), and total number of active compounds (Ha) (Halgren et al., 2004). The decoy set database (D) contained 325 compounds, including 25 BCR-ABL inhibitors and 300 inactive compounds. Our results indicated that 26 compounds were identified as being active based on the Hypogen1 model and 21 compounds were known inhibitors of the BCR-ABL protein. The relevant parameters for the decoy set are shown in Table 3. An EF of 10.5 and goodness of GF of 0.8022 verifies the high efficiency of Hypogen1, indicating the high predictive ability of the model.
Table 3. Statistical parameters of the decoy set for Hypogen1.
Hiphop Pharmacophore Model
Construction of Hiphop Model
Eight compounds were selected as the training set to generate qualitative Hiphop models for common feature pharmacophore generation. The structures and bioactivities of the training set are shown in Supplementary Figure 3. The top 10 common feature hypotheses were produced and are shown in Supplementary Table 2, according to their ranking scores. Based on the ranking scores and feature similarities of hypotheses, ten hypotheses were clustered into two groups: the first group models contained 2RA, 3H, 1HBD, and 2HBA chemical features and the second group models contained 1RA, 4H, 1HBD, and 2HBA chemical features.
Hiphop1 and Hiphop9 were the top ranked models for each cluster and therefore, were selected for further analysis. The bestfit values of Hiphop1 and Hiphop9 model with all training set compounds are shown in Supplementary Table 3, and the average of the bestfit value for active and inactive compounds were calculated for each hypothesis. Based on the results in Supplementary Table 3, we can see that compared with Hiphop9, Hiphop1 has a more pronounced difference in the active and inactive compounds, and thus, Hiphop1 was chosen as the optimal Hiphop model (Figure 5A). Figure 5B shows the most active compound in the training set aligned with Hiphip1.
Figure 5. (A) The Hiphop1 pharmacophore model. (B) Hiphop1 mapped with the most active compound.
Validation of Hiphop1
To validate the Hiphop1 model, 21 active compounds, with different scaffold and efficacies (0.0005–41 μM) (Supplementary Figure 4) and 220 compounds with unknown efficacies from Maybridge, were combined to create a new database (D1). Hiphop1 was used to select the active compounds. As a result, 22 compounds were screened from the database by Hiphop1 and 21 of these compounds have been previously proven to be known effective compounds. The validation results (Table 4) indicated an EF of 10.95 and a goodness of GF of 0.91, thus validating the high efficiency of Hiphop1 and demonstrating that Hiphop1 can discriminate active compounds from inactive compounds.
Table 4. Statistical parameters of Hiphop-1.
Virtual Screening
The best validated Hypogen1 and Hiphop1 pharmacophore models were used as 3D queries to screen the drug-like subset of the ZINC Database. ZINC is a free database of purchasable molecules, including multiple subsets such as drug-like and lead-like and it offers similarity searching on its website (Irwin and Shoichet, 2005; Awale et al., 2015). The work flow for the virtual screening is shown in Figure 1. Overall, 31,966 compounds matched all of the chemical features of Hypogen1 and 4,188 compounds had fit values > 10. By using Hiphop1, 896 compounds were selected out for their high magnitude of matching with all of the chemical features of Hiphop and 192 compounds with fit values > 1.8 were selected. The hit compounds were further subjected to Lipinski’s rule of five and ADMET to assess their drug-like and pharmacokinetic properties (Rampogu et al., 2018), and 3,865 and 162 compounds meet these standards, respectively. Subsequently, molecular docking was used to determine if these compounds interact with the BCR-ABL protein. In total, 3,856 and 161 compounds were screened using Lipinski’s rule of five and ADMET, respectively, interacted with the active site of protein. Finally, identical compounds screened from the two pharmacophore models were selected for the in vitro efficacy assay. The docking scores for these compounds are shown in Supplementary Table 4.
Compounds Inhibit the Proliferation of K562 Cells
There were 13 identical compounds (Supplementary Table 4) and we purchased six of these for the proliferation experiments: ZINC36617838, ZINC30201139, ZINC65008391, ZINC45895251, ZINC36617852, and ZINC36617849.
MTT assays were performed to determine the inhibitory efficacy of compounds on the proliferation of K562 leukemia cells. The cells were incubated with various concentrations of the six compounds for 72 h and the results are shown in Figure 6. The compound, ZINC36617838, had a greater efficacy in inhibiting the proliferation of K562 leukemia cells (IC50 10.002 μM) compared to the other ZINC compounds. However, the efficacy of ZINC36617838 was approximately 200 times lower than that of imatinib (IC50 0.047 μM). This finding led us to select ZINC36617838 as a potential hit compound for structural optimization. The structure of ZINC36617838 and the interaction between ZINC36617838 and 1IEP is shown in Supplementary Figure 5.
Figure 6. Cytotoxic efficacy of compounds in K562 leukemia cells. Cells were incubated with various concentrations of IM, ZINC36617838, ZINC30201139, ZINC65008391, ZINC45895251, ZINC36617852, or ZINC36617849 for 72 h. Cell viability was determined using the MTT assay. The results are represented as the mean ± SD.
Structure Optimization
In order to find new compounds with higher inhibitory efficacy, the skeleton structure of ZINC36617838 (Figure 7), was used as a substructure to search for its derivatives in the ZINC database. The compounds that had 40% similarity with ZINC36617838 were selected (circled in red, green, and blue in Figure 7). As shown in Figure 7, the analysis identified 290 compounds. The validated Hypogen1 and Hiphop1 pharmacophores were used as 3D queries to screen the database and four and 153 compounds matched all of the pharmacophore features of Hiphop1 and Hypogen1, respectively. Subsequently, four and 148 compounds meet the standards of the Lipinski’s rules of five and ADMET filters, respectively, and were docked into active site of protein to remove the detected false positives (Rampogu et al., 2018).
Figure 7. The work flow chart of structure optimization and its screening process.
After screening, four compounds from the Hiphop screening were also identified from the Hypogen screening. However, the predicted IC50 values or efficacy of 29 out of the 148 compounds were less than 0.1 μM in the Hypogen pharmacophore screening pathway, and the CDOCKER INTERACTION ENERGY of six compounds were greater than 60. Therefore, these ten compounds were selected for further experimental evaluation. We purchased ZINC21710815, ZINC36617889, and ZINC20617585. The binding modes of these three compounds are shown in Figures 8A–C.
Figure 8. The binding modes of three compounds and imatinib with the BCR-ABL protein (PDB: 1IEP): hydrogen bond (green), pi-pi (yellow). (A) ZINC20617585; (B) ZINC36617889; (C) ZINC21710815; and (D) Imatinib (IM). The effect of these compounds on the proliferation of (E) K562, (F) BaF3, (G) BaF3/T315I leukemia cells, and (H) CCC-HEL-1 cells. The cells were incubated with various concentrations of imatinib (IM), ZINC21710815, ZINC36617889, or ZINC20617585 for 72 h. The viability of the cells was assessed using the MTT assay. The results are represented as the mean ± S.D.
We used the MTT assay to determine the efficacy of these compounds to inhibit the proliferation of K562, BaF3/WT, and BaF3/T315I leukemia cells (Figure 8). All cells were incubated with different concentrations of each compound for 72 h. As shown in Figure 8E, imatinib and ZINC21710815 inhibited the proliferation of K562 cells with IC50 values of 0.047 and 0.531 μM, respectively. The IC50 value for ZINC21710815 was approximately 10 times lower than that of imatinib. Imatinib and ZINC21710815 decreased the viability of BaF3/WT leukemia cells, with IC50 values of 0.196 and 0.512 μM, respectively (Figure 8F). In the BaF3/T315I leukemia cells, which have the T315I mutation that produces resistance to imatinib, ZINC21710815 significantly inhibited the proliferation of BaF3/T315I leukemia cells, (IC50 = 0.88 μM), and its IC50 value was significantly lower than that for imatinib (IC50 = 10.11 μM), and its IC50 for normal CCC-HEL-1 cells was 89.587 μM (Figure 8H). The selectivity indexes (SI) were 1, 168.71, 174.97, and 101.80, respectively, for the CCC-DEL-1, K562, BaF3/WT, and BaF3/T315I cells. These results suggest that ZINC21710815 was not cytotoxic in normal CCC-HEL-1 cells and it significantly inhibited the proliferation of wild-type BCR-ABL and T315I mutated BCR-ABL leukemia cells.
ZINC21710815 Inhibits the Growth of CML Cells
Based on the proliferation experiments, we next determined the effect of ZINC21710815 on the cell cycle of the leukemia cells. The K562, BaF3/WT, and BaF3/T315I leukemia cells were incubated with various concentrations of ZINC21710815 for 24 h. These cells were incubated and stained with PI, which stains DNA, allowing for the analysis of cellular DNA content using flow cytometry. ZINC21710815 significantly increased the accumulation of K562 leukemia cells in the G2 phase, suggesting an arrest of the cells in the G2 phase (Figures 9A,D). In contrast, ZINC217108155 significantly increased the accumulation of BaF3/WT and BaF3/T315I cells leukemia cells in the G1 phase (Figures 9B,C,D). These findings suggested that ZINC21710815 may inhibit the proliferation of leukemia cells by interrupting the cell cycle in CML cells.
Figure 9. The effect of ZINC21710815 on the cell cycle in K562, BaF3/WT and BaF3/T315I leukemia cells. (A) K562 leukemia cells were incubated with various concentrations of ZINC21710815 for 24 h at a density of 1 × 106 cells/well. (B) BaF3 leukemia cells were incubated with various concentrations of ZINC21710815 for 24 h at a density of 5 × 105 cells/well. (C) BaF3/T315I leukemia cells were incubated with various concentrations of ZINC21710815 for 24 h at a density of 5 × 105 cells/well. (D) The relative percentages of the cell cycles for K562, BaF3/WT, and BaF3/T315I leukemia cells were calculated using ModFit LT software. The results are represented as the mean ± SD.
ZINC 21710815 Induces Apoptosis in CML Cells
We conducted experiments to determine if ZINC21710815 inhibits the proliferation of K562, BaF3/WT, and BaF3/T315I leukemia cells by inducing apoptosis. The cells were cultured with various concentrations of ZINC21710815 for 24 h and the cells were stained using Annexin V-FITC/PI and analyzed using flow cytometry. As shown in Figure 10A, ZINC21710815 produce a concentration-dependent increase in the apoptosis of K562, BaF3/WT, and BaF3/T315I leukemia cells. Apoptosis is regulated by specific caspases, and cleaved caspase-3 has been shown to be a reliable marker for detecting apoptotic cells (Crowley and Waterhouse, 2016). Therefore, we measured the expression levels of caspase-3 using western blotting (Figure 10B). ZINC21710815 significantly increased the levels of cleaved caspase-3 in BaF3/WT and BaF3/T315I leukemia cells, suggesting that this compound induces apoptosis.
Figure 10. ZINC21710815 induces apoptosis and autophagy in CML cells. (A) K562, BaF3/WT, and BaF3/T315I cells were incubated with various concentration of ZINC21710815 for 24 h at a density of 5 × 105 cells/well, stained using Annexin V-FITC and PT and analyzed by FlowJo. (B) Caspase-3 protein expression in BaF3/WT and BaF3/T315I cells incubated with 0 (vehicle), 0.1, or 1 μM of ZINC21710815 for 48 h. LC3 (C) and Beclin1 (D) protein expression in BaF3/WT and BaF3/T315I cells incubated with 0 (vehicle), 0.1, or 1 μM of ZINC21710815 for 48 h.
ZINC21710815 Induces Autophagy in BaF3/WT and BaF3/T315I Leukemia Cells
To determine whether ZINC21710815 induced autophagy, we determined the expression of the autophagy-related proteins, LC3-II, and Beclin1. LC3-II is a hallmark protein of autophagy and the level of LC3-II expression is significantly correlated with the number of autophagosomes (Mizushima and Yoshimori, 2007; Yoshii and Mizushima, 2017; Lu et al., 2018). In the early stages of autophagosome formation, Beclin1 is considered to be an essential component for the initiation of autophagy (Hamurcu et al., 2018). As shown in Figures 10C,D, the level of LC3-II/LC3-I and Beclin1 are significantly increased in BaF3/WT and BaF3/T315I leukemia cells following incubation with ZINC21710815 for 48 h. These data suggest that ZINC21710815 induces autophagy, in part, by increasing the expression of the levels LC3-II and Beclin1 in BaF3/WT and BaF3/T315I leukemia cells.
ZINC21710815 Inhibits Tyrosine Phosphorylation of the BCR-ABL Protein and Its Downstream Protein Targets, Signal Transducer and Activator of Transcription 5 (STAT5) and CRK Like Proto-Oncogene, Adaptor Protein (Crkl)
BCR-ABL activates multiple downstream signaling pathways, including STAT5 and p-Crkl (Gupta et al., 2016). The activation of the transcription factor STAT5 is crucial for the progression of CML (Warsch et al., 2012). STAT5 was necessary for cellular proliferation and overexpression of constitutively active STAT5 could stimulate cell proliferation (de Groot et al., 1999). Crkl is an adaptor protein and a major substrate of BCR-ABL in CML cells (Oda et al., 1996). Crkl is phosphorylated by BCR-ABL and it subsequently activates other pathways that play a role in leukemic cell transformation (Senechal et al., 1996). To confirm the connection between the compound and BCR-ABL and downstream proteins, the effect of ZINC21710815 on the tyrosine phosphorylation of BCR-ABL and its downstream protein, STAT5 and Crkl, was determined using western blotting. BaF3/WT leukemia cells were incubated with various concentrations of imatinib and ZINC21710815 for 48 h. ZINC21710815 significantly reduced the phosphorylation of BCR-ABL at a concentration of 0.1 μM, decreased the phosphorylation of STAT5 at a concentration of 0.1 μM and reduced the phosphorylation of Crkl at a concentration of 1 μM in BaF3/WT leukemia cells (Figure 11).
Figure 11. BaF3/WT leukemia cells were incubated with 0 (vehicle), 0.1, or 1 μM for 48 h with either imatinib and ZINC21710815 the levels of BCR-ABL and phosphorylated BCR-ABL (A), and the downstream proteins p-Crkl (B), p-STAT5 (C) were analyzed using western blotting. The results are represented as the mean ± SD.
Discussion
In this study, Hypogen pharmacophore models were constructed based on 21 BCR-ABL TKIs and the best quantitative pharmacophore model consisted of five features. The correlation coefficient of Hypogen1 with the training and test set were 0.973 and 0.847, respectively. Hypogen1 was further validated by Fischer’s randomization method and by using a decoy set. Hiphop pharmacophore models were constructed based on five highly activity BCR-ABL inhibitors and three inhibitors with low efficacy. The best Hiphop pharmacophore model consists of eight features. The best pharmacophore models were selected to validate the test set and were highly efficient in distinguishing active from inactive compounds. The best Hypogen and Hiphop pharmacophore models were used as 3D queries to screen the ZINC database in order to find novel potential BCR-ABL inhibitors. The hit compounds were further screened by Lipinski’s rule of five, ADMET and molecular docking. Finally, three compounds, ZINC21710815, ZINC36617889, and ZINC20617585, were subjected to further analysis, of which compound ZINC21710815 had the greatest efficacy in inhibiting the proliferation of K562, BaF3/WT, and BaF3/T315I leukemia cells.
ZINC21710815, ZINC36617889, and ZINC20617585, as well as imatinib, interact with the BCR-ABL protein via a carbonyl of amide, which forms an O-NH hypogen interaction with Asp381, and all of these compounds form π-π interactions with residue Lys271. Given that imatinib binds to the ABL kinase domain known as the DFG-out conformation (Asp381-Phe382-Gly383) (Pereira et al., 2012), we hypothesize that ZINC21710815, ZINC36617889, and ZINC20617585 can induce the DFG-out conformation. Our results indicated that ZINC21710815 and imatinib have a higher –CDOCKER ENERGY compared to ZINC36617889 and ZINC 20617585 (Figure 8) and this could potentially explain why ZINC21710815 and imatinib have similar efficacies in inhibiting the proliferation of the leukemia cell lines used in this study. Imatinib inhibited the proliferation of K562 (IC50 = 0.047 μM), BaF3/WT (IC50 = 0.196 μM), and BaF3/T315I (IC50 = 10.11 μM) leukemia cells. ZINC21710815 was shown to be a novel compound that inhibited the proliferation of K562 (IC50 = 0.531 μM), BaF3/WT (IC50 = 0.512 μM), and BaF3/T315I leukemia cells (IC50 = 0.88 μM) in a concentration–dependent manner. ZINC21710815 and imatinib decreased the viability of K562 and BaF3/WT leukemia cells and ZINC21710815 produced a greater inhibition of BaF3/T315I leukemia cell viability compared to imatinib (Figure 8F). Although ZINC21710815 was the lead compound, it produced a lower magnitude of inhibition of BaF3/T315I leukemia cell proliferation compared to dasatinib in BaF3/WT and BaF3/T315T leukemia cells (IC50 = 1 nM/300 nM) and ponatinib in BaF3/WT and BaF3/T315I leukemia cells (IC50 = 1 nM/8 nM) (Cassuto et al., 2012), we will continue to optimize the structure of ZINC21710815 to improve its efficacy, particularly in the leukemia cells with the T315I mutation.
Our mechanistic experiments suggested that the ZINC21710815-induced decrease in the K562, BaF3/WT, and BaF3/T315I leukemia is due to it producing cell cycle arrest and apoptosis. ZINC21710815 increased the accumulation of K562 leukemia cells the G2 phase in and increased accumulation of BaF3/WT and BaF3/T315I leukemia cells in the G1 phase. As reported, imatinib induced an accumulation in the S phase, and dasatinib caused accumulation at the G1 phase and were efficacious in inhibiting leukemia cell proliferation (Song et al., 2013; Wu et al., 2014). In contrast, ZINC21710815 significantly increased the accumulation of cells in the G1 or G2 phases. These results suggest that ZINC21710815, in part, decreases the proliferation of the leukemia cell lines used in this study by producing cell cycle arrest in either the of the G1 or G2 phases. The anti-proliferative efficacy of ZINC21710815 may also be due to it inducing the apoptosis of CML cells as this compound increased the expression of the apoptosis-inducing protein, cleaved capase-3 in BaF3/WT and BaF3/T315I. It has been reported that the activation of the enzyme caspase-3 plays a critical role in inducing cancer cell apoptosis, and imatinib, ponatinib, and dasatinib can increase the levels of the caspase-3 by increasing the proteolysis of its inactive zymogen, pro-caspase-3, thereby increasing apoptosis (Quintás-Cardama et al., 2008; Okabe et al., 2014; Wu et al., 2014). Similarly, ZINC21710815, at 0.1 uM, increased the level of cleaved caspase-3, producing leukemia cell apoptosis. Finally, our results indicated that ZINC21710815 induces autophagy in BaF3/WT and BaF3/T315I leukemia cells, based on the increase in the expression of the autophagy-related proteins LC3-II and Beclin1, in BaF3/WT, and BaF3/T315I leukemia cells. Numerous studies have reported that the induction of autophagy by various TKIs, including imatinib and dasatinib, may protect CML cells from death, thereby promoting their survival (Bellodi et al., 2009; Helgason et al., 2011). Indeed, the inhibition of the BCR-ABL kinase activity is a main trigger of autophagy that allows the cells to survive in a stressful environment (Bellodi et al., 2009; Calabretta and Salomoni, 2012). In vitro studies indicate that the inhibition of autophagy increases the efficacy of certain TKIs in CML cells (Bellodi et al., 2009; Salomoni and Calabretta, 2009; Crowley et al., 2011; Yu et al., 2012). In the current study, we do not know if ZINC21710815’s induction of autophagy increases or decreases its anti-proliferative efficacy. Therefore, future experiments, where autophagy of the three leukemia cell lines is inhibited, must be done to ascertain if ZINC21710815’s efficacy is increased or decreased. However, it is important to note that the induction of autophagy in leukemia cells by imatinib increases the sequestration of the BCR-ABL protein in autophagosomes, decreasing its levels, thereby attenuating its oncogenic efficacy (Elzinga et al., 2013). Furthermore, the inhibition of autophagy in wild type and T315I leukemia cells decreases the biodegradation of the BCR-ABL protein, thus increasing its stability (Shinohara et al., 2019).
In order to further determine the efficacy ZINC21710815 on the BCR-ABL kinase, we measured the effect of ZINC21710815 on levels of the proteins STAT5 and Crkl, which are substrates for the BCR-ABL kinase. BCR-ABL kinase phosphorylates and activates the protein transcription factor, STAT5, which is translocated to the nucleus, where it increases the transcription of genes that produce proteins that promote cell survival and proliferation (Huang et al., 2002). Our results indicated that ZINC21710815 and imatinib, at 0.1 μM, significantly decreased p-STAT5 level in BaF3/WT leukemia cells. Our imatinib results are congruent with previous studies in leukemia cells (Dong et al., 2018; Tu et al., 2018) and similarly, dasitinib also decreases p-STAT5 levels in BaF3/WT leukemia cells (Fiskus et al., 2006). We also determined the effect of ZINC21710815 on the levels of phosphorylated Crkl, a protein that is phosphorylated by the BCR-ABL (Bhat et al., 1997). Phosphorylated Crkl is one of main tyrosyl-phosphoproteins present in peripheral blood cells of CML patients, where it functions as a nuclear adaptor and transcriptional activator in BCR-ABL expressing cells (Nichols et al., 1994; Rhodes et al., 2000). Furthermore, the levels of phosphorylated Crkl can be used as an indirect marker of BCR-ABL function and its levels can be used to determine the status of BCR-ABL kinase activity due to the fact that Crkl is only phosphorylated by BCR-ABL (Nichols et al., 1994; Hamilton et al., 2006). Our in vitro results indicated ZINC21710815 significantly decreased the levels of phosphorylated Crkl in BaF3/WT leukemia cell lines, the magnitude of inhibition of phosphorylated Crkl levels by ZINC21710815 was similar to that of imatinib. Similarly, dasatinib has been reported to decrease p-Crkl (Fiskus et al., 2006) levels with an efficacy comparable to that of ZINC21710815. Based on these findings, we hypothesize that ZINC21710815 may decrease leukemia cell proliferation by decreasing the phosphorylation of BCR-ABL, STAT5, and Crkl.
Our experimental results showed that the compound ZINC21710815 inhibits the proliferation of leukemia cells, inhibits cell cycle, induces autophagy and apoptosis and inhibits the tyrosine phosphorylation of BCR-ABL target protein and downstream proteins. However, ZINC21710815 had no significant effect on the expression of BCR-ABL kinase in the T315I leukemia cells. Therefore, in the future, we will conduct experiments to optimize the structure of ZINC21710815 to increase its inhibitory efficacy in T315I leukemia cells.
The BCR-ABL fusion gene can be seen in CML patients, it is also found in other types of leukemia patients, such as acute myeloid leukemia (AML) patients (Bucur et al., 2013), B-cell acute lymphoblastic leukemia (B-ALL) patients (Zhou et al., 2018), and mixed phenotype acute leukemia (MPAL) patients (Alexander et al., 2018). ZINC21710815 could decrease the expression of BCR-ABL, which may indicate that the compound may also have an inhibitory effect in other cancers.
In conclusion, ZINC21710815 significantly decreased the proliferation of K562, BaF3/WT, and BaF3/T315I leukemia cells due to it (1) producing cell cycle arrest; (2) inducing apoptosis; and (3) inhibiting the phosphorylation of BCR-ABL kinase and the downstream targets, STAT5 and Crkl, in BaF3/WT leukemia cells. As stated above, it remains to be determined if inducing autophagy increases or decreases the anti-proliferative efficacy of ZINC21710815. Overall, our results suggest that ZINC21710815 is an inhibitor of BCR-ABL tyrosine kinase activity. Additional studies must be done to determine the in vivo efficacy and toxicity of ZINC21710815.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
Author Contributions
T-TH and XW designed research. T-TH, S-JQ, and Z-XW performed research. T-TH and Z-NZ analyzed data. T-TH and CA wrote the manuscript. CA and Z-SC revised the manuscript. J-ZL and Z-SC supervised the research. All authors contributed to the article and approved the submitted version.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We would like to thank Brian J. Druker for kindly providing the plasmids of BCR-ABL (The Knight Cancer Institute, Oregon Health and Science University, Portland, OR, United States).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.649434/full#supplementary-material
References
Abbas, R., and Hsyu, P. H. (2016). Clinical pharmacokinetics and pharmacodynamics of Bosutinib. Clin. Pharmacokinet. 55, 1191–1204. doi: 10.1007/s40262-016-0391-6
Al-Balas, Q. A., Amawi, H. A., Hassan, M. A., Qandil, A. M., Almaaytah, A. M., and Mhaidat, N. M. (2013). Virtual lead identification of farnesyltransferase inhibitors based on ligand and structure-based pharmacophore techniques. Pharmaceuticals (Basel) 6, 700–715. doi: 10.3390/ph6060700
Alexander, T. B., Gu, Z., Iacobucci, I., Dickerson, K., Choi, J. K., Xu, B., et al. (2018). The genetic basis and cell of origin of mixed phenotype acute leukaemia. Nature 562, 373–379.
Aparoy, P., Kumar Reddy, K., Kalangi, S. K., Chandramohan Reddy, T., and Reddanna, P. (2010). Pharmacophore modeling and virtual screening for designing potential 5-lipoxygenase inhibitors. Bioorg. Med. Chem. Lett. 20, 1013–1018. doi: 10.1016/j.bmcl.2009.12.047
Apperley, J. F. (2007). Part I: mechanisms of resistance to imatinib in chronic myeloid leukaemia. Lancet Oncol. 8, 1018–1029. doi: 10.1016/s1470-2045(07)70342-x
Arana-Trejo, R. M., Ruíz Sánchez, E., Ignacio-Ibarra, G., Báez De La Fuente, E., Garces, O., Gómez Morales, E., et al. (2002). BCR/ABL p210, p190 and p230 fusion genes in 250 Mexican patients with chronic myeloid leukaemia (CML). Clin Lab Haematol 24, 145–150. doi: 10.1046/j.1365-2257.2002.00413.x
Awale, M., Jin, X., and Reymond, J. L. (2015). Stereoselective virtual screening of the ZINC database using atom pair 3D-fingerprints. J. Cheminform. 7:3.
Azam, M., Powers, J. T., Einhorn, W., Huang, W. S., Shakespeare, W. C., Zhu, X., et al. (2010). AP24163 inhibits the gatekeeper mutant of BCR-ABL and suppresses in vitro resistance. Chem. Biol. Drug Des. 75, 223–227. doi: 10.1111/j.1747-0285.2009.00911.x
Baccarani, M., Cortes, J., Pane, F., Niederwieser, D., Saglio, G., Apperley, J., et al. (2009). Chronic myeloid leukemia: an update of concepts and management recommendations of European LeukemiaNet. J. Clin. Oncol. 27, 6041–6051. doi: 10.1200/jco.2009.25.0779
Barnes, D. J., De, S., Van Hensbergen, P., Moravcsik, E., and Melo, J. V. (2007). Different target range and cytotoxic specificity of adaphostin and 17-allylamino-17-demethoxygeldanamycin in imatinib-resistant and sensitive cell lines. Leukemia 21, 421–426. doi: 10.1038/sj.leu.2404533
Bellodi, C., Lidonnici, M. R., Hamilton, A., Helgason, G. V., Soliera, A. R., Ronchetti, M., et al. (2009). Targeting autophagy potentiates tyrosine kinase inhibitor-induced cell death in Philadelphia chromosome-positive cells, including primary CML stem cells. J. Clin. Invest. 119, 1109–1123. doi: 10.1172/jci35660
Bhadauriya, A., Dhoke, G. V., Gangwal, R. P., Damre, M. V., and Sangamwar, A. T. (2013). Identification of dual Acetyl-CoA carboxylases 1 and 2 inhibitors by pharmacophore based virtual screening and molecular docking approach. Mol. Divers. 17, 139–149. doi: 10.1007/s11030-013-9425-2
Bhat, A., Kolibaba, K., Oda, T., Ohno-Jones, S., Heaney, C., and Druker, B. J. (1997). Interactions of CBL with BCR-ABL and CRKL in BCR-ABL-transformed myeloid cells. J. Biol. Chem. 272, 16170–16175. doi: 10.1074/jbc.272.26.16170
Bu, Q., Cui, L., Li, J., Du, X., Zou, W., Ding, K., et al. (2014). SAHA and S116836, a novel tyrosine kinase inhibitor, synergistically induce apoptosis in imatinib-resistant chronic myelogenous leukemia cells. Cancer Biol. Ther. 15, 951–962. doi: 10.4161/cbt.28931
Buchdunger, E., Zimmermann, J., Mett, H., Meyer, T., Müller, M., Druker, B. J., et al. (1996). Inhibition of the Abl protein-tyrosine kinase in vitro and in vivo by a 2-phenylaminopyrimidine derivative. Cancer Res. 56, 100–104.
Bucur, O., Stancu, A. L., Goganau, I., Petrescu, S. M., Pennarun, B., Bertomeu, T., et al. (2013). Combination of bortezomib and mitotic inhibitors down-modulate Bcr-Abl and efficiently eliminates tyrosine-kinase inhibitor sensitive and resistant Bcr-Abl-positive leukemic cells. PLoS One 8:e77390. doi: 10.1371/journal.pone.0077390
Calabretta, B., and Salomoni, P. (2012). Suppression of autophagy by BCR/ABL. Front. Biosci. (Schol Ed) 4:453–460. doi: 10.2741/278
Cao, R., Wang, Y., and Huang, N. (2015). Discovery of 2-Acylaminothiophene-3-Carboxamides as multitarget inhibitors for BCR-ABL kinase and microtubules. J. Chem. Inf. Model. 55, 2435–2442. doi: 10.1021/acs.jcim.5b00540
Cassuto, O., Dufies, M., Jacquel, A., Robert, G., Ginet, C., Dubois, A., et al. (2012). All tyrosine kinase inhibitor-resistant chronic myelogenous cells are highly sensitive to ponatinib. Oncotarget 3, 1557–1565. doi: 10.18632/oncotarget.692
Coutre, P. L., Tassi, E., Varella-Garcia, M., Barni, R., Mologni, L., Cabrita, G. A., et al. (2000). Induction of resistance to the Abelson inhibitor STI571 in human leukemic cells through gene amplification. Blood 95, 1758–1766. doi: 10.1182/blood.v95.5.1758.005a41_1758_1766
Crowley, L. C., Elzinga, B. M., O’sullivan, G. C., and Mckenna, S. L. (2011). Autophagy induction by Bcr-Abl-expressing cells facilitates their recovery from a targeted or nontargeted treatment. Am. J. Hematol. 86, 38–47. doi: 10.1002/ajh.21914
Crowley, L. C., and Waterhouse, N. J. (2016). Detecting Cleaved Caspase-3 in Apoptotic Cells by Flow Cytometry. Cold Spring Harbor, NY: Cold Spring Harbor Protocols.
de Groot, R. P., Raaijmakers, J. A., Lammers, J. W., Jove, R., and Koenderman, L. (1999). STAT5 activation by BCR-Abl contributes to transformation of K562 leukemia cells. Blood 94, 1108–1112. doi: 10.1182/blood.v94.3.1108.415k07_1108_1112
Deininger, M., Buchdunger, E., and Druker, B. J. (2005). The development of imatinib as a therapeutic agent for chronic myeloid leukemia. Blood 105, 2640–2653. doi: 10.1182/blood-2004-08-3097
Deininger, M. W. N., Goldman, J. M., and Melo, J. V. (2000). the molecular biology of chronic myeloid leukemia. Blood 96, 3343–3354.
Deng, M., Guan, X., Wen, X., Xiao, J., An, X., and Yu, J. (2020). Clinical efficacy and safety of imatinib treatment in children and adolescents with chronic myeloid leukemia: a single-center experience in China. Medicine (Baltimore) 99:e19150. doi: 10.1097/md.0000000000019150
Dong, B., Liang, Z., Chen, Z., Li, B., Zheng, L., Yang, J., et al. (2018). Cryptotanshinone suppresses key onco-proliferative and drug-resistant pathways of chronic myeloid leukemia by targeting STAT5 and STAT3 phosphorylation. Sci. China Life Sci. 61, 999–1009. doi: 10.1007/s11427-018-9324-y
Dube, D., Periwal, V., Kumar, M., Sharma, S., Singh, T. P., and Kaur, P. (2012). 3D-QSAR based pharmacophore modeling and virtual screening for identification of novel pteridine reductase inhibitors. J. Mol. Model. 18, 1701–1711. doi: 10.1007/s00894-011-1187-0
Elzinga, B. M., Nyhan, M. J., Crowley, L. C., O’donovan, T. R., Cahill, M. R., and Mckenna, S. L. (2013). Induction of autophagy by Imatinib sequesters Bcr-Abl in autophagosomes and down-regulates Bcr-Abl protein. Am. J. Hematol. 88, 455–462. doi: 10.1002/ajh.23428
Fiskus, W., Pranpat, M., Balasis, M., Bali, P., Estrella, V., Kumaraswamy, S., et al. (2006). Cotreatment with vorinostat (suberoylanilide hydroxamic acid) enhances activity of dasatinib (BMS-354825) against imatinib mesylate-sensitive or imatinib mesylate-resistant chronic myelogenous leukemia cells. Clin. Cancer Res. 12, 5869–5878. doi: 10.1158/1078-0432.ccr-06-0980
Gontarewicz, A., Balabanov, S., Keller, G., Colombo, R., Graziano, A., Pesenti, E., et al. (2008). Simultaneous targeting of Aurora kinases and Bcr-Abl kinase by the small molecule inhibitor PHA-739358 is effective against imatinib-resistant BCR-ABL mutations including T315I. Blood 111, 4355–4364. doi: 10.1182/blood-2007-09-113175
Gora-Tybor, J., and Robak, T. (2008). Targeted drugs in chronic myeloid leukemia. Curr. Med. Chem. 15, 3036–3051. doi: 10.2174/092986708786848578
Gumireddy, K., Baker, S. J., Cosenza, S. C., John, P., Kang, A. D., Robell, K. A., et al. (2005). A non-ATP-competitive inhibitor of BCR-ABL overrides imatinib resistance. Proc. Natl. Acad. Sci. U.S.A. 102, 1992–1997. doi: 10.1073/pnas.0408283102
Gupta, P., Kathawala, R. J., Wei, L., Wang, F., Wang, X., Druker, B. J., et al. (2016). PBA2, a novel inhibitor of imatinib-resistant BCR-ABL T315I mutation in chronic myeloid leukemia. Cancer Lett. 383, 220–229. doi: 10.1016/j.canlet.2016.09.025
Gupta, P., Zhang, G. N., Barbuti, A. M., Zhang, X., Karadkhelkar, N., Zhou, J., et al. (2020). Preclinical development of a novel BCR-ABL T315I inhibitor against chronic myeloid leukemia. Cancer Lett. 472, 132–141. doi: 10.1016/j.canlet.2019.11.040
Halbach, S., Hu, Z., Gretzmeier, C., Ellermann, J., Wöhrle, F. U., Dengjel, J., et al. (2016). Axitinib and sorafenib are potent in tyrosine kinase inhibitor resistant chronic myeloid leukemia cells. Cell Commun. Signal. 14:6.
Halgren, T. A., Murphy, R. B., Friesner, R. A., Beard, H. S., Frye, L. L., Pollard, W. T., et al. (2004). Glide: a new approach for rapid, accurate docking and scoring. 2. Enrichment factors in database screening. J. Med. Chem. 47, 1750–1759. doi: 10.1021/jm030644s
Hamilton, A., Elrick, L., Myssina, S., Copland, M., Jorgensen, H., Melo, J. V., et al. (2006). BCR-ABL activity and its response to drugs can be determined in CD34+ CML stem cells by CrkL phosphorylation status using flow cytometry. Leukemia 20, 1035–1039. doi: 10.1038/sj.leu.2404189
Hamurcu, Z., Delibasi, N., Gecene, S., Sener, E. F., Donmez-Altuntas, H., Ozkul, Y., et al. (2018). Targeting LC3 and Beclin-1 autophagy genes suppresses proliferation, survival, migration and invasion by inhibition of Cyclin-D1 and uPAR/Integrin beta1/Src signaling in triple negative breast cancer cells. J. Cancer Res. Clin. Oncol. 144, 415–430. doi: 10.1007/s00432-017-2557-5
Helgason, G. V., Karvela, M., and Holyoake, T. L. (2011). Kill one bird with two stones: potential efficacy of BCR-ABL and autophagy inhibition in CML. Blood 118, 2035–2043. doi: 10.1182/blood-2011-01-330621
Huang, D., Zhu, X., Tang, C., Mei, Y., Chen, W., Yang, B., et al. (2012). 3D QSAR pharmacophore modeling for c-Met kinase inhibitors. Med. Chem. 8, 1117–1125. doi: 10.2174/1573406411208061117
Huang, M., Dorsey, J. F., Epling-Burnette, P. K., Nimmanapalli, R., Landowski, T. H., Mora, L. B., et al. (2002). Inhibition of Bcr-Abl kinase activity by PD180970 blocks constitutive activation of Stat5 and growth of CML cells. Oncogene 21, 8804–8816. doi: 10.1038/sj.onc.1206028
Irwin, J. J., and Shoichet, B. K. (2005). ZINC–a free database of commercially available compounds for virtual screening. J. Chem. Inf. Model. 45, 177–182. doi: 10.1021/ci049714%2B
Jiang, Y., Gao, H., Liu, M., and Mao, Q. (2016). Sorting and biological characteristics analysis for side population cells in human primary hepatocellular carcinoma. Am. J. Cancer Res. 6, 1890–1905.
Kandakatla, N., and Ramakrishnan, G. (2014). Ligand based pharmacophore modeling and virtual screening studies to design novel HDAC2 inhibitors. Adv. Bioinform. 2014, 812148.
Katsoulas, A., Rachid, Z., Mcnamee, J. P., Williams, C., and Jean-Claude, B. J. (2008). Combi-targeting concept: an optimized single-molecule dual-targeting model for the treatment of chronic myelogenous leukemia. Mol. Cancer Ther. 7, 1033–1043. doi: 10.1158/1535-7163.mct-07-0179
Kavitha, R., Karunagaran, S., Chandrabose, S. S., Lee, K. W., and Meganathan, C. (2015). Pharmacophore modeling, virtual screening, molecular docking studies and density functional theory approaches to identify novel ketohexokinase (KHK) inhibitors. Biosystems 138, 39–52. doi: 10.1016/j.biosystems.2015.10.005
Khorashad, J. S., Kelley, T. W., Szankasi, P., Mason, C. C., Soverini, S., Adrian, L. T., et al. (2013). BCR-ABL1 compound mutations in tyrosine kinase inhibitor-resistant CML: frequency and clonal relationships. Blood 121, 489–498. doi: 10.1182/blood-2012-05-431379
Kim, S. J., Jung, K. H., Yan, H. H., Son, M. K., Fang, Z., Ryu, Y. L., et al. (2015). HS-543 induces apoptosis of Imatinib-resistant chronic myelogenous leukemia with T315I mutation. Oncotarget 6, 1507–1518. doi: 10.18632/oncotarget.2837
Kimura, S. (2006). Second generation Abl kinase inhibitors and novel compounds to eliminate the Bcr-Abl/T315I clone. Recent Pat. Anticancer Drug Discov. 1, 347–355. doi: 10.2174/157489206778776907
Kumar, R., Son, M., Bavi, R., Lee, Y., Park, C., Arulalapperumal, V., et al. (2015). Novel chemical scaffolds of the tumor marker AKR1B10 inhibitors discovered by 3D QSAR pharmacophore modeling. Acta Pharmacol. Sin. 36, 998–1012. doi: 10.1038/aps.2015.17
Lu, N., Li, X., Tan, R., An, J., Cai, Z., Hu, X., et al. (2018). HIF-1alpha/Beclin1-mediated autophagy is involved in neuroprotection induced by hypoxic preconditioning. J. Mol. Neurosci. 66, 238–250. doi: 10.1007/s12031-018-1162-7
Machova Polakova, K., Kulvait, V., Benesova, A., Linhartova, J., Klamova, H., Jaruskova, M., et al. (2015). Next-generation deep sequencing improves detection of BCR-ABL1 kinase domain mutations emerging under tyrosine kinase inhibitor treatment of chronic myeloid leukemia patients in chronic phase. J. Cancer Res. Clin. Oncol. 141, 887–899. doi: 10.1007/s00432-014-1845-6
Mizushima, N., and Yoshimori, T. (2007). How to interpret LC3 immunoblotting. Autophagy 3, 542–545. doi: 10.4161/auto.4600
Musoev, A., Numonov, S., You, Z., and Gao, H. (2019). Discovery of novel DPP-IV inhibitors as potential candidates for the treatment of type 2 Diabetes mellitus predicted by 3D QSAR pharmacophore models, molecular docking and de novo evolution. Molecules 24:2870. doi: 10.3390/molecules24162870
Nagar, B., Bornmann, W. G., Pellicena, P., Schindler, T., Veach, D. R., Miller, W. T., et al. (2002). Crystal structures of the kinase domain of c-Abl in complex with the small molecule inhibitors PD173955 and imatinib (STI-571). Cancer Res. 62, 4236–4243.
Nichols, G. L., Raines, M. A., Vera, J. C., Lacomis, L., Tempst, P., and Golde, D. W. (1994). Identification of CRKL as the constitutively phosphorylated 39-kD tyrosine phosphoprotein in chronic myelogenous leukemia cells. Blood 84, 2912–2918. doi: 10.1182/blood.v84.9.2912.2912
Oda, A., Miyakawa, Y., Druker, B. J., Ishida, A., Ozaki, K., Ohashi, H., et al. (1996). Crkl is constitutively tyrosine phosphorylated in platelets from chronic myelogenous leukemia patients and inducibly phosphorylated in normal platelets stimulated by thrombopoietin. Blood 88, 4304–4313. doi: 10.1182/blood.v88.11.4304.bloodjournal88114304
O’Hare, T., Walters, D. K., Stoffregen, E. P., Jia, T., Manley, P. W., Mestan, J., et al. (2005). In vitro activity of Bcr-Abl inhibitors AMN107 and BMS-354825 against clinically relevant imatinib-resistant Abl kinase domain mutants. Cancer Res. 65, 4500–4505. doi: 10.1158/0008-5472.can-05-0259
Okabe, S., Tauchi, T., Kimura, S., Maekawa, T., Kitahara, T., Tanaka, Y., et al. (2014). Combining the ABL1 kinase inhibitor ponatinib and the histone deacetylase inhibitor vorinostat: a potential treatment for BCR-ABL-positive leukemia. PLoS One 9:e89080. doi: 10.1371/journal.pone.0089080
Pal, S., Kumar, V., Kundu, B., Bhattacharya, D., Preethy, N., Reddy, M. P., et al. (2019). Ligand-based pharmacophore modeling, virtual screening and molecular docking studies for discovery of potential topoisomerase I inhibitors. Comput. Struct. Biotechnol. J. 17, 291–310. doi: 10.1016/j.csbj.2019.02.006
Pan, X., Dong, J., Shao, R., Su, P., Shi, Y., Wang, J., et al. (2015). Expanding the structural diversity of Bcr-Abl inhibitors: hybrid molecules based on GNF-2 and Imatinib. Bioorg. Med. Chem. Lett. 25, 4164–4168. doi: 10.1016/j.bmcl.2015.08.013
Peng, L., Jiang, H., and Bradley, C. (2001). [Annexin V for flow cytometric detection of phosphatidylserine expression on lymphoma cells undergoing apoptosis]. Hua Xi Yi Ke Da Xue Xue Bao 32, 602–604620.
Pereira, E. G., Moreira, M. A., and Caffarena, E. R. (2012). Molecular interactions of c-ABL mutants in complex with imatinib/nilotinib: a computational study using linear interaction energy (LIE) calculations. J. Mol. Model. 18, 4333–4341. doi: 10.1007/s00894-012-1436-x
Protein Data Bank. (1971). Available online at: http://www.rcsb.org/ (accessed May 4, 2018).
Quintás-Cardama, A., Tong, W., Manshouri, T., Vega, F., Lennon, P. A., Cools, J., et al. (2008). Activity of tyrosine kinase inhibitors against human NUP214-ABL1-positive T cell malignancies. Leukemia 22, 1117–1124. doi: 10.1038/leu.2008.80
Rampogu, S., Baek, A., Son, M., Park, C., Yoon, S., Parate, S., et al. (2020). Discovery of Lonafarnib-like compounds: pharmacophore modeling and molecular dynamics studies. ACS Omega 5, 1773–1781. doi: 10.1021/acsomega.9b02263
Rampogu, S., Son, M., Park, C., Kim, H. H., Suh, J. K., and Lee, K. W. (2017). Sulfonanilide derivatives in identifying novel aromatase inhibitors by applying docking, virtual screening, and MD simulations studies. Biomed. Res. Int. 2017:2105610.
Rampogu, S., Zeb, A., Baek, A., Park, C., Son, M., and Lee, K. W. (2018). Discovery of potential plant-derived peptide deformylase (PDF) inhibitors for multidrug-resistant bacteria using computational studies. J. Clin. Med. 7:563. doi: 10.3390/jcm7120563
Ren, R. (2005). Mechanisms of BCR-ABL in the pathogenesis of chronic myelogenous leukaemia. Nat. Rev. Cancer 5, 172–183. doi: 10.1038/nrc1567
Ren, X., Pan, X., Zhang, Z., Wang, D., Lu, X., Li, Y., et al. (2013). Identification of GZD824 as an orally bioavailable inhibitor that targets phosphorylated and nonphosphorylated breakpoint cluster region-Abelson (Bcr-Abl) kinase and overcomes clinically acquired mutation-induced resistance against imatinib. J. Med. Chem. 56, 879–894. doi: 10.1021/jm301581y
Rhodes, J., York, R. D., Tara, D., Tajinda, K., and Druker, B. J. (2000). CrkL functions as a nuclear adaptor and transcriptional activator in Bcr-Abl-expressing cells. Exp. Hematol. 28, 305–310. doi: 10.1016/s0301-472x(99)00148-4
Rossari, F., Minutolo, F., and Orciuolo, E. (2018). Past, present, and future of Bcr-Abl inhibitors: from chemical development to clinical efficacy. J. Hematol. Oncol. 11:84.
Sakkiah, S., Thangapandian, S., and Lee, K. W. (2012). Ligand-based virtual screening and molecular docking studies to identify the critical chemical features of potent cathepsin D inhibitors. Chem. Biol. Drug Des. 80, 64–79. doi: 10.1111/j.1747-0285.2012.01339.x
Salerno, L., Romeo, G., Modica, M. N., Amata, E., Sorrenti, V., Barbagallo, I., et al. (2017). Heme oxygenase-1: a new druggable target in the management of chronic and acute myeloid leukemia. Eur. J. Med. Chem. 142, 163–178. doi: 10.1016/j.ejmech.2017.07.031
Salomoni, P., and Calabretta, B. (2009). Targeted therapies and autophagy: new insights from chronic myeloid leukemia. Autophagy 5, 1050–1051. doi: 10.4161/auto.5.7.9509
Senechal, K., Halpern, J., and Sawyers, C. L. (1996). The CRKL adaptor protein transforms fibroblasts and functions in transformation by the BCR-ABL oncogene. J. Biol. Chem. 271, 23255–23261. doi: 10.1074/jbc.271.38.23255
Shinohara, H., Minami, Y., Naoe, T., and Akao, Y. (2019). Autophagic degradation determines the fate of T315I-mutated BCR-ABL protein. Haematologica 104, e191–e194.
Song, Y., Sun, X., Bai, W. L., and Ji, W. Y. (2013). Antitumor effects of Dasatinib on laryngeal squamous cell carcinoma in vivo and in vitro. Eur. Arch. Otorhinolaryngol. 270, 1397–1404. doi: 10.1007/s00405-013-2394-3
Tanaka, R., and Kimura, S. (2008). Abl tyrosine kinase inhibitors for overriding Bcr-Abl/T315I: from the second to third generation. Expert Rev. Anticancer Ther. 8, 1387–1398. doi: 10.1586/14737140.8.9.1387
Tu, Y. X., Wang, S. B., Fu, L. Q., Li, S. S., Guo, Q. P., Wu, Y., et al. (2018). Ovatodiolide targets chronic myeloid leukemia stem cells by epigenetically upregulating hsa-miR-155, suppressing the BCR-ABL fusion gene and dysregulating the PI3K/AKT/mTOR pathway. Oncotarget 9, 3267–3277. doi: 10.18632/oncotarget.23231
von Bubnoff, N., Schneller, F., Peschel, C., and Duyster, J. (2002). BCR-ABL gene mutations in relation to clinical resistance of Philadelphia-chromosome-positive leukaemia to STI571: a prospective study. Lancet 359, 487–491. doi: 10.1016/s0140-6736(02)07679-1
Wang, L. X., Wang, J. D., Chen, J. J., Long, B., Liu, L. L., Tu, X. X., et al. (2016). Aurora A kinase inhibitor AKI603 induces cellular senescence in chronic myeloid leukemia cells harboring T315I mutation. Sci. Rep. 6:35533.
Warsch, W., Grundschober, E., Berger, A., Gille, L., Cerny-Reiterer, S., Tigan, A.-S., et al. (2012). STAT5 triggers BCR-ABL1 mutation by mediating ROS production in chronic myeloid leukaemia. Oncotarget 3, 1669–1687.
Winter, S. S., Greene, J. M., Mcconnell, T. S., and Willman, C. L. (1999). Pre-B acute lymphoblastic leukemia with b3a2 (p210) and e1a2 (p190) BCR-ABL fusion transcripts relapsing as chronic myelogenous leukemia with a less differentiated b3a2 (p210) clone. Leukemia 13, 2007–2011. doi: 10.1038/sj.leu.2401598
Wu, L. X., Wu, Y., Chen, R. J., Liu, Y., Huang, L. S., Lou, L. G., et al. (2014). Curcumin derivative C817 inhibits proliferation of imatinib-resistant chronic myeloid leukemia cells with wild-type or mutant Bcr-Abl in vitro. Acta Pharmacol. Sin. 35, 401–409. doi: 10.1038/aps.2013.180
Ye, Y., Liao, Q., Wei, J., and Gao, Q. (2010). 3D-QSAR study of corticotropin-releasing factor 1 antagonists and pharmacophore-based drug design. Neurochem. Int. 56, 107–117. doi: 10.1016/j.neuint.2009.09.008
Yoshii, S. R., and Mizushima, N. (2017). Monitoring and measuring autophagy. Int. J. Mol. Sci. 18:1865. doi: 10.3390/ijms18091865
Yu, Y., Yang, L., Zhao, M., Zhu, S., Kang, R., Vernon, P., et al. (2012). Targeting microRNA-30a-mediated autophagy enhances imatinib activity against human chronic myeloid leukemia cells. Leukemia 26, 1752–1760. doi: 10.1038/leu.2012.65
Zabriskie, M. S., Eide, C. A., Tantravahi, S. K., Vellore, N. A., Estrada, J., Nicolini, F. E., et al. (2014). BCR-ABL1 compound mutations combining key kinase domain positions confer clinical resistance to ponatinib in Ph chromosome-positive leukemia. Cancer Cell 26, 428–442. doi: 10.1016/j.ccr.2014.07.006
Zabriskie, M. S., Vellore, N. A., Gantz, K. C., Deininger, M. W., and O’hare, T. (2015). Radotinib is an effective inhibitor of native and kinase domain-mutant BCR-ABL1. Leukemia 29, 1939–1942. doi: 10.1038/leu.2015.42
Zheng, X., Ekins, S., Raufman, J. P., and Polli, J. E. (2009). Computational models for drug inhibition of the human apical sodium-dependent bile acid transporter. Mol. Pharm. 6, 1591–1603. doi: 10.1021/mp900163d
Keywords: CML, BCR-ABL, pharmacophore model, ZINC21710815, apoptosis, autophagy
Citation: Huang T-T, Wang X, Qiang S-J, Zhao Z-N, Wu Z-X, Ashby CR Jr, Li J-Z and Chen Z-S (2021) The Discovery of Novel BCR-ABL Tyrosine Kinase Inhibitors Using a Pharmacophore Modeling and Virtual Screening Approach. Front. Cell Dev. Biol. 9:649434. doi: 10.3389/fcell.2021.649434
Received: 04 January 2021; Accepted: 10 February 2021;
Published: 04 March 2021.
Edited by:
Ri Cui, Wenzhou Medical University, ChinaReviewed by:
Andaleeb Sajid, Yale University, United StatesPeng Zou, Wenzhou Medical University, China
Jin-Jian Lu, University of Macau, China
Copyright © 2021 Huang, Wang, Qiang, Zhao, Wu, Ashby, Li and Chen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jia-Zhong Li, lijiazhong@lzu.edu.cn; Zhe-Sheng Chen, chenz@stjohns.edu