In silico identification of novel kinase inhibitors targeting wild-type and T315I mutant ABL1 from FDA-approved drugs

Anticancer potential of mebendazole against chronic myeloid leukemia: in silico and in vitro studies revealed new insights about the mechanism of action

Chronic myeloid leukemia (CML) is caused by constitutively active fusion protein BCR-ABL1, and targeting ABL1 is a promising therapy option. Imatinib, dasatinib, and nilotinib have all been shown to work effectively in clinical trials. ABL1 mutations, particularly the T315I gate-keeper mutation, cause resistance in patients. As a result, broad-spectrum ABL1 medicines are desperately needed. In order to screen potential drugs targeting CML, mebendazole (MBZ) was subjected to the in vitro test against CML cell lines (K562 and FEPS) and computational assays. The antiproliferative effect of MBZ and the combination with tyrosine kinase inhibitors (TKIs) was tested using end-point viability assays, cell cycle distribution analysis, cell membrane, and mitochondrial dyes. By interrupting the cell cycle and causing cell death, MBZ and its combination with imatinib and dasatinib have a significant antiproliferative effect. We identified MBZ as a promising “new use” drug targeting wild-type and mutant ABL1 using molecular docking. Meanwhile, we determined which residues in the allosteric site are important in ABL1 drug development. These findings may not only serve as a model for repositioning current authorized medications but may also provide ABL1-targeted anti-CML treatments a fresh lease of life.

Computational and Biological Investigations on Abl1 Tyrosine Kinase: A Review

Abl1 tyrosine kinase is a validated target for the treatment of chronic myeloid leukemia. It is a form of cancer that is difficult to treat and much research is being done to identify new molecular entities and to tackle drug resistance issues. In recent years, drug resistance of Abl1 tyrosine kinase has become a major healthcare concern. Second and third-generation TKI reported better responses against the resistant forms; still they had no impact on long-term survival prolongation. New compounds derived from natural products and organic small molecule inhibitors can lay the foundation for better clinical therapies in the future. Computational methods, experimental and biological studies can help us understand the mechanism of drug resistance and identify novel molecule inhibitors. ADMET parameters analysis of reported drugs and novel small molecule inhibitors can also provide valuable insights. In this review, available therapies, point mutations, structure-activity relationship and ADMET parameters of reported series of Abl1 tyrosine kinase inhibitors and drugs are summarised. We summarise in detail recent computational and molecular biology studies that focus on designing drug molecules, investigation of natural product compounds and organic new chemical entities. Current ongoing research suggests that selective targeting of Abl1 tyrosine kinase at the molecular level to combat drug resistance in chronic myeloid leukemia is promising.

Secondary Resistant Mutations to Small Molecule Inhibitors in Cancer Cells

Secondary resistant mutations in cancer cells arise in response to certain small molecule inhibitors. These mutations inevitably cause recurrence and often progression to a more aggressive form. Resistant mutations may manifest in various forms. For example, some mutations decrease or abrogate the affinity of the drug for the protein. Others restore the function of the enzyme even in the presence of the inhibitor. In some cases, resistance is acquired through activation of a parallel pathway which bypasses the function of the drug targeted pathway. The Catalogue of Somatic Mutations in Cancer (COSMIC) produced a compendium of resistant mutations to small molecule inhibitors reported in the literature. Here, we build on these data and provide a comprehensive review of resistant mutations in cancers. We also discuss mechanistic parallels of resistance.

New drugs are not enough‑drug repositioning in oncology: An update

Drug repositioning refers to the concept of discovering novel clinical benefits of drugs that are already known for use treating other diseases. The advantages of this are that several important drug characteristics are already established (including efficacy, pharmacokinetics, pharmacodynamics and toxicity), making the process of research for a putative drug quicker and less costly. Drug repositioning in oncology has received extensive focus. The present review summarizes the most prominent examples of drug repositioning for the treatment of cancer, taking into consideration their primary use, proposed anticancer mechanisms and current development status.

Computer-aided screening for suppressor of variegation 4-20 homolog 1 inhibitors and their preliminary activity validation in human osteosarcoma

Histone methylation/demethylation facilitate to maintain balanced histone methylation levels and underpin gene regulation, playing the key roles in epigenetic regulation. Suppressor of variegation 4-20 homolog 1 (SUV420H1), a member of class Histone Lysine Methyltransferase and a key enzyme in the epigenetic regulation of the pathways controlling metabolism and tumorigenesis, is crucial to maintain cell homeostasis. The inhibition of SUV420H1 has emerged as a promising candidate for drug development and cancer therapy. Herein, two potential and potent SUV420H1 inhibitors (ZINC08398384, ZINC08439608) were identified through in silico approach and in vitro biological experiments. In vitro biological tests demonstrated that these compounds can inhibit the proliferation of U2OS cells and restrict its migration ability. And the level of dimethylation of lysine 20 on histone H4 (H4K20me2) was markedly decreased by these compounds-treatment in a dose-dependent manner. These results indicated that ZINC08398384 and ZINC08439608 are potential SUV420H1 inhibitors and could be developed as promising drug candidates applied to cancer epigenetic therapy.Communicated by Ramaswamy H. Sarma.

Computer-aided identification of a novel pyruvate kinase M2 activator compound

OBJECTIVES

The aim of this study was to obtain antitumour molecules targeting to activate PKM2 through adequate computational methods combined with biological activity experiments.

METHODS

The structure-based virtual screening was utilized to screen effective activator targeting PKM2 from ZINC database. Molecular dynamics simulations were performed to evaluate the stability of the small molecule-binding PKM2 complex systems. Then, cell survival experiments, glutaraldehyde crosslinking reaction, western blot, and qPCR experiments were used to detect the effects of top hits on various cancer cells and the targeting specificity of PKM2.

RESULTS

Two small molecules in 1,5-2H-pyrrole-dione were obtained after virtual screening. In vitro experiments demonstrated that ZINC08383544 specifically activated PKM2 and affected the expression of upstream and downstream genes of PKM2 during glycolysis, leading to the inhibition of tumour cell growth. These results indicate that ZINC08383544 conforms to the characteristics of PKM2 activator and is potential to be a novel PKM2 activator as antitumour drug.

DISCUSSION

This work proves that ZINC08383544 promotes the formation of PKM2 tetramer, effectively blocks PKM2 nuclear translocation, and inhibits the growth of tumour, and ZINC08383544 may be a novel activator of PKM2. This work may provide a good choice of drug or molecular fragments for the antitumour strategy targeting PKM2. Screening of targeted drugs by combination of virtual screening and bioactivity experiments is a rapid method for drug discovery.

Identification of two potential glycogen synthase kinase 3β inhibitors for the treatment of osteosarcoma

Osteosarcoma is the most common primary malignant bone tumor among adolescents worldwide with high mortality rate. Glycogen synthase kinase 3β (GSK3β) is a serine/threonine kinase and is considered as a validated target in osteosarcoma therapy. Therefore, the study of GSK3β inhibitors is one of the most popular fields in anti-osteosarcoma drug development. Here, the tools of bioinformatics were used to screen novel effective inhibitors of GSK3β from ZINC Drug Database. The molecular docking, molecular dynamic simulations, MM/GBSA, and energy decomposition analysis were performed to identify the inhibitors. Finally, ZINC08383479 and ZINC08441251 were selected as potential GSK3β inhibitors. These two inhibitors were evaluated by GSK3β kinase inhibition assay in vitro. The inhibition of cell proliferation was tested in osteosarcoma cell lines U2OS and MG63 in vitro. The result showed that ZINC08383479 and ZINC08441251 had high inhibition activity against GSK3β. We found that CHIR99021 (a known GSK3β inhibitor), ZINC08383479, and ZINC08441251 had significant inhibition activity in U2OS cells and MG63 cells. These findings may provide new ideas for the design of more potent GSK3β inhibitors and therapeutic targets for osteosarcoma.

In silico identification of small molecules as novel LXR agonists for the treatment of cardiovascular disease and cancer

Liver X receptor (LXR), a member of the nuclear receptor superfamily, mainly serves as a reverse cholesterol transporter in lipid metabolism. It has been demonstrated that LXR is a promising target for the treatment of cardiovascular diseases. LXR is also involved in cancer metabolism, glucose homeostasis, immunity, and various physiological processes. The antitumor function of LXR has become of great interest to researchers in recent years. However, while it is believed that activating LXR with small molecules could be a promising approach to cancer treatment, effective drugs that target LXR are yet to be reported. To find compounds that are potentially capable of activating LXR, we utilized a high-throughput screening method to search the MolMall database for suitable compounds. Seven candidates with lower GB/SA Hawkins scores than the reference ligand T0901317 were identified. Based on the results of molecular dynamics (MD) simulations, binding free energy analysis, and an analysis of the agonism mechanism, ZINC90512020 and ZINC3845032 were predicted to have high affinities for LXR and high relative stabilization, and were therefore selected as potential LXR agonists. Both of these compounds will undergo further development with a view to utilizing them for the treatment of LXR-related cardiovascular diseases or cancers.

Structure and function analysis of Polygonatum cyrtonema lectin by site-directed mutagenesis

The crystal structure of mature Polygonatum cyrtonema lectin (PCL) showed three similar carbohydrate-binding sites (CBS I, CBS II, and CBS III). The Gln58 and Asp60 residues of CBS II are substituted with His58 and Asn60. To establish the relationship between the key amino acid residues and structure or activity of PCL, we constructed four recombinant mutants in CBS I, CBS II, and CBS III. The experimental results indicate that CBS I, CBS III and the disulfide bond play vital roles in the binding with mannose. Furthermore, molecular dynamics simulations and binding free energy calculation illustrate that CBS I has a direct and strong relationship with the activity of PCL. CBS II does not play a critical role in the model for mannose binding by PCL. Although CBS III does not enhance the activity, it helps to maintain the activity and 3D structure. These results suggest that the carbohydrate-binding site of PCL may be in a hydrophilic environment, and Asn and Tyr are the key amino acids involved in its binding with sugar, but Gln and Asp are not necessary to maintain its activity.

In silico drug repurposing of FDA-approved drugs to predict new inhibitors for drug resistant T315I mutant and wild-type BCR-ABL1: A virtual screening and molecular dynamics study

The BCR-ABL fusion gene is one of the major causes of 95% of Chronic Myeloid Leukemia (CML). While, BCR-ABL protein is currently being used as a major target to treat CML. Although, current FDA-approved drugs such as; Imatinib and Nilotinib have stupendously improved the patients 5-year's survival rates, the drug resistance has dramatically reduced their effects. So, more accurate and effective alternative treatments are crucially needed. To address this issue, we screened the FDA-approved drugs by virtual screening and binding free energy calculations to identify new inhibitors for the wild-type and T315I gatekeeper mutant ABL1. It was invigorating to identify that chlorohexidine, paromomycin and deferoxamine could inhibit the wild-type ABL1, while chlorohexidine and ritonavir could inhibit the T315I mutant ABL1. The applications of these newly identified drugs are not just an effortless hypothesis in drug discovery. These drugs can be evaluated in phase 2 clinical trials after a simple kinase selectivity assay.

Insight into drug resistance mechanisms and discovery of potential inhibitors against wild-type and L1196M mutant ALK from FDA-approved drugs

Anaplastic lymphoma kinase (ALK) plays a crucial role in multiple malignant cancers. It is known as a well-established target for the treatment of ALK-dependent cancers. Even though substantial efforts have been made to develop ALK inhibitors, only crizotinib, ceritinib, and alectinib had been approved by the U.S. Food and Drug Administration for patients with ALK-positive non-small cell lung cancer (NSCLC). The secondary mutations with drug-resistance bring up difficulties to develop effective drugs for ALK-positive cancers. To give a comprehensive understanding of molecular mechanism underlying inhibitor response to ALK tyrosine kinase mutations, we established an accurate assessment for the extensive profile of drug against ALK mutations by means of computational approaches. The molecular mechanics-generalized Born surface area (MM-GBSA) method based on molecular dynamics (MD) simulation was carried out to calculate relative binding free energies for receptor-drug systems. In addition, the structure-based virtual screening was utilized to screen effective inhibitors targeting wild-type ALK and the gatekeeper mutation L1196M from 3180 approved drugs. Finally, the mechanism of drug resistance was discussed, several novel potential wild-type and L1196M mutant ALK inhibitors were successfully identified.

Discovery of a Potential HER2 Inhibitor from Natural Products for the Treatment of HER2-Positive Breast Cancer

Breast cancer is one of the most lethal types of cancer in women worldwide due to the late stage detection and resistance to traditional chemotherapy. The human epidermal growth factor receptor 2 (HER2) is considered as a validated target in breast cancer therapy. Even though a substantial effort has been made to develop HER2 inhibitors, only lapatinib has been approved by the U.S. Food and Drug Administration (FDA). Side effects were observed in a majority of the patients within one year of treatment initiation. Here, we took advantage of bioinformatics tools to identify novel effective HER2 inhibitors. The structure-based virtual screening combined with ADMET (absorption, distribution, metabolism, excretion and toxicity) prediction was explored. In total, 11,247 natural compounds were screened. The top hits were evaluated by an in vitro HER2 kinase inhibition assay. The cell proliferation inhibition effect of identified inhibitors was evaluated in HER2-overexpressing SKBR3 and BT474 cell lines. We found that ZINC15122021 showed favorable ADMET properties and attained high binding affinity against HER2. Moreover, ZINC15122021 showed high kinase inhibition activity against HER2 and presented outstanding cell proliferation inhibition activity against both SKBR3 and BT474 cell lines. Results reveal that ZINC15122021 can be a potential HER2 inhibitor.

Identification of Novel Pathways in Plant Lectin-Induced Cancer Cell Apoptosis

Plant lectins have been investigated to elucidate their complicated mechanisms due to their remarkable anticancer activities. Although plant lectins seems promising as a potential anticancer agent for further preclinical and clinical uses, further research is still urgently needed and should include more focus on molecular mechanisms. Herein, a Naïve Bayesian model was developed to predict the protein-protein interaction (PPI), and thus construct the global human PPI network. Moreover, multiple sources of biological data, such as smallest shared biological process (SSBP), domain-domain interaction (DDI), gene co-expression profiles and cross-species interolog mapping were integrated to build the core apoptotic PPI network. In addition, we further modified it into a plant lectin-induced apoptotic cell death context. Then, we identified 22 apoptotic hub proteins in mesothelioma cells according to their different microarray expressions. Subsequently, we used combinational methods to predict microRNAs (miRNAs) which could negatively regulate the abovementioned hub proteins. Together, we demonstrated the ability of our Naïve Bayesian model-based network for identifying novel plant lectin-treated cancer cell apoptotic pathways. These findings may provide new clues concerning plant lectins as potential apoptotic inducers for cancer drug discovery.

Screening of novel inhibitors targeting lactate dehydrogenase A via four molecular docking strategies and dynamics simulations

Lactate dehydrogenase A (LDHA) is a metabolic enzyme which catalyzes the interconversion of lactate and pyruvate in the glycolysis pathway, thus playing key roles in aerobic glycolysis. The inhibition of LDHA by small molecules has become an attractive strategy for anticancer therapy in recent years. However, very few LDHA inhibitors have been reported, even though a great deal of effort has directed into identifying LDHA inhibitors using structure-based approaches. Therefore, high-throughput and high-accuracy screening approaches are still urgently needed in order to target LDHA effectively. In the present work, after establishing that our docking strategies performed well using test datasets, we screened 32791 Specs products for their docking scores with the substrate-binding pocket and, separately, the cofactor-binding pocket of LDHA. We subsequently identified 76 hits (i.e., ligands that show low docking scores) for the cofactor-binding pocket and 27 hits for the substrate-binding pocket. Two representative compounds, ZINC20036549 and ZINC19369718, were then chosen for further MD simulation analysis, and we found that these compounds maintained their inhibitory activity during the MD simulations. Meanwhile, we found that ZINC19369718 interacts with a novel binding site close to the active site, and that this interaction may inhibit the catalytic activity of LDHA. Together, these results offer not only a new paradigm for identifying Specs drug-like products for novel therapeutic use but they also provide further opportunity to adopt LDHA inhibition as a strategy for cancer therapy.

Induction of heme oxygenase-1 by Na+-H+ exchanger 1 protein plays a crucial role in imatinib-resistant chronic myeloid leukemia cells

Resistance toward imatinib (IM) and other BCR/ABL tyrosine kinase inhibitors remains troublesome in the treatment of advanced stage chronic myeloid leukemia (CML). The aim of this study was to estimate the reversal effects of down-regulation of Na(+)/H(+) exchanger 1 (NHE1) on the chemoresistance of BCR-ABL-positive leukemia patients' cells and cell lines. After treatment with the specific NHE1 inhibitor cariporide to decrease intracellular pH (pHi), the heme oxygenase-1 (HO-1) levels of the K562R cell line and cells from IM-insensitive CML patients decreased. HO-1, as a Bcr/Abl-dependent survival molecule in CML cells, is important for the resistance to tyrosine kinase inhibitors in patients with newly diagnosed CML or IM-resistant CML. Silencing PKC-β and Nrf-2 or treatment with inhibitors of p38 pathways obviously blocked NHE1-induced HO-1 expression. Furthermore, treatment with HO-1 or p38 inhibitor plus IM increased the apoptosis of the K562R cell line and IM-insensitive CML patients' cells. Inhibiting HO-1 enhanced the activation of caspase-3 and poly(ADP-ribose) polymerase-1. Hence, the results support the anti-apoptotic role of HO-1 induced by NHE1 in the K562R cell line and IM-insensitive CML patients and provide a mechanism by which inducing HO-1 expression via the PKC-β/p38-MAPK pathway may promote tumor resistance to oxidative stress.

Discovery of potential drugs for human-infecting H7N9 virus containing R294K mutation

Background

After the first epidemic wave from February through May 2013, the influenza A (H7N9) virus emerged and has followed a second epidemic wave since June 2013. As of June 27, 2014, the outbreak of H7N9 had caused 450 confirmed cases of human infection, with 165 deaths included. The case-fatality rate of all confirmed cases is about 36%, making the H7N9 virus a significant threat to people’s health. At present, neuraminidase inhibitors are the only licensed antiviral medications available to treat H7N9 infections in humans. Oseltamivir is the most commonly used inhibitor, and it is also a front-line drug for the threatening H7N9. Unfortunately, it has been reported that patients treated with oseltamivir can induce R294K (Arg294Lys) substitution in the H7N9 virus, which is a rare mutation and can reduce the antiviral efficacy of inhibitors. Even worse, deaths caused by such mutation after oseltamivir treatment have already been reported, indicating that the need to find substitutive neuraminidase inhibitors for currently available drugs to treat drug-resistant H7N9 is really pressing.

Materials and methods

First, the structure of H7N9 containing the R294K substitution was downloaded from the Protein Data Bank, and structural information of approved drugs was downloaded from the ZINC (ZINC Is Not Commercial) database. Taking oseltamivir carboxylate as a reference drug, we then filtered these molecules through virtual screening to find out potential inhibitors targeting the mutated H7N9 virus. For further evaluation, we carried out a 14 ns molecular dynamic simulation for each H7N9–drug complex and calculated the binding energy for each candidate drug.

Results

We found five inhibitors that could be candidate drugs for treating the mutated H7N9 virus. Docking poses showed these drugs could bind to the virus effectively, with the contribution of hydrogen bonds and hydrophobic interactions. With regard to the molecular dynamic simulations, receptor–ligand complexes formed by these candidate drugs were more stable than the one formed by oseltamivir carboxylate. The binding energy of oseltamivir carboxylate was −122.4 kJ/mol, while those for these potential inhibitors were −417.5, −404.7, −372.2, −304.3, and −289.9 kJ/mol, much better than the reference drug.

Conclusion

Given the current and future threat of the mutated H7N9 virus, it is urgent that potent drugs and effective antiviral therapeutics be found. Our study therefore is able to complement currently available drugs for influenza A infectors and helps to prevent the ongoing threat of H7N9 virus.

In silico discovery of potential VEGFR-2 inhibitors from natural derivatives for anti-angiogenesis therapy

Angiogenesis is the growth of new capillaries from existing blood vessels that supply oxygen and nutrients and provide gateways for immune surveillance. Abnormal vessel growth in term of excessive angiogenesis is a hallmark of cancer, inflammatory and eye diseases. VEGFR-2 (vascular endothelial growth factor receptor 2) dominating the process of angiogenesis has led to approval of therapeutic inhibitors and is becoming a promising target for anti-angiogenic drugs. Notwithstanding these successes, the clinical use of current VEGFR-2 blockers is more challenging than anticipated. Taking axitinib as a reference drug, in our study we found three potent VEGFR-2 inhibitors (ZINC08254217, ZINC08254138, and ZINC03838680) from natural derivatives. Each of the three inhibitors acquired a better grid score than axitinib (−62.11) when docked to VEGFR-2. Molecular dynamics simulations demonstrated that ZINC08254217– and ZINC08254138–VEGFR-2 complexes were more stable than axitinib. Similar to bind free energy for axitinib (−54.68 kcal/mol), such for ZINC03838680, ZINC08254217, and ZINC08254138 was −49.37, −43.32, and −32.73 kcal/mol respectively. These results suggested these three compounds could be candidate drugs against angiogenesis, with comparable VEGFR-2 binding affinity of axitinib. Hence findings in our study are able to provide valuable information on discovery of effective anti-angiogenesis therapy.