Combined pharmacophore modeling, docking, and 3D-QSAR studies of PLK1 inhibitors

Polo-like kinase 1, an important enzyme with diverse biological actions in cell mitosis, is a promising target for developing novel anticancer drugs. A combined molecular docking, structure-based pharmacophore modeling and three-dimensional quantitative structure-activity relationship (3D-QSAR) study was performed on a set of 4,5-dihydro-1H-pyrazolo[4,3-h]quinazoline derivatives as PLK1 inhibitors. The common substructure, molecular docking and pharmacophore-based alignment were used to develop different 3D-QSAR models. The comparative molecular field analysis (CoMFA) and comparative molecule similarity indices analysis (CoMSIA) models gave statistically significant results. These models showed good q(2) and r(2) (pred) values and revealed a good response to test set validation. All of the structural insights obtained from the 3D-QSAR contour maps are consistent with the available crystal structure of PLK1. The contour maps obtained from the 3D-QSAR models in combination with the structure based pharmacophore model help to better interpret the structure-activity relationship. These satisfactory results may aid the design of novel PLK1 inhibitors. This is the first report on 3D-QSAR study of PLK1 inhibitors.     

Identification of Green Tea Catechins as Potent Inhibitors of the Polo‐Box Domain of Polo‐Like kinase 1

Polo‐like kinase 1 (PLK1) plays crucial functions in multiple stages of mitosis and is considered to be a potential drug target for cancer therapy. The functions of PLK1 are mediated by its N‐terminal kinase domain and C‐terminal polo‐box domain (PBD). Most inhibitors targeting the kinase domain of PLK1 have a selectivity issue because of a high degree of structural conservation within kinase domains of all protein kinases. Here, we combined virtual and experimental screenings to identify green tea catechins as potent inhibitors of the PLK1 PBD. Initially, (?)‐epigallocatechin, one of the main components of green tea polyphenols, was found to significantly block the binding of fluorescein‐labeled phosphopeptide to the PBD at a concentration of 10 μm. Next, additional catechins were evaluated for their dose‐dependent inhibition of the PBD and preliminary structure–activity relationships were derived. Cellular analysis further showed that catechins interfere with the proper subcellular localization of PLK1, lead to cell‐cycle arrest in the S and G2M phases, and induce growth inhibition of several human cancer cell types, such as breast adenocarcinoma (MCF7), lung adenocarcinoma (A549), and cervical adenocarcinoma (HeLa). Our data provides new insight into understanding the anticancer activities of green tea catechins.     

Peptidomimetic Polo-Box-Targeted Inhibitors that Engage PLK1 in Tumor Cells and Are Selective against the PLK3 Tumor Suppressor

The polo-box domain (PBD) of PLK1 determines mitotic substrate recognition and subcellular localization. Compounds that target PLK1 selectively are required due to the tumor-suppressor roles of PLK3. A structure-activity analysis of the PBD phosphopeptide binding motif has identified potent peptides that delineate the determinants required for mimicry by nonpeptidic inhibitors and provide insights into the structural basis for the selectivity of inhibitors for the PLK1 PBD. Fragment-ligated inhibitory peptides (FLIPs) obtained through REPLACE have been optimized to enhance in vitro binding and a systematic analysis of selectivity for PLK1 vs PLK3 has been carried out for peptides and peptidomimetics. Furthermore, these more drug-like non-ATP-competitive inhibitors had on-target engagement in a cellular context, as evidenced by stabilization of PLK1 in a thermal-shift assay and by inhibition of the phosphorylation of TCTP, a target of PLK1. Investigation in cells expressing a mutant PLK1 showed that these cells are sensitive to PBD inhibitors but dramatically resistant to clinically investigated ATP-competitive compounds. These results further validate targeting the PBD binding site in the move towards PLK1 inhibitors that are active against tumors resistant to ATP inhibitors.     

3D-QSAR Study of Potent Inhibitors of Phosphodiesterase-4 Using a CoMFA Approach

Phosphodiesterase-4 (PDE4) plays an important role in treatment of asthma and chronic obstructive pulmonary disease. Thirty-one analogs displaying variable inhibition of PDE4 were selected to develop models for establishing three-dimensional quantitative structure-activity relationships (3D-QSAR). Comparative molecular field analysis (CoMFA) was conducted on the group of analogs to determine the structural requirements for potency in inhibiting PDE4. The resulting model exhibited good q² and r² values up to 0.741 and 0.954 for CoMFA. The contributions from the steric and electrostatic fields were 0.915 and 0.085 respectively. The 3D-QSAR model should be very useful for design of novel PDE 4 inhibitors.     

Combinatorial Pharmacophore-Based 3D-QSAR Analysis and Virtual Screening of FGFR1 Inhibitors

The fibroblast growth factor/fibroblast growth factor receptor (FGF/FGFR) signaling pathway plays crucial roles in cell proliferation, angiogenesis, migration, and survival. Aberration in FGFRs correlates with several malignancies and disorders. FGFRs have proved to be attractive targets for therapeutic intervention in cancer, and it is of high interest to find FGFR inhibitors with novel scaffolds. In this study, a combinatorial three-dimensional quantitative structure-activity relationship (3D-QSAR) model was developed based on previously reported FGFR1 inhibitors with diverse structural skeletons. This model was evaluated for its prediction performance on a diverse test set containing 232 FGFR inhibitors, and it yielded a SD value of 0.75 pIC₅₀ units from measured inhibition affinities and a Pearson’s correlation coefficient R² of 0.53. This result suggests that the combinatorial 3D-QSAR model could be used to search for new FGFR1 hit structures and predict their potential activity. To further evaluate the performance of the model, a decoy set validation was used to measure the efficiency of the model by calculating EF (enrichment factor). Based on the combinatorial pharmacophore model, a virtual screening against SPECS database was performed. Nineteen novel active compounds were successfully identified, which provide new chemical starting points for further structural optimization of FGFR1 inhibitors.     

A Multi-step Virtual Screening Protocol for the Identification of Novel Non-acidic Microsomal Prostaglandin E2 Synthase-1 (mPGES-1) Inhibitors

Microsomal prostaglandin E₂ synthase-1 (mPGES-1) is a potential therapeutic target for the treatment of inflammatory diseases and certain types of cancer. To identify novel scaffolds for mPGES-1 inhibition, we applied a virtual screening (VS) protocol that comprises molecular docking, fingerprints-based clustering with diversity-based selection, protein–ligand interactions fingerprints, and molecular dynamics (MD) simulations with molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations. The hits identified were carefully analyzed to ensure the selection of novel scaffolds that establish stable interactions with key residues in the mPGES-1 binding pocket and inhibit the catalytic activity of the enzyme. As a result, we discovered two promising chemotypes, 4-(2-chlorophenyl)-N-[(2-{[(propan-2-yl)sulfamoyl]methyl}phenyl)methyl]piperazine-1-carboxamide ( 6 ) and N-(4-methoxy-3-{[4-(6-methyl-1,3-benzothiazol-2-yl)phenyl]sulfamoyl}phenyl)acetamide ( 8 ), as non-acidic mPGES-1 inhibitors with IC₅₀ values of 1.2 and 1.3 μm , respectively. Minimal structural optimization of 8 resulted in three more compounds with promising improvements in inhibitory activity (IC₅₀: 0.3–0.6 μm ). The unprecedented chemical structures of 6 and 8 , which are amenable to further derivatization, reveal a new and attractive approach for the development of mPGES-1 inhibitors with potential anti-inflammatory and anticancer properties.     

Combined 3D-QSAR Modeling and Molecular Docking Studies on Pyrrole-Indolin-2-ones as Aurora A Kinase Inhibitors

Aurora kinases have emerged as attractive targets for the design of anticancer drugs. 3D-QSAR (comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA)) and Surflex-docking studies were performed on a series of pyrrole-indoline-2-ones as Aurora A inhibitors. The CoMFA and CoMSIA models using 25 inhibitors in the training set gave r(2) (cv) values of 0.726 and 0.566, and r(2) values of 0.972 and 0.984, respectively. The adapted alignment method with the suitable parameters resulted in reliable models. The contour maps produced by the CoMFA and CoMSIA models were employed to rationalize the key structural requirements responsible for the activity. Surflex-docking studies revealed that the sulfo group, secondary amine group on indolin-2-one, and carbonyl of 6,7-dihydro-1H-indol-4(5H)-one groups were significant for binding to the receptor, and some essential features were also identified. Based on the 3D-QSAR and docking results, a set of new molecules with high predicted activities were designed.     

Molecular Targeting of the Oncoprotein PLK1 in Pediatric Acute Myeloid Leukemia: RO3280, a Novel PLK1 Inhibitor,Induces Apoptosis in Leukemia Cells

Polo-like kinase 1 (PLK1) is highly expressed in many cancers and therefore a biomarker of transformation and potential target for the development of cancer-specific small molecule drugs. RO3280 was recently identified as a novel PLK1 inhibitor; however its therapeutic effects in leukemia treatment are still unknown. We found that the PLK1 protein was highly expressed in leukemia cell lines as well as 73.3% (11/15) of pediatric acute myeloid leukemia (AML) samples. PLK1 mRNA expression was significantly higher in AML samples compared with control samples (82.95 ± 110.28 vs. 6.36 ± 6.35; p < 0.001). Kaplan-Meier survival analysis revealed that shorter survival time correlated with high tumor PLK1 expression (p = 0.002). The 50% inhibitory concentration (IC₅₀) of RO3280 for acute leukemia cells was between 74 and 797 nM. The IC₅₀ of RO3280 in primary acute lymphocytic leukemia (ALL) and AML cells was between 35.49 and 110.76 nM and 52.80 and 147.50 nM, respectively. RO3280 induced apoptosis and cell cycle disorder in leukemia cells. RO3280 treatment regulated several apoptosis-associated genes. The regulation of DCC, CDKN1A, BTK, and SOCS2 was verified by western blot. These results provide insights into the potential use of RO3280 for AML therapy; however, the underlying mechanisms remain to be determined.     

Plumbagin,a Natural Product with Potent Anticancer Activities,Binds to and Inhibits Dihydroorotase,a Key Enzyme in Pyrimidine Biosynthesis

Dihydroorotase (DHOase) is the third enzyme in the de novo biosynthesis pathway for pyrimidine nucleotides, and an attractive target for potential anticancer chemotherapy. By screening plant extracts and performing GC–MS analysis, we identified and characterized that the potent anticancer drug plumbagin (PLU), isolated from the carnivorous plant Nepenthes miranda, was a competitive inhibitor of DHOase. We also solved the complexed crystal structure of yeast DHOase with PLU (PDB entry 7CA1), to determine the binding interactions and investigate the binding modes. Mutational and structural analyses indicated the binding of PLU to DHOase through loop-in mode, and this dynamic loop may serve as a drug target. PLU exhibited cytotoxicity on the survival, migration, and proliferation of 4T1 cells and induced apoptosis. These results provide structural insights that may facilitate the development of new inhibitors targeting DHOase, for further clinical anticancer chemotherapies.     

Evaluation of Interactions of Selected Olivacine Derivatives with DNA and Topoisomerase II

Olivacine and ellipticine are model anticancer drugs acting as topoisomerase II inhibitors. Here, we present investigations performed on four olivacine derivatives in light of their antitumor activity. The aim of this study was to identify the best antitumor compound among the four tested olivacine derivatives. The study was performed using CCRF/CEM and MCF-7 cell lines. Comet assay, polarography, inhibition of topoisomerase II activity, histone acetylation, and molecular docking studies were performed. Each tested compound displayed interaction with DNA and topoisomerase II, but did not cause histone acetylation. Compound 2 (9-methoxy-5,6-dimethyl-1-({[1-hydroxy-2-(hydroxymethyl)butan-2-yl]amino}methyl)-6H-pyrido[4,3-b]carbazole) was found to be the best candidate as an anticancer drug because it had the highest affinity for topoisomerase II and caused the least genotoxic damage in cells.     

Identification of KPNB1 as a Cellular Target of Aminothiazole Derivatives with Anticancer Activity

We found that aminothiazole derivative (E)‐N‐(5‐benzylthiazol‐2‐yl)‐3‐(furan‐2‐yl)acrylamide ( 1 ) has strong anticancer activity, and undertook proteomics approaches to identify the target protein of compound 1 , importin β1 (KPNB1). A competitive binding assay using fluorescein‐labeled 1 showed that 1 has strong binding affinity for KPNB1 (K_d: ～20 nm ). Furthermore, through western blotting assays for KPNB1, KPNA2, EGFR, ErbB2, and STAT3, we confirmed that 1 has inhibitory effects on the importin pathway. KPBN1 appears to be overexpressed in several cancer cells, and siRNA‐induced inhibition of KPNB1 shows significant inhibition of cancer cell proliferation, while leaving non‐cancerous cells unaffected. Therefore, compound 1 is a promising new lead for the development of KPNB1‐targeted anticancer agents. Fluorescein‐labeled 1 could be a useful quantitative probe for the development of novel KPNB1 inhibitors.     

Studies on [5,6]-Fused Bicyclic Scaffolds Derivatives as Potent Dual B-RafV600E/KDR Inhibitors Using Docking and 3D-QSAR Approaches

Research and development of multi-target inhibitors has attracted increasing attention as anticancer therapeutics. B-Raf^V600E synergistically works with vascular endothelial growth factor receptor 2 (KDR) to promote the occurrence and progression of cancers, and the development of dual-target drugs simultaneously against these two kinds of kinase may offer a better treatment advantage. In this paper, docking and three-dimensional quantitative structure activity relationship (3D-QSAR) studies were performed on a series of dual B-Raf/KDR inhibitors with a novel hinge-binding group, [5,6]-fused bicyclic scaffold. Docking studies revealed optimal binding conformations of these compounds interacting with both B-Raf and KDR. Based on these conformations, comparative molecular field analysis (CoMFA) and comparative molecular similarity indices analysis (CoMSIA) 3D-QSAR models were constructed, and the best CoMFA (q² = 0.542, r² = 0.989 for B-Raf; q² = 0.768, r² = 0.991 for KDR) and CoMSIA models (q² = 0.519, r² = 0.992 for B-Raf; q² = 0.849, r² = 0.993 for KDR) were generated. Further external validations confirmed their predictability, yielding satisfactory correlation coefficients (r²_pred = 0.764 (CoMFA), r²_pred = 0.841 (CoMSIA) for B-Raf, r²_pred = 0.912 (CoMFA), r²_pred = 0.846 (CoMSIA) for KDR, respectively). Through graphical analysis and comparison on docking results and 3D-QSAR contour maps, key amino acids that affect the ligand-receptor interactions were identified and structural features influencing the activities were discussed. New potent derivatives were designed, and subjected to preliminary pharmacological evaluation. The study may offer useful references for the modification and development of novel dual B-Raf/KDR inhibitors.     

Insight into the Structural Determinants of Imidazole Scaffold-Based Derivatives as TNF-α Release Inhibitors by in Silico Explorations

Presently, 151 widely-diverse pyridinylimidazole-based compounds that show inhibitory activities at the TNF-α release were investigated. By using the distance comparison technique (DISCOtech), comparative molecular field analysis (CoMFA), and comparative molecular similarity index analysis (CoMSIA) methods, the pharmacophore models and the three-dimensional quantitative structure-activity relationships (3D-QSAR) of the compounds were explored. The proposed pharmacophore model, including two hydrophobic sites, two aromatic centers, two H-bond donor atoms, two H-bond acceptor atoms, and two H-bond donor sites characterizes the necessary structural features of TNF-α release inhibitors. Both the resultant CoMFA and CoMSIA models exhibited satisfactory predictability (with Q² (cross-validated correlation coefficient) = 0.557, R²_ncv (non-cross-validated correlation coefficient) = 0.740, R²_pre (predicted correlation coefficient) = 0.749 and Q² = 0.598, R²_ncv = 0.767, R²_pre = 0.860, respectively). Good consistency was observed between the 3D-QSAR models and the pharmacophore model that the hydrophobic interaction and hydrogen bonds play crucial roles in the mechanism of actions. The corresponding contour maps generated by these models provide more diverse information about the key intermolecular interactions of inhibitors with the surrounding environment. All these models have extended the understanding of imidazole-based compounds in the structure-activity relationship, and are useful for rational design and screening of novel 2-thioimidazole-based TNF-α release inhibitors.     

Combination of Inhibitors of USP7 and PLK1 has a Strong Synergism against Paclitaxel Resistance

USP7 is a promising target for the development of cancer treatments because of its high expression and the critical functions of its substrates in carcinogenesis of several different carcinomas. Here, we demonstrated the effectiveness of targeting USP7 in advanced malignant cells showing high levels of USP7, especially in taxane-resistant cancer. USP7 knockdown effectively induced cell death in several cancer cells of lung, prostate, and cervix. Depletion of USP7 induced multiple spindle pole formation in mitosis, and, consequently, resulted in mitotic catastrophe. When USP7 was blocked in the paclitaxel-resistant lung cancer NCI-H460^TXR cells, which has resistance to mitotic catastrophe, NCI-H460^TXR cells underwent apoptosis effectively. Furthermore, combination treatment with the mitotic kinase PLK1 inhibitor volasertib and the USP7 inhibitor {"type":"entrez-protein","attrs":{"text":"P22077","term_id":"134707"}}P22077 showed a strong synergism through down-regulation of MDR1/ABCB1 in paclitaxel-resistant lung cancer. Therefore, we suggest USP7 is a promising target for cancer therapy, and combination therapy with inhibitors of PLK1 and USP7 may be valuable for treating paclitaxel-resistant cancers, because of their strong synergism.     

靛玉红类衍生物的三维定量构效关系研究

采用比较分子场方法(CoMFA)对16个靛玉红类衍生物进行三维定量构效关系(3D-QSAR)研究,探讨靛玉红类衍生物的抗癌活性与其结构的关系。建立了3D-QSAR的CoMFA模型,得到了非交叉验证相关系数r2=0.987,标准偏差SD=2.141,统计方差比F=117.136,交叉验证相关系数q2=0.496。该模型合理、可信,具有一定的预测能力,表明改变取代基R1基团的体积才是提高靛玉红类衍生物抗癌活性的主要途径,为开发活性更好的抗癌药物提供一定的理论参考。     

Synthetic Chalcones with Potent Antioxidant Ability on H2O2-Induced Apoptosis in PC12 Cells

Chalcone derivatives (E)-3-(4-hydroxy-3-methoxyphenyl)-1-(4-methoxyphenyl) prop-2-en-1-one and (E)-3-(4-hydroxyphenyl)-1-(4-methoxyphenyl) prop-2-en-1-one (Compounds 1 and 2) have been demonstrated to be potent anti-inflammatory agents in our previous study. In light of the relationship of intracellular mechanisms between anti-inflammatories and antioxidants, we further designed and synthesized a series of chalcone derivatives based on 1 and 2, to explore their antioxidant efficacy. The majority of the derivatives exhibited strong protective effects on PC12 (PC12 rat pheochromocytoma) cells exposed to H₂O₂, and all compounds were nontoxic. A preliminary structure-activity relationship was proposed. Compounds 1 and 1d ((E)-2-methoxy-4-(3-(4-methoxyphenyl)-3-oxoprop-1-en-1-yl) phenyl acrylate) exerted the action in a good dose-dependent manner. Quantitative RT-PCR (qRT-PCR) and western blot analysis showed that 1 and 1d significantly improve the expression of nuclear factor erythroid 2 p45-related factor 2 (Nrf2)-dependent antioxidant genes g-Glutamylcysteine Ligase Catalytic Subunit (GCLC) and heme oxygenase-1 (HO-1) and their corresponding proteins (γ-glutamyl cysteine synthase (γ-GCS) and HO-1) in PC12 cells. Inhibition of GCLC and HO-1 by specific inhibitors, l-buthionine-S-sulfoximine (BSO) and zinc protoporphyrin (ZnPP), respectively, partially reduce the protective effect of 1 and 1d. These data present a series of novel chalcone analogs, especially compounds 1 and 1d, as candidates for treating oxidative stress-related disease by activating the Nrf2-antioxidant responsive element (ARE) pathway.     

Pyrrolo[3,2-d]pyrimidine Derivatives as Type II Kinase Insert Domain Receptor (KDR) Inhibitors: CoMFA and CoMSIA Studies

Kinase insert domain receptor (KDR) inhibitors have been proved to be very effective anticancer agents. Molecular docking, 3D-QSAR methods, CoMFA and CoMSIA were performed on pyrrolo[3,2-d]pyrimidine derivatives as non-ATP competitive KDR inhibitors (type II). The bioactive conformation was explored by docking one potent compound 20 into the active site of KDR in its DFG-out inactive conformation. The constructed CoMFA and CoMSIA models produced statistically significant results with the cross-validated correlation coefficients q(2) of 0.542 and 0.552, non-cross-validated correlation coefficients r(2) of 0.912 and 0.955, and predicted correction coefficients r(2) (pred) of 0.913 and 0.897, respectively. These results ensure the CoMFA and CoMSIA models as a tool to guide the design of a series of new potent KDR inhibitors.     

CCNE1 Is a Putative Therapeutic Target for ARID1A-Mutated Ovarian Clear Cell Carcinoma

Background: Ovarian clear cell carcinoma (OCCC) is resistant to platinum chemotherapy and is characterized by poor prognosis. Today, the use of poly (ADP-ribose) polymerase (PARP) inhibitor, which is based on synthetic lethality strategy and characterized by cancer selectivity, is widely used for new types of molecular-targeted treatment of relapsed platinum-sensitive ovarian cancer. However, it is less effective against OCCC. Methods: We conducted siRNA screening to identify synthetic lethal candidates for the ARID1A mutation; as a result, we identified Cyclin-E1 (CCNE1) as a potential target that affects cell viability. To further clarify the effects of CCNE1, human OCCC cell lines, namely TOV-21G and KOC7c (ARID1A mutant lines), and RMG-I and ES2 (ARID1A wild type lines) were transfected with siRNA targeting CCNE1 or a control vector. Results: Loss of CCNE1 reduced proliferation of the TOV-21G and KOC7c cells but not of the RMG-I and ES2 cells. Furthermore, in vivo interference of CCNE1 effectively inhibited tumor cell proliferation in a xenograft mouse model. Conclusion: This study showed for the first time that CCNE1 is a synthetic lethal target gene to ARID1A-mutated OCCC. Targeting this gene may represent a putative, novel, anticancer strategy in OCCC treatment.     

Design,Synthesis, and Biological Evaluation of Orally Bioavailable CHK1 Inhibitors Active against Acute Myeloid Leukemia

Checkpoint kinase 1 (CHK1) is a central component in DNA damage response and has emerged as a target for antitumor therapeutics. Herein, we describe the design, synthesis, and biological evaluation of a novel series of potent diaminopyrimidine CHK1 inhibitors. The compounds exhibited moderate to potent CHK1 inhibition and could suppress the proliferation of malignant hematological cell lines. The optimized compound 13 had a CHK1 IC₅₀ value of 7.73±0.74 nM, and MV-4-11 cells were sensitive to it (IC₅₀=0.035±0.007 μM). Furthermore, compound 13 was metabolically stable in mouse liver microsomes in vitro and displayed moderate oral bioavailability in vivo. Moreover, treatment of MV-4-11 cells with compound 13 for 2 h led to robust inhibition of CHK1 autophosphorylation on serine 296. Based on these biochemical results, we consider compound 13 to be a promising CHK1 inhibitor and potential anticancer therapeutic agent.