Skip to content
Open access

Identification of novel PKR inhibitory chemotypes via integrated molecular modelling and machine learning

Aug 2026 · RSC Advances · 0 citations · 60 references
Medicine

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.

Read PDF

Similar papers

Open access Jul 2026

Machine learning–driven identification of PIM2 kinase inhibitors through QSAR modeling and molecular dynamics simulations

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. · 0 citations
Open access Aug 2026

Machine Learning-Driven Drug Repositioning Identifies Putative IRAK4 Inhibitors Through Structure-Based Computational Evaluation

Interleukin-1 receptor-associated kinase 4 (IRAK4) is one of the IRAK family proteins and plays an important role in the regulation of innate and inflammatory responses. In particular, IRAK4 acts as a key regulator of the Toll-like receptor (TLR) and interleukin-1 receptor (IL-1R) signaling pathways and has attracted attention as a therapeutic target for immune and inflammatory diseases. In this study, an integrated computational approach combining machine learning, molecular docking, and molecular dynamics simulations was applied to identify putative IRAK4 inhibitor candidates. Bioactivity data of IRAK4 were obtained from the ChEMBL and PubChem databases and evaluated for multiple binary classification models. The optimized XGBoost model based on ECFP4 and PubChem fingerprints achieved an ROC-AUC of 0.996 and an average precision (AP) of 0.991 on the independent test set. After that, 20 candidate compounds with high predictive probability score were finally selected through subsequent screening of the DrugBank database. Among them, DB12168 (MK-0557), DB15040 (TP-271), and DB18152 (Zilurgisertib) exhibited favorable binding free energies and stable complex formation with IRAK4 through molecular dynamics simulations and MM-PBSA calculations. Overall, these 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 · 0 citations
Open access Aug 2026

An Integrated Consensus Machine Learning and Structure-Based Workflow for the Discovery of Novel Tankyrase 1 Inhibitors

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. · 0 citations
Open access Aug 2026

From Descriptor Learning to Binding Stability: An Explainable Machine Learning Pipeline for EGFR Double-Mutant Inhibitor Discovery

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 · 0 citations