A Back-to-Front Fragment-Based Drug Design Search Strategy Targeting the DFG-Out Pocket of Protein Tyrosine Kinases

X-ray crystallographic analyses of 14 IPMK inhibitor complexes

Inositol polyphosphate multikinase (IPMK) is a ubiquitously expressed kinase that has been linked to several cancers. Here, we report 14 new co-crystal structures (1.7Å - 2.0Å resolution) of human IPMK complexed with various IPMK inhibitors developed by another group. The new structures reveal two ordered water molecules that participate in hydrogen-bonding networks, and an unoccupied pocket in the ATP-binding site of human IPMK. New Protein Data Bank (PDB) codes of these IPMK crystal structures are: 8V6W (1.95Å), 8V6X (1.75Å), 8V6Y (1.70Å), 8V6Z (1.85Å), 8V70 (1.85Å), 8V71 (1.70Å), 8V72 (2.0Å), 8V73 (1.90Å), 8V74 (1.85Å), 8V75 (1.85Å), 8V76 (1.95Å), 8V77 (1.95Å), 8V78 (1.95Å), 8V79 (1.95Å).

Novel VEGFR-2 inhibitors as antiangiogenic and apoptotic agents via paracrine and autocrine cascades: Design, synthesis, and biological evaluation

Appertaining to its paracrine and autocrine signaling loops, VEGFR-2 succeeded in grabbing attention as one of the leading targets in cancer treatment. Based on the foregoing and our comprehensive studies regarding pharmacophoric features and activity of sorafenib, novel phenylpyridazinone based VEGFR-2 inhibitors 4, 6a-e, 7a,b, 9a,b, 12a-c, 13a,b, 14a,b, 15a,b, and 17a-d were optimized. An assortment of biological assays was conducted to assess the antiangiogenic and apoptotic activities of the synthesized derivatives. In vitro VEGFR-2 kinase assay verified the inhibitory activity of the synthesized derivatives with IC50 values from 49.1 to 418.0 nM relative to the reference drug sorafenib (IC50 = 81.8 nM). Antiproliferative activity against HUVECs revealed that compounds 2-{2-[2-(6-oxo-3-phenylpyridazin-1(6H)-yl)acetyl]hydrazineyl}-N-(p-tolyl)acetamide (12c) and 2-[(5-mercapto-4-methyl-4H-1,2,4-triazol-3-yl)methyl]-6-phenylpyridazin-3(2H)-one (13a) possessed superior activity (IC50 values = 11.5 and 12.3 nM, respectively) in comparison to sorafenib (IC50 = 23.2 nM). For the purpose of appraising their antiproliferative effect, derivatives 12c and 13a were exposed to cell cycle analysis, apoptotic, cell invasion and migration assays in addition to determination of VEGFR-2 in protein level. Moreover, cytotoxicity as well as selectivity index against WI-38 cell line was measured to examine safety of derivatives 12c and 13a. After that, molecular docking study was executed on the top five compounds in the in vitro VEGFR-2 kinase assay 6d, 12c, 13a, 14a and 17c to get a deep perception on binding mode of the synthesized compounds and correlate the design strategy with biological results. Finally, physicochemical, pharmacokinetic properties, and drug-likeness studies were performed on the top five derivative in in vitro VEGFR-2 kinase assay.

Pharmacophore Modeling of Targets Infested with Activity Cliffs via Molecular Dynamics Simulation Coupled with QSAR and Comparison with other Pharmacophore Generation Methods: KDR as Case Study

Activity cliffs (ACs) are defined as pairs of structurally similar compounds with large difference in their potencies against certain biotarget. We recently proposed that potent AC members induce significant entropically-driven conformational modifications of the target that unveil additional binding interactions, while their weakly-potent counterparts are enthalpically-driven binders with little influence on the protein target. We herein propose to extract pharmacophores for ACs-infested target(s) from molecular dynamics (MD) frames of purely "enthalpic" potent binder(s) complexed within the particular target. Genetic function algorithm/machine learning (GFA/ML) can then be employed to search for the best possible combination of MD pharmacophore(s) capable of explaining bioactivity variations within a list of inhibitors. We compared the performance of this approach with established ligand-based and structure-based methods. Kinase inserts domain receptor (KDR) was used as a case study. KDR plays a crucial role in angiogenic signaling and its inhibitors have been approved in cancer treatment. Interestingly, GFA/ML selected, MD-based, pharmacophores were of comparable performances to ligand-based and structure-based pharmacophores. The resulting pharmacophores and QSAR models were used to capture hits from the national cancer institute list of compounds. The most active hit showed anti-KDR IC50 of 2.76 µM.

Novel 1,3-Thiazole Analogues with Potent Activity against Breast Cancer: A Design, Synthesis, In Vitro, and In Silico Study

Breast cancer is the most common cancer in women, responsible for over half a million deaths in 2020. Almost 75% of FDA-approved drugs are mainly nitrogen- and sulfur-containing heterocyclic compounds, implying the importance of such compounds in drug discovery. Among heterocycles, thiazole-based heterocyclic compounds have demonstrated a broad range of pharmacological activities. In the present study, a novel set of 1,3-thiazole derivatives was designed and synthesized based on the coupling of acetophenone derivatives, and phenacyl bromide was substituted as a key reaction step. The activity of synthesized compounds was screened against the proliferation of two breast cancer cell lines (MCF-7 and MDA-MB-231). Almost all compounds exhibited a considerable antiproliferative activity toward the breast cancer cells as compared to staurosporine, with no significant cytotoxicity toward the epithelial cells. Among the synthesized compounds, compound 4 exhibited the most potent antiproliferative activity, with an IC₅₀ of 5.73 and 12.15 µM toward MCF-7 and MDA-MB-231 cells, respectively, compared to staurosporine (IC₅₀ = 6.77 and 7.03 µM, respectively). Exploring the mechanistic insights responsible for the antiproliferative activity of compound 4 revealed that compound 4 possesses a significant inhibitory activity toward the vascular endothelial growth factor receptor-2 (VEGFR-2) with (IC₅₀ = 0.093 µM) compared to Sorafenib (IC₅₀ = 0.059 µM). Further, compound 4 showed the ability to induce programmed cell death by triggering apoptosis and necrosis in MCF-7 cells and to induce cell cycle arrest on MCF-7 cells at the G1 stage while decreasing the cellular population in the G2/M phase. Finally, detailed in silico molecular docking studies affirmed that this class of compounds possesses a considerable binding affinity toward VEGFR2 proteins. Overall, these results indicate that compound 4 could be a promising lead compound for developing potent anti-breast cancer compounds.

Fragment Hit Discovery and Binding Site Characterization by Indirect Affinity Capillary Electrophoresis: Application to Factor XIIa

Fragment-based lead discovery is a usual strategy in drug discovery to identify innovative lead compounds. The success of this approach strongly relies on the capacity to detect weak binders and characterize their binding site. NMR and X-ray crystallography are the conventional technologies used to tackle this challenge. However, their large protein consumption and the cost of equipment reduce their accessibility. Here, an affinity capillary electrophoresis methodology was developed that enables the detection of mM binders, the determination of dissociation constants, and the characterization of the fragment binding site. On the basis of multiple equilibrium theory, dissociation constants in the μM-mM range were determined, and a new methodology is proposed to establish graphically if two fragments bind the same protein pocket. The applicability of this methodology was demonstrated experimentally on coagulation factor XIIa by evaluating pairs of fragments with expected behavior. This study reinforces the significance of using affinity capillary electrophoresis to gather valuable information for medicinal chemistry projects.

Design and synthesis of novel furan, furo[2,3-d]pyrimidine and furo[3,2-e][1,2,4]triazolo[1,5-c]pyrimidine derivatives as potential VEGFR-2 inhibitors

Novel furan 6a-c, furo[2,3-d]pyrimidine 7a-f, 9, 10a-f, 12a,b, 14a-d and furo[3,2-e][1,2,4]triazolo[1,5-c]pyrimidine 8a-f derivatives were designed based on their structural similarity to a previously described oxazole VEGFR-2 back pocket binding fragment. The designed compounds were synthesized and screened for their in vitro VEGFR-2 inhibitory activity where they exhibited good to moderate nanomolar inhibition with improved ligand efficiencies. 8b and 10c (IC₅₀ = 38.72 ± 1.7 and 41.40 ± 1.8 nM, respectively) were equipotent to sorafenib and 6a, 6c, 7f, 8a, 8c, 10b, 10f, 12b, 14c and 14d showed good activity (IC₅₀ = 43.31-98.31 nM). The furotriazolopyrimidines 8a-c and furopyrimidine derivative 10c were further evaluated for their in vitro antiproliferative activity against human umbilical vein endothelial cells (HUVECs) where 8b showed higher potency than sorafenib and resulted in cell cycle arrest at G2/M phase whereas 8c revealed good antiproliferative activity with cell cycle arrest at G1 phase. Moreover, 8a-c and 10c showed significant inhibitory effects on the invasion and migration of HUVECs. Molecular docking study was conducted to gain insight about the potential binding mode. The furo[3,2-e][1,2,4]triazolo[1,5-c]pyrimidine derivatives 8b and 8c represent interesting starting point for antiangiogenic compounds based on their activity and favorable drug likeness profiles.

Synthesis and biological evaluation of 3,5-diaryl-pyrazole derivatives as potential antiprostate cancer agents

Prostate cancer is the most frequently diagnosed tumor in men and the second leading cause of cancer-associated mortality in most developed countries. 3,5-Diaryl substituted pyrazole derivatives (20-28) were prepared starting from related chalcones and biologically evaluated for in vitro growth inhibition activity against PC3 and DU145 human prostate cancer cell lines. Compounds 23, 26, and 28 were found to be more potent as compared to the other halogen-substituted derivatives. Especially, the 2-bromo-substituted pyrazole derivative (26) was found to be more potent against PC3 and DU145 cells. Epidermal growth factor receptor (EGFR) and vascular endothelial growth factor receptor 2 (VEGFR2) are known to be expressed in DU145 and PC3 cancer cells. The binding mode of the most selective compound 26 toward EGFR and VEGFR2 was investigated by employing docking simulations based on GLIDE standard precision (-5.912 and -6.949 kcal/mol, respectively).

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.

Glechomanamides A-C, Germacrane Sesquiterpenoids with an Unusual Δ8-7,12-Lactam Moiety from Salvia scapiformis and Their Antiangiogenic Activity

Three new germacrane sesquiterpenoid-type alkaloids with an unusual Δ⁸-7,12-lactam moiety, glechomanamides A-C (1-3), and two pairs of 7,12-hemiketal sesquiterpenoid epimers (4a/b, 5a/b) were isolated from Salvia scapiformis. Their structures were elucidated by spectroscopic methods including HRESIMS, IR, UV, and 1D and 2D NMR and also confirmed by single-crystal X-ray diffraction analysis. The chemical transformation of compounds 1-5 in a solution environment was analyzed by 2D NMR spectroscopy. The aza acetallactams (1-3) were stable in organic solvent, while single crystals of the hemiacetal esters (4a/b, 5a/b) underwent a tautomeric equilibrium after being dissolved. Single crystals of 4a, 4b, and 5a were obtained for the first time as their naturally occurring forms. Glechomanamide B (2) exhibited antiangiogenic activity by suppression of vascular endothelial growth factor (VEGF)-induced tube formation through modulation of VEGF receptor 2 (VEGFR2)-mediated signaling pathways in human umbilical vascular endothelial cells (HUVECs). In addition, compound 2 also showed the significant suppression of mRNA expression associated with glycolysis and angiogenesis biomarkers in high glucose (30 mM)-induced HUVECs. These findings suggest that compound 2 might be a potential lead compound candidate for the management of diabetic retinopathy.

Discovery of Small Molecules that Target Vascular Endothelial Growth Factor Receptor-2 Signalling Pathway Employing Molecular Modelling Studies

Angiogenesis is defined as the formation of new blood vessels and is a key phenomenon manifested in a host of cancers during which tyrosine kinases play a crucial role. Vascular endothelial growth factor receptor-2 (VEGFR-2) is pivotal in cancer angiogenesis, which warrants the urgency of discovering new anti-angiogenic inhibitors that target the signalling pathways. To obtain this objective, a structure-based pharmacophore model was built from the drug target VEGFR-2 (PDB code: 4AG8), complexed with axitinib and was subsequently validated and employed as a 3D query to retrieve the candidate compounds with the key inhibitory features. The model was escalated to molecular docking studies resulting in seven candidate compounds. The molecular docking studies revealed that the seven compounds displayed a higher dock score than the reference-cocrystallised compound. The GROningen MAchine for Chemical Simulations (GROMACS) package guided molecular dynamics (MD) results determined their binding mode and affirmed stable root mean square deviation. Furthermore, these compounds have preserved their key interactions with the residues Glu885, Glu917, Cys919 and Asp1046. The obtained findings deem that the seven compounds could act as novel anti-angiogenic inhibitors and may further assist as the prototype in designing and developing new inhibitors.

Design, Synthesis and Biological Activity Evaluation of S-Substituted 1H-5-Mercapto-1,2,4-Triazole Derivatives as Antiproliferative Agents in Colorectal Cancer

Colon cancer is a widespread pathology with complex biochemical etiology based on a significant number of intracellular signaling pathways that play important roles in carcinogenesis, tumor proliferation and metastasis. These pathways function due to the action of key enzymes that can be used as targets for new anticancer drug development. Herein we report the synthesis and biological antiproliferative evaluation of a series of novel S-substituted 1H-3-R-5-mercapto-1,2,4-triazoles, on a colorectal cancer cell line, HT-29. Synthesized compounds were designed by docking based virtual screening (DBVS) of a previous constructed compound library against protein targets, known for their important role in colorectal cancer signaling: MEK1, ERK2, PDK1, VEGFR2. Among all synthesized structures, TZ55.7, which was retained as a possible PDK1 (phospholipid-dependent kinase 1) inhibitor, exhibited the most significant cytotoxic activity against HT-29 tumor cell line. The same compound alongside other two, TZ53.7 and TZ3a.7, led to a significant cell cycle arrest in both sub G0/G1 and G0/G1 phase. This study provides future perspectives for the development of new agents containing the 1,2,4-mercapto triazole scaffold with antiproliferative activities in colorectal cancer.

Exploration for novel inhibitors showing back-to-front approach against VEGFR-2 kinase domain (4AG8) employing molecular docking mechanism and molecular dynamics simulations

Background

Angiogenesis is a process of formation of new blood vessels and is an important criteria demonstrated by cancer cells. Over a period of time, these cancer cells infect the other parts of the healthy body by a process called progression. The objective of the present article is to identify a drug molecule that inhibits angiogenesis and progression.

Methods

In this pursuit, ligand based pharmacophore virtual screening was employed, generating a pharmacophore model, Hypo1 consisting of four features. Furthermore, this Hypo1 was validated recruiting, Fischer’s randomization, test set method and decoy set method. Later, Hypo1 was allowed to screen databases such as Maybridge, Chembridge, Asinex and NCI and were further filtered by ADMET filters and Lipinski’s Rule of Five. A total of 699 molecules that passed the above criteria, were challenged against 4AG8, an angiogenic drug target employing GOLD v5.2.2.

Results

The results rendered by molecular docking, DFT and the MD simulations showed only one molecule (Hit) obeyed the back-to-front approach. This molecule displayed a dock score of 89.77, involving the amino acids, Glu885 and Cys919, Asp1046, respectively and additionally formed several important hydrophobic interactions. Furthermore, the identified lead molecule showed interactions with key residues when challenged with CDK2 protein, 1URW.

Conclusion

The lead candidate showed several interactions with the crucial residues of both the targets. Furthermore, we speculate that the residues Cys919 and Leu83 are important in the development of dual inhibitor. Therefore, the identified lead molecule can act as a potential inhibitor for angiogenesis and progression.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4050-1) contains supplementary material, which is available to authorized users.

Eriocalyxin B, a natural diterpenoid, inhibited VEGF-induced angiogenesis and diminished angiogenesis-dependent breast tumor growth by suppressing VEGFR-2 signaling

Eriocalyxin B (EriB), a natural ent-kaurane diterpenoid isolated from the plant Isodon eriocalyx var. laxiflora, has emerged as a promising anticancer agent. The effects of EriB on angiogenesis were explored in the present study. Here we demonstrated that the subintestinal vein formation was significantly inhibited by EriB treatment (10, 15 μM) in zebrafish embryos, which was resulted from the alteration of various angiogenic genes as shown in transcriptome profiling. In human umbilical vein endothelial cells, EriB treatment (50, 100 nM) could significantly block vascular endothelial growth factors (VEGF)-induced cell proliferation, tube formation, cell migration and cell invasion. Furthermore, EriB also caused G1 phase cell cycle arrest which was correlated with the down-regulation of the cyclin D1 and CDK4 leading to the inhibition of phosphorylated retinoblastoma protein expression. Investigation of the signal transduction revealed that EriB inhibited VEGF-induced phosphorylation of VEGF receptor-2 via the interaction with the ATP-binding sites according to the molecular docking simulations. The suppression of VEGFR-2 downstream signal transduction cascades was also observed. EriB was showed to inhibit new blood vessel formation in Matrigel plug model and mouse 4T1 breast tumor model. EriB (5 mg/kg/day) treatment was able to decrease tumor vascularization and suppress tumor growth and angiogenesis. Taken together, our findings suggested that EriB is a novel inhibitor of angiogenesis through modulating VEGFR-2 signaling pathway, which could be developed as a promising anti-angiogenic agent for treatment of angiogenesis-related human diseases, such as cancer.

Fragment-Based Discovery of Potent and Selective DDR1/2 Inhibitors

The DDR1 and DDR2 receptor tyrosine kinases are activated by extracellular collagen and have been implicated in a number of human diseases including cancer. We performed a fragment-based screen against DDR1 and identified fragments that bound either at the hinge or in the back pocket associated with the DFG-out conformation of the kinase. Modeling based on crystal structures of potent kinase inhibitors facilitated the "back-to-front" design of potent DDR1/2 inhibitors that incorporated one of the DFG-out fragments. Further optimization led to low nanomolar, orally bioavailable inhibitors that were selective for DDR1 and DDR2. The inhibitors were shown to potently inhibit DDR2 activity in cells but in contrast to unselective inhibitors such as dasatinib, they did not inhibit proliferation of mutant DDR2 lung SCC cell lines.

Novel approaches for targeting kinases: allosteric inhibition, allosteric activation and pseudokinases

Protein kinases are involved in many essential cellular processes and their deregulation can lead to a variety of diseases, including cancer. The pharmaceutical industry has invested heavily in the identification of kinase inhibitors to modulate these disease-promoting pathways, resulting in several successful drugs. However, the field is challenging as it is difficult to identify novel selective inhibitors with good pharmacokinetic/pharmacodynamic properties. In addition, resistance to kinase inhibitor treatment frequently arises. The identification of non-ATP site targeting ('allosteric') inhibitors, the identification of kinase activators and the expansion of kinase target space to include the less studied members of the family, including atypical- and pseudo-kinases, are potential avenues to overcome these challenges. In this perspective, the opportunities and challenges of following these approaches and others will be discussed.

Optimization of potent DFG-in inhibitors of platelet derived growth factor receptorβ (PDGF-Rβ) guided by water thermodynamics

In this study we report on the hit optimization of substituted 3,5-diaryl-pyrazin-2(1H)-ones toward potent and effective platelet-derived growth factor receptor (PDGF-R) β-inhibitors. Originally, the 3,5-diaryl-pyrazin-2-one core was derived from the marine sponge alkaloid family of hamacanthins. In our first series compound 2 was discovered as a promising hit showing strong activity against PDGF-Rβ in the kinase assay (IC50 = 0.5 μM). Furthermore, 2 was shown to be selective for PDGF-Rβ in a panel of 24 therapeutically relevant protein kinases. Molecular modeling studies on a PDGF-Rβ homology model using prediction of water thermodynamics suggested an optimization strategy for the 3,5-diaryl-pyrazin-2-ones as DFG-in binders by using a phenolic OH function to replace a structural water molecule in the ATP binding site. Indeed, we identified compound 38 as a highly potent inhibitor with an IC50 value of 0.02 μM in a PDGF-Rβ enzymatic assay also showing activity against PDGF-R dependent cancer cells.

Fragment-based approaches to the discovery of kinase inhibitors

Protein kinases are one of the most important families of drug targets, and aberrant kinase activity has been linked to a large number of disease areas. Although eminently targetable using small molecules, kinases present a number of challenges as drug targets, not least obtaining selectivity across such a large and relatively closely related target family. Fragment-based drug discovery involves screening simple, low-molecular weight compounds to generate initial hits against a target. These hits are then optimized to more potent compounds via medicinal chemistry, usually facilitated by structural biology. Here, we will present a number of recent examples of fragment-based approaches to the discovery of kinase inhibitors, detailing the construction of fragment-screening libraries, the identification and validation of fragment hits, and their optimization into potent and selective lead compounds. The advantages of fragment-based methodologies will be discussed, along with some of the challenges associated with using this route. Finally, we will present a number of key lessons derived both from our own experience running fragment screens against kinases and from a large number of published studies.

Novel VEGFR-2 kinase inhibitors identified by the back-to-front approach

We report a novel VEGFR-2 inhibitor, developed by the back-to-front approach. Docking experiments indicated that the 3-chloromethylphenylurea motif of the lead compound occupied the back pocket of VEGFR-2 kinase. An attempt was made to enhance the binding affinity of 1 by expanding the structure to access the front pocket using a triazole linker. A library of 1,4-(disubstituted)-1H-1,2,3-triazoles were screened in silico, and one compound (VH02) was identified with an IC50 against VEGFR-2 of 0.56μM. VH02 showed antiangiogenic effects, inhibiting tube formation in HUVEC cells (EA.hy926) at 0.3μM, 13 times lower than its cytotoxic dose. These enzymatic and cellular activities suggest that VH02 has potential as a lead for further optimization.