Innovative computer-aided methods for the discovery of new kinase ligands

Ligand-based pharmacophore modeling and machine learning for the discovery of potent aurora A kinase inhibitory leads of novel chemotypes

Aurora-A (AURKA) is serine/threonine protein kinase involved in the regulation of numerous processes of cell division. Numerous studies have demonstrated strong association between AURKA and cancer. AURKA is overexpressed in many cancers, such as colon, breast and prostate cancers. Consequently, AURKA has emerged as promising target for therapeutic intervention in cancer management. Herein, we describe a computational workflow for the discovery of novel anti-AURKA inhibitory leads starting with ligand-based assessment of the pharmacophoric space of six diverse sets of inhibitors. Subsequently, machine learning/QSAR modeling was coupled with genetic function algorithm to search for the best possible combination of machine learner, ligand-based pharmacophore(s) and molecular descriptors capable of explaining variation in anti-AURKA bioactivities within a collected list of inhibitors. Two learners succeeded in achieving acceptable structure/activity correlations, namely, random forests and extreme gradient boosting (XGBoost). Three pharmacophores emerged in the successful ML models. These were then used as 3D search queries to mine the National Cancer Institute database for novel anti-AURKA leads. Top-ranking 38 hits were assessed in vitro for their anti-AURKA bioactivities. Among them, three compounds exhibited promising dose-response curves, demonstrating experimental IC50 values ranging from sub-micromolar to low micromolar values. Remarkably, two of these compounds are of novel chemotypes.

Structure‐Based Approaches for the Prediction of Alzheimer's Disease Inhibitors: Comparative Interactions Analysis, Pharmacophore Modeling and Molecular Dynamics Simulations

Due to its significant role in neurodegeneration, Cyclin‐dependent kinase 5 (CDK5) has emerged as a potential target for addressing neuropathological disorders, including Alzheimer's disease (AD). The application of CDK5 inhibitors has demonstrated promise in the treatment of AD. This prompted us to model this interesting target using a computational workflow named Docking‐based Comparative Intermolecular Contacts Analysis (dbCICA). Approaches including 3D‐QSAR, genetic algorithm, and pharmacophore modeling were employed to discover new CDK inhibitors.

Discovery of new PKN2 inhibitory chemotypes via QSAR-guided selection of docking-based pharmacophores

Serine/threonine-protein kinase N2 (PKN2) plays an important role in cell cycle progression, cell migration, cell adhesion and transcription activation signaling processes. In cancer, however, it plays important roles in tumor cell migration, invasion and apoptosis. PKN2 inhibitors have been shown to be promising in treating cancer. This prompted us to model this interesting target using our QSAR-guided selection of docking-based pharmacophores approach where numerous pharmacophores are extracted from docked ligand poses and allowed to compete within the context of QSAR. The optimal pharmacophore was sterically-refined, validated by receiver operating characteristic (ROC) curve analysis and used as virtual search query to screen the National Cancer Institute (NCI) database for new promising anti-PKN2 leads of novel chemotypes. Three low micromolar hits were identified with IC₅₀ values ranging between 9.9 and 18.6 µM. Pharmacological assays showed promising cytotoxic properties for active hits in MTT and wound healing assays against MCF-7 and PANC-1 cancer cells.

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.

“MedChemVR”: A Virtual Reality Game to Enhance Medicinal Chemistry Education

Medicinal chemistry (MC) is an indispensable component of the pharmacy curriculum. The pharmacists’ unique knowledge of a medicine’s chemistry enhances their understanding of the pharmacological activity, manufacturing, storage, use, supply, and handling of drugs. However, chemistry is a challenging subject for both teaching and learning. These challenges are typically caused by the inability of students to construct a mental image of the three-dimensional (3D) structure of a drug molecule from its two-dimensional presentations. This study explores a prototype virtual reality (VR) gamification option, as an educational tool developed to aid the learning process and to improve the delivery of the MC subject to students. The developed system is evaluated by a cohort of 41 students. The analysis of the results was encouraging and provided invaluable feedback for the future development of the proposed system.

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.

Discovery of New Phosphoinositide 3-kinase Delta (PI3Kδ) Inhibitors via Virtual Screening using Crystallography-derived Pharmacophore Modelling and QSAR Analysis

BACKGROUND

PI3Kδ is predominantly expressed in hematopoietic cells and participates in the activation of leukocytes. PI3Kδ inhibition is a promising approach for treating inflammatory diseases and leukocyte malignancies. Accordingly, we decided to model PI3Kδ binding.

METHODS

Seventeen PI3Kδ crystallographic complexes were used to extract 94 pharmacophore models. QSAR modelling was subsequently used to select the superior pharmacophore(s) that best explain bioactivity variation within a list of 79 diverse inhibitors (i.e., upon combination with other physicochemical descriptors).

RESULTS

The best QSAR model (r2 = 0.71, r2 LOO = 0.70, r2 press against external testing list of 15 compounds = 0.80) included a single crystallographic pharmacophore of optimal explanatory qualities. The resulting pharmacophore and QSAR model were used to screen the National Cancer Institute (NCI) database for new PI3Kδ inhibitors. Two hits showed low micromolar IC50 values.

CONCLUSION

Crystallography-based pharmacophores were successfully combined with QSAR analysis for the identification of novel PI3Kδ inhibitors.

Combining Stochastic Deformation/Relaxation and Intermolecular Contacts Analysis for Extracting Pharmacophores from Ligand-Receptor Complexes

We previously combined molecular dynamics (classical or simulated annealing) with ligand-receptor contacts analysis as a means to extract valid pharmacophore model(s) from single ligand-receptor complexes. However, molecular dynamics methods are computationally expensive and time-consuming. Here we describe a novel method for extracting valid pharmacophore model(s) from a single crystallographic structure within a reasonable time scale. The new method is based on ligand-receptor contacts analysis following energy relaxation of a predetermined set of randomly deformed complexes generated from the targeted crystallographic structure. Ligand-receptor contacts maintained across many deformed/relaxed structures are assumed to be critical and used to guide pharmacophore development. This methodology was implemented to develop valid pharmacophore models for PI3K-γ, RENIN, and JAK1. The resulting pharmacophore models were validated by receiver operating characteristic (ROC) analysis against inhibitors extracted from the CHEMBL database. Additionally, we implemented pharmacophores extracted from PI3K-γ to search for new inhibitors from the National Cancer Institute list of compounds. The process culminated in new PI3K-γ/mTOR inhibitory leads of low micromolar IC₅₀s.

Ligand-based modeling of Akt3 lead to potent dual Akt1/Akt3 inhibitor

Akt1 and Akt3 are important serine/threonine-specific protein kinases involved in G2 phase required by cancer cells to maintain cell cycle and to prevent cell death. Accordingly, inhibitors of these kinases should have potent anti-cancer properties. This prompted us to use pharmacophore/QSAR modeling to identify optimal binding models and physicochemical descriptors that explain bioactivity variation within a set of 74 diverse Akt3 inhibitors. Two successful orthogonal pharmacophores were identified and further validated using receiver operating characteristic (ROC) curve analyses. The pharmacophoric models and associated QSAR equation were applied to screen the national cancer institute (NCI) list of compounds for new Akt3 inhibitors. Six hits showed significant experimental anti-Akt3 IC₅₀ values, out of which one compound exhibited dual low micromolar anti-Akt1 and anti-Akt3 inhibitory profiles.

Simulated annealing molecular dynamics and ligand-receptor contacts analysis for pharmacophore modeling

AIM

Ligand-based pharmacophore modeling requires long list of inhibitors, while pharmacophores based on single ligand-receptor crystallographic structure can be too restricted or promiscuous.

METHODOLOGY

This prompted us to combine simulated annealing molecular dynamics (SAMD) with ligand-receptor contacts analysis as means to construct pharmacophore model(s) from single ligand-receptor complex. Ligand-receptor contacts that survive numerous heating-cooling SAMD cycles are considered critical and are used to guide pharmacophore development.

RESULTS

This methodology was implemented to develop pharmacophores for acetylcholinesterase and protein kinase C-θ. The resulting models were validated by receiver-operating characteristic analysis and in vitro bioassay. Assay identified four new protein kinase C-θ inhibitors among captured hits, two of which exhibited nanomolar potencies.

CONCLUSION

The results illustrate the ability of the new method to extract valid pharmacophores from single ligand-protein complex.

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

AIM

FGFR-1 is an oncogenic kinase involved in several cancers. FGFR1-specific inhibitors have shown promising results against several human cancers prompting us to model this interesting target. Toward the end, we implemented elaborate ligand-based and structure-based computational workflows to explore the pharmacophoric requirements for potent FGFR-1 inhibitors. Results & methodology: Structure-based and ligand-based modeling applied on 59 diverse FGFR-1 inhibitors yielded novel pharmacophore and quantitative structure-activity relationship models that were used to scan the National Cancer Institute's structural database for novel leads. Four potent hits were captured, with the most active having IC₅₀ of 426 nM. Identities and purities of active hits were established using nuclear magnetic resonance and mass spectroscopy.

CONCLUSION

Elaborate ligand-based (pharmacophore/quantitaive structure-activity relationship) and structure-based (docking-based comparative intermolecular contacts analysis) modeling provided deep understanding of ligand binding within FGFR-1 as evidenced by the virtually captured new potent leads.