Integrated Metabolomics and Network Pharmacology Approach to Explain Possible Action Mechanisms of Xin-Sheng-Hua Granule for Treating Anemia

As a well-known traditional Chinese medicine (TCM) prescription, Xin-Sheng-Hua Granule (XSHG) has been applied in China for more than 30 years to treat postpartum diseases, especially anemia. However, underlying therapeutic mechanisms of XSHG for anemia were still unclear. In this study, plasma metabolomics profiling with UHPLC-QTOF/MS and multivariate data method was firstly analyzed to discover the potential regulation mechanisms of XSHG on anemia rats induced by bleeding from the orbit. Afterward, the compound-target-pathway network of XSHG was constructed by the use of network pharmacology, thus anemia-relevant signaling pathways were dissected. Finally, the crucial targets in the shared pathways of metabolomics and network pharmacology were experimentally validated by ELISA and Western Blot analysis. The results showed that XSHG could exert excellent effects on anemia probably through regulating coenzyme A biosynthesis, sphingolipids metabolism and HIF-1α pathways, which was reflected by the increased levels of EPOR, F2, COASY, as well as the reduced protein expression of HIF-1α, SPHK1, and S1PR1. Our work successfully explained the polypharmcological mechanisms underlying the efficiency of XSHG on treating anemia, and meanwhile, it probed into the potential treatment strategies for anemia from TCM prescription.

Introduction

Anemia, one of the most common diseases in obstetrics, could cause severe maternal and fetal complications, including preterm birth and placental mesenchymal dysplasia (Mousa et al., 2014). On a global scale, more than 500,000 maternal deaths occurred during pregnancy each year, and 20% of them are resulted from anemia and postpartum hemorrhage (Balarajan et al., 2011; Milman, 2011). Researchers have demonstrated that postpartum anemia was commonly caused by blood clots, coagulation abnormalities, retained placenta or genital tract trauma (Froessler et al., 2014; Markova et al., 2015). However, the available therapeutic drugs and their therapeutic effects were limited, and some of them even caused obvious side effects (Milman, 2012; Tunçalp et al., 2013). For these reasons, the discoveries of potential strategies for treating anemia could have a positive effect on maternal as well as fetal outcomes.

Traditional Chinese medicine (TCM) prescriptions are usually comprised of some specific herbs in order to achieve the synergistic effects in a holistic way (Yang et al., 2017). Comparing with modern medicine, they had obvious advantages of fewer side effects for treating multi-factorial diseases (Wang et al., 2016b). In the past few decades, Xin-Sheng-Hua granule (XSHG) has been widely applied for the treatment of anemia-related diseases in China, especially for the anemia caused by postpartum hemorrhage, whose prescription is the boiled water extraction followed by seven common herbs, Angelica sinensis (Oliv.) Diels (ASD; Danggui), Leonurus artemisia (Laur.) S. Y. Hu F (LA; Yimucao), Ligusticum chuanxiong Hort. (LC; Chuanxiong), Prunus persica (L.) Batsch (PP; Taoren), roasted rhizome of Zingiber officinale Rosc. (RZR; Jiangtan), Radix Glycyrrhizae preparata (RGP; Zhigancao) and Carthamus tinctorius L. (CT; Honghua), at a proportion of 80:90:30:8:5:5:5 (Pang et al., 2016b). According to the National Bureau of Statistics of China¹, its annual sales are expected to exceed 100 million dollars by 2020, and it ranked the forefront among all of the obstetrical medications. A growing body of evidence indicated that XSHG could exhibit excellent therapeutic effects on anemia via promoting hematopoietic stem cell proliferation and differentiation (Pang et al., 2017; Zhou et al., 2017). However, up to date, the action mechanism for treating anemia of XSHG remained poorly understood.

Metabolomics aims at constructing the metabolic profile of endogenous metabolites with low molecular weight in biological systems through modern analytical techniques (Naz et al., 2014). It could be applied to monitor and diagnose disease progression by evaluating the content variations of metabolic biomarkers between control and model groups (Courant et al., 2014). Due to these advantages, increasing attention has been focused on metabolomics to reveal the mechanisms of TCM prescriptions, such as Huangqin Decoction (Cui et al., 2017) and Liu Wei Di Huang Wan (Wang et al., 2010). Network pharmacology is now regarded as a holistic and efficient technique to study the role of TCM prescriptions (Yue et al., 2017b). It can help to understand the active compounds and therapeutic targets of TCMs by combining pharmacology and pharmacodynamics analysis (Wang et al., 2017). Thus this approach is quite important for interpreting the underlying mechanisms of TCMs (Chen et al., 2016; Guo et al., 2016). Recently, some researchers have successfully used the integrated metabolomics and network pharmacology strategy to explore the interactions between organisms and drugs (Kell and Goodacre, 2014; Li S. et al., 2014), bringing great inspiration to the mechanism research of XSHG.

In current study, in order to explain possible action mechanisms of XSHG for treating anemia, we performed the integrated strategy of metabolomics and network pharmacology (Figure 1). Firstly, the anemia model was induced by bleeding from the orbit. Secondly, the metabolic pathways were performed by studying plasma metabolomics through ultra-high-performance liquid chromatography combined with quadrupole time-of-flight mass spectrometry (UHPLC-QTOF/MS) and multivariate data method. Thirdly, the anemia-relevant pathways were profiled and potential mechanisms of XSHG were revealed through network pharmacology. Finally, the crucial targets of the shared pathways of metabolomics and network pharmacology were experimentally validated to account for the therapeutic effects of XSHG. This approach could offer a significant guidance for elucidating mechanism of XSHG on anemia, and then promote its clinical application.

FIGURE 1

FIGURE 1. The integrated strategy of metabolomics and network pharmacology. C-T-P network: compound-target-pathway network.

Materials and Methods

Chemicals and Reagents

The prescription XSHG was generously provided by Jiangsu Revolence Pharmaceutical Co., Ltd. (Huai’an, China). Acetonitrile, methanol and formic acid of LC-MS grade were purchased from Merck KGaA (Darmstadt, Germany). Ethylene diamine tetraacetic acid-2Na (EDTA-2Na) was offered by Shanghai Danning Chemical Co., Ltd. (Shanghai China). The deionized water purified by Milli-Q water purification system (Millipore, Billerica, MA, United States) was applied to prepare and extract plasma samples. Other reagents and chemicals were all analytical grade.

Animals and Drug Administration

Sprague-Dawley (SD) female rats (220 - 250 g) were purchased from Shanghai Slac Laboratory Animal Co., Ltd., (Shanghai, China). Before the experiments, the animals were housed in an environmentally controlled breeding room for 12 h light-dark cycle, at a relative humidity of 55 ± 5°C and a temperature for 20 ± 2°C. All experimental animals were free access to food and tab water during this study. The protocol was strictly in accordance with the Guide for the Care and Use of Laboratory Animals (US National Research Council, 1996), and was approved by the Animal Experimental Ethics Committee of Nanjing University of Chinese Medicine.

After 10 days for acclimatization, 30 SD rats were randomly divided into 3 groups with 10 rats in each: the control, model and XSHG groups. The method of inducing anemia model was according to our previous research (Shi et al., 2014a). Except the control group, rats in other groups were performed by bleeding from the orbit at the account of 5.0 mL/kg each day, and the experimental period was 12 successive days.

The XSHG extract was dispersed and dissolved homogeneously by using pure water, and its dosage of 9.86 g/kg in XSHG group was equivalent to two times of the adult daily dose of its prescription (36 g), which was inferred from the following formula: human dose of crude herbs in clinic × 0.018/200 × 1000 × the multiple of clinical equivalency dose (Shi et al., 2014b). Model and control groups were intragastrically administered with the same volume of saline solution, and all animals were given orally one time each day for 12 days. After 1 h of the last administration, 1.0 mL blood samples were collected into 1.5 mL centrifuge tubes via the orbit. 200 μL whole blood was applied to evaluate peripheral blood routine, including hemoglobin (HGB/g.L^-1), red blood cell count (RBC/10¹².L^-1), white blood cell count (WBC/10⁹.L^-1), platelet count (PLT/10⁹.L^-1), lymphocyte ratio (LY/%), monocyte ratio (MO/%), neutrophils ratio (GR/%) and haematocrit (HCT/%). The rest of whole blood samples were immediately centrifuged at 3000 rpm for 10 min, then the supernatant was transferred into another 1.5 mL centrifuge tube. Tissue samples from right lobe of liver were grinded with saline solution (1:3, v/v). All of the plasma and liver samples were kept at -80°C until analysis.

Metabolomics Study

Sample Preparation

Plasma samples were thawed at room temperature before preparation. To precipitate protein, 200 μL plasma samples were extracted with 600 μL methanol, and the mixture was then vortexed for 1 min and centrifuged at 3000 rpm for 10 min. Afterward, 400 μL supernatants were transferred into new centrifuge tubes, and were evaporated to dryness at 37°C by using Labconco CentriVap concentrator (Kansas City, MO, United States). The residues were re-dissolved with 100 μL of initial mobile phase, vortexed for 3 min and centrifuged at 13000 rpm for 10 min at 4°C. Finally, an aliquot of 5 μL supernatant was injected for UHPLC-MS analysis.

Furthermore, 15 plasma samples were randomly selected from these three groups, which were mixed together as quality control (QC) samples. As QC samples contained the most data of each group, it was applied to validate the stability of UHPLC-MS system. Before analysis, QC samples were analyzed for six times to equilibrate the UHPLC system, then every ten samples were used to monitor the stability of this method (Godzien et al., 2011; Prosser et al., 2014). All samples above were maintained at 4°C during analysis.

Liquid Chromatography

Chromatographic separation was performed on an UHPLC^TM system (Waters, Corp., Milford, MA, United States). Considering the complicated analytes, a Thermo Scientific Hypersil GOLD column (100 mm × 3 mm, 1.9 μm) was firstly applied for all the analytes. And the column temperature was maintained at 35°C. Chromatographic analysis was performed with gradient elution using water/0.1% formic acid (solvent A) and acetonitrile (solvent B), which programmed as follows: 5–40% B from 0 to 3 min; 40–52% B from 3 to 3.5 min; 52% B from 3.5 to 6 min; 52–95% B from 6 to 17 min; 95% B from 17 to 19 min; 5% B from 19 to 19.5 min, held at 5% B for 3 min to equilibrate the column. The mobile phase was directly delivered into electrospray ion (ESI) source at 0.4 mL/min.

Mass Spectrometry

In the experiment, a Waters Synapt^TM QTOF/MS spectrometer (Waters Corp., Milford, MA, United States) was combined with UHPLC system via the ESI source in both positive and negative ionization modes. The detailed parameters were as follows: extraction cone voltage of 2.0 V, sample cone voltage of 30.0 V, capillary voltage of 3.0 kV, desolvation temperature of 400°C, and source temperature of 120°C. Nitrogen was used as the cone and desolvation gas at the flow rate of 50 and 800 L/h accordingly. The metabolomics in centroid mode were set from 100 to 1000 Da. To ensure accuracy, leucine encephalin was selected as lock-mass solution for the acquisition of accurate mass.

Pattern Recognition Analysis and Data Processing

The acquired raw data were introduced to MakerLynx within MassLynx software for peak alignment and detection, then m/z data and retention time of each analyte could be obtained. The primary parameters were: mass range 100–1000 amu, retention time range 1–20 min, mass tolerance 0.1 Da, and noise elimination level 5. After normalization, the resultant data matrices were imported to EZinfo 2.0 for principal component analysis (PCA), partial least-squares-discriminant analysis (PLS-DA), and orthogonal projection to latent structures (OPLS) analysis (Ciborowski et al., 2012). Prior to multivariate statistical analysis, raw data acquired from UHPLC–MS were scaled to Pareto variance. In PLS-DA score plots, the variables with VIP > 1 were usually regarded as candidate biomarkers and were subjected to further identification of the chemical formula. Pathway analysis was conducted with MetaboAnalyst, a web-based tool for visualization of metabolomics results (Park et al., 2015; Takayama et al., 2015).

Network Pharmacology Study

Compound Data Preparation and ADME Screening

All of the compounds data of seven herbs in XSHG were obtained from traditional Chinese medicine systems pharmacology database and analysis platform,² (TCMSP^TM). Afterward, three important in silico ADME indexes including OB (oral bioavailability), Caco-2 (cell permeability), and DL (drug-likeness) were employed to screen the candidate active ingredients (Cao et al., 2015). The threshold values of these three indexes were set as follows: OB ≥ 30%, Coca-2 ≥-0.4, and DL ≥ 0.18, respectively (Liu et al., 2016). These ingredients, which met all of three criteria above, were selected as candidate molecules for further analysis. Moreover, some ingredients were also supplemented manually by a wide-scale data-mining way.

Targets Fishing and Network Construction

A comprehensive approach of information integration, chemometric method and text-mining was introduced to discover the targets of XSHG. Firstly, the most likely biological targets were retrieved from PharmMapper sever (Yue et al., 2017a), STITCH³, and SEA⁴ (Keiser et al., 2007; Liu et al., 2010; Kuhn et al., 2012). Afterward, the candidate bioactive ingredients were imported to Google Scholar, DrugBank⁵, TTD⁶, and Herbal Ingredients Targets database (HIT) (Yue et al., 2017a) to discover the relevant targets supported by literature. Besides, to better understand the mechanism of XHSG for the treatment of anemia, all obtained targets were also sent to TTD, PharmGKB⁷, and Comparative Toxicogenomics Database (CTD⁸) (Ye et al., 2011; Zhu et al., 2012). And targets which were implicated in the clinical manifestations of anemia were retained, while others were eliminated. Finally, the signaling pathways related to anemia were extracted from KEGG database, and compound-target-pathway (C-T-P) network was generated by using Cytoscape 3.0 software.

Experimental Validation

Analysis of Key Targets

Based on the results of metabolomics and network pharmacology, we focused on their shared pathways to discover the potential mechanisms. The metabolites of these pathways could be found, then metabolite-protein interactions from HMDB, protein–protein interactions from KEGG, and the relevant interaction genes in STRING could be further induced. As a result, the network containing relationships among metabolites, enzymes, pathway-related proteins and relevant interaction genes was established. By analyzing this network, the key targets in this interaction network were selected for experimental validation. Furthermore, other significant targets closely related to anemia were also added based on the previous reports. Collectively, these targets obtained above were validated by ELISA and Western blot analysis.

ELISA Analysis

Plasma levels of TNF-α, IL6, F2, and EPOR in different groups were measured by commercial ELISA kits (R&D Systems, United States) according to the protocols provided from the manufacture. Absorbance was detected at 450 nm by the EnSpire^TM microplate reader (PerkingElmer, United States).

Western Blot

To extract the protein of liver, the liver samples were homogenized with whole lysis buffer (1 mmol/L phenylmethylsulfonyl fluoride, 10 mmol/L Tris-HCl, 30 mmol/L sodium pyrophosphate, 50 mmol/L sodium fluoride, 250 mmol/L sodium chloride, 0.5% Triton X-100, 10% glycerol, 1 × proteinase inhibitor mixture, 2 mmol/L iodoacetic acid, and 5 mmol/L ZnCl₂). Western Blot assays of ACSS1, COASY, CBS, AHCY, SPHK1, S1PR1, and HIF-1α (Proteintech, United States) were carried out by standard protocols based on the manufacturer’s instructions.

Statistics

Statistical analysis was performed by using SPSS 19.0 software, all data were expressed as mean ± standard deviation. The statistical results were conducted with two-tailed unpaired Student’s t-test, and value of P ≤ 0.05 was regarded significant difference.

Results

Peripheral Blood Routine Analysis

The peripheral blood routine analysis was firstly evaluated to test whether the anemia model was constructed successfully. As shown in Table 1, compared with control group, HGB, RBC, and LY% in model group were obviously reduced (P < 0.01), while PLT, WBC, HCT, MO%, and GR% in model group were significantly increased (P < 0.01). The results indicated that anemia model was successfully induced. The levels of PLT, WBC, HCT, MO, and GR in XSHG were significantly reduced compared with model group (P < 0.05), whereas the decreased HGB, RBC, and LY% in model group were obviously elevated after oral administration of XSHG (P < 0.01). The MS data were successfully submitted to the online MassIVE datasets⁹, and the accession number was MSV000081992.

TABLE 1

TABLE 1. Effects of XSHG after 12 days administration on peripheral blood parameters of anemia rats (-X ± SD, n = 10).

Results of Plasma Metabolomics

Metabolomics Profiling of Plasma Samples

Under the optimized chromatography conditions, the typical based peak intensity (BPI) chromatograms of plasma samples were analyzed in both positive and negative modes (Figure 2). And the endogenous markers obtained excellent separation within 20 min. Among different groups, the subtle changes from complex MS data could be discovered through multivariate data analysis methods, including PCA, PLS-DA, and OPLS-DA.

FIGURE 2

FIGURE 2. Representative base peak intensity chromatograms of plasma samples of different groups derived from UHPLC-QTOF/MS. (A) Control group, (B) model group, and (C) XSHG group.

Multivariate Data Analysis

Specifically, to observe the subtle changes, raw data from UHPLC-MS were imported into EZinfo 2.0 software. Then an obvious separation in the PCA score plots between the model and control groups was clearly observed, indicating that biochemical perturbation occurred in the anemia group (Figure 3A). Then the S-plot of PLS-DA (Figure 3B) and VIP-value plot (Figure 3C) were drawn to unveil the metabolic biomarkers of anemia. Generally, the variables of VIP plots exhibiting values larger than 1 could be considered as the potential biomarkers, which were responsible for the differences between model and control groups (Wang et al., 2016a). The potential endogenous markers could also be discovered from the loading plot of the PLS-DA (Figure 3D). R²Y of the PLS-DA model in positive and negative modes was 0.996 and 0.978, respectively; Q² was 0.938 and 0.855 accordingly, which showed that this PLS-DA model was excellent for prediction and fitness.

FIGURE 3

FIGURE 3. The plasma PCA score plot (A); S-plot of OPLS-DA (B); and VIP-plot of OPLS-DA (C) between control and model groups and loading plot of PLS-DA (D) of control, model and XSHG groups in negative mode (n = 10).

Moreover, the relative distances from treatment groups (model and XSHG groups) to control group in the PLS-DA score plot were applied for quantitation (Supplementary Figure S1). And the mean metabolic patterns of normal rats were defined as the jumping-off point of other groups (Duan et al., 2016). In comparison with model group, the relative distance of XSHG group decreased obviously in both positive and negative modes, demonstrating that XSHG could adjust the anemia to the normal (Table 2).

TABLE 2

TABLE 2. Relative distance between drug groups (model and XSHG) and control group from the PLS-DA score plot of the plasma samples (± SD, n = 10).

Identification of Endogenous Metabolites

The potential endogenous biomarkers could be identified based on MS/MS fragments, retention behavior and online databases query. Firstly, candidate biomarkers were extracted from S-plots and loading plots performed, following analysis of OPLS-DA from MarkerLynx (Wettersten and Weiss, 2013). The accurate molecular mass and m/z values of metabolites were then determined. Subsequently, we could predict these metabolites initially through querying some on-line databases (HMDB, KEGG, and METLIN). And the potential markers were analyzed within an appropriate degree of measurement error (<10 ppm). Moreover, the MS/MS spectrum of different fragment ions were obtained in a targeted MS/MS mode. At last, these endogenous metabolites were identified by comparing with standard references and mass assignments in online databases (Xu et al., 2015).

Based on the protocols described above, a total of 18 endogenous biomarkers in plasma samples were tentatively identified (Table 3). LysoPC (20:4) and deoxycholic acid were identified by chemical standards, whereas other endogenous markers were identified tentatively by comparing accuracy mass information and corresponding assignments of product ions. To illustrate the identification process, one endogenous biomarker (Rt = 9.78, m/z = 543.6729) with VIP value of 8.03 was detailed as an example. Firstly, the exact mass of this biomarker ([M + H]⁺ at 544.0729) was obtained based on the mass spectrum in total ion chromatogram. Afterward, specific MS/MS fragments of this biomarker were determined by QTOF system. From the positive spectrum, main fragment ions of the candidate metabolite were observed at m/z 526.33, 484.29, 258.10, 184.07, 125.00, and 104.12, which could be the [M + H]⁺ of lost H₂O, C₃H₁₀N, C₂₃H₃₇O₂, C₂₆H₄₄O₃, and C₂₄H₄₂O₄P, respectively. Molecular formula C₂₈H₅₀NO₇P was located as the potential marker by the use of measurement error of 5 ppm. Finally, this biomarker was identified as LysoPC (20:4) through querying HMDB and METLIN databases.

TABLE 3

TABLE 3. The disturbed plasma endogenous metabolites between the control and anaemia rats as well as their identification results.

Among the 18 endogenous metabolites, the levels of eicosapentaenoic acid, docosahexaenoic acid and 16-dehydroprogesterone in model group were up-regulated significantly (P < 0.05), whereas other biomarkers in model group were down-regulated obviously (P < 0.01). Furthermore, the relative intensities of 18 endogenous markers were determined for semi-quantitation, and it was discovered that the contents of these metabolites in XSHG group tended to the normal group (Supplementary Figure S2).

Metabolic Pathway Analysis

To explore the potential metabolic pathways of anemia, 18 endogenous biomarkers were then imported into MetaboAnalyst. Ten corresponding pathways were constructed, which could be highly related to the anemia (Supplementary Table S1). Among these metabolic pathways, pantothenate and CoA biosynthesis with impact-value 0.33, cysteine and methionine metabolism with impact-value 0.22, and sphingolipid metabolism with impact-value 0.17 were filtered out as significant pathways (Figure 4), for the reason that the threshold of the impact-value over 0.10 was usually considered as the most important metabolic pathways (Hoffmann et al., 2014).

FIGURE 4

FIGURE 4. Summary of pathway analysis of plasma samples collected from anaemia rats. (a) pantothenate and CoA biosynthesis, (b) cysteine and methionine metabolism, (c) sphingolipid metabolism, (d) steroid hormone biosynthesis, (e) glycerophospholipid metabolism, (f) biosynthesis of unsaturated fatty acids, (g) nitrogen metabolism, (h) retinol metabolism, (i) glycine, serine and threonine metabolism;, (j) purine metabolism.

Of the 18 potential metabolites identified from these pathways, several played important roles in the development of anemia. 4-Phosphopantothenoylcysteine was a key intermediate metabolite of coenzyme A (CoA), which could provide 90% of the body energy. Pantothenoylcysteine could be combined with triphosphate and adenosine to form CoA, could promote hematopoietic function through improving energy metabolism, activating the immune system, and increasing the synthesis of erythrocyte. In anemia patients, the activities of ATPase decreased obviously, which might result from the low levels of 4-phosphopantothenoylcysteine (Pietrocola et al., 2015). S-Adenosylhomocysteinase and L-Cystathionine were mainly involved in cysteine and methionine metabolism, and they could play key roles in cardiovascular disease and biosynthesis of folic acid (Pacana et al., 2015). The decreased contents suggested that cysteine and methionine metabolism was inhibited, besides, several decreased biomarkers like sphinganine, and sphingosine 1-phosphate could be found and applied to explain sphingolipid metabolism. Researches have demonstrated that sphingolipids played crucial roles in metabolic diseases by regulating cell proliferation, differentiation and death (Harvald et al., 2015). The contents of sphinganine and sphingosine 1-phosphate reduced in anemia, which could be caused by the injured membrane of erythrocyte and disordered peripheral blood parameters (Li et al., 2015).

Correlation Analysis Between Biochemical Indicators and Biomarkers

To test whether biochemical indicators have relations with biomarkers, the Pearson correlation matrix was applied for the correlation analysis. As was shown in Figure 5, there was a certain correlation between peripheral blood routine and endogenous metabolites in anemia rats. The markers 4-phosphopantothenoylcysteine, dIMP, 8-isoprostaglandin F2a, ubiquinone-2, sphinganine, sphingosine 1-phosphate, 4-oxo-retinoic acid and leucyl-phenylalanine were strong positive associated with WBC, PLT and MO (r > 0.5); 4-phosphopantothenoylcysteine, dIMP, 8-isoprostaglandin F2a, sphinganine, sphingosine 1-phosphate and 4-oxo-retinoic acid had strong negative association with RBC, HGB, HCT, LY, and NE (r < -0.5); additionally, 15-deoxy-d-12,14-PGJ2 had strong positive association with WBC and MO (r = 0.776, 0.996), as well as strong negative association with RBC, HGB, HCT, and LY (r = -0.737, -0.725, -0.509, -0.683). Deoxycholic acid was significantly positive associated with WBC and MO (r = 0.702, 0.780), and negative associated with RBC, HGB, and HCT (r = -0.803, -0.674, -0.643). The level of PLT was positive associated with LysoPC (20:4) (r = 0.503). Eicosapentaenoic acid showed strong positive association with WBC and LY (r = 0.560, 0.612). An obvious positive correlation was discovered between docosahexaenoic acid and MO (r = 0.745), and it was strong negative associated with LY (r = -0.535). These correlationships might be valuable for understanding the pathological process of anemia disease.

FIGURE 5

FIGURE 5. Correlation analysis between potential biomarkers and efficacy indicators of periphery blood parameters based on their Pearson correlation coefficient. Pm1, 4-phosphopantothenoylcysteine; Pm2, dIMP; Pm3, L-cystathionine; Pm4, 8-isoprostaglandin F2a; Pm5, 15-deoxy-d-12,14-PGJ2; Pm6, ubiquinone-2; Pm7, sphinganine; Pm8, sphingosine 1-phosphate; Pm9, deoxycholic acid; Pm10, pregnenolone; Pm11, 4-oxo-retinoic acid; Pm12, LysoPC (20:4); Pm13, 9-cis-retinoic acid; Pm14, Leucyl-phenylalanine; Pm15, Eicosapentaenoic acid; Pm16, docosahexaenoic acid; Pm17, 16-dehydroprogesterone; Pm18, S-adenosylhomocysteinase; WBC, white blood cell count; RBC, red blood cell count; HGB, hemoglobin; HCT, haematocrit; PLT, platelet count; LY, lymphocyte ratio; MO, monocyte ratio; NE, neutrophils ratio.

Results of Network Pharmacology

Screening of XSHG Active Compounds

Although the TCM prescriptions or the single herbs usually contained thousands of constituents, only a few with satisfying pharmacokinetic and pharmacodynamics features contributed to its curative effects. Our previous work has shown that major bioactive compounds of XSHG were aromatic acids, phthalides, alkaloids, flavonoids, and gingerols (Pang et al., 2016a). In present study, three key ADME parameters, including OB, DL and Caco-2, were applied to mine the major active components in XSHG. Additionally, several active constituents which did not satisfy these three criteria above were also introduced due to relative high contents or excellent bioactivities.

First of all, a gross of 250, 102, 378, 132, 378, 296, and 560 candidate ingredients were obtained from Danggui, Yimucao, Chuanxiong, Taoren, Honghua, Jiangtan, and Zhigancao, respectively. Combined with the ADME profiles and literature confirmation, some bioactive ingredients could be screened from these compounds. The general filtering criteria were: OB ≥ 30%, DL ≥ 0.18, and Caco-2 ≥-0.4. Consequently, a total of 58 potential constituents were selected for further analysis (Supplementary Table S2). To evaluate whether these constituents could be detected in the rat plasma after oral administration of XSHG prescriptions, 43 ingredients were characterized by UHPLC-QTOF-MS, which suggested that these screening compounds were believable to be applied in the network pharmacology analysis. The numbers of bioactive constituents were 22 for ASD and LC, 15 for LA, 7 for PP, 5 for CT, 4 for RZR, and 5 for RGP, respectively. Among the 58 ingredients, most of them had desirable pharmacokinetic properties. For instance, ferulic acid (M1, OB = 39.56%, Caco-2 = 0.47, and DL = 0.18) has antibacterial, anti-inflammatory, antioxidant and radio-protective effects (Mancuso and Santangelo, 2014). It should be pointed out that stachydrine held low OB values (OB = 0.27), but the content of stachydrine was found to be about 10 mg/g in XSHG (Pang et al., 2016a), which possessed potent anti-inflammatory, antiplatelet aggregation and antioxidant activities (Li B.H. et al., 2013). Similarly, the phthalides, the marker constituents in ASD and LC, showed low DL values, but they exhibited remarkable pharmacological activities (Wei et al., 2013). Thus, these compounds were also selected for further analysis.

Target Proteins of XSHG

Subsequently, a comprehensive in silico method was intended to construct the protein targets of the 58 compounds. Integrated models contained PharmMapper Sever, STITCH and SEA. Other databases were also used, including Google Scholar, DrugBank, Therapeutic Targets Database, and Herbal Ingredients’ Targets database. As a result, 46 candidate targets related to anemia were tentatively fished out, generating 373 target-ligand interactions: 168 for ASD and LC, 100 for LA, 39 for PP, 22 for CT, 18 for RZR, and 26 for RGP. The average number of candidate targets for each constituent was 6.4. Among the 58 active compounds, 7 had high degrees, including ferulic acid (degree = 17), ligustilide (degree = 16), stachydrine (degree = 15), leonurine (degree = 14), amygdalin (degree = 13), hydroxysafflor yellow A (degree = 11), glycyrrhizic acid (degree = 10), which showed the poly-pharmacology and multi-target properties of this prescription. At the same time, our previous studies exhibited that these compounds had larger areas under the curve compared to others (Pang et al., 2017). Based on the above analysis, these compounds undoubtedly played crucial roles in the treatment of anemia. Target information of the 58 bioactive constituents was listed in Table 4.

TABLE 4

TABLE 4. The protein targets information of 58 active constituents of XHSG.

Remarkably, multiple molecular targets were interacted with the bioactive constituents, such as TNF-α, IL-1β, MAPK1, MAPK10, IL-6, PTGS2, F2, HIF-1α, and EPOR. And the majority of 46 targets were directly related to anemia and vascular systems. For examples, TNF-α, IL-1β, and IL-6 could regulate the vascular systems to improve the vascular injury, and promote the hematopoietic function after postpartum hemorrhage (Janssens et al., 2015); MAPK1 and MAPK10, the mitogen-activated protein kinases, were involved in cell biological responses and immune defense by mediating multiple cellular processes (Mohanta et al., 2015); PTGS2 contributed to the development of thrombosis and atherosclerosis in response to the production of eicosanoids (Deeb et al., 2006). Moreover, stachydrine, leonurine, ferulic acid and Z-ligustilide might mediate the expressions of F2, HIF1α, and EPOR to improve blood circulation, thus exerting the effects of immune enhancement, erythrocyte synthesis and blood-vessel dilation (Bi et al., 2012). Senkyunolide I may interact with MAPK1, IL-1β, and TNF-α, which were also involved in anemia and vascular systems (Hu et al., 2016).

Except for β-sitosterol, nicotinic acid and vanillin, other bioactive constituents also interacted with protein targets relevant to pain, inflammation, microcirculation disorders, immune and nervous system diseases (Yue et al., 2017b). Five potential targets containing F2, FGG, MAPK14, TBXAS1, and ALOX5 were closely related to micro-circulation disorders (Yue et al., 2017b), which could explain why XSHG had the function of improving microcirculation. Besides, several targets such as CA2, HMOX1, GSTA1, and NOS2 were closely associated with anemia (Lin et al., 2013), which could exhibit significant influences on the therapeutic effects of XSHG.

C-T-P Network Analysis

To clarify the underlying mechanism of XSHG, 24 anemia-relevant signaling pathways were obtained from KEGG. Among these pathways, seven were consistent with the results of metabolomics. Then we mapped all potential ingredients with their corresponding targets onto these predictable pathways. And a global view of C-T-P network was generated by Cytoscape 3.4.0, which consisted of 128 nodes (58 ingredients, 46 genes and 24 pathways) and 453 edges (Figure 6A). Besides, the relationships between betweenness centrality and degree in the C-T-P network (Figure 6B) demonstrated that the distribution of betweenness centrality and degree was strongly correlated and the most highly connected vertices had high centrality scores. Most pathways were interacted with a few targets, but about 20% candidate targets were involved in multi-pathways (≥5), reflecting the crucial roles in treating anemia. For example, NOS3, which was closely associated with immune enhancement, vasodilation and neuroprotective effects (Shu et al., 2015), was influenced by six constituents through six pathways.

FIGURE 6

FIGURE 6. Compound-Target-Pathway network of active ingredients in XSHG for anemia disease (A); the correlation plots of the degree and betweenness centrality for this network (B). Fifty-eight active ingredients map 46 potential targets and 24 signaling pathways. Green nodes represent the ingredients, brown nodes represent the targets, and red nodes are pathways.

Moreover, we discovered that four biological pathways, i.e., HIF-1α, TNF, MAPK, and sphingolipid signaling pathways were remarkably enriched for the mined targets, implying the significant roles of these potential pathways. Among them, HIF-1α signaling pathway had the highest degree and the crucial targets in this pathway such as HMOX1, HIF-1α, F2, and IL6 were highly related to the hematopoietic function (Haase, 2013). Also, it could be observed that most compounds of XHSG were involved in these signaling pathways, providing a scientific basis for prevention and treatment of anemia. For example, 33 compounds such as Z-ligustilide, ferulic acid, and stachydrine were referred to mediating key targets in HIF-1α signaling pathway.

Furthermore, five biological pathways containing VEGF, Estrogen, Toll-like receptor, NF-kB, and hematopoietic cell lineage signaling pathways were the second major pathways able to exhibit anti-oxidative, anti-inflammatory, immune responses, and neuroprotective effects (DiDonato et al., 2012; Cabas et al., 2013; Maglione et al., 2015). It should be noted that different types of active ingredients could interact with various pathways. For alkaloids, VEGF, Estrogen and Toll-like receptor, signaling pathways exhibited higher correlations with their corresponding target proteins, whereas phthalides were mainly involved in HIF-1α and hematopoietic cell lineage pathways. These results indicated that various constituents of XHSG might act on different biological pathways.

Results of Experimental Validation

Selection of Crucial Targets

Integrating the results of metabolomics and network pharmacology, the metabolite-protein interaction network, containing five metabolites, six enzymes, six pathway-related proteins and eight interaction genes, was constructed (Figure 7). Then the crucial targets of shared pathways were firstly selected from the constructed network: coenzyme A biosynthesis (ACSS1 and COASY), cysteine and methionine metabolism (CBS and AHCY), and sphingolipid signaling pathway (SPHK1 and S1PR1). Moreover, other key targets, which were closely associated with the erythrocyte function (EPOR and HIF-1α) and blood microcirculation (TNF-α, IL6, and F2), were also applied for the experiment validation. As shown in Supplementary Figure S3, the correlation plots of the degree and betweenness centrality indicated that targets AHCY, CBS, EPOR, IL6, and TNF had the relatively high connected vertices compared with the average values, suggesting these targets had the significant roles.

FIGURE 7

FIGURE 7. Five candidate metabolites for anemia associated with relevant protein targets. Metabolites (green), enzymes (brown), pathway-related proteins (purple), and interaction genes (red).

ELISA Analysis

As shown in Figure 8, the levels of TNF-α, IL6, F2, and EPOR were significantly reduced (P < 0.01) in the model, presenting the inhibition of these cytokines. After oral administration of XSHG, the decreased levels of TNF-α, IL6, F2, and EPOR were obviously increased (P < 0.05). TNF-α could inhibit erythropoiesis through inducing dendritic cell proliferation (Mizrahi et al., 2013), and the plasma TNF-α level might show a compensatory decrease in anemia rats. IL6 had the function of hematopoiesis, immune response and inflammation (Janssens et al., 2015). XSHG could raise the level of IL6 and promote the hematopoietic and immune function. Moreover, XHSG could improve the erythrocyte immune functions and wound healing by enhancing the expression of EPOR and F2 in anemia.

FIGURE 8

FIGURE 8. Expression of EPOR, F2, TNF-α, and IL-6 protein levels in rat plasma using ELISA analysis. Data were presented as the mean ± SD (n = 6), ^ΔP < 0.05, ^ΔΔP < 0.01 vs. normal group; ▲P < 0.05, ▲▲P < 0.01 vs. model group.

Western Blot Analysis

As we all know, the liver serves as the main hematopoietic organ in mammals, so certain enzymes undoubtedly play significant roles in hematopoietic function. ACSS1 and COASY levels in liver were firstly investigated to explore the effects of XSHG on energy metabolism, and Western blot results were exhibited in Figure 9. As expected, anemia rats showed obvious inhibition of ACSS1 and COASY. However, XSHG enhanced the expression of these targets, indicating that it could improve energy metabolism through activating coenzyme A biosynthesis (Tomita et al., 2013). We also explored whether XSHG could mediate the cysteine and methionine pathway in this model. Clearly, the reduced levels of AHCY and CBS were observed in anemia. After oral administration of XSHG, the levels of these two protein targets were markedly increased. CBS and AHCY mainly participated in DNA methylation, which could promote the biosynthesis of erythrocyte (Amabile et al., 2015). Consequently, XSHG might up-regulate cysteine and methionine pathway to improve erythrocyte function.

FIGURE 9

FIGURE 9. Western blot analysis of ACSS1, COASY, AHCY, CBS, S1PR1, SPHK1, and HIF-1α protein levels in rat liver. Data were presented as the mean ± SD (n = 6), ^ΔP < 0.05, ^ΔΔP < 0.01 compared to normal group; while ▲P < 0.05, ▲▲P < 0.01 compared to the model group.

Furthermore, crucial targets of sphingolipid signaling pathway were investigated for the liver-protective effects of XSHG, and the expression of SPHK1 and S1PR1 was evaluated. The results indicated that XSHG could significantly reduce the expression of disturbed SPHK1 and S1PR1 to improve signal transduction (Harvald et al., 2015). Besides, to evaluate whether XSHG had the anti-anoxia effects, the expression of HIF-1α was determined in anemia rats. As expected, XSHG indeed down-regulated the protein level of HIF-1α, which could explain why HIF-1α pathway was blocked in XSHG group.

Discussion

According to the above results, we assumed that bioactive compounds of XSHG might simultaneously interact with multiple pathways like coenzyme A biosynthesis, sphingolipids and HIF-1α biological pathways, thereby exhibiting synergistic effects in anemia disease. In fact, lots of researches have indicated the significant roles of HIF-1α, coenzyme A biosynthesis and sphingolipids biological pathways in treating anemia and related diseases.

As is well known, in anemia patients, the balance of internal environment is weakened, thus the energy and substance metabolism is disturbed, and it will lead to the aggravation of anemia disease. Coenzyme A, a significant coenzyme participated in the biosynthesis of fatty acid and pyruvate, had the effects of immune enhancement and tissue repair, thus exerting a vital role in the energy metabolism of body (Tomita et al., 2013). Several coenzyme A biosynthesis target genes, including ACSS1, COASY, and PKLR, were particularly important for improving the disordered energy metabolism. By activating coenzyme A biosynthesis, disturbed energy metabolism could be improved, thus to enhance the synthesis of erythrocyte. Recent studies have demonstrated that Ge-Gen-Qin-Lian Decoction could regulate energy metabolism to treat metabolic diseases through activation of coenzyme A biosynthesis (Zhang et al., 2016), as confirmed by our results.

Moreover, after hemorrhage, the numbers of erythrocyte decreased significantly. Correspondingly, the levels of several cytokines such as EPOR, F2, and IL6 were reduced, leading to the occurrence of hypoxia. HIF-1α, a mainly mediator of maintaining the oxygen balance, could mediate the generation and survival of circulating erythrocyte (Rankin et al., 2012). Due to the increased expression of EPOR, the hypoxia could be improved, which could explain why the protein level of HIF-1α in XSHG group was down-regulated. By blocking HIF-1α signaling pathway, XSHG could maintain the balance of active oxygen metabolism, thereby promoting cell differentiation and enhancing the generation of EPOR. Some researchers have revealed that Herba Leonuri obviously increased the numbers of erythrocyte, attenuated the levels of EPOR and IL6, inhibited HIF-1α expression, and consequently improved the anemia, which supported its crucial role in the improvement of anemia (Li X. et al., 2013).

Sphingolipids, the critical composition of cell membrane, could regulate the cellular growth, differentiation and aging (Testai et al., 2014). They included sphingomyelin, cerebroside, and ganglioside, among which, sphingomyelin might play significant roles in the acute leukemia disease (Li et al., 2015). Some studies have indicated that sphingolipids were involved in inflammation, apoptosis, cellular immunity response (Testai et al., 2014). During the process of anemia, the membrane of erythrocyte was injured, resulting in erythrocyte apoptosis. So far, more and more evidence has shown that organisms could initiate erythrocyte apoptosis to maintain homeostasis of anemia patients (Knapp et al., 2012; Egler and Lang, 2017). Our research group have confirmed that herb pair Danggui-Chuanxiong could suppress erythrocyte apoptosis by inhibiting sphingolipids dependent signaling molecules (Li W.X. et al., 2014). Thus, the levels of sphingolipids decreased in anemia, which was consistent with the injured membrane of erythrocyte and disordered peripheral blood indicators.

In this work, an integrative approach of metabolomics and network pharmacology was performed to investigate the biological mechanisms of XSHG for treating anemia. Since the main ingredients of XSHG were aromatic acids, phthalides and alkaloids, we assumed that the treatments of XSHG might be ascribed to the integrate effects of multiple compounds rather than single constituent. The main action mechanisms of its therapeutic effects on anemia were shown in Figure 10. Collectively, the action mechanisms for treating anemia of XSHG were attributed to the activation of coenzyme A biosynthesis, inhibition of sphingolipid metabolism and suppression of HIF-1α pathway. And XSHG could play critical roles in anemia via increasing the expression of erythrocyte, improving the substance metabolism of liver, and promoting blood circulation. Further studies should evaluate the efficiency of the main active constituents such as stachydrine, ligustilide, and ferulic acid on other pathways and their interactions. This investigation was valuable for performing a systematical investigation of herb medicines, as well as for efficiently predicting the therapeutic targets of TCM formulae.

FIGURE 10

FIGURE 10. Hypothetical molecular mechanisms of XSHG for treating postpartum anaemia. Lines with arrow-heads represent activation, and lines with bars at the end denote inhibition.

Author Contributions

Y-PT and J-AD proposed the idea and designed this experiment. H-QP, Y-YC, Y-JT, and Y-JC performed the experiments. H-QP and S-JY designed the study. H-QP, X-QS, G-SZ, and AK participated in data analysis. S-JY, S-LH, Y-JS, JS, and Z-ST contributed to writing, revision and proof-reading this manuscript. All authors read and approved the final manuscript.

Funding

This work was supported by National Natural Science Foundation of China (81773882 and 81603258), and State Project of TCM standardization “standardization construction of Xin-Sheng-Hua granule” (ZYBZH-C-JS-34). This work was also supported by the Open Project Program of Jiangsu Key Laboratory for High Technology Research of TCM Formulae and Jiangsu Collaborative Innovation Center of Chinese Medicinal Resources Industrialization (No. FJGJS-2015-11), 333 High-level Talents Training Project funded by Jiangsu Province (No. BRA2016387), Technology Innovation Venture fund by Nanjing University of Chinese Medicine (CX201503). This research was also financially supported by a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD). We are also pleased to thank Waters Ltd. for technical support.

Conflict of Interest Statement

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.

The reviewer P-MA and handling Editor declared their shared affiliation.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2018.00165/full#supplementary-material

Footnotes

1. ^ http://www.stats.gov.cn/
2. ^ http://lsp.nwu.edu.cn/tcmsp.php
3. ^ http://stitch.embl.de/
4. ^ http://sea.bkslab.org/
5. ^ http://www.drugbank.ca/
6. ^ http://bidd.nus.edu.sg/group/ttd/
7. ^ http://www.pharmgkb.org
8. ^ http://ctdbase.org/
9. ^ https://massive.ucsd.edu/ProteoSAFe/static/massive.jsp

References

Amabile, G., Ruscio, D. A., Müller, F., Welner, R. S., Yang, H., Ebralidze, A. K., et al. (2015). Dissecting the role of aberrant DNA methylation in human leukemia. Nat. Commun. 6:7091. doi: 10.1038/ncomms8091

PubMed Abstract | CrossRef Full Text | Google Scholar

Balarajan, Y., Ramakrishnan, U., Ozaltin, E., Shankar, A. H., and Subramanian, S. V. (2011). Anaemia in low-income and middle-income countries. Lancet 378, 2123–2135. doi: 10.1016/S0140-6736(10)62304-5

CrossRef Full Text | Google Scholar

Bi, C. W., Xu, L., Zhang, W. L., Zhan, J. Y., Fu, Q., Zheng, K. Y., et al. (2012). Fo shou san, an ancient herbal decoction prepared from angelicae sinensis radix and chuanxiong rhizoma, induces erythropoietin expression: a signaling mediated by the reduced degradation of hypoxia-inducible factor in cultured liver cells. Planta Med. 78, 122–127. doi: 10.1055/s-00311-280356

PubMed Abstract | CrossRef Full Text

Cabas, I., Rodenas, M. C., Abellán, E., Meseguer, J., Mulero, V., and García-Ayala, A. (2013). Estrogen signaling through the G protein-coupled estrogen receptor regulates granulocyte activation in fish. J. Immunol. 191, 4628–4639. doi: 10.4049/jimmunol.1301613

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, C., Chen, J., Lyu, C., Yu, J., Zhao, W., Wang, Y., et al. (2015). Bioinformatics analysis of the effects of tobacco smoke on Gene expression. PLoS One 10:e0143377. doi: 10.1371/journal.pone.0143377

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, M. M., Yang, F. F., Yang, X. M., Lai, X. M., and Gao, Y. X. (2016). Systematic understanding of mechanism of a Chinese herbal formula in treatment of metabolic syndrome by an integrated pharmacology approach. Int. J. Mol. Sci. 17:E2114. doi: 10.3390/ijms17122114

PubMed Abstract | CrossRef Full Text | Google Scholar

Ciborowski, M., Lipska, A., Godzien, J., Ferrarini, A., Korsak, J., Radziwon, P., et al. (2012). Combination of LC-MS- and GC-MS-based metabolomics to study the effect of ozonated autohemotherapy on human blood. J. Proteome. Res. 11, 6231–6241. doi: 10.1021/pr3008946

PubMed Abstract | CrossRef Full Text | Google Scholar

Courant, F., Antignac, J. P., Dervilly-Pinel, G., and Le, B. B. (2014). Basics of mass spectrometry based metabolomics. Proteomics 14, 2369–2388. doi: 10.1002/pmic.201400255

PubMed Abstract | CrossRef Full Text | Google Scholar

Cui, D. N., Wang, X., Chen, J. Q., Lv, B., Zhang, P., Zhang, W., et al. (2017). Quantitative evaluation of the compatibility effects of Huangqin Decoction on the treatment of irinotecan-induced gastrointestinal toxicity using untargeted metabolomics. Front. Pharmacol. 8:211. doi: 10.3389/fphar.2017.00211

PubMed Abstract | CrossRef Full Text | Google Scholar

Deeb, R. S., Shen, H., Gamss, C., Gavrilova, T., Summers, B. D., Kraemer, R., et al. (2006). Inducible nitric oxide synthase mediates prostaglandin H2 synthase nitration and suppresses eicosanoid production. Am. J. Pathol. 168, 349–362. doi: 10.2353/ajpath.2006.050090

PubMed Abstract | CrossRef Full Text | Google Scholar

DiDonato, J. A., Mercurio, F., and Karin, M. (2012). NF-κB and the link between inflammation and cancer. Immunol. Rev. 246, 379–400. doi: 10.1111/j.1600-065X.2012.01099.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Duan, Y., Pei, K., Cai, H., Tu, S. C., Cheng, X. W., Zhang, Z. W., et al. (2016). Strategy of integrated evaluation on treatment of traditional Chinese medicine as interaction of system to system and establishment of novel fuzzy target contribution recognition with herb pairs, a case study on Astragali Radix-Fructus Corni. Mol. Cell. Endocrinol. 434, 219–237. doi: 10.1016/j.mce.2016.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Egler, J., and Lang, F. (2017). Triggering of eryptosis, the suicidal erythrocyte death, by perifosine. Cell Physiol. Biochem. 41, 2534–2544. doi: 10.1159/000475977

PubMed Abstract | CrossRef Full Text | Google Scholar

Froessler, B., Collingwood, J., Hodyl, N. A., and Dekker, G. (2014). Intravenous ferric carboxymaltose for anaemia in pregnancy. BMC Pregnancy Childbirth 14:115. doi: 10.1186/1471-2393-14-115

PubMed Abstract | CrossRef Full Text | Google Scholar

Godzien, J., Ciborowski, M., Angulo, S., Ruperez, F. J., Martínez, M. P., Señorans, F. J., et al. (2011). Metabolomic approach with LC-QTOF to study the effect of a nutraceutical treatment on urine of diabetic rats. J. Proteome. Res. 10, 837–844. doi: 10.1021/pr100993x

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Q. Y., Mao, X., Zhang, Y. Q., Meng, S. Q., Xi, Y., Ding, Y., et al. (2016). Guizhi-Shaoyao-Zhimu decoction attenuates rheumatoid arthritis partially by reversing inflammation-immune system imbalance. J. Transl. Med. 14:165. doi: 10.1186/s12967-016-0921-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Haase, V. H. (2013). Regulation of erythropoiesis by hypoxia-inducible factors. Blood Rev. 27, 41–53. doi: 10.1016/j.blre.2012.12.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Harvald, E. B., Olsen, A. S., and Fargeman, N. J. (2015). Autophagy in the light of sphingolipid metabolism. Apoptosis 20, 658–670. doi: 10.1007/s10495-015-1108-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffmann, T., Krug, D., Hüttel, S., and Müller, R. (2014). Improving natural products identification through targeted LC-MS/MS in an untargeted secondary metabolomics workflow. Anal. Chem. 86, 10780–10788. doi: 10.1021/ac502805w

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, Y. Y., Wang, Y., Liang, S., Yu, X. L., Zhang, L., Feng, L. Y., et al. (2016). Senkyunolide I attenuates oxygen-glucose deprivation/reoxygenation-induced inflammation in microglial cells. Brain Res. 1649, 123–131. doi: 10.1016/j.brainres.2016.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Janssens, K., Slaets, H., and Hellings, N. (2015). Immunomodulatory properties of the IL-6 cytokine family in multiple sclerosis. Ann. N. Y. Acad. Sci. 1351, 52–60. doi: 10.1111/nyas.12821

PubMed Abstract | CrossRef Full Text | Google Scholar

Keiser, M. J., Roth, B. L., Armbruster, B. N., Ernsberger, P., Irwin, J. J., and Shoichet, B. K. (2007). Relating protein pharmacology by ligand chemistry. Nat. Biotechnol. 25, 197–206. doi: 10.1038/nbt1284

PubMed Abstract | CrossRef Full Text | Google Scholar

Kell, D. B., and Goodacre, R. (2014). Metabolomics and systems pharmacology: why and how to model the human metabolic network for drug discovery. Drug Discov. Today 19, 171–182. doi: 10.1016/j.drudis.2013.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Knapp, M., Zendzian-Piotrowska, M., Błachnio-Zabielska, A., Zabielski, P., Kurek, K., and Górski, J. (2012). Myocardial infarction differentially alters sphingolipid levels in plasma, erythrocytes and platelets of the rat. Basic Res. Cardiol. 107:294. doi: 10.1007/s00395-012-0294-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuhn, M., Szklarczyk, D., Franceschini, A., Mering, C. V., Jensen, L. J., and Bork, P. (2012). STITCH 3: zooming in on protein-chemical interactions. Nucleic Acids Res. 40, D876–D880. doi: 10.1093/nar/gkr1011

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, B. H., Wu, J. D., and Li, X. L. (2013). Simultaneous determination and pharmacokinetic study of stachydrine and leonurine in rat plasma after oral administration of Herba Leonuri extract by LC–MS/MS. J. Pharmaceut. Biomed. 76, 192–199. doi: 10.1016/j.jpba.2012.12.029

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, L., Wang, L. L., Shangguan, D. H., Wei, Y. B., Han, J. J., Xiong, S. X., et al. (2015). Ultra-high-performance liquid chromatography electrospray ionization tandem mass spectrometry for accurate analysis of glycerophospholipids and sphingolipids in drug resistance tumor cells. J. Chromatogr. A 1381, 140–148. doi: 10.1016/j.chroma.2015.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S., Fan, T. P., Jia, W., Lu, A., and Zhang, W. (2014). Network pharmacology in traditional Chinese medicine. Evid. Based Complement. Alternat. Med. 2014:138460. doi: 10.1155/2014/138460

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, W. X., Tang, Y. P., Guo, J. M., Shang, E. X., Qian, Y. F., Wang, L. Y., et al. (2014). Comparative metabolomics analysis on hematopoietic functions of herb pair Gui-Xiong by ultra-high-performance liquid chromatography coupled to quadrupole time-of-flight mass spectrometry and pattern recognition approach. J. Chromatogr. A 1346, 49–56. doi: 10.1016/j.chroma.2014.04.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, X., Wang, B., Li, Y., Wang, L., Zhao, X., Zhou, X., et al. (2013). The Th1/Th2/Th17/Treg paradigm induced by stachydrine hydrochloride reduces uterine bleeding in RU486-induced abortion mice. J. Ethnopharmacol. 145, 241–253. doi: 10.1016/j.jep.2012.10.059

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C. Y., Hsiao, W. C., Huang, C. J., Kao, C. F., and Hsu, G. S. (2013). Heme oxygenase-1 induction by the ROS-JNK pathway plays a role in aluminum-induced anaemia. J. Inorg. Biochem. 128, 221–228. doi: 10.1016/j.jinorgbio.2013.07.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H. P., Zeng, L. T., Yang, K. L., and Zhang, G. M. (2016). A network pharmacology approach to explore the pharmacological mechanism of Xiaoyao Powder on an ovulatory infertility. Evid. Based Complement. Alternat. Med. 2016:2960372. doi: 10.1155/2016/2960372

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X. F., Ouyang, S. S., Yu, B., Liu, Y. B., Huang, K., Gong, J. Y., et al. (2010). PharmMapper server: a web server for potential drug target identification using pharmacophore mapping approach. Nucleic Acids Res. 38, 609–614. doi: 10.1093/nar/gkq300

PubMed Abstract | CrossRef Full Text | Google Scholar

Maglione, P. J., Simchoni, N., and Cunningham-Rundles, C. (2015). Toll-like receptor signaling in primary immune deficiencies. Ann. N. Y. Acad. Sci. 1356, 1–21. doi: 10.1111/nyas.12763

PubMed Abstract | CrossRef Full Text | Google Scholar

Mancuso, C., and Santangelo, R. (2014). Ferulic acid: pharmacological and toxicological aspects. Food Chem. Toxicol. 65, 185–195. doi: 10.1016/j.fct.2013.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Markova, V., Norgaard, A., Jørgensen, K. J., and Langhoff-Roos, J. (2015). Treatment for women with postpartum iron deficiency anaemia. Cochrane Database Syst. Rev. 8:CD010861. doi: 10.1002/14651858.CD010861.pub2

PubMed Abstract | CrossRef Full Text | Google Scholar

Milman, N. (2011). Postpartum anemia I: definition, prevalence, causes, and consequences. Ann. Hematol. 90, 1247–1253. doi: 10.1007/s00277-011-1279-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Milman, N. (2012). Postpartum anemia II: prevention and treatment. Ann. Hematol. 91, 143–154. doi: 10.1007/s00277-011-1381-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Mizrahi, K., Stein, J., Yaniv, I., Kaplan, O., and Askenasy, N. (2013). TNF-α has tropic rather than apoptotic activity in human hematopoietic progenitors: involvement of TNF receptor-1 and caspase-8. Stem Cells 31, 156–166. doi: 10.1002/stem.1259

PubMed Abstract | CrossRef Full Text | Google Scholar

Mohanta, T. K., Arora, P. K., Mohanta, N., Parida, P., and Bae, H. (2015). Identification of new members of the MAPK gene family in plants shows diverse conserved domains and novel activation loop variants. BMC Genomics 16:58. doi: 10.1186/s12864-015-1244-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Mousa, H. A., Blum, J., Abou, E. I., Senoun, G., Shakur, H., and Alfirevic, Z. (2014). Treatment for primary postpartum haemorrhage. Cochrane Database Syst. Rev. 2:CD003249 doi: 10.1002/14651858.CD003249.pub3

PubMed Abstract | CrossRef Full Text | Google Scholar

Naz, S., Vallejo, M., García, A., and Barbas, C. (2014). Method validation strategies involved in non-targeted metabolomics. J. Chromatogr. A 1353, 99–105. doi: 10.1016/j.chroma.2014.04.071

PubMed Abstract | CrossRef Full Text | Google Scholar

Pacana, T., Cazanave, S., Verdianelli, A., Patel, V., Min, H. K., Mirshahi, F., et al. (2015). Deregulated hepatic methionine metabolism drives homocysteine elevation in diet-induced nonalcoholic fatty liver disease. PLoS One 10:e0136822. doi: 10.1371/journal.pone.0136822

PubMed Abstract | CrossRef Full Text

Pang, H. Q., Tang, Y. P., Cao, Y. J., Tan, Y. J., Jin, Y., Shi, X. Q., et al. (2017). Comparatively evaluating the pharmacokinetic of fifteen constituents in normal and anaemia rats after oral administration of Xin-Sheng-Hua Granule by UPLC-MS/MS. J. Chromatogr. B 106, 372–381. doi: 10.1016/j.jchromb.2017.07.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, H. Q., Wang, J., Tang, Y. P., Wu, L., Xu, H. Q., Ji, Y., et al. (2016a). Comparatively evaluating effect contribution of promoting blood circulation of herb pairs containing Angelicae Sinensis Radix on Xin-Sheng-Hua granule by withdrawal analysis. Chin. J. Chin. Mater. Med. 41, 4006–4014. doi: 10.4268/cjcmm20162120

PubMed Abstract | CrossRef Full Text | Google Scholar

Pang, H. Q., Wang, J., Tang, Y. P., Xu, H. Q., Wu, L., Jin, Y., et al. (2016b). Comparative analysis of the main bioactive components of Xin-Sheng-Hua granule and its single herbs by ultrahigh performance liquid chromatography with tandem mass spectrometry. J. Sep. Sci. 39, 4096–4106. doi: 10.1002/jssc.201600606

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, S., Sadanala, K. C., and Kim, E. K. (2015). A Metabolomic approach to understanding the metabolic link between obesity and diabetes. Mol. Cells 38, 587–596. doi: 10.14348/molcells.2015.0126

PubMed Abstract | CrossRef Full Text | Google Scholar

Pietrocola, F., Galluzzi, L., Bravo-San, P. J. M., Madeo, F., and Kroemer, G. (2015). Acetyl coenzyme A: a central metabolite and second messenger. Cell Metab. 21, 805–821. doi: 10.1016/j.cmet.2015.05.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Prosser, G. A., Larrouy-Maumus, G., and de Carvalho, L. P. (2014). Metabolomic strategies for the identification of new enzyme functions and metabolic pathways. EMBO Rep. 15, 657–669. doi: 10.15252/embr.201338283

PubMed Abstract | CrossRef Full Text | Google Scholar

Rankin, E. B., Wu, C., Khatri, R., Wilson, T. L., Andersen, R., Araldi, E., et al. (2012). The HIF signaling pathway in osteoblasts directly modulates erythropoiesis through the production of EPO. Cell 149, 63–74. doi: 10.1016/j.cell.2012.01.051

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, X. Q., Tang, Y. P., Zhu, H. X., Li, W. X., Li, Z. H., Li, W., et al. (2014a). Comparative tissue distribution profiles of five major bio-active components in normal and anaemia rats after oral administration of Danggui Buxue Decoction by UPLC-TQ/MS. J. Pharm. Biomed. Anal. 88, 207–215. doi: 10.1016/j.jpba.2013.08.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Shi, X. Q., Tang, Y. P., Zhu, H. X., Li, W. X., Li, Z. H., Li, W., et al. (2014b). Pharmacokinetic comparison of seven major bio-active components in normal and anaemia rats after oral administration of Danggui Buxue decoction by UPLC-TQ/MS. J. Ethnopharmacol. 153, 169–177. doi: 10.1016/j.jep.2014.02.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Shu, X., Keller, T. C., Begandt, D., Butcher, J. T., Biwer, L., Keller, A. S., et al. (2015). Endothelial nitric oxide synthase in the microcirculation. Cell Mol. Life Sci. 72, 4561–4575. doi: 10.1007/s00018-015-2021-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Takayama, T., Mochizuki, T., Todoroki, K., Min, J. Z., Mizuno, H., Inoue, K., et al. (2015). A novel approach for LC-MS/MS-based chiral metabolomics fingerprinting and chiral metabolomics extraction using a pair of enantiomers of chiral derivatization reagents. Anal. Chim. Acta 898, 73–84. doi: 10.1016/j.aca.2015.10.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Testai, F. D., Kilkus, J. P., Berdyshev, E., Gorshkova, I., Natarajan, V., and Dawson, G. (2014). Multiple sphingolipid abnormalities following cerebral micro endothelial hypoxia. J. Neurochem. 131, 530–540. doi: 10.1111/jnc.12836

PubMed Abstract | CrossRef Full Text | Google Scholar

Tomita, H., Imanaka, T., and Atomi, H. (2013). Identification and characterization of an archaeal ketopantoate reductase and its involvement in regulation of coenzyme A biosynthesis. Mol. Microbiol. 90, 307–321. doi: 10.1111/mmi.12363

PubMed Abstract | CrossRef Full Text | Google Scholar

Tunçalp, O., Souza, J. P., Gülmezoglu, M., and World Health Organization (2013). New WHO recommendations on prevention and treatment of postpartum hemorrhage. Int. J. Gynaecol. Obstet. 123, 254–256. doi: 10.1016/j.ijgo.2013.06.024

PubMed Abstract | CrossRef Full Text | Google Scholar

US National Research Council (1996). Guide for the Care and Use of Laboratory Animals. Washington, DC: The National Academies Press. doi: 10.17226/5140

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, P., Sun, H., Lv, H., Sun, W., Yuan, Y., Han, Y., et al. (2010). Thyroxine and reserpine-induced changes in metabolic profiles of rat urine and the therapeutic effect of Liu Wei Di Huang Wan detected by UPLC-HDMS. J. Pharm. Biomed. Anal. 53, 631–645. doi: 10.1016/j.jpba.2010.04.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Yu, S. Y., Jia, Q., Chen, L. C., Zhong, J. Q., Pan, Y. H., et al. (2017). Niao Du Qing granules relieve chronic kidney disease symptoms by decreasing renal fibrosis and anaemia. Oncotarget 8, 55920–55937. doi: 10.18632/oncotarget.18473

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. J., Zhang, A. H., Sun, H., Han, Y., and Yan, G. L. (2016a). Discovery and development of innovative drug from traditional medicine by integrated chinmedomics strategies in the post-genomic era. Trac Trend Anal. Chem. 76, 88–94. doi: 10.1016/j.trac.2015.11.010

CrossRef Full Text | Google Scholar

Wang, X. J., Zhang, A. H., Zhou, X. H., Liu, Q., Nan, Y., Guan, Y., et al. (2016b). An integrated chinmedomics strategy for discovery of effective constituents from traditional herbal medicine. Sci. Rep. 6:18997. doi: 10.1038/srep18997

PubMed Abstract | CrossRef Full Text | Google Scholar

Wei, Y., Huang, W. W., and Gu, Y. X. (2013). Online isolation and purification of four phthalide compounds from Chuanxiong rhizoma using high-speed counter-current chromatography coupled with semi-preparative liquid chromatography. J. Chromatogr. A 1284, 53–58. doi: 10.1016/j.chroma.2013.01.103

PubMed Abstract | CrossRef Full Text | Google Scholar

Wettersten, H. I., and Weiss, R. H. (2013). Applications of metabolomics for kidney disease research: from biomarkers to therapeutic targets. Organogenesis 9, 11–18. doi: 10.4161/org.24322

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, W., Chen, D. Y., Wang, N., Zhang, T., Zhou, R. K., Huan, T., et al. (2015). Development of high-performance chemical isotope labeling LC–MS for profiling the human fecal metabolome. Anal. Chem. 87, 829–836. doi: 10.1021/ac503619q

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Y. F., Li, Y., Wang, J. H., Sun, K., Tao, W. Y., Wang, Z. Z., et al. (2017). Systematic investigation of Ginkgo Biloba Leaves for treating cardio-cerebrovascular diseases in an animal model. ACS Chem. Biol. 12, 1363–1372. doi: 10.1021/acschembio.6b00762

PubMed Abstract | CrossRef Full Text | Google Scholar

Ye, H., Ye, L., Kang, H., Zhang, D. F., Tao, L., Tang, K. L., et al. (2011). HIT: linking herbal active ingredients to targets. Nucleic Acids Res. 39, 1055–1059. doi: 10.1093/nar/gkq1165

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, S.-J., Liu, J., Feng, W.-W., Zhang, F.-L., Chen, J.-X. Xin, L.-T., et al. (2017a). System pharmacology-based dissection of the synergistic mechanism of Huangqi and Huanglian for diabetes mellitus. Front. Pharmacol. 8:694. doi: 10.3389/fphar.2017.00694

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, S. J., Xin, L. T., Fan, Y. C., Li, S. J., Tang, Y. P., Duan, J. A., et al. (2017b). Herb pair Danggui-Honghua: mechanism underlying blood stasis syndrome by system pharmacology approach. Sci. Rep. 7:40318. doi: 10.1038/srep40318

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q. Y., Xu, G. L., Li, J., Guo, X. F., Wang, H., Li, B. T., et al. (2016). Metabonomic study on the plasma of streptozotocin-induced diabetic rats treated with Ge Gen Qin Lian Decoction by ultra high performance liquid chromatography–mass spectrometry. J. Pharm. Biomed. Anal. 120, 175–180. doi: 10.1016/j.jpba.2015.12.031

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, J. W., Li, X., Deng, P. Y., Wei, Y., Liu, J., Chen, M., et al. (2017). Chinese herbal formula, modified Danggui Buxue Tang, attenuates apoptosis of hematopoietic stem cells in immune-mediated aplastic anemia mouse model. J. Immunol. Res. 2017:9786972. doi: 10.1155/2017/9786972

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, F., Shi, Z., Qin, C., Tao, L., Liu, X., Xu, F., et al. (2012). Therapeutic target database update 2012: a resource for facilitating target-oriented drug discovery. Nucleic Acids Res. 40, 1128–1136. doi: 10.1093/nar/gkr797

PubMed Abstract | CrossRef Full Text | Google Scholar