Literature retrieval results
A total of 7714 articles were retrieved, including 2585 articles on CNKI, 3689 articles on Wanfang, 1142 articles on VIP, 60 articles on pubmed, and 238 articles on web of science. After removing the duplicate literature, the remaining 5310 articles were collected, 298 prescriptions were collected, finally, 269 prescriptions were remained after deleting the repeated prescriptions.
Medication rule analysis
A total of 295 traditional Chinese medicines were used from 269 prescriptions, and among them, the five most frequently used were HQ, FL, BZ, Ophiopogon japonicus (MD, in Chinese), Glycyrrhiza uralensis (GC, in Chinese) (Fig. 2A); the three drug pairs of BZ-HQ, Hedyotis diffusa(BHSSC, in Chinese)-HQ and BZ-FL were ranked the highest in support, and they were also the most frequently used in the analysis of association rules (Table1). After clustering analysis of TCMs, it was found that Codonopsis pilosula(DS, in Chinese), HQ, BZ and FL were clustered into the same category, which indicated that these TCMs have similar properties and functions (Fig. 2B). Due to the highest frequency and correlation of HQ, BZ and FL based on the above analysis results, they were selected as the core drugs for subsequent analysis.
Table 1 Correlated herb pairs (top 10)
couplet drugs
|
support
|
confidence
|
degree
|
BZ-HQ
|
0.29
|
0.74
|
79
|
BHSSC-HQ
|
0.29
|
0.81
|
78
|
BZ-FL
|
0.29
|
0.72
|
77
|
FL-HQ
|
0.29
|
0.66
|
77
|
MD-HQ
|
0.28
|
0.68
|
74
|
DS-HQ
|
0.27
|
0.71
|
72
|
GC-HQ
|
0.26
|
0.66
|
71
|
DS-BZ
|
0.25
|
0.65
|
66
|
DS-FL
|
0.23
|
0.62
|
63
|
CP-FL
|
0.18
|
0.72
|
47
|
* CP (Tangerine Peel, in Chinese)
Network Pharmacology Analysis
Core targets and component screening
A total of 40 active ingredients of core drugs that met the screening conditions were collected from TCMSP, and 9 were supplemented through literature review (Xu et al. 2018; Chen et al. 2019; Zhang et al. 2019; Dou et al. 2021; Xu et al. 2021; Acharya et al. 2022; Mao et al. 2022; Zhang et al. 2022; Jiang et al. 2023). The collected components were imported into the SwissTargetPrediction database, then 156 potential targets that met the screening conditions were obtained. A total of 1309,1490,3926 and 168 NSCLC targets were collected from GeneCards, OMIM, DisGeNET and TTD databases, respectively. After removing duplicate targets, 5549 NSCLC targets were obtained. After intersecting the potential targets of active ingredients and NSCLC targets, 116 intersection targets were obtained, that is, the potential targets of core TCM in the treatment of NSCLC (Fig. 2C). According to the intersection targets obtained above, three key components were screened by sorting the degree values of the components corresponding to these targets: Quercetin, Kaempferol and Isorhamnetin.
PPI network construction
Table 2 Key targets of herbal medicine against NSCLC (top 10)
Target
|
Betweenness
|
Closeness
|
Degree
|
AKT1
|
1955.66
|
0.01
|
65.00
|
EGFR
|
1547.19
|
0.01
|
60.00
|
ESR1
|
1281.93
|
0.01
|
52.00
|
PTGS2
|
1144.71
|
0.01
|
46.00
|
SRC
|
792.72
|
0.01
|
46.00
|
MMP9
|
356.01
|
0.01
|
41.00
|
GSK3B
|
393.78
|
0.01
|
37.00
|
PARP1
|
498.49
|
0.01
|
38.00
|
MCL1
|
257.82
|
0.01
|
32.00
|
KDR
|
214.27
|
0.01
|
32.00
|
The 116 intersection targets were imported into the STRING database to obtain the PPI network. The obtained PPI network was analyzed by Cytoscape and the screening condition was greater than the mean value, that is, degree > 14.85, closeness centrality > 0.0042, betweenness centrality > 129.01. The key targets of TCMs components in the treatment of NSCLC were screened (Table 2). Because the degree, betweenness centrality and closeness centrality of EGFR and AKT1 were significantly higher than those of other targets, they were identified as the core targets of Astragalus, Atractylodes and Poria in the treatment of NSCLC (Fig. 2D). Then, they were molecularly docked with the three core components (Quercetin, Kaempferol and Isorhamnetin) screened above.
Enrichment analysis
A total of 110 KEGG signaling pathways were enriched, including PI3K-Akt signaling pathway, EGFR tyrosine kinase inhibitor resistance, FoxO signaling pathway and other signaling pathways closely related to malignant tumors (Fig. 2E). 1982 GO entries was enriched by GO functional enrichment, involving 1703 biological processes, 90 cellular components, and 189 molecular functions (Fig. 2F). According to the results of KEGG enrichment, both the significance and the number of enriched intersection targets are higher in the PI3K-Akt signaling pathway, and existing studies have also shown that it is closely related to tumors (He et al. 2021; Khezri et al. 2022; et al. 2022). Therefore, we will futher verify it by Western blot experiments. In addition, we plotted the component-target-PI3K/Akt signaling pathway according to the corresponding relationship between PI3K-Akt pathway, targets and components (Fig. 2G).
Molecular docking
Table 3 Docking binding energy
Components
|
AKT1
kcal/mol
|
EGFR
kcal/mol
|
Quercetin
|
-6.53
|
-6.64
|
Isorhamnetin
|
-6.58
|
-6.79
|
Kaempferol
|
-6.25
|
-6.50
|
The results of molecular docking showed that the binding energies of the two core targets and the three key components were all less than -1.2 kcal/mol, which indicated that they could be spontaneously and stably combined, and that Quercetin, Kaempferol and Isorhamnetin could act on AKT1 and EGFR targets to play an anti-NSCLC role. Among them, the binding energy of EGFR and Isorhamnetin was the lowest (-6.79 kcal/mol) (Fig. 2H, Table 3). From the diagram, it could be found that the main forces between isorhamnetin and EGFR are Carbon Hydrogen Bond, Pi-Sigma, Pi-Sulfur and Pi-Alkyl. Isorhamnetin can form Pi-Alkyl bond with Lysa745, Alaa743 and Vala726 of EGFR, Pi-Sigma bond with Meta790, Leua718 and Leua844, Pi-Sulfur bond with Meta766, and Carbon Hydrogen Bond with meta793, Glna791, Glua792, Lysa745 and Leua718.
Cell experiment
Determined by cell viability assay
In this experiment, the effects of Quercetin, Kaempferol and Isorhamnetin on the viability of A549 cells at 6 different concentrations and 3 time points were investigated. The results showed that the viability of A549 cells and the three components showed the dependence between dose and time. The difference between each treatment group and the blank group was statistically significant (P < 0.05), and the inhibitory effect of Isorhamnetin was the strongest. Therefore, the subsequent experiments focused on the effect of Isorhamnetin on A549 cells (Fig. 3A-D). The uncontrolled proliferation of cells is a major feature of tumors. When the tumor proliferates to a certain extent, it will transfer from the original lesion to other parts by means of tumor blood vessels. Therefore, the inhibition of tumor cell proliferation can effectively control tumor progression. It can be seen from the result that the viability of A549 cells was significantly inhibited after treatment with 100 μmol/L Isorhamnetin for 48 h.
Cell morphology observation
From the inverted microscope and confocal microscope, it could be found that the cell morphology of the blank group was full and irregular paving stone, while the cell morphology of the treatment group became elongated and slender. In the laser confocal microscope image, it could also be observed that the PI dye could not penetrate the membrane due to the integrity of the plasma membrane, , while the Dio dye could penetrate the membrane due to its lipophilicity in the normal group, so the cells in the normal group only have green fluorescence in the cytoplasmic region, while there is no red fluorescence in the nuclear region; however, for the treatment group, due to the influence of isorhamnetin, the cells began to apoptosis and the cell membrane was damaged. Therefore, most cells not only had green fluorescence in the cytoplasmic region, but also had red fluorescence in the nuclear region (Fig. 3E-H). The above results suggest that isorhamnetin will shrink the cytoplasm of A549 cells and recover the pseudopods, thus their morphology will change. At the same time, it will also promote the apoptosis of the cells and lead to the damage of the cell membrane.
Cell migration
The results of scratch test showed that the cell migration rates were 39.30%, 32.40%, 21.02% and 15.17% after 48 h treatment with Isorhamnetin at concentrations of 0, 90, 180 and 270 μmol / L, respectively. Although there was no significant difference between the 90 μmol/L group and the control group, the inhibition of cell migration ability became stronger with the increase of drug concentration. When 180 and 270 μmol/L were compared with the control group, the difference was statistically significant (P < 0.05), indicating that the higher the concentration of Isorhamnetin, the more significant the inhibition of A549 cell migration ability (Fig. 4A-H, M). NSCLC is easy to metastasize, which is closely related to poor prognosis. It can be seen from this experiment that although the low concentration group has no obvious inhibitory effect on the migration of A549 cells, with the increase of administration dose, the viability and migration ability of A549 cells are gradually inhibited. Its migration ability is concentration-dependent, indicating that isorhamnetin can effectively inhibit the migration of A549 cells.
Cell invasion
The results of this experiment showed that with the increase of dosage, the fewer cells penetrated the membrane. The number of A549 cells in the 0, 90, 180 and 270 μmol/L groups was (451 ± 23.52), (210 ± 17.52), (119.67 ± 4.93) and (71.67± 13.65), respectively. The difference between the administration group and the control group was statistically significant (P < 0.05) (Fig. 4I-L, N). Therefore, Isorhamnetin can effectively inhibit the invasion of A549 cells, and the higher the concentration, the stronger the inhibitory effect on its invasion ability.
Apoptosis and cell cycle
The effect of Isorhamnetin on the apoptosis of A549 cells was detected by FITC and PI staining. The results showed that the apoptosis rates of 0,90,180 and 270 μmol/L groups were 7.29%, 15.04%, 23.87% and 30.74 %, respectively. The early apoptosis rate of 90 μmol/L treatment group was 11.94%, while the late apoptosis rate was 3.10%. The early apoptosis rate of 180 μmol/L treatment group was 18.67%, and the late apoptosis rate was 5.20%. The early apoptosis rate of 270 μmol/L treatment group was 20.51%, and the late apoptosis rate was 10.23%. It can be seen that with the increase of drug concentration, the number of apoptotic cells also increased, and the number of late apoptotic cells also increased synchronously. It shows that with the increase of Isorhamnetin concentration, its effect on A549 cells gradually changes from proliferation inhibition to direct killing, and the higher the concentration, the more obvious the effect (Fig. 5A-D). This is also consistent with the results of cell morphology observation, indicating that isorhamnetin can induce apoptosis of A549 cells.
PI staining was used to detect the effect of isorhamnetin on the cell cycle of A549 cells. The results showed that the proportion of G1 phase cells in the blank group was 63.17%, while the proportion of G1 phase cells in the 90, 180 and 270 μmol/L groups was 77.73%, 83.49% and 90.92%, respectively. It can be seen that with the increase of drug concentration, the G1 phase arrest is more obvious (Fig. 5E-H). This result further indicated thatA549 cells were arrested in G1 phase, and the process of DNA synthesis and cell division became slower with the increase of drug concentration. It is suggested that isorhamnetin can inhibit DNA replication and promote cell apoptosis, thereby reduce the rate of cell proliferation and controlling tumor diseases.
Western blot
The results of Western blot showed that there was no significant difference in the expression of AKT and PI3K between the blank group and each treatment group (P > 0.05), while the expression of P-AKT and P-PI3K decreased with the increase of drug concentration. Compared with the blank group, there was no statistical significance in the 90 μmol/L group,but there was statistical significance in the 180 and 270 μmol/L groups (P < 0.05) (Fig. 5I-K). From the previous theoretical analysis, it can be seen that the PI3K-Akt signaling pathway has a regulatory effect on a series of biological processes such as proliferation, apoptosis and migration of tumor cells, and the two targets of PI3K and AKT are the core targets of the pathway. It can be seen from the results that the expression of P-PI3K and P-Akt in each treatment group decreased. Although there was no significant difference in the expression between the 90 μmol/L group and the blank group, but the expression decreased with the increase of drug concentration, showing a dose-dependent manner. It was further proved that isorhamnetin can regulate the PI3K-Akt signaling pathway by inhibiting AKT1 and PI3K phosphorylation, thereby inhibit cell growth and proliferation and promoting apoptosis.