# The Pharmacological Mechanism of the Effect of Plant Extract Compound Drugs on Cancer Pain Based on Network Pharmacology

Authors: Yuanyuan Shen; Jun Wang; Pengpeng Yan; Tiantian Chen; Xingrui Li; Ming Jiang
Journal: Journal of Healthcare Engineering
 (2022)

---

## Abstract

Objective. We systematically analyzed the mechanism of plant-derived drugs alleviating cancer pain in our hospital through network pharmacology, so as to provide the possibility of further application of traditional Chinese medicine in the treatment of cancer pain. Methods. We used TCMSP, ETCM, and TCMID databases to mine the active ingredients of plant-derived drugs. We combined OMIM, GeneCards, and DrugBank databases to mine and match the common targets of plant-derived drugs for cancer pain. We used the STRING platform and Cytoscape software to analyze and screen out the core targets. We used GO and KEGG methods to analyze the biological processes, molecular functions, cellular composition, and signaling pathways involved in the reduction of cancer pain by plant-derived drugs. Results. We found 153 active ingredients from botanical drugs by TCMSP (Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform, TCMSP), ETCM (The Encyclopedia of Traditional Chinese Medicine), and TCMID (Traditional Chinese Medicine Integrated Database) databases, covering 341 protein targets in human body. Combined with OMIM (Online Mendelian Inheritance in Man), GeneCards, and DrugBank databases, we excavated and matched 141 targets of plant-derived drugs and cancerous pain diseases. Through the analysis of the STRING platform and Cytoscape software, 19 core targets including TNF, MAPK1, JUN, and IL-6 were screened out. Go and KEGG enrichment showed that plant-derived drugs alleviated cancer pain processes involving 193 biological processes, 47 molecular functions, 22 cell components, and 118 signaling pathways. By screening genes involved in KEGG signaling pathway, it was found that plant-derived drugs were mainly associated with PI3K-Akt signaling pathway, tumor necrosis factor signaling pathway, MAPK signaling pathway, Toll-like receptor signaling pathway, and HIF-1 signaling pathway in alleviating cancer pain. Conclusion. These results indicate that botanical drugs can positively affect the expression of inflammatory factors and apoptotic factors in the process of treatment and relief of cancer pain, which is expected to have a potential therapeutic effect on the relief of cancer pain.


---

## Body

## 1. Introduction

Cancer pain is a chronic pain with a unique and complex mechanism. At the early stage, cancer pain is mainly sensitized by primary afferent nerves caused by a series of pain-causing substances released by tumor cells and inflammatory cells and continuously activated osteoclasts, while at the later stage, cancer pain is mainly caused by nerve compression and injury caused by tumor growth [1]. In recent years, with the aging of the society and the increasing incidence of cancer, about half of the patients with malignant tumors have different degrees of cancer pain, especially the incidence of cancer pain in advanced patients is as high as 60–80%, which seriously affects the patients’ disease recovery and quality of life [2]. Cancerous pain is a kind of pain with areas in need of repair or adjustment after the information transmission to the nerve center of feeling, although the clinical use of “three steps” analgesia method can partly control or relieve patients’ pain, but as the disease progresses, the tumor load increases, and the increase of the analgesic drug resistance and adverse drug reactions are led to the decrease of the tolerance [3–5].The understanding of pain in Huangdi Neijing is as follows: “The meridians are not only popular, but also endless around; the cold enters the meridians but delays, and the weeping fails; the blood is little outside the veins, and the qi is blocked in the middle, so the death is painful.” Chinese medicine research believes that cancer pain is cancer toxin cohesion, meridian block blockage, and not general pain. If the human body is in a state of qi deficiency, it can not nourish the viscera and meridians, so pain will occur. Pain and its treatment are the diseases that are first felt and recognized by human beings. In other words, after the body is stimulated by the internal and external environment, pain-causing substances are generated and released from the tissues, pain receptors are sensitized, pain information is transmitted, sensory center perceives, and finally, pain enters the consciousness stage, which leads to pain [6]. Among them, pain-causing substances include bradykinin, histamine, 5-HT, and prostaglandin, which can excite or sensitize the receptors when they are stimulated by various stimuli [7]. At present, western medicine is widely used in the treatment of various kinds of pain, which has the effect of addiction and tolerance, which limits the clinical application and reduces the therapeutic effect [8, 9]. The application of traditional Chinese medicine in the treatment of pain has a history of thousands of years. In ancient books and modern scientific research, it has been proved that more than 800 kinds of traditional Chinese medicine can effectively relieve pain [10, 11]. Traditional Chinese medicine compound has the characteristics of multiherbs, multicomponents, and multitarget coordination, but the current research methods cannot fully explain its effects. Network pharmacology, from a macroscopic and systematic perspective, explores the overall effect of drug treatment on diseases, breaking the previous research mode of single compound and single target, and providing a new idea for the study of the TCM compound. There are 11 herbal medicines of plant origin, including Astragalus membranaceus, Poria cocos, honeysuckle, Cyperus, Paeonia lactiflora, bupleurum, Hedyotis diffusa, Scutellaria barbata, Trichosanthes kirilowii, Arisaema, and Scutellaria baicalensis. Among them, Astragalus membranaceus, Poria cocos, Cyperus, Paeonia lactiflora, and Scutellaria baicalensis have been used many times in the clinical research of pain treatment . However, the specific mechanism of pain relief by botanical drugs is still unclear. In order to optimize the clinical therapeutic effect of botanical drugs and clarify its mechanism of action, we obtained the effective chemical components and corresponding targets through TCMSP, ETCM, TCMID, and other databases; screened the pain-related targets from GeneCards, OMIM, DrugBank, and other databases; and constructed the “humanized” model by using Cytoscape software and STRING database. To screen the key components and key targets of therapeutic effect, GO and KEGG pathway enrichment analysis was carried out on the key targets through the DAVID platform, and the potential mechanism of 11 herbal extracts was obtained [12–18].


## 2. Materials and Methods

### 2.1. Identification and ADME Screening of Candidate Compounds of Traditional Chinese Medicine

In this study, all the active ingredients of 11 herbal medicines of botanical origin were searched by TCMSP, ETCM, and TCMID [19, 20]. The names of 11 traditional Chinese medicines, Astragalus membranaceus, Poria cocos, Honeysuicerae lonicerae, Radix Paeoniae lactiflora, Radix Bupleurum, Hedyotis diffusa, Scutellaria scleroides, Fructus Trichosanthis, Araceae aratidae, Radix Scutellariae, were input, and the main active components were further screened according to the oral availability (OB) ≥30% and drug-like properties (DL) ≥0.18. OB and DL are the key indicators to evaluate the effective availability of drugs. Generally speaking, active ingredients with OB ≥ 30% and DL ≥ 0.18 can be regarded as the main active ingredients of drugs [21].


### 2.2. Prediction of Related Targets of Active Compounds in Traditional Chinese Medicine

TCMSP, PharmMapper, SwissTarget, and other databases were used to screen the protein targets corresponding to the main active components of the abovementioned traditional Chinese medicine [22]. By using the retrieval function of UniproKBT in the UniProt database [23], all target gene names were corrected into official gene symbol names by inputting protein names, and active ingredients without targets were removed to obtain the information of active ingredients and related targets.


### 2.3. Cancer Pain Disease Targets and Protein Collection

According to DrugBank, GeneCards, and OMIM databases, the English “cancer pain” corresponding to pain was used as the key word to search and screen the related targets, and the targets of three disease databases were combined to remove the repetitive base, so as to obtain the potential targets of pain, and by using the retrieval function of UniProKBT in the UniProt database and by inputting the name of protein and limiting the species to human, all the potential targets of pain were identified. The name of target gene was corrected to official gene symbol, and the active components without target were eliminated to obtain the final disease target [24, 25].


### 2.4. Construction of “Compound-Target-Cancerous Pain” Signaling Pathway

The protein-protein interaction (PPI) network of potential targets and disease targets was constructed by using Cytoscape 3.7.2, and the two network graphs were fused and the intersection network was extracted by using the correlation function in the software. The direct and indirect target regulatory networks of 14 herbal medicines for cancer pain relief were obtained.


### 2.5. GO and KEGG Analysis

In order to clarify the potential action targets of active ingredients in 14 herbal medicines of plant origin, as well as the role of cancer-related pain targets in gene function and signaling pathway, GO and KEGG enrichment analysis was conducted in the DAVID database and Metascape platform based on annotation and visualization modules [26]. In this study, R language programming was used for GO and KEGG pathway analysis, and set thresholds P < 0.05. The possible mechanism of analgesic effects of 14 herbal herbs derived from plants was predicted through gene enrichment analysis, and the results were imported into the mapping software for visual mapping [27] (Table 1).Table 1
Website of database and network platform.
TCMSPhttps://tcmspw.com/tcmsp.phpETCMhttp://www.tcmip.cn/ETCM/index.php/TCMIDhttp://www.Megabionet.org/tcmid/PharmMapperhttp://lilab-ecust.cn/pharmmapper/SwissTargethttp://new.swisstargetprediction.ch/UniProthttps://www.uniprot.org/DrugBankhttps://omim.org/GeneCardshttps://omim.org/OMIMhttps://omim.org/DAVIDhttps://metascape.org/gp/Metascapehttps://metascape.org/gp/STRINGhttps://string-db.org/


## 2.1. Identification and ADME Screening of Candidate Compounds of Traditional Chinese Medicine

In this study, all the active ingredients of 11 herbal medicines of botanical origin were searched by TCMSP, ETCM, and TCMID [19, 20]. The names of 11 traditional Chinese medicines, Astragalus membranaceus, Poria cocos, Honeysuicerae lonicerae, Radix Paeoniae lactiflora, Radix Bupleurum, Hedyotis diffusa, Scutellaria scleroides, Fructus Trichosanthis, Araceae aratidae, Radix Scutellariae, were input, and the main active components were further screened according to the oral availability (OB) ≥30% and drug-like properties (DL) ≥0.18. OB and DL are the key indicators to evaluate the effective availability of drugs. Generally speaking, active ingredients with OB ≥ 30% and DL ≥ 0.18 can be regarded as the main active ingredients of drugs [21].


## 2.2. Prediction of Related Targets of Active Compounds in Traditional Chinese Medicine

TCMSP, PharmMapper, SwissTarget, and other databases were used to screen the protein targets corresponding to the main active components of the abovementioned traditional Chinese medicine [22]. By using the retrieval function of UniproKBT in the UniProt database [23], all target gene names were corrected into official gene symbol names by inputting protein names, and active ingredients without targets were removed to obtain the information of active ingredients and related targets.


## 2.3. Cancer Pain Disease Targets and Protein Collection

According to DrugBank, GeneCards, and OMIM databases, the English “cancer pain” corresponding to pain was used as the key word to search and screen the related targets, and the targets of three disease databases were combined to remove the repetitive base, so as to obtain the potential targets of pain, and by using the retrieval function of UniProKBT in the UniProt database and by inputting the name of protein and limiting the species to human, all the potential targets of pain were identified. The name of target gene was corrected to official gene symbol, and the active components without target were eliminated to obtain the final disease target [24, 25].


## 2.4. Construction of “Compound-Target-Cancerous Pain” Signaling Pathway

The protein-protein interaction (PPI) network of potential targets and disease targets was constructed by using Cytoscape 3.7.2, and the two network graphs were fused and the intersection network was extracted by using the correlation function in the software. The direct and indirect target regulatory networks of 14 herbal medicines for cancer pain relief were obtained.


## 2.5. GO and KEGG Analysis

In order to clarify the potential action targets of active ingredients in 14 herbal medicines of plant origin, as well as the role of cancer-related pain targets in gene function and signaling pathway, GO and KEGG enrichment analysis was conducted in the DAVID database and Metascape platform based on annotation and visualization modules [26]. In this study, R language programming was used for GO and KEGG pathway analysis, and set thresholds P < 0.05. The possible mechanism of analgesic effects of 14 herbal herbs derived from plants was predicted through gene enrichment analysis, and the results were imported into the mapping software for visual mapping [27] (Table 1).Table 1
Website of database and network platform.
TCMSPhttps://tcmspw.com/tcmsp.phpETCMhttp://www.tcmip.cn/ETCM/index.php/TCMIDhttp://www.Megabionet.org/tcmid/PharmMapperhttp://lilab-ecust.cn/pharmmapper/SwissTargethttp://new.swisstargetprediction.ch/UniProthttps://www.uniprot.org/DrugBankhttps://omim.org/GeneCardshttps://omim.org/OMIMhttps://omim.org/DAVIDhttps://metascape.org/gp/Metascapehttps://metascape.org/gp/STRINGhttps://string-db.org/


## 3. Results

### 3.1. Analysis of Components and ADME of Botanical Drugs

By searching TCMSP, ETCM, and TCMID databases and analyzing ADME, based on OB ≥ 30% and DL ≥ 0.18, a total of 154 active components were obtained from 11 kinds of herbal medicines, including 23 Astragalus membranaceus, 20 Poria cocos, 26 Honeysuckle, 18 Cyperus, 13 Paeonia alba, 17 Bupleurum, 7 Hedyotis diffusa, 29 Scutellaria barbata, 11 Trichosanthes kirilowii, 7 Arisaema, and 36 Scutellaria baicalensis. The compound information of 11 Traditional Chinese medicines is shown in Figure1 and Table 2–12 [28].Figure 1
Network of botanical drugs, traditional Chinese medicine, and compounds.Table 2
Related information of active components of Astragalus membranaceus.
Mol IDMolecule nameOB (%)DLMOL000211Mairin55.380.78MOL000239Jaranol50.830.29MOL000296Hederagenin36.910.75MOL00003324-Propylcholesterol36.230.78MOL000354Isorhamnetin49.600.31MOL0003713,9-Di-O-methylnissolin53.740.48MOL0003745′-Hydroxyiso-muronulatol-2′,5′-di-O-glucoside41.720.69MOL0003787-O-Methylisomucronulatol74.690.30MOL0003799,10-Dimethoxypterocarpan-3-O-β-D-glucoside36.740.92MOL000380(6aR,11aR)-9,10-Dimethoxy-6a,11a-dihydro-6H-benzofurano[ [2, 3]-c]chromen-3-ol64.260.42MOL000387Bifendate31.100.67MOL000392Formononetin69.670.21MOL000398Isoflavanone109.990.30MOL000417Calycosin47.750.24MOL000433FA68.960.71MOL000438(3R)-3-(2-Hydroxy-3,4-dimethoxyphenyl)chroman-7-ol67.670.26MOL000439Isomucronulatol-7,2′-di-O-glucosiole49.280.62MOL0004421,7-Dihydroxy-3,9-dimethoxy pterocarpene39.050.48MOL000098Quercetin46.430.28MOL000359Sitosterol36.910.75MOL000433Folic acid68.960.71MOL000239Kumatakenin50.830.29MOL002565Medicarpin49.220.34Table 3
Related information of active components of Poria cocos.
Mol IDMolecule nameOB (%)DLMOL000273(2R)-2-[(3S,5R,10S,13R,14R,16R,17R)-3,16-Dihydroxy-4,4,10,13,14-pentamethyl-2,3,5,6,12,15,16,17-octahydro-1H-cyclopenta[a]phenanthren-17-yl]-6-methylhept-5-enoic acid30.930.81MOL000275Trametenolic acid38.710.80MOL0002767,9(11)-Dehydropachymic acid35.110.81MOL000279Cerevisterol37.960.77MOL000280Dehydrotumulosic acid31.070.82MOL000282Ergosta-7,22E-dien-3beta-ol43.510.72MOL000283Ergosterol peroxide40.360.81MOL000285Polyporenic acid C38.260.82MOL0002873β-Hydroxy-24-methylene-8-lanostene-21-oic acid38.700.81MOL000289Pachymic acid33.630.81MOL000290Poricoic acid A30.610.76MOL000291Poricoic acid B30.520.75MOL000292Poricoic acid C38.150.75MOL000296Hederagenin36.910.75MOL000300Dehydroeburicoic acid44.170.83MOL000289Pachymic acid33.630.81MOL000290Poricoic acid A30.610.76MOL000291Poricoic acid B30.520.75MOL000292Poricoic acid C38.150.75MOL000275Trametenolic acid38.710.80Table 4
Related information of active components of honeysuckle.
Mol IDMolecule nameOB (%)DLMOL001494Mandenol42.000.19MOL001495Ethyl linolenate46.100.20MOL002707Phytofluene43.180.50MOL002914Eriodyctiol (flavanone)41.350.24MOL003006(-)-(3R,8S,9R,9aS,10aS)-9-Ethenyl-8-(beta-D-glucopyranosyloxy)-2,3,9,9a,10,10a-hexahydro-5-oxo-5H,8H-pyrano[ [3, 4]-d]oxazolo[ [2, 3]-a]pyridine-3-carboxylic acid_qt87.470.23MOL003014Secologanic dibutylacetal_qt53.650.29MOL002773Beta-carotene37.180.58MOL003036ZINC0397878143.830.76MOL003044Chryseriol35.850.27MOL003059Kryptoxanthin47.250.57MOL0030624,5′-Retro-.beta.,.beta.-Carotene-3,3′-dione, 4′,5′-didehydro-31.220.55MOL0030955-Hydroxy-7-methoxy-2-(3,4,5-trimethoxyphenyl)chromone51.960.41MOL0031017-Epi-Vogeloside46.130.58MOL003108Caeruloside C55.640.73MOL003111Centauroside_qt55.790.50MOL003117Ioniceracetalides B_qt61.190.19MOL003124Xylostosidine43.170.64MOL003128Dinethylsecologanoside48.460.48MOL000358Beta-sitosterol36.910.75MOL000422Kaempferol41.880.24MOL000449Stigmasterol43.830.76MOL000006Luteolin36.160.25MOL000098Quercetin46.430.28MOL000359Sitosterol36.910.75MOL000296Hederagenin36.910.75MOL003044Chrysoeriol35.850.27Table 5
Related information of active components of cyperus.
Mol IDMolecule nameOB (%)DLMOL003044Chryseriol35.850.27MOL000354Isorhamnetin49.600.31MOL0035428-Isopentenyl-kaempferol38.040.39MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL0040271,4-Epoxy-16-hydroxyheneicos-1,3,12,14,18-pentaene45.100.24MOL004053Isodalbergin35.450.20MOL004058Khell33.190.19MOL004059Khellol glucoside74.960.72MOL010489Resivit30.840.27MOL004068Rosenonolactone79.840.37MOL004071Hyndarin73.940.64MOL004074Stigmasterol glucoside_qt43.830.76MOL004077Sugeonyl acetate45.080.20MOL000422Kaempferol41.880.24MOL000449Stigmasterol43.830.76MOL000006Luteolin36.160.25MOL000098Quercetin46.430.28Table 6
Related information of active components ofPaeonia lactiflora.
Mol IDMolecule nameOB (%)DLMOL00191011-Alpha,12-alpha-epoxy-3beta-23-dihydroxy-30-norolean-20-en-28,12beta-olide64.770.38MOL001918Paeoniflorgenone87.590.37MOL001919(3S,5R,8R,9R,10S,14S)-3,17-Dihydroxy-4,4,8,10,14-pentamethyl-2,3,5,6,7,9-hexahydro-1H-cyclopenta[a]phenanthrene-15,16-dione43.560.53MOL001921Lactiflorin49.120.80MOL001924Paeoniflorin53.870.79MOL001925paeoniflorin_qt68.180.40MOL001928albiflorin_qt66.640.33MOL001930Benzoyl paeoniflorin31.270.75MOL000211Mairin55.380.78MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL000422Kaempferol41.880.24MOL000492(+)-Catechin54.830.24Table 7
Related information of active components of bupleurum.
Mol IDMolecule nameOB (%)DLMOL001645Linoleyl acetate42.100.20MOL002776Baicalin40.120.75MOL000449Stigmasterol43.830.76MOL000354Isorhamnetin49.600.31MOL000422Kaempferol41.880.24MOL0045983,5,6,7-Tetramethoxy-2-(3,4,5-trimethoxyphenyl)chromone31.970.59MOL004609Areapillin48.960.41MOL013187Cubebin57.130.64MOL004624Longikaurin A47.720.53MOL004628Octalupine47.820.28MOL004644Sainfuran79.910.23MOL004648Troxerutin31.600.28MOL004653(+)-Anomalin46.060.66MOL004702Saikosaponin c_qt30.500.63MOL004718α-Spinasterol42.980.76MOL000490Petunidin30.050.31MOL000098Quercetin46.430.28Table 8
Related information of active components ofHedyotis diffusa.
Mol IDMolecule nameOB (%)DLMOL0016462,3-Dimethoxy-6-methyanthraquinone34.860.26MOL001659Poriferasterol43.830.76MOL001663Oleanolic acid32.030.76MOL0016702-Methoxy-3-methyl-9,10-anthraquinone37.830.21MOL000449Stigmasterol43.830.76MOL000358Beta-sitosterol36.910.75MOL000098Quercetin46.430.28Table 9
Related information of active components of Scutellaria barbata.
Mol IDMolecule nameOB (%)DLMOL001040(2R)-5,7-Dihydroxy-2-(4-hydroxyphenyl)chroman-4-one42.360.21MOL0122455,7,4′-Trihydroxy-6-methoxyflavanone36.630.27MOL0122465,7,4′-Trihydroxy-8-methoxyflavanone74.240.26MOL0122485-Hydroxy-7,8-dimethoxy-2-(4-methoxyphenyl)chromone65.820.33MOL0122507-Hydroxy-5,8-dimethoxy-2-phenyl-chromone43.720.25MOL012251Chrysin-5-methylether37.270.20MOL0122529,19-Cyclolanost-24-en-3-ol38.690.78MOL002776Baicalin40.120.75MOL012254Campesterol37.580.71MOL000953CLR37.870.68MOL000358Beta-sitosterol36.910.75MOL012266Rivularin37.940.37MOL001973Sitosteryl acetate40.390.85MOL012269Stigmasta-5,22-dien-3-ol-acetate46.440.86MOL012270Stigmastan-3,5,22-triene45.030.71MOL000449Stigmasterol43.830.76MOL000173Wogonin30.680.23MOL001735Dinatin30.970.27MOL00175524-Ethylcholest-4-en-3-one36.080.76MOL002714Baicalein33.520.21MOL0027196-Hydroxynaringenin33.230.24MOL002915Salvigenin49.070.33MOL000351Rhamnazin47.140.34MOL000359Sitosterol36.910.75MOL005190Eriodictyol71.790.24MOL005869Daucostero_qt36.910.75MOL000006Luteolin36.160.25MOL008206Moslosooflavone44.090.25MOL000098Quercetin46.430.28Table 10
Related information of active components of Trichosanthes kirilowii.
Mol IDMolecule nameOB (%)DLMOL001494Mandenol42.000.19MOL002881Diosmetin31.140.27MOL004355Spinasterol42.980.76MOL005530Hydroxygenkwanin36.470.27MOL006756Schottenol37.420.75MOL00716510α-Cucurbita-5,24-diene-3β-ol44.020.74MOL0071715-Dehydrokarounidiol30.230.77MOL0071727-Oxo-dihydrokaro-unidiol36.850.75MOL007175Karounidiol 3-o-benzoate43.990.50MOL007179Linolenic acid ethyl ester46.100.20MOL007180Vitamin-e32.290.70MOL0131468,11,14-Docosatrienoic acid, methyl ester43.230.30MOL013156[(2R)-2-[[[(2R)-2-(Benzoylamino)-3-phenylpropanoyl]amino]methyl]-3-phenylpropyl] acetate38.880.56MOL00151024-Epicampesterol37.580.71MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL000449Stigmasterol43.830.76MOL000953CLR37.870.68Table 11
Related information of active components of Rhizoma Arisaema.
Mol IDMolecule nameOB (%)DLMOL0131468,11,14-Docosatrienoic acid, methyl ester43.230.30MOL013156[(2R)-2-[[[(2R)-2-(Benzoylamino)-3-phenylpropanoyl]amino]methyl]-3-phenylpropyl] acetate38.880.56MOL00151024-Epicampesterol37.580.71MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL000449Stigmasterol43.830.76MOL000953CLR37.870.68Table 12
Related information of active components ofScutellaria baicalensis Georgi.
Mol IDMolecule nameOB (%)DLMOL001689Acacetin34.970.24MOL000173Wogonin30.680.23MOL000228(2R)-7-Hydroxy-5-methoxy-2-phenylchroman-4-one55.230.20MOL002714Baicalein33.520.21MOL0029085,8,2′-Trihydroxy-7-methoxyflavone37.010.27MOL0029095,7,2,5-Tetrahydroxy-8,6-dimethoxyflavone33.820.45MOL002910Carthamidin41.150.24MOL0029112,6,2′,4′-Tetrahydroxy-6′-methoxychaleone69.040.22MOL002913Dihydrobaicalin_qt40.040.21MOL002914Eriodyctiol (flavanone)41.350.24MOL002915Salvigenin49.070.33MOL0029175,2′,6′-Trihydroxy-7,8-dimethoxyflavone45.050.33MOL0029255,7,2′,6′-Tetrahydroxyflavone37.010.24MOL002926Dihydrooroxylin A38.720.23MOL002927Skullcapflavone II69.510.44MOL002928Oroxylin a41.370.23MOL002932Panicolin76.260.29MOL0029335,7,4′-Trihydroxy-8-methoxyflavone36.560.27MOL002934Neobaicalein104.340.44MOL002937Dihydrooroxylin66.060.23MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL000525Norwogonin39.400.21MOL0005525,2′-Dihydroxy-6,7,8-trimethoxyflavone31.710.35MOL000073Ent-Epicatechin48.960.24MOL000449Stigmasterol43.830.76MOL001458Coptisine30.670.86MOL001490Bis[(2S)-2-ethylhexyl] benzene-1,2-dicarboxylate43.590.35MOL001506Supraene33.550.42MOL002879Diop43.590.39MOL002897Epiberberine43.090.78MOL008206Moslosooflavone44.090.25MOL01041511,13-Eicosadienoic acid, methyl ester39.280.23MOL0122455,7,4′-Trihydroxy-6-methoxyflavanone36.630.27MOL0122465,7,4′-Trihydroxy-8-methoxyflavanone74.240.26MOL012266Rivularin37.940.37


### 3.2. Target Prediction of Active Components in Traditional Chinese Medicine

Through TCMSP, PharmMapper, SwissTarget database retrieval, and comparison with the UniProt database, a total of 341 protein targets of 153 active components of botanical drugs were obtained, including 247 Astragalus, 94 Poria cocos, 262 Honeysuckle, 178 Cyperus, 72 Paeonia lactiflora, 148 Hupleurum, 153 Hedyotis diffusa, 178 Scutellaria barbata, 9 Trichosanthes kirilowii, 39 Arisaema, and 90 Scutellaria baicalensis. The protein target information of 11 active components of traditional Chinese medicine is shown in Figure2.Figure 2
Network of traditional Chinese medicine active ingredients target.


### 3.3. Screening of Pain-Related Targets

Based on OMIM, GeneCards, and DrugBank databases, pain disease-related targets were screened with “cancer pain” as the key word. The results showed that there were 670 pain disease-related targets in the OMIM database, 1086 GeneCards database, and 215 DrugBank database. A total of 1695 targets were obtained from the three databases. Compared with the UniProt database, 1110 verified human pain disease-related targets were obtained. By matching 341 protein targets of active components of plant-derived drugs, 141 common targets of plant-derived drugs and cancer pain diseases were obtained. The Venn diagram and detailed target information are shown in Figure3 and Table 13, respectively.Figure 3
Venn map of target. ((a) is the Venn map of cancer pain-related targets in OMIM, GeneCards, and DrugBank databases; (b) is the Venn map of common targets of botanical drugs and cancer pain diseases).
(a)(b)Table 13
Common targets of botanical drugs and cancer pain diseases.
141 common target proteinsACHECD40LGEGFRHRH1MAPK14NR1I2RAF1ACTBCES1ERBB2HSD3B2MAPK8NR3C1RASA1ADRA1ACHEK2ESR1HSP90AA1MDM2NR3C2RB1ADRA2ACHRM2ESR2HSPB1METODC1RELAADRA2CCHRM3F3HTR3AMMP1OPRD1SCN5AADRB1COMTF7ICAM1MMP2OPRM1SELEADRB2COX5AFASLGIFNGMMP3PCNASERPINE1AKR1C1COX6B1FECHIGF2MMP9PGRSLC6A2AKT1CRPFN1IGFBP3MPOPIK3CGSLC6A3ALOX5CXCL10FOSIKBKBMT-CO1PLATSLC6A4ARCXCL8GABRA2IL10MTHFRPLAUSPP1BAXCYCSGABRA3IL1AMTRPON1STAT1BCL2CYP19A1GRIN1IL1BMYCPPARATHBDBIRC5CYP1A1GRIN2AIL2NFE2L2PPARGTNFCACNA2D1CYP1A2GRIN2BIL4NFKB1PRKCATOP2ACASP8CYP27B1GRIN2DIL6NFKBIAPRSS1TP53CASP9CYP2C9GSTM1IRF1NOS2PTGESTYRCATCYP3A4GSTP1JUNNOS3PTGS1VCAM1CAV1DNMT1HIF1AKDRNPPBPTGS2VDRCCL2EGFHMOX1MAPK1NR0B1PYGMVEGFACCND1


### 3.4. PPI Network Construction

The interaction network of 141 potential targets was constructed by using the STRING database platform. The minimum interaction threshold was set to “0.9,” and the results are shown in Figure4. The nodes in the figure are intersection genes, and the edges represent the association degree of the intersection genes, and the thickness represents the binding degree. As shown in the figure, there are 140 nodes and 538 edges. The PPI network diagram is imported into the software Cytoscape 3.7.2, and the network pictures are analyzed by tools. We selected those targets with greater than average degree and betweenness centrality, a total of 19 targets, as potential key targets for plant-derived drugs to exert analgesic effects The main proteins closely related to analgesic effect of botanical drugs are TNF, mapk1, Jun, and IL-6, as shown in Figure 5 and Table 14.Figure 4
PPI network of botanical drugs cancer pain diseases.Figure 5
Core target map of different parameters of botanical drug cancer pain disease (selectivity and betweenness centrality are (a) 30; (b) 25; (c) 20; (d) 15).
(a)(b)(c)(d)Table 14
Core target map of different parameters of botanical drugs cancer pain disease.
Degree value 30Degree value 25Degree value 20Degree value 15TNFTNFTNFCXCL8NR3C1MAPK1MAPK14MAPK1IL6AKT1MAPK1IL1BNFKB1IL4JUNMAPK1RELAEGFRTNFRELAESR1AKT1JUNMAPK14IL6EGFRHSP90AA1RELAVEGFAESR1MAPK8JUNTP53JUNTP53HSP90AA1CXCL8VEGFAHSP90AA1IL6AKT1FOSTP53


### 3.5. Construction of “Compound Traditional Chinese Medicine Compound Target Disease” Network Diagram of Analgesic Effect of Botanical Drugs

Nineteen core targets of botanical drug analgesia were obtained from 2.4, and the corresponding compounds and traditional Chinese medicine were screened and matched. The network diagram of “compound traditional Chinese medicine compound target disease” of botanical drug analgesia was predicted by using Cytoscape 3.7.2 software. The results are shown in Figures6 and 7.Figure 6
Traditional Chinese medicine target disease network.Figure 7
Compound Chinese medicine compound target disease network.


### 3.6. GO and KEGG Enrichment Analysis

Through DAVID and STRING databases, 141 targets of botanical drug pain disease were enriched and analyzed, including biological process (BP), molecular function (MF), and cell composition (CC). There are 193 biological processes, 47 molecular functions, and 22 cell components; 20 enriched genes are selected and mapped through the mapping software, as shown in Figure8; 118 signaling pathways are involved in KEGG pathway abdominal muscle analysis, and 20 enriched genes are screened according to the number of enriched genes, as shown in Figure 9.Figure 8
Enrichment analysis of GO in pain target of botanical drugs.Figure 9
Enrichment analysis of the KEGG pathway, a pain target of botanical drugs.


## 3.1. Analysis of Components and ADME of Botanical Drugs

By searching TCMSP, ETCM, and TCMID databases and analyzing ADME, based on OB ≥ 30% and DL ≥ 0.18, a total of 154 active components were obtained from 11 kinds of herbal medicines, including 23 Astragalus membranaceus, 20 Poria cocos, 26 Honeysuckle, 18 Cyperus, 13 Paeonia alba, 17 Bupleurum, 7 Hedyotis diffusa, 29 Scutellaria barbata, 11 Trichosanthes kirilowii, 7 Arisaema, and 36 Scutellaria baicalensis. The compound information of 11 Traditional Chinese medicines is shown in Figure1 and Table 2–12 [28].Figure 1
Network of botanical drugs, traditional Chinese medicine, and compounds.Table 2
Related information of active components of Astragalus membranaceus.
Mol IDMolecule nameOB (%)DLMOL000211Mairin55.380.78MOL000239Jaranol50.830.29MOL000296Hederagenin36.910.75MOL00003324-Propylcholesterol36.230.78MOL000354Isorhamnetin49.600.31MOL0003713,9-Di-O-methylnissolin53.740.48MOL0003745′-Hydroxyiso-muronulatol-2′,5′-di-O-glucoside41.720.69MOL0003787-O-Methylisomucronulatol74.690.30MOL0003799,10-Dimethoxypterocarpan-3-O-β-D-glucoside36.740.92MOL000380(6aR,11aR)-9,10-Dimethoxy-6a,11a-dihydro-6H-benzofurano[ [2, 3]-c]chromen-3-ol64.260.42MOL000387Bifendate31.100.67MOL000392Formononetin69.670.21MOL000398Isoflavanone109.990.30MOL000417Calycosin47.750.24MOL000433FA68.960.71MOL000438(3R)-3-(2-Hydroxy-3,4-dimethoxyphenyl)chroman-7-ol67.670.26MOL000439Isomucronulatol-7,2′-di-O-glucosiole49.280.62MOL0004421,7-Dihydroxy-3,9-dimethoxy pterocarpene39.050.48MOL000098Quercetin46.430.28MOL000359Sitosterol36.910.75MOL000433Folic acid68.960.71MOL000239Kumatakenin50.830.29MOL002565Medicarpin49.220.34Table 3
Related information of active components of Poria cocos.
Mol IDMolecule nameOB (%)DLMOL000273(2R)-2-[(3S,5R,10S,13R,14R,16R,17R)-3,16-Dihydroxy-4,4,10,13,14-pentamethyl-2,3,5,6,12,15,16,17-octahydro-1H-cyclopenta[a]phenanthren-17-yl]-6-methylhept-5-enoic acid30.930.81MOL000275Trametenolic acid38.710.80MOL0002767,9(11)-Dehydropachymic acid35.110.81MOL000279Cerevisterol37.960.77MOL000280Dehydrotumulosic acid31.070.82MOL000282Ergosta-7,22E-dien-3beta-ol43.510.72MOL000283Ergosterol peroxide40.360.81MOL000285Polyporenic acid C38.260.82MOL0002873β-Hydroxy-24-methylene-8-lanostene-21-oic acid38.700.81MOL000289Pachymic acid33.630.81MOL000290Poricoic acid A30.610.76MOL000291Poricoic acid B30.520.75MOL000292Poricoic acid C38.150.75MOL000296Hederagenin36.910.75MOL000300Dehydroeburicoic acid44.170.83MOL000289Pachymic acid33.630.81MOL000290Poricoic acid A30.610.76MOL000291Poricoic acid B30.520.75MOL000292Poricoic acid C38.150.75MOL000275Trametenolic acid38.710.80Table 4
Related information of active components of honeysuckle.
Mol IDMolecule nameOB (%)DLMOL001494Mandenol42.000.19MOL001495Ethyl linolenate46.100.20MOL002707Phytofluene43.180.50MOL002914Eriodyctiol (flavanone)41.350.24MOL003006(-)-(3R,8S,9R,9aS,10aS)-9-Ethenyl-8-(beta-D-glucopyranosyloxy)-2,3,9,9a,10,10a-hexahydro-5-oxo-5H,8H-pyrano[ [3, 4]-d]oxazolo[ [2, 3]-a]pyridine-3-carboxylic acid_qt87.470.23MOL003014Secologanic dibutylacetal_qt53.650.29MOL002773Beta-carotene37.180.58MOL003036ZINC0397878143.830.76MOL003044Chryseriol35.850.27MOL003059Kryptoxanthin47.250.57MOL0030624,5′-Retro-.beta.,.beta.-Carotene-3,3′-dione, 4′,5′-didehydro-31.220.55MOL0030955-Hydroxy-7-methoxy-2-(3,4,5-trimethoxyphenyl)chromone51.960.41MOL0031017-Epi-Vogeloside46.130.58MOL003108Caeruloside C55.640.73MOL003111Centauroside_qt55.790.50MOL003117Ioniceracetalides B_qt61.190.19MOL003124Xylostosidine43.170.64MOL003128Dinethylsecologanoside48.460.48MOL000358Beta-sitosterol36.910.75MOL000422Kaempferol41.880.24MOL000449Stigmasterol43.830.76MOL000006Luteolin36.160.25MOL000098Quercetin46.430.28MOL000359Sitosterol36.910.75MOL000296Hederagenin36.910.75MOL003044Chrysoeriol35.850.27Table 5
Related information of active components of cyperus.
Mol IDMolecule nameOB (%)DLMOL003044Chryseriol35.850.27MOL000354Isorhamnetin49.600.31MOL0035428-Isopentenyl-kaempferol38.040.39MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL0040271,4-Epoxy-16-hydroxyheneicos-1,3,12,14,18-pentaene45.100.24MOL004053Isodalbergin35.450.20MOL004058Khell33.190.19MOL004059Khellol glucoside74.960.72MOL010489Resivit30.840.27MOL004068Rosenonolactone79.840.37MOL004071Hyndarin73.940.64MOL004074Stigmasterol glucoside_qt43.830.76MOL004077Sugeonyl acetate45.080.20MOL000422Kaempferol41.880.24MOL000449Stigmasterol43.830.76MOL000006Luteolin36.160.25MOL000098Quercetin46.430.28Table 6
Related information of active components ofPaeonia lactiflora.
Mol IDMolecule nameOB (%)DLMOL00191011-Alpha,12-alpha-epoxy-3beta-23-dihydroxy-30-norolean-20-en-28,12beta-olide64.770.38MOL001918Paeoniflorgenone87.590.37MOL001919(3S,5R,8R,9R,10S,14S)-3,17-Dihydroxy-4,4,8,10,14-pentamethyl-2,3,5,6,7,9-hexahydro-1H-cyclopenta[a]phenanthrene-15,16-dione43.560.53MOL001921Lactiflorin49.120.80MOL001924Paeoniflorin53.870.79MOL001925paeoniflorin_qt68.180.40MOL001928albiflorin_qt66.640.33MOL001930Benzoyl paeoniflorin31.270.75MOL000211Mairin55.380.78MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL000422Kaempferol41.880.24MOL000492(+)-Catechin54.830.24Table 7
Related information of active components of bupleurum.
Mol IDMolecule nameOB (%)DLMOL001645Linoleyl acetate42.100.20MOL002776Baicalin40.120.75MOL000449Stigmasterol43.830.76MOL000354Isorhamnetin49.600.31MOL000422Kaempferol41.880.24MOL0045983,5,6,7-Tetramethoxy-2-(3,4,5-trimethoxyphenyl)chromone31.970.59MOL004609Areapillin48.960.41MOL013187Cubebin57.130.64MOL004624Longikaurin A47.720.53MOL004628Octalupine47.820.28MOL004644Sainfuran79.910.23MOL004648Troxerutin31.600.28MOL004653(+)-Anomalin46.060.66MOL004702Saikosaponin c_qt30.500.63MOL004718α-Spinasterol42.980.76MOL000490Petunidin30.050.31MOL000098Quercetin46.430.28Table 8
Related information of active components ofHedyotis diffusa.
Mol IDMolecule nameOB (%)DLMOL0016462,3-Dimethoxy-6-methyanthraquinone34.860.26MOL001659Poriferasterol43.830.76MOL001663Oleanolic acid32.030.76MOL0016702-Methoxy-3-methyl-9,10-anthraquinone37.830.21MOL000449Stigmasterol43.830.76MOL000358Beta-sitosterol36.910.75MOL000098Quercetin46.430.28Table 9
Related information of active components of Scutellaria barbata.
Mol IDMolecule nameOB (%)DLMOL001040(2R)-5,7-Dihydroxy-2-(4-hydroxyphenyl)chroman-4-one42.360.21MOL0122455,7,4′-Trihydroxy-6-methoxyflavanone36.630.27MOL0122465,7,4′-Trihydroxy-8-methoxyflavanone74.240.26MOL0122485-Hydroxy-7,8-dimethoxy-2-(4-methoxyphenyl)chromone65.820.33MOL0122507-Hydroxy-5,8-dimethoxy-2-phenyl-chromone43.720.25MOL012251Chrysin-5-methylether37.270.20MOL0122529,19-Cyclolanost-24-en-3-ol38.690.78MOL002776Baicalin40.120.75MOL012254Campesterol37.580.71MOL000953CLR37.870.68MOL000358Beta-sitosterol36.910.75MOL012266Rivularin37.940.37MOL001973Sitosteryl acetate40.390.85MOL012269Stigmasta-5,22-dien-3-ol-acetate46.440.86MOL012270Stigmastan-3,5,22-triene45.030.71MOL000449Stigmasterol43.830.76MOL000173Wogonin30.680.23MOL001735Dinatin30.970.27MOL00175524-Ethylcholest-4-en-3-one36.080.76MOL002714Baicalein33.520.21MOL0027196-Hydroxynaringenin33.230.24MOL002915Salvigenin49.070.33MOL000351Rhamnazin47.140.34MOL000359Sitosterol36.910.75MOL005190Eriodictyol71.790.24MOL005869Daucostero_qt36.910.75MOL000006Luteolin36.160.25MOL008206Moslosooflavone44.090.25MOL000098Quercetin46.430.28Table 10
Related information of active components of Trichosanthes kirilowii.
Mol IDMolecule nameOB (%)DLMOL001494Mandenol42.000.19MOL002881Diosmetin31.140.27MOL004355Spinasterol42.980.76MOL005530Hydroxygenkwanin36.470.27MOL006756Schottenol37.420.75MOL00716510α-Cucurbita-5,24-diene-3β-ol44.020.74MOL0071715-Dehydrokarounidiol30.230.77MOL0071727-Oxo-dihydrokaro-unidiol36.850.75MOL007175Karounidiol 3-o-benzoate43.990.50MOL007179Linolenic acid ethyl ester46.100.20MOL007180Vitamin-e32.290.70MOL0131468,11,14-Docosatrienoic acid, methyl ester43.230.30MOL013156[(2R)-2-[[[(2R)-2-(Benzoylamino)-3-phenylpropanoyl]amino]methyl]-3-phenylpropyl] acetate38.880.56MOL00151024-Epicampesterol37.580.71MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL000449Stigmasterol43.830.76MOL000953CLR37.870.68Table 11
Related information of active components of Rhizoma Arisaema.
Mol IDMolecule nameOB (%)DLMOL0131468,11,14-Docosatrienoic acid, methyl ester43.230.30MOL013156[(2R)-2-[[[(2R)-2-(Benzoylamino)-3-phenylpropanoyl]amino]methyl]-3-phenylpropyl] acetate38.880.56MOL00151024-Epicampesterol37.580.71MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL000449Stigmasterol43.830.76MOL000953CLR37.870.68Table 12
Related information of active components ofScutellaria baicalensis Georgi.
Mol IDMolecule nameOB (%)DLMOL001689Acacetin34.970.24MOL000173Wogonin30.680.23MOL000228(2R)-7-Hydroxy-5-methoxy-2-phenylchroman-4-one55.230.20MOL002714Baicalein33.520.21MOL0029085,8,2′-Trihydroxy-7-methoxyflavone37.010.27MOL0029095,7,2,5-Tetrahydroxy-8,6-dimethoxyflavone33.820.45MOL002910Carthamidin41.150.24MOL0029112,6,2′,4′-Tetrahydroxy-6′-methoxychaleone69.040.22MOL002913Dihydrobaicalin_qt40.040.21MOL002914Eriodyctiol (flavanone)41.350.24MOL002915Salvigenin49.070.33MOL0029175,2′,6′-Trihydroxy-7,8-dimethoxyflavone45.050.33MOL0029255,7,2′,6′-Tetrahydroxyflavone37.010.24MOL002926Dihydrooroxylin A38.720.23MOL002927Skullcapflavone II69.510.44MOL002928Oroxylin a41.370.23MOL002932Panicolin76.260.29MOL0029335,7,4′-Trihydroxy-8-methoxyflavone36.560.27MOL002934Neobaicalein104.340.44MOL002937Dihydrooroxylin66.060.23MOL000358Beta-sitosterol36.910.75MOL000359Sitosterol36.910.75MOL000525Norwogonin39.400.21MOL0005525,2′-Dihydroxy-6,7,8-trimethoxyflavone31.710.35MOL000073Ent-Epicatechin48.960.24MOL000449Stigmasterol43.830.76MOL001458Coptisine30.670.86MOL001490Bis[(2S)-2-ethylhexyl] benzene-1,2-dicarboxylate43.590.35MOL001506Supraene33.550.42MOL002879Diop43.590.39MOL002897Epiberberine43.090.78MOL008206Moslosooflavone44.090.25MOL01041511,13-Eicosadienoic acid, methyl ester39.280.23MOL0122455,7,4′-Trihydroxy-6-methoxyflavanone36.630.27MOL0122465,7,4′-Trihydroxy-8-methoxyflavanone74.240.26MOL012266Rivularin37.940.37


## 3.2. Target Prediction of Active Components in Traditional Chinese Medicine

Through TCMSP, PharmMapper, SwissTarget database retrieval, and comparison with the UniProt database, a total of 341 protein targets of 153 active components of botanical drugs were obtained, including 247 Astragalus, 94 Poria cocos, 262 Honeysuckle, 178 Cyperus, 72 Paeonia lactiflora, 148 Hupleurum, 153 Hedyotis diffusa, 178 Scutellaria barbata, 9 Trichosanthes kirilowii, 39 Arisaema, and 90 Scutellaria baicalensis. The protein target information of 11 active components of traditional Chinese medicine is shown in Figure2.Figure 2
Network of traditional Chinese medicine active ingredients target.


## 3.3. Screening of Pain-Related Targets

Based on OMIM, GeneCards, and DrugBank databases, pain disease-related targets were screened with “cancer pain” as the key word. The results showed that there were 670 pain disease-related targets in the OMIM database, 1086 GeneCards database, and 215 DrugBank database. A total of 1695 targets were obtained from the three databases. Compared with the UniProt database, 1110 verified human pain disease-related targets were obtained. By matching 341 protein targets of active components of plant-derived drugs, 141 common targets of plant-derived drugs and cancer pain diseases were obtained. The Venn diagram and detailed target information are shown in Figure3 and Table 13, respectively.Figure 3
Venn map of target. ((a) is the Venn map of cancer pain-related targets in OMIM, GeneCards, and DrugBank databases; (b) is the Venn map of common targets of botanical drugs and cancer pain diseases).
(a)(b)Table 13
Common targets of botanical drugs and cancer pain diseases.
141 common target proteinsACHECD40LGEGFRHRH1MAPK14NR1I2RAF1ACTBCES1ERBB2HSD3B2MAPK8NR3C1RASA1ADRA1ACHEK2ESR1HSP90AA1MDM2NR3C2RB1ADRA2ACHRM2ESR2HSPB1METODC1RELAADRA2CCHRM3F3HTR3AMMP1OPRD1SCN5AADRB1COMTF7ICAM1MMP2OPRM1SELEADRB2COX5AFASLGIFNGMMP3PCNASERPINE1AKR1C1COX6B1FECHIGF2MMP9PGRSLC6A2AKT1CRPFN1IGFBP3MPOPIK3CGSLC6A3ALOX5CXCL10FOSIKBKBMT-CO1PLATSLC6A4ARCXCL8GABRA2IL10MTHFRPLAUSPP1BAXCYCSGABRA3IL1AMTRPON1STAT1BCL2CYP19A1GRIN1IL1BMYCPPARATHBDBIRC5CYP1A1GRIN2AIL2NFE2L2PPARGTNFCACNA2D1CYP1A2GRIN2BIL4NFKB1PRKCATOP2ACASP8CYP27B1GRIN2DIL6NFKBIAPRSS1TP53CASP9CYP2C9GSTM1IRF1NOS2PTGESTYRCATCYP3A4GSTP1JUNNOS3PTGS1VCAM1CAV1DNMT1HIF1AKDRNPPBPTGS2VDRCCL2EGFHMOX1MAPK1NR0B1PYGMVEGFACCND1


## 3.4. PPI Network Construction

The interaction network of 141 potential targets was constructed by using the STRING database platform. The minimum interaction threshold was set to “0.9,” and the results are shown in Figure4. The nodes in the figure are intersection genes, and the edges represent the association degree of the intersection genes, and the thickness represents the binding degree. As shown in the figure, there are 140 nodes and 538 edges. The PPI network diagram is imported into the software Cytoscape 3.7.2, and the network pictures are analyzed by tools. We selected those targets with greater than average degree and betweenness centrality, a total of 19 targets, as potential key targets for plant-derived drugs to exert analgesic effects The main proteins closely related to analgesic effect of botanical drugs are TNF, mapk1, Jun, and IL-6, as shown in Figure 5 and Table 14.Figure 4
PPI network of botanical drugs cancer pain diseases.Figure 5
Core target map of different parameters of botanical drug cancer pain disease (selectivity and betweenness centrality are (a) 30; (b) 25; (c) 20; (d) 15).
(a)(b)(c)(d)Table 14
Core target map of different parameters of botanical drugs cancer pain disease.
Degree value 30Degree value 25Degree value 20Degree value 15TNFTNFTNFCXCL8NR3C1MAPK1MAPK14MAPK1IL6AKT1MAPK1IL1BNFKB1IL4JUNMAPK1RELAEGFRTNFRELAESR1AKT1JUNMAPK14IL6EGFRHSP90AA1RELAVEGFAESR1MAPK8JUNTP53JUNTP53HSP90AA1CXCL8VEGFAHSP90AA1IL6AKT1FOSTP53


## 3.5. Construction of “Compound Traditional Chinese Medicine Compound Target Disease” Network Diagram of Analgesic Effect of Botanical Drugs

Nineteen core targets of botanical drug analgesia were obtained from 2.4, and the corresponding compounds and traditional Chinese medicine were screened and matched. The network diagram of “compound traditional Chinese medicine compound target disease” of botanical drug analgesia was predicted by using Cytoscape 3.7.2 software. The results are shown in Figures6 and 7.Figure 6
Traditional Chinese medicine target disease network.Figure 7
Compound Chinese medicine compound target disease network.


## 3.6. GO and KEGG Enrichment Analysis

Through DAVID and STRING databases, 141 targets of botanical drug pain disease were enriched and analyzed, including biological process (BP), molecular function (MF), and cell composition (CC). There are 193 biological processes, 47 molecular functions, and 22 cell components; 20 enriched genes are selected and mapped through the mapping software, as shown in Figure8; 118 signaling pathways are involved in KEGG pathway abdominal muscle analysis, and 20 enriched genes are screened according to the number of enriched genes, as shown in Figure 9.Figure 8
Enrichment analysis of GO in pain target of botanical drugs.Figure 9
Enrichment analysis of the KEGG pathway, a pain target of botanical drugs.


## 4. Discussion

Modern medicine believes that cancer pain is mainly caused by the disease itself or the pain caused in the treatment process. Tumors infiltrate bone tissue, and bone tissue destruction results in the release of prostaglandins. It invades the viscera, causing vasospasm, occlusion, and eventually necrosis of the viscera. Or it may compress peripheral nerves, nerve roots, spinal cord, etc. The patient’s own mental tension and psychological pressure can cause pain disorders. There is no name for cancer pain in Chinese medicine, but according to the pain of different parts of cancer, it can be assigned to the pain syndrome of the corresponding parts. For example, pain caused by brain tumor, nasopharyngeal cancer, and cancer brain metastasis is classified as “headache.” Cancer of the esophagus and lung can be classified as “chest pain.” Pain caused by gastric cancer is classified as “stomachache.” Pain caused by liver cancer is classified as “flank pain.” Pain caused by pancreatic and colorectal cancer is classified as “abdominal pain.” The pain caused by bone cancer and cancer bone metastasis is classified as “Bi syndrome” or “Gubi syndrome.”Pain is a kind of feeling of human disease caused by tissue injury, which stimulates the sensory system [29]. It is of great significance and research value in clinical diagnosis and treatment. Usually, the clinical pain mainly comes from the noxious stimulation of internal and external conditions caused by disease or surgery, which not only includes the defense response of the human body to noxious stimulation, but also a clinical manifestation of a variety of diseases and postoperative reactions [30]. Most of the patients have severe pain or chronic long-term pain and need to be treated with analgesic drugs. Cancer pain usually accompanies the whole course of cancer patients, so that patients are in pain, anxiety and even depression, and other negative emotions for a long time, causing some patients to lose confidence in life, which seriously affects the quality of life of patients. Therefore, how to effectively control the degree of cancer pain of cancer patients has great significance for the treatment of cancer and the improvement of life quality of cancer patients [31, 32]. We need to understand the meaning of this. Therefore, the research of analgesic drugs has become one of the hot directions of modern drug research, but the clinical use of chemical drugs for acute or chronic pain, such as opioids and nonopioids: morphine, codeine, pethidine, aspirin, indomethacin, ibuprofen, and so on [33]. Although the above methods have achieved good clinical efficacy, the side effects are often more serious due to the continuous and large dose of use. Traditional Chinese medicine believes that pain is due to body deficiency, deficiency of healthy qi, and invasion of external pathogens, which eventually leads to heat toxin and internal stagnation, obstruction of meridians and collaterals, qi stagnation, and blood stasis. Herbal medicine is one of the most important methods of traditional Chinese medicine in treating cancer pain, and studies have shown that bitter cold herbal medicine has exact antitumor effect, with less toxic and side effects and obvious effect [34, 35]. The clinical analgesic effect of chemical drugs has been fully recognized, but for patients with chronic pain, long-term use will produce dependence, peptic ulcer, and other adverse reactions. Accordingly, with the development of traditional Chinese medicine industry, the research of natural plant extracts for pain treatment is gradually deepening [36, 37].Botanical medicine is an important method for the treatment of cancer pain. This method should be highly valued by medical staff. External application of traditional Chinese medicine is directly administered to patients on the body surface of medical staff. Under the condition of absorption on the skin or mucous membrane surface of the patient, it can reach the pain site of the patient, produce good analgesic effect, and gradually reduce the side effects caused by oral drugs [38]. In this study, 11 plant-derived drugs widely used in clinical pain treatment were selected as follows: Astragalus membranaceus, Poria cocos, Honeysuckle, Cyperus, Paeonia lactiflora, Bupleurum, Hedyotis diffusa, Scutellaria barbata, Trichosanthes kirilowii, Arisaema, and Scutellaria baicalensis. Among them, Astragalus is sweet in taste and mild in nature. It is mainly used for blood numbness of limbs and hemiplegia. Poria cocos is sweet and light, and has a mild nature. Honeysuckle has the characteristics of sweet-cold clearing, light aroma, and drooping, is good for clearing away heat and detoxification, and treats carbuncle sores and boils. Cyperus is pungent, relieves bitterness, and tastes sweet, and it is harmonious, smooth, and unbiased, and works well for soothing the liver, regulating qi, and relieving pain. Paeonia lactiflora has sour, sweet, and bitter taste, it is mild cold in nature, and it has the effect of softening the liver and relieving pain. Bupleurum is bitter and mildly cold, and aromatic, relieves diarrhea, and mainly treats flank pain, irregular menstruation, and dysmenorrhea. Hedyotis diffusa is bitter and cold to clear diarrhea, and sweet and cold to infiltrate and benefit, clear away heat and toxins, dissipate carbuncle, treat sore carbuncle, sore throat, and intestinal carbuncle. Scutellaria barbata has a pungent flavor that can disperse, relieve bitter cold, and clear diarrhea, not only dissipate blood stasis and cool blood to stop bleeding, but also clear heat and water to relieve drenching and swelling, and it is good at treating sores and cancers. Trichosanthes kirilowii is sweet and cold lubricating, and has the effect of clearing lungs, moisturizing dryness and resolving phlegm, reducing swelling, and relieving lumps. Arisaema has bitterness, dryness, and strong medicinal power, and it is especially good at dispelling meridian wind and phlegm, not only drying dampness and reducing phlegm, but also dispelling wind and relieving spasm. Scutellaria baicalensis is bitter cold, clears venting and dryness, and can cool blood and stop bleeding [39, 40].Based on the database of TCMSP, ETCM and TCMID, and OB ≥ 30% and DL ≥ 0.18, 23 Astragalus membranaceus, 20 Poria cocos, 26 Honeysuckle, 18 Cyperus, 13 Paeonia lactiflora, 17 Bupleurum, 7 Hedyotis diffusa, 29 Scutellaria barbata, 11 Trichosanthes kirilowii, 7 Arisaema, and 36 Scutellaria baicalensis were obtained, with a total of 154 active components. Through TCMSP, PharmMapper, SwissTarget database retrieval, and comparison with the UniProt database, 341 protein targets of plant-derived drugs were obtained. Through 1110 targets related to human pain diseases, 141 common targets of plant-derived drugs pain diseases were obtained. The interaction network of 141 potential targets was constructed by using the STRING database platform, and the PPI network diagram was imported into Cytoscape 3.7.2 in the software, the network images were analyzed by tools, and the 19 targets with degree and betweenness centrality >15 were the potential key targets for the analgesic effect of botanical drugs. The main proteins closely related to the analgesic effect of botanical drugs were TNF, MAPK1, JUN, IL-6, and so on. GO enrichment analysis by the DAVID platform showed that the mechanism of pain relief by botanical drugs mainly involved 193 biological processes, 47 molecular functions, and 22 cell components; 118 signaling pathways were obtained by the abdominal muscle analysis of the KEGG pathway. The main biological processes are as follows: RNA polymerase II promoter, transcriptional DNA template, inflammatory reaction, transcriptional positive regulatory DNA template, apoptosis process, immune reaction, active regulation of angiogenesis, regulation of cell proliferation, activation of cysteine-type endopeptidase activity, regulation of apoptosis process, steroid hormone-mediated signaling pathway, lipopolysaccharide-mediated signaling pathway, and sequence-specific DNA binding transcription, positive regulation of NF-kappa B transcription factor activity, negative regulation of cell proliferation, negative regulation of transcription, and so on. Cell components include nucleus, cytoplasm, extracellular space, extracellular exosomes, cytosol, organic components of plasma membrane, intracellular mitochondria, external plasma membrane, endoplasmic reticulum, cell surface, extracellular matrix, postsynaptic membrane, complex receptors, membrane rafts, neuronal projection, cell junction, nuclear chromatin, spindle, RNA polymerase II transcription factor complexThings. Molecular functions include zinc ion binding, ATP binding transcription factor activity, DNA binding, heme node, cytokine activity, protein dimer activity, RNA polymerase II core promoter proximal region-specific sequence DNA binding, chromatin binding, steroid hormone receptor activity, growth factor activity, iron ion binding, transcriptional activator activity, steroid binding, transcriptional regulatory region DNA binding, amino acid endopeptidase activity, adrenaline binding, monooxygenase activity metalloendopeptidase activity, norepinephrine binding, aromatase activity, tumor necrosis factor receptor binding, MAPK activity, peroxidase activity, transmembrane receptor protein tyrosine kinase activity, core promoter sequence-specific DNA binding, drug binding, and chemokine activity. KEGG signaling pathways mainly include cancer, hepatitis B, PI3K-Akt signaling pathway, proteoglycan in cancer, influenza, tumor necrosis factor signaling pathway, pulmonary tuberculosis, MAPK signaling pathway, HTLV-I, toxoplasmosis infection, prostate cancer, Toll-like receptor signaling pathway, adhesive plaque, HIF-1 signaling pathway, hepatitis C, Ras signaling pathway, pancreatic cancer, small cell lung cancer, T-cell receptor signaling pathway, osteoclast differentiation, nonalcoholic fatty liver disease (NAFLD), herpes simplex infection, and cAMP signaling pathway ssc05206: microRNAs in tumors. It can be seen from the above signaling pathways that plant-derived drugs not only regulate the expression of inflammatory factors through inflammatory signaling pathways to relieve pain, but also actively participate in prostate cancer, pancreatic cancer, small cell lung cancer, and other cancer pathways, in order to achieve the effect of relieving cancer pain.


## 5. Conclusion

In conclusion, the material basis and the possible mechanism of action of the active components of plant-derived drugs were preliminarily discussed through network pharmacology. In this article, the effective chemical components and corresponding targets of plant-derived drugs were obtained by TCMSP, ETCM, TCMID, and other databases. GeneCards, OMIM, DrugBank, and other databases were used to screen out pain-related targets. Cytoscape software and STRING database were used to construct the “compound target” interaction network and protein interaction (PPI) network to screen the key components and key targets for therapeutic effects. GO and KEGG pathway enrichment analysis was performed on key targets through the DAVID platform, and the potential action mechanism of 11 herbal extracts of plant origin drugs was learned. We found 153 active ingredients from botanical drugs by TCMSP, ETCM, and TCMID databases, covering 341 protein targets in the human body. Combined with OMIM, GeneCards, and DrugBank databases, we excavated and matched 141 targets of plant-derived drugs and cancerous pain diseases. Through the analysis of the STRING platform and Cytoscape software, 19 core targets including TNF, MAPK1, JUN, and IL-6 were screened out. Go and KEGG enrichment showed that plant-derived drugs alleviated cancer pain processes involving 193 biological processes, 47 molecular functions, 22 cell components, and 118 signaling pathways. By screening genes involved in the KEGG signaling pathway, it was found that plant-derived drugs were mainly associated with PI3K-Akt signaling pathway, tumor necrosis factor signaling pathway, MAPK signaling pathway, Toll-like receptor signaling pathway, and HIF-1 signaling pathway in alleviating cancer pain. These results indicate that botanical drugs can positively affect the expression of inflammatory factors and apoptotic factors in the process of treatment and relief of cancer pain, which is expected to have a potential therapeutic effect on the relief of cancer pain. Although TCMSP and other databases were used to obtain the chemical components of plant-derived drugs, these data are not complete. In the later stage, UPLC, UPLC-Q-TOF-MS, and other methods should be used to detect the chemical components of plant-derived drugs, in order to obtain more comprehensive chemical information to study the mechanism of action of traditional Chinese medicine monomer.



---

*Source: 9326373-2022-02-27.xml*