Computer-aided discovery of new FGFR-1 inhibitors followed by in vitro validation

Synthesis and Biological Evaluation of Thiazole-based Fibroblast Growth Factor Receptor-1 Inhibitors

BACKGROUND

The Fibroblast Growth Factor Receptor-1 (FGFR-1) is a tyrosine kinase and a validated target for the treatment of different cancer types.

OBJECTIVE

Design and synthesis of novel thiazole-based analogues of anticancer agents.

METHODS

Series of 2-aryl-5-methylthiazole analogues linked to structurally variable basic heads were synthesized as novel anticancer agents. Developed compounds were tested for their cytotoxic activities against several cancer cell lines.

RESULTS

Many analogues exhibited strong antiproliferative activities against breast cancer cell lines, with higher potency towards the highly metastatic form (MDA-MB-231). Pharmacophoric profiling using an in-house pharmacophore database identified FGFR-1 as a molecular target of active analogues. Synthesized compounds were bioassayed for their FGFR-1 inhibitory activities and many hits exhibited IC50 values in the low micromolar to nanomolar range.

CONCLUSION

The 2-aryl-5-methylthiazole linked to a basic head is a novel chemical scaffold of ATP-competitive inhibitor of FGFR-1 with potential therapeutic activities against different types of cancer.

Discovery of new potent lysine specific histone demythelase-1 inhibitors (LSD-1) using structure based and ligand based molecular modelling and machine learning

Lysine-specific histone demethylase 1 (LSD-1) is an epigenetic enzyme that oxidatively cleaves methyl groups from monomethyl and dimethyl Lys4 of histone H3 and is highly overexpressed in different types of cancer. Therefore, it has been widely recognized as a promising therapeutic target for cancer therapy. Towards this end, we employed various Computer Aided Drug Design (CADD) approaches including pharmacophore modelling and machine learning. Pharmacophores generated by structure-based (SB) (either crystallographic-based or docking-based) and ligand-based (LB) (either supervised or unsupervised) modelling methods were allowed to compete within the context of genetic algorithm/machine learning and were assessed by Shapley additive explanation values (SHAP) to end up with three successful pharmacophores that were used to screen the National Cancer Institute (NCI) database. Seventy-five NCI hits were tested for their LSD-1 inhibitory properties against neuroblastoma SH-SY5Y cells, pancreatic carcinoma Panc-1 cells, glioblastoma U-87 MG cells and in vitro enzymatic assay, culminating in 3 nanomolar LSD-1 inhibitors of novel chemotypes.

Biological Evaluation and Reverse Pharmacophore Mapping of Innovative Bis-Triazoles as Promising Anticancer Agents

Here, we describe further cytotoxic studies and reverse pharmacophore mapping (pharmacophore profiling) for bis-triazoles MS44-53, which were designed and synthesized previously to stabilize the G-quadruplex nucleic acids capable of being formed at the telomeric region and promoter sequences of genes involved in cellular proliferation and oncogenes. Pharmacophore-based activity profiling screen demonstrated some biological targets that MS44-53 may modulate their biological response, and thus can be considered as potential drugs to treat different kinds of diseases, such as carcinoma, diabetes type II, bacterial infection and cardiovascular diseases. Potent cell growth inhibitory properties were shown by ligands MS47 and MS49 against human melanoma MDA-MB-435, colon cancer HCT-116 and COLO 205, and pancreatic cancer MIA PaCa-2 cell lines, as evidenced by MTT assay. Both ligands were more potent against cancer cells than in skin normal CCD-1064Sk fibroblasts.

Aim:

The aim of this study is to identify the molecular target and mechanism of action of our promising anticancer bis-triazoles MS44-53, focusing specifically on the G-quadruplex stabilizers MS47 and MS49.

Background:

In molecular biology, G-quadruplexes (also known as G4-DNA), one of the higher-order structures of polynucleotides, are four stranded structures formed by nucleic acid sequences which are rich in guanine. They are formed mainly at the single-stranded G-overhang of telomeric DNA and within promoter sequences of genes involved in cellular proliferation and oncogenes such as c-myc, c-kit, and Hsp90. Stabilization of DNA G-quadruplexes is one of the anticancer strategies that has the potential to treat all cancers regardless of the type. A new series of bis-triazoles MS44-53 were developed to stabilize G-quadruplex structures selectively, as G4 ligands and experimental antitumour agents. FRET assay showed that MS47 and MS49 were only the best binders towards the Hsp90 promoter G-guadruplexes. While all bis-triazoles MS44-53 exhibited potent cell growth inhibitory activity against human carcinoma cell lines, suggesting that the ligands perturb molecular targets and mechanisms of action, other than stabilizing G-quadruplexes, contributing to antitumor activity. Therefore, the molecular targets and mechanisms of action of bis-triazoles MS44-53 in different types of human cancer cell lines should be determined by performing further computational studies to MS44-53 and in vitro evaluations for the G-quadruplex stabilizers MS47 and MS49.

Objectives:

1- Determining the exact IC50 for bis-triazoles MS47 & MS49 against four different types of human cancer cell lines; melanoma MDA-MB-435, pancreatic cancer MIA PaCa-2, and colon cancer HCT-116 and COLO 205 cell lines.

2- Predicting the biological targets that bis-triazoles MS44-53 may interact with to trigger or block their biological response.

Methods:

1- MTT assay was used for in vitro evaluation of the antiproliferative activities of MS47 and MS49, and determination of IC50 values.

2- Reverse pharmacophore mapping (pharmacophore profiling) was used for predicting the biological targets of bis-triazoles MS44-53, and determining the % binding probabilities.

Results:

MS49 exhibited more potent proliferation inhibitory activity than MS47 and higher IC50 value against skin normal fibroblasts. Pharmacophore profiling demonstrated FGFR1, PDGFR2, FLT3, mTOR, PPAR-gamma, MUR-F and CETP as biological targets for bis-triazoles MS44-53.

Conclusion:

Bis-triazoles MS47 and MS49 are promising selective innovative compounds with wide spectrum cytotoxic activities against distinct cancer types. Bis-triazoles MS44-53 can be considered as potential drugs to treat different types of carcinoma, in addition to diabetes type II, bacterial infection and cardiovascular diseases.

Other:

Further in vitro evaluations will be performed for bis-triazoles MS44-53 in order to identify their molecular targets and mechanisms of action in different types of human cancer cell lines.

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.

Structure-Based Discovery and Bioactivity Evaluation of Novel Aurora-A Kinase Inhibitors as Anticancer Agents via Docking-Based Comparative Intermolecular Contacts Analysis (dbCICA)

Aurora-A kinase plays a central role in mitosis, where aberrant activation contributes to cancer by promoting cell cycle progression, genomic instability, epithelial-mesenchymal transition, and cancer stemness. Aurora-A kinase inhibitors have shown encouraging results in clinical trials but have not gained Food and Drug Administration (FDA) approval. An innovative computational workflow named Docking-based Comparative Intermolecular Contacts Analysis (dbCICA) was applied—aiming to identify novel Aurora-A kinase inhibitors—using seventy-nine reported Aurora-A kinase inhibitors to specify the best possible docking settings needed to fit into the active-site binding pocket of Aurora-A kinase crystal structure, in a process that only potent ligands contact critical binding-site spots, distinct from those occupied by less-active ligands. Optimal dbCICA models were transformed into two corresponding pharmacophores. The optimal one, in capturing active hits and discarding inactive ones, validated by receiver operating characteristic analysis, was used as a virtual in-silico search query for screening new molecules from the National Cancer Institute database. A fluorescence resonance energy transfer (FRET)-based assay was used to assess the activity of captured molecules and five promising Aurora-A kinase inhibitors were identified. The activity was next validated using a cell culture anti-proliferative assay (MTT) and revealed a most potent lead 85(NCI 14040) molecule after 72 h of incubation, scoring IC₅₀ values of 3.5–11.0 μM against PANC1 (pancreas), PC-3 (prostate), T-47D and MDA-MB-231 (breast)cancer cells, and showing favorable safety profiles (27.5 μM IC₅₀ on fibroblasts). Our results provide new clues for further development of Aurora-A kinase inhibitors as anticancer molecules.

Pharmacophore modeling of JAK1: A target infested with activity-cliffs

Janus kinase 1 (JAK1) is protein kinase involved in autoimmune diseases (AIDs). JAK1 inhibitors have shown promising results in treating AIDs. JAK1 inhibitors are known to exhibit regions of SAR discontinuity or activity cliffs (ACs). ACs represent fundamental challenge to successful QSAR/pharmacophore modeling because QSAR modeling rely on the basic premise that activity is a smooth continuous function of structure. We propose that ACs exist because active ACs members exhibit subtle, albeit critical, enthalpic features absent from their inactive twins. In this context we compared the performances of two computational modeling workflows in extracting valid pharmacophores from 151 diverse JAK1 inhibitors that include ACs: QSAR-guided pharmacophore selection versus docking-based comparative intermolecular contacts analysis (db-CICA). The two methods were judged based on the receiver operating characteristic (ROC) curves of their corresponding pharmacophore models and their abilities to distinguish active members among established JAK1 ACs. db-CICA modeling significantly outperformed ligand-based pharmacophore modeling. The resulting optimal db-CICA pharmacophore was used as virtual search query to scan the National Cancer Institute (NCI) database for novel JAK1 inhibitory leads. The most active hit showed IC₅₀ of 1.04 μM. This study proposes the use of db-CICA modeling as means to extract valid pharmacophores from SAR data infested with ACs.

Deciphering role of FGFR signalling pathway in pancreatic cancer

Recently, fibroblast growth factors are identified to play a vital role in the development and progression of human pancreatic cancer. FGF pathway is critical involved in numerous cellular processes through regulation of its downstream targets, including proliferation, apoptosis, migration, invasion, angiogenesis and metastasis. In this review article, we describe recent advances of FGFR signalling pathway in pancreatic carcinogenesis and progression. Moreover, we highlight the available chemical inhibitors of FGFR pathway for potential treatment of pancreatic cancer. Furthermore, we discuss whether targeting FGFR pathway is a novel therapeutic strategy for pancreatic cancer clinical management.

Recent developments and advances of FGFR as a potential target in cancer

FGFs and their receptors (FGFRs) are critical for many biologic processes, including angiogenesis, wound healing and tissue regeneration. Aberrations in FGFR signaling are common in cancer, making FGFRs a promising target in antitumor studies. To date, many FGFR inhibitors are being detected in clinical studies, and resistance to some inhibitors has emerged. Understanding the mechanisms of resistance is a fundamental step for further implementation of targeted therapies. In this review, we will describe the basic knowledge regarding FGF/FGFR signaling and categorize the clinical FGFR inhibitors. The mechanisms of resistance to FGFR inhibitors and corresponding strategies of overcoming drug resistance will also be discussed.

Ligand-based computational modelling of platelet-derived growth factor beta receptor leading to new angiogenesis inhibitory leads

Platelet derived growth factor beta receptor (PDGFR- β) plays an important role in angiogenesis. PDGFR-β expression is correlated with increased vascularity and maturation of blood vessels in cancer. Pharmacophore modeling and quantitative structure-activity relationship (QSAR) analysis were combined to explore the structural requirements for ligand-PDGFR-β recognition using 107 known PDGFR-β inhibitors. Genetic function algorithm (GFA) coupled to k nearest neighbor (kNN) and multiple linear regression (MLR) analysis were employed to generate predictive QSAR models based on optimal combinations of pharmacophores and physicochemical descriptors. The successful pharmacophores were complemented with exclusion spheres to optimize their receiver operating characteristic curve (ROC) profiles. The QSAR models and their associated pharmacophore hypotheses were validated by identification and experimental evaluation of new angiogenesis inhibitory leads retrieved from the National Cancer Institute (NCI) structural database. Two hits illustrated low micromolar IC₅₀ values in two distinct anti-angiogenesis bioassays.

Computational modeling of the bat HKU4 coronavirus 3CLpro inhibitors as a tool for the development of antivirals against the emerging Middle East respiratory syndrome (MERS) coronavirus

The Middle East respiratory syndrome coronavirus (MERS‐CoV) is an emerging virus that poses a major challenge to clinical management.

The 3C‐like protease (3CL^pro) is essential for viral replication and thus represents a potential target for antiviral drug development. Presently, very few data are available on MERS‐CoV 3CL^pro inhibition by small molecules. We conducted extensive exploration of the pharmacophoric space of a recently identified set of peptidomimetic inhibitors of the bat HKU4‐CoV 3CL^pro. HKU4‐CoV 3CL^pro shares high sequence identity (81%) with the MERS‐CoV enzyme and thus represents a potential surrogate model for anti‐MERS drug discovery. We used 2 well‐established methods: Quantitative structure‐activity relationship (QSAR)‐guided modeling and docking‐based comparative intermolecular contacts analysis. The established pharmacophore models highlight structural features needed for ligand recognition and revealed important binding‐pocket regions involved in 3CL^pro‐ligand interactions. The best models were used as 3D queries to screen the National Cancer Institute database for novel nonpeptidomimetic 3CL^pro inhibitors. The identified hits were tested for HKU4‐CoV and MERS‐CoV 3CL^pro inhibition. Two hits, which share the phenylsulfonamide fragment, showed moderate inhibitory activity against the MERS‐CoV 3CL^pro and represent a potential starting point for the development of novel anti‐MERS agents. To the best of our knowledge, this is the first pharmacophore modeling study supported by in vitro validation on the MERS‐CoV 3CL^pro.

Highlights

MERS‐CoV is an emerging virus that is closely related to the bat HKU4‐CoV.

3CL^pro is a potential drug target for coronavirus infection.

HKU4‐CoV 3CL^pro is a useful surrogate model for the identification of MERS‐CoV 3CL^pro enzyme inhibitors.

dbCICA is a very robust modeling method for hit identification.

The phenylsulfonamide scaffold represents a potential starting point for MERS coronavirus 3CL^pro inhibitors development.

Mechanistic Insights into the Binding of Class IIa HDAC Inhibitors toward Spinocerebellar Ataxia Type-2: A 3D-QSAR and Pharmacophore Modeling Approach

Spinocerebellar ataxia (SCA-2) type-2 is a rare neurological disorder among the nine polyglutamine disorders, mainly caused by polyQ (CAG) trinucleotide repeats expansion within gene coding ataxin-2 protein. The expanded trinucleotide repeats within the ataxin-2 protein sequesters transcriptional cofactors i.e., CREB-binding protein (CBP), Ataxin-2 binding protein 1 (A2BP1) leading to a state of hypo-acetylation and transcriptional repression. Histone de-acetylases inhibitors (HDACi) have been reported to restore transcriptional balance through inhibition of class IIa HDAC's, that leads to an increased acetylation and transcription as demonstrated through in-vivo studies on mouse models of Huntington's. In this study, 61 di-aryl cyclo-propanehydroxamic acid derivatives were used for developing three dimensional (3D) QSAR and pharmacophore models. These models were then employed for screening and selection of anti-ataxia compounds. The chosen QSAR model was observed to be statistically robust with correlation coefficient (r²) value of 0.6774, cross-validated correlation coefficient (q²) of 0.6157 and co-relation coefficient for external test set (pred_r²) of 0.7570. A high F-test value of 77.7093 signified the robustness of the model. Two potential drug leads ZINC 00608101 (SEI) and ZINC 00329110 (ACI) were selected after a coalesce procedure of pharmacophore based screening using the pharmacophore model ADDRR.20 and structural analysis using molecular docking and dynamics simulations. The pharmacophore and the 3D-QSAR model generated were further validated for their screening and prediction ability using the enrichment factor (EF), goodness of hit (GH), and receiver operating characteristics (ROC) curve analysis. The compounds SEI and ACI exhibited a docking score of −10.097 and −9.182 kcal/mol, respectively. An evaluation of binding conformation of ligand-bound protein complexes was performed with MD simulations for a time period of 30 ns along with free energy binding calculations using the g_mmpbsa technique. Prediction of inhibitory activities of the two lead compounds SEI (7.53) and ACI (6.84) using the 3D-QSAR model reaffirmed their inhibitory characteristics as potential anti-ataxia compounds.