Integrating Network Pharmacology and Experimental Validation Deciphers the Mechanism of Guizhi Fuling Wan against Adenomyosis

Wang, Haoxian; Zhang, Jihong; Zhu, Qinqin; Fu, Xianyun; Li, Chenjie

doi:https://doi.org/10.1155/2021/6034147

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results and Discussion Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2021 | Article ID 6034147 | https://doi.org/10.1155/2021/6034147

Integrating Network Pharmacology and Experimental Validation Deciphers the Mechanism of Guizhi Fuling Wan against Adenomyosis

Haoxian Wang,^1,2Jihong Zhang,¹Qinqin Zhu,²Xianyun Fu,²and Chenjie Li¹

Academic Editor: Samuel Martins Silvestre

Received06 Aug 2021

Revised28 Sept 2021

Accepted01 Oct 2021

Published26 Oct 2021

Abstract

Aim. This study aimed to predict the key targets and endocrine mechanisms of Guizhi Fuling Wan (GZFLW) in treating adenomyosis (AM) through network pharmacology, molecular docking, and animal experiment verification. Methods. The related ingredients and targets of GZFLW in treating AM were screened out using TCMSP, BATMAN-TCM, SwissTargetPrediction, and PubChem Database. Then, the protein-protein interaction (PPI) analysis and the network of compound-hub targets were constructed. At the same time, the key targets were uploaded to the Metascape Database for KEGG pathway enrichment analysis. After that, the molecular docking technology of the main active components and hub targets was performed. Furthermore, animal experiments were used to verify the results of network pharmacology analysis. Results. A total of 55 active ingredients of GZFLW and 44 overlapping targets of GZFLW in treating AM were obtained. After screening, 25 hub targets were collected, including ESR1, EGF, and EGFR. Then, the KEGG pathway enrichment analysis results indicated that the endocrine therapeutic mechanism of GZFLW against AM is mainly associated with the estrogen signaling pathway, endocrine resistance, and an EGFR tyrosine kinase signaling pathway. Then, molecular docking showed that the significant compounds of GZFLW had a strong binding ability with ERα and EGFR. More importantly, the animal experiments confirmed that the GZFLW could downregulate the abnormal infiltration of the endometrial epithelium into the myometrium and had no interference with the normal sexual cycle. This effect may be directly related to intervening the local estrogen signaling pathway of the endometrial myometrial interface (EMI). It may also be associated with the myometrium cells’ estrogen resistance via GPER/EGFR signaling pathway. Conclusion. The endocrine mechanism of GZFLW in treating AM was explored based on network pharmacology, molecular docking, and animal experiments, which provided a theoretical basis for the clinical application of GZFLW.

1. Introduction

In 2020, 1.03% of the women were suffering from adenomyosis (AM), and approximately 90% of them were accompanied by menorrhagia, progressive dysmenorrhea, pelvic pain, dysmenorrhea, and infertility [1]. AM is characterized by the progressive invasion of endometrial cells into the myometrium, attributed to endocrine disorders, chronic hypoxia, and inflammation in the endometrial myometrial interface (EMI) [2, 3]. Current therapeutic approaches, including operative treatment, nonsteroidal anti-inflammatory drugs, and endocrine therapy, have been used to suppress inflammation, relieve pain, and improve fertility. However, several limitations including poor patients’ response, high prices, and long-term side effects [4, 5]. Therefore, it is of great clinical significance to develop alternative therapeutic strategies to improve the outcome of AM patients.

Traditional Chinese medicine (TCM), which has been widely used in clinical practice in China for thousands of years, has made considerable contributions to the relief of AM [6]. Guizhi Fuling Wan (GZFLW) is an effective Chinese medicinal formula for AM which is used in TCM [7, 8]. This medication consists of the following five herbs: Cinnamon Twig (CT), Poria cocos (PC), Cortex Moutan (CM), Radix Paeoniae Rubra (RPR), and Peach Kernel (PK). The complex composition of GZFLW makes it not easy to conduct mechanistic studies on it. Thus, the “multicomponent and multitarget” therapeutic mechanism of GZFLW in treating AM is still obscure, hindering its clinical use. Hence, for a comprehensive understanding of the mechanism of action of GZFLW, a novel analysis method is needed [9].

In this study, a network pharmacology strategy was adopted for investigating the mechanism of GZFLW against AM. We screened representative compounds in five herbs of GZFLW. Their potential targets in treating AM were identified. Furthermore, molecular docking technology was used to verify the binding ability of active ingredients to several potential targets. In addition, animal experiments have further confirmed some of the hypotheses. The results revealed that GZFLW might modulate the estrogen signaling pathway, endocrine resistance, and EGFR tyrosine kinase signaling pathway to exert an effect. These results may provide novel sights into the anti-AM effect of GZFLW.

2. Materials and Methods

2.1. Chemical Ingredients Database Building

The active ingredients of the five herbs contained in GZFLW were obtained using the Traditional Chinese Medicine Systems Pharmacology Database (TCMSP, http://lsp.nwu.edu.cn/tcmsp.php) and the Bioinformatics Analysis Tool for Molecular Mechanism of Traditional Chinese Medicine (BATMAN-TCM, http://bionet.ncpsb.org/batman-tcm/index.php), which are both bioinformatics analysis tools for the main components of TCM.

2.2. Active Ingredients Screening

The oral bioavailability (OB) and drug-likeness (DL) were used as screening indexes. OB is the percentage of oral drugs absorbed into the systemic circulation, which is a frequently used pharmacokinetic parameter. DL refers to the property of measuring the absorption, distribution, metabolism, and excretion (ADME) of drug molecules, which could help optimize pharmacokinetic and pharmaceutical properties. The OB and DL indices of all the related ingredients are presented in the TCMSP. We set the indexes of OB ≥ 30% and DL ≥ 0.18 as the screening criteria of active ingredients. Briefly, components that satisfied the screening threshold of OB and DL were regarded as the candidate compounds.

2.3. GZFLW Targets Prediction

The targets of active ingredients were searched in TCMSP, BATMAN-TCM, SwissTargetPrediction Database (http://www.swisstargetprediction.ch/index.php), and PubChem Database (https://pubchem.ncbi.nlm.nih.gov/), with the species limited as ‘‘Homo sapiens.’’

2.4. Adenomyosis Therapeutic Targets Prediction

Targets associated with AM were collected by GeneCards Database (https://www.genecards.org/) and Online Mendelian Inheritance in Man Database (OMIM, https://www.omim.org). We searched these databases with the keyword “adenomyosis” to obtain related targets. After deletion of duplications, candidate targets against AM were identified by matching the targets of the active ingredients of GZFLW and the genes related to AM. These targets were then uploaded to STRING Database (https://string-db.rog/) to obtain the protein-protein interaction (PPI) information between two proteins or multiple proteins with the confidence ≥0.7.

2.5. Network Topological Feature Set Definitions

For each node in the interaction network, we selected two parameters to evaluate its topological features: “Degree” reflects the number of other nodes interacting with this node; the percentage of all shortest paths measures “Betweenness Centrality” (BC). Nodes greater than or equal to the median are screened as hub nodes for GZFLW against AM.

2.6. KEGG Enrichment Analysis

We used the Metascape Database (http://metascape.org/) to obtain the potential pathways of GZFLW in treating AM. was regarded as the cutoff value.

2.7. Molecular Docking of Hub Targets and Active Ingredients

Molecular docking is a computational tool to predict the binding ability and connection type of proteins and ligands. It can calculate and predict the conformation and directions of ligands at active protein sites. AutoDock 1.5.6 (http://autodock.scripps.edu/) was a molecular docking software for network pharmacology based prediction and analysis, which could be used to dock the hub targets and active ingredients. In this process, structures files in the formats of 3D of hub targets we chose were obtained from RCSB Protein Data Bank (PDB, http://www.rcsb.org/): EGFR (PDB ID: 2GS2) and ERα (PDB ID: 6CHW); structures files in the formats of 3D of active compounds were provided from ZINC Database (http://zinc.docking.org/) and PubChem Database (https://pubchem.ncbi.nlm.nih.gov/). During the docking process, the PDBQT files of these structures were used in AutoDock software, as these files contained the information of potential ligands and receptors required for the process. After that, the Lamarckian genetic algorithm was used to obtain docking results and the binding energy of ligands and receptors. At last, the results were output in .dlg files and uploaded to PyMOL software for establishing the 3D results binding models. Binding energy was used as a docking score to evaluate the protein-ligand binding potential of active ingredients to hub targets. Among them, those results with value ≤ −5 were selected and considered to have moderate binding potential and tight combination.

2.8. HPLC and Mid-Infrared Spectroscopy Analysis

The HPLC finger printing of CT, RPR, PK, and CM in GZFLW was constructed using Series LC-2010A Liquid Chromatograph (Shimadzu, Japan). Octadecyl silane bonded silica was used as the stationary phase (column: hypersil ODS C18 4.6 mm × 250 mm, 5 μm). The mobile phase was ammonium acetate (A) and water (B). The flow rate was kept constant at 1 mL/min. The system operated at 30°C and the injection volume was 10 μL (CT, RPR, and PK) and 5 μL (CM). The detection wavelength was, respectively, kept at 285 nm (CT), 230 nm (RPR), 274 nm (CM), and 210 nm (PK). The reference substances (cinnamic acid, cinnamaldehyde, gallic acid, catechin, paeonol, paeoniflorin, and amygdalin) were prepared by Chengdu Ruifensi Biotechnology Co., Ltd. (Chengdu, China).

The active ingredients of PC were analyzed by mid-infrared spectroscopy. The ingredient was weighed, adding about 100 mg KBr powder which was dried at 120°C for 4 h. Then, grind and mix well in the mortar, transfer to the mold, and press the samples into transparent slices at a pressure of about 30 MPa for 1∼2 min. The sample slice thickness was about 1∼2 mm. The Fourier transform infrared spectrometer (deuterated triglycine sulfate detector) was used with a spectral range of 4000∼400 cm⁻¹. Each sample was repeated twice.

2.9. In Vivo Experiment

Virgin female ICR mice (26–29 g, 7 weeks old) were provided by the Experimental Animals Center of China Three Gorges University. All the animals were housed at a temperature of 22 ± 2°C, humidity of 55–65%, and light controlled (12 h light/12 h dark cycle) vivarium, with food and water ad libitum. Animal welfare and experimental procedures were carried out in accordance with the Guide for the Care and Use of Laboratory Animals and were approved by the China Three Gorges University Medical Animal Ethics Commission (Approval No. 20190801).

Twenty-five mice were implanted with a single pituitary gland each under the capsule of the right uterus at 7 weeks of age and then were divided randomly into 2 groups: model group (n = 14) and GZFLW group (n = 11). After intraperitoneally anesthetizing with emulsion (Xian Nippon Pharmaceutical Co., Ltd., Xian, China, No. H19990282, 100 mg/kg), pituitary obtained from age-matched male mice was grafted into the right uterine lumen in female mice. Then, the gentamycin solution (0.25 ml, 20000 units/20 g) was dropped intraperitoneally to prevent infection. Mice receiving no pituitary grafts served as controls (n = 8). One day after the pituitary grafting, 11 mice were intragastriced with decocting-free granules of GZFLW (4.25 mg/g of body weight, consisting of CT, PC, CM, RPR, and PK at the ratio of 1 : 1:1 : 1:1, Beijing Kangrentang Pharmaceutical Co., Ltd., Beijing, China, Nos. 19026231, 19027291, 19019661, 19004891, and 19006412) and dissolved in aqueous solution every day for 8 weeks (GZFLW group). The remaining 22 mice (model group 14, control group 8) were intragastriced with the same dosage of saline as a control vehicle. All mice were killed at 15 weeks of age by decapitation, after the estrous cycle had been checked by vaginal smear.

Haematoxylin and eosin (HE) staining was used to observe the depth of invasion. Immunofluorescence was adopted to identify the expressions of ERα, GPER, and EGFR in the glandular, stromal, and muscular parts in the endometrial myometrial interface (EMI) (magnification, 200×, every microscope field taking three different parts for analysis).

2.10. Data and Statistical Analysis

Data were analyzed using SPSS 20 software and presented as mean ± standard deviation (SD). Statistical differences between groups were evaluated by one-way analysis of variance (ANOVA). The Chi-square test was applied to analyze the categorical variables. Pearson analysis was used for correlation. was regarded as a statistically significant difference.

3. Results and Discussion

3.1. Active Ingredients and Targets of GZFLW

By retrieving from the databases, there were 498 related ingredients of GZFLW in total: 222 ingredients in CT, 33 ingredients in PC, 56 ingredients in CM, 121 ingredients in RPR, and 66 ingredients in PK. After eliminating the redundancy and screening with OB and DL, 55 active components and 364 corresponding targets of GZFLW were obtained. It is worth noting that cinnamaldehyde (OB = 31.99%, DL = 0.02), cinnamic acid (OB = 19.68%, DL = 0.03), benzoic acid (OB = 31.55%, DL = 0.02), paeonol (OB = 28.79%, DL = 0.04), gallic acid (OB = 31.69, DL = 0.04), and amygdalin (OB = 4.42%, DL = 0,61) have a relatively low OB or DL value, but they were included in this study because they were the vital pharmacological compounds identified, which were confirmed to the previous study (Table 1) (Table S1).

3.2. Target Analysis

Via the keyword of “adenomyosis,” the therapeutic targets of AM were obtained from GeneCards and OMIM Database with a total of 211 (Table S2). A total of 48 overlapping targets were then obtained as the related targets of GZFLW in treating AM (Figure 1).

From the STRING Database of “multiple proteins” interaction of targets, the data of the protein-protein interaction (PPI) network of those 48 targets were subsequently obtained (Table S3). With confidence ≥0.7, we gained 268 edges of 48 nodes. With betweenness ≥14 and degree ≥0.01, the topological analysis of targets mentioned above was performed. The results were 25 hub nodes with 183 edges, in which 25 nodes represented 25 central targets, including ESR1, EGF, and EGFR. When the 25 significant nodes and the other 23 nodes were distributed with “degree” and “betweenness,” the network of 48 nodes was constructed as in Figure 2. The comprehensive network of “Herbs-Compounds-Hub Targets” had 25 nodes related to 26 ingredients, as determined by the Cytoscape 3.7.1 software (Figure 3).

Figure 2

The network of 48 nodes. The different color represents different nodes from each ingredient: The blue nodes represent the targets from CT, the pink nodes represent the targets from CM, the purple nodes represent the targets from RPR, the orange nodes represent the targets from PC, and the yellow nodes represent the targets that are targeted by more than one ingredient. The node size is proportional to the target degree in the network. The edge color changes from light to dark reflect that the betweenness value changes from low to high in the network.

3.3. KEGG Pathway Enrichment Analysis

In the Metascape Database, we performed the KEGG pathway enrichment analysis on 48 targets and obtained 114 pathways (Table S4). We ranked these pathways according to the P-value of each enriched pathway in ascending order, and the top 20 pathways were selected, which mainly related to endocrine resistance, EGFR tyrosine kinase inhibitor resistance, and estrogen signaling pathway (Figure 4).

3.4. Molecular Docking

To further explore the binding ability of active components and hub targets in GZFLW, molecular docking technology was used. Due to the fact that ERα and EGFR were the core targets of both estrogen pathway and endocrine resistance pathway, the two were selected to bind to active ingredients. 71.43% of the compounds (10/17) had binding energies lower than −5, which verified that those active components may be well binding with ERα and EGFR (Table 2). The 3D binding mode of 2 major active ingredients in the active site of ERα and EGFR, respectively, is represented in Figure 5.

(a)

(b)

(c)

(d)

3.5. HPLC Analysis of CT, CM, RPR, and PK

There were six main peaks in the fingerprint of CT at 285 nm, and two active components were determined after comparison with standards; the ten main peaks in the fingerprint of CM at 274 nm were obtained, and three vital components were confirmed; there were five main peaks in the fingerprint of RPR at 230 nm, and, among these, three active ingredients were determined; the results in the fingerprint of PK at 210 nm were four main peaks, defining one significant ingredient (Figure 6). All of the main peaks had good separation.

(a)

(b)

(c)

(d)

3.6. Mid-Infrared Spectroscopy

Around 3400 cm⁻¹, these frequencies were associated with O–H stretching. Around 2800 cm⁻¹, the peaks of C-H antisymmetric stretching of methylene group were observed. Frequencies around 1600 cm⁻¹ and 1450 cm⁻¹ were associated with O-H and C-H bends. Below 1300 cm⁻¹, the frequencies were associated with C-O stretching, with frequencies between 1200 and 1000 cm⁻¹ showing the presence of polysaccharide. In addition, frequencies below 950 cm⁻¹ were related to sugar ring skeleton stretching (Figure 7).

3.7. Effect on Infiltration Depth in Mice

AM model was successfully established after pituitary grafts, and HE staining indicated that endometrial stroma and glandular epithelial cells infiltrated into the myometrium. According to the infiltration depth, the models can be divided into mild, moderate, and severe infiltration (Figure 8). The results showed that the proportion of moderate and severe infiltration decreased after GZFLW treatment.

3.8. Effect on the Sexual Cycle in Mice

The estrous cycle of mice can be divided into four stages: proestrus (mainly nucleated epithelial cells), estrus (mainly nonnucleated keratinocytes), metestrus (mainly a small number of keratinocytes and more leukocytes), and diestrus (mainly a large number of leukocytes). We divided the estrous cycles into two phases, combining metestrus with diestrus because of the leucocytic predominance [10]. There is no significant difference in vaginal smears among the three groups by Chi-square test, which confirmed that GZFLW treatment had no critical interference in the sexual cycle of mice (Figure 9).

3.9. Immunofluorescence Assay

The immunofluorescence assay results showed that ERα and GPER were expressed in glandular epithelial cells, stromal cells, and muscularis cells. EGFR was expressed in stromal cells, muscularis cells, and some glandular epithelial cells (Figure 10). Compared with the control group, the expression of ERα in the gland, stroma, and muscularis of the model group was significantly enhanced; after GZFLW treatment, the expression of ERα in the gland and muscularis was obviously decreased and there was no significant difference between the control group and GZFLW group. Compared with the control group, the expression of GPER in the gland was markedly enhanced in the model group; after GZFLW treatment, the expression was considerably decreased. The expression of EGFR in the muscularis was significantly higher in the model group when compared to the control group, and the expression of EGFR in the muscularis of GZFLW group was significantly lower compared to that in the control group. Moreover, the expressions of GPER and EGFR in the muscle layer were positively correlated (Figure 11).

Figure 10

Immunofluorescence staining of ERα, GPER1, and EGFR at the EMI (×200). (a) ERα expression in glandular, stromal, and muscle layer of EMI. (b) GPER expression in glandular, stromal, and muscle layer of EMI. (c) EGFR expression in the stromal and muscle layer of EMI. The cells indicated by yellow arrow represent gland, the cells indicated by white arrow represent stroma, and the cells indicated by green arrow represent muscularis. Compared with the control group, ; compared with the model group, .

4. Discussion

Due to AM’s complex pathophysiology, a totally resolving cure is yet to be found. The efficacy of TCM in treating AM has been confirmed through several hundred years’ practices and experimental research [11]. Previous studies have found that TCM treatment can alleviate dysmenorrhea, shorten the menstrual time, decrease menstrual volume, reduce uterine volume, and lower CA125 levels in AM patients [12–14]. GZFLW is one of the classic traditional Chinese herbal formulas widely used in the clinical treatment of AM [8], and studies have confirmed that the side effect of GZFLW on AM is only 1.99% [15]. We found that the oral treatment of GZFLW exhibited a therapy effect on AM mice by suppressing the depth of invasion without disruption to the normal sexual cycle.

As a hormone-dependent disease, the disorder of estrogen and its receptor is known to promote and exacerbate AM [16, 17]. It has been reported that AM is usually accompanied by a high level of estrogen [17], and the critical enzyme expression of estrogen metabolism regulation in ectopic endometrial stromal cells is significantly increased [18]. Targeting estrogen has excellent potential in the therapy of AM [19, 20]. In our study, 9 out of 55 compounds from GZFLW may act as estrogen regulators, including beta-sitosterol, hederagenin, quercetin, paeoniflorin, kaempferol, taxifolin, pachymic acid, ellagic acid, and baicalein. Quercetin and baicalin are flavonoids derived from GZFLW component herbs RPR. Both of them can inhibit the proliferation and migration, consequently depress the invasion of AM ectopic endometrial stromal cells, and significantly alleviate pain allergy symptoms [21–23]. Kaempferol and beta-sitosterol, derived from GZFLW component herbs including CM, CT, RPR, and PK, have been shown to have antiproliferative effects on mice models of endometriosis diseases [24–26]. The effect of reducing estrogen-dependent abnormal proliferation has been demonstrated to be mediated via a downregulation of ERα [21, 27–29]. Paeoniflorin can downregulate the ERα level of AM primary cells [30]. Several active components, including taxifolin, beta-sitosterol, baicalein, and hederagenin, were validated as the potential antiestrogen compounds in GZFLW by molecular docking because of their high binding ability with ERα. The animal experiments found that ERα expression in AM mice model was significantly upregulated in glandular epithelial cells, stromal cells, and myometrial cells. It is acknowledged that endocrine regulatory treatments, such as oral contraceptives, oral and nonoral progesterone, danazol, and gonadotropin-releasing hormone (GnRH) analogs, have been used to control menstrual pain and menorrhagia by inhibiting uterine ERα expression. Still, they interfere with the menstrual cycle of AM women [31]. Those therapies lead to prolonged diestrus stages and even persistent diestrus in AM mice [32], which was associated with the cycle lengthening and disorder [33]. Our data showed that the ERα expression decreased and the pathological infiltration depth reduced in GZFLW group, while the diestrus/metestrus stages time was not prolonged. The result suggested that the treatment of TCM did not interfere with the estrous cycle in mice, which may be attributed to the selective action of the GZFLW on the local lesion.

Of note, the endocrine resistance pathway and EGFR tyrosine kinase signaling pathway may also contribute to the therapeutic mechanism of GZFLW against AM predicted by the network pharmacology analysis. Indeed, endocrine resistance is a complex pathophysiological process, known to lead to unresponsiveness of the regular endocrine therapy in a particular proportion of AM patients. It was reported that the EGFR gene polymorphism is involved in the endometriosis’s susceptibility [34]. EGF/EGFR is of importance to the migration and invasion ability of endometrial stromal cells (ESC) [35], the reduction of which would inhibit the proliferation and vitality of AM ESC [36]. Subsequently, the results based on the molecular docking showed that the main compounds of GZFLW had a good binding ability with EGFR. In line with this observation, studies have shown that GZFLW can reduce the level of EGF in rats, thereby inhibiting the abnormal proliferation [37]. It has been proved that quercetin [38] can lower blood concentrations of EGF. Baicalein [39] and taxifolin [40] have also been reported to inhibit EGFR signaling by repressing EGFR-Tyr1068 autophosphorylation and then interfere with EGFR associated downstream signaling transduction. The results of our animal experiments showed that EGFR was upregulated in the model group compared to the control group and remarkably decreased after GZFLW treatment, which confirmed that EGF/EGFR is a potential target of GZFLW in treating AM. Besides, we proved there was a positive correlation between the expressions of GPER and EGFR, which may contribute to endocrine resistance. Indeed, endocrine resistance is a complex pathophysiological process, known to lead to unresponsiveness of the regular endocrine therapy in a particular proportion of AM patients [41]. In addition to ER deletion mutation, the EGF tyrosine kinase signaling pathway activated via GPER is another important cause of acquired estrogen tolerance. Previous studies have found that the gene polymorphism of GPER affects the risk of AM [42], and the expression of GPER is significantly increased in AM [43]. GPER can promote heparin-binding EGF release from the cell surface via activating Src family tyrosine kinases, and the activation of the EGF/EGFR signaling pathway eventually leads to epithelial-mesenchymal transition in AM [44]. In addition to the epithelial cells of the uterine gland, a similar trend was observed in the myometrium. Myometrium has the ability of stem cell-like abnormal differentiation [45]. Considering that the metaplasia of uterine smooth muscle cells was regarded as one crucial reason for the occurrence and development of AM [46], we speculate that GZFLW leads to the downregulation of the GPER/EGFR signaling pathway, which subsequently suppresses the estrogen resistance of smooth muscle cells and inhibits the abnormal differentiation of them.

In conclusion, based on the results of network pharmacology and animal experiments, we think that GZFLW can synergistically inhibit the abnormal proliferation and differentiation of AM via downregulating the ERα estrogen signaling pathway and suppress estrogen resistance of myocytes. This effect has no apparent interference in the normal endocrine.

5. Conclusion

There are several limitations in our study. Firstly, the databases we used can help predict the potential mechanism of GZFLW in treating AM. Nevertheless, most of them are based on the current research results, limiting the discovery of new therapeutic targets. Moreover, the primary detection indexes of the animal experiment only include ERα, GPER1, and EGFR, and this is not enough and comprehensive for data analysis. Finally, the number of animals is relatively small and may lead to study biases. Although network pharmacology strategy could help researchers simplify the complex system of TCM formulas and offer a novel research strategy for studying the mechanism of other TCM formulas, the mechanism remains confirmed by further experiments.

Data Availability

The underlying data used to support the findings of this study are included within the supplementary information file.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This study was supported by grants of the National Natural Science Foundation of China (no. 81973897), Scientific Research Foundation of Hubei Provincial Health Department (no. WJ2019F079), and Yichang Medical and Health Research Foundation (no. A19-301-44).

Supplementary Materials

Table S1: detailed information of GZFLW compounds. Table S2: detailed information on therapeutic targets of AM. Table S3: detailed information of PPI network of 48 overlapped targets. Table S4: KEGG enrichment of overlapped targets from Metascape Database. (Supplementary Materials)

References

O. Yu, R. Schulze-Rath, J. Grafton, K. Hansen, D. Scholes, and S. D. Reed, “Adenomyosis incidence, prevalence and treatment: United States population-based study 2006–2015,” American Journal of Obstetrics and Gynecology, vol. 223, no. 1, p. 94.e10, 2020.
View at: Publisher Site | Google Scholar
M. Kossaï and F. Penault-Llorca, “Role of hormones in common benign uterine lesions: endometrial polyps, leiomyomas, and adenomyosis,” Advances in Experimental Medicine and Biology, vol. 1242, pp. 37–58, 2020.
View at: Publisher Site | Google Scholar
Y. Y. Zeng, Y. G. Guan, and K. Y. Li, “Role of estrogen, estrogen receptors, and aromatase in the pathogenesis of uterine adenomyosis,” Journal of Southern Medical University, vol. 37, no. 3, pp. 383–387, 2017.
View at: Google Scholar
K. D. Baboo, Z. Chen, and X. Zhang, “Progress on medical treatment in the management of adenomyosis,” Journal of Zhejiang University, vol. 48, no. 2, pp. 142–147, 2019.
View at: Google Scholar
A. Suardika and T. G. Astawa Pemayun, “New insights on the pathogenesis of endometriosis and novel non-surgical therapies,” Journal of the Turkish-German Gynecological Association, vol. 19, no. 3, pp. 158–164, 2018.
View at: Publisher Site | Google Scholar
Z. Xiong, L. Lang, X. Gao, W. Xiao, Z. Wang, and L. Zhao, “An integrative urinary metabolomic study of the therapeutic effect of Guizhi Fuling capsule on primary dysmenorrheal rats based 1H NMR and UPLC-MS,” Journal of Pharmaceutical and Biomedical Analysis, vol. 164, pp. 750–758, 2019.
View at: Publisher Site | Google Scholar
W. Zheng, M. Li, Y. Wang et al., “Guizhi fuling capsule exhibits antidysmenorrhea activity by inhibition of cyclooxygenase activity,” Evidence-based Complementary and Alternative Medicine: eCAM, vol. 2020, Article ID 8607931, 10 pages, 2020.
View at: Publisher Site | Google Scholar
S. Hu, “Effect of Guizhi Fuling pill combined with Marvelon on patients of qi stagnation and blood stasis type of Adenomyosis,” Hubei University of Traditional Chinese Medicine, vol. 2015, 2015.
View at: Google Scholar
H. Wang, G. Zhou, M. Zhuang, W. Wang, and X. Fu, “Utilizing network pharmacology and molecular docking to explore the underlying mechanism of Guizhi Fuling Wan in treating endometriosis,” PeerJ, vol. 9, p. e11087, 2021.
View at: Publisher Site | Google Scholar
T. Bahougne, E. Angelopoulou, N. Jeandidier, and V. Simonneaux, “Individual evaluation of luteinizing hormone in aged C57BL/6 J female mice,” GeroScience, vol. 42, no. 1, pp. 323–331, 2020.
View at: Publisher Site | Google Scholar
J. He and G. H. Hu, “Progress on TCM treatment for dysmenorrhea due to endometriosis (adenomyosis),” Chin Arch of Traditional Chinese Medicine, vol. 35, no. 3, pp. 692–695, 2017.
View at: Google Scholar
L. Du, D.-H. Du, B. Chen, Y. Ding, T. Zhang, and W. Xiao, “Anti-inflammatory activity of sanjie zhentong capsule assessed by network pharmacology analysis of adenomyosis treatment,” Drug Design, Development and Therapy, vol. 14, pp. 697–713, 2020.
View at: Publisher Site | Google Scholar
H. Feng, L. Zang, and D. Z. Zhang, “Clinical study of modified shaofu zhuyu decoction in treating adenomyosis dysmenorrheal,” Journal of Nanjing University Traditional Chinese Medicine, vol. 32, no. 2, pp. 114–117, 2016.
View at: Google Scholar
Y. H. Wang, “Study on regulating the expression of ER α, ER β, GPER, PR-A and PR-B in mice with adenomyosis by clearing dampness and removing blood stasis,” Chengdu University Traditional Chinese Medicine, 2015.
View at: Google Scholar
G. Q. Wang, Y. Gao, F. Liu, R. L. Wei, and Y. M. Xie, “Post-marketing surveillance on Guizhi Fuling Jiaonang based on literature review,” China Journal of Chinese Materia Medica, vol. 43, no. 4, pp. 820–832, 2018.
View at: Google Scholar
M. Galvankar, N. Singh, and D. Modi, “Estrogen is essential but not sufficient to induce endometriosis,” Journal of biosciences, vol. 42, no. 2, pp. 251–263, 2017.
View at: Publisher Site | Google Scholar
P. Chen, D.-B. Wang, and Y.-M. Liang, “Evaluation of estrogen in endometriosis patients: regulation of GATA-3 in endometrial cells and effects on Th2 cytokines,” Journal of Obstetrics and Gynaecology Research, vol. 42, no. 6, pp. 669–677, 2016.
View at: Publisher Site | Google Scholar
C. A. Piccinato, R. M. Neme, N. Torres et al., “Effects of steroid hormone on estrogen sulfotransferase and on steroid sulfatase expression in endometriosis tissue and stromal cells,” The Journal of Steroid Biochemistry and Molecular Biology, vol. 158, pp. 117–126, 2016.
View at: Publisher Site | Google Scholar
Z. Jin, H. Liu, and C. Xu, “Estrogen degrades Scribble in endometrial epithelial cells through E3 ubiquitin ligase HECW1 in the development of diffuse adenomyosis†,” Biology of Reproduction, vol. 102, no. 2, pp. 376–387, 2020.
View at: Google Scholar
R. Hu, G.-Q. Peng, D.-Y. Ban, C. Zhang, X.-Q. Zhang, and Y.-P. Li, “High-expression of neuropilin 1 correlates to estrogen-induced epithelial-mesenchymal transition of endometrial cells in adenomyosis,” Reproductive Sciences, vol. 27, no. 1, pp. 395–403, 2020.
View at: Publisher Site | Google Scholar
Y. Chen, J. Wang, D.-Y. Hong et al., “βBaicalein has protective effects on the 17β-estradiol-induced transformation of breast epithelial cells,” Oncotarget, vol. 8, no. 6, pp. 10470–10484, 2017.
View at: Publisher Site | Google Scholar
J. Nie and X. Liu, “Quercetin alleviates generalized hyperalgesia in mice with induced adenomyosis,” Molecular Medicine Reports, vol. 16, no. 4, pp. 5370–5376, 2017.
View at: Publisher Site | Google Scholar
W. Xu, Y. Song, K. Li, B. Zhang, and X. Zhu, “Quercetin inhibits adenomyosis by attenuating cell proliferation, migration and invasion of ectopic endometrial stromal cells,” Drug Design, Development and Therapy, vol. 14, pp. 3815–3826, 2020.
View at: Publisher Site | Google Scholar
M. Ilhan, Z. Ali, I. A. Khan, H. Taştan, and E. Küpeli Akkol, “The regression of endometriosis with glycosylated flavonoids isolated from Melilotus officinalis (L.) Pall. in an endometriosis rat model,” Taiwanese Journal of Obstetrics & Gynecology, vol. 59, no. 2, pp. 211–219, 2020.
View at: Publisher Site | Google Scholar
M. Ilhan, Z. Ali, I. A. Khan, H. Taştan, and E. Küpeli Akkol, “Promising activity of Anthemis austriaca Jacq. on the endometriosis rat model and isolation of its active constituents,” Saudi Pharmaceutical Journal, vol. 27, no. 6, pp. 889–899, 2019.
View at: Publisher Site | Google Scholar
S. Park, W. Lim, F. W. Bazer, K.-Y. Whang, and G. Song, “Quercetin inhibits proliferation of endometriosis regulating cyclin D1 and its target microRNAs in vitro and in vivo,” The Journal of Nutritional Biochemistry, vol. 63, pp. 87–100, 2019.
View at: Publisher Site | Google Scholar
L. Schröder, P. Marahrens, and J. G. Koch, “Effects of green tea, matcha tea and their components epigallocatechin gallate and quercetin on MCF-7 and MDA-MB-231 breast carcinoma cells,” Oncology, vol. 41, no. 1, pp. 387–396, 2019.
View at: Google Scholar
G.-A. Lee, K.-C. Choi, and K.-A. Hwang, “Kaempferol, a phytoestrogen, suppressed triclosan-induced epithelial-mesenchymal transition and metastatic-related behaviors of MCF-7 breast cancer cells,” Environmental Toxicology and Pharmacology, vol. 49, pp. 48–57, 2017.
View at: Publisher Site | Google Scholar
Z. Xu, D. Zhao, X. Zheng, B. Huang, X. Xia, and X. Pan, “Quercetin exerts bidirectional regulation effects on the efficacy of tamoxifen in estrogen receptor‐positive breast cancer therapy: an in vitro study,” Environmental Toxicology, vol. 35, no. 11, pp. 1179–1193, 2020.
View at: Publisher Site | Google Scholar
Y. G. Guan, J. B. Liao, K. Y. Li et al., “Potential mechanisms of an antiadenomyosis Chinese herbal formula shaoyao-gancao decoction in primary cell culture model,” Evidence-based Complementary and Alternative Medicine: eCAM, vol. 2014, Article ID 982913, 11 pages, 2014.
View at: Publisher Site | Google Scholar
M. Szubert, E. Koziróg, O. Olszak, K. Krygier-Kurz, J. Kazmierczak, and J. Wilczynski, “Adenomyosis and infertility-review of medical and surgical approaches,” International Journal of Environmental Research and Public Health, vol. 18, no. 3, 2021.
View at: Google Scholar
T. Singtripop, T. Mori, S. Sakamoto, S. Sassa, M. K. Park, and S. Kawashima, “Suppression of the development of uterine adenomyosis by danazol treatment in mice,” Life Sciences, vol. 51, no. 14, pp. 1119–1125, 1992.
View at: Publisher Site | Google Scholar
J. F. Nelson, L. S. Felicio, P. K. Randall, C. Sims, and C. E. Finch, “A longitudinal study of estrous cyclicity in aging C57bl/6J mice: I. Cycle frequency, length and vaginal Cytology1,” Biology of Reproduction, vol. 27, no. 2, pp. 327–339, 1982.
View at: Publisher Site | Google Scholar
Y.-M. Wang, M.-J. Wu, Y.-H. Lin, and J. Chen, “Association of epidermal growth factor receptor (EGFR) gene polymorphisms with endometriosis,” Medicine, vol. 98, no. 17, p. e15137, 2019.
View at: Publisher Site | Google Scholar
H. Zhan, B. Peng, J. Ma et al., “Epidermal growth factor promotes stromal cells migration and invasion via up-regulation of hyaluronate synthase 2 and hyaluronan in endometriosis,” Fertility and Sterility, vol. 114, no. 4, pp. 888–898, 2020.
View at: Publisher Site | Google Scholar
L. Liu, N. Luo, J. Guo, Y. Xie, L. Chen, and Z. Cheng, “Berberine inhibits growth and inflammatory invasive phenotypes of ectopic stromal cells: imply the possible treatment of adenomyosis,” Journal of Pharmacological Sciences, vol. 137, no. 1, pp. 5–11, 2018.
View at: Publisher Site | Google Scholar
T. Wang, Y. G. Cao, and Y. Liu, “Study on the mechanism of Guizhi Fuling Wan combined with Buzhong Yiqi Wan on benign prostatic hyperplasia in ovariectomized rats,” Chin ExpTradit Med Form, vol. 16, no. 17, pp. 154–157, 2010.
View at: Google Scholar
P. Wang, S. M. Henning, C. E. Magyar, Y. Elshimali, D. Heber, and J. V. Vadgama, “Green tea and quercetin sensitize PC-3 xenograft prostate tumors to docetaxel chemotherapy,” Journal of Experimental & Clinical Cancer Research, vol. 35, no. 1, p. 73, 2016.
View at: Publisher Site | Google Scholar
A. Carazo Fernández, T. Smutny, L. Hyrsová, K. Berka, and P. Pavek, “Chrysin, baicalein and galangin are indirect activators of the human constitutive androstane receptor (CAR),” Toxicology letters, vol. 233, no. 2, pp. 68–77, 2015.
View at: Publisher Site | Google Scholar
N. Oi, H. Chen, M. Ok Kim, R. A. Lubet, A. M. Bode, and Z. Dong, “Taxifolin suppresses UV-induced skin carcinogenesis by targeting EGFR and PI3K,” Cancer Prevention Research, vol. 5, no. 9, pp. 1103–1114, 2012.
View at: Publisher Site | Google Scholar
K. N. Khan, A. Fujishita, A. Koshiba et al., “Biological differences between focal and diffuse adenomyosis and response to hormonal treatment,” Reproductive BioMedicine Online, vol. 38, no. 4, pp. 634–646, 2019.
View at: Publisher Site | Google Scholar
D. G. Hong, J. Y. Park, G. O. Chong et al., “Transmembrane G protein-coupled receptor 30 gene polymorphisms and uterine adenomyosis in Korean women,” Gynecological Endocrinology, vol. 35, no. 6, pp. 498–501, 2019.
View at: Publisher Site | Google Scholar
J. J. Li, H. Duan, S Wang et al., “Expression pattern of G-protein-coupled estrogen receptor in myometrium of uteri with and without adenomyosis,” BioMed Research International, vol. 2017, Article ID 5974693, 6 pages, 2017.
View at: Publisher Site | Google Scholar
K. Chatterjee, S. Jana, P. DasMahapatra, and S. Swarnakar, “EGFR‐mediated matrix metalloproteinase‐7 up‐regulation promotes epithelial‐mesenchymal transition via ERK1‐AP1 axis during ovarian endometriosis progression,” The FASEB Journal, vol. 32, no. 8, pp. 4560–4572, 2018.
View at: Publisher Site | Google Scholar
Y. Kishi, K. Shimada, T. Fujii et al., “Phenotypic characterization of adenomyosis occurring at the inner and outer myometrium,” PLoS One, vol. 12, no. 12, p. e0189522, 2017.
View at: Publisher Site | Google Scholar
J. García-Solares, J. Donnez, O. Donnez, and M.-M. Dolmans, “Pathogenesis of uterine adenomyosis: invagination or metaplasia?” Fertility and Sterility, vol. 109, no. 3, pp. 371–379, 2018.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2021 Haoxian Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies