Unraveling tissue-specific molecular targets of dihydroartemisinin in non-small cell lung cancer: an integrative machine learning and network pharmacology approach
摘要
Non-small cell lung cancer (NSCLC) presents significant therapeutic challenges due to resistance and immune evasion. Dihydroartemisinin (DHA), a derivative of artemisinin, exhibits broad anti-tumor activity, but its molecular targets and mechanisms in NSCLC remain unclear. To identify the core therapeutic targets and elucidate the mechanism of action of DHA against NSCLC using an integrated computational and bioinformatics approach. Potential DHA targets were predicted using PharmMapper, SEA, SwissTargetPrediction, SuperPred, and TargetNet. NSCLC-associated targets were retrieved from OMIM, GeneCards, and CTD. Transcriptomic datasets (GSE101929, GSE118370, GSE116959, GSE159857) were integrated and analyzed for differential expression (limma) and co-expression networks (WGCNA). KEGG pathway enrichment identified key pathways. Protein–protein interaction networks, machine learning (Lasso regression, Random Forest), nomogram construction, immune infiltration analysis (ssGSEA), miRNA-mRNA network analysis (miRTarBase), and molecular docking (CB-Dock2) were performed to identify and validate core targets. We identified 1277 potential DHA targets and 44 consensus NSCLC targets. Integration of DEGs (1240 genes) and WGCNA modules (3 key modules, 2860 genes) yielded 1128 overlapping genes. KEGG enrichment revealed 15 key pathways. Machine learning on 196 pathway-enriched DHA targets identified 12 candidate genes. Validation confirmed 6 core targets: AR, CASP3, CDK1, CDK4, PTK2, MMP9. A nomogram based on the 12 targets showed excellent predictive power (AUC = 0.987). Immune profiling revealed significant alterations in 21 immune cell types in NSCLC, and correlation analysis linked core targets (e.g., CDK1/CDK4 with T cell subsets, MMP9 with myeloid cells) to immune dysregulation. Molecular docking confirmed strong binding affinities between DHA and all 6 core targets, with CDK1 exhibiting the highest affinity (− 8.8 kcal/mol). miRNA networks identified key regulators like hsa-miR-15b-5p and hsa-miR-302a-3p. This study delineates AR, CASP3, CDK1, CDK4, PTK2, and MMP9 as core therapeutic targets of DHA in NSCLC. DHA exerts its anti-NSCLC effects through direct inhibition of these targets (particularly high-affinity binding to CDK1) and modulation of the tumor immune microenvironment, including T-cell memory, cytotoxic function, myeloid-mediated remodeling, and immunosuppressive cell subsets. These findings provide a mechanistic foundation for developing DHA as a therapeutic agent or adjuvant for NSCLC, especially in combination with immunotherapy.
Graphical abstract