Findings highlight Z4P as a promising mutation-resilient IRE1 inhibitor and validate the effectiveness of the integrated computational pipeline for identifying potential anti-cancer therapeutics.
Abstract
Breast cancer is one of the most prevalent and lethal malignancies affecting women globally. The increasing resistance to current therapeutic strategies highlights the need for novel molecular targets. Inositol-requiring enzyme 1 alpha (IRE1α), a key sensor in the unfolded protein response (UPR), has emerged as a promising therapeutic target due to its role in tumour progression and survival. This study employed an integrative in silico approach combining machine learning, molecular docking, and molecular dynamics simulations to identify potent, non-toxic IRE1α inhibitors for breast cancer treatment. An initial library of 115 compounds retrieved from ChEMBL and MedChemExpress was used for machine learning-based toxicity modelling. Literature curation identified 44 reported IRE1α inhibitors, which were reduced to 38 unique compounds following duplicate removal. Drug-likeness and ADMET screening using SwissADME and ProTox retained 22 compounds for further evaluation. Molecular docking was performed using AutoDock, followed by Dynamics simulations in GROMACS to assess stability. Machine Learning (ML) models were developed for both toxicity regression and binary toxicity classification analyses. Toxicity prediction models were developed using twenty physicochemical and pharmacokinetic descriptors. In the regression analysis, Random Forest demonstrated the strongest cross-validation performance (R² = 0.5998 ± 0.3439), while the stacking ensemble achieved the highest test-set performance (R² = 0.9765), although differences among ensemble methods were not statistically significant. In the complementary classification analysis, the Support Vector Machine (SVM) achieved the highest discriminative performance with an ROC-AUC value of 0.98. Docking studies revealed that Z4P exhibited the strongest binding affinity (- 7.93 kcal/mol) to the wild-type IRE1, compared with the control drug MKC8866 (- 6.7 kcal/mol). Additionally, Z4P exhibited a higher binding energy of - 9.5 kcal/mol, whereas MKC8866 had a binding energy of - 6.94 kcal/mol. MD simulations over 200 ns confirmed the stability of the IRE1-Z4P complex, with favourable RMSD, RMSF, Rg, and SASA profiles relative to the control. These findings highlight Z4P as a promising mutation-resilient IRE1 inhibitor and validate the effectiveness of the integrated computational pipeline for identifying potential anti-cancer therapeutics.
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
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
Tyrosyl-DNA phosphodiesterase I (TDP1) repairs topoisomerase I (TOP1)–mediated DNA damage and is a promising anticancer target, particularly in combination with TOP1 inhibitors. However, the discovery of potent and drug-like TDP1 inhibitors remains challenging due to the limited structural diversity of known active compounds. Here, we developed an integrated computational framework combining machine learning (ML), deep learning (DL), and structure-based docking with experimental validation. A curated dataset of 2040 compounds (857 active, 1183 inactive) was assembled and analyzed by scaffold composition. A total of 40 binary classification models were constructed using six ML algorithms and a deep neural network (DNN), each paired with five molecular fingerprint representations, along with five graph neural network architectures (GCN, GAT, MPNN, AttentiveFP, and FPGNN). The SVM::RDKitDes model performed best (AUC = 0.89, F1 = 0.78, BA = 0.80), with robustness confirmed by Y-scrambling and randomized-split analyses, and SHAP analysis identified 20 key descriptors of TDP1 inhibition. The model was deployed as a web application (http://drugpred.top:5050) and standalone desktop applications (.exe) are available at https://github.com/zenghuang8006/TDP1-inhibitor-prediction. The validated model was applied to screen 201 231 compounds, followed by drug-likeness filtering and hierarchical docking, yielding 16 candidates. Biological evaluation identified compound AO65 as a potent TDP1 inhibitor (IC50 = 0.80 ± 0.02 µM), and quantum chemical calculations and docking elucidated its electronic properties and binding within the catalytic domain. This work demonstrates the value of integrating ML-driven prediction with structure-based approaches and identifies AO65 as a promising lead for further TDP1-focused investigation.
Huang Zeng, Manyi Zhang, Bo Qiu et al.· RSC Advances· 0 citations
Background: HER2 is a key oncogenic gene in breast cancer, involved in tumor progression, metastasis, and therapeutic resistance. This study aimed to find new HER2 inhibitors using a hybrid of machine learning (ML) and structure-based virtual screening (VS), combined with molecular dynamics (MD) simulations on various scaffolds. Methods: Four supervised molecular fingerprint classification models were trained on a dataset of 10,000 validated compounds from ChEMBL. Random Forest was the top model for screening a large compound library. Selected compounds underwent molecular docking in the HER2 ATP binding site, ADMET, drug likeness, toxicity analysis, and 200 ns MD simulations. Methods like PCA, FEL, hydrogen-bond analysis, DCCM, RDF, salt-bridge analysis, and MM/PBSA were used to assess binding stability. Results: Virtual screening identified three compounds, CHMEBL193865 (Lead-1), CHMEBL46740 (Lead-2), and CHMEBL151318 (Lead-3)—with better binding affinity and interaction profiles than the reference inhibitor. MD simulations showed stable protein–ligand complexes with RMSD values of 2.32–2.76 Å. Among these, Lead-2 was the most structurally stable, and Lead-1 had the most favorable binding free energy. All three compounds showed good drug likeness, ADMET properties, and low predicted toxicity. Conclusions: These findings support further in vitro and in vivo testing for developing new therapeutics against HER2-overexpressing breast cancer, highlighting two scaffolds with promising lead optimization potential.
Alhumaidi B. Alabbas, Safar M. Alqahtani· Pharmaceuticals· 0 citations