Abstract

Atrophic rhinitis (AR) is a chronic disease that causes severe structural changes to the nasal mucosa leading to squamous epithelial metaplasia. However, treatment regarding AR remains a major challenge. We used network pharmacology and molecular docking methods to explore the potential mechanisms of the Yiqi Qingre Ziyin method to modulate neuropeptides in the treatment of AR. The active ingredients of the Yiqi Qingre Ziyin method and their targets of action were obtained from the Traditional Chinese Medicine Systematic Pharmacology Database Analysis Platform (TCMSP). Disease targets for AR were obtained from four databases: GeneCards, PharmGKB, DrugBank, and Online Mendelian Inheritance in Man (OMIM). A total of 59 active ingredients, 39 potential targets, and 76 relevant neuropeptides were obtained after deduplication. We constructed target interaction networks with the STRING database. Gene Ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed on the 14 potential target proteins. We used Cytoscape software to construct the “drug-active ingredient-potential target” and “ingredient-target-pathway” networks of the Yiqi Qingre Ziyin method for treating AR. Molecular docking results suggest that dipeptidyl peptidase 4 (DPP4), opioid receptor gene d1 (OPRD1), and opioid receptor m1 (OPRM1) are key targets for the Yiqi Qingre Ziyin method. Therefore, this study proposed a potential mechanism for the treatment of AR by affecting the expression of neuropeptide-related genes (including DPP4, OPRD1, and OPRM1), which may potentially improve the immune microenvironment of the nasal mucosa.

1. Introduction

Atrophic rhinitis (AR) is a chronic disease that causes severe structural changes to the nasal mucosa leading to squamous epithelial metaplasia, loss of cupped cells, and occlusion of the vascular cavity [1, 2]. Chronic infections, autoimmune diseases, endocrine disorders, nutritional deficiencies, such as fat-soluble vitamin deficiencies, and hereditary or iron deficiency anemia are suspects as causative agents of this disease. However, a consensus has not been reached [3, 4], leading to unsatisfactory results in modern medical treatment, which is mostly based on conservative treatment with nasal rinses, drops, antibiotics, antibacterial drugs, immunosuppressants, and vasodilators. Our previous studies have found that AR may be associated with septoplasty [5]. Surgical treatments are available to reduce the size, promote regeneration of normal tissue, and improve blood circulation in the nasal cavity. However, the benefits of surgical treatments are uncertain and can be painful. Besides, symptomatic drugs such as immunosuppressive drugs have certain toxic side effects [3]. These factors make AR a progressively difficult disease to treat. The internal conditioning methods of traditional Chinese medicines (TCM) have certain advantages in treating such diseases. Patients are opting for them increasingly [6].

In Chinese medicine, AR, also known as “Bigao,” is commonly associated with the deficiency of the “lung and spleen.” The identification and typing of diseases is the core of diagnosis in the theory of Chinese medicine. The treatment should benefit the Qi, clear heat, and nourish Yin. However, that theory of Chinese medicine was performed at the subjective level, lacking conclusive evidence. In previous studies, scholars found that neurogenic inflammation tests are more specific than traditional eosinophil counts or sIgE tests. And there is a stronger correspondence between these tests and TCM evidence [79]. With the development of network pharmacology and structural biology in recent years, a reliable approach has been provided to explore the mechanism of AR progression.

The sympathetic and parasympathetic nervous systems regulate the complex structure of the nasal mucosa. However, it is also controlled by the nonadrenergic noncholinergic nervous system (NANC) [1012]. NANC mainly consists of sensory nerves like nerve fibers. Generally, the nerve fiber receptors are isolated because of the integrity of the nasal mucosa epithelium. AR disrupts the epithelial integrity and exposes the receptors, resulting in the production of large amounts of sense neuropeptides [13]. Therefore, we hypothesized that neuropeptides may play an important role in the development and progression of AR.

This study explored the mechanism of the Yiqi Qingre Ziyin method for the treatment of AR through network pharmacology and molecular docking.

2. Methods

2.1. Drug Formulations

Based on the literature review, an herbal prescription was summarized [9]. The prescription consisted of radix astragali (20 g), Poria (15 g), Atractylodes macrocephala Koidz (AMK) (12 g), angelica (12 g), scrophulariae (12 g), Ligustrum lucidum Ait (LIA) (15 g), kudzu (12 g), and red peony root (RPR) (9 g). Codonopsis pilosula (20 g) was added to the prescription when patients presented with significant “qi” deficiency syndrome. In previous research, ingredients and targets of these herbs were downloaded from the traditional Chinese medicine systems pharmacology database [1417]. An ingredient was defined as active when a combination of the following conditions was true: and .

2.2. Screening of Chinese Herbal Ingredients and Target Proteins and Construction of a Chinese Herbal-Component-Target Visualization Network Map

Drug targets were converted to gene names using the UniProt (https://www.uniprot.org/) database. Proteins related to AR were gathered from four datasets, including GeneCards (https://www.genecards.org/), OMIM (https://omim.org/), PharmGKB (https://www.pharmgkb.org), and DrugBank (https://www.drugbank.com). The aggregation and repetition cut-out work was finished using R software [18]. Then, Cytoscape 3.7.1 was used to construct the herbal-component-target visualization network diagram [19].

2.3. Extraction of Key Targets and Construction of Protein-Protein Interaction Networks

Potential targets were imported into the STRING 11.0 database (https://string-db.org/) as previous researches [20, 21]. Species were set to Homo sapiens. Protein interaction relationships were obtained using a filter with a and removing free nodes. Download the result file in tsv format, and import it into Cytoscape 3.7.1. Use the network analyzer in this software to perform network topology analysis and build a PPI network [22].

2.4. Enrichment Analysis of Pathways

Enrichment analysis was performed for the 14 proteins screened using the Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO). The top-ranked genes and those with values less than 0.05 were visualized with R software as previous researches [2326].

2.5. Molecular Docking Analysis

The molecular docking software AUTODOCK 4 was used to validate the molecular docking of the target with the highest correspondence value for each herbal medicine in the PPI network [27]. The RCSB PDB database (http://www.rcsb.org/) and the ZINC database (http://zinc.docking.org/) were used to obtain the structures of target proteins and ligands. Docking results are presented in the form of binding energy, smaller scores indicating more stable binding [28, 29].

2.6. Molecular Dynamics Simulation (MDS)

The complex system of compound and target protein obtained by molecular docking was used as the initial structure for all-atom molecular dynamics simulations (MDS). Create small molecule force fields using ACPYPE Server [30, 31]. The protein force field is described using CHARMM36 [32]. The truncated tetrahedral TIP3P was added as a solvent cartridge in the system. In this case, Na+/Cl- is used to balance the system charge.

GROMACS software was used in this study to perform the system simulations [33, 34]. Energy optimisation is carried out first. The temperature of the system is then slowly increased from 0 K to 307 K at a fixed volume and a constant rate of temperature increase. To further homogenize the distribution of solvent molecules in the solvent box, the NVT and NPT simulations were run in order at a system maintenance temperature of 307 K. Finally, to perform molecular dynamics simulations, the composite system was simulated under periodic boundary conditions with a system pressure of 1 atm, a 2 fs integration step, and trajectories stored at 1 ps intervals.

3. Results

3.1. Screening of Active Ingredients for Yiqi Qingre Ziyin Method

The research flowchart for this study is shown in Figure 1. We first obtained the active ingredients of the Yiqi Qingre Ziyin method. The proteins corresponding to 147 genes were found to be potential targets for the active ingredients of this method. The proteins corresponding to 76 genes were found to be potential targets of Codonopsis pilosula. A total of 843 disease-associated genes were screened (Figure 2(a)), and 39 genes were identified as potential targets (Figure 2(b)). All targets of Codonopsis pilosula, also named as Dangshen (DS), were found to be associated with neuropeptide-related genes (Figure 2(c)).

3.2. Construction of Component-Target Gene Network and Protein-Protein Interaction (PPI) Networks

For regular AR patients, Yiqi Qingre Ziyin herbal soup is often given orally in clinical practice. A disease-drug-component-potential-target network map was constructed using Cytoscape (Figure 3(a)). As shown in the figure, red, blue, pink, cyan, yellow, green, and reseda indicate Atractylodes macrocephala Koidz (AMK), red peony root (RPR), Poria, kudzu, radix astragali, scrophulariae, and angelica, respectively (Figure 3(a)). Besides, these genes were imported into STRING to build the PPI network (Figure 3(b)). Based on this data, the network topology analysis was performed using Cytoscape 3.7.1 [35]. The network has 14 nodes and 89 edges. A total of 14 key genes were acquired: AKT1, CASP3, CXCL8, HIF1A, HMOX1, JUN, MAPK14, MAPK8, MMP2, MMP9, PPARG, PTGS2, TP53, and VCAM1 (Figure 3(c)).

3.3. GO and KEGG Enrichment Analysis

Fourteen key genes (AKT1, CASP3, CXCL8, HIF1A, HMOX1, JUN, MAPK14, MAPK8, MMP2, MMP9, PPARG, PTGS2, TP53, and VCAM1) were used for enrichment analysis. The enriched biological processes are responses to lipopolysaccharide, molecules of bacterial origin, chemical stress (cellular response), oxygen levels, hypoxia, metal ion, nutrient levels, and oxidative stress (cellular response). The enriched cellular components are membrane raft, membrane microdomain, membrane region, nuclear chromatin, transcription regulator complex, RNA polymerase II, caveola, plasma membrane raft, ficolin-1-rich granule, and ficolin-1-rich granule lumen. The enriched molecular functions are nuclear receptor activity, ligand-activated transcription factor activity, RNA polymerase II general transcription initiation factor binding, steroid hormone receptor activity, heme binding, tetrapyrrole binding, Hsp90 protein binding, general transcription initiation factor binding, RNA polymerase II-specific DNA-binding transcription factor binding, and DNA-binding transcription factor binding. The above results suggest that the Yiqi Qingre Ziyin method presents a multitarget and multipathway action for the treatment of AR (Figure 4(a)) (Supplementary Table 1). In the KEGG metabolic pathway enrichment, a total of 78 relevant metabolic pathways were obtained at . The key KEGG enrichment pathway maps were extracted and visualized in descending order of (Figure 4(b)) (Supplementary Table 2). The main pathways of the Yiqi Qingre Ziyin method for the treatment of AR are AGE-RAGE in diabetic complications, TNF, IL-17, relaxin, endocrine resistance, Th17 cell differentiation, NF-kappa B, neurotrophin, osteoclast differentiation, apoptosis, Toll-like receptor, C-type lectin receptor, nonalcoholic fatty liver disease, Th1 and Th2 cell differentiation, estrogen, alcoholic liver disease, T cell receptor, HIF-1, RIG-I-like receptor, prolactin, leukocyte transendothelial migration, MAPK, sphingolipid, NOD-like receptor, VEGF, Fc epsilon RI, adipocytokine, mitophagy–animal, platinum drug resistance, and B cell receptor.

3.4. Analysis of Neuropeptide-Related Target Gene-Pathway Networks and Molecular Docking Results

Codonopsis pilosula was added to the prescription when patients presented with significant “qi” deficiency syndrome. Of the 76 targeted neuropeptide genes in Codonopsis pilosula (Figure 2(c)), DPP4, OPRD1, and OPRM1 are the target genes for AR (Figure 5(a)). Luteolin, frutinone A, 7-methoxy-2-methylisoflavone, 3-beta-hydroxymethyllenetanshiquinone, and 11-hydroxyrankinidine are effective components of Codonopsis pilosula and may have a potential role in AR therapy (Figures 5(b)5(i)) (Table 1). Drug components (including luteolin, frutinone A, 7-methoxy-2-methyl isoflavone, and 3-beta-hydroxymethyllenetanshiquinone) from DS bound to DPP4 with binding energies of -7.6 kcal/mol, -8.6 kcal/mol, -7.8 kcal/mol, and -8.9 kcal/mol, respectively. 7-Methoxy-2-methylisoflavone, 3-beta-hydroxymethyllenetanshiquinone, and 11-hydroxyrankinidine bound to OPRM1 with binding energies of -8.9 kcal/mol, -8.9 kcal/mol, and -7.4 kcal/mol, respectively. 3-beta-Hydroxymethyllenetanshiquinone also binds to OPRD1 with a binding energy of -8.3 kcal/mol. It is generally accepted that a docking fraction below -4.5 kcal/mol is a binding possibility, and a docking fraction below -6 kcal/mol is strong binding power. The binding docking fractions in this analysis were all below -6 kcal/mol, indicating stable and strong binding. This suggests that the effective components in Codonopsis pilosula have a high binding capacity towards neuropeptide-related proteins DPP4, OPRD1, and OPRM1.

3.5. Molecular Dynamics Simulations Indicate That the Tiny Chemical 3-beta-Hydroxymethyllenetanshiquinone Binds to DPP4 in a Stable Manner

The small chemical 3-beta-hydroxymethyllenetanshiquinone was shown to have the maximum binding stability to DPP4 in the current anticipated results. The purple color indicates the pre-MDS conformation, and the orange color indicates the post-MDS conformation. It can be seen that both molecules before and after MDS maintain the same binding site (Figure 6(a)). 3-beta-Hydroxymethyllenetanshiquinone and DPP4 binding sites were virtually unaltered over MDS, showing that 3-beta-hydroxymethyllenetanshiquinone was always securely bound inside the active pocket (Figure 6(b)). And, after MDS, the small molecule seems to be closer to the protein, indicating that 3-beta-hydroxymethyllenetanshiquinone can bind to DPP4 at this position in a stable manner. As a result, we hypothesize that 3-beta-hydroxymethyllenetanshiquinone, the active component in the Yiqi Qingre Ziyin method, can suppress the development of AR by specifically targeting DPP4.

4. Discussion

The early symptoms of primary AR are mild with no obvious dysfunction. Chinese medicine believes that AR is caused by “wind,” “phlegm,” and “dampness.” Patients present with nasal congestion and large amounts of purulent discharge. This leads to a high rate of missed and delayed diagnoses. The etiology of this disease is complex; however, it is generally believed to be related to bacterial infections, plant nervous disorder, and vitamin B2, vitamin A, and vitamin D deficiency [36]. The disease occurs in adolescence with higher incidences in women and gets aggravated during menstruation. Some scholars believe that this disease is related to endocrine immune dysfunction. Considering the complex etiology of the disease, clinical treatment is often symptomatic. However, the therapeutic efficacy is poor.

In the present study, new neuropeptide-related target gene-pathway networks were proposed using network pharmacology and molecular docking techniques. Through MDS research, the active component 3-beta-hydroxymethyllenetanshiquinone in Yiqi Qingre Ziyin method was demonstrated to suppress the development of AR by target binding to DPP4.

Luteolin, frutinone A, 7-methoxy-2-methylisoflavone, 3-beta-hydroxymethyllenetanshiquinone, and 11-hydroxyrankinidine may serve as key players in the treatment of AR. These drug components have been found in previous studies to play a key role in several diseases. Luteolin has many properties, including antioxidant, antimicrobial, anti-inflammatory, chemopreventive, chemotherapeutic, cardioprotective, antidiabetic, antidiabetic neuroprotective, and antiallergic effects [3741]. Luteolin inhibits interleukin- (IL-) 1β, IL-2, IL-6, IL-8, IL-12, IL-17, TNF-α, interferon- (IFN-) β, and granulocyte-macrophage colony-stimulating factor. It also increases the levels of IL-10, an anti-inflammatory cytokine. 7-Methoxy-2-methylisoflavone significantly reduces depression in patients, and this antidepressant affects both menopausal women and men [4244]. GSEA analysis suggests that the Yiqi Qingre Ziyin method targets the IL-17 signaling pathway, Th17 cell differentiation, which suggests that the Yiqi Qingre Ziyin method may regulate the expression of DPP4 in Th17 cells, thereby alleviating the symptoms of AR.

In this study, DPP4 was found to be associated with AR. It is also a target gene of Codonopsis pilosula. Indeed, as a dipeptidyl peptidase, DPP4 could cleave dipeptides from the n-terminus of dipeptides, including enterotrophic hormone, neuropeptides, and some chemokines. DPP4 was found abundantly expressed in IL-17-producing CD4 T cells (Th17 cells) [45, 46]. DPP4 is an enzymatically active molecule expressed by Th17 cells and has a biologically active enzymatic role that may contribute to the regulation of inflammatory T cell chemotaxis [46]. Studies suggest that DPP4-mediated modifications of enzymatic chemokines may regulate T cell migration. Therefore, upregulation of DPP4 is also associated with disease activity in human autoimmune diseases, such as rheumatoid arthritis and multiple sclerosis [4750]. In the present study, we verified for the first time that the active ingredient 3-beta-hydroxymethyllenetanshiquinone stably binds to DPP4 protein in a computer-simulated system by MDS technique. This suggests that 3-beta-hydroxymethyllenetanshiquinone, the active ingredient in Yiqi Qingre Ziyin method, may be able to inhibit the progression of AR by targeting binding to DPP4.

The m-opioid receptor is encoded by the opioid receptor m1 (OPRM1) gene and is the primary site of action for the most used opioids. The delta-opioid receptor (DOR) is encoded by the opioid receptor gene (OPRD1). Opioid receptors are the targets of a large number of drugs, including heroin and many analgesics. Genetic diversity and DOR expression levels of OPRD1 are associated with a variety of diseases, including substance abuse and addiction, anorexia nervosa (AN), obesity, and Alzheimer’s disease (AD) [5153]. Under normal conditions, DOR is localized within the cell. Inflammation can cause receptors on dorsal spinal cord neurons to migrate to the membrane surface, where they can be agonized by delta opioid agonists [5456]. However, further studies on the role of OPRD1 and OPRM1 in AR are needed, and this study found that they may serve as potential targets for the Yiqi Qingre Ziyin method in the treatment of AR.

Network pharmacology and molecular docking were used to explore the potential mechanisms of Yiqi Qingre Ziyin method-mediated neurogenic peptides for the treatment of AR. Three neuropeptide-related proteins, DPP4, OPRD1, and OPRM1, were found to be potentially associated with the development of AR. Nonetheless, the particular mechanism of these proteins’ action on AR has to be explored further. The active ingredients luteolin, frutinone A, 7-methoxy-2-methylisoflavone, 3-beta-hydroxymethyllenetanshiquinone, and 11-hydroxylancidin from Codonopsis pilosula bind with DPP4, OPRM1, and OPRD1. These could be key targets for the Yiqi Qingre Ziyin method. The results of this study suggest that the Yiqi Qingre Ziyin method has good potential efficacy against AR. In this study, the pathogenesis of AR was found to be associated with OPRD1 and OPRM1. However, the exact mechanism remains to be further investigated.

5. Conclusion

This study revealed that the Yiqi Qingre Ziyin method may improve the immune microenvironment of the nasal mucosa by affecting the expression of neuropeptide-related genes (including DPP4, OPRD1, and OPRM1), which provides a theoretical basis for further elucidating the molecular mechanisms of this method in the treatment of AR.

Abbreviations

AMK:Atractylodes macrocephala Koidz
AKT1:AKT serine/threonine kinase 1
AR:Atrophic rhinitis
BP:Biological processes
CASP3:Caspase 3
CC:Cellular components
CXCL8:C-X-C motif chemokine ligand 8
DS:Dangshen
DPP4:Dipeptidyl peptidase 4
GO:Gene Ontology
HMOX1:Heme oxygenase 1
HIF1A:Hypoxia-inducible factor 1-alpha
KEGG:Kyoto Encyclopedia of Genes and Genomes
LIA:Ligustrum lucidum Ait
MMP2:Matrix metallopeptidase 2
MMP9:Matrix metallopeptidase 9
MDS:Molecular dynamics simulation
MF:Molecular functions
NANC:Nonadrenergic noncholinergic nervous system
OPRD1:Opioid receptor delta 1
OPRM1:Opioid receptor mu 1
PPARG:Peroxisome proliferator-activated receptor gamma
PTGS2:Prostaglandin-endoperoxide synthase 2
MAPK14:Recombinant mitogen-activated protein kinase 14
MAPK8:Recombinant mitogen-activated protein kinase 8
RPR:Red peony root
TCMSP:Traditional Chinese Medicine Systematic Pharmacology Database Analysis Platform
TCM:Traditional Chinese medicines
TP53:Tumor protein P53
VCAM1:Vascular cell adhesion molecule 1.

Data Availability

Genes related to atrophic rhinitis were gathered from 4 datasets, including GeneCards (https://www.genecards.org/), OMIM (https://omim.org/), PharmGKB (https://www.pharmgkb.org), and DrugBank (https://www.drugbank.com). The RCSB PDB database (http://www.rcsb.org/) and the ZINC database (http://zinc.docking.org/) were used to obtain the structures of target proteins and ligands.

Conflicts of Interest

The authors declare there are no competing interests.

Authors’ Contributions

Lixing Lu, Xueran Kang, and Yisheng Chen contributed equally to this work.

Acknowledgments

This work was supported by the Medical and Industrial Project, Shanghai Jiao Tong University School of Medicine (YG2021ZD16), and the Clinical Research Program of 9th People´s Hospital, Shanghai Jiao Tong University School of Medicine (JYLJ028).

Supplementary Materials

Supplementary Table 1: enrichment analysis of GO suggests that the Yiqi Qingre Ziyin method presents a multitarget and multipathway action for the treatment of AR. Supplementary Table 2: in the KEGG metabolic pathway enrichment, a total of 78 relevant metabolic pathways were obtained at . (Supplementary Materials)