Structure-based virtual screening approach to identify novel classes of Cdc25B phosphatase inhibitors

A Comprehensive Overview of the Developments of Cdc25 Phosphatase Inhibitors

Cdc25 phosphatases have been considered promising targets for anticancer development due to the correlation of their overexpression with a wide variety of cancers. In the last two decades, the interest in this subject has considerably increased and many publications have been launched concerning this issue. An overview is constructed based on data analysis of the results of the previous publications covering the years from 1992 to 2021. Thus, the main objective of the current review is to report the chemical structures of Cdc25s inhibitors and answer the question, how to design an inhibitor with better efficacy and lower toxicity?

Medicinal chemistry insights into novel CDC25 inhibitors

Cell division cycle 25 (CDC25) phosphatases, a kind of cell cycle regulators, have become an attractive target for drug discovery, as they have been found to be over-expressed in various human cancer cells. Several CDC25 inhibitors have achieved significant attention in clinical trials with possible mechanistic actions. Prompted by the significance of CDC25 inhibitors with medicinal chemistry prospect, it is an apt time to review the various drug discovery methods involved in CDC25 drug discovery including high throughput screening (HTS), virtual screening (VS), fragment-based drug design, substitution decorating approach, structural simplification approach and scaffold hopping method to seek trends and identify promising new avenues of CDC25 drug discovery.

CDC25 Inhibition in Acute Myeloid Leukemia-A Study of Patient Heterogeneity and the Effects of Different Inhibitors

Cell division cycle 25 (CDC25) protein phosphatases regulate cell cycle progression through the activation of cyclin-dependent kinases (CDKs), but they are also involved in chromatin modulation and transcriptional regulation. CDC25 inhibition is regarded as a possible therapeutic strategy for the treatment of human malignancies, including acute myeloid leukemia (AML). We investigated the in vitro effects of CDC25 inhibitors on primary human AML cells derived from 79 unselected patients in suspension cultures. Both the previously well-characterized CDC25 inhibitor NSC95397, as well as five other inhibitors (BN82002 and the novel small molecular compounds ALX1, ALX2, ALX3, and ALX4), only exhibited antiproliferative effects for a subset of patients when tested alone. These antiproliferative effects showed associations with differences in genetic abnormalities and/or AML cell differentiation. However, the responders to CDC25 inhibition could be identified by analysis of global gene expression profiles. The differentially expressed genes were associated with the cytoskeleton, microtubules, and cell signaling. The constitutive release of 28 soluble mediators showed a wide variation among patients and this variation was maintained in the presence of CDC25 inhibition. Finally, NSC95397 had no or only minimal effects on AML cell viability. In conclusion, CDC25 inhibition has antiproliferative effects on primary human AML cells for a subset of patients, and these patients can be identified by gene expression profiling.

Insights into the interaction of high potency inhibitor IRC-083864 with phosphatase CDC25

CDC25 phosphatases play a crucial role in cell cycle regulation. They have been found to be over-expressed in various human tumours and to be valuable targets for cancer treatment. Here, we report the first model of binding of the most potent CDC25 inhibitor to date, the bis-quinone IRC-083864, into CDC25B obtained by combining molecular modeling and NMR studies. Our study provides new insights into key interactions of the catalytic site inhibitor and CDC25B in the absence of any available experimental structure of CDC25 with a bound catalytic site inhibitor. The docking model reveals that IRC-083864 occupies both the active site and the inhibitor binding pocket of the CDC25B catalytic domain. NMR saturation transfer difference and WaterLOGSY data indicate the binding zones of the inhibitor and support the docking model. Probing interactions of analogues of the two quinone units of IRC-083864 with CDC25B demonstrate that IRC-083864 competes with each monomer. Proteins 2017; 85:593-601. © 2016 Wiley Periodicals, Inc.

Druggability analysis and classification of protein tyrosine phosphatase active sites

Protein tyrosine phosphatases (PTP) play important roles in the pathogenesis of many diseases. The fact that no PTP inhibitors have reached the market so far has raised many questions about their druggability. In this study, the active sites of 17 PTPs were characterized and assessed for its ability to bind drug-like molecules. Consequently, PTPs were classified according to their druggability scores into four main categories. Only four members showed intermediate to very druggable pocket; interestingly, the rest of them exhibited poor druggability. Particularly focusing on PTP1B, we also demonstrated the influence of several factors on the druggability of PTP active site. For instance, the open conformation showed better druggability than the closed conformation, while the tight-bound water molecules appeared to have minimal effect on the PTP1B druggability. Finally, the allosteric site of PTP1B was found to exhibit superior druggability compared to the catalytic pocket. This analysis can prove useful in the discovery of new PTP inhibitors by assisting researchers in predicting hit rates from high throughput or virtual screening and saving unnecessary cost, time, and efforts via prioritizing PTP targets according to their predicted druggability.

Ligand-based chemoinformatic discovery of a novel small molecule inhibitor targeting CDC25 dual specificity phosphatases and displaying in vitro efficacy against melanoma cells

CDC25 phosphatases are important regulators of the cell cycle and represent promising targets for anticancer drug discovery. We recently identified NSC 119915 as a new quinonoid CDC25 inhibitor with potent anticancer activity. In order to discover more active analogs of NSC 119915, we performed a range of ligand-based chemoinformatic methods against the full ZINC drug-like subset and the NCI lead-like set. Nine compounds (3, 5-9, 21, 24, and 25) were identified with Ki values for CDC25A, -B and -C ranging from 0.01 to 4.4 μM. One of these analogs, 7, showed a high antiproliferative effect on human melanoma cell lines, A2058 and SAN. Compound 7 arrested melanoma cells in G2/M, causing a reduction of the protein levels of CDC25A and, more consistently, of CDC25C. Furthermore, an intrinsic apoptotic pathway was induced, which was mediated by ROS, because it was reverted in the presence of antioxidant N-acetyl-cysteine (NAC). Finally, 7 decreased the protein levels of phosphorylated Akt and increased those of p53, thus contributing to the regulation of chemosensitivity through the control of downstream Akt pathways in melanoma cells. Taken together, our data emphasize that CDC25 could be considered as a possible oncotarget in melanoma cells and that compound 7 is a small molecule CDC25 inhibitor that merits to be further evaluated as a chemotherapeutic agent for melanoma, likely in combination with other therapeutic compounds.

GRID and docking analyses reveal a molecular basis for flavonoid inhibition of Src family kinase activity

Flavonoids reduce cardiovascular disease risk through anti-inflammatory, anti-coagulant and anti-platelet actions. One key flavonoid inhibitory mechanism is blocking kinase activity that drives these processes. Flavonoids attenuate activities of kinases including phosphoinositide-3-kinase, Fyn, Lyn, Src, Syk, PKC, PIM1/2, ERK, JNK and PKA. X-ray crystallographic analyses of kinase-flavonoid complexes show that flavonoid ring systems and their hydroxyl substitutions are important structural features for their binding to kinases. A clearer understanding of structural interactions of flavonoids with kinases is necessary to allow construction of more potent and selective counterparts. We examined flavonoid (quercetin, apigenin and catechin) interactions with Src family kinases (Lyn, Fyn and Hck) applying the Sybyl docking algorithm and GRID. A homology model (Lyn) was used in our analyses to demonstrate that high-quality predicted kinase structures are suitable for flavonoid computational studies. Our docking results revealed potential hydrogen bond contacts between flavonoid hydroxyls and kinase catalytic site residues. Identification of plausible contacts indicated that quercetin formed the most energetically stable interactions, apigenin lacked hydroxyl groups necessary for important contacts and the non-planar structure of catechin could not support predicted hydrogen bonding patterns. GRID analysis using a hydroxyl functional group supported docking results. Based on these findings, we predicted that quercetin would inhibit activities of Src family kinases with greater potency than apigenin and catechin. We validated this prediction using in vitro kinase assays. We conclude that our study can be used as a basis to construct virtual flavonoid interaction libraries to guide drug discovery using these compounds as molecular templates.

Therapeutic targeting the cell division cycle 25 (CDC25) phosphatases in human acute myeloid leukemia--the possibility to target several kinases through inhibition of the various CDC25 isoforms

The cell division cycle 25 (CDC25) phosphatases include CDC25A, CDC25B and CDC25C. These three molecules are important regulators of several steps in the cell cycle, including the activation of various cyclin-dependent kinases (CDKs). CDC25s seem to have a role in the development of several human malignancies, including acute myeloid leukemia (AML); and CDC25 inhibition is therefore considered as a possible anticancer strategy. Firstly, upregulation of CDC25A can enhance cell proliferation and the expression seems to be controlled through PI3K-Akt-mTOR signaling, a pathway possibly mediating chemoresistance in human AML. Loss of CDC25A is also important for the cell cycle arrest caused by differentiation induction of malignant hematopoietic cells. Secondly, high CDC25B expression is associated with resistance against the antiproliferative effect of PI3K-Akt-mTOR inhibitors in primary human AML cells, and inhibition of this isoform seems to reduce AML cell line proliferation through effects on NFκB and p300. Finally, CDC25C seems important for the phenotype of AML cells at least for a subset of patients. Many of the identified CDC25 inhibitors show cross-reactivity among the three CDC25 isoforms. Thus, by using such cross-reactive inhibitors it may become possible to inhibit several molecular events in the regulation of cell cycle progression and even cytoplasmic signaling, including activation of several CDKs, through the use of a single drug. Such combined strategies will probably be an advantage in human cancer treatment.

Protein tyrosine phosphatases: Ligand interaction analysis and optimisation of virtual screening

Docking-based virtual screening is an established component of structure-based drug discovery. Nevertheless, scoring and ranking of computationally docked ligand libraries still suffer from many false positives. Identifying optimal docking parameters for a target protein prior to virtual screening can improve experimental hit rates. Here, we examine protocols for virtual screening against the important but challenging class of drug target, protein tyrosine phosphatases. In this study, common interaction features were identified from analysis of protein-ligand binding geometries of more than 50 complexed phosphatase crystal structures. It was found that two interactions were consistently formed across all phosphatase inhibitors: (1) a polar contact with the conserved arginine residue, and (2) at least one interaction with the P-loop backbone amide. In order to investigate the significance of these features on phosphatase-ligand binding, a series of seeded virtual screening experiments were conducted on three phosphatase enzymes, PTP1B, Cdc25b and IF2. It was observed that when the conserved arginine and P-loop amide interactions were used as pharmacophoric constraints during docking, enrichment of the virtual screen significantly increased in the three studied phosphatases, by up to a factor of two in some cases. Additionally, the use of such pharmacophoric constraints considerably improved the ability of docking to predict the inhibitor's bound pose, decreasing RMSD to the crystallographic geometry by 43% on average. Constrained docking improved enrichment of screens against both open and closed conformations of PTP1B. Incorporation of an ordered water molecule in PTP1B screening was also found to generally improve enrichment. The knowledge-based computational strategies explored here can potentially inform structure-based design of new phosphatase inhibitors using docking-based virtual screening.

Ligand pose and orientational sampling in molecular docking

Molecular docking remains an important tool for structure-based screening to find new ligands and chemical probes. As docking ambitions grow to include new scoring function terms, and to address ever more targets, the reliability and extendability of the orientation sampling, and the throughput of the method, become pressing. Here we explore sampling techniques that eliminate stochastic behavior in DOCK3.6, allowing us to optimize the method for regularly variable sampling of orientations. This also enabled a focused effort to optimize the code for efficiency, with a three-fold increase in the speed of the program. This, in turn, facilitated extensive testing of the method on the 102 targets, 22,805 ligands and 1,411,214 decoys of the Directory of Useful Decoys - Enhanced (DUD-E) benchmarking set, at multiple levels of sampling. Encouragingly, we observe that as sampling increases from 50 to 500 to 2000 to 5000 to 20000 molecular orientations in the binding site (and so from about 1×10¹⁰ to 4×10¹⁰ to 1×10¹¹ to 2×10¹¹ to 5×10¹¹ mean atoms scored per target, since multiple conformations are sampled per orientation), the enrichment of ligands over decoys monotonically increases for most DUD-E targets. Meanwhile, including internal electrostatics in the evaluation ligand conformational energies, and restricting aromatic hydroxyls to low energy rotamers, further improved enrichment values. Several of the strategies used here to improve the efficiency of the code are broadly applicable in the field.