Qingfei Jiedu decoction inhibits PD-L1 expression in lung adenocarcinoma based on network pharmacology analysis, molecular docking and experimental verification

Objective: We aim at investigating the molecular mechanisms through which the Qingfei Jiedu decoction (QFJDD) regulates PD-L1 expression in lung adenocarcinoma (LUAD). Methods: Bioactive compounds and targets of QFJDD were screened from TCMSP, BATMAN-TCM, and literature. Then, GeneCard, OMIM, PharmGKB, Therapeutic Target, and DrugBank databases were used to identify LUAD-related genes. The protein-protein interaction (PPI) network was constructed using overlapping targets of bioactive compounds in LUAD with the Cytoscape software and STRING database. The potential functions and pathways in which the hub genes were enriched by GO, KEGG, and DAVID pathway analyses. Molecular docking of bioactive compounds and key genes was executed via AutoDock Vina. Qualitative and quantitative analyses of QFJDD were performed using UPLC-Q-TOF-MS and UPLC. Expressions of key genes were determined by qRT-PCR, immunoreactivity score (IRS) of PD-L1 was assessed by immunohistochemistry (IHC), while the CD8+PD-1+T% derived from spleen tissues of Lewis lung cancer (LLC) bearing-mice was calculated using flow cytometry (FCM). Results: A total of 53 bioactive compounds and 288 targets of QFJDD as well as 8151 LUAD associated genes were obtained. Further, six bioactive compounds, including quercetin, luteolin, kaempferol, wogonin, baicalein, and acacetin, and 22 hub genes were identified. The GO analysis showed that the hub genes were mainly enriched in DNA or RNA transcription. KEGG and DAVID pathway analyses revealed that 20 hub genes were primarily enriched in virus, cancer, immune, endocrine, and cardiovascular pathways. The EGFR, JUN, RELA, HIF1A, NFKBIA, AKT1, MAPK1, and MAPK14 hub genes were identified as key genes in PD-L1 expression and PD-1 checkpoint pathway. Moreover, ideal affinity and regions were identified between core compounds and key genes. Notably, QFJDD downregulated EGFR, JUN, RELA, HIF1A, NFKBIA, and CD274 expressions (p < 0.05), while it upregulated AKT1 and MAPK1 (p < 0.05) levels in A549 cells. The PD-L1 IRS of LLC tissue in the QFJDD high dose (Hd) group was lower than model group (p < 0.01). CD8+PD-1+T% was higher in the QFJDD Hd group than in normal and model groups (p < 0.05). Conclusion: QFJDD downregulates PD-L1 expression and increases CD8+PD-1+T% via regulating HIF-1, EGFR, JUN and NFκB signaling pathways. Therefore, QFJDD is a potential treatment option for LUAD.


Introduction
According to the global cancer data 2020, lung cancer, especially small cell lung cancer and non-small cell lung cancer (NSCLC), is the most lethal cancer (Sung et al., 2021). In China, among all the histological subtypes of NSCLC, lung adenocarcinoma (LUAD) is the most prevalent. Lung cancer is often diagnosed in the late stages and therefore has a poor overall survival rate (Duma et al., 2019). Anti-tumor drugs have been shown to improve the prognosis of lung cancer . Moreover, the advent of programmed cell death-1 (PD-1)/programmed cell death ligand-1 (PD-L1) has resulted in better prognostic outcomes of NSCLC patients (Jain et al., 2018). However, due to the high heterogeneity of NSCLC and the varying degrees of T cell infiltrations into the tumor microenvironment, the prognostic outcomes for patients remains poor. Therefore, there is a need to identify new therapeutic agents for NSCLC.
For thousands of years, Traditional Chinese medicine (TCM) has been used as a complementary and alternative medicine (Cyranoski, 2018) to treat various diseases, including malignancies . In recent years, TCM are combined with chemotherapy, targeted therapies, or immune checkpoint inhibitors to treat various cancers . However, the mechanisms of action of TCM prescriptions have not been fully established. Over the years, network pharmacology have been used to study the pharmacological mechanisms of famous prescriptions recorded in ancient TCM books, but, the mechanisms of empirical prescriptions have rarely been explored. Qingfei Jiedu Decoction (QFJDD) is an empirical prescription that has been used for the complementary intervention of lung cancer by our medical team (Pan et al., 2020a). This prescription is prepared from Scutellariae Barbatae Herba (Banzhilian, BZL), Lobeliae Chinensis Herba (Banbianlian, BBL), Hedyotis Diffusae Herba (Baihuasheshecao, BHSSC), Herba Solani Lyrati (Baimaoteng, BMT), Solanum Nigrum (Longkui, LK), and Coicis Semen (Yiyiren, YYR) ( Table 1). Previously, we showed that QFJDD inhibits Lewis lung cancer (LLC) cell proliferations and downregulates PD-L1 expression in tumor tissues (Pan et al., 2020b). However, the active ingredients and pharmacological mechanisms of QFJDD have not been established.
In this study, network pharmacology and experimental verification were performed to investigate the mechanisms through which QFJDD inhibits PD-L1 expression in LUAD ( Figure 1). First, bioactive compounds and targets of QFJDD against LUAD were searched and screened. Then, hub genes, biological functions, and the key signaling pathways associated with QFJDD were identified by network construction and analysis. Finally, the results were validated through cellular and animal experiments.
2 Materials and methods 2.1 Data preparation and collation 2.1.1 Bioactive compounds and action targets of Qingfei Jiedu decoction At first, the Traditional Chinese Medicine Systems Pharmacology database and Analysis Platform (TCMSP) (Ru et al., 2014), Bioinformatics Analysis Tool for Molecular  (Liu et al., 2016a), and the published literature were used to screen bioactive compounds in BHSSC, BZL, BBL, BMT, YYR, and LK. Bioactive compounds from BATMAN-TCM database were numbered successively by BATMAN001, BATMAN002, etc. Then, the SwissADME web tool (Daina et al., 2017) was used to detect the oral bioavailability (OB), drug-likeness (DL) and pharmacokinetics of each bioactive compounds. The screening criteria were OB ≥ 30% and DL ≥ 0.18 (Tu et al., 2021). The screening criterion of core bioactive compounds was degree >24 . Moreover, the potential targets of bioactive compounds in QFJDD were searched from TCMSP, BATMAN-TCM, and SwissTargetPrediction (Daina et al., 2019). Finally, the duplicate bioactive compounds and targets were removed to obtain the herb-gene text file.

Network construction and analysis 2.2.1 Herb-compound-target network
The herb-gene text file was processed using Strawberry Perl programming language to obtain files of the all-net-list, col-herblist, col-gene-list, col-id-list, id-pie-input-list and item-typesinput-list. Afterward, the Cytoscape 3.9.1 software (Shannon et al., 2003) was used to establish the herb-compound-target (H-C-T) network.

Construction of the protein-protein interaction and hub genes network
A total of 244 herb-LUAD genes were submitted to the Cytoscape 3.9.1 software to identify hub genes. Then, the STRING 11.5 database (https://www.string-db.org/) (Szklarczyk et al., 2021) was used to construct the PPI network with a confidence score >0.9. Targets were limited to the Homo sapiens species.
nih.gov/) website (Kim et al., 2021). Subsequently, Uniprot IDs of seven protein receptor conformations corresponding to key genes were retrieved from the Protein Data Bank (PDB, https://www.rcsb. org/) database (Velankar et al., 2021). Screening parameters were: 1) X-crystal diffraction and crystal resolution<3Å; 2) Protein structures of Homo sapiens. Based on PyMOL version 2.5 software, the following options were carried out: assignation of bond orders, addition of hydrogens, creation of zero-order bonds to metals, creation of disulfide bonds, deletion of waters beyond 5Å from het groups and calculation of the molecular activity pockets. Molecular docking simulation was executed using AutoDock Vina version 1.1.2 version software (Seeliger and de Groot, 2010). Binding energy was determined from affinity.

Qualitative and quantitative analyses of six bioactive compounds in Qingfei Jiedu decoction
Ultra-Performance Liquid Chromatography Quadrupole Timeof-Flight Mass Spectrometry (UPLC-Q-TOF-MS) was used to isolate and identify six bioactive compounds of QFJDD which were extracted and characterized using UPLC-Q-TOF-MS on the Waters ® SYNAPT ® G2-Si system with a waters CORTECS UPLC T3 Column (2.1 mm × 100 mm, 1.67 μm). Solvent A (acetonitrile) and solvent B (0.1% formic acid-H 2 O) were used for gradient elution as shown in Supplementary Table S1. The MS, which was run in both positive and negative electrospray ionization, parameters are shown in Supplementary Table S2. Quercetin, luteolin, kaempferol, wogonin, baicalein, and acacetin were purchased from Chengdu herbpurify Co., Ltd, and the purity of all standards has met the analytical requirements. The MS instrumentation and data acquisition were conducted using the Masslynx ® V4.1 software.
Moreover, a volume of 2.8 ml QFJDD was accurately taken into a volumetric flask and fixed the volume to 25 ml with ultrapure water. After that, the supernatant was obtained by centrifugation at 1847 (×g) for 5 min and filtered via a 0.45 µm filter membrane. Quercetin, luteolin, kaempferol, wogonin, baicalein and acacetin solutions were prepared in methanol at final concentrations of 0.1088 mg/ml, 0.0856 mg/ ml, 0.1072 mg/ml, 0.1440 mg/ml, 0.1108 mg/ml, and 0.1164 mg/ ml, separately. A further 100 μl of each standard solution was taken into a volumetric flask and diluted with 900 μl methanol. UPLC was used to quantify six bioactive compounds of QFJDD which were analyzed using UPLC on the Waters ACQUITY UPLC ® H-Class system with a ACQUITY UPLC BEH C18 Column (2.1 mm × 100 mm, 1.7 μm) (Supplementary Figure S1). Solvent A (0.2% phosphoric acid-distilled water) and solvent B (methanol) were used for gradient elution as shown in Supplementary Table S3. Standard curves were constructed as shown in Supplementary Table S4 and the amount of six bioactive compounds in QFJDD was calculated. The setting parameters were as follows: column temperature = 40°C; injection volume = 5 μl; wavelength = 350 nm.

Cell and animal cultures
The A549 and LLC cell lines were purchased from the Stem Cell Bank, Chinese Academy of Sciences (Shanghai, China). The A549 cells were cultured in F-12K (Invitrogen) with 10% fetal bovine serum (FBS, Gibco) and 1% Glutamax (Invitrogen). Incubation was done at 37°C in a 5% CO 2 atmosphere. The LLC cells were cultured in DMEM (Gibco) with 10% FBS and incubated at 37°C in a 5% CO 2 atmosphere. Six-week-old male Sprague-Dawley (SD) rats (n = 20, 200 ± 10 g) and five-week-old male C57BL/6 mice (n = 25, 14-16 g) of SPF grade were purchased from the Shanghai SLAC Laboratory Animal Co., Ltd. [license No. SCXK (Shanghai) 2017-0005]. Animals were maintained at the Zhejiang Chinese Medical University Laboratory Animal Center (22 ± 2°C and 40%-70% relative humidity with a 12-h light/12-h dark cycle) with ad libitum access to food and water.

Preparation of Qingfei Jiedu decoctioncontaining serum
After 5 days of adaptive feeding, SD rats were randomized into two groups (n = 10). Equivalent doses between humans and rats were converted according to body surface area (Nair and Jacob, 2016). Thus, the control group was administered with 0.9% saline (1 ml/100 g) while the QFJDD group received 8.64 g/kg of QFJDD. Rats were intragastrically administered with drugs once daily for seven consecutive days ). An hour after the final drug administration, rats were anesthetized by 3% sodium pentobarbital via peritoneal injection after which whole blood was collected by cardiac puncture. Blood samples were allowed to stand at room temperature for 2 h, and afterwards the serum was prepared by centrifugation at 1847 (×g) for 10 min at 4°C. After that, the complement was inactivated at 56°C for 60 min. Serum was passed through a 0.22 μm filter to remove bacteria and stored at −80°C.
Frontiers in Pharmacology frontiersin.org 04 2.4.5 Establishment of Lewis lung cancerbearing mouse model and drug intervention After 7 days of adaptive feeding, C57BL/6 mice were randomized into six groups, each containing five mice. After that, 0.2 ml of 1×10 7 / ml LLC cell suspension was subcutaneously inoculated into the right axilla of C57BL/6 mice, apart from the normal group. When the tumor diameter reached 0.5 cm, drugs were initiated and administered for 14 days. The high dose (H d , 57.6 g/kg/d), medium dose (M d , 28.8 g/kg/d) and low dose (L d , 14.4 g/kg/d) groups were administered with once-daily doses of QFJDD via oral gavage. The PD-L1 inhibitor group (BMS-202, MedChemExpress, United States, 17.2 mg/kg) received onceweekly injections via the caudal vein. In addition, model and normal groups received a once-daily oral gavage of 0.9% sterile saline. At 24 h after the last dose, tumor and the spleen tissues were aseptically harvested. The equivalent dose between human and mouse was determined based on body surface area (Nair and Jacob, 2016). The doses of QFJDD in different groups were based on our previous study (Pan et al., 2020b).

Flow cytometry
Connective tissues of the spleen were aseptically removed from mice and washed in precooled PBS. The spleen was punctured repeatedly with a sterile needle after which PBS was injected to Frontiers in Pharmacology frontiersin.org drain the lymphocytes. After that, tissues were filtered using a 200mesh screen, washed twice using PBS, and centrifuged at 1,200 rpm for 5 min to obtain a single cell suspension of 1×10 6 /ml. Subsequently, 100 μl of the single cell suspension was added into a flow tube. 2 μl of CD3 (eBioscience, United States), 2 μl of CD8 (eBioscience, United States), and 2 μl PD-1 (eBioscience, United States) labeled antibody were added to the experimental tube, while the corresponding same type of control antibody was added to the control tube, and incubated at room temperature in dark for 15 min. 500 μl red blood cell lysis buffer (Sigma Aldrich, United States) was then added and incubated at room temperature for 10 min. 5 ml precooled PBS was added to cells and centrifuged at 1,200 rpm for 5 min, after which the supernatant was discarded. After centrifugation and re-suspension, the percentage of CD8 + PD-1 + T was determined by flow cytometry (Agilent Novocyte, United States).

Statistical analysis
Data are presented as mean ± SD. Statistical analyses were performed using SPSS 25.0 software. Expressions of key genes in A549 cells were compared by one-way analysis of variance (ANOVA) (Bewick et al., 2004). Differences in PD-L1 IRS in subcutaneous tumor tissues between multiple groups were determined by Kruskal-Wallis nonparametric analysis of variance with multiple group comparisons (Theodorsson-Norheim, 1986). Comparison of CD8 + PD-1 + T% among multiple groups were performed by Kruskal-Wallis nonparametric analysis of variance. p-value<0.05 was considered statistically significant.

Bioactive compounds and candidate targets network of Qingfei Jiedu decoction
Based on the filter criteria of OB ≥ 30% and DL ≥ 0.18, preliminary screening generated 4 compounds in BHSSC, 27 compounds in BZL, 18 compounds in BBL, 1 compound in BMT, 6 compounds in YYR and 7 compounds in LK. After removing duplicate values, 53 bioactive compounds and 288 targets were obtained (Figure 2, Supplementary Table S6). Moreover, six core bioactive compounds, including quercetin, luteolin, kaempferol, wogonin, baicalein, and acacetin were identified with the screening criterion of degree>24 (Table 2). Pharmacokinetic parameters are presented in Table 3.

The gene ontology and kyoto encyclopedia of genes and genomes pathway enrichment analyses
To investigate the potential functions of the hub genes, GO enrichment analysis of enriched biological processes (BP), cellular components (CC) and molecular functions (MF) was performed ( Figure 4A). With regards to BPs, hub genes were enriched in regulation of DNA-binding transcription factor activity, positive regulation of pri-miRNA transcription by RNA polymerase II, response to radiation, cellular response to chemical stress, and response to drug among others. The top five enriched CCs were transcription regulator complex, RNA polymerase II transcription regulator complex, spindle, transferase complex, transferring phosphorus-containing groups and vesicle lumen. The significantly enriched MFs were DNA-binding transcription factor binding, RNA polymerase II-specific DNA-binding transcription factor binding, ubiquitin protein ligase binding, ubiquitin-like protein ligase binding, DNA-binding transcription GI, absorption, Gatrointestinal absorption; BBB, permeant, according to the yolk of the BOILED-Egg. All pharmacokinetic parameters were detected using SwissADME web tool.
Frontiers in Pharmacology frontiersin.org 07 activator activity, RNA polymerase II-specific and DNA-binding transcription activator activity.
The KEGG analysis of the top 30 signaling pathways revealed that 20 of the 22 hub genes were primarily involved in viral or bacterial infections (Kaposi sarcoma-associated herpesvirus infection, hepatitis B, human T-cell leukemia virus 1 infection, Epstein-Barr virus infection, human cytomegalovirus infection, hepatitis C, Measles, and salmonella infection), various cancers (prostate cancer, breast cancer, colorectal cancer, chronic myeloid leukemia, small cell lung cancer, gastric cancer, endometrial cancer, pancreatic cancer, bladder cancer, and proteoglycans in cancer), immune system (PD-L1 expression and PD-1 checkpoint pathway in cancer, Th17 cell differentiation, IL-17 signaling pathway, chemical carcinogenesis-receptor activation, MAPK signaling pathway, apoptosis, and TNF signaling pathway), endocrine system (endocrine resistance, relaxin signaling pathway, and prolactin signaling pathway) and cardiovascular system (lipid and atherosclerosis and fluid shear stress and atherosclerosis) ( Figure 4B). Noticeably, the top three genes enriched in the above-mentioned signaling pathways were AKT1, MAPK1, and RELA ( Figure 4C, Supplementary Table S7).

The database for annotation, visualization, and integrated discovery pathway analysis
The PD-1/PD-L1 inhibitors are vital in cancer immunotherapy (Pan et al., 2021;Sezer et al., 2021). Based on the results shown in Figures 4B,C, we explored the regulation of tumor cells by the PD-1/PD-L1 pathway.

FIGURE 3 (A) A Venn diagram of the intersecting LUAD-related genes constructed using DrugBank, GeneCard, OMIM, PharmGKB, and TTD databases. (B)
The overlapping genes between QFJDD and LUAD. (C-E) The process of identifying twenty two hub genes. (F) The PPI network of twenty two hub genes targeted by QFJDD in LUAD. In the PPI diagram, one gene is represented by a circle, and the protein structure is displayed in the center of the circle.

Molecular docking simulation
AutoDock Vina version 1.1.2 software was used to calculate binding active pockets and binding energy between receptor proteins and ligands. Using PyMOL version 2.5, the six core bioactive compounds were observed to enter the active pockets of eight receptors, respectively. The receptor-ligand affinities are shown in Table 5. Three-dimensional structures of several receptorligand binding regions are shown in (Figure 6). Taking the HIF1Aquercetin pair in affinity (kcal/mol) as an example for analysis, quercetin formed hydrogen bonds with residues with SER-274, ASP-249, LEU-248, and THR-288 ( Figure 6C).

Effects of Qingfei Jiedu decoction-Containing serum on PD-L1 expression and PD-1 checkpoint pathway in a A549 cell line
Cellular experiments were conducted to validate the effects of QFJDD-containing serum on key genes of the hsa05235 pathway. mRNA expressions of EGFR, HIF1A and CD274 in QFJDDcontaining serum group and Osimertinib group were significantly reduced relative to control group (p < 0.01; p < 0.001, Figures 8C,E,H). Transcription levels of JUN and RELA (p < 0.01, Figures 8A,B) as well as NFKBIA (p < 0.001, Figure 8D) were significantly decreased in the QFJDD-containing group relative to the control group. On the contrary, expression of JUN and NFKBIA (p < 0.001, Figures 8A,D) as well as RELA (p < 0.01, Figure 8B) were markedly increased in the Osimertinib group relative to QFJDD-containing serum group. Moreover, compared with the control group, the mRNA expression of AKT1 (p < 0.001, Figure 8F) and MAPK1 (p < 0.05, Figure 8G) in the QFJDD- containing serum group and Osimertinib group was significantly increased. In conclusion, QFJDD-containing serum inhibits PD-L1 expression by regulating the above key genes.

Effects of Qingfei Jiedu decoction on PD-L1 expression in tumor tissues and CD8 + PD-1 + % in spleen tissues from Lewis lung cancer-bearing mice
Animal assays to assess the effects of QFJDD on PD-L1 expression levels in LLC tissues (Figures 9A-E) and CD8 + PD-1 + % spleen tissues (Figures 9G-L). The results showed that PD-L1 IRS of the H d group was lower than Model group (p < 0.01). However, differences between the other groups were insignificant ( Figure 9F). Further, analysis of CD8 + PD-1 + T percentage in the spleen tissue of the Normal group, Model group, L d group, M d group, H d group, and BMS-202 group by FCM revealed that CD8 + PD-1 + T% in H d and BMS-202 groups were significantly elevated relative to Normal or Model groups (p < 0.05, p < 0.01, Figure 9M). Therefore, QFJDD reduces the PD-L1 IRS in tumor tissues and facilitates the differentiations of CD8 + T into toxic CD8 + PD-1 + T cells.

Discussion
TCM is safe, feasible and improves the clinical efficacies of chemotherapy, radiotherapy, targeted therapy and immunotherapy (Zhang et al., 2018;Su et al., 2020;Zhang FIGURE 4 (A,B) The GO and KEGG enrichment analyses of hub genes of QFJDD. (C) Network of twenty hub genes and thirty signaling pathways. The pink square represents hub gene, and the green triangle indicates signaling pathway. The square and triangle size reflects node degree. The higher the degree value, the bigger the node size.
Frontiers in Pharmacology frontiersin.org Frontiers in Pharmacology frontiersin.org 11 Pan et al. 10.3389/fphar.2022.897966    Frontiers in Pharmacology frontiersin.org 13   , 2021). Therefore, it has been widely used in adjuvant lung cancer treatment. The PD-1/PD-L1 inhibitors are effective for malignant tumor treatment, and have become a focus of immunotherapy (Steven et al., 2016). Primarily, PD-L1 is expressed on tumor cells. This protein inhibits T-cell activation and proliferation, particularly cytotoxic T lymphocytes (CTLs) by binding to PD-1, which is predominantly expressed on activated T cells (Poggio et al., 2019). In healthy cells, PD-1 and PD-L1 regulate T cell response amplitudes and maintains self-tolerance (Ai et al., 2020). However, cancer cells hijack the PD-1/PD-L1 pathway to evade immune surveillance by overexpressing PD-L1, resulting in cancer cell proliferation and metastasis. QFJDD suppresses PD-L1 levels in Lewis lung cancer cells and increases the proportions and activities of spleen-derived CD8 + T cells (Pan et al., 2020b). However, the underlying mechanism of QFJDD was unclear. From an immunoadjuvant therapy perspective, we reveal six core bioactive compounds of QFJDD and eight key genes involved in regulating PD-L1 expression and the proportion of CD8 + PD-1 + T cells in LUAD. HIF1A (hypoxia inducible factor 1, alpha subunit) is a basic helix-loop-helix transcription factor that decreases the survival and cytolytic activities of CD8 + CTLs, and promotes the expression of immune checkpoint inhibition molecules. In LUAD, PD-L1, a direct target of HIF1A, is positively correlated with HIF1A expression . Blocking HIF1A and PD-L1 enhance T cell activities (Noman et al., 2014). UPLC-Q-TOF-MS and molecular docking revealed that bioactive components of QFJDD acting on HIF1A were quercetin, luteolin, kaempferol, wogonin, baicalein and acacetin (Table 5, Table 6). According to previous studies, quercetin, luteolin and baicalein inhibit the expression of PD-L1 and restore the destruction of tumor cells by T cells (Ke et al., 2019;Jiang et al., 2021;Jing et al., 2021). Moreover, kaempferol significantly inhibits PD-1/PD-L1 interaction (Kim et al., 2020). We found that A549 cells exposed to QFJDD-containing serum had significantly reduced mRNA expressions of HIF1A and CD274 ( Figures 8C,H). KEGG analysis showed that HIF1A was a key gene in HIF-1 signaling pathway involved in regulating PD-L1 expression ( Figure 5B). These findings suggest that QFJDD suppresses PD-L1 expression in LUAD by inhibiting the HIF-1 signaling pathway and restoring CD8 + T cell activities.
EGFR (epidermal growth factor receptor) is an important oncogenic signaling pathway in NSCLC. It can directly or indirectly drive PD-L1 overexpression (Li et al., 2018) and affect the abundance of CD8 + T infiltration in tumor tissues (Zhao et al., 2020). In EGFR mutant NSCLC, activated EGFR induced PD-L1 expression through PI3K/AKT1 and MAPK signaling pathways (Luo et al., 2021), which is in accord with KEGG pathway analysis ( Figure 5). In recent years, the inhibitory effects of quercetin, luteolin, kaempferol, wogonin, baicalein and acacetin on EGFR have been reported in literature (Wenzel et al., 2001;Hong et al., 2014;Liu et al., 2016b;Yao et al., 2016;Tian et al., 2021;Ganthala et al., 2022). Furthermore, JUN (jun proto-oncogene, c-Jun) is the key component of dimeric transcription factor AP-1 involved in regulation of PD-L1 expression and activation of CD8 + T cells (Atsaves et al., 2019;Papavassiliou and Musti, 2020). It is also a target of the MAPK signaling cascade (Sumimoto et al., 2016;Zerdes et al., 2018). According to this study, the QFJDDcontaining serum significantly suppressed EGFR, JUN and CD274 mRNA expressions in A549 cells ( Figures 8A,E,H). Expressions of MAPK1 in the QFJDD group were Frontiers in Pharmacology frontiersin.org significantly high relative to the control group ( Figure 8G). Green et al. (2012) reported that the PD-L1 enhancer can bind AP-1 components and increase PD-L1 promoter activities in cHL Reed-Sternberg cells. The IHC assay showed that IRS of PD-L1 in the H d group was lower than in the Model group ( Figure 9F). The percentage of CD8 + PD-1 + T derived from spleen tissues of LLC bearing mice were apparently increased following treatment with high-dose QFJDD or BMS-202 ( Figure 9M). In accordance with previous studies, quercetin, luteolin, kaempferol, wogonin, baicalein, and acacetin inhibit the activation of c-Jun N-terminal kinase and attenuate the activation of the AP-1 transcription factor (Jang et al., 2008;Chen et al., 2008;Fong et al., 2010;Chen et al., 2012;Kim et al., 2012;Chen et al., 2013;Wang et al., 2014). Furthermore, kaempferol blocks the interactions between PD-1 and PD-L1 (Kim et al., 2020) whereas wogonin can down-regulate mRNA expression of ERK2 in melanoma cells (Chen et al., 2017). Thus, QFJDD suppresses PD-L1 expression in LUAD by inhibiting mRNA transcriptions of EGFR and JUN, and increasing the proportions of CD8 + PD-1 + T. NF-κB (Nuclear factor kappa-light-chain-enhancer of activated B cells) is synthesized in the cytoplasm and binds to IκB to form an inactive complex (Zhang et al., 2017a). The heterodimer p50/p65 is a classic representative of the NF-κB family encoded by NFKB1 and RELA genes, respectively. IκBα is encoded by NFKBIA and is phosphorylated and degraded by IκB kinase (IKK), leading to activation of free p-p50/p65 and nuclear translocation (Antonangeli et al., 2020). The in vitro cell assay revealed that following treatment with QFJDD-containing serum, mRNA expressions of RELA and NFKBIA in A549 cells were suppressed while AKT1 levels were increased ( Figures 8B,D,F). Quercetin, luteolin, kaempferol, wogonin, baicalein, and acacetin have been shown to suppress the NF-κB pathway (Chen et al., 2012;Kim et al., 2014;Sikder et al., 2014;Li et al., 2016;Tuorkey, 2016;Zhang et al., 2017b), consistent with molecular docking and qRT-PCR outcomes. The KEGG pathway analysis showed that NFκB is associated with PD-L1, and was also related to T cells apoptosis ( Figure 5B). The PI3K/Akt/ mTOR signaling pathway has been shown to affect immunity by regulating PD-L1 expression. Besides, IKK can be activated by AKT1 . The NF-κB signal is also associated with PD-L1 and CD8 + T (Antonangeli et al., 2020). These results suggest that QFJDD downregulates mRNA expression level of PD-L1 by suppressing p65 synthesis and inhibiting IκBα phosphorylation. In addition, the increase in CD8 + PD-1 + T % was also related to inhibition of the NFκB signaling pathway by QFJDD.

Conclusion
This study identified six core bioactive compounds in QFJDD. Then, eight key genes, including EGFR, JUN, RELA, HIF1A, NFKBIA, MAPK1, AKT1, and MAPK14, were identified by mapping the bioactive compounds of QFJDD to the targets of hsa05235 pathway. Meanwhile, EGFR, HIF-1, JUN, and NFκB signaling pathways were shown to be involved in regulating PD-L1 expression and CD8 + PD-1 + T% in LUAD. Moreover, the therapeutic potential of QFJDD holds tremendous promise for five areas, including virus, cancer, immunity, endocrine system, and cardiovascular system. Our findings provide a scientific basis for clinical applications of QFJDD in LUAD treatment.

Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author.