This work demonstrates that integrating docking-derived pharmacophores with conformational ensemble-based machine learning provides an effective approach for discovering novel inhibitors against underexplored kinase targets.
Abstract
Interferon-inducible RNA-dependent protein kinase (PKR) is an emerging therapeutic target involved in cancer, neurodegeneration, and inflammatory disorders; however, the discovery of potent and structurally diverse PKR inhibitors remains limited. In this study, we report an integrated computational–experimental strategy for the identification of novel PKR inhibitory chemotypes. Pharmacophore models were generated from flexible docking of structurally diverse reference inhibitors and validated using receiver operating characteristic (ROC) analysis. To better capture ligand flexibility and address data scarcity, multiple conformations of 223 PKR inhibitors were employed as a data augmentation strategy in machine learning-based QSAR modelling. Among several algorithms evaluated, a Naïve Bayes classifier combined with genetic function algorithm (GFA) feature selection provided the most predictive model. The optimized model, incorporating a single pharmacophore hypothesis and key physicochemical descriptors, was applied to virtual screening of the Boehringer Ingelheim opnMe compound library. Experimental validation using a PKR kinase assay identified the pre-synthesized compound BI-8128 as a potent PKR inhibitor, with an IC50 value of 174.8 nM, demonstrating higher potency than the reference inhibitor C16 under identical conditions. Notably, this compound represents a structurally distinct scaffold compared to reported PKR inhibitors, indicating effective scaffold hopping. To the best of our knowledge, this is the first report of PKR inhibitory activity for BI-8128. Overall, this work demonstrates that integrating docking-derived pharmacophores with conformational ensemble-based machine learning provides an effective approach for discovering novel inhibitors against underexplored kinase targets.
The proto-oncogene serine/threonine kinase PIM2 is a critical regulator of cell proliferation, survival, and tumor progression and represents an attractive therapeutic target for several cancers. In this study, an integrated machine learning–guided computational pipeline was developed to identify potential PIM2 inhibitors by combining quantitative structure–activity relationship (QSAR) modeling, virtual screening, molecular docking, molecular dynamics (MD) simulations, and pharmacokinetic prediction. Bioactivity data for PIM2 inhibitors were retrieved from the ChEMBL database, yielding 5953 compounds. After data cleaning, structural standardization, and removal of duplicates and invalid entries, a curated dataset of 1584 compounds was obtained for QSAR modeling. To address dataset imbalance, the Synthetic Minority Oversampling Technique (SMOTE) was applied before model development. Twelve molecular fingerprint descriptors were generated and used to construct 180 QSAR models using five machine learning algorithms, including Random Forest (RF), Extreme Gradient Boosting (XGBoost), Support Vector Regression (SVR), k-Nearest Neighbors (KNN), and Multilayer Perceptron (MLP). Among these models, the Random Forest–fingerprint model demonstrated the best predictive performance, achieving a mean R2 of 0.971 with low prediction errors (RMSE = 0.271; MAE = 0.125) across training, testing, and cross-validation datasets. The optimized model was subsequently applied to virtual screening of multiple chemical libraries, including FDA-approved drugs, natural product databases, and commercial compound collections. Several promising candidates were identified, including TCMBANKIN000009 (emetine), Amb28533044 (4,6′-Anhydrooxysporidinone), NPC170963 (Lysophosphatidylcholine (15:0)), NPC262768 (Endosulfan), and NPC469603 (8-hydroxyircinialactam A). Molecular docking showed that these compounds bind within the ATP-binding pocket of PIM2 kinase, forming interactions with key residues such as Lys62, Asp125, Asp128, and Glu168. Subsequent molecular dynamics simulations confirmed the stability of selected complexes, demonstrating reduced residue fluctuations, stable protein compactness, and persistent intermolecular interactions during the simulation. Furthermore, ADMET prediction suggested favorable pharmacokinetic and toxicity profiles for several compounds. Collectively, these findings highlight the potential of the identified molecules as promising PIM2 inhibitor candidates, providing valuable leads for future experimental validation and anticancer drug development.
A. Fahira, M. Shahab, Zaheer Ud Din et al.· Journal of Genetic Engineeri...· 0 citations
Results demonstrate that approaches incorporating machine learning and structure-based computational analysis can be useful for discovering and prioritizing potential IRAK4 inhibitor candidates.
H. Na, Juwon Park, Jiwon Choi· Current Issues in Molecular...· 0 citations
The findings from molecular docking, 500-ns molecular dynamics simulations, MM-GBSA calculations, and alanine scanning analyses collectively corroborate a stable binding mode of BTB11556 within the MPO active site, support further investigation of BTB11556 as a candidate compound associated with MPO-targeted therapeutic strategies.
Maysoon Raed Alnajdawi, H. Wahab, Belal Alnajjar et al.· Journal of Computer-Aided Mo...· 0 citations
The proposed workflow efficiently reduced a large chemical space to a focused set of TNKS1 inhibitor candidates while substantially reducing the experimental screening burden, highlighting the value of integrating consensus ML, SBVS, and experimental validation to accelerate early-stage hit discovery for TNKS1 and other therapeutic targets.
M. Bilotta, Adriana Gargano, R. Rocca et al.· Pharmaceuticals· 0 citations
An integrated computational workflow combining explainable machine learning, virtual screening, molecular dynamics simulations, and binding free-energy calculations to identify novel inhibitors of this drug-resistant EGFR variant may support the development of new therapeutic strategies for overcoming resistance in EGFR-driven cancers.
Jurica Novak· International Journal of Mol...· 0 citations