A real-world study and network pharmacology analysis of EGFR-TKIs combined with ZLJT to delay drug resistance in advanced lung adenocarcinoma
Authors: Han Xue , Liang Lan , He Chenming , Ren Qinyou , Su Jialin , Cao Liang , Zheng Jin 
Published: November 2023

AbstractObjectiveThis study aimed to explore the efficacy and safety of combining epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) with ZiLongJin Tablet (ZLJT) in delaying acquired resistance in advanced EGFR-mutant lung adenocarcinoma (LUAD) patients. Furthermore, we employed network pharmacology and molecular docking techniques to investigate the underlying mechanisms.MethodsA retrospective comparative study was conducted on stage IIIc/IV LUAD patients treated with EGFR-TKIs alone or in combination with ZLJT at the Second Affiliated Hospital of the Air Force Medical University between January 1, 2017, and May 1, 2023. The study evaluated the onset of TKI resistance, adverse reaction rates, safety indicators (such as aspartate aminotransferase, alanine aminotransferase, and creatinine), and inflammatory markers (neutrophil-to-lymphocyte ratio and platelet-to-lymphocyte ratio) to investigate the impact of EGFR-TKI combined with ZLJT on acquired resistance and prognostic indicators. Additionally, we utilized the Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform, the Bioinformatics Analysis Tool for Molecular Mechanism of Traditional Chinese Medicine, PubChem, UniProt, and Swiss Target Prediction databases to identify the active ingredients and targets of ZLJT. We obtained differentially expressed genes related to EGFR-TKI sensitivity and resistance from the Gene Expression Omnibus database using the GSE34228 dataset, which included sensitive (n = 26) and resistant (n = 26) PC9 cell lines. The "limma" package in R software was employed to detect DEGs. Based on this, we constructed a protein‒protein interaction network, performed gene ontology and KEGG enrichment analyses, and conducted pathway network analysis to elucidate the correlation between the active ingredients in ZLJT and signaling pathways. Finally, molecular docking was performed using AutoDockVina, PYMOL 2.2.0, and Discovery Studio Client v19.1.0 software to simulate spatial and energy matching during the recognition process between predicted targets and their corresponding compounds.Results(1) A total of 89 patients were included, with 40 patients in the EGFR-TKI combined with ZLJT group (combination group) and 49 patients in the EGFR-TKI alone group (monotherapy group). The baseline characteristics of the two groups were comparable. There was a significant difference in the onset of resistance between the combination group and the monotherapy group (P < 0.01). Compared to the monotherapy group, the combination group showed a prolongation of 3.27 months in delayed acquired resistance. There was also a statistically significant difference in the onset of resistance to first-generation TKIs between the two groups (P < 0.05). (2) In terms of safety analysis, the incidence of adverse reactions related to EGFR-TKIs was 12.5% in the combination group and 14.3% in the monotherapy group, but this difference was not statistically significant (P > 0.05). There were no statistically significant differences in serum AST, ALT, CREA, TBIL, ALB and BUN levels between the two groups after medication (P > 0.05). (3) Regarding inflammatory markers, there were no statistically significant differences in the changes in neutrophil-to-lymphocyte Ratio(NLR) and Platelet-to-lymphocyte Ratio(PLR) values before and after treatment between the two groups (P > 0.05). (4) Network pharmacology analysis identified 112 active ingredients and 290 target genes for ZLJT. From the GEO database, 2035 differentially expressed genes related to resistant LUAD were selected, and 39 target genes were obtained by taking the intersection. A "ZLJT-compound-target-disease" network was successfully constructed using Cytoscape 3.7.0. GO enrichment analysis revealed that ZLJT mainly affected biological processes such as adenylate cyclase-modulating G protein-coupled receptor. In terms of cellular components, ZLJT was associated with the cell projection membrane. The molecular function primarily focused on protein heterodimerization activity. KEGG enrichment analysis indicated that ZLJT exerted its antitumor and anti-drug resistance effects through pathways such as the PI3K-Akt pathway. Molecular docking showed that luteolin had good binding activity with FOS (-9.8 kJ/mol), as did tanshinone IIA with FOS (-9.8 kJ/mol) and quercetin with FOS (-8.7 kJ/mol).ConclusionZLJT has potential antitumor progression effects. For patients with EGFR gene-mutated non-small cell LUAD, combining ZLJT with EGFR-TKI treatment can delay the occurrence of acquired resistance. The underlying mechanisms may involve altering signal transduction pathways, blocking the tumor cell cycle, inhibiting tumor activity, enhancing cellular vitality, and improving the bioavailability of combination therapy. The combination of EGFR-TKI and ZLJT represents an effective approach for the treatment of tumors using both Chinese and Western medicine.