Comprehensive Metabolomics and Network Pharmacology to Explore the Mechanism of 5-Hydroxymethyl Furfural in the Treatment of Blood Deficiency Syndrome

Radix Rehmanniae (RR, from Radix Rehmanniae (Gaertn.) DC.) is a natural medicine used in traditional Chinese medicine (TCM) since ancient times for the treatment of blood disorders. RR is steamed to get Rehmanniae Radix Praeparata (RP), which has a tonic effect on blood; the content of 5-hydromethylfurfural (5-HMF) increases more than four times after steaming. Studies have shown that 5-HMF has positive pharmacological effects on cardiovascular and hematological disorders. This study aimed to explore and verify the impact of 5-HMF on rats with chemotherapy-induced blood deficiency syndrome (BDS). Rats were given cyclophosphamide (CP) and acetophenhydrazine (APH) to induce BDS, the coefficients of some organs (liver, spleen, and kidney) were measured, and a routine blood test examined the coefficients of several peripheral blood cells. Metabolomics and network pharmacology were combined to find important biomarkers, targets, and pathways. Western blot was used to detect the expression of CYP17A1 and HSD3B1 proteins in the spleen. All these findings suggested that the 5-HMF significantly increased the number of peripheral blood cells and reversed splenomegaly in rats. In addition, 5-HMF upregulated CYP17A1 and HSD3B1 protein expression in splenic tissues. Also, 5-HMF ameliorated chemotherapy-induced BDS in rats, and its therapeutic mechanism might depend on steroid hormone biosynthesis and other pathways. It acts on blood deficiency via multiple targets and pathways, which is unique to Chinese medicine.

1 Introduction

According to the theory of traditional Chinese medicine (TCM), blood deficiency is a common disease and a pathological state of blood dysfunction and organ malnutrition (Li et al., 2015). Blood deficiency syndrome (BDS) is mainly related to excessive blood loss (Shi et al., 2014), deficient spleen and stomach function, insufficient hematogenesis, and blood stasis (Yong et al., 2012). The main diagnostic index is the decrease in blood cells or hemoglobin content, similar to anemia (Wang et al., 2013). Modern medicine holds that patients with BDS often have clinical manifestations such as impaired hematopoietic function, decreased visceral function, malnutrition, and bone marrow suppression (Zhang et al., 2014). BDS can cause spleen enlargement (Zhang et al., 2020). After chemotherapy, most patients have a hemoglobin synthesis disorder, which causes BDS and severely interferes with blood system circulation (Jiang et al., 2013). However, anemia is defined only as a decrease in the concentration of hemoglobin in the body, which is usually related to a reduction in red blood cells (Li et al., 2017).During cancer treatment, the standard treatment for BDS is iron and vitamin B12 supplementation to promote the synthesis of erythropoietin, which can accelerate the recovery of the hematopoietic system (Brown 2020). However, the application of these regimens is limited because of their unstable efficacy, high cost, and side effects (Milano and Schneider 2007).

A main component of TCM is 5-HMF, among which the most representative herb is Rehmanniae Radix Praeparata (RP) (Won et al., 2014). Rehmanniae Radix (RR) is steamed to obtain RP, which is a plant of the Scrophulariaceae family. RP has been used to treat BDS for thousands of years. It is a commonly used TCM in Asia (Li et al., 2021). The Compendium of Materia Medica records that RP is used as a blood tonic (Guo et al., 2016), which regulates and fortifies blood and treats BDS (Zee-Cheng 1992). The main components of RR are monosaccharides, oligosaccharides (Kitagawa et al., 1995), iridoids (Xia et al., 2020), and glycosides (Xu et al., 2012). The decomposition products of these compounds constitute RP, of which 5-HMF is the most prominent (Won et al., 2014). During the processing of RR, stachyose is reduced to hexose, and these components finally form 5-HMF, which increases the content of 5-HMF in RP after steaming by more than four times compared with RR (Zhu et al., 2007). Therefore, 5-HMF has become the main active component of RP and is used as a marker compound for quality control (Li et al., 2005; Korean Food and Drug Administration, 2008). This transformation of active ingredients in raw materials and processed products leads to differences in their efficacy. For example, RP is better at supplementing blood (Xia et al., 2020). According to current research results, 5-HMF has positive pharmacological effects on cardiovascular and blood diseases. For example, 5-HMF can increase hemoglobin’s oxygen affinity to support cardiac function during severe hypoxia (Lucas et al., 2019). It also has sound biological effects, such as antioxidant activity (Li et al., 2011) and strengthening the resistance of sickle red blood cells to injury caused by hypoxia (Qiang et al., 2021). Therefore, we speculate that 5-HMF can be used as a potential therapeutic drug for BDS. However, the treatment of BDS by 5-HMF has not been documented, and its therapeutic mechanism is unclear.

Network pharmacology and metabolomics are practical tools for elucidating the underlying mechanisms of TCM (Pan et al., 2020). First, a blood deficiency rat model induced by chemotherapy was used to evaluate the effect of 5-HMF on BDS. Second, a network pharmacology method was established to explore the mechanism of 5-HMF in treating BDS from the perspective of targets and pathways. At the same time, plasma metabolomics revealed the synergistic metabolic mechanism of metabolites and metabolic pathways. Finally, Western blot validated the metabolic pathway. This study was the first to use a combination of network pharmacology and plasma metabolomics to determine the mechanism of 5-HMF in the treatment of BDS. Hopefully, these results will provide a theoretical basis for elucidating the mechanism of 5-HMF in BDS treatment.

2 Materials and Methods

2.1 Animals and Ethic Statement

We used 40 healthy male SD rats weighing 220–260 g, batch number SCXK (Liao) 2020–0001, from Liaoning Changsheng Biotechnology Co., Ltd. The animal experiment was approved by the Ethics Committee of Heilongjiang University of Traditional Chinese Medicine (approval number 2019121101). Animals were kept in a room at 25 C and 40–60% humidity. Rats were randomly divided into five groups with 8 rats in each group: control group, model group, 5-HMF-L group, 5-HMF-M group, and 5-HMF-H group.

2.2 Chemicals and Reagents

The 5-HMF was purchased from Aladdin (batch number H2020020, >95% (Gas chromatography), containing 3–5% water as a stabilizer). N-acetophenazine was purchased from Aladdin (batch number L1909182), and cyclophosphamide (CP) was purchased from MACKLIN (batch number C11147187).

2.3 Establishment of the Blood Deficiency Syndrome Model

The establishment of the BDS animal model was consistent with previous studies (Li et al., 2015). APH (20 mg/kg, 10 mg/kg) was injected subcutaneously on the first and fourth days. On the 4th day, 2 h after subcutaneous injection of APH, CP (20 mg/kg) was injected intraperitoneally. CP was injected for 4 consecutive days after that. On the 7th day, the indexes of red blood cells (RBC), white blood cells (WBC), hemoglobin (HGB), and hematocrit (HCT) in peripheral blood were measured. The control group was injected with the same volume of saline in the same way. From the first day, the treatment group was given 5-HMF (5.19 mg/kg, 2.595 mg/kg, 1.2975 mg/kg) by gavage (Duan 2021), respectively, and the control group and model group were given the same volume of saline by gavage for 15 days.

2.4 Routine Blood Test

After the last administration, blood was taken from the abdominal aorta and collected in a sterile vacuum blood collection tube containing ethylenediaminetetraacetic acid (EDTA). The blood was analyzed by the HEMAVET 950 automatic hematology analyzer (Drew Scientific Group, Dallas, Texas, United States) to quantify RBC, WBC, HGB, and HCT.

2.5 Network Analysis

2.5.1 Target Prediction

The targets of 5-HMF were predicted in the following ways: Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP) (Ru et al., 2014), Swiss Target Prediction (Gfeller et al., 2014), Targetnet (Yao et al., 2016), PharmMapper (Wang et al., 2017). We used GeneCards (Stelzer et al., 2016), Online Mendelian Inheritance in Man (OMIM) (Amberger et al., 2015), Drugbank (Wishart et al., 2018) to identify targets related to BDS. In the OMIM database, we opened the Gene Map and entered “blood deficiency syndrome” to obtain 532 essential genes. We entered “blood deficiency syndrome” in the GeneCard database, selected genes with a relevance score >20, and obtained a total of 1,542 critical genes. We entered “anemia” in the Drugbank database and got a total of 11 essential genes and then removed the duplicated genes.

2.5.2 Kyoto Encyclopedia of Genes and Genomes and Gene Ontology Enrichment Analysis

In analyzing the common targets of the 5-HMF action target and BDS-related targets, the server Metascape (Zhou et al., 2019) was used for Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) enrichment analysis.

2.5.3 Network Construction

Two visualization networks were constructed: 1) Target point network, a network composed of 5-HMF, BDS, and common targets; 2) Targetpath network, in which the targets and their pathways were used to generate spatiotemporal networks. All networks were constructed by Cytoscape 3.8.0 (Liu et al., 2020).

2.6 LC-MS Metabolomics Analysis

2.6.1 Sample Preparation

After the final administration, the blood of rats was collected from the abdominal aorta, left to stand for 30 min at 4 C, centrifuged at 3,500 rpm/min for 15 min. Then the plasma supernatant sample was mixed 1:3 with acetonitrile, vortexed for 30 s, and left to stand for 30 min at −20°C. The blood was centrifuged at 13,500 rpm/min for 15 min at 4 C, 300 μL of supernatant was evaporated, and 150 μL of 50% acetonitrile was added to redissolve the sample. After centrifugation, the sample was injected into the injection cup for testing (Gong et al., 2019).

2.6.2 UPLC-QTOF/MS Analysis

The G2-Si (Waters Q-TOF SYNAPT from Manchester Waters Company) has a continuous mode and is equipped with an electrospray ion source. The plasma samples were chromatographed on an Acquity UPLC HSS T3 column with a gradient elution program at the flow rate of 0.2 μL/min. The mobile phase was comprised of water and 0.1% formic acid (phase A) and acetonitrile and 0.1% formic acid (phase B). The elution gradient of phase B was as follows: 2% (0–0.5 min), 2% up to 40% (0.5–9 min), 40% up to 98% (9–16 min), then maintained at 98% for 1 min. Optimized parameters for the mass spectrometer were as follows: ion spray voltage, 5 or −4 kV (for positive and negative mode, respectively); declustering potential, 60 V; curtain gas, 25 psi; nebulizer gas of 40 psi; interface heater temperature, 650 C; scan range (m/z), 50–1,200. The stability of the analysis was continuously monitored by analyzing QC samples at intervals of every 10 samples.

2.6.3 Multivariate Data Analysis

The advanced Progressive QI platform (Waters, United States) was used to obtain potential and critical annotations from the original centroid data file obtained by SYNAPT G2-Si. Metabolic data completed by principal component analysis and dimensionality reduction had the most significant impact on health and disease. The potential differences in metabolites were identified from the orthogonal partial least squares discriminant analysis (OPLS-DA). The components that significantly impacted the grouping were screened out through objective numerical values. The t-test value was lower than 0.05.

2.6.4 Biomarker Identification and Metabolic Pathway Analysis

LC-MS/MS data were imported into Ezinfo software (Waters Progressional QI) for principal component analysis and OPLS-DA data analysis. The variable projection importance value obtained by OPLS-DA was used to evaluate the impact of each metabolite difference on the discrimination and interpretation ability of each sample. The greater the VIP value, the more significant the contribution of metabolites to sample differentiation. It is generally believed that there are significant differences between variables with VIP >1. The potential biomarkers were the S-plot constructed by Ezinfo, and the V-plot was labeled.

2.7 Statistical Analysis

All data were expressed in terms of mean ± standard deviation. The two-tailed unpaired t-test was performed using SPSS 21.0 software (from SPSS Corporation, United States). A value of p < 0.05 was considered to show a significant difference. A value of p < 0.01 was considered to indicate a highly significant difference. The histogram was drawn by GraphPad Prism 7 software (GraphPad Software, United States).

2.7.1 Pathway Analysis

The targets from network pharmacology and the metabolites from plasma metabolomics were jointly analyzed by MetaboAnalyst to select crucial metabolism pathways.

2.8 Western Blot Was Applied to Verify the Expression of Common Key Proteins in Metabolomics and Network Pharmacology

The bincinchoninic acid (BCA) method was used to determine the total protein of conventionally extracted splenocytes. The extract was separated by 10% SDS polyacrylamide gel electrophoresis and transferred to the PVDF membrane. It was blocked with 5% skim milk at room temperature for 1 h, and antibodies of CYP17A1, HSD3B1 (all 1:1,000), and GAPDH (all 1:2,000) were added. These were incubated overnight at 4 C and rinsed with tris buffered saline with Tween 20 (TBST) three times for 10 min each. The membrane was immersed in the horseradish peroxidase-labeled secondary antibody (1:1,000), diluted in 2% skimmed milk, incubated at room temperature for 1 h, and rinsed with TBST 3 times, each for 10 min. Image Quant LAS 500 imaging equipment acquired the signal using the ECL chemiluminescence method.

3 Results

3.1 Behavioral Analysis of Rats

Changes in the general behavior of rats could reflect the occurrence of and recovery from BDS in rats. After BDS induction, rats in the model group developed fatigue and lethargy, accompanied by weight loss, thinning hair, pale ears and tails, and loss of appetite. These symptoms were consistent with the description of BDS in Chinese medicine. In contrast, the control and BDS-induced rats treated with 5-HMF were relatively robust, showing thick and shiny hair, pink and moist nose and lips, round and pink tails, and stable body weight along with increased appetite.

3.2 Routine Blood Testing

After 15 days of administration of normal saline and 5-HMF, the RBC, WBC, HGB, and HCT levels in the peripheral blood of the rats were measured (Figure 1). Compared with the control group, the RBC, WBC, HGB, and HCT levels in the model group were significantly reduced (p < 0.01), indicating that the blood deficiency model was successfully established. The levels of WBC, RBC, HGB, and HCT in the Model group were significantly increased after 5-HMF treatment (p < 0.05 or p < 0.01).

FIGURE 1

FIGURE 1. The blood parameters of rats in the control and model groups treated with saline and different doses of 5-HMF. (A) The level of WBC. (B) The level of RBC. (C) The level of HGB. (D) The level of HCT. Each value represents the mean SD (n = 8); ^#p < 0.05 and ^##p < 0.01, compared with the control group; *p < 0.05 and **p < 0.01, compared with the model group.

3.3 Organ Index

The results of the organ index are shown in Figure 2 Compared with the control group, the liver and kidney function in the model group were relatively unchanged, but the spleen was abnormally enlarged. Except for the spleen coefficient, the experimental group (model group, 5-HMF-H group, 5-HMF-M group, and 5-HMF-L group) and the control group had no significant differences in other organ coefficients. Compared with the control group, the spleen coefficient in the model group increased significantly (p < 0.01). Compared with the model group, the spleen coefficient in the 5-HMF-L group was significantly reduced (p < 0.01), the spleen coefficient in the 5-HMF-M group was decreased significantly (p < 0.05), and the spleen coefficient in the 5-HMF-H group was not very different. The results showed that 5-HMF had a dose-dependent reversal effect on spleen enlargement.

FIGURE 2

FIGURE 2. Organ index of mice in the control group and model group treated with saline (model) and different doses of 5-HMF. (A) Liver index. (B) Kidney index. (C) Spleen index. Each value represents the mean SD (n = 8); ^##p < 0.01, compared with the control group; *p < 0.05 and **p < 0.01, compared with the model group.

3.4 Results of Network Pharmacology Analysis

3.4.1 5-Hydromethylfurfural and Blood Deficiency Syndrome Common Target Interaction

In this study, four online websites were used to predict 5-HMF targets, and three were used to predict BDS disease targets. Finally, it was determined that 5-HMF had 494 targets, BDS had a total of 1887 targets, and a total of 84 common targets. To further explore the mechanism of 5-HMF treatment of BDS, KEGG enrichment analysis was performed on 84 common targets, and 14 metabolic pathways were obtained (Figure 3), and the common targets and pathways are illustrated in a network diagram (Figure 4). The results show that the 5-HMF treatment of BDS has the characteristics of multiple pathways and multiple targets.

FIGURE 3

FIGURE 3. The blue triangle node represents the common pathway of 5-HMF and BDS, while the red quadrilateral node represents the common target of 5-HMF and BDS.

FIGURE 4

FIGURE 4. (A) 5-HMF and BDS common target KEGG enrichment analysis; (B) 5-HMF and BDS common target KEGG enrichment analysis bubble chart; (C) 5-HMF and BDS common target GO enrichment analysis of dots; (D) Bubble chart of GO enrichment analysis of the common target of 5-HMF and BDS.

3.5 Results of Metabolomics Analysis

3.5.1 Multivariate Data Analysis

The samples were analyzed with the UPLC/Q-TOF metabolic spectrum. No abnormalities were found in the total ion chromatogram, and the retention time of chromatographic peaks and their ionic strength were observed (Figure 5) Masslynx and Progenesis QI software were used to analyze many raw files from G2-Si total quantification for data preprocessing, such as peak picking, denoising, and normalization, to generate a bioinformatics matrix. Then we exported the data set to Simca 14.1 software for chemometric analysis. The principal component analysis score chart showed that the model group was separated from the control group (Figure 6A). A significant separation between the model group and the control group was observed in the OPLS-DA chart (Figure 6B). The S-plot of OPLS-DA showed multiple metabolites (Figure 6C). The V-plot diagram of OPLS-DA showed a remarkable diversity of various metabolites (Figure 6D). The replacement test was used to verify the OPLS-DA mode (Figure 6E). Figure 6F shows that all experimental groups were separate in the liquid chromatography-mass spectrometry metabolic map. Among them, the metabolic profiles of the 5-HMF-H group, 5-HMF-M group, and 5-HMF-L group were closer to the control group than the model group, suggesting that the model’s metabolic disorder that was induced was reversed after drug treatment, and the 5-HMF-L group had a better separation effect. The four points in the QC group were in clustered, which proved that the machine was stable.

FIGURE 5

FIGURE 5. UPLC-MS/MS TIC diagram: (A) Control group, (B) Model group (C) 5-HMF-H group, (D) 5-HMF-M group, (E) 5-HMF- L group.

FIGURE 6

FIGURE 6. Multivariate data of UPLC-MS/MS: (A) Principal Component Analysis Score Chart, (B) OPLS-DA Analysis Score Chart (C) OPLS-DA S-plot Chart, (D) OPLS-DA V-plot Chart (E) Partial Least Squares-Data Analysis Model Validation Diagram, (F) PLS-DA Analysis Score Diagram. Plasma samples were collected from the different groups.

3.5.2 Identification of Endogenous Metabolites

The difference metabolites between the control group and the model group were screened according to the VIP value (VIP >1.0) and t-test (p < 0.05). As a result, a total of 65 different metabolites were screened (Table 1). Compared with the control group, 35 different metabolites (epitestosterone, etc.) were significantly increased, and 30 different metabolites (17a-hydroxypregnenolone, etc.) significantly reduced. The 5-HMF-H group could regulate 48 metabolites (p-xanthine, etc.). The 5-HMF-M group could regulate 50 kinds of metabolites (5-KETE, etc.). The 5-HMF-L group could regulate 57 metabolites (5-KETE, etc.). The 5-HMF-L group had significantly more metabolic differences than the other two groups.

TABLE 1

TABLE 1. Detection of BDS-related Metabolites in Plasma by LC-MS. Trend 1 is Control group compared with Model group; Trend 2 is Model group compared with 5-HML-L group.

There were 9 metabolites involved in metabolism in the control group, model group, and 5-HMF-L group (Figure 7). To further clarify the distribution of 9 different metabolites in different groups, unsupervised clustering was performed by using a hierarchical clustering analysis heat map. As shown in Figure 8, there were significant differences between the model group and the control group. These metabolites could be significantly regulated, and metabolic disorders in rats after 5-HMF treatment were improved.

FIGURE 7

FIGURE 7. Comparison of the relative content of different metabolites in plasma of the different administration groups. All data are expressed as mean ± standard deviation (n = 8). ^#p < 0.05, ^##p < 0.01, ^###p < 0.001 compared with the control group; *p < 0.05, ***p < 0.001 compared with the model group.

FIGURE 8

FIGURE 8. Cluster analysis heatmap of different metabolites in BDS rat plasma. Brown and blue represent higher and lower than average levels of different metabolites, respectively. The rows represent different metabolites, and the columns represent rat plasma samples (n = 8).

3.5.3 Analysis of Metabolic Pathways

The metabolites of the 5-HMF-L group were input into MetaboAnalyst to explore the metabolic pathways of 5-HMF in the treatment of BDS. The results showed that the 5-HMF-L group may play a role through metabolism of cysteine, methionine, glycerophospholipid, sphingolipid, caffeine, vitamin B6, glycerol, and tyrosine and steroid hormone biosynthesis (Figure 9).

FIGURE 9

FIGURE 9. Pathway analysis of plasma samples from BDS rats.

3.5.4 Pathway Analysis of Targets and Metabolites

To explore the crucial metabolic pathways, joint pathway analysis by MetaboAnalyst was conducted with the 57 differential metabolites and 84 targets. The results showed that only steroid hormone biosynthesis was simultaneously enriched with the targets from network pharmacology and the differential metabolites from metabolomics. According to the number of targets from the network pharmacology in the pathway, steroid hormones biosynthesis was selected as the most crucial metabolic pathway (Figure 10).

FIGURE 10

FIGURE 10. Pathway analysis of plasma samples in 5-HMF-L group.

3.6 Effects of 5-Hydromethylfurfural on the Expression of Key Proteins in the Spleen of Blood Deficiency Rats

Compared with the control group, the expression of CYP17A1 and HSD3B1 protein in the spleen of the model group was decreased (p < 0.05). Compared with the model group, the 5-HMF-H group had increased expression of CYP17A1 protein (p < 0.05), and the 5-HMF-M and 5-HMF-L groups had significantly increased expression of CYP17A1 protein (p < 0.01). Compared with the model group, the 5-HMF-L group had increased expression of HSD3B1 protein (p < 0.05) (Figure 11).

FIGURE 11

FIGURE 11. Effect of each group on protein expression. (A) Control group, (B) Model group, (C) 5-HMF-H group, (D) 5-HMF-M group, (E) 5-HMF-L group. ^#p < 0.05, compared with the control group; *p < 0.05 and **p < 0.01, compared with the model group.

4 Discussion

According to the theory of TCM, the pathogenesis of BDS might be qi-blood discord, the imbalance between yin and yang, blood stasis, and visceral dysfunction. Therefore, Chinese herbal medicines that could nourish yin and blood, fill the essence, and nourish the marrow were often used clinically to treat BDS and anemia (Bailly et al., 1978). Another medicine, RP, has been used for more than 1,000 years to nourish yin and blood, but the possible mechanism of the effect of “tonifying blood” is still unclear. The appearance of rats in the model group changed significantly, including hair loss, slow movement, and weight loss. From the behavior of BDS rats, 5-HMF was effective in treating BDS. Peripheral RBC, WBC, HGB, and HCT were significantly reduced, a typical clinical manifestation of BDS and anemia (Safeukui et al., 2015). At the same time, the 5-HMF-H group did not have increases in WBC, RBC, and HGB, but the 5-HMF-M group and the 5-HMF-L group had more obvious therapeutic effects. Another indication of anemia is splenomegaly, which is one of the most common and earliest pathological manifestations in different types of anemia (Altamura et al., 2019). It was generally believed that hematopoietic dysfunction was related to the discharge of bone marrow cells from bone marrow through blood, resulting in accumulation in the spleen (Koldkjaer et al., 2013). The collection of RBCs in the enlarged spleen could exacerbate anemia. Therefore, the reduction of RBCs in circulation and anemia are related to the degree of spleen enlargement (Gilep et al., 2011). According to Figure 2, the spleen mass in the model group was significantly greater than that in the control group. At the same time, compared with the model group, the spleen of the rats treated with 5-HMF was not enlarged. This showed that 5-HMF could reduce the side effects of chemotherapy on the spleen.

First, using the network pharmacology method, the therapeutic effect of 5-HMF on BDS was analyzed from the two perspectives of target and pathway. From the results of network pharmacology, the regulation of metabolic pathways might be one of the main mechanisms in 5-HMF treatment of BDS. Second, from the perspective of metabolomics, 5-HMF did have a curative effect on BDS, and the curative effect in the 5-HMF-L group was greater than that in the other two groups. Therefore, it showed the efficacy of 5-HMF in BDS treatment, and the most effective dose was in the 5-HMF-L group. Combined with network pharmacology and metabolomics analysis, the steroid hormone biosynthesis may be the key pathway of 5-HMF in the treatment of BDS.

Steroids are commonly used drugs to treat autoimmune hemolytic anemia, but long-term use of steroids can lead to drug resistance (Schiffer et al., 2019). Steroid hormones play a vital role in regulating water and salt balance, metabolism and stress response, and initiating and maintaining sexual differentiation and reproduction (Porcu et al., 2016). The production of new steroids produces steroid hormones in the adrenal cortex, gonads, and placenta. In addition, a series of neurosteroids are produced in the brain (Christakoudi et al., 2018; Sackett et al., 2018). They can promote the body’s hematopoietic function, increase HGB, accelerate the circulation of oxygen and CO₂, and enhance physical strength (Wu et al., 2016).

The metabolic pathway of steroid hormone biosynthesis was entered into the KEGG database and compared with the critical metabolic pathway genes predicted by network pharmacology. We found that CYP17A1 and HSD3B1 are essential proteins for steroid synthesis (Almassi et al., 2018). Western blot verified that in the 5-HMF-L group, the expression of CYP17A1 and HSD3B1 proteins was significantly up-regulated, which was consistent with the changing trends of metabolic results, verifying v the reliability of the metabolomics and network pharmacology prediction. Therefore, we have reason to believe that the steroid hormone biosynthesis pathway may be related to BDS and may be the potential mechanism of 5-HMF in the treatment of BDS. Our results suggest that 5-HMF provides a new therapeutic target by regulating related signaling pathways.

5 Conclusion

Our data showed that 5-HMF could ameliorate blood deficiency in rats induced by chemotherapy. Its mechanism might be regulating the expression of CYP17A1 and HSD3B1 in the steroid hormone biosynthesis pathway. These findings suggest that 5-HMF may be an effective alternative drug for BDS treatment.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics Statement

The animal study was reviewed and approved by Ethics Committee of Heilongjiang University of Traditional Chinese Medicine.

Author Contributions

WZ and NC have contributed equally to this work. HK and QW are responsible for this article.

Funding

Chief Scientist of Qi-Huang Project of National Traditional Chinese Medicine Inheritance and Innovation “One Hundred Million” Talent Project (2021). Qi-Huang Scholar of National Traditional Chinese Medicine Leading Talents Support Program (2018). Heilongjiang Touyan Innovation Team Program.

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.

Acknowledgments

The authors thank Charlesworth Author Services (www.dashboard.cwauthors.com.cn) for English language editing and review services.

Abbreviations

SD, Standard deviation; EDTA, Ethylenediaminetetraacetic acid; TCMSP, Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform; OMIM, Online Mendelian Inheritance in Man; KEGG, Kyoto Encyclopedia of Genes and Genomes; GO, Gene Ontology; VIP, Variable importance in the projection.

References

Almassi, N., Reichard, C., Li, J., Russell, C., Perry, J., Ryan, C. J., et al. (2018). HSD3B1 and Response to a Nonsteroidal CYP17A1 Inhibitor in Castration-Resistant Prostate Cancer. JAMA Oncol. 4 (4), 554–557. doi:10.1001/jamaoncol.2017.3159

PubMed Abstract | CrossRef Full Text | Google Scholar

Altamura, S., Vegi, N. M., Hoppe, P. S., Schroeder, T., Aichler, M., Walch, A., et al. (2019). Glutathione Peroxidase 4 and Vitamin E Control Reticulocyte Maturation, Stress Erythropoiesis and Iron Homeostasis. Haematologica 105 (4), 937–950. doi:10.3324/haematol.2018.212977

PubMed Abstract | CrossRef Full Text | Google Scholar

Amberger, J. S., Bocchini, C. A., Schiettecatte, F., Scott, A. F., and Hamosh, A. (2015). OMIM.org: Online Mendelian Inheritance in Man (OMIM®), an Online Catalog of Human Genes and Genetic Disorders. Nucleic Acids Res. 43 (Database issue), D789–D798. doi:10.1093/nar/gku1205

PubMed Abstract | CrossRef Full Text | Google Scholar

Bailly, Y., Sergant, M., and Raynaud, G. (1978). Hepato- and Thyreotoxicity of Tifemoxone. Arch. Toxicol. Suppl. 1, 229–233. doi:10.1007/978-3-642-66896-8_42

PubMed Abstract | CrossRef Full Text | Google Scholar

Brown, G. (2020). Towards a New Understanding of Decision-Making by Hematopoietic Stem Cells. Int. J. Mol. Sci. 21 (7), 2362. doi:10.3390/ijms21072362

PubMed Abstract | CrossRef Full Text | Google Scholar

Christakoudi, S., Runglall, M., Mobillo, P., Rebollo-Mesa, I., Tsui, T. L., Nova-Lamperti, E., et al. (2018). Steroid Regulation: An Overlooked Aspect of Tolerance and Chronic Rejection in Kidney Transplantation. Mol. Cel Endocrinol 473, 205–216. doi:10.1016/j.mce.2018.01.021

CrossRef Full Text | Google Scholar

Duan, S. (2021). Effects of Angelica on Na∼+-K∼+-ATPase and Ca∼(2+)-Mg∼(2+)-ATPase in Blood Deficiency Rats[J]. Gansu Anim. Husbandry Vet. Med. 51 (04), 50–53.

Google Scholar

Gfeller, D., Grosdidier, A., Wirth, M., Daina, A., Michielin, O., and Zoete, V. (2014). SwissTargetPrediction: a Web Server for Target Prediction of Bioactive Small Molecules. Nucleic Acids Res. 42 (Web Server issue), W32–W38. doi:10.1093/nar/gku293

PubMed Abstract | CrossRef Full Text | Google Scholar

Gilep, A. A., Sushko, T. A., and Usanov, S. A. (2011). At the Crossroads of Steroid Hormone Biosynthesis: the Role, Substrate Specificity and Evolutionary Development of CYP17. Biochim. Biophys. Acta 1814 (1), 200–209. doi:10.1016/j.bbapap.2010.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, W., Zhu, S., Chen, C., Yin, Q., Li, X., Du, G., et al. (2019). The Anti-depression Effect of Angelicae Sinensis Radix Is Related to the Pharmacological Activity of Modulating the Hematological Anomalies. Front. Pharmacol. 10, 192. doi:10.3389/fphar.2019.00192

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, Z. J., Zhai, H. Q., Yuan, Y. P., Zhang, T., Li, R., Kong, L. T., et al. (2016). National Physician Master Jin Shiyuan's Dispensing Technology of Rehmanniae Radix Praeparata Based on Li Shizhen's Pharmaceutical Academic Thought. Zhongguo Zhong Yao Za Zhi 41 (9), 1759–1762. Chinese. doi:10.4268/cjcmm20160930

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C., Lin, S. Y., and Zhao, J. L. (2013). Zhongguo Zhong Xi Yi Jie He Za Zhi Zhongguo Zhongxiyi Jiehe Zazhi = Chinese Journal of Integrated Traditional and Western Medicine. Chin. J. integrated traditional West. Med. 33 (1), 44–46.

Google Scholar

Kitagawa, I., Fukuda, Y., Taniyama, T., and Yoshikawa, M. (1995). Chemical Studies on Crude Drug Processing. X. On the Constituents of Rehmanniae Radix (4) : Comparison of the Constituents of Various Rehmanniae Radixes Originating in China, Korea, and Japan. Yakugaku Zasshi 115 (12), 992–1003. doi:10.1248/yakushi1947.115.12_992

PubMed Abstract | CrossRef Full Text | Google Scholar

Koldkjaer, P., McDonald, M. D., Prior, I., and Berenbrink, M. (2013). Pronounced In Vivo Hemoglobin Polymerization in Red Blood Cells of Gulf Toadfish: a General Role for Hemoglobin Aggregation in Vertebrate Hemoparasite Defense? Am. J. Physiol. Regul. Integr. Comp. Physiol. 305 (10), R1190–R1199. doi:10.1152/ajpregu.00246.2013

CrossRef Full Text | Google Scholar

Li, J., Zhang, L. P., Liu, W., Zhang, Z. L., and Liu, L. J. (2005). Studies on the Changes of 5-hydroxymethyl Furfuraldehyde Content in Radix Rehmanniae Steamed for Different Time. Zhongguo Zhong Yao Za Zhi 30 (18), 1438–1440. Chinese. PMID: 16381466.

PubMed Abstract | Google Scholar

Li, M., Jiang, H., Hao, Y., Du, K., Du, H., Ma, C., et al. (2021). A Systematic Review on Botany, Processing, Application, Phytochemistry and Pharmacological Action of Radix Rehmnniae. J. ethnopharmacology 2021, 114820. Advance online publication. doi:10.1016/j.jep.2021.114820

CrossRef Full Text | Google Scholar

Li, M. M., Wu, L. Y., Zhao, T., Xiong, L., Huang, X., Liu, Z. H., et al. (2011). The Protective Role of 5-HMF against Hypoxic Injury. Cell Stress Chaperones 16 (3), 267–273. doi:10.1007/s12192-010-0238-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, P. L., Sun, H. G., Hua, Y. L., Ji, P., Zhang, L., Li, J. X., et al. (2015). Metabolomics Study of Hematopoietic Function of Angelica Sinensis on Blood Deficiency Mice Model. J. Ethnopharmacol 166, 261–269. doi:10.1016/j.jep.2015.03.010

CrossRef Full Text | Google Scholar

Li, X., Dao, M., Lykotrafitis, G., and Karniadakis, G. E. (2017). Biomechanics and Biorheology of Red Blood Cells in Sickle Cell Anemia. J. Biomech. 50, 34–41. doi:10.1016/j.jbiomech.2016.11.022

CrossRef Full Text | Google Scholar

Liu, W., Fan, Y., Tian, C., Jin, Y., Du, S., Zeng, P., et al. (2020). Deciphering the Molecular Targets and Mechanisms of HGWD in the Treatment of Rheumatoid Arthritis via Network Pharmacology and Molecular Docking. Evid. Based Complement. Alternat Med. 2020, 7151634. doi:10.1155/2020/7151634

PubMed Abstract | CrossRef Full Text | Google Scholar

Lucas, A., Ao-Ieong, E. S. Y., Williams, A. T., Jani, V. P., Muller, C. R., Yalcin, O., et al. (2019). Increased Hemoglobin Oxygen Affinity with 5-Hydroxymethylfurfural Supports Cardiac Function during Severe Hypoxia. Front. Physiol. 10, 1350. doi:10.3389/fphys.2019.01350

PubMed Abstract | CrossRef Full Text | Google Scholar

Milano, M., and Schneider, M. (2007). EPO in Cancer Anemia: Benefits and Potential Risks. Crit. Rev. Oncol. Hematol. 62 (2), 119–125. doi:10.1016/j.critrevonc.2006.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, L., Li, Z., Wang, Y., Zhang, B., Liu, G., and Liu, J. (2020). Network Pharmacology and Metabolomics Study on the Intervention of Traditional Chinese Medicine Huanglian Decoction in Rats with Type 2 Diabetes Mellitus. J. Ethnopharmacol 258, 112842. doi:10.1016/j.jep.2020.112842

CrossRef Full Text | Google Scholar

Porcu, P., Barron, A. M., Frye, C. A., Walf, A. A., Yang, S. Y., He, X. Y., et al. (2016). Neurosteroidogenesis Today: Novel Targets for Neuroactive Steroid Synthesis and Action and Their Relevance for Translational Research. J. Neuroendocrinol 28 (2), 12351. doi:10.1111/jne.12351

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiang, Y., Liu, J., Dao, M., and Du, E. (2021). In Vitro assay for Single-Cell Characterization of Impaired Deformability in Red Blood Cells under Recurrent Episodes of Hypoxia. Lab. Chip 21 (18), 3458–3470. doi:10.1039/d1lc00598g

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (1), 13. doi:10.1186/1758-2946-6-13

CrossRef Full Text | Google Scholar

Sackett, K., Cohn, C. S., Fahey-Ahrndt, K., Smith, A. R., and Johnson, A. D. (2018). Successful Treatment of Pure Red Cell Aplasia Because of ABO Major Mismatched Stem Cell Transplant. J. Clin. Apher. 33 (1), 108–112. doi:10.1002/jca.21553

CrossRef Full Text | Google Scholar

Safeukui, I., Gomez, N. D., Adelani, A. A., Burte, F., Afolabi, N. K., Akondy, R., et al. (2015). Malaria Induces Anemia through CD8+ T Cell-dependent Parasite Clearance and Erythrocyte Removal in the Spleen. mBio 6 (1), e02493–14. doi:10.1128/mBio.02493-14

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiffer, L., Barnard, L., Baranowski, E. S., Gilligan, L. C., Taylor, A. E., Arlt, W., et al. (2019). Human Steroid Biosynthesis, Metabolism and Excretion Are Differentially Reflected by Serum and Urine Steroid Metabolomes: A Comprehensive Review. J. Steroid Biochem. Mol. Biol. 194, 105439. doi:10.1016/j.jsbmb.2019.105439

CrossRef Full Text | Google Scholar

Shi, X., Tang, Y., Zhu, H., Li, W., Li, Z., Li, W., et al. (2014). Comparative Tissue Distribution Profiles of Five Major Bio-Active Components in normal and Blood Deficiency 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

CrossRef Full Text | Google Scholar

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–33. doi:10.1002/cpbi.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Shen, Y., Wang, S., Li, S., Zhang, W., Liu, X., et al. (2017). PharmMapper 2017 Update: a Web Server for Potential Drug Target Identification with a Comprehensive Target Pharmacophore Database. Nucleic Acids Res. 45, W356–W360. doi:10.1093/nar/gkx374

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Li, C., Chuo, W., Liu, Z., Ouyang, Y., Li, D., et al. (2013). Integrated Proteomic and Metabolomic Analysis Reveals the NADH-Mediated TCA Cycle and Energy Metabolism Disorders Based on a New Model of Chronic Progressive Heart Failure. Mol. Biosyst. 9 (12), 3135–3145. doi:10.1039/c3mb70263d

PubMed Abstract | CrossRef Full Text | Google Scholar

Wishart, D. S., Feunang, Y. D., Guo, A. C., Lo, E. J., Marcu, A., Grant, J. R., et al. (2018). DrugBank 5.0: a Major Update to the DrugBank Database for 2018. Nucleic Acids Res. 46 (D1), D1074–D1082. doi:10.1093/nar/gkx1037

PubMed Abstract | CrossRef Full Text | Google Scholar

Won, T. H., Liao, L., Kang, S. S., and Shin, J. (2014). Simultaneous Analysis of Furfural Metabolites from Rehmanniae Radix Preparata by HPLC-DAD-ESI-MS. Food Chem. 142, 107–113. doi:10.1016/j.foodchem.2013.07.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y. S., Chang, C. H., and Nan, F. H. (2016). Steroid Hormone "cortisone" and "20-hydroxyecdysone" Involved in the Non-specific Immune Responses of white Shrimp (Litopenaeus Vannamei). Fish. Shellfish Immunol. 56, 272–277. doi:10.1016/j.fsi.2016.07.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Xia, F., Liu, C., and Wan, J. B. (2020). Characterization of the Cold and Hot Natures of Raw and Processed Rehmanniae Radix by Integrated Metabolomics and Network Pharmacology. Phytomedicine 74, 153071. doi:10.1016/j.phymed.2019.153071

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, J., Wu, J., Zhu, L. Y., Shen, H., Xu, J. D., Jensen, S. R., et al. (2012). Simultaneous Determination of Iridoid Glycosides, Phenethylalcohol Glycosides and Furfural Derivatives in Rehmanniae Radix by High Performance Liquid Chromatography Coupled with Triple-Quadrupole Mass Spectrometry. Food Chem. 135 (4), 2277–2286. doi:10.1016/j.foodchem.2012.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

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. Computer-Aided Mol. Des. 30, 413. doi:10.1007/s10822-016-9915-2

CrossRef Full Text | Google Scholar

Yong, W., Chun, L., Wenjing, C., Zhongyang, L., Yulin, O., Dong, L., et al. (2012). Integrated Proteomic and Metabolomic Analysis Reveals Nadh-Mediated Tca Energy Metabolism Disorder in Chronic Progressive Heart Failure. Heart 98 (Suppl. 2), E115. doi:10.1136/heartjnl-2012-302920b.12

CrossRef Full Text | Google Scholar

Zee-Cheng, R. K. (1992). Shi-quan-da-bu-tang (Ten Significant Tonic Decoction), SQT. A Potent Chinese Biological Response Modifier in Cancer Immunotherapy, Potentiation and Detoxification of Anticancer Drugs. Methods Find Exp. Clin. Pharmacol. 14 (9), 725–736.

PubMed Abstract | Google Scholar

Zhang, H., Wang, H. F., Liu, Y., Huang, L. J., Wang, Z. F., and Li, Y. (2014). The Haematopoietic Effect of Panax Japonicus on Blood Deficiency Model Mice. J. Ethnopharmacol 154 (3), 818–824. doi:10.1016/j.jep.2014.05.008

CrossRef Full Text | Google Scholar

Zhang, Z., Zhang, Y., Gao, M., Cui, X., Yang, Y., van Duijn, B., et al. (2020). Steamed Panax Notoginseng Attenuates Anemia in Mice with Blood Deficiency Syndrome via Regulating Hematopoietic Factors and JAK-STAT Pathway. Front. Pharmacol. 10, 1578. doi:10.3389/fphar.2019.01578

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat. Commun. 10 (1), 1523. doi:10.1038/s41467-019-09234-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, M. F., Liu, X. Q., Ju-hee, O., Chang-soo, Y., and Je-Hyun, L. (2007). Zhongguo Zhong Yao Za Zhi = Zhongguo Zhongyao Zazhi = China Journal of Chinese Materia Medica. China J. Chin. materia Med. 32 (12), 1155–1157.

Google Scholar