# Network Pharmacology-Based Study on the Active Ingredients and Mechanism of Pan Ji Sheng Traditional Chinese Medicine Formula in the Treatment of Inflammation

Authors: Shiji Wu; Hongliang Jiang; Zongwen Chen; Weining Lu; Qin Chen
Journal: Evidence-Based Complementary and Alternative Medicine
 (2022)

---

## Abstract

Background. Pan Ji Sheng Formula is a Chinese medicine formula that enables heat-free detoxification as well as anti-inflammatory and immune-boosting properties. This formula contains eight herbs. Its underlying mechanism is unknown. The bioactive ingredients were screened in our work, and the mechanism of this formula was investigated. Methods. Using traditional Chinese medicine systems pharmacology database and analysis platform (TCMSP), ingredients in Pan Ji Sheng Chinese medicine formula were screened, and we selected the main bioactive ingredients for web-based research. The targets of bioactive ingredients are primarily obtained from the SwissTargetPrediction and TCMSP databases, and the text mining method is used. STRING and Cytoscape were then used to examine the protein-protein interaction (PPI) networks. To explore the biological function and related pathways, functional annotation and pathway analysis were performed. Results. This research discovered 96 bioactive ingredients. Then, 215 potential targets of bioactive ingredients were screened. Through the analysis of the PPI network, we discovered 25 key target genes, which can be described as hub target genes regulated by bioactive ingredients. Bioactive ingredients primarily regulate CASP3, AKT1, JUN, and other proteins. The formula works synergistically to enhance immune response and antiinfection by regulating immune-related pathways, TNF signaling pathways, and apoptosis. Conclusions. A variety of bioactive ingredients in the formula could play roles in regulating CASP3, AKT1, and other genes in immune, infection, apoptosis, and tumor-related signaling pathways. Our data point the way forward for future studies on the mechanism of action of this formula.


---

## Body

## 1. Introduction

The climate in China’s Lingnan region is standard subtropical. Summers are hot, rainy, as well as wet [1]. Furthermore, Cantonese people prefer to eat fried, dry, and hot foods. It is easy to make people “heat” and “dampness” due to the hot and humid climate, poor diet, and insufficient sleep [2, 3]. The symptoms of “heat” contain fever, thirst, sweating, fatigue, yellow urine, and yellow tongue. The common symptoms of “dampness” contain head pain, chest tightness, sluggishness, and sore or swollen joints. “Heat” and “dampness” are considered to be the cause of many inflammatory disease, cancer, and metabolic disorders [2].Inflammation is a pathological defense response and it is also the most important protective response [4]. In modern western medicine, clinical experimental data show that the current conventional treatment for inflammation is anti-inflammatory drugs and antibiotic drugs [5, 6]. Nonsteroidal anti-inflammatory drugs (NSAIDs) are extensively used to reduce inflammation [7]. NSAIDs, such as aspirin and ibuprofen, are effective by inhibiting cyclooxygenase (COX) activity, thereby suppressing inflammatory responses [8]. Although it is effective, some anti-inflammatory drugs can lead to some side effects, such as gastrointestinal damage, gastrointestinal bleeding, and cardiovascular risk [9, 10]. The long-term use of antibiotic drugs can also lead to drug-resistance and seriously affect the treatment effect [11]. Traditional Chinese medicine (TCM) has the advantages of long efficacy and safety, so it is necessary to excavate the TCM compound formulas for treating inflammation.The ancestors attempted to collect herbs for clearing heat and detoxification, and boiling water for drinking to eliminate the “heat” in order to get rid of dampness and heat and adapt to the environment. Since this type of herbal medicine was safe to drink, it gradually spread among the people [12, 13]. People gradually dig up various therapeutic properties of traditional Chinese medicine substances under the research of ancient and modern science, and make formulas with heat-clearing and detoxification features with honeysuckle, Scutellaria baicalensis, chrysanthemum, isatis root, and other traditional Chinese medicines, so as to enhance immune response and alleviate problems such as getting angry and heavy moisture caused by improper diet and lack of sleep [14, 15]. TCM (traditional Chinese medicine) is a type of traditional medicine. TCM is still a vital resource with such a long history. TCM can still influence the advancement of modern medicine [16, 17]. The Pan Ji Sheng formula, which contains eight different herbs, is the subject of this research: Microctis Folium (the leaves of Microcos paniculata), Polygonum chinense (creeping smartweed), Ecliptae Herba (false daisy), Perilla Frutescens (the leaves of Beefsteak Plant), Isatidis Radix (the dried roots of the plant Isatis indigotica Fort or Isatis tinctoria L.), Chrysanthemi Flos (the flower of Chrysanthemum indicum Linne or Chrysanthemum morifolium Ramatuelle), Glycyrrhiza uralensis (Chinese liquorice, the root of Glycyrrhiza uralensis), and Chimonanthus salicifolius (wintersweet). All of these herbs are commonly used to treat diseases by clinicians. According to published research, these Chinese herbal medicines can prevent and treat diseases by utilizing a wide range of chemical components and multiple targets [18–21]. For example, isatis root lectin can directly kill influenza viruses by blocking the expression of nuclear proteins of new influenza viruses [22]; at the same time, nucleoside components such as uridine, guanosine, and adenosine can interfere with the synthesis of viral nucleic acid and perform critical roles for influenza virus defense [23], and polysaccharides have immunomodulatory effects and play indirect roles for influenza virus defense [24].There is, however, no systematic research report on the specific formula and network mechanism of the formula’s effects of clearing heat, detoxifying, anti-inflammatory, and enhancing immune response. Now, researchers have realized the “one key, one lock” model is insufficient for deciphering drug effects, particularly in complex diseases [25]. Network pharmacology is a new technology that uses the receptor theory and biological network technology to elucidate drug action mechanisms [26]. Its research mode of “multicomponent network target action” opens up a new research field and its compound prescriptions with multicomponent and multitarget synergy [27]. Furthermore, the rapid development of biomedical data, such as the TCMSP (traditional Chinese medicine system pharmacology database and analysis platform), has facilitated such research [28]. As a result, web-based pharmacological analysis can provide us with a thorough understanding of the significance of each component, target, and pathway. Based on the research concept of traditional Chinese medicine’s multicomponent and multitarget effect, this study explains the biological mechanism of clearing heat, detoxifying, anti-inflammatory, and enhancing immune response by using the network pharmacology technology and analyzing the target characteristics, biological function, and pathway of the Pan Ji Sheng formula. Our research provides a scientific basis for experimental research and product development.


## 2. Methods

### 2.1. Screening of Bioactive Ingredients

Through TCMSP, we search the relevant information about the bioactive ingredients in eight herbals in Pan Ji Sheng formula and screen the qualified compounds as the formula’s active ingredients. The screening conditions are oral bioavailability (OB) ≥ 30%, number of hydrogen bond donors (Hdon) < 5, lipid water partition coefficient (Alogp) < 5, number of hydrogen bond receptors (HACC) < 10, intestinal epithelial permeability (Caco-2) > 0, drug class (DL) ≥ 0.18, and drug half-life (HL) ≥ 4. We obtained bioactive ingredients of six herbals (Microctis Folium, Ecliptae Herba, Perilla Frutescens, Isatidis Radix, Chrysanthemi Flos, and Glycyrrhiza uralensis) from the TCMSP database. There is no information about Polygonum chinense and Chimonanthus salicifolia in the TCMSP database, so we search the literature for bioactive ingredients of these two herbals, then test OB ≥ 30% and DL ≥ 0.18 in TCMSP to determine the active ingredients.


### 2.2. Target Prediction of Bioactive Ingredients

The formula’s bioactive ingredients were imported to TCMSP to obtain information on ingredient-target interaction. Second, we use the Swiss Target Prediction online analysis tool to predict the active ingredient’s targets, screen potential targets, extract the names of the target genes, and build the chemical ingredient-target interaction network. The specific method is to convert all ingredients into standard smiles format and import the smiles format file into the Swiss Target Prediction online analysis platform [29], set the species to “Homo sapiens,” and set Probability ≥0.7, and export the target data in the CSV format.The target genes were then imported to the UniProt database to confirm their gene names. Through computer research, this study obtained the list of target genes for the traditional Chinese medicine Pan Ji Sheng formula.


### 2.3. Construction of the Protein-Protein Interaction (PPI) Network

We import target genes into STRING [30] and set the species to “Homo sapiens (human)” and use a confidence level of 0.9 to build the target interaction network (PPI). We hide the discrete points in the network, then export the results to a TSV file and import it to Cytoscape 3.9.1 [31]. Cytoscape was then used to construct the target’s PPI network.Then, in Cytoscape, the MCODE and Cytohubba plug-ins were used to extract the functional modules and top 25 hub genes of the PPI network, respectively.


### 2.4. Gene Ontology (GO) Functional Annotation and KEGG Pathway Analysis

All screened target genes were entered into the Metascape platform for enrichment analysis [32].The hub targets were imported into the David database to clarify their function and role in signal transduction. GO biological process enrichment analysis and KEGG signal pathway analysis are carried out. The enrichment analysis results are enhanced with the R program package and displayed in the form of a bubble diagram.


### 2.5. Construction of the Bioactive Ingredients-Hub Target Network

Cytoscape 3.9.1 software was used to build the bioactive ingredients-hub target network. In this network, nodes represent bioactive ingredients and hub targets.


#### 2.5.1. Hub Target-GO BP/Pathway/Disease Network

Use Cytoscape 3.9.1 to build the network model. Nodes represent hub targets, pathways, and diseases, and edges represent interactions between these nodes.


## 2.1. Screening of Bioactive Ingredients

Through TCMSP, we search the relevant information about the bioactive ingredients in eight herbals in Pan Ji Sheng formula and screen the qualified compounds as the formula’s active ingredients. The screening conditions are oral bioavailability (OB) ≥ 30%, number of hydrogen bond donors (Hdon) < 5, lipid water partition coefficient (Alogp) < 5, number of hydrogen bond receptors (HACC) < 10, intestinal epithelial permeability (Caco-2) > 0, drug class (DL) ≥ 0.18, and drug half-life (HL) ≥ 4. We obtained bioactive ingredients of six herbals (Microctis Folium, Ecliptae Herba, Perilla Frutescens, Isatidis Radix, Chrysanthemi Flos, and Glycyrrhiza uralensis) from the TCMSP database. There is no information about Polygonum chinense and Chimonanthus salicifolia in the TCMSP database, so we search the literature for bioactive ingredients of these two herbals, then test OB ≥ 30% and DL ≥ 0.18 in TCMSP to determine the active ingredients.


## 2.2. Target Prediction of Bioactive Ingredients

The formula’s bioactive ingredients were imported to TCMSP to obtain information on ingredient-target interaction. Second, we use the Swiss Target Prediction online analysis tool to predict the active ingredient’s targets, screen potential targets, extract the names of the target genes, and build the chemical ingredient-target interaction network. The specific method is to convert all ingredients into standard smiles format and import the smiles format file into the Swiss Target Prediction online analysis platform [29], set the species to “Homo sapiens,” and set Probability ≥0.7, and export the target data in the CSV format.The target genes were then imported to the UniProt database to confirm their gene names. Through computer research, this study obtained the list of target genes for the traditional Chinese medicine Pan Ji Sheng formula.


## 2.3. Construction of the Protein-Protein Interaction (PPI) Network

We import target genes into STRING [30] and set the species to “Homo sapiens (human)” and use a confidence level of 0.9 to build the target interaction network (PPI). We hide the discrete points in the network, then export the results to a TSV file and import it to Cytoscape 3.9.1 [31]. Cytoscape was then used to construct the target’s PPI network.Then, in Cytoscape, the MCODE and Cytohubba plug-ins were used to extract the functional modules and top 25 hub genes of the PPI network, respectively.


## 2.4. Gene Ontology (GO) Functional Annotation and KEGG Pathway Analysis

All screened target genes were entered into the Metascape platform for enrichment analysis [32].The hub targets were imported into the David database to clarify their function and role in signal transduction. GO biological process enrichment analysis and KEGG signal pathway analysis are carried out. The enrichment analysis results are enhanced with the R program package and displayed in the form of a bubble diagram.


## 2.5. Construction of the Bioactive Ingredients-Hub Target Network

Cytoscape 3.9.1 software was used to build the bioactive ingredients-hub target network. In this network, nodes represent bioactive ingredients and hub targets.


### 2.5.1. Hub Target-GO BP/Pathway/Disease Network

Use Cytoscape 3.9.1 to build the network model. Nodes represent hub targets, pathways, and diseases, and edges represent interactions between these nodes.


## 2.5.1. Hub Target-GO BP/Pathway/Disease Network

Use Cytoscape 3.9.1 to build the network model. Nodes represent hub targets, pathways, and diseases, and edges represent interactions between these nodes.


## 3. Results

### 3.1. Screening of Bioactive Ingredients of the Pan Ji Sheng Formula

The bioactive ingredients of eight Chinese herbal medicines from the Pan Ji Sheng formula were screened from the TCMSP platform in this study. Because there is no relevant information on the TCMSP platform forPolygonum chinense and Chimonanthus salicifolia, we obtained the active components of these two herbals through literature retrieval and then tested whether they meet the standards of oral bioavailability (OB) ≥ 30 percent and drug class (DL) ≥ 0.18 in TCMSP. We obtained the active components of the other six herbals from TCMSP. In total, this study screened 96 active ingredients from eight herbals in the Pan Ji Sheng formula (Table 1).Table 1
Herbal and bioactive ingredients of Pan Ji Sheng formula.
HerbalsMolecule namesMicroctis foliumIsorhamnetinKaempferol4′,5-DihydroxyflavoneKaempferolQuercetinPolygonum chinense3-O-Methylellagic acidKaempferol-7-O-glucoside3,3′-Di-O-methylellagic acidProtocatechuic acidIsorhamnetinLuteolinAcacetinEcliptae herbaButin1,3,8,9-Tetrahydroxybenzofurano [3,2-c] chromen-6-one3′-O-MethylorobolPratenseinDemethylwedelolactoneWedelolactoneLuteolinPerilla frutescensLuteolinAcacetinEupatorinDinatinQuindolineHydroxyindirubinIndigo(2Z)-2-(2-Oxoindolin-3-ylidene) indolin-3-oneIsatidis radix2-(9-((3-Methyl-2-oxopent-3-en-1-yl) oxy)-2-oxo-1,2,8,9-tetrahydrofuro [2,3-h] quinolin-8-yl) propan-2-yl acetateDFV(E)-2-[(3-Indole) cyanomethylene-]-3-indolinoneneohesperidin_qtSinensetin6-(3-Oxoindolin-2-ylidene) indolo[2,1-b]quinazolin-12-one(E)-3-(3,5-Dimethoxy-4-hydroxy-benzylidene)-2-indolinone(E)-3-(3,5-Dimethoxy-4-hydroxyb-enzylidene)-2-indolinone3-[(3,5-Dimethoxy-4-oxo-1-cyclohexa-2,5-dienylidene)methyl]-2,4-dihydro-1H-pyrrolo[2,1-b] quinazolin-9-one[(1S,5S,7S)-7-Acetoxy-5-isopropenyl-2,8-dimethylene-cyclodecyl] acetateAcacetinChryseriolIsorhamnetinChrysanthemi flosKaempferol5,7-Dihydroxy-2-(3-hydroxy-4-methoxyphenyl) chroman-4-oneLuteolinEupatorinDiosmetinNaringeninArtemetinJaranolIsorhamnetinFormononetinLicoriceCalycosinKaempferolLicochalcone aInermineDFVGlycyrolMedicarpinLupiwighteone7-Methoxy-2-methyl isoflavoneNaringeninGlyasperin BGlyasperin FIsotrifoliol(E)-1-(2,4-Dihydroxyphenyl)-3-(2,2-dimethylchromen-6-yl) prop-2-en-1-one(2S)-6-(2,4-Dihydroxyphenyl)-2-(2-hydroxypropan-2-yl)-4-methoxy-2,3-dihydrofuro [3,2-g] chromen-7-oneSemilicoisoflavone BGlepidotin AGlepidotin BGlypallichalcone8-(6-Hydroxy-2-benzofuranyl)-2,2-dimethyl-5-chromenolLicochalcone BLicochalcone GLicoriconeGancaonin AGancaonin B3-(3,4-Dihydroxyphenyl)-5,7-dihydroxy-8-(3-methylbut-2-enyl) chromone5,7-Dihydroxy-3-(4-methoxyphenyl)-8-(3-methylbut-2-enyl) chromone2-(3,4-Dihydroxyphenyl)-5,7-dihydroxy-6-(3-methylbut-2-enyl) chromoneLicocoumaroneLicoisoflavoneLicoisoflavone BLicoisoflavanoneShinpterocarpin(E)-3-[3,4-Dihydroxy-5-(3-methylbut-2-enyl)phenyl]-1-(2,4-dihydroxyphenyl) prop-2-en-1-oneGlyzaglabrinGlabraninGlabrone1,3-Dihydroxy-9-methoxy-6-benzofurano[3,2-c] chromenone1,3-Dihydroxy-8,9-dimethoxy-6-benzofurano[3,2-c] chromenoneEurycarpin ASigmoidin-B(2R)-7-Hydroxy-2-(4-hydroxyphenyl) chroman-4-one(2S)-7-Hydroxy-2-(4-hydroxyphenyl)-8-(3-methylbut-2-enyl) chroman-4-oneIsoglycyrolIsolicoflavonolHMO1-MethoxyphaseollidinQuercetin der.6-Prenylated eriodictyol7-Acetoxy-2-methylisoflavone8-Prenylated eriodictyolGancaonin GGancaonin HLicoagrocarpinGlyasperins MLicoagroisoflavoneOdoratinPhaseolXambioonaChimonanthus salicifoliusLuteolin-5-O-glucosideQuercetinKaempferol


### 3.2. Screening of Target Genes

Target genes of bioactive components were obtained using the TCMSP platform and Swiss target prediction screening. After removing the repeated target genes, we obtained a total of 214 target genes in this study (Table2). For details of target genes, see Table S1.Table 2
Potential target genes of bioactive ingredients of Pan Ji Sheng formula.
No.Target gene namesString Id1NOS29606.ENSP000003272512PTGS19606.ENSP000003546123ESR19606.ENSP000004053304AR9606.ENSP000003638225PPARG9606.ENSP000002878206PTGS29606.ENSP000003564387PTPN19606.ENSP000003606838ESR29606.ENSP000003439259DPP49606.ENSP0000035373110MAPK149606.ENSP0000022979511GSK3B9606.ENSP0000032480612HSP90AA19606.ENSP0000033515313CDK29606.ENSP0000026697014PIK3CG9606.ENSP0000035212115PKIA9606.ENSP0000037969616PRSS19606.ENSP0000030872017PIM19606.ENSP0000036260818CCNA29606.ENSP0000027402619NCOA29606.ENSP0000039996820CALM29606.ENSP0000027229821PYGM9606.ENSP0000016413922PPARD9606.ENSP0000031092823CHEK19606.ENSP0000038864824AKR1B19606.ENSP0000028593025NCOA19606.ENSP0000038521626F79606.ENSP0000036473127F29606.ENSP0000030854128NOS39606.ENSP0000029749429ACHE9606.ENSP0000030321130GABRA19606.ENSP0000039309731MAOB9606.ENSP0000036730932GRIA29606.ENSP0000029652633RELA9606.ENSP0000038427334XDH9606.ENSP0000036872735NCF19606.ENSP0000028947336OLR19606.ENSP0000030912437PGR9606.ENSP0000032512038CHRM19606.ENSP0000030649039GABRA29606.ENSP0000042182840SLC6A29606.ENSP0000021983341CHRM29606.ENSP0000039974542ADRA1B9606.ENSP0000030666243TOP2A9606.ENSP0000041153244IKBKB9606.ENSP0000043068445AKT19606.ENSP0000045182846BCL29606.ENSP0000038118547BAX9606.ENSP0000029328848CD40LG9606.ENSP0000035966349JUN9606.ENSP0000036026650AHSA19606.ENSP0000021647951CASP39606.ENSP0000031103252MAPK89606.ENSP0000037897453MMP19606.ENSP0000032278854STAT19606.ENSP0000035439455CDK19606.ENSP0000037869956HMOX19606.ENSP0000021611757CYP3A49606.ENSP0000033791558CYP1A19606.ENSP0000036905059ICAM19606.ENSP0000026483260SELE9606.ENSP0000033173661VCAM19606.ENSP0000029472862NR1I29606.ENSP0000033652863CYP1B19606.ENSP0000047856164ALOX59606.ENSP0000036351265HAS29606.ENSP0000030699166AHR9606.ENSP0000024205767PSMD39606.ENSP0000026463968SLC2A49606.ENSP0000032093569NR1I39606.ENSP0000035695970INSR9606.ENSP0000030383071DIO19606.ENSP0000035464372GSTM19606.ENSP0000031146973GSTM29606.ENSP0000024133774AKR1C39606.ENSP0000036992775SLPI9606.ENSP0000034208276NOX49606.ENSP0000026331777AVPR29606.ENSP0000035180578MAOA9606.ENSP0000034068479IGF1R9606.ENSP0000026803580FLT39606.ENSP0000024145381CYP19A19606.ENSP0000037968382EGFR9606.ENSP0000027549383CA29606.ENSP0000028537984AURKB9606.ENSP0000031395085DRD49606.ENSP0000017618386ADORA19606.ENSP0000035620587CA79606.ENSP0000034565988GLO19606.ENSP0000036246389MPO9606.ENSP0000022527590PIK3R19606.ENSP0000042805691ADORA2A9606.ENSP0000033663092DAPK19606.ENSP0000038613593PYGL9606.ENSP0000021639294CA19606.ENSP0000043065695SRC9606.ENSP0000036268096PTK29606.ENSP0000034118997HSD17B29606.ENSP0000019993698KDR9606.ENSP0000026392399MMP139606.ENSP00000260302100CA129606.ENSP00000178638101CA139606.ENSP00000318912102CA99606.ENSP00000367608103GPR359606.ENSP00000411788104ERBB29606.ENSP00000269571105CCND19606.ENSP00000227507106CDK49606.ENSP00000257904107PDGFRB9606.ENSP00000261799108FLT49606.ENSP00000261937109CCNA19606.ENSP00000255465110PLK19606.ENSP00000300093111CA69606.ENSP00000366654112CA149606.ENSP00000358107113CSNK2A19606.ENSP00000217244114MET9606.ENSP00000317272115CA49606.ENSP00000300900116PLK49606.ENSP00000270861117TEK9606.ENSP00000369375118TNF9606.ENSP00000398698119IL29606.ENSP00000226730120RPS6KA39606.ENSP00000368884121CD389606.ENSP00000226279122PDE5A9606.ENSP00000347046123NQO29606.ENSP00000369822124ADRA2C9606.ENSP00000386069125ALDH29606.ENSP00000261733126NMUR29606.ENSP00000255262127ADRA2A9606.ENSP00000280155128SLC29A19606.ENSP00000377424129AURKA9606.ENSP00000216911130CA5A9606.ENSP00000309649131BACE19606.ENSP00000318585132MAP3K89606.ENSP00000263056133BRAF9606.ENSP00000288602134BCL2L19606.ENSP00000302564135CDKN1A9606.ENSP00000384849136CASP99606.ENSP00000330237137MMP29606.ENSP00000219070138MMP99606.ENSP00000361405139MAPK19606.ENSP00000215832140IL109606.ENSP00000412237141RB19606.ENSP00000267163142CDK49606.ENSP00000257904143IL69606.ENSP00000385675144TP539606.ENSP00000269305145NFKBIA9606.ENSP00000216797146TOP19606.ENSP00000354522147MDM29606.ENSP00000258149148APP9606.ENSP00000284981149PCNA9606.ENSP00000368458150CASP79606.ENSP00000358327151MCL19606.ENSP00000358022152BIRC59606.ENSP00000301633153CCNB19606.ENSP00000256442154TYR9606.ENSP00000263321155IFNG9606.ENSP00000229135156IL49606.ENSP00000231449157XIAP9606.ENSP00000360242158PTGES9606.ENSP00000342385159NUF29606.ENSP00000271452160ADCY29606.ENSP00000342952161ADRB29606.ENSP00000305372162PDE3A9606.ENSP00000351957163CASP89606.ENSP00000351273164FASN9606.ENSP00000304592165FASLG9606.ENSP00000356694166RXRA9606.ENSP00000419692167LACTBL19606.ENSP00000402297168SCN5A9606.ENSP00000410257169F109606.ENSP00000364709170RHO9606.ENSP00000296271171KCNH29606.ENSP00000262186172KCNMA19606.ENSP00000286628173SLC6A49606.ENSP00000261707174CHRNA79606.ENSP00000407546175PPP3CA9606.ENSP00000378323176MAPK39606.ENSP00000263025177LDLR9606.ENSP00000454071178BAD9606.ENSP00000378040179SOD19606.ENSP00000270142180MTTP9606.ENSP00000427679181APOB9606.ENSP00000233242182PLB19606.ENSP00000330442183HMGCR9606.ENSP00000287936184UGT1A89606.ENSP00000304845185PPARA9606.ENSP00000385523186SREBF19606.ENSP00000348069187GSR9606.ENSP00000221130188ABCC19606.ENSP00000382342189ADIPOQ9606.ENSP00000389814190SOAT29606.ENSP00000301466191AKR1C19606.ENSP00000370254192GOT19606.ENSP00000359539193ABAT9606.ENSP00000379845194CES19606.ENSP00000353720195SOAT19606.ENSP00000356591196ADRA1D9606.ENSP00000368766197SLC6A39606.ENSP00000270349198SIRT19606.ENSP00000212015199ATP5B9606.ENSP00000262030200MT-ND69606.ENSP00000354665201HSD3B29606.ENSP00000445122202HSD3B19606.ENSP00000358421203STAT39606.ENSP00000264657204EIF69606.ENSP00000363574205FOSL29606.ENSP00000264716206CHRM39606.ENSP00000255380207OPRM19606.ENSP00000394624208DRD19606.ENSP00000377353209CHRM59606.ENSP00000372750210CHRM49606.ENSP00000409378211HTR2A9606.ENSP00000437737212MAPK109606.ENSP00000352157213OPRD19606.ENSP00000234961214ADRB19606.ENSP00000358301215LTA4H9606.ENSP00000228740


### 3.3. Enrichment Analysis of All Target Genes

Using the Metascape website, this study firstly discovered relevant significantly enriched GO/KEGG terms for all target genes. Figure1 depicts the findings of the analysis. Many target genes are enriched in cancer and lipid metabolism-related pathways (Figures 1(a) and 1(b)). A subset of enriched terms was chosen and rendered as a network plot to further capture the relationships between the terms (Figure 1(c)).Figure 1
Enrichment analysis for bioactive ingredient targets by Metascape website. (a, b) Top 20 clusters with their representative enriched terms. (c) :Enrichment heatmap of the selected GO parents.
(a)(b)(c)We also analyzed related diseases and expression patterns of all target genes through Metascape, as shown in Figure2. Diabetes, reperfusion injury, and fatty liver disease are the three most common diseases associated with target genes. The tissues that expressed the target genes were the lung and liver. According to preliminary findings, the target gene may be linked to lung and liver diseases.Figure 2
Related diseases and expression patterns of all target genes. (a) The summary of enrichment analysis in Disgenet. (b) The summary of enrichment analysis in PaGenBase.
(a)(b)


#### 3.3.1. PPI Network for All Targets

We upload the names of all target genes to STRING. According to network statistics, the number of nodes is 214, the number of edges is 3057, and the average node degree is 28.6. The expected number of edges is 1173, and the local clustering coefficient is 0.583. We discovered that the network had far more interactions than expected. This suggests that the target proteins as a group are at least partially biologically connected.Using Cytoscape 3.9.1, we constructed a PPI network (Figure3(a)). Then, using the Cytoscape plug-in “cytohubba,” we analyzed hub targets and chose the top 25 target genes as hub genes (Figure 3(b)). CASP3, AKT1, Jun, STAT3, TP53, MMP9, BCL2l1, SRC, and other proteins. The higher the rank, the more important these target genes are in disease treatment. Hub targets are painted red and located at the center of the network for further analysis and research.Figure 3
PPI network of all target genes. (a) PPI network, colored and in the middle are 25 hub genes. (b) Top 25 genes in the network ranked by the MCC method in “Cytohubba”.
(a)(b)We also used the Cytoscape plug-in “MCODE” to examine the PPI network clusters and modules of all target genes (Figure4). The PPI network is divided into six clusters, with 25 hub target genes located in Cluster 1, indicating that hub genes have biological function relevance and may play a synergistic role.Figure 4
Clusters 1–6 in the PPI network. Among them, 25 hub genes are painted red and orange.


### 3.4. Herbal-Key Bioactive Ingredient-Hub Target Network

After obtaining the hub target genes, we analyzed the active ingredients corresponding to these 25 hub genes, which are named as key bioactive ingredients. For more information, see Table S2. The network of herbal-key bioactive ingredient-hub targets was constructed using Cytoscape 3.9.1 (Figure5). In addition to Perilla frutescens, the other seven Chinese herbal medicines have three or more corresponding key bioactive ingredients. Some hub genes are affected by multiple bioactive ingredients at the same time. The primary targets of the active ingredients are MAPK14, HSP90AA1, PTGS2, and ESR1. These genes may be the primary targets of the formula.Figure 5
Herbal-key bioactive ingredient-hub target network.


### 3.5. GO Functional Annotation and KEGG Pathway Analysis

To investigate the biological processes engaged in hub targets, GO enrichment analysis and KEGG enrichment analysis on 25 hub genes were analyzed in the David website. The mechanism of action of the formula can be researched, based on the biological process regulated by the hub target.Beautify the enrichment analysis results withR (Figure 6). In total, 226 GO biological process enrichment results were obtained. Negative regulation of the apoptotic process, positive regulation of the nitric oxide biosynthetic process, and positive regulation of transcription from the RNA polymerase II promoter are the top three enrichment biological processes. As shown in Figure 6(a), the top 20 GO biological processes are represented in the form of a bubble diagram, where the size of the circle represents the enrichment of relevant targets in the pathway, and the darker the color of the circle represents the degree of enrichment of targets, indicating that the formula could have physiological effects by regulating these biological processes.Figure 6
GO and KEGG enrichment analysis of hub genes.
(a)(b)(c)(d)For KEGG pathway enrichment analysis, 25 hub targets were mapped into the David database. The species was defined as “human,” and a total of 94 pathways were obtained. As shown in Figure6(d), the top 20 pathways with high significance of KEGG enrichment results are closely related to the mechanism of the Pan Ji Sheng formula. The top five pathways include hepatitis B, pathways in cancer, TNF signaling pathway, toxoplasmosis, and toll-like receptor signaling pathway. The majority of these pathways are linked to the genes TP53, JUN, AKT1, MAPK14, HSP90AA1, and PTGS2.We also performed disease enrichment analysis to investigate diseases associated with hub targets. Figure7 shows the classification of diseases enriched in hub targets. The three major categories are cancer, infection, and immune system. Our findings indicate that the formula studied in this study may primarily target these diseases.Figure 7
Disease and disease class enrichment analysis of hub genes.
(a)(b)


#### 3.5.1. Hub Target-GO BP/Pathway/Disease Class Network

In order to demonstrate the biological process of the hub target and the relationship between the hub target and the pathway more clearly, the hub target-GO BP/pathway/disease class network was built with Cytoscape 3.9.1 software (Figure8).Figure 8
Hub target-GO BP/pathway network.The hub target is represented by the circle in the center of Figure8. The left and right sides of Figure 8 show the top 20 enriched biological processes and pathways, respectively. We can clearly understand the relationship between the targets and biological processes or pathways. MAPK14, hSP90AA1, and PTGS2 genes are associated with apoptotic biological processes, TNF signaling pathways, toll-like receptor signaling pathways, and cancer pathways. The formula could play a significant role by regulating these pathways.In order to demonstrate the link between the hub targets and diseases more clearly, Cytoscape 3.9.1 software was used to create a network of hub targets and diseases (Figure9). The genes MAPK14, HSP90AA1, PTGS2, and ESR1 have been linked to cancer, infection, and immune disease.Figure 9
Hub target-disease class network.


## 3.1. Screening of Bioactive Ingredients of the Pan Ji Sheng Formula

The bioactive ingredients of eight Chinese herbal medicines from the Pan Ji Sheng formula were screened from the TCMSP platform in this study. Because there is no relevant information on the TCMSP platform forPolygonum chinense and Chimonanthus salicifolia, we obtained the active components of these two herbals through literature retrieval and then tested whether they meet the standards of oral bioavailability (OB) ≥ 30 percent and drug class (DL) ≥ 0.18 in TCMSP. We obtained the active components of the other six herbals from TCMSP. In total, this study screened 96 active ingredients from eight herbals in the Pan Ji Sheng formula (Table 1).Table 1
Herbal and bioactive ingredients of Pan Ji Sheng formula.
HerbalsMolecule namesMicroctis foliumIsorhamnetinKaempferol4′,5-DihydroxyflavoneKaempferolQuercetinPolygonum chinense3-O-Methylellagic acidKaempferol-7-O-glucoside3,3′-Di-O-methylellagic acidProtocatechuic acidIsorhamnetinLuteolinAcacetinEcliptae herbaButin1,3,8,9-Tetrahydroxybenzofurano [3,2-c] chromen-6-one3′-O-MethylorobolPratenseinDemethylwedelolactoneWedelolactoneLuteolinPerilla frutescensLuteolinAcacetinEupatorinDinatinQuindolineHydroxyindirubinIndigo(2Z)-2-(2-Oxoindolin-3-ylidene) indolin-3-oneIsatidis radix2-(9-((3-Methyl-2-oxopent-3-en-1-yl) oxy)-2-oxo-1,2,8,9-tetrahydrofuro [2,3-h] quinolin-8-yl) propan-2-yl acetateDFV(E)-2-[(3-Indole) cyanomethylene-]-3-indolinoneneohesperidin_qtSinensetin6-(3-Oxoindolin-2-ylidene) indolo[2,1-b]quinazolin-12-one(E)-3-(3,5-Dimethoxy-4-hydroxy-benzylidene)-2-indolinone(E)-3-(3,5-Dimethoxy-4-hydroxyb-enzylidene)-2-indolinone3-[(3,5-Dimethoxy-4-oxo-1-cyclohexa-2,5-dienylidene)methyl]-2,4-dihydro-1H-pyrrolo[2,1-b] quinazolin-9-one[(1S,5S,7S)-7-Acetoxy-5-isopropenyl-2,8-dimethylene-cyclodecyl] acetateAcacetinChryseriolIsorhamnetinChrysanthemi flosKaempferol5,7-Dihydroxy-2-(3-hydroxy-4-methoxyphenyl) chroman-4-oneLuteolinEupatorinDiosmetinNaringeninArtemetinJaranolIsorhamnetinFormononetinLicoriceCalycosinKaempferolLicochalcone aInermineDFVGlycyrolMedicarpinLupiwighteone7-Methoxy-2-methyl isoflavoneNaringeninGlyasperin BGlyasperin FIsotrifoliol(E)-1-(2,4-Dihydroxyphenyl)-3-(2,2-dimethylchromen-6-yl) prop-2-en-1-one(2S)-6-(2,4-Dihydroxyphenyl)-2-(2-hydroxypropan-2-yl)-4-methoxy-2,3-dihydrofuro [3,2-g] chromen-7-oneSemilicoisoflavone BGlepidotin AGlepidotin BGlypallichalcone8-(6-Hydroxy-2-benzofuranyl)-2,2-dimethyl-5-chromenolLicochalcone BLicochalcone GLicoriconeGancaonin AGancaonin B3-(3,4-Dihydroxyphenyl)-5,7-dihydroxy-8-(3-methylbut-2-enyl) chromone5,7-Dihydroxy-3-(4-methoxyphenyl)-8-(3-methylbut-2-enyl) chromone2-(3,4-Dihydroxyphenyl)-5,7-dihydroxy-6-(3-methylbut-2-enyl) chromoneLicocoumaroneLicoisoflavoneLicoisoflavone BLicoisoflavanoneShinpterocarpin(E)-3-[3,4-Dihydroxy-5-(3-methylbut-2-enyl)phenyl]-1-(2,4-dihydroxyphenyl) prop-2-en-1-oneGlyzaglabrinGlabraninGlabrone1,3-Dihydroxy-9-methoxy-6-benzofurano[3,2-c] chromenone1,3-Dihydroxy-8,9-dimethoxy-6-benzofurano[3,2-c] chromenoneEurycarpin ASigmoidin-B(2R)-7-Hydroxy-2-(4-hydroxyphenyl) chroman-4-one(2S)-7-Hydroxy-2-(4-hydroxyphenyl)-8-(3-methylbut-2-enyl) chroman-4-oneIsoglycyrolIsolicoflavonolHMO1-MethoxyphaseollidinQuercetin der.6-Prenylated eriodictyol7-Acetoxy-2-methylisoflavone8-Prenylated eriodictyolGancaonin GGancaonin HLicoagrocarpinGlyasperins MLicoagroisoflavoneOdoratinPhaseolXambioonaChimonanthus salicifoliusLuteolin-5-O-glucosideQuercetinKaempferol


## 3.2. Screening of Target Genes

Target genes of bioactive components were obtained using the TCMSP platform and Swiss target prediction screening. After removing the repeated target genes, we obtained a total of 214 target genes in this study (Table2). For details of target genes, see Table S1.Table 2
Potential target genes of bioactive ingredients of Pan Ji Sheng formula.
No.Target gene namesString Id1NOS29606.ENSP000003272512PTGS19606.ENSP000003546123ESR19606.ENSP000004053304AR9606.ENSP000003638225PPARG9606.ENSP000002878206PTGS29606.ENSP000003564387PTPN19606.ENSP000003606838ESR29606.ENSP000003439259DPP49606.ENSP0000035373110MAPK149606.ENSP0000022979511GSK3B9606.ENSP0000032480612HSP90AA19606.ENSP0000033515313CDK29606.ENSP0000026697014PIK3CG9606.ENSP0000035212115PKIA9606.ENSP0000037969616PRSS19606.ENSP0000030872017PIM19606.ENSP0000036260818CCNA29606.ENSP0000027402619NCOA29606.ENSP0000039996820CALM29606.ENSP0000027229821PYGM9606.ENSP0000016413922PPARD9606.ENSP0000031092823CHEK19606.ENSP0000038864824AKR1B19606.ENSP0000028593025NCOA19606.ENSP0000038521626F79606.ENSP0000036473127F29606.ENSP0000030854128NOS39606.ENSP0000029749429ACHE9606.ENSP0000030321130GABRA19606.ENSP0000039309731MAOB9606.ENSP0000036730932GRIA29606.ENSP0000029652633RELA9606.ENSP0000038427334XDH9606.ENSP0000036872735NCF19606.ENSP0000028947336OLR19606.ENSP0000030912437PGR9606.ENSP0000032512038CHRM19606.ENSP0000030649039GABRA29606.ENSP0000042182840SLC6A29606.ENSP0000021983341CHRM29606.ENSP0000039974542ADRA1B9606.ENSP0000030666243TOP2A9606.ENSP0000041153244IKBKB9606.ENSP0000043068445AKT19606.ENSP0000045182846BCL29606.ENSP0000038118547BAX9606.ENSP0000029328848CD40LG9606.ENSP0000035966349JUN9606.ENSP0000036026650AHSA19606.ENSP0000021647951CASP39606.ENSP0000031103252MAPK89606.ENSP0000037897453MMP19606.ENSP0000032278854STAT19606.ENSP0000035439455CDK19606.ENSP0000037869956HMOX19606.ENSP0000021611757CYP3A49606.ENSP0000033791558CYP1A19606.ENSP0000036905059ICAM19606.ENSP0000026483260SELE9606.ENSP0000033173661VCAM19606.ENSP0000029472862NR1I29606.ENSP0000033652863CYP1B19606.ENSP0000047856164ALOX59606.ENSP0000036351265HAS29606.ENSP0000030699166AHR9606.ENSP0000024205767PSMD39606.ENSP0000026463968SLC2A49606.ENSP0000032093569NR1I39606.ENSP0000035695970INSR9606.ENSP0000030383071DIO19606.ENSP0000035464372GSTM19606.ENSP0000031146973GSTM29606.ENSP0000024133774AKR1C39606.ENSP0000036992775SLPI9606.ENSP0000034208276NOX49606.ENSP0000026331777AVPR29606.ENSP0000035180578MAOA9606.ENSP0000034068479IGF1R9606.ENSP0000026803580FLT39606.ENSP0000024145381CYP19A19606.ENSP0000037968382EGFR9606.ENSP0000027549383CA29606.ENSP0000028537984AURKB9606.ENSP0000031395085DRD49606.ENSP0000017618386ADORA19606.ENSP0000035620587CA79606.ENSP0000034565988GLO19606.ENSP0000036246389MPO9606.ENSP0000022527590PIK3R19606.ENSP0000042805691ADORA2A9606.ENSP0000033663092DAPK19606.ENSP0000038613593PYGL9606.ENSP0000021639294CA19606.ENSP0000043065695SRC9606.ENSP0000036268096PTK29606.ENSP0000034118997HSD17B29606.ENSP0000019993698KDR9606.ENSP0000026392399MMP139606.ENSP00000260302100CA129606.ENSP00000178638101CA139606.ENSP00000318912102CA99606.ENSP00000367608103GPR359606.ENSP00000411788104ERBB29606.ENSP00000269571105CCND19606.ENSP00000227507106CDK49606.ENSP00000257904107PDGFRB9606.ENSP00000261799108FLT49606.ENSP00000261937109CCNA19606.ENSP00000255465110PLK19606.ENSP00000300093111CA69606.ENSP00000366654112CA149606.ENSP00000358107113CSNK2A19606.ENSP00000217244114MET9606.ENSP00000317272115CA49606.ENSP00000300900116PLK49606.ENSP00000270861117TEK9606.ENSP00000369375118TNF9606.ENSP00000398698119IL29606.ENSP00000226730120RPS6KA39606.ENSP00000368884121CD389606.ENSP00000226279122PDE5A9606.ENSP00000347046123NQO29606.ENSP00000369822124ADRA2C9606.ENSP00000386069125ALDH29606.ENSP00000261733126NMUR29606.ENSP00000255262127ADRA2A9606.ENSP00000280155128SLC29A19606.ENSP00000377424129AURKA9606.ENSP00000216911130CA5A9606.ENSP00000309649131BACE19606.ENSP00000318585132MAP3K89606.ENSP00000263056133BRAF9606.ENSP00000288602134BCL2L19606.ENSP00000302564135CDKN1A9606.ENSP00000384849136CASP99606.ENSP00000330237137MMP29606.ENSP00000219070138MMP99606.ENSP00000361405139MAPK19606.ENSP00000215832140IL109606.ENSP00000412237141RB19606.ENSP00000267163142CDK49606.ENSP00000257904143IL69606.ENSP00000385675144TP539606.ENSP00000269305145NFKBIA9606.ENSP00000216797146TOP19606.ENSP00000354522147MDM29606.ENSP00000258149148APP9606.ENSP00000284981149PCNA9606.ENSP00000368458150CASP79606.ENSP00000358327151MCL19606.ENSP00000358022152BIRC59606.ENSP00000301633153CCNB19606.ENSP00000256442154TYR9606.ENSP00000263321155IFNG9606.ENSP00000229135156IL49606.ENSP00000231449157XIAP9606.ENSP00000360242158PTGES9606.ENSP00000342385159NUF29606.ENSP00000271452160ADCY29606.ENSP00000342952161ADRB29606.ENSP00000305372162PDE3A9606.ENSP00000351957163CASP89606.ENSP00000351273164FASN9606.ENSP00000304592165FASLG9606.ENSP00000356694166RXRA9606.ENSP00000419692167LACTBL19606.ENSP00000402297168SCN5A9606.ENSP00000410257169F109606.ENSP00000364709170RHO9606.ENSP00000296271171KCNH29606.ENSP00000262186172KCNMA19606.ENSP00000286628173SLC6A49606.ENSP00000261707174CHRNA79606.ENSP00000407546175PPP3CA9606.ENSP00000378323176MAPK39606.ENSP00000263025177LDLR9606.ENSP00000454071178BAD9606.ENSP00000378040179SOD19606.ENSP00000270142180MTTP9606.ENSP00000427679181APOB9606.ENSP00000233242182PLB19606.ENSP00000330442183HMGCR9606.ENSP00000287936184UGT1A89606.ENSP00000304845185PPARA9606.ENSP00000385523186SREBF19606.ENSP00000348069187GSR9606.ENSP00000221130188ABCC19606.ENSP00000382342189ADIPOQ9606.ENSP00000389814190SOAT29606.ENSP00000301466191AKR1C19606.ENSP00000370254192GOT19606.ENSP00000359539193ABAT9606.ENSP00000379845194CES19606.ENSP00000353720195SOAT19606.ENSP00000356591196ADRA1D9606.ENSP00000368766197SLC6A39606.ENSP00000270349198SIRT19606.ENSP00000212015199ATP5B9606.ENSP00000262030200MT-ND69606.ENSP00000354665201HSD3B29606.ENSP00000445122202HSD3B19606.ENSP00000358421203STAT39606.ENSP00000264657204EIF69606.ENSP00000363574205FOSL29606.ENSP00000264716206CHRM39606.ENSP00000255380207OPRM19606.ENSP00000394624208DRD19606.ENSP00000377353209CHRM59606.ENSP00000372750210CHRM49606.ENSP00000409378211HTR2A9606.ENSP00000437737212MAPK109606.ENSP00000352157213OPRD19606.ENSP00000234961214ADRB19606.ENSP00000358301215LTA4H9606.ENSP00000228740


## 3.3. Enrichment Analysis of All Target Genes

Using the Metascape website, this study firstly discovered relevant significantly enriched GO/KEGG terms for all target genes. Figure1 depicts the findings of the analysis. Many target genes are enriched in cancer and lipid metabolism-related pathways (Figures 1(a) and 1(b)). A subset of enriched terms was chosen and rendered as a network plot to further capture the relationships between the terms (Figure 1(c)).Figure 1
Enrichment analysis for bioactive ingredient targets by Metascape website. (a, b) Top 20 clusters with their representative enriched terms. (c) :Enrichment heatmap of the selected GO parents.
(a)(b)(c)We also analyzed related diseases and expression patterns of all target genes through Metascape, as shown in Figure2. Diabetes, reperfusion injury, and fatty liver disease are the three most common diseases associated with target genes. The tissues that expressed the target genes were the lung and liver. According to preliminary findings, the target gene may be linked to lung and liver diseases.Figure 2
Related diseases and expression patterns of all target genes. (a) The summary of enrichment analysis in Disgenet. (b) The summary of enrichment analysis in PaGenBase.
(a)(b)


### 3.3.1. PPI Network for All Targets

We upload the names of all target genes to STRING. According to network statistics, the number of nodes is 214, the number of edges is 3057, and the average node degree is 28.6. The expected number of edges is 1173, and the local clustering coefficient is 0.583. We discovered that the network had far more interactions than expected. This suggests that the target proteins as a group are at least partially biologically connected.Using Cytoscape 3.9.1, we constructed a PPI network (Figure3(a)). Then, using the Cytoscape plug-in “cytohubba,” we analyzed hub targets and chose the top 25 target genes as hub genes (Figure 3(b)). CASP3, AKT1, Jun, STAT3, TP53, MMP9, BCL2l1, SRC, and other proteins. The higher the rank, the more important these target genes are in disease treatment. Hub targets are painted red and located at the center of the network for further analysis and research.Figure 3
PPI network of all target genes. (a) PPI network, colored and in the middle are 25 hub genes. (b) Top 25 genes in the network ranked by the MCC method in “Cytohubba”.
(a)(b)We also used the Cytoscape plug-in “MCODE” to examine the PPI network clusters and modules of all target genes (Figure4). The PPI network is divided into six clusters, with 25 hub target genes located in Cluster 1, indicating that hub genes have biological function relevance and may play a synergistic role.Figure 4
Clusters 1–6 in the PPI network. Among them, 25 hub genes are painted red and orange.


## 3.3.1. PPI Network for All Targets

We upload the names of all target genes to STRING. According to network statistics, the number of nodes is 214, the number of edges is 3057, and the average node degree is 28.6. The expected number of edges is 1173, and the local clustering coefficient is 0.583. We discovered that the network had far more interactions than expected. This suggests that the target proteins as a group are at least partially biologically connected.Using Cytoscape 3.9.1, we constructed a PPI network (Figure3(a)). Then, using the Cytoscape plug-in “cytohubba,” we analyzed hub targets and chose the top 25 target genes as hub genes (Figure 3(b)). CASP3, AKT1, Jun, STAT3, TP53, MMP9, BCL2l1, SRC, and other proteins. The higher the rank, the more important these target genes are in disease treatment. Hub targets are painted red and located at the center of the network for further analysis and research.Figure 3
PPI network of all target genes. (a) PPI network, colored and in the middle are 25 hub genes. (b) Top 25 genes in the network ranked by the MCC method in “Cytohubba”.
(a)(b)We also used the Cytoscape plug-in “MCODE” to examine the PPI network clusters and modules of all target genes (Figure4). The PPI network is divided into six clusters, with 25 hub target genes located in Cluster 1, indicating that hub genes have biological function relevance and may play a synergistic role.Figure 4
Clusters 1–6 in the PPI network. Among them, 25 hub genes are painted red and orange.


## 3.4. Herbal-Key Bioactive Ingredient-Hub Target Network

After obtaining the hub target genes, we analyzed the active ingredients corresponding to these 25 hub genes, which are named as key bioactive ingredients. For more information, see Table S2. The network of herbal-key bioactive ingredient-hub targets was constructed using Cytoscape 3.9.1 (Figure5). In addition to Perilla frutescens, the other seven Chinese herbal medicines have three or more corresponding key bioactive ingredients. Some hub genes are affected by multiple bioactive ingredients at the same time. The primary targets of the active ingredients are MAPK14, HSP90AA1, PTGS2, and ESR1. These genes may be the primary targets of the formula.Figure 5
Herbal-key bioactive ingredient-hub target network.


## 3.5. GO Functional Annotation and KEGG Pathway Analysis

To investigate the biological processes engaged in hub targets, GO enrichment analysis and KEGG enrichment analysis on 25 hub genes were analyzed in the David website. The mechanism of action of the formula can be researched, based on the biological process regulated by the hub target.Beautify the enrichment analysis results withR (Figure 6). In total, 226 GO biological process enrichment results were obtained. Negative regulation of the apoptotic process, positive regulation of the nitric oxide biosynthetic process, and positive regulation of transcription from the RNA polymerase II promoter are the top three enrichment biological processes. As shown in Figure 6(a), the top 20 GO biological processes are represented in the form of a bubble diagram, where the size of the circle represents the enrichment of relevant targets in the pathway, and the darker the color of the circle represents the degree of enrichment of targets, indicating that the formula could have physiological effects by regulating these biological processes.Figure 6
GO and KEGG enrichment analysis of hub genes.
(a)(b)(c)(d)For KEGG pathway enrichment analysis, 25 hub targets were mapped into the David database. The species was defined as “human,” and a total of 94 pathways were obtained. As shown in Figure6(d), the top 20 pathways with high significance of KEGG enrichment results are closely related to the mechanism of the Pan Ji Sheng formula. The top five pathways include hepatitis B, pathways in cancer, TNF signaling pathway, toxoplasmosis, and toll-like receptor signaling pathway. The majority of these pathways are linked to the genes TP53, JUN, AKT1, MAPK14, HSP90AA1, and PTGS2.We also performed disease enrichment analysis to investigate diseases associated with hub targets. Figure7 shows the classification of diseases enriched in hub targets. The three major categories are cancer, infection, and immune system. Our findings indicate that the formula studied in this study may primarily target these diseases.Figure 7
Disease and disease class enrichment analysis of hub genes.
(a)(b)


### 3.5.1. Hub Target-GO BP/Pathway/Disease Class Network

In order to demonstrate the biological process of the hub target and the relationship between the hub target and the pathway more clearly, the hub target-GO BP/pathway/disease class network was built with Cytoscape 3.9.1 software (Figure8).Figure 8
Hub target-GO BP/pathway network.The hub target is represented by the circle in the center of Figure8. The left and right sides of Figure 8 show the top 20 enriched biological processes and pathways, respectively. We can clearly understand the relationship between the targets and biological processes or pathways. MAPK14, hSP90AA1, and PTGS2 genes are associated with apoptotic biological processes, TNF signaling pathways, toll-like receptor signaling pathways, and cancer pathways. The formula could play a significant role by regulating these pathways.In order to demonstrate the link between the hub targets and diseases more clearly, Cytoscape 3.9.1 software was used to create a network of hub targets and diseases (Figure9). The genes MAPK14, HSP90AA1, PTGS2, and ESR1 have been linked to cancer, infection, and immune disease.Figure 9
Hub target-disease class network.


## 3.5.1. Hub Target-GO BP/Pathway/Disease Class Network

In order to demonstrate the biological process of the hub target and the relationship between the hub target and the pathway more clearly, the hub target-GO BP/pathway/disease class network was built with Cytoscape 3.9.1 software (Figure8).Figure 8
Hub target-GO BP/pathway network.The hub target is represented by the circle in the center of Figure8. The left and right sides of Figure 8 show the top 20 enriched biological processes and pathways, respectively. We can clearly understand the relationship between the targets and biological processes or pathways. MAPK14, hSP90AA1, and PTGS2 genes are associated with apoptotic biological processes, TNF signaling pathways, toll-like receptor signaling pathways, and cancer pathways. The formula could play a significant role by regulating these pathways.In order to demonstrate the link between the hub targets and diseases more clearly, Cytoscape 3.9.1 software was used to create a network of hub targets and diseases (Figure9). The genes MAPK14, HSP90AA1, PTGS2, and ESR1 have been linked to cancer, infection, and immune disease.Figure 9
Hub target-disease class network.


## 4. Discussion

Traditional Chinese medicine formulas are typically difficult to decipher due to the action mode of traditional Chinese medicine formulas [33]. Using network pharmacology, this study explains the action mechanism of the Pan Ji Sheng Chinese medicine formula. According to the findings of this study, CASP3, AKT1, JUN, and other genes are the hub targets of the formula to enhance immune response and anti-inflammatory.According to the active ingredient-target network, HSP90AA1, PTGS2, ESR1, and MAPK14 are the four key genes regulated by the active ingredient of the Pan Ji Sheng formula. HSP90AA1 is an inflammation-related protein that can be significantly upregulated with some inflammation-related genes in the inflammatory response [34, 35]; PTGS2 is involved in inflammation, immunity, and other processes [36, 37]; ESR1 is also involved in inflammation and immunity and is one of the key targets for the treatment of pneumonia [38, 39]; and MAPK14 is related to autophagy and plays an important role in immune response [40].As shown in the results, 19 of the 25 hub targets were discovered to be involved in the pathways in cancer, with the pathways in cancer being the most significant pathway. This could be due to the fact that respiratory inflammation and lung disease are risk factors for cancer [41, 42]. Other top KEGG enrichment pathways include hepatitis B, the TNF signaling pathway, toxoplasmosis, and the toll-like receptor signaling pathway. A key target gene is tumor necrosis factor (TNF), a cytokine secreted by macrophages and adipocytes. It can cause IR by suppressing the activity of the PI3K/Akt signaling pathway. TNF has been shown to activate MAPK and NF-B signaling pathways, which regulate inflammatory response, oxidative stress, and apoptosis [43, 44].The network pharmacological analysis reveals that the Pan Ji Sheng formula could regulate HSP90AA1, PTGS2, ESR1, MAPK14, and other genes, modulating pathways such as cancer pathways, TNF signaling pathways, and toll-like receptor signaling pathways to regulate inflammatory response and immune processes.This study investigated the anti-inflammatory and immune mechanisms of Pan Ji Sheng formula. However, in vivo and in vitro experiments are needed to provide more information on the mechanism of action of the formula.


## 5. Conclusions

The active components of the Pan Ji Sheng formula could regulate certain proteins, including HSP90AA1, PTGS2, ESR1, and MAPK14. The Chinese herbs in the Pan Ji Sheng formula have a synergistic therapeutic effect, primarily by acting on inflammation and immune-related signal pathways. Pan Ji Sheng formula plays the functions through multicomponents, multitargets (HSP90AA1, PTGS2, ESR1, MAPK14, and other hub targets), and multipathways (inflammation and immune-related signal pathways). These findings could serve as guidelines for future research into this formula. Based on the present study, functional experiments can be performed on animal models or human cells to validate the pharmacological mechanisms of the Pan Ji Sheng formula in the future. This research has theoretical significance for the TCM pharmacology and has application value for the development and utilization of TCMs.



---

*Source: 5340933-2022-09-28.xml*