Abstract
Background. Qi She Pill (QSP) is a traditional prescription for the treatment of neuropathic pain (NP) that is widely used in China. However, no network pharmacology studies of QSP in the treatment of NP have been conducted to date. Objective. To verify the potential pharmacological effects of QSP on NP, its components were analyzed via target docking and network analysis, and network pharmacology methods were used to study the interactions of its components. Materials and Methods. Information on pharmaceutically active compounds in QSP and gene information related to NP were obtained from public databases, and a compound-target network and protein-protein interaction network were constructed to study the mechanism of action of QSP in the treatment of NP. The mechanism of action of QSP in the treatment of NP was analyzed via Gene Ontology (GO) biological process annotation and Kyoto Gene and Genomics Encyclopedia (KEGG) pathway enrichment, and the drug-like component-target-pathway network was constructed. Results. The compound-target network contained 60 compounds and 444 corresponding targets. The key active compounds included quercetin and beta-sitosterol. Key targets included PTGS2 and PTGS1. The protein-protein interaction network of the active ingredients of QSP in the treatment of NP featured 48 proteins, including DRD2, CHRM, β2-adrenergic receptor, HTR2A, and calcitonin gene-related peptide. In total, 53 GO entries, including 35 biological process items, 7 molecular function items, and 11 cell related items, were identified. In addition, eight relevant (KEGG) pathways were identified, including calcium, neuroactive ligand-receptor interaction, and cAMP signaling pathways. Conclusion. Network pharmacology can help clarify the role and mechanism of QSP in the treatment of NP and provide a foundation for further research.
1. Introduction
Neuropathic pain (NP) is a common condition globally [1]. With a complex pathogenesis, reliable treatments are currently unavailable; consequently, NP has a substantial negative impact on quality of life, and it causes a heavy social burden. NP was defined by the International Pain Society in 2011 as “damaged by the somatosensory nervous system or pain caused by disease” [2]. The clinical symptoms of NP mainly include hyperalgesia, allodynia, and spontaneous pain [3]. NP has a high incidence, accounting for approximately one-fifth of all cases of chronic pain, thereby causing significant damage to patients’ physical and mental health [4]. In addition, NP remains a major focus of both clinical medicine and basic research.
Qi She Pill (QSP) has been approved by the China Food and Drug Administration (approval number: Z20090978), and it is widely available in China [5]. QSP is composed of six Chinese herbs, namely Caulis Sinomenii (CS, Qing Feng Teng), Calculus Bovis (BC, Niu Huang), Chuanxiong Rhizoma (CR, Chuan Xiong), Stephaniae Tetrandrae Radix (STR, Fang Ji), Hedysarum Multijugum Maxim (HMM, Huang Qi), and Moschus (MO, She Xiang). QSP has several effects, such as improving Qi, dissipating phlegm, eliminating stagnation, and relieving pain [6].
Preclinical pharmacological studies reported that QSP has analgesic effects and obvious anti-inflammatory effects against acute and chronic inflammation [7]. Pathology studies found that QSP can suppress NP. Furthermore, clinical trials proved that QSP effectively relieves pain and improves clinical symptoms and that the therapy is safe for treating NP. Phase II and III clinical trials of QSP also verified its clinical effect in patients with NP, and no obvious adverse effects were identified. QSP is widely used to treat NP in China.
The development of network pharmacology [8] has been bolstered by the rapid development of multidisciplinary fields such as pharmacy, computer science, biology, and systems biology [9]. The comprehensive observation of drugs using network pharmacology systems is achieved through network analysis covering a wide range of diseases, genes, targets, and drug networks. The results of this analysis illustrate that a disease can be regulated by a variety of different proteins/genes, and some proteins can also be regulated by a variety of diseases. Such results fully reflect the theory of traditional Chinese medicine (TCM) that the same disease can be treated with different methods and that different diseases can be treated with the same methods. In-depth studies of TCM compounds revealed the material basis and molecular mechanism of Chinese medicine pharmacodynamics [10].
Although the analgesic effects and chemical profiling of QSP have been reported, the active compounds in QSP and the specific molecular mechanisms of QSP in treating NP are unclear, and network pharmacology-based prediction of the bioactive components and their target pathways has not been performed. Therefore, the present study explored the pharmacological mechanisms by which QSP alleviates NP and investigated the relationships among herbs, compounds, and target genes by employing QSP-based network pharmacology. First, we collected information on the bioactive compounds in QSP and retrieved candidate target genes from public databases. Network analysis was then performed to connect compounds in QSP with their target genes and related signaling pathways. Finally, based on the results of network analysis, multicompound, multitarget, and multipathway mechanisms related to QSP were predicted. The workflow is shown in Figure 1.

2. Methods
2.1. Collection and Treatment of Active Compounds and Targets in QSP
The chemical components of CS, BC, CR, STR, and HMM meeting the criteria of oral bioavailability (OB) >30% and drug likeness (DL) >0.18 were selected as candidate active components using the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database [11], and their potential targets were retrieved [12, 13]. The chemical components of MO were determined using the BATMAN-TCM server [14], and the principle of “drug-target similarity” was used to set the threshold “score cutoff = 30” to calculate the potential targets corresponding to chemical components.
2.2. Prediction of the Targets of Active Ingredients and Construction of the Active Ingredient-Target Network Map of QSP
Potential protein targets of QSP were screened using the TCMSP database, and all retrieved target proteins were corrected to the official symbols using the UniProt database. The network diagram of “active component-biological target” was constructed using Cytoscape 3.2.1 software [15]. Nodes in the network represent active components and biological targets. Active components and biological targets were connected by edges.
2.3. Retrieval of Relevant Biological Targets for NP
In this study, “neuropathic pain” was used as the search term to retrieve recognized biological targets using various databases (TTD [16], DrugBank [17], DisGeNET [18], GAD [19], and PharmGKB [20]), and the search results were used to identify currently known biological targets related to NP. Because of the poor uniformity of biological target information sources, the UniProt ID of each target was obtained by correcting the target protein and gene information screened via protein database retrieval [21].
2.4. Construction of the Protein-Protein Interaction (PPI) Network Diagram and Screening of Key Biological Targets
The network diagram of PPIs is an important method for identifying direct and indirect interactions between targets [22]. The PPI network was constructed using the STRING platform, and the species was set as “Homo sapiens [23].” The minimum threshold of interaction was set as “medium confidence = 0.4,” which was used to draw the PPI maps of the active components of QSP and PPI maps of NP-related targets. The “Network Analyzer” function of Cytoscape 3.2.1 [24] was utilized to analyze the topological properties of PPI network and calculate two important topological parameters: “node degree distribution” and “betweenness centrality” of the overall network. The node degree reflects the number of connections with other nodes, whereas betweenness reflects a network of all of the shortest paths through a node in the path. These measures indicate the importance of a node in the network and determine whether a target protein is a “core target.” To study the possible pharmacological effects of a core target, the mean degree and betweenness are taken as the “card value,” and the target between degree and betweenness on the card value is selected as the “core target [25].”
2.5. Gene Ontology (GO) Enrichment Analysis of NP-Related Targets following QSP Treatment
To illustrate the QSP treatment target NP protein in the role of gene function, this study uses David [26] analysis QSP PPI for the treatment of NP protein involved in the network GO function analysis of enrichment. According to the number of targets involves sorting, the first few biological processes establish GO bioprocess enrichment and analyze the possible bioprocesses of core targets in the body.
2.6. Kyoto Gene and Genomics Encyclopedia (KEGG) Pathway Enrichment Analysis of the Targets of QSP in the Treatment of NP
Using DAVID, the protein targets of QSP in the treatment of NP were examined using KEGG signaling pathway enrichment [27]. The correlation between core targets and NP in the KEGG signal pathway and clear QSP possible pharmacological mechanism is researched. By using Cytoscape 3.2.1 software, the KEGG pathway of “drug—active ingredient—target—pathway” network is constructed. The results of the enrichment analysis of visualization were processed by ImageGP software.
3. Results
3.1. Screening of the Active Ingredients in QSP
QSP consists of CS, BC, CR, STR, HMM, and MO. Using the TCMSP database, a total of 361 compounds were retrieved, including 16 in CS, 19 in BC, 189 in CR, 50 in STR, and 87 in HMM. Of the 16 compounds in CS, six satisfied the criteria of OB ≥ 30% and DL ≥ 0.18. Of the 19 compounds in BC, five satisfied the criteria of OB ≥ 30% and DL ≥ 0.18. Of the 189 compounds in CR, 115 satisfied the criteria of OB ≥ 30%, whereas seven satisfied the criteria of OB ≥ 30% and DL ≥ 0.18. Of the 50 compounds in STR, all satisfied the criteria of OB ≥ 30% and four satisfied the criteria of OB ≥ 30% and DL ≥ 0.18. Of the 87 compounds in HMM, 44 satisfied the criteria of OB ≥ 30% and 20 satisfied the criteria of OB ≥ 30% and DL ≥ 0.15. Thus, 42 of the total 361 compounds satisfied the criteria, and 37 compounds were finally obtained after excluding duplicates. BATMAN-TCM was used to excavate 30 active components of MO, and 25 active components were obtained after excluding five components not meeting the score cutoff of 30. The remaining 60 active ingredients were removed after deleting duplicates. The resulting network of the final selected compounds from the six herbal medicines is shown in Table 1.
3.2. Prediction of the Targets of QSP and Construction of an Active Component-Target Network Diagram
Using the TCMSP database, 706 potential targets of CS, BC, CR, STR, and HMM were screened, and 734 potential targets of MO in astragalus musk deer pills were screened using the BATMAN-TCM database. The remaining 444 targets were eliminated after deleting duplicates. Compounds were linked to biological targets, and networks of active-component-target interactions were plotted using Cytoscape 3.2.1 software. Network diagrams more intuitively reflect the interactions of these compounds with targets. It was found that the network diagram contained 504 nodes, and the interaction network between drugs and targets contained 1425 edges. The 504 nodes included 60 active components of QSP. Quercetin acted on 151 targets, and beta-sitosterol acted on 74 targets. PTGS2 had a total of 28 potential targets, followed by 24 targets for PTGS1. This fully reflects the complex network relationship between multiple targets, verifying that QSP functions in a multichannel, multilink, multitarget, and integrated synergistic manner. These genes, proteins, compounds, and herbs are listed in Supplementary Table 1, and the resulting herb-compound-target gene network is shown in Figure 2. The key nodes and topological properties of the active component-target network are shown in Table 2.

3.3. Retrieval of NP-Related Targets
By searching databases such as TTD, DrugBank, DisGeNET, GAD, and PharmGKB, a total of 70 QSP targets related to the occurrence and development of NP were identified after deleting repeated targets, as shown in Supplementary Table 2.
3.4. Construction of the PPI Network for the Treatment of NP by the Active Ingredients of QSP and the Screening and Analysis of Key Targets
The PPI network of the active components of QSP was constructed using STRING (Figure 3), and the results revealed 434 direct or indirect targets and 7070 interactions between the targets. At the same time, the PPI network graph of NP-related targets identified 66 targets directly or indirectly related to the selected targets (Figure 4). In total, 324 interactions were identified between the targets. Cytoscape was used to extract the intersection network in the PPI network graph and delete targets without intersections. The screening principle was that when the active components in the vine acted on an NP target, the active components were associated with NP. In total, 48 potential targets were analyzed. There were 191 interactions between targets. Figure 5 presents the PPI network of the NP-related targets of the active ingredients of QSP. Table 3 presents information about the targets in the PPI network in the treatment of NP by QSP. Cytoscape 3.2.1 software was used to create an active component-target interaction network. The network map more intuitively reflects the interactions of the active components of QSP with NP-related targets.



Analysis of the network of active ingredient-target interactions revealed that the network map contained 90 nodes and 296 edges. The 90 nodes included 42 active components of QSP in blood and 48 targets of QSP. The active ingredients of QSP and NP target network are shown in Figure 6. In the subsequent experiment, Network Analyzer topology properties were used to filter key targets. After a target is screened according to the network topology theory, it can affect more nodes, increasing the efficiency of network information transmission. This study screened 10 key targets: DRD2, CHRM2, β2-adrenergic receptor (ADRB2), HTR7, OPRM1, HTR2A, SLC6A4, HTR2C, CHRM1, and calcitonin gene-related peptide (CALCA). The interaction between the target and the target is 18 sides. Table 4 shows the 10 key nodes and the topological properties of QSP in the NP-related target PPI network.

3.5. GO Biological Process Annotation for NP-Related Targets of QSP
Using the GO bioprocess enrichment analysis function of the DAVID platform, the functions of the 48 proteins in the PPI network of QSP in the treatment of NP were studied, and 53 GO entries were determined using a false discovery rate (FDR) of <0.05. The 53 GO entries are listed in Supplementary Table 3. GO is an analytical method for classifying genes according to their function, which can be divided into three parts: biological process, molecular function, and cell composition. GO analysis identifies important functions that cause a trait to change, and the gene corresponding to the function is revealed. Among them, 35 items were related to biological processes, including the regulation of systemic processes, multicellular biological processes, signal transduction, cell proliferation, and response to stimuli. Meanwhile, seven items related to molecular function were found, including items linked to receptor activity, neurotransmitter binding, drug binding, and adrenaline binding, and 11 related items were related to cell composition, including items related to the plasma membrane, synapses, and dendrites.
3.6. KEGG Pathway Enrichment Analysis of QSP in the Treatment of NP
KEGG pathway enrichment analysis was performed using DAVID. The pathways of 48 proteins involved in PPI network of NP-related targets of QSP were analyzed via pathway enrichment, and 24 signaling pathways were obtained. According to the principle of FDR <0.05, eight signaling pathways, including calcium signaling and neural active ligand-receptor interaction signaling pathways, were identified, indicating that the active components of QSP can act on these signaling pathways in the treatment of NP. Table 5 shows the pathways for QSP in the treatment of NP target. The enrichment analysis results were visualized using Image GP software. The bubble size in the figure represents the number of genes enriched in the pathway. The difference in bubble color represents the degree of enrichment of the target gene in the pathway. The bubble diagram is shown in Figure 7. In this study, 32 genes were enriched in the neuroactive ligand-receptor interaction signaling pathway, 19 genes were enriched in the calcium signaling pathway, and 11 genes were enriched in the cAMP signaling pathway.

3.7. Construction of the Drug-Active Ingredient-Target-Pathway Network Diagram of QSP in the Treatment of NP
The drug-active ingredient-target-pathway network model map of QSP was constructed using Cytoscape 3.2.1 software. The map included 103 nodes, including 5 drugs, 42 active components, and 48 target proteins. The metabolic pathways were connected by 444 edges to form a complete network. Figure 8 presents the drug-active ingredient-target-pathway network diagram.

4. Discussion
TCMSP is a system-level pharmacological database of TCM that provides systematic TCM information, and it is a platform for the development of new drugs for the treatment of human diseases. After a critical review of pharmacological and clinical knowledge, the TCMSP collected all 499 herbs registered in the 2010 edition of the Chinese Pharmacopoeia, totaling 12,144 compounds [11]. BATMAN-TCM is the first online bioinformatics analysis tool of TCM molecular mechanisms, and it is used to clarify the “multicomponent, multitarget, and multipathway” combinatorial therapeutic mechanism of TCM and provide valuable clues for subsequent experimental validation, thus accelerating the elucidation of the molecular mechanisms of TCMs [14]. It is the first step of network pharmacology analysis for identifying chemical components and their target information as completely as possible. The collection of the chemical constituents of QSP is accurate and complete using the two databases currently recognized for network pharmacology research.
QSP can relieve pain and promote blood circulation. The drug has been widely used in China. Clinical studies revealed that QSP can significantly improve mild-to-moderate NP and improve the quality of life of patients with mild-to-moderate NP, especially concerning physiological function, social function, emotional function, vitality, and other aspects. In a number of randomized, double-blind, placebo-controlled studies, improvements in dysfunction associated with mild-to-moderate NP were observed in both the QSP and placebo group, although the improvement was significantly better in the QSP group [5]. To study the mechanism of action of QSP, this study relied on the TCMSP and BATMAN-TCM platforms to study the active constituents of six Chinese herbs (CS, BC, CR, STR, HMM, and MO) and constructed a compound-target network to analyze the direct interactions between compounds and targets. These findings provided a reference for the multicomponent, multitarget, and multipath treatment mechanism of QSP. The results of this study suggest that QSP most likely exerts its therapeutic effects on NP through interactions among MOL021, MOL001, PTGS2, and PTGS1, thereby laying a foundation for subsequent confirmatory studies. The key compounds of QSP in the target network included quercetin and beta-sitosterol, which had much higher degrees than other compounds, and beta-sitosterol was detected in both CS and STR. It can thus be speculated that CS and STR play key roles in the effects of QSP [28]. Prior research indicated that beta-sitosterol has a wide range of pharmacological effects including anti-inflammatory, immunomodulatory, antioxidant, and analgesic effects [29]. Quercetin can lower blood pressure, enhance capillary resistance, and reduce capillary fragility [30]. The key targets of QSP in the key target network included PTGS1 and PTGS2. PTGS is a key rate-limiting enzyme for the synthesis of prostaglandins from arachidonic acid [31], and two isomerases, namely, PTGS1 and PTGS2, have been identified. PTGS and its products are involved in various physiological and pathological processes such as inflammatory responses and blood coagulation balance [32]. Ten key target proteins were screened using the Network Analyzer topology properties in the PPI network for the treatment of NP by the active ingredients of QSP. It was previously found that DRD2 is the major receptor for NP [33], and KLF15 is induced in NP in a TNF-α-dependent manner. KLF15 induces NP by promoting the expression of the DRD2 receptor. NP is associated with persistent changes in gene expression in primary sensory neurons, but the underlying epigenetic mechanisms leading to these changes are unclear. In addition, CHRM plays a key role in the regulation of spinal cord nociceptive transmission [34]. Norepinephrine and the adrenergic system control the transmission of damage at the level of the spinal cord and play a key role in this process. Stimulation of ADRB2 downregulates the phosphorylation of p38 MAPK in microglia and JNK in astrocytes. The aforementioned studies indicated that spinal ADRB2 regulates MAPK phosphorylation and downregulates glial activity, thereby inhibiting neuropathic damage conduction at the spinal cord level [35]. The HTR family is involved in the central nociceptive mechanism, which plays a key role in suppressing analgesic pathways, and its stimulation can regulate neuronal signaling induced by neuropathy [36]. Animal studies revealed that CALCA expressed in dorsal root ganglia is involved in the generation and maintenance of NP [37]. Based on the aforementioned information, it is speculated that the testes directly or indirectly act on these key targets to play a pharmacological role in the treatment of NP.
PPI networks are composed of individual proteins interacting with each other to participate in biological signal transmission, gene expression regulation, energy and substance metabolism, cell cycle regulation, and other aspects of life processes. It is of great significance to understand the functional connection between proteins [22]. Enrichment analysis is helpful for studying gene and expression information as a whole network. In organisms, different gene products coordinate with each other to perform biological functions, and pathway annotation analysis of differentially expressed genes is helpful for further interpreting gene functions [26, 27]. In this study, we preliminarily explored the primary pathway of QSP in treating NP using PPI network and enrichment analyses. In the present study, the analysis yielded the following results. Of the 5 herbs and 42 compounds examined, a total of 48 target genes were associated with NP. To illustrate the role of the active targets of QSP in gene function and signaling pathways, GO functional enrichment and KEGG pathway enrichment analyses of the 48 genes in the PPI network were performed. GO function enrichment analysis found that the effects of QSP in the treatment of NP are mainly reflected in the regulation of receptor activity, plasma membrane, and systemic processes. The main KEGG pathways involved in the treatment of NP by QSP included calcium signaling and neuroactive ligand-receptor interaction signaling pathways. From the aforementioned studies, it is concluded that the active components of QSP alleviate NP by acting on these signaling pathways. A neuroactive ligand-receptor interaction signaling pathway is a collection of receptor ligands on the plasma membrane that are associated with intracellular and extracellular signaling pathways and are closely related physiologically to neurological function [38]. Key molecules in the calcium signaling pathway have a variety of biological functions, among which calcium ions play an important role in maintaining biopotential [39], maintaining normal nerve conduction function [40], and regulating physiological actions including intracellular enzyme activity, muscle contraction and relaxation [41], nervous system function [42], stimulation and secretion-coupled participation [43], glucose metabolism [44], and prostaglandin [45] and insulin synthesis [46].
Network pharmacology combines the concepts of systems biology and multidirectional pharmacology [47]. It explores the relationship between drugs and diseases from a holistic perspective, emphasizing the holistic and systematic interactions among drugs, targets, and diseases [48]. The Chinese herb system is a multicomponent complex system. Bioinformatic and network pharmacology research helps to both clarify the mechanism of action of QSP in the treatment of NP and provide an important theoretical basis for further discussion of subsequent molecular mechanisms and optimization of experimental design. This research also provides new ideas and methods for the research and development of new Chinese herbs and elucidation of the related mechanisms of drugs and diseases [49]. Meanwhile, some limitations of this study must be addressed. First, because of limitations regarding screening conditions, only the main compounds in QSP could be analyzed, and thus the research results are limited to some extent. Second, although a large number of targets and pathways can be screened using network pharmacology, these results need to be verified via pharmacological experiments. Third, further validation experiments are needed to increase the drug concentration to accurately reveal the effective ingredients and mechanism of QSP in the treatment of NP.
5. Conclusion
The network pharmacology analysis of QSP identified 444 compounds in six herbs, and 48 of the compounds in the six herbs were screened to identify NP-related target genes. This study demonstrated that network pharmacology studies can reveal close interactions between multiple components and multiple targets. The network pharmacology approach of this study provides another research strategy for a comprehensive understanding of the mechanism by which QSP treats NP.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (81873317) and Clinical Medical Center Construction Project in Shanghai (no. 2017ZZ01010).
Supplementary Materials
Supplementary Table 1: list of the genes, proteins, compounds, and herbs of QSP. Supplementary Table 2: list of a total of 70 QSP targets related to the occurrence and development of NP. Supplementary Table 3: list of the 53 GO entries. (Supplementary Materials)