Discovery of Novel AURKA Inhibitors for Triple-Negative Breast Cancer (TNBC) Therapy via a Hybrid Virtual Screening Pipeline, Biological Evaluation and Molecular Dynamics Simulation
A cascaded AI-driven virtual screening pipeline is developed, integrating sequence-based affinity prediction, equivariant deep learning docking (KarmaDock), and geometric rescoring (DeepDock) to identify novel AURKA inhibitor candidates.
Abstract
Aurora kinase A (AURKA) is a pivotal driver of malignant progression and poor prognosis in triple-negative breast cancer (TNBC). In this study, we developed a cascaded AI-driven virtual screening pipeline, integrating sequence-based affinity prediction (PSICHIC), equivariant deep learning docking (KarmaDock), and geometric rescoring (DeepDock) to identify novel AURKA inhibitor candidates. From an in-house 160,000-compound screening library assembled from commercially available collections, three leads (compounds 3, 5, and 8) were selected and subsequently validated via HTRF biochemical assays, exhibiting potent enzymatic inhibition with IC50 values of 157 nM, 21.64 nM, and 46.03 nM, respectively. Cell-based assays demonstrated that compound 3 produced stronger short-term cell-growth inhibition in MDA-MB-231 (TNBC) cells compared to clinical benchmarks MLN8237 and CCT241736, whereas compounds 3 and 5 showed cell-growth inhibition in NIH/3T3 cells within the same concentration range as the reference inhibitors. Triplicate 500 ns molecular dynamics simulations supported stable binding modes of the identified leads in the AURKA binding pocket. Additional computational analyses further provided supportive information for subsequent lead optimization. This study provides a transparent and open-source workflow for AI-assisted identification of AURKA-active chemotypes.
Neuroblastoma (NB) is the most common extracranial malignant solid tumor in children, and aberrant activation of anaplastic lymphoma kinase (ALK)—driven primarily by hotspot mutations F1174L and R1275Q—represents a core oncogenic driver of high-risk and relapsed NB. Clinical ALK inhibitors are frequently hampered by acquired drug resistance and off-target toxicity, highlighting an unmet need for novel selective ALK inhibitors for NB treatment. Here, we employed a drug repurposing strategy integrating virtual screening, multi-dimensional molecular dynamics (MD) simulations, steered MD, and umbrella sampling to screen 7,322 compounds for novel ALK inhibitory activity, followed by systematic
in vitro
and
in vivo
validation. We identified CDK12-IN-3 (C18) as a high-affinity ligand for wild-type ALK and both NB-associated ALK mutants (F1174L, R1275Q), with favorable binding selectivity over 10 homologous tyrosine kinases.
In vitro
, CDK12-IN-3 exerted nanomolar-level antiproliferative, anti-clonogenic, and anti-migratory effects across three genetically distinct NB cell lines (SH-SY5Y, SK-N-BE (2), IMR-32), while showing no significant cytotoxicity in ALK-negative HEK293T normal cells even at 20 μM, demonstrating a wide
in vitro
safety window. Mechanistic studies confirmed that CDK12-IN-3 directly bound ALK (KD = 5.33 μM by surface plasmon resonance), downregulated ALK expression, suppressed the downstream PI3K/AKT/mTOR signaling pathway, and induced NB cell apoptosis.
In vivo
, oral administration of CDK12-IN-3 (50 mg/kg) significantly suppressed NB xenograft tumor growth in nude mice with no observable systemic toxicity. Collectively, our findings validate CDK12-IN-3 as a novel selective ALK inhibitor with potent anti-NB activity, providing a promising lead compound and an efficient dynamic simulation-guided drug repurposing paradigm for overcoming ALK inhibitor resistance in high-risk NB.
Tianyi Liu, Xuejiao Hu, Wenxin Yan et al.· Frontiers in Chemistry· 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
The natural product compounds CNP0456830 and CNP0467494 exhibited the lowest binding free energies for both EGFR and PIK3CA, identifying them as the most promising dual-target inhibitors.
Si-miao Lu, Yi Zhu, Yong-tao Han et al.· Diseases of the esophagus· 0 citations
This LBVS-SBVS-ADMET-MD pipeline effectively identified three promising TNIK inhibitors, providing a solid foundation for future experimental validation and potential development of targeted therapies for Wnt-driven malignancies.
D. Mishra, Rajnish Kumar, Anurag T. K. Baidya et al.· Talanta: The International J...· 0 citations
In prostate cancer (PCa), the TMPRSS2-ERG fusion gene drives aberrant activation of androgen receptor (AR) signaling, while TMPRSS2 protease activity contributes to remodeling the tumor microenvironment. Targeting TMPRSS2 can suppress cancer metastasis and represents a potential therapeutic strategy for castration-resistant prostate cancer (CRPC). This study aimed to screen and identify promising TMPRSS2-binding lead candidates via integrated virtual screening and experimental validation. Screening of the ChemDiv, ChemBridge, and TargetMol compound libraries yielded 53 hit compounds. Surface plasmon resonance (SPR) assays were subsequently performed to evaluate binding affinities. Six candidate compounds (15, 26, 29, 31, 41, and 42) exhibited higher affinity for TMPRSS2 compared to the reference inhibitor Nafamostat (KD = 6.09 × 10-5 M). Among them, compound 26 showed the highest affinity (KD = 3.88 × 10-6 M). Subsequent in vitro inhibition assays using LNCaP clone FGC cells demonstrated that compound 26 possessed the strongest anti-proliferative activity (IC50 = 10.95 μM), indicating its therapeutic potential against PCa. Confocal microscopy showed apparent reduced TMPRSS2 staining after 24 h compound 26 treatment, implying potential TMPRSS2 modulation. Further mechanistic insights were gained through molecular dynamics (MD) simulations, alanine scanning mutagenesis, quantum mechanics/molecular mechanics (QM/MM) calculations, and dynamical cross-correlation matrix (DCCM) analysis, which revealed specific ligand-protein interactions in the compound 26-TMPRSS2 complex and reinforced the reliability of our findings. Collectively, these results identify compound 26 as a promising TMPRSS2-binding lead candidate, providing a solid foundation for further biochemical verification and structural optimization.
Huiru Xie, Qibo Hu, Qiuhong Zhang et al.· Journal of Molecular Graphic...· 0 citations