Ziqing Gao1
Rui Mi2
Zhaoxi Cheng2
Xiaofeng Li3
Huawu Zeng2
Gaosong Wu1
Jing Zhao1*
Weidong Zhang1,2*
Ji Ye2*- 1Institute of Interdisciplinary Integrative Medicine Research, Shanghai University of Traditional Chinese Medicine, Shanghai, China
- 2School of Pharmacy, Second Military Medical University, Shanghai, China
- 3School of Pharmacy, Fudan University, Shanghai, China
Hong-Hua-Xiao-Yao Tablet (HHXYT) is a traditional Chinese medicine (TCM) formula that has been approved for the treatment of mammary gland hyperplasia (MGH), but its mechanism of action is unclear. In this study, a strategy that integrated metabolomics and network pharmacology was applied to systemically reveal the mechanism of HHXYT in the treatment of MGH. Our pharmacodynamic study indicated that the proliferation of mammary gland was inhibited in rats, and serum-level disorder of estradiol and progesterone was reversed after HHXYT treatment. 54 compounds absorbed in rat plasma were identified after administration of HHXYT. The serum metabolome revealed 58 endogenous differential metabolites, of which 31% were steroid lipids metabolites, with steroid hormone biosynthesis being the most significant metabolic module. 7 targets, 6 herbs, and 17 ingredients were found to play key roles in HHXYT’s treatment of MGH. 3 of the 7 key targets (CYP11A1, HSD3B2, and CYP17A1) were directly involved in androgen synthesis, while 2 targets (AR and ESR1) were receptors for the direct action of androgens and estrogens. Molecular docking was utilized to confirm the bindings between the 5 targets and their corresponding compounds. In an in vitro test, HHXYT (50 µg/ml) and its ingredient formononetin (3.2, 6.3, and 12.5 µM) were found to significantly reduce the increase of testosterone level induced by dexamethasone (10 µM) in thecal cells. In summary, this study illustrated that the mechanism of HHXYT’s treatment of MGH was to regulate hormone disorder. HHXYT could reduce estrogen-stimulated hyperplasia by inhibiting the production of its precursor androgen.
1 Introduction
The most common breast disease, mammary gland hyperplasia (MGH), is a proliferative lesion of mammary tissue accompanying by the lumps and pain symptoms of the breast. It has overtly interfered with the patient’s daily life (You et al., 2017; Li et al., 2018). In the past few years, it has been increasingly reported that the incidence of breast cancer in MGH patients is remarkably higher than that in average person (Hartmann et al., 2005), since the pathogenesis of MGH is similar to that of breast cancer, both of which are closely related to endocrine and hormone level disorders (Coussens and Pollard, 2011; Arendt and Kuperwasser, 2015). Hormone drugs are currently applied to the therapy of MGH to regulate internal hormone levels. However, due to its obvious side effects, such as gastrointestinal and cardiovascular problems, the safety and effectiveness of this approach are not optimistic (Wibowo et al., 2016).
As a traditional Chinese medicine (TCM) formula on MGH treatment in the clinic for years, Hong-Hua-Xiao-Yao Tablet (HHXYT) was composed of nine herbal medicines, including Carthamus tinctorius L. (Honghua), Angelica sinensis (Oliv.) Diels (Danggui), Paeonia lactiflora Pall. (Baishao), Gleditsia sinensis Lam. (Zaojiaoci), Atractylodes macrocephala Koidz. (Baizhu), Wolfiporia cocos (F.A. Wolf) Ryvarden & Gilb (Fuling), Mentha canadensis L. (Bohe), Bupleurum marginatum Wall. ex DC. (Zhuye Chaihu), and Glycyrrhiza uralensis Fisch. ex DC. (Gancao). Chemical profiling of HHXYT has been previously elucidated by high-performance liquid chromatography combined with mass spectrometry technique, and 55 constituents were characterized. Moreover, in our earlier study, 14 components that were identified as potential HHXYT quality markers were quantified using high-performance liquid chromatography-tandem triple quadrupole mass spectrometry (HPLC-QQQ-MS) technique (Mi et al., 2020). In clinical trials, HHXYT was found to be effective in the treatment of perimenopausal syndrome, ovarian syndrome, menstrual disorders, as well as mammary gland hyperplasia (Wang, 2012; Pei et al., 2016; Chen et al., 2020). However, the mechanism of its therapy was still unelucidated.
Metabolomics could integrally focus on alterations of endogenous metabolites to monitor and analyze metabolic disturbances for elucidating therapeutic mechanisms (Nicholson and Wilson, 2003). Moreover, compared with other omics techniques, it is more in line with the change of disease phenotype (Newgard, 2017). Network pharmacology performs network analysis of complex biological systems by constructing network models (Hopkins, 2007; Li and Zhang, 2013). The combined analysis of metabolomics and network pharmacology can characterize the overall changes of biological systems, which is in accordance with the holistic therapeutic concept of TCM. This strategy has been increasingly applied in the study of TCM pharmacology and has contributed to elucidating the mechanism of multi-pathway and multi-target of TCM (Xu et al., 2018; Wu, Zhang and Li, 2019; Zhao et al., 2019).
In this study, the integrated metabolomics and network pharmacology strategy were used to reveal the mechanism of HHXYT in the treatment of MGH on a systemic level. Rat model of MGH was established and the metabolomics approach based on the HPLC-QTOF-MS/MS analysis was operated to search for serum endogenous differential metabolites. The ingredients and putative targets of HHXYT were then analyzed using the network pharmacology method. Next, an integrated analysis of differential metabolites and targets identified HHXYT’s key targets and substances against MGH. Ultimately, the cell-based assay on the synthesis of sex hormones was performed to validate HHXYT’s effects in vitro (Figure 1).
FIGURE 1. Integrated workflow for exploring the potential mechanisms of HHXYT against MGH.
2 Materials and Methods
2.1 Chemicals and Reagents
Reagents of the LC-MS purity grade were acquired from Fisher Scientific (Fair Lawn, NJ, United Staes), including formic acid (purity >98%), acetonitrile, and methanol. Estradiol benzoate and progesterone injection were purchased from Zhejiang Xianju Pharmaceutical Co., Ltd. (Taizhou, China). ELISA kits of estradiol (E2), testosterone (T), and progesterone (P) were obtained from Boster Biological Technology Co., Ltd. (Shanghai, China). McCoy’s 5A medium was provided by Gibco® (Grand Island, NY, United States). Follicle-stimulating hormone (FSH) and luteinizing hormone (LH) were acquired from Shanghai Linc-Bio Science Co. Ltd. (Shanghai, China). Dexamethasone (DEX) and formononetin (FT) were supplied by Chengdu Must Biotechnology Co., Ltd. (Chengdu, China). The purity of FT was above 98%. Cell Counting Kit-8 (CCK-8) was purchased from Shanghai Beyotime Biotechnology Co., Ltd. (Shanghai, China). HHXYT (batch number: 190106), the same batch named as “PZ-5” in our previous study (Mi et al., 2020), was supplied by Puzheng Pharmaceutical Co., Ltd. (Jiangxi, China). The contents of fourteen quality marker compounds of HHXYT were quantified by HPLC-QQQ-MS (Table 1). The chromatograms and chemical structure were shown in Supplementary Figure S1.
TABLE 1. The contents of fourteen quality marker components of HHXYT (batch number: 190106, mg/g, Mean ± SD, n = 3).
2.2 Animal Experiment and Sample Preparation
All sexually mature and non-pregnant female Sprague-Dawley rats were obtained from Shanghai Xipuer-Bikai Experimental Animal Co. Ltd., Shanghai, China (License No.: SCXK-Shanghai-2018-0006). They were bred in a stabilized environment for 1 week with room temperature 25 ± 2°C, relative humidity 55 ± 5%, 12 h each of light and dark, and free access to drinking water and feed. All procedures were approved by the Animal Care and Use Committee of Second Military Medical University and abided by the Guide for the Care and Use of Laboratory Animals.
2.2.1 Prototypes of HHXYT and Related Metabolites in Plasma
The animal experiment used to study the chemicals absorbed in rat plasma was previously described (Mi et al., 2020). Briefly, the HHXYT sample was suspended in 0.5% carboxymethylcellulose sodium (CMC-Na) and administered orally to twenty-five female rats at five time-set points of 0.5, 1, 2, 4, and 6 h, at a dose of 8.82 g/kg (the optimal dose among 4 different doses in the amount and intensity of target compounds, Supplementary Figure S2), respectively (n = 5). Five other rats received oral administration of an equivalent volume of 0.5% CMC-Na. The blood samples were gathered from the hepatic portal vein, and all the samples were centrifuged at 3,500 rpm for 10 min at 4°C. A four-fold volume of methanol was used to extract the supernatants by the protein precipitation approach. After centrifuging at 13,000 rpm for 10 min at 4°C, the supernatants were blow-dried and reconstituted with 100 µL of 75% methanol for further analysis.
2.2.2 Pharmacodynamics and Metabolomics Study
Twenty-four female Rats were haphazardly and equally divided into four groups: sham operation group (Sham), the mammary gland hyperplasia model group (MGH), low dose group of HHXYT (HHXYT-L), and high dose group of HHXYT (HHXYT-H). The MGH model was achieved by injecting intramuscularly with estradiol benzoate (0.5 mg/kg/d) for 25 days and then with progesterone (5 mg/kg/d) for 5 days (Wang et al., 2011). Meanwhile, the Sham group was injected with saline injection (0.5 mg/kg/d). After modeling, the HHXYT-L group and HHXYT-H group were administrated HHXYT orally with suspending in 0.5% CMC-Na aqueous solution, at a dosage of 0.5 g/kg/d (clinical equivalent dose) and 4.5 g/kg/d for 3 weeks, respectively. The Sham group and MGH group were given oral administration of 0.5% CMC-Na.
Rats were injected intraperitoneally with 10% chloral hydrate for anesthesia after 3 weeks of administration, and the diameter of the second pair of nipples was measured. The blood samples were obtained from the abdominal aorta, centrifuged at 3,500 rpm at 4°C for 15 min. The serum was separated from the upper layer and stored at −80°C. The left mammary gland of the second pair of breasts was collected and fixed in tissue fixative for 24 h, then sectioned in paraffin and stained with Hematoxylin-Eosin (HE) for optical microscopy. The contents of estradiol (E2) and progesterone (P) in serum were determined by ELISA kits.
300 µL of methanol was put in 100 µL of serum, which was eddied for 1 min to ensure complete mingling before centrifugation at 14,000 rpm for 15 min at 4°C. The supernatant was aspirated and stored overnight at −80°C, followed by centrifugation under the same conditions for 15 min to remove proteins once again. The final supernatant was pipetted into the inner tube of vials and allowed to stand for LC-MC analysis. The quality control samples (QC) were prepared by pooling 20 µL of 6 random samples (containing all groups) from all samples.
2.3 Instrument Parameters
2.3.1 UPLC-QTOF-MS/MS for Identifying Prototypes of HHXYT and Related Metabolites in Plasma
The compounds absorbed in rat plasma were detected using a Waters Acquity UPLC I-class system (Waters Corp., Milford, MA, United States) and Q-tof Synapt G2-Si mass spectrometry system (Waters, Manchester, United Kingdom). A Waters Acquity UPLC BEH reverse-phase C18 analytical column (1.7 µm, 3.0 × 100 mm, Waters) was used to separate the compounds. The column temperature was set at 35°C and the injection volume was 5 µL. The mobile phase was made up of 0.1% formic acid in water (v/v) (A) and acetonitrile (B). At a flow rate of 0.4 ml/min, the elution gradient was processed in 37 min. The base peak chromatograms of both positive and negative ion modes were acquired under multiple detection modes of MSE and Fast DDA. The parameters were set following the previous report (Mi et al., 2020).
2.3.2 HPLC-QTOF-MS/MS for Metabolomics Analysis
LC system was considered to apply with a Shimadzu HPLC system (Nexera XR LC-20AD, Japan) and use an ACQUITY UPLC® HSS T3 (2.1 × 100 mm, 1.8 µm, Waters, Milford, MA, United States) as the chromatographic column. Mass spectrometric detection was carried out on a SCIEX Triple TOF 5600+, using the electrospray ion (ESI) source and information-dependent acquisition (IDA) mode. The column temperature was set at 40°C and the injection volume was 2 µL. The mobile phase was made up of mass spectrometric grade 0.1% formic acid in water (v/v) (A) and acetonitrile (B). The calibration deliver system (CDS) was applied to calibrate the mass spectrometric acquisitions automatically. Additional parameter settings are provided in the Supplementary Materials.
2.4 Data Processing for Metabolomics Results
The raw data files were imported into the Progenesis QI 2.0 software (Waters, Milford, MA, United States) for data preprocessing, including alignment, peak picking, and compound identification. StatTarget was utilized to remove system errors as well as signal drift (Luan et al., 2018). The data was then imported into SICMA 14.1 (Umetric, Umeå, Sweden) for multivariable statistical analysis. Variates with p-values < 0.05, fold change (FC) ≥1.2 or FC ≤ 0.8, and VIP >1 (variable importance in the projection), were considered as endogenous differential metabolites, and further identified with online databases of HMDB (Wishart et al., 2018), LIPID MAPS (Liebisch et al., 2020) and KEGG (Kanehisa et al., 2017).
2.5 Identification and Analysis of Metabolic Modules
We applied the diffusion kernel algorithm to detect metabolites that locate at the proximity of the endogenous differential metabolites in the genome-scale metabolic network of rats (Luan et al., 2018). The diffusion parameter was set as 0.1. All endogenous differential metabolites of the MGH, HHXYT-L, and HHXYT-H groups were input as seed nodes to perform the diffusion kernel algorithm. The nodes with high diffusion scores were taken as compounds close to the endogenous differential metabolites in the genome-scale metabolic network. Then, the metabolic network affected by mammary glands hyperplasia and HHXYT treatment was constructed by extracting these nodes and links between them from the background network.
2.6 Putative Targets for HHXYT’s Plasma Absorbed Compounds
For prototypes and their metabolites in the plasma of rats after intragastric administration of HHXYT, we collected or predicted their putative targets from 3 databases and web servers: ETCM (Encyclopaedia of Traditional Chinese Medicine) (Xu et al., 2019), HERB (High-throughput Experiment- and Reference-guided dataBase of traditional Chinese medicine) (Fang et al., 2021), and TargetNet (Yao et al., 2016). Target information in HERB was collected by manual reference mining and from multiple TCM databases, including SymMap (Wu, Zhang, Yang et al., 2019), TCMID (Xue et al., 2013), TCMSP (Ru et al., 2014), and TCM-ID (Chen et al., 2006). The threshold of confidence score was set as 0.9 to filter the compound-protein interactions predicted by TargetNet.
2.7 Collection of MGH Associated Genes and Functional Enrichment Analysis
Disease genes of MGH were collected from the GeneCards database (Stelzer et al., 2016). Inputting the keyword “mammary gland hyperplasia” generated 4,251 results with relevance scores in the area [0.53, 248.35]. We set the threshold of relevance score as 20 to filter the results and obtained 420 top-ranked genes with high relevance to MGH, taking about 10% of the results. These 420 genes were considered as MGH associated genes for further analysis.
We used the STRING platform to conduct functional enrichment analysis for gene groups under study (Szklarczyk et al., 2019).
2.8 Network Construction, Analysis, and Visualization
The genomic scale protein-protein interaction (PPI) network of human beings was downloaded from the HINT (High-quality INTeractomes) database (Das and Yu, 2012), whose PPI data was integrated from 8 databases (BioGRID, MINT, iRefWeb, DIP, IntAct, HPRD, MIPS, and the PDB) and filtered both systematically and manually to keep high-confidence and -quality interactions. It was used as a background network for the construction of HHXYT’s target network. We used Cytoscape version 3.6.0 to construct, analysis and visualize networks (Shannon et al., 2003). The Cytohubba app of Cytoscape was applied to compute topological measures of nodes (Chin et al., 2014).
Metabolite-gene network was constructed by the Network Analysis module of MetaboAnalyst platform version 5.0 (Chong et al., 2019).
2.9 Molecular Docking
The SDF file for a ligand’s structure was downloaded from the PubChem database (Kim et al., 2021). The PDB file for a protein receptor’s crystal structure was obtained from RCSB Protein Data Bank (PDB) database (Burley et al., 2021). For proteins without known crystal structures, the predicted crystal structures were downloaded from AlphaFold Protein Structure Database (Jumper et al., 2021). The preprocessing and preparing of ligands and receptors were conducted by AutoDock Tools (Tanchuk et al., 2016) following the tutorial and manual (http://vina.scripps.edu/manual.html). Each ligand was docked with corresponding receptors by AutoDock Vina (Trott and Olson, 2010), and a compound-target pair with the docking score of less than −5 kcal/mol was considered as a binding pair. PLIP platform was used to analyze the binding sites of the binding compound-target pair (Adasme et al., 2021). PyMol was applied for visualizing the results of AutoDock Vina and PLIP (Alexander et al., 2011).
2.10 Cell-Based Assay Evaluation of HHXYT In Vitro
Thecal cells were extracted and purified from the ovaries of female rats at 3–4 weeks old (Li and Hearn, 2000). 5 × 103 thecal cells per well were seeded into 96-well plates and set as the control group and experimental groups. The control group was McCoy’s 5A medium containing thecal cell. The experimental groups were based on the control group supplemented with different concentrations of luteinizing hormone (LH), follicle-stimulating hormone (FSH), HHXYT, and its ingredient formononetin (FT). After 36 h of adherent culture at 37°C and 5% CO2, these drugs were supplemented according to different groups and cultured for another 24 h. Then 10 µL of CCK-8 reagent was added to each well in the dark, and the absorbance was measured at a wavelength of 450 nm 2 h later to calculate the cell survival rate (Chen et al., 2021).
4 × 104 Cells were added to each well of 48-well plates containing McCoy’s 5A medium, which consisted of 10% FBS and 1% penicillin-streptomycin. After culturing at 37°C and 5% CO2 for 36 h, different concentrations of LH, dexamethasone (DEX), HHXYT, and FT were added into different groups of wells and cultured for another 24 h (Ong et al., 2019). ELISA kits were used to detect the levels of estradiol (E2) or testosterone (T) in cell media.
3 Results
3.1 Pharmacodynamic Results of HHXYT on MGH Model
To evaluate the therapeutic effect of HHXYT on MGH, we detected variations of nipple diameter and serum sex hormone level. After establishing the MGH model, the nipple diameter grew significantly (p < 0.01, MGH vs. Sham) from 0.858 ± 0.068 mm to 1.230 ± 0.086 mm, compared to the Sham group. Different doses of HHXYT, on the other hand, could significantly (p < 0.01, HHXYT vs. MGH) reduce nipple enlargement (Figure 2A). Meanwhile, in the MGH group, the estradiol level in serum increased prominently (p < 0.01, MGH vs. Sham) while the progesterone level decreased notably (p < 0.01, MGH vs. Sham). After treatment with different doses of HHXYT, the serum level of estradiol and progesterone was reversed, respectively (Figures 2B,C).
FIGURE 2. The effect of HHXYT on pharmacodynamic results of the MGH model. (A) The diameter of nipples. The serum levels of sex hormone estradiol (B) and progesterone (C). n = 6 per group, ##p < 0.01, compared to Sham; *p < 0.05, **p < 0.01, compared to MGH. HE staining results for mammary gland sections in Sham (D), MGH (E), HHXYT-H (F) group (10 × 10: D-1, E-1, and F-1; 10 × 40: D-2, E-2, and F-2).
According to Hematoxylin-Eosin (HE) staining, mammary tissue in the Sham group had no aberrant proliferation, a tiny number of the acinus, nearly no dilatation, and no secretion of mammary ducts and lumens (Figure 2D). The mammary tissue in the MGH group showed distinct proliferation; the volume of mammary lobules was increased; the number of glandular bubbles was also increased and became close to fusion; the acinus and lumens were revealed dilatation and excessive secretions (Figure 2E). However, after HHXYT treatment, mammary gland hyperplasia was prominently alleviated due to a considerable reduction in the account of the volume of mammary lobules, the number of glandular vesicles, and the secretion of lumens (Figure 2F).
3.2 Characterization of Prototypes and Metabolites in HHXYT, and Functional Enrichment Analysis of Putative Targets
It is widely accepted that prototypes and their metabolites in blood after oral TCM intake have a therapeutic effect (Ye et al., 2019). In order to determine the active components of HHXYT, we comprehensively analyzed the prototypes and their metabolites of HHXYT in rat plasma. A total of 54 components were identified, comprising 22 prototypes (Figures 3A,B; Supplementary Table S1) and 32 metabolites (Supplementary Table S2). After HHXYT treatment, three prototype components, named desbenzoylpaeoniflorin, 1-O-β-d–glucopyranosylpaeonisuffrone, and bohenoside A, as well as 32 metabolites, were identified in rat plasma for the first time. These 32 metabolites include 2 monoterpenoids, 8 triterpenes, 11 organic acids, 10 flavonoids, and 1 other metabolite (Figure 3C; Supplementary Table S4-S7).
FIGURE 3. The prototypes and their metabolites in the plasma of rats after intragastric administration of HHXYT. Total ion chromatograms of rat plasma samples after intragastric administration of HHXYT and the distribution of 22 prototypes, positive mode (A) and negative mode (B). (C) Categories of prototypes and metabolites.
Based on 54 active ingredients identified in the rat plasma, we predicted their potential target through network pharmacology (hydroxysafflower yellow A was also added to the analysis because it was reported to be an active ingredient of Carthami Flos (Honghua) (Li et al., 2021)). Using 3 databases and web servers (ETCM, HERB, and TargetNet), 373 distinct protein targets were predicted as HHXYT’s putative targets regulated by the 55 active ingredients (Supplementary Table S3).
The STRING platform was utilized to analyze the 373 targets for functional enrichment. Setting FDR<10−4 as the threshold, we identified 93 statistically enriched KEGG pathways. The top enriched pathways were presented in Table 2, and these pathways were mainly involved in the biological processes of the endocrine, immune, and nervous system.
TABLE 2. Enriched KEGG pathways of HHXYT’s targets. The pathways with bold type are associated with metabolic modules affected by HHXYT.
Then, based on the KEGG pathways in Table 2, we performed enrichment analysis for targets of HHXYT’s herbs and prototypes. When conducting enrichment analysis, we combined targets of metabolites derived from its prototype with targets of the prototype. The threshold for statistical significance was set as FDR<10−4. As shown in Figure 4, targets of 6 herbs and 9 prototypes from HHXYT were enriched in at least 3 of the 20 KEGG pathways in Table 2, respectively. These results suggested that the 6 herbs and 9 prototypes could play key roles in HHXYT’s treatment of MGH.
FIGURE 4. KEGG pathway enrichment analysis for the targets of herbs and prototypes of HHXYT. Targets of 6 herbs (A) and 9 prototypes (B) from HHXYT were enriched in at least 3 of the top 20 KEGG pathways regulated by HHXYT, respectively.
3.3 HHXYT Alters the Serum Metabolome
On the basis of research on prototypes and their metabolites of HHXYT, the metabolomics analysis was used to further identify the serum endogenous differential metabolites of rats intervened by HHXYT. The PLS-DA model was applied to characterize the metabolic disturbances. As shown in Figures 5A,B, the serum samples of Sham, MGH, HHXYT-L, and HHXYT-H groups were distinctly separated in both positive and negative modes. Based on the criteria of p-values < 0.05, FC ≥ 1.2 or FC ≤ 0.8, and VIP >1, 58 endogenous differential metabolites in serum were screened and identified, including 29 metabolites in the Sham/MGH group, 19 metabolites in the MGH/HHXYT-L group, and 18 metabolites in the MGH/HHXYT-H group (Figure 5C). There were 5 same endogenous differential metabolites between HHXYT-H and MGH (Figure 5D, red compound names), as well as 2 overlaps between HHXYT-L and MGH (blue compound names). Furthermore, the regulatory trends of HHXYT-H on the 5 same differential metabolites were all opposite to that of MGH. Meanwhile, only one differential metabolite had been oppositely regulated between HHXYT-L and MGH groups. These results indicated that the HHXYT-H group had more pronounced alterations in metabolites endogenous on treating MGH. The majority of these endogenous differential metabolites are lipid metabolites, accounting for 68%, including 31% steroids lipid (12% bile acids, 8% pregnane steroids, 7% hydroxysteroids, etc.), 8% fatty acids, 7% fatty acid esters, 8% prenol lipids, 8% glycerophospholipids and so on (Figure 5E). The above results suggested that HHXYT could prominently regulate lipid metabolism, especially hormone disorder in rats with hyperplasia of mammary glands.
FIGURE 5. Endogenous differential metabolites were identified from the metabolomics experiment. PLS-DA score plots of Sham, MGH HHXYT-L, and HHXYT-H, in positive mode (A) and negative mode (B). (C) The intersection of endogenous differential metabolites in the 3 groups. (D) Heatmap of 58 endogenous differential metabolites. n = 6, #p < 0.05, ##p < 0.01, compared to Sham; *p < 0.05, **p < 0.01, compared to MGH. The red compound name means the common differential metabolites between HHXYT-H and MGH; Blue: the common differential metabolites of HHXYT-L/MGH; Green: the common differential metabolites of HHXYT-H/HHXYT-L. (E) Categories of endogenous differential metabolites. 31% of steroids were marked with the orange border.
3.4 Metabolic Modules Affected by the Treatment of HHXYT
To reveal the underlying significance of metabolites variation in vivo, we proceeded with metabolic module establishment and pathway analysis of endogenous differential metabolites. Inputting all the endogenous differential metabolites as seed nodes, we ran the diffusion kernel algorithm to score all nodes in the genome-scale metabolic network of rats. The metabolic network influenced by MGH and HHXYT therapy was then constructed by extracting the nodes with diffusion scores at least 0.05 and the links between them (Figure 6). This network consists of 19 isolated modules, in which only the first three large modules include endogenous differential metabolites affected by both MGH and HHXYT treatment, suggesting that they could be key metabolic modules by which HHXYT exerts its therapeutic effect on MGH.
FIGURE 6. Metabolic modules affected by HHXYT treatment and their functions. Circle nodes are endogenous differential metabolites identified in metabolomics, and hexagonal nodes are metabolites close to the endogenous differential metabolites in the genomic metabolic network which were identified by the diffusion kernel algorithm. Yellow circles are endogenous differential metabolites in MGH, while circles with blue and red colors or borders were in HHXYT-L and HHXYT-H, respectively.
We searched the KEGG database (Kanehisa et al., 2017) to annotate the first 3 modules in Figure 6. It was found that metabolites in the first module were mainly involved in pathways of steroid hormone biosynthesis, cortisol synthesis and secretion, and ovarian steroidogenesis; the compounds in the second module participated in primary bile acid biosynthesis and bile secretion; the third module was associated with glycerophospholipid metabolism and biosynthesis of unsaturated fatty acids. Notably, all the 3 pathways associated with the first metabolic module and 1 pathway of the second metabolic module, were also found in the pathways listed in Table 2, which were enriched by 373 HHXYT’s predicted targets. It demonstrated that the pathways involved by HHXYT’s putative targets discovered by network pharmacology were consistent with the analyzed regulation of metabolic disturbances in rat serum by metabolomics.
3.5 Network Pharmacological Analysis for HHXYT’s Putative Targets and Metabolic Modules
The first step was to conduct a network analysis of HHXYT’s potential targets. We calculated the intersections of HHXYT’s 373 target genes with the 420 diseases genes for MGH (Supplementary Table S4; Figure 7A). It could be seen that one-fourth of HHXYT’s targets were disease genes. Then we constructed the protein-protein interaction (PPI) network between the targets based on the human genomic-scale PPI network from the HINT database. As shown in Figure 7B, 114 of the 373 targets interacted with each other to form HHXYT’s target network, with 53 of them being disease genes associated with MGH.
FIGURE 7. Network pharmacological analysis for HHXYT’s putative targets and integrated analysis of metabolic modules. (A) The intersection between HHXYT’s putative targets and MGH disease genes from GeneCards. (B) HHXYT’s putative target network, where the nodes with red borders are MGH disease genes and the green nodes are central nodes with the highest global centrality. (C) Global (Closeness, EcCentricity, Radiality) and local centrality measures (Degree) of the four central nodes in the target network. Metabolite-gene networks for the first (D) and second metabolic module (E) affected by HHXYT, respectively. Hexagons were endogenous differential metabolites identified from the metabolomics analysis and circles were HHXYT’s target and genes. Yellow represented in MGH disease conditions, while red represented in HHXYT-H. Pink circles were the hub genes of the network.
The centrality measures for nodes in the target network were computed using the Cytohubba App of the Cytoscape software. The centrality of nodes was then ranked applying 3 global parameters based on shortest paths (Closeness, EcCentricity, and Radiality). As shown in Figure 7C, the 4 nodes with the highest global centrality were AR, ESR1, AKT1, and JUN, suggesting that they were central nodes whose impacts were easy to reach out into the whole target network.
The Network Analysis module of the MetaboAnalyst platform version 5.0 was then utilized to construct metabolite-gene networks corresponding to the metabolic modules affected by HHXYT. We input the endogenous differential metabolites in the first and second modules of Figure 6, together with HHXYT’s targets enriched in the corresponding pathways of the modules, respectively. Then we got the metabolite-gene interaction modules for the first and second metabolite module affected by HHXYT (Figures 7D,E). As shown in Figure 7D, the first module included 3 hub genes, i.e., CYP11A1, HSD3B2, and CYP17A1, which linked to all the metabolites in the network, implying that they performed a key role in HHXYT’s treatment to MGH.
The key targets of HHXYT for treating MGH were then obtained by combining the 4 central genes of the target network and the 3 hub genes in the first metabolite-gene module. We integrated the 7 key targets with the key herbs, prototypes (Figure 4), and their metabolites on HHXYT’s treatment to MGH, to construct the herb-ingredient-target network as shown in Figure 8A. There were 7 targets, 6 herbs, and 17 ingredients in this network. 8 of 17 ingredients were HHXYT’s prototypes absorbed in plasma, and 9 were metabolites derived from them. These herbs, ingredients, and targets constitute key substances for HHXYT to exert therapeutic effects on MGH.
FIGURE 8. The key substances for HHXYT to exert therapeutic effects on MGH. (A) Herb-ingredient-target relationships for key targets, herbs, and ingredients of HHXYT. The integrated name of prototypes and their metabolite ID means that the metabolite is derived from the corresponding prototypes. (B–F) Molecular docking of HHXYT’s five key targets with corresponding compounds. The binding poses and binding sites of formononetin complexed with CYP17A1 (B) and HSD3B2 (C); The binding poses and binding sites of glycyrrhetinic acid complexed with AR (D), CYP11A1 (E), and ER (F).
At last, to verify the binding between 5 of HHXYT’s key targets and corresponding prototypes, we utilized AutoDock Vina software to conduct molecular docking. The crystal structures of 4 targets, i.e., AR (PDB ID:2q7i), CYP11A1 (PDB ID:3n9y), CYP17A1 (PDB ID:6ciz), ESR1 (PDB ID: 2qgw), which have higher resolution, were downloaded from RCSB Protein Data Bank. The predicted structure of HSD3B2 (ID: AF-P26439-F1) was downloaded from AlphaFold Protein Structure Database. The SDF files for chemical structures of the compounds (ferulic acid, formononetin, glycyrrhetinic acid, glycyrrhizic acid, isoliquiritigenin, liquiritigenin, and protocatechuic acid) were downloaded from the PubChem database. The docking scores between 14 pairs of compound-target are less than −5 kcal/mol (Table 3), suggesting that the corresponding compounds and targets have a higher binding affinity. Binding poses and binding sites for 5 pairs of compound-target interactions were shown in Figures 8B–F.
TABLE 3. The docking scores between HHXYT’s key targets and compounds (kcal/mol).
3.6 Effects of HHXYT and its Ingredient on the Testosterone Synthesis in Thecal Cells
Firstly, the cell-based assay was performed to demonstrate that the extracted and purified cells in vitro were thecal cells. After LH (50 ng/ml) was added, the cell proliferation rate (maximum ratio 147%) and testosterone (T) level (from 1965.15 ± 96.65 to 2,329.43 ± 87.93 pg/ml) both increased prominently (Figures 9A,B). However, after the addition of FSH, no significant changes in cell quantity were observed, and estradiol (E2) was not detected. It indicated that the purified cells were primarily thecal cells, with no obvious interference from granulosa cells, which were also presented in follicles of the ovary. The cell survival rate results of HHXYT and its ingredient formononetin (FT, Figure 9C) were shown in Supplementary Figure S9.
FIGURE 9. Cell-based assay to evaluate the HHXYT’s effects on testosterone synthesis in thecal cells. (A) The proliferation effect of luteinizing hormone (LH) and follicle-stimulating hormone (FSH) on thecal cells. (B) Effects of LH on the production of testosterone (T) and estradiol (E2) in thecal cells. (C) The chemical structures of formononetin (FT). (D) Effects of dexamethasone (DEX), HHXYT, and FT on the production of testosterone by thecal cell. n = 4, ##p < 0.01, to the Control group; *p < 0.05, **p < 0.01, compared to the DEX (10 µM) group. (E) Pathways for the synthesis of testosterone and estradiol by thecal cells and granulosa cells in follicles. The compounds and targets involved in this study were marked by orange-red colors.
To evaluate the effect of HHXYT and its ingredient formononetin (FT) on the production of testosterone in thecal cells, thecal cells were stimulated with different concentrations of dexamethasone (DEX) and then treated with different concentrations of HHXYT and FT. The results showed that testosterone levels were markedly elevated at the DEX concentration of 10, 20, 25 µM (Figure 9D). Whereas, HHXYT (50 µg/ml) and FT (3.2, 6.3, and 12.5 µM) could distinctly inhibit the increase of testosterone level caused by DEX (10 µM) in thecal cells, suggesting that HHXYT and FT have an inhibitory effect on testosterone production in thecal cells.
4 Discussion
MGH is the most prevalent breast disease and the incidence of breast cancer in MGH patients is remarkably higher than that in the average person (Basree et al., 2019). However, in the current clinical treatment of MGH, hormone therapy such as tamoxifen is not suitable for long-term use with apparent ill effects and complications (Wibowo et al., 2016); operative treatment is unacceptable for patients because of its complexity. Fortunately, HHXYT has been applied in the treatment of MGH for years and its efficacy has been validated clinically (Meng et al., 2014; Pei et al., 2016).
In this study, the integrated metabolomics and network pharmacology strategy were applied to systemically reveal the mechanism of HHXYT in the treatment of MGH. Firstly, the rat model of MGH was established to reconfirm the therapeutic effect of HHXYT. HHXYT effectively inhibited the abnormal proliferation of mammary lobules, glandular bubbles, and lumens in mammary glands, according to HE staining of mammary tissue. Meanwhile, the abnormal levels of estradiol (E2) and progesterone (P) in serum were reversed after treatment with HHXYT. In fact, the mammary gland is closely regulated by sex hormones (Li et al., 2017). Furthermore, it has been reported that the disturbance of sex hormones, particularly estrogen and progesterone, is the primary cause and defining feature of MGH: estrogen stimulation can accelerate breast proliferation, whereas progesterone can facilitate the growth and maturation of breast acinous cells and resist the proliferation impact of estrogen (Kleinberg et al., 2011; Chang et al., 2014). In other words, high levels of E2 or a lack of P cause estrogen to stimulate breast epithelium for an extended period without the protection of progesterone, resulting in significant mammary tissue hyperplasia (Liu et al., 2020). As a result, the above results confirmed that HHXYT could inhibit the symptoms of MGH visibly.
To figure out active compounds which exert therapeutic effects in HHXYT, we studied the plasma of rats after oral administration of HHXYT and identified 22 prototypes and their 32 metabolites. Following that, these compounds were used for HHXYT’s putative target prediction and functional enrichment analysis by network pharmacology. 373 putative targets, 93 pathways, 6 key herbs, and 9 key prototypes of HHXYT were sought out and the most relevant of these pathways was biological processes of the endocrine, especially the synthesis of hormones (Table 2). To further explore the therapeutic mechanism of HHXYT, the analysis of serum metabolomics was performed to intuitively characterize changes in endogenous metabolites. 58 endogenous different metabolites were identified, of which 31% were steroid lipids. In addition, the most significant metabolic modules were primarily involved in steroid hormone biosynthesis, cortisol synthesis and secretion, and ovarian steroidogenesis pathways. The compounds in the second module participated in primary bile acid biosynthesis and bile secretion. Notably, all these 4 pathways were also found to be enriched by HHXYT’s putative targets. It demonstrated that metabolomic confirmed the regulation of metabolic disturbances similar to network pharmacology’s prediction results, both of which were related to hormones and bile acid biosynthesis.
Actually, bile acid was closely related to the synthesis, absorption, and metabolism of cholesterol, which was the precursor of all hormone synthesis (Capponi, 2002). As shown in Figure 9E, in the pathway of sex hormone synthesis, cholesterol was applied as a raw material to synthesize androgens, and estrogens and androgens were precursors for the production of estrogen (Williams et al., 2001). The superiority of network pharmacology is that biological networks can be constructed to assist in elucidating drug treatment mechanisms from a holistic perspective (Zhao et al., 2019; Zhao et al., 2021). Ultimately, the herb-ingredient-target network was established as the substance basis for HHXYT to exert therapeutic effects on mammary gland hyperplasia. Among the 7 key targets of HHXYT, CYP11A1, HSD3B2, and CYP17A1 were directly involved in androgen synthesis and played vital roles, meanwhile, AR and ESR1 were receptors for the direct action of androgens and estrogens, respectively. In detail, estrogen was principally generated by the synergistic effects of thecal and granulosa cells of the ovary in females. Cholesterol was converted to pregnenolone in thecal cells by CYP11A1, which was turned into androgen under the catalysis of HSD3B2 and CYP17A1, especially CYP17A1. Then androgen was generated to estrogen in granulosa cells by CYP19A1 (Imamichi et al., 2017; Burris-Hiday and Scott, 2021). According to the above results, the synthesis of sex hormones, especially androgens, was significantly regulated by HHXYT.
Therefore, it was speculated that HHXYT may reduce the level of estrogen in vivo by inhibiting androgen synthesis, consequently protecting the mammary gland from excessive proliferation stimulated by estrogen. There were thecal cells and granulosa cells in the follicles of the ovary. Thecal cells possessed the luteinizing hormone (LH) receptor and synthesized testosterone (T), while granulosa cells owned follicle-stimulating hormone (FSH) receptors and produced E2. In addition, thecal cells would promote cell proliferation and testosterone production after being stimulated by LH (Li and Hearn, 2000). To evaluate the HHXYT’s effects on the synthesis of testosterone, we extracted thecal cells from the ovaries and experimented in vitro. Since dexamethasone (DEX) would promote the synthesis of T in thecal cells (Ong et al., 2019), DEX was added as the model group to observe the effects of HHXYT on testosterone synthesis. The results displayed that HHXYT (50 µg/ml) could distinctly inhibit the increase of T level caused by DEX (10 µM) in thecal cells. Meanwhile, formononetin (3.2, 6.3, and 12.5 µM), a prototype of HHXYT, also exhibited the inhibiting effect on testosterone levels rising. It demonstrated that they have an inhibitory effect on testosterone production in thecal cells. In summary, these findings contribute to elucidating the treatment mechanism of HHXYT on MGH and provide data and theoretical support for the in-depth study of its mechanism.
5 Conclusion
In this study, through integrating metabolomics and network pharmacology analysis, our results revealed the key substances for HHXYT to exert therapeutic effects on MGH, which included 7 targets, 6 herbs, and 17 ingredients. This study illustrated that HHXYT could regulate the hormone disorder prominently and the results suggested that HHXYT might reduce estrogen-stimulated mammary gland hyperplasia by inhibiting the production of its precursor androgen. The data and theoretical support were provided for the in-depth study of HHXYT’s mechanism. In addition, this study provides new ideas for the research of the mechanism of the multi-component TCM formula.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The animal study was reviewed and approved by the Animal Care and Use Committee of Second Military Medical University.
Author Contributions
WZ, JZ and ZG designed the study. ZG and RM performed the experiments. ZG, ZC and JZ analyzed the data. HZ and XL directed the design of animal and in vitro experiments. ZG and JZ wrote the manuscript. GW, JZ, JY, and WZ revised the manuscript. All the authors approved its final version and agreed to be accountable for all aspects of the work.
Funding
The work was supported by the National Key Research and Development Program of China (No. 2017YFC1700200), Professor of Chang Jiang Scholars Program, National Natural Science Foundation of China (Nos. 81520108030, 81572378, and 81973706), Shanghai Engineering Research Center for the Preparation of Bioactive Natural Products (No. 16DZ2280200), the Scientific Foundation of Shanghai China (No. 13401900103 and 13401900101).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.788019/full#supplementary-material
References
Adasme, M. F., Linnemann, K. L., Bolz, S. N., Kaiser, F., Salentin, S., Haupt, V. J., et al. (2021). PLIP 2021: Expanding the Scope of the Protein-Ligand Interaction Profiler to DNA and RNA. Nucleic Acids Res. 49, W530–W534. doi:10.1093/nar/gkab294
Alexander, N., Woetzel, N., and Meiler, J. (2011). “bcl: Cluster: A Method for Clustering Biological Molecules Coupled with Visualization in the Pymol Molecular Graphics System,” in IEEE Int. Conf. Comput. Adv. Bio Med. Sci. 2011, Orlando, FL, USA, 3-5 Feb. 2011, 13–18. doi:10.1109/iccabs.2011.5729867
Arendt, L. M., and Kuperwasser, C. (2015). Form and Function: How Estrogen and Progesterone Regulate the Mammary Epithelial Hierarchy. J. Mammary Gland Biol. Neoplasia 20, 9–25. doi:10.1007/s10911-015-9337-0
Basree, M. M., Shinde, N., Koivisto, C., Cuitino, M., Kladney, R., Zhang, J., et al. (2019). Abrupt Involution Induces Inflammation, Estrogenic Signaling, and Hyperplasia Linking Lack of Breastfeeding with Increased Risk of Breast Cancer. Breast Cancer Res. 21, 80. doi:10.1186/s13058-019-1163-7
Burley, S. K., Bhikadiya, C., Bi, C., Bittrich, S., Chen, L., Crichlow, G. V., et al. (2021). RCSB Protein Data Bank: Powerful New Tools for Exploring 3D Structures of Biological Macromolecules for Basic and Applied Research and Education in Fundamental Biology, Biomedicine, Biotechnology, Bioengineering and Energy Sciences. Nucleic Acids Res. 49, D437–D451. doi:10.1093/nar/gkaa1038
Burris-Hiday, S. D., and Scott, E. E. (2021). Steroidogenic Cytochrome P450 17A1 Structure and Function. Mol. Cel Endocrinol 528, 111261. doi:10.1016/j.mce.2021.111261
Capponi, A. M. (2002). Regulation of Cholesterol Supply for Mineralocorticoid Biosynthesis. Trends Endocrinol. Metab. 13, 118–121. doi:10.1016/s1043-2760(01)00538-0
Chang, X. J., Zhou, J., Zhang, S., Chen, J., Chen, C. M., Wang, Z. Z., et al. (2014). Effects of Guizhi Fuling Capsule on Sex Hormone Levels and Breast Issue Morphology of Mammary Gland Hyperplasia Model Rats. Chin. J. Chin. Mater. Med. 39, 4139–4142. doi:10.4268/cjcmm20142111
Chen, X., Zhou, H., Liu, Y. B., Wang, J. F., Li, H., Ung, C. Y., et al. (2006). Database of Traditional Chinese Medicine and its Application to Studies of Mechanism and to Prescription Validation. Br. J. Pharmacol. 149, 1092–1103. doi:10.1038/sj.bjp.0706945
Chen, Y., Yang, X., Qi, M., Yu, H., Du, Q., and Zhao, H. (2020). Observation on the Effect of Honghua Xiaoyao Tablets Assisted Western Medicine in the Treatment of Mastitis Guided by Real-Time Monitoring of Color Doppler Ultrasound. Mod. J. Integr. Tradit. West Med. 29, 2227–2229+2233. doi:10.3969/j.issn.1008-8849.2020.20.01
Chen, Z. M., Hu, J., Xu, Y. M., He, W., Meng, L., Huang, T., et al. (2021). Cryptotanshinone Inhibits Cytotoxin-Associated Gene A-Associated Development of Gastric Cancer and Mucosal Erosions. World J. Gastrointest. Oncol. 13, 693–705. doi:10.4251/wjgo.v13.i7.693
Chin, C. H., Chen, S. H., Wu, H. H., Ho, C. W., Ko, M. T., and Lin, C. Y. (2014). cytoHubba: Identifying Hub Objects and Sub-networks from Complex Interactome. BMC Syst. Biol. 8 (Suppl. 4), S11. doi:10.1186/1752-0509-8-s4-s11
Chong, J., Wishart, D. S., and Xia, J. (2019). Using MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Curr. Protoc. Bioinformatics 68, e86. doi:10.1002/cpbi.86
Coussens, L. M., and Pollard, J. W. (2011). Leukocytes in Mammary Development and Cancer. Cold Spring Harb Perspect. Biol. 3, a003285. doi:10.1101/cshperspect.a003285
Das, J., and Yu, H. (2012). HINT: High-Quality Protein Interactomes and Their Applications in Understanding Human Disease. BMC Syst. Biol. 6, 92. doi:10.1186/1752-0509-6-92
Fang, S., Dong, L., Liu, L., Guo, J., Zhao, L., Zhang, J., et al. (2021). HERB: a High-Throughput experiment- and Reference-Guided Database of Traditional Chinese Medicine. Nucleic Acids Res. 49, D1197–D1206. doi:10.1093/nar/gkaa1063
Hartmann, L. C., Sellers, T. A., Frost, M. H., Lingle, W. L., Degnim, A. C., Ghosh, K., et al. (2005). Benign Breast Disease and the Risk of Breast Cancer. N. Engl. J. Med. 353, 229–237. doi:10.1056/NEJMoa044383
Hopkins, A. L. (2007). Network Pharmacology. Nat. Biotechnol. 25, 1110–1111. doi:10.1038/nbt1007-1110
Imamichi, Y., Sekiguchi, T., Kitano, T., Kajitani, T., Okada, R., Inaoka, Y., et al. (2017). Diethylstilbestrol Administration Inhibits Theca Cell Androgen and Granulosa Cell Estrogen Production in Immature Rat Ovary. Sci. Rep. 7, 8374. doi:10.1038/s41598-017-08780-7
Jumper, J., Evans, R., Pritzel, A., Green, T., Figurnov, M., Ronneberger, O., et al. (2021). Highly Accurate Protein Structure Prediction with AlphaFold. Nature 596, 583–589. doi:10.1038/s41586-021-03819-2
Kanehisa, M., Furumichi, M., Tanabe, M., Sato, Y., and Morishima, K. (2017). KEGG: New Perspectives on Genomes, Pathways, Diseases and Drugs. Nucleic Acids Res. 45, D353–D361. doi:10.1093/nar/gkw1092
Kim, S., Chen, J., Cheng, T., Gindulyte, A., He, J., He, S., et al. (2021). PubChem in 2021: New Data Content and Improved Web Interfaces. Nucleic Acids Res. 49, D1388–D1395. doi:10.1093/nar/gkaa971
Kleinberg, D. L., Ameri, P., and Singh, B. (2011). Pasireotide, an IGF-I Action Inhibitor, Prevents Growth Hormone and Estradiol-Induced Mammary Hyperplasia. Pituitary 14, 44–52. doi:10.1007/s11102-010-0257-0
Li, S., and Zhang, B. (2013). Traditional Chinese Medicine Network Pharmacology: Theory, Methodology and Application. Chin. J. Nat. Med. 11, 110–120. doi:10.1016/s1875-5364(13)60037-0
Li, S. K., and Hearn, M. T. (2000). Isolation of Thecal Cells: an Assessment of Purity and Steroidogenic Potential. J. Biochem. Biophys. Methods 45, 169–181. doi:10.1016/s0165-022x(00)00107-x
Li, T., Zhang, W., Hu, E., Sun, Z., Li, P., Yu, Z., et al. (2021). Integrated Metabolomics and Network Pharmacology to Reveal the Mechanisms of Hydroxysafflor Yellow A against Acute Traumatic Brain Injury. Comput. Struct. Biotechnol. J. 19, 1002–1013. doi:10.1016/j.csbj.2021.01.033
Li, H. T., Liu, H. H., Yang, Y. X., Wang, T., Zhou, X. L., Yu, Y., et al. (2018). Therapeutic Effects of a Traditional Chinese Medicine Formula Plus Tamoxifen vs. Tamoxifen for the Treatment of Mammary Gland Hyperplasia: A Meta-Analysis of Randomized Trials. Front. Pharmacol. 9, 45. doi:10.3389/fphar.2018.00045
Li, X., Xin, P., Wang, C., Wang, Z., Wang, Q., and Kuang, H. (2017). Mechanisms of Traditional Chinese Medicine in the Treatment of Mammary Gland Hyperplasia. Am. J. Chin. Med. 45, 443–458. doi:10.1142/s0192415x17500276
Liebisch, G., Fahy, E., Aoki, J., Dennis, E. A., Durand, T., Ejsing, C. S., et al. (2020). Update on LIPID MAPS Classification, Nomenclature, and Shorthand Notation for MS-derived Lipid Structures. J. Lipid Res. 61, 1539–1555. doi:10.1194/jlr.S120001025
Liu, J., Zhang, J. F., Gong, G. H., Zhang, B., and Wei, C. X. (2020). Therapeutic Effect of Mongolian Medicine RuXian-I on Hyperplasia of Mammary Gland Induced by Estrogen/Progesterone through CRYAB-Promoted Apoptosis. Evid. Based Complement. Alternat Med. 2020, 5707106. doi:10.1155/2020/5707106
Luan, H., Ji, F., Chen, Y., and Cai, Z. (2018). statTarget: A Streamlined Tool for Signal Drift Correction and Interpretations of Quantitative Mass Spectrometry-Based Omics Data. Anal. Chim. Acta 1036, 66–72. doi:10.1016/j.aca.2018.08.002
Meng, Q., Liu, S., and Wu, X. (2014). Clinical Observation on 215 Cases of Mammary Gland Hyperplasia Treated with Honghua Xiaoyao Tablets. Hebei J. TCM 36, 1536–1537.
Mi, R., Li, X., Zhang, Z., Cheng, T., Tian, S., Xu, X., et al. (2020). Chemical Profiling of Honghua Xiaoyao Tablet and Simultaneous Determination of its Quality Markers by Liquid Chromatography-Tandem Mass Spectrometry Combined with Chemometrics Methods. J. Sep. Sci. 43, 4263–4280. doi:10.1002/jssc.202000689
Newgard, C. B. (2017). Metabolomics and Metabolic Diseases: Where Do We Stand?. Cell Metab 25, 43–56. doi:10.1016/j.cmet.2016.09.018
Nicholson, J. K., and Wilson, I. D. (2003). Opinion: Understanding 'global' Systems Biology: Metabonomics and the Continuum of Metabolism. Nat. Rev. Drug Discov. 2, 668–676. doi:10.1038/nrd1157
Ong, M., Cheng, J., Jin, X., Lao, W., Johnson, M., Tan, Y., et al. (2019). Paeoniflorin Extract Reverses Dexamethasone-Induced Testosterone Over-secretion through Downregulation of Cytochrome P450 17A1 Expression in Primary Murine Theca Cells. J. Ethnopharmacol 229, 97–103. doi:10.1016/j.jep.2018.09.006
Pei, S., Zhang, Q., and Yang, J. (2016). Clinical Observation of Honghua Xiaoyao Tables in the Treatment of Delayed Menstrual Cycle of Liver-Depression and Qi-Blood Stagnation. China Pharm. 27, 245–247. doi:10.6039/j.issn.1001-0408.2016.02.3210.1353/jod.2016.0057
Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: a Database of Systems Pharmacology for Drug Discovery from Herbal Medicines. J. Cheminform 6, 13. doi:10.1186/1758-2946-6-13
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303
Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analyses. Curr. Protoc. Bioinformatics 54, 1.30.1–1.30.33. doi:10.1002/cpbi.5
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING V11: Protein-Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-wide Experimental Datasets. Nucleic Acids Res. 47, D607–D613. doi:10.1093/nar/gky1131
Tanchuk, V. Y., Tanin, V. O., Vovk, A. I., and Poda, G. (2016). A New, Improved Hybrid Scoring Function for Molecular Docking and Scoring Based on AutoDock and AutoDock Vina. Chem. Biol. Drug Des. 87, 618–625. doi:10.1111/cbdd.12697
Trott, O., and Olson, A. J. (2010). AutoDock Vina: Improving the Speed and Accuracy of Docking with a New Scoring Function, Efficient Optimization, and Multithreading. J. Comput. Chem. 31, 455–461. doi:10.1002/jcc.21334
Wang, J. (2012). 35 Cases of Perimenopausal Syndrome Treated with Honghua Xiaoyao Tablets. Chin. J. Exp. Tradit Med. Form 18, 312–313. doi:10.13422/j.cnki.syfjx.2012.18.01510.1007/s11655-012-1064-4
Wang, L., Zhao, D., Di, L., Cheng, D., Zhou, X., Yang, X., et al. (2011). The Anti-hyperplasia of Mammary Gland Effect of Thladiantha Dubia Root Ethanol Extract in Rats Reduced by Estrogen and Progestogen. J. Ethnopharmacol 134, 136–140. doi:10.1016/j.jep.2010.11.071
Wibowo, E., Pollock, P. A., Hollis, N., and Wassersug, R. J. (2016). Tamoxifen in Men: a Review of Adverse Events. Andrology 4, 776–788. doi:10.1111/andr.12197
Williams, K., Mckinnell, C., Saunders, P. T., Walker, M., Fisher, J. S., Turner, K. J., et al. (2001). Neonatal Exposure to Potent and Environmental Oestrogens and Abnormalities of the Male Reproductive System in the Rat: Evidence for Importance of the Androgen-Oestrogen Balance and Assessment of the Relevance to Man. Hum. Reprod. Update 7, 236–247. doi:10.1093/humupd/7.3.236
Wishart, D. S., Feunang, Y. D., Marcu, A., Guo, A. C., Liang, K., Vázquez-Fresno, R., et al. (2018). HMDB 4.0: the Human Metabolome Database for 2018. Nucleic Acids Res. 46, D608–D617. doi:10.1093/nar/gkx1089
Wu, G., Zhang, W., and Li, H. (2019). Application of Metabolomics for Unveiling the Therapeutic Role of Traditional Chinese Medicine in Metabolic Diseases. J. Ethnopharmacol 242, 112057. doi:10.1016/j.jep.2019.112057
Wu, Y., Zhang, F., Yang, K., Fang, S., Bu, D., Li, H., et al. (2019). SymMap: an Integrative Database of Traditional Chinese Medicine Enhanced by Symptom Mapping. Nucleic Acids Res. 47, D1110–D1117. doi:10.1093/nar/gky1021
Xu, H. Y., Zhang, Y. Q., Liu, Z. M., Chen, T., Lv, C. Y., Tang, S. H., et al. (2019). ETCM: an Encyclopaedia of Traditional Chinese Medicine. Nucleic Acids Res. 47, D976–D982. doi:10.1093/nar/gky987
Xu, W., Zhang, Y., Yu, Y., Li, B., Liu, J., Wang, P., et al. (2018). Dose-dependent Target Diversion of Danhong Injection on the Glu-GLT-1/Gly-GlyRα Dynamic Balance Module of Cerebral Ischemia. Pharmacol. Res. 135, 80–88. doi:10.1016/j.phrs.2018.07.020
Xue, R., Fang, Z., Zhang, M., Yi, Z., Wen, C., and Shi, T. (2013). TCMID: Traditional Chinese Medicine Integrative Database for Herb Molecular Mechanism Analysis. Nucleic Acids Res. 41, D1089–D1095. doi:10.1093/nar/gks1100
Yao, Z. J., Dong, J., Che, Y. J., Zhu, M. F., Wen, M., Wang, N. N., et al. (2016). TargetNet: a Web Service for Predicting Potential Drug-Target Interaction Profiling via Multi-Target SAR Models. J. Comput. Aided Mol. Des. 30, 413–424. doi:10.1007/s10822-016-9915-2
Ye, J., Li, R., Zeng, H., Liu, R., and Zhang, W. (2019). Discovery and Research Advances in Quality Markers of Chinese Material Medica Based on Holistic Characteristics. Chin. Tradit Herbal Drugs 50, 4529–4537. doi:10.7501/j.issn.0253-2670.2019.19.002
You, Z., Sun, J., Xie, F., Chen, Z., Zhang, S., Chen, H., et al. (2017). Modulatory Effect of Fermented Papaya Extracts on Mammary Gland Hyperplasia Induced by Estrogen and Progestin in Female Rats. Oxid. Med. Cell. Longev. 2017, 8235069. doi:10.1155/2017/8235069
Zhao, J., Lv, C., Wu, Q., Zeng, H., Guo, X., Yang, J., et al. (2019). Computational Systems Pharmacology Reveals an Antiplatelet and Neuroprotective Mechanism of Deng-Zhan-Xi-Xin Injection in the Treatment of Ischemic Stroke. Pharmacol. Res. 147, 104365. doi:10.1016/j.phrs.2019.104365
Zhao, J., Tian, S., Lu, D., Yang, J., Zeng, H., Zhang, F., et al. (2021). Systems Pharmacological Study Illustrates the Immune Regulation, Anti-infection, Anti-inflammation, and Multi-Organ protection Mechanism of Qing-Fei-Pai-Du Decoction in the Treatment of COVID-19. Phytomedicine 85, 153315. doi:10.1016/j.phymed.2020.153315
Keywords: Hong-Hua-Xiao-Yao tablet, mammary gland hyperplasia, traditional medicine, metabolomics, network pharmacology, pharmacological effects
Citation: Gao Z, Mi R, Cheng Z, Li X, Zeng H, Wu G, Zhao J, Zhang W and Ye J (2022) Integrated Metabolomics and Network Pharmacology Revealed Hong-Hua-Xiao-Yao Tablet’s Effect of Mediating Hormone Synthesis in the Treatment of Mammary Gland Hyperplasia. Front. Pharmacol. 13:788019. doi: 10.3389/fphar.2022.788019
Received: 01 October 2021; Accepted: 07 January 2022;
Published: 01 February 2022.
Edited by:
Yuanjia Hu, University of Macau, ChinaReviewed by:
Jianxin Chen, Beijing University of Chinese Medicine, ChinaMarta C. Romano, Instituto Politécnico Nacional de México (CINVESTAV), Mexico
Copyright © 2022 Gao, Mi, Cheng, Li, Zeng, Wu, Zhao, Zhang and Ye. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jing Zhao, Zhaojanne@gmail.com; Weidong Zhang, wdzhangy@hotmail.com; Ji Ye, catheline620@163.com