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

Computational workflow for discovering small molecular binders for shallow binding sites by integrating molecular dynamics simulation, pharmacophore modeling, and machine learning: STAT3 as case study

STAT3 belongs to a family of seven transcription factors. It plays an important role in activating the transcription of various genes involved in a variety of cellular processes. High levels of STAT3 are detected in several types of cancer. Hence, STAT3 inhibition is considered a promising therapeutic anti-cancer strategy. However, since STAT3 inhibitors bind to the shallow SH2 domain of the protein, it is expected that hydration water molecules play significant role in ligand-binding complicating the discovery of potent binders. To remedy this issue, we herein propose to extract pharmacophores from molecular dynamics (MD) frames of a potent co-crystallized ligand complexed within STAT3 SH2 domain. Subsequently, we employ genetic function algorithm coupled with machine learning (GFA-ML) to explore the optimal combination of MD-derived pharmacophores that can account for the variations in bioactivity among a list of inhibitors. To enhance the dataset, the training and testing lists were augmented nearly a 100-fold by considering multiple conformers of the ligands. A single significant pharmacophore emerged after 188 ns of MD simulation to represent STAT3-ligand binding. Screening the National Cancer Institute (NCI) database with this model identified one low micromolar inhibitor most likely binds to the SH2 domain of STAT3 and inhibits this pathway.

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.

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.

AKT Inhibitors: The Road Ahead to Computational Modeling-Guided Discovery

AKT, is a serine/threonine protein kinase comprising three isoforms-namely: AKT1, AKT2 and AKT3, whose inhibitors have been recognized as promising therapeutic targets for various human disorders, especially cancer. In this work, we report a systematic evaluation of multi-target Quantitative Structure-Activity Relationship (mt-QSAR) models to probe AKT' inhibitory activity, based on different feature selection algorithms and machine learning tools. The best predictive linear and non-linear mt-QSAR models were found by the genetic algorithm-based linear discriminant analysis (GA-LDA) and gradient boosting (Xgboost) techniques, respectively, using a dataset containing 5523 inhibitors of the AKT isoforms assayed under various experimental conditions. The linear model highlighted the key structural attributes responsible for higher inhibitory activity whereas the non-linear model displayed an overall accuracy higher than 90%. Both these predictive models, generated through internal and external validation methods, were then used for screening the Asinex kinase inhibitor library to identify the most potential virtual hits as pan-AKT inhibitors. The virtual hits identified were then filtered by stepwise analyses based on reverse pharmacophore-mapping based prediction. Finally, results of molecular dynamics simulations were used to estimate the theoretical binding affinity of the selected virtual hits towards the three isoforms of enzyme AKT. Our computational findings thus provide important guidelines to facilitate the discovery of novel AKT inhibitors.

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.

miR-664b-5p Inhibits Hepatocellular Cancer Cell Proliferation Through Targeting Oncogene AKT2

Background: miR-664b-5p accelerates the development of certain cancers, but the role of miR-664b-5p in hepatocellular carcinoma (HCC) has been less reported. Therefore, the authors aimed to study the role of miR-664b-5p in HCC progression. Materials and Methods: miR-664b-5p expression in liver cancer and adjacent tissues, and in HepG2 and SUN-475 cells, was measured by quantitative real-time polymerase chain reaction (qRT-PCR). Relationship between miR-664b-5p and AKT2 was predicted by TargetScan and confirmed by dual-luciferase reporter assay, and gene or protein expressions were determined by performing qRT-PCR and Western blotting. The viability and apoptosis, and the migration and invasion of HepG2 and SUN-475 cells were determined by CCK-8 assay and flow cytometry, and transwell assay, respectively. Results: Downregulated miR-664b-5p was observed in hepatocellular cancer tissues. Functional analyses revealed that miR-664b-5p mimic suppressed viability, migration, and invasion, but promoted apoptosis in HepG2 and SUN-475 cells. AKT2 was a target of miR-664b-5p, whose mimics inhibited the expression of AKT2. However, upregulated AKT2 promoted viability, migration, and invasion, but inhibited apoptosis in HepG2 and SUN-475 cells, and such effects were reversed by miR-664b-5p mimics. Conclusions: miR-664b-5p acts as a cancer suppressor through negatively regulating AKT2 expression in HepG2 and SUN-475 cells, suggesting that miR-664b-5p could be a protective target for HCC patients.

Investigation of Binding Characteristics of Phosphoinositide-dependent Kinase-1 (PDK1) Co-crystallized Ligands Through Virtual Pharmacophore Modeling Leading to Novel Anti-PDK1 Hits

BACKGROUND

3-Phosphoinositide Dependent Protein Kinase-1 (PDK1) is being lately considered as an attractive and forthcoming anticancer target. A Protein Data Bank (PDB) cocrystallized crystal provides not only rigid theoretical data but also a realistic molecular recognition data that can be explored and used to discover new hits.

OBJECTIVE

This incited us to investigate the co-crystallized ligands' contacts inside the PDK1 binding pocket via a structure-based receptor-ligand pharmacophore generation technique in Discovery Studio 4.5 (DS 4.5).

METHODS

Accordingly, 35 crystals for PDK1 were collected and studied. Every single receptorligand interaction was validated and the significant ones were converted into their corresponding pharmacophoric features. The generated pharmacophores were scored by the Receiver Operating Characteristic (ROC) curve analysis.

RESULTS

Consequently, 169 pharmacophores were generated and sorted, 11 pharmacophores acquired good ROC-AUC results of 0.8 and a selectivity value above 8. Pharmacophore 1UU3_2_01 was used in particular as a searching filter to screen NCI database because of its acceptable validity criteria and its distinctive positive ionizable feature. Several low micromolar PDK1 inhibitors were revealed. The most potent hit illustrated anti-PDK1 IC50 values of 200 nM with 70% inhibition against SW480 cell lines.

CONCLUSION

Eventually, the active hits were docked inside the PDK1 binding pocket and the recognition points between the active hits and the receptor were analyzed that led to the discovery of new scaffolds as potential PDK1 inhibitors.