The Medicinal Plant Pair Bupleurum chinense-Scutellaria baicalensis – Metabolomics and Metallomics Analysis in a Model for Alcoholic Liver Injury

Traditional Chinese Medicine (TCM), a complex natural herbal medicine system, has increasingly attracted attention from all over the world. Most research has illustrated the mechanism of TCM based on the active components or single herbs. It was fruitful and effective but far from satisfactory as it failed to gain insights into the interactivity and combined effects of TCM. In this work, we used Bupleurum chinense (B. chinense DC, a species in the genus Bupleurum, family Apiaceae) and Scutellaria baicalensis (S. baicalensis Georgi, a species in the genus Scutellaria, family Lamiaceae), an herbal pair in TCM, to illustrate the combined effect. We compared the diverse effects between the B. chinense-S. baicalensis herbal pair and its compositions in an animal model of Alcoholic Liver Injury to highlight the advantages of the formula. Biochemical and histological indicators revealed that the effect of B. chinense-S. baicalensis was better than its individual parts. Furthermore, metabolite profiling of the serum, liver tissue, and feces were conducted to reveal that the herbal pair largely presented its effects through enhanced tissue penetration to maintain liver-located intervention with less global and symbiotic disturbance. Furthermore, we analyzed the distribution of the metal elements in extracts of the serum and liver tissue and found that the herbal pair significantly regulated the distribution of endogenous selenium in liver tissue. As selenium plays an important role in the anti-oxidative and hepatoprotective effects, it may be the reason for combined effects in BS formula. This research could open new perspectives for exploring the material basis of combined effects in natural herbal medicine.

Introduction

Traditional Chinese Medicine, a complex medical system with a long history of usage, has increasingly been examined. Intervention in the human body by compound herbal formulas is one of the major characters of TCM. So far, many endeavors have been made to unveil the mechanisms of TCM by way of the separation strategy, which is to study the pharmacodynamic material basis by tapering the entire formula to single herbs, even to single active components. However, study of the complexity of TCM requires a different strategy as the separation strategy may lose sight of the interactivity and combined effects (Zheng et al., 2013). In fact, combination therapies have been increasingly utilized during the past few years (Yohannes et al., 2010). A multi-drug strategy usually yields more therapeutic efficacy with less adverse effects, which may lead to a triumph for modern medicine and provide fertile ground for modern drug development (Lee et al., 2011; Wang et al., 2013). In particular, the herbal pair, the minor formula composed of two herbs with distinct properties, is used by clinical TCM practitioners to increase the therapeutic effect or reduce toxic effects (Li et al., 2010). To investigate the integration mechanism of the formula, technologies with holistic concepts are needed (Cappuccio et al., 2015). Recently, systems biology combined with physiology, biochemistry, molecular biology and tongue image digitization were employed to evaluate the treatment effect and explore the mechanisms of the TCM formula (Wu et al., 2010; Dai et al., 2013a; Sun et al., 2013). The focus has been on the constituent changes and chemical transformations to explain the potential combined effect of the TCM formula (Li et al., 2010; Ye et al., 2011; Cao et al., 2014; Zhang Y. et al., 2016). Furthermore, metal or semi-metal elements have attracted our attention for the numerous physiological and biochemical functions, thus helping researchers to better illustrate the complex pharmacodynamic mechanism of TCM (Dong et al., 2006; Yan et al., 2015; Deng et al., 2018).

Bupleurum chinense-Scutellaria baicalensis (BS), the simultaneous usage of B. chinense and S. baicalensis, was originally recorded in the ancient Chinese medical book Shanghan Lun (Treatise on Cold Damage Diseases). It is one of the most common herbal pairs used in clinical trials to treat hepatic disease (Lin et al., 2011; Wang et al., 2014). It is considered as being able to achieve better effects with less toxicity, yet the in-depth mechanism is unclear. Physiological and biochemistry studies were conducted to highlight the advantages of the formula in an ALI animal model. Then, systematic metabolomics combined with metallomics were further employed to illustrate the precise regulation of disease-intervened metabolites and trace element distribution in liver tissue, which may explain the combined effects present in the BS formula. Our research not only presented the systematic comparison of the pharmacodynamic effects between herbal pairs and its compositions, but also provided new perspective for exploring the material basis of the formula’s combined effect.

Materials and Methods

Materials and Reagents

Analytical grade methoxyamine hydrochloride, methanol, ethanol, chloroform, and pyridine were obtained from China National Pharmaceutical Group Corporation in Shanghai of China. N,O-Bis(trimethylsilyl)trifluoroacetamide (BSTFA:TMCS, 99 : 1), heptadecanoic acid and urease were purchased from Sigma-Aldrich. Liquor (56 %, v/v, 00215) was obtained from Beijing Red Star, China. Double-distilled water was produced by a Milli-Q Ultra-pure water system (Millipore Corporation, United States). Dry S. baicalensis and B. chinense were purchased from Inner Mongolia and An’hui was identified as S. baicalensis Georgi and B. chinense DC by Dr. Lihong Wu of Shanghai University of Traditional Chinese Medicine. Baicalin and Saikosaponin were purchased from Ciyuan Biotechnology (Shanxi, China).

Animal Experiments and Samples Collection

All protocols for animal experimentation and maintenance were approved by the Animal Ethics Committee of Shanghai University of Traditional Chinese Medicine. Male Wistar rats (weighing 200 ± 10g) were provided with unrestricted amounts of food and water, housed in temperature- and humidity-controlled rooms, and kept on a 12-h light/dark cycle. One week later, all rats were divided into 8 groups of 8 rats each, randomized as follows: normal group, ALI modeling group and other six treatment groups. The ALI modeling groups were fed a high-glucose-fat diet (Hapon et al., 2005) (prescription: 10% cane sugar, 2% cholesterol, 10 % lard, 5 % yolk powder, 0.2% propylthiouracil, and 72.8% basal diet, purchased by Suzhou Shuangshi Laboratory Animal Feed Science, Suzhou, China), combined with 10 ml/kg liquor (Niemela et al., 1995) every 2 days for 30 days. The biochemistry indexes and histological results were employed to illustrate the successful modeling of ALI. All treatment groups were performed by intragastric administration, respectively, with decoctions of B. chinense (BU, amounting to 4 g/kg dried medicinal herbs), decoctions of S. baicalensis (SC, amount to 4 g/kg raw dried medicinal herbs), decoctions of B. chinense-S. baicalensis (BS, 1:1, amounting to 4 g/kg dried medicinal herbs), saikosaponin solution in water (the major active component of B. chinense referred to in Chinese Pharmacopeia (Lu et al., 2012) as SA, amounting to 4 g/kg dried medicinal herbs), baicalin solution in water [the major active component of S. baicalensis referred to in Chinese Pharmacopeia (Xi et al., 2015) as BA, amounting to 4 g/kg dried medicinal herbs], as well as the mixture of BA and SA (1:1, MI, amounting to 4 g/kg dried medicinal herbs) for 20 days. Additionally, normal and modeling rats were given equivalent double-distilled water.

After the modeling and treatment, the rats were placed in individual metabolic cages with fasting for 24 h to collect the urine and feces samples, which were kept at -80°C until analysis. Then, all the animals were anesthetized, and the serum samples were collected via abdominal aorta. All serum samples were centrifuged at 3,000 rpm at 4°C for 10 min. The supernatants were stored at -80°C for metabolomics and biochemical analysis. Liver tissues, dissected from the same lobe of each animal, were fixed in 4% paraformaldehyde overnight at room temperature, dehydrated through ethanol gradients, infiltrated in xylene, and embedded in paraffin. Serial sections of 5 μm thickness were made from paraffin-embedded tissue for hematoxylin and eosin (H&E) staining. Another set of liver samples, taken from the same lobe of each animal, were snap-frozen in liquid N₂, and stored at -80°C.

Biochemistry Indexes Analysis

Alanine aminotransferase, AST, GLU, CHOL, and TG were assayed by a Hitachi 7080 automatic biochemistry analyzer. All results were presented as the mean ± SD and analyzed using SPSS 17.0 (SPSS, Chicago, IL, United States). Multiple comparisons for two groups were performed by independent sample T-test. P-values below 0.05 were considered as statistically significant.

Metabolomic Analysis

All serum, feces, and tissues were prepared for metabolomics analysis. For serum samples, the preparation was the same as in our previous reports (Sun et al., 2012; Dai et al., 2013b,c). For liver tissue, 50 mg of tissue was homogenized by 150 μL of normal saline. And for feces, 80 mg of fecal sample was transferred into a screw tube with 320 μL of normal saline and extracted with ultrasonic treatment for 20 min. Then, the derivatization process was performed as in our previous publications (Sun et al., 2012; Dai et al., 2013b,c).

All GC-MS analyses were performed by a mass spectrometer (5975B, Agilent, United States) coupled with a gas chromatography instrument (6890, Agilent, United States). The GC-MS operating condition was adapted from our previous reports (Sun et al., 2012; Dai et al., 2013b,c). The programmed column temperature was optimized to acquire a well separation; the column temperature was initially maintained at 70°Ñ for 2 min and then increased from 70 to 120°Ñ at a rate of 5°Ñ/min for 2 min. Then, the column temperature was increased to 190°Ñ at a rate of 3°Ñ/min, to 210°Ñ at a rate of 5°Ñ/min, to 260°Ñ at a rate of 10°Ñ/min. After that, the column temperature was increased to 290°Ñ at a rate of 5°Ñ/min and held for 3 min. The method of data analysis was also adapted from our previous reports (Sun et al., 2012; Dai et al., 2013b,c), optimizing the VIP values to 1.0 and the threshold of P value to 0.1. The variables that were discovered were identified by searching the NIST 2011 database and verified by standards. The Kyoto Encyclopedia of Genes and Genomes (KEGG)¹ and the Metabolites Biological Role (MBRole)² were employed to highlight the most related pathway.

Metallomic Analysis

Five hundred (500) μL of herbal extracts, 150 μL of serum or 100 mg of liver sample were removed into digestion tubes which were soaked with in 30% HNO₃ at least 24 h followed by another 24 h in high purity water, and then rinsed again with high purity water. Three mL of 65% HNO₃ was added into the tubes and these tubes were then put into an airing chamber for 12 h. The samples carried out microwave digestion with the different temperature program (Supplementary Data Sheet 1 for details). Subsequently, these solutions were transferred into 50 mL centrifuge tubes and the tubes were washed with high ultrapure water three times to ensure that the concentration of HNO₃ was less than 5%. All these solutions were stored at 4°C (Mesko et al., 2011), and awaited analysis by ICP-MS.

Environmental calibration standard solutions—10 mg/L^-1 of Al, As, Ba, Be, Co, Cr, Cu, Mn, Ni, Pb, Se, V, and Zn, in addition to 1000 mg/L^-1 of Ca, Fe, K, Mg, and Na—were diluted into a series of concentrations with 5% HNO₃.

Before ICP-MS analysis, it was necessary to be purged by helium (99.99%) at 2 L/min. When the forward power reached 1500 W and the reflected power reached 1 W, it indicated that the system could ignite to start up. The experimental parameters used for ICP-MS analysis are summarized in Supplementary Data Sheet 1. Furthermore, calibration curve, LLOQ (the lower limit of quantification), precision test and repeated experiment were carried out to validate the method reliability.

All results were presented as the mean ± SD and analyzed using SPSS 17.0 (SPSS, IBM, United States). Multiple comparisons for two groups were performed by independent sample T-test. A P value less than 0.05 was considered statistically significant.

Results

Pharmacodynamic Evaluation of Bupleurum-Scutellaria and Its Compositions

First, biochemistry indexes and histological results were employed to illustrate pharmacodynamic effects of B. chinense-S. baicalensis and its compositions. The biomarkers, including ALT, AST, GLU, CHOL, and TG, which reflect the liver-associated disorder, were used to evaluate the physiological changes with treatment of BS and its compositions in ALI (detail in Supplementary Table S1). We found that decoctions of B. chinense-S. baicalensis and S. baicalensis, alongside a water solution of Baicalin and in addition to the mixture of Baicalin and Saikosaponin, could reverse the biochemical disturbances in ALI (Figures 1A,B). However, B. chinense and its major active component, Saikosaponin, present little pharmacodynamic effects (Figures 1A,B). Then, the histological results evaluated the improvement of pathological conditions. Some liver injury characteristics, such as liver cell swelling, lipid droplets accumulation, and inflammatory infiltration occurred in the ALI model group (Figure 1C). It is worth pointing out that B. chinense-S. baicalensis nearly reversed all ALI-induced pathological conditions to normal (Figure 1C). The decoctions of S. baicalensis showed slightly worse reversion than B. chinense-S. baicalensis in all ALI-induced pathological changes. Treatment of Baicalin and the mixture of Baicalin and Saikosaponin could significantly reduce lipid droplets accumulation (Figure 1C). Furthermore, we found that B. chinense and Saikosaponin did not present improvement of liver injury characteristics (Figure 1C), which was consistent with the results of biochemistry indexes. So, it seems that SC and BA may, respectively, play an important role in BS and MI’s treatment efficacy for ALI. More importantly, these results provided evidence of the synergy, more exact curative effect, in BS treatment in the ALI group at both the biochemical and histopathological levels. Why it happened? How were these emergent characteristics brought out? The above questions prompted us to investigate the potential mechanism by which the efficacy of B. chinense-S. baicalensis, rather than the single herbs and their active ingredients, was enhanced.

FIGURE 1

Figure 1. Biochemical and histological analysis. (A) Biochemical analysis of major biomarkers for liver functions in blood of rates. The biochemical indicators of Bupleurum chinense-Scutellaria baicalensis and their composition were compared to the normal control group; (B) liver tissue histopathology (H&E staining) of every group; (C) the quantized results of liver tissue histopathology. The results were based on the area of intercellular space and lipid droplets to present liver cell swelling and lipid droplets accumulation. N: Normal control; M: Alcoholic Liver Injury Modeling; ALI: Alcoholic Liver Injury; BS: decoctions of B. chinense-S. baicalensis; BU: decoctions of B. chinense; SC: decoctions of S. baicalensis; SA: water solution of Saikosaponin; BA: water solution of Baicalin; MI: the mixture water solution of Baicalin and Saikosaponin; Alanine aminotransferase (ALT), Aspartate aminotransferase (AST), Glucose (GLU), Cholesterol (CHOL), Triglyceride (TG).

Systematic Analysis of Drug-Induced Metabolic Perturbations

Compared with the ALI group, the metabolites of serum, liver tissue, and fecal extracts, disturbed by BS or its components, were employed to construct global, local, and symbiotic metabolic profiling (Figure 2A), respectively. First, we analyzed the metabolites that only changed in BS treatment (Table 1 and Supplementary Table S2). After preliminary structural and functional analysis, we found that unique BS-intervened metabolites were enriched in alpha amino acids and derivatives (serine, phenylalanine and ornithine) and aliphatic acyclic compounds (urea and methanamine), which are involved in ATP-binding cassette transporter (urea, phenylalanine, ornithine and serine), Microbial metabolism in diverse environments (methanamine, urea and serine) and biosynthesis of amino acids (phenylalanine, ornithine and serine).

FIGURE 2

Figure 2. Metabolite and pathway analysis of drug-induced perturbations. (A) Numbers of drug-intervened metabolites in ALI; (B) numbers of drug-intervened pathways in ALI. The summarization of drug-intervened global, local and symbiotic pathways is presented in the center. BS: decoctions of B. chinense-S. baicalensis; SC: decoctions of S. baicalensis; BU: decoctions of B. chinense; BA: water solution of baicalin; SA: water solution of saikosaponin; MI: the mixture of baicalin and saikosaponin.

TABLE 1

Table 1. The special metabolites detected in BS group of ALI compared with control group.

However, these metabolites could not provide the mechanisms needed to unveil BS’s unique efficacy. All drug-disturbed metabolites were, respectively, enriched into pathways by MBRole, which were further categorized into functional modules by KEGG (Figure 2B). Surprisingly, we found that the numbers of changed serum metabolites in the BS group were the least among all groups (Figure 2A). According to the results of biochemical indicators and pathology, we presumed that BS’s unique efficiency may be derived from the precise regulation of disease-intervened metabolites with less bypass perturbations. Further metabolic analysis could provide the evidence for this hypothesis. The symbiotic bacteria metabolism was disturbed by BS group most slightly but amino acid metabolism, carbohydrate metabolism and signal transduction. Furthermore, it was observed BS in serum-metabolite profiling only intervened in lipid metabolism, and SC may play a major role in this pathway. However, SC could not intervene in liver metabolism but BU could. We hypothesized that BS largely exerted its prefect effects through tissue penetration with the effect of BU. Integrating the global, local and symbiotic perspective, striking differences were observed in the BS group compared with other groups, such as more effects in amino acid metabolism, nucleotide metabolism and signal transduction but less disturbances in other disease-related metabolism (such as endocrine and metabolic diseases as well as neurodegenerative diseases). Finally, we were eager to find out the reasons for which usage of the herbal pair could provide beneficial pharmacodynamic effects.

Diversity of Metal Elements in Serum and Liver Tissue May Derive From Better Pharmacodynamic Effects of Herbal Pair

Preliminarily, we wanted to explain the mechanism of the formation of new compounds, but previous research proved that there were not more constituents from chemical transformation (Dai Peihui et al., 2017). As even the major active components (Baicalin and Saikosaponin) could not bring better efficacy (Figure 1), we changed our focus to the metal elements with many physiological and biochemical effects related to the pharmacodynamic mechanism (Dong et al., 2006; Yan et al., 2015; Deng et al., 2018) of TCM. Metallomics was employed to explore diverse metal elements in the herbal extraction, serum and liver tissue. This approach could help us to illustrate the combined effects of the B. chinense-S. baicalensis herbal pair.

We found that there was no significant difference at the comparation of metal contents between the average of BU and SC (white column) and BS (black column) (Figure 3A). The different distributions of metal elements in the serum and liver tissue were not because of the distinguished input of mental elements. Interestingly, the distributions of elements were obviously different between the serum and liver tissue. In the serum, most elements presented no differences between the average of BU-treated and SC-treated rats (white column) and BS-treated rats (black column), except arsenic (Figure 3B). However, in liver tissue, calcium, vanadium, manganese, iron, and arsenic were significantly reduced by treatment of BS, while the content of selenium was clearly increased (Figure 3C). Selenium particularly caught our attention because it was not detected in the extraction of BU, SC, and BS (Figure 3A). These results suggested that BS significantly regulated the distribution of endogenous selenium in liver tissue. As selenium plays an important role in the anti-oxidative and hepatoprotective effects, it may be the reason behind the combined effects of the BS formula.

FIGURE 3

Figure 3. B. chinense-S. baicalensis herbal pair altered the distribution of metal contents in the serum and liver tissue. (A) Comparison of metal contents between average of BU and SC (white column) and BS (black column). (B) Comparison of serum metal contents between average of BU-treated and SC-treated rats (white column) and BS-treated rats (black column). (C) Comparison of hepatic metal contents between average of BU-treated and SC-treated rats (white column) and BS-treated rats (black column). BU: decoctions of B. chinense; SC: decoctions of S. baicalensis; BS: decoctions of Bupleurum chinense-S. baicalensis.

Discussion

The TCM formula, which is composed of related herbs under the theories of TCM, has been applied in clinical treatment to provide more efficiency and safety. There is a trend of combination therapy to achieve combined effects in complex diseases, such as breast cancers (Hyung et al., 2011), liver fibrosis (Zhang H. et al., 2016). However, it is challenging to explore the combined effect of more complex TCM with deeper insight. To elucidate the combined mechanisms of the TCM formula, researchers have undertaken much research. The “Realgar-Indigo naturalis formula,” which was composed of Realgar, Indigo naturalis and Salvia miltiorrhiza, has been proven to be very effective in treating human acute promyelocytic leukemia promyelocytic leukemia. Three major active ingredients (tetraarsenic tetrasulfide, indirubin, and tanshinone IIA are the major active ingredients in Realgar, Indigo naturalis and Salvia miltiorrhiza, respectively) yielded synergy in the treatment of acute promyelocytic leukemia in vivo and in vitro. They dissected the mechanisms of TCM formulae at the molecular, cellular and organismal levels (Wang et al., 2008). There have also been other examples, such as Shuang-Huang-Lian (Feng et al., 2005; Di et al., 2006), Yin-Chen-Hao-Tang (Hyung et al., 2011; Wang et al., 2011) and Wen-Xin-Formula (Dai et al., 2013b), yet systematic evaluation for the pharmacodynamic effects and the material basis of the combined effect have not been well illustrated.

In this study, we employed biochemical and pathological analyses combined with metabolomics to provide the system evaluation for the pharmacodynamic effects of the B. chinense-S. baicalensis herbal pair, the single herbs and their active ingredients. In our study, we found that BS clearly intensified the therapeutic efficacy on both biochemical and histological levels. These observations provided insight into the combined effect. As analysis of biochemical and pathological results, we could claim that S. baicalensis and its major active compound, baicalin, may play important roles to ameliorate lipid droplet accumulation, as observed in the B. chinense-S. baicalensis herbal pair group. Some research that showed that baicalin presented antisteatotic effects also support our results (Guo et al., 2009; Kim and Lee, 2012; Xi et al., 2015). Though lipid metabolism disorder is one of the most notable characteristics in ALI, improvement of cell swelling and inflammatory infiltration gradually diminished form the herbal pair to single herb and baicalin (Figure 1). This may suggest that other components in S. baicalensis such as baicalein (Fan et al., 2013; Patwardhan et al., 2016), wogonin (Lin and Shieh, 1996; Lee et al., 2015) and oroxylin A (Shimizu et al., 2018), showed the main anti-inflammatory and hepatoprotective activity. In contrast to our expectation, B. chinense and saikosaponin may not be able to cure ALI (Figure 1). However, the metabolomic analysis of liver indicated that S. baicalensis hardly intervened in the metabolomic pathways in liver tissue (Figure 2B). Yet, B. chinense-S. baicalensis could significantly enhance tissue penetration and influenced the pathways of intestinal flora. This may address one of the common problems in combination therapy, diarrhea. All these results suggested that the B. chinense-S. baicalensis herbal pair presented better effects in the treatment of ALI than single herbs and their active compounds.

To reveal the connotation of the formula, it is necessary to understand not only the main ingredients but also the combined application of multilevel analytical and synthetic research approaches at the molecular, cellular, and organismal levels. In our case, we paid attention to metal elements, as previous research had proved that only the changes of baicalin and saikosaponin (Dai Peihui et al., 2017) could not contribute to the better efficacy in our results (Figure 1). After metallomic analysis of BS herbal pair treatment, we found that arsenic was increased in the serum, but calcium, vanadium, manganese, iron and arsenic were significantly decreased in liver tissue, while contents of selenium were increased in liver tissue (Figures 3B,C). The results revealed that BS herbal pair treatment could reduce toxic metal accumulation in the liver, such as vanadium (Hosseini et al., 2013), manganese (Fujimura et al., 2012), iron (Zhang et al., 2014) and arsenic (Tan et al., 2011; Massey et al., 2015). The arsenic in the serum could intensify the expression of aquaglyceroporin and enhance PML-RAR degradation and anti-inflammatory efficacy (Wang et al., 2008). At the same time, selenium in the liver could significantly inhibit adipocyte hypertrophy and abdominal fat accumulation via induction of fatty acid beta-oxidation (de Castro Barra et al., 2015) and alter lipid metabolism and protein synthesis in the liver (Zhao et al., 2016). Furthermore, selenium could be incorporated into selenoproteins which play important roles in the anti-oxidative stress and hepatoprotective effect via apoptosis inhibition (Tang et al., 2018). The selenoproteins could participate in pyrimidine metabolism (Carlson et al., 2007), which is the unique pathway of intervention by the BS herbal pair in liver tissue (as shown in Table 2). Thus, the transition from pyrimidine to uridine was increased (as shown in Table 2). The enhanced uridine homeostasis-linked liver pyrimidine metabolism could ameliorate lipid accumulation (Le et al., 2013). Furthermore, simultaneously increased contents of selenium and uridine could show a cytoprotection effect through enhancing the availability of precursors and cofactors in the circulation (Rijpma et al., 2015). Finally, selenium could reverse drug-induced immunotoxic effects (Latorre et al., 2011) to relieve drug pressure. In summary, the elevated selenium may be one of the potential causes that lead to enhancing efficacy and reducing toxicity in the B. chinense-S. baicalensis herbal pair.

TABLE 2

Table 2. Contributions of the single herbs and their active ingredients to the disturbed pathways of B. chinense-S. baicalensis herbal pair.

Conclusion

In this study, we conducted an animal model of ALI to study the characteristics of the B. chinense-S. baicalensis herbal pair, which was controlled by B. chinense, S. baicalensis, and their main active ingredients baicalin and saikosaponin. Biochemical and histological indicators revealed that the effect of B. chinense-S. baicalensis is the best in all treatment groups. Furthermore, metabolite profiling of the serum, liver tissue and feces was conducting with GC-MS to comprehensively reveal the global, local and symbiotic metabolic modulation. Our results indicated that BS largely developed its effects through enhanced tissue penetration to maintain liver-located intervention with less global and symbiotic disturbance. Furthermore, we analyzed the distribution of the metal elements in the serum and liver tissue in B. chinense-S. baicalensis, and found that selenium was significantly increased in liver tissue. The increased selenium may show multi-beneficial effects for the treatment efficacy in B. chinense-S. baicalensis. The schematic diagram of the combined effects of the B. chinense-S. baicalensis herbal pair in ALI is shown in Figure 4. This observation preliminarily showed the combined effects in the B. chinense-S. baicalensis herbal pair, which may provide a new perspective for exploring the material basis of combined effects in natural herbal medicine, especially TCM.

FIGURE 4

Figure 4. Schematic diagram of the combined effects of the B. chinense-S. baicalensis herbal pair in Alcoholic Liver Injury. Alcoholic Liver Injury was used to research the characteristics of B. chinense- S. baicalensis, which is controlled by B. chinense and S. baicalensis, and their main active ingredients baicalin and saikosaponin. The biochemical indicators and histology indicated that B. chinense-S. baicalensis presented the best efficacy, which could be further illustrated by global, local, and symbiotic metabolic analysis. Finally, metallomic analysis was employed to uncover the distribution changes of selenium which, unlike other metal elements, was significantly increased by B. chinense-S. baicalensis in liver tissue to demonstrate multiple benefits such as anti-oxidative stress and hepatoprotective effect, to relieve drug pressure.

Author Contributions

YZ and JD contributed to conception and design of the study. NZ, QL, SS, YW, ZD, and JD organized the database. NZ, QL, SS, YW, ZD, and JD performed the statistical analysis. JD wrote the first draft of the manuscript. NZ, QL, SS, YW, ZD, and YZ wrote sections of the manuscript. All authors contributed to manuscript revision, read and approved the submitted version. All authors have given approval to the final version of the manuscript.

Funding

This work is financially supported by the National Natural Science Foundation Project (81373519), Shanghai Interdisciplinary Cultivation Platform of Outstanding and Innovative Postgraduates, and Shanghai “085” Science.

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.

Supplementary Material

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

Abbreviations

ALI, alcoholic liver injury; ALT, alanine aminotransferase; AST, aspartate aminotransferase; BA, baicalin; BS, Bupleurum chinense-Scutellaria baicalensis; BU, Bupleurum chinense; CHOL, cholesterol; GLU, glucose; OPLS, orthogonal partial least squares; PCA, principle component analysis; PLS-DA, partial least squares-discriminant analysis; SA, saikosaponin; SC, Scutellaria baicalensis; TCM, traditional chinese medicine; TG, triglyceride.

Footnotes

1. ^www.genome.ad.jp/kegg/
2. ^http://csbg.cnb.csic.es/mbrole2/

References

Cao, G., Li, Q., Cai, H., Tu, S., and Cai, B. (2014). Investigation of the chemical changes from crude and processed paeoniae radix alba-atractylodis macrocephalae rhizoma herbal pair extracts by using Q exactive high-performance benchtop quadrupole-orbitrap LC-MS/MS. Evid. Based Complement. Alternat. Med. 2014:170959. doi: 10.1155/2014/170959

PubMed Abstract | CrossRef Full Text | Google Scholar

Cappuccio, A., Zollinger, R., Schenk, M., Walczak, A., Servant, N., Barillot, E., et al. (2015). Combinatorial code governing cellular responses to complex stimuli. Nat. Commun. 6:6847. doi: 10.1038/ncomms7847

PubMed Abstract | CrossRef Full Text | Google Scholar

Carlson, B. A., Moustafa, M. E., Sengupta, A., Schweizer, U., Shrimali, R., Rao, M., et al. (2007). Selective restoration of the selenoprotein population in a mouse hepatocyte selenoproteinless background with different mutant selenocysteine tRNAs lacking Um34. J. Biol. Chem. 282, 32591–32602. doi: 10.1074/jbc.M707036200

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, J., Fang, J., Sun, S., Chen, Q., Cao, H., Zheng, N., et al. (2013a). ZHENG-omics application in ZHENG classification and treatment: chinese personalized medicine. Evid. Based Complement. Alternat. Med. 2013:235969. doi: 10.1155/2013/235969

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, J., Sun, S., Cao, J., Zhao, Y., Cao, H., Zheng, N., et al. (2013b). Similar connotation in chronic hepatitis B and nonalcoholic fatty liver patients with dampness-heat syndrome. Evid. Based Complement. Alternat. Med. 2013:793820. doi: 10.1155/2013/793820

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai, J., Sun, S., Peng, J., Cao, H., Zheng, N., Fang, J., et al. (2013c). Exploration of macro-micro biomarkers for dampness-heat syndrome differentiation in different diseases. Evid. Based Complement. Alternat. Med. 2013:706762. doi: 10.1155/2013/706762

PubMed Abstract | CrossRef Full Text | Google Scholar

Dai Peihui, M. R., Gao, J., and Liu, T. (2017). Study on main chemical components from compatibility of chaihu and huangqin by HPLC. World Chin. Med. 12:4.

de Castro Barra, P. M., Sabarense, C. M., Alvarenga, M. B., de Sousa, R. A., and de Oliveira, M. A. (2015). Selenium content in the liver of wistar rats fed diets of different fatty acid quality. Biol. Trace Elem. Res. 168, 441–446. doi: 10.1007/s12011-015-0359-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Deng, Z., Yang, Z., Ma, X., Tian, X., Bi, L., Guo, B., et al. (2018). Urinary metal and metalloid biomarker study of henoch-schonlein purpura nephritis using inductively coupled plasma orthogonal acceleration time-of-flight mass spectrometry. Talanta 178, 728–735. doi: 10.1016/j.talanta.2017.10.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Di, B., Feng, N., and Liu, W. (2006). Pharmacokinetic comparisons of shuang-huang-lian with the different combinations of its constitutional herbs. J. Ethnopharmacol. 107, 401–405. doi: 10.1016/j.jep.2006.03.030

PubMed Abstract | CrossRef Full Text | Google Scholar

Dong, S. F., Liu, J., Dong, H. B., Zhang, W. P., and Liu, J. H. (2006). Study on the determination and blurred classification of metal elements content in danshen etc 11 kinds of Chinese traditional medicine. Guang Pu Xue Yu Guang Pu Fen Xi 26, 2150–2153.

PubMed Abstract | Google Scholar

Fan, G. W., Zhang, Y., Jiang, X., Zhu, Y., Wang, B., Su, L., et al. (2013). Anti-inflammatory activity of baicalein in LPS-stimulated RAW264.7 macrophages via estrogen receptor and NF-kappaB-dependent pathways. Inflammation 36, 1584–1591. doi: 10.1007/s10753-013-9703-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, N. P., Di, B., and Liu, W. Y. (2005). Comparison of the metabolism of baicalin in rats orally administered with radix scutellariae extract and shuang-huang-lian extract. Chem. Pharm. Bull. 53, 978–983. doi: 10.1248/cpb.53.978

PubMed Abstract | CrossRef Full Text | Google Scholar

Fujimura, T., Terachi, T., Funaba, M., and Matsui, T. (2012). Reduction of liver manganese concentration in response to the ingestion of excess zinc: identification using metallomic analyses. Metallomics 4, 847–850. doi: 10.1039/c2mt20100c

PubMed Abstract | CrossRef Full Text | Google Scholar

Guo, H. X., Liu, D. H., Ma, Y., Liu, J. F., Wang, Y., Du, Z. Y., et al. (2009). Long-term baicalin administration ameliorates metabolic disorders and hepatic steatosis in rats given a high-fat diet. Acta Pharmacol. Sin. 30, 1505–1512. doi: 10.1038/aps.2009.150

PubMed Abstract | CrossRef Full Text | Google Scholar

Hapon, M. B., Varas, S. M., Jahn, G. A., and Gimenez, M. S. (2005). Effects of hypothyroidism on mammary and liver lipid metabolism in virgin and late-pregnant rats. J. Lipid Res. 46, 1320–1330. doi: 10.1194/jlr.M400325-JLR200

PubMed Abstract | CrossRef Full Text | Google Scholar

Hosseini, M. J., Shaki, F., Ghazi-Khansari, M., and Pourahmad, J. (2013). Toxicity of vanadium on isolated rat liver mitochondria: a new mechanistic approach. Metallomics 5, 152–166. doi: 10.1039/c2mt20198d

PubMed Abstract | CrossRef Full Text | Google Scholar

Hyung, S. W., Lee, M. Y., Yu, J. H., Shin, B., Jung, H. J., Park, J. M., et al. (2011). A serum protein profile predictive of the resistance to neoadjuvant chemotherapy in advanced breast cancers. Mol. Cell Proteomics 10:M111011023. doi: 10.1074/mcp.M111.011023

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, S. J., and Lee, S. M. (2012). Effect of baicalin on toll-like receptor 4-mediated ischemia/reperfusion inflammatory responses in alcoholic fatty liver condition. Toxicol. Appl. Pharmacol. 258, 43–50. doi: 10.1016/j.taap.2011.10.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Latorre, A. O., Caniceiro, B. D., Wysocki, H. L. Jr., Haraguchi, M., Gardner, D. R., and Górniak, S. L. (2011). Selenium reverses pteridium aquilinum-induced immunotoxic effects. Food Chem. Toxicol. 49, 464–470. doi: 10.1016/j.fct.2010.11.026

PubMed Abstract | CrossRef Full Text | Google Scholar

Le, T. T., Ziemba, A., Urasaki, Y., Hayes, E., Brotman, S., and Pizzorno, G. (2013). Disruption of uridine homeostasis links liver pyrimidine metabolism to lipid accumulation. J. Lipid Res. 54, 1044–1057. doi: 10.1194/jlr.M034249

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, J. Y., Kim, J. Y., Park, G. W., Cheon, M. H., Kwon, K. H., Ahn, Y. H., et al. (2011). Targeted mass spectrometric approach for biomarker discovery and validation with nonglycosylated tryptic peptides from N-linked glycoproteins in human plasma. Mol. Cell Proteomics 10:M111009290. doi: 10.1074/mcp.M111.009290

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, W., Ku, S. K., and Bae, J. S. (2015). Anti-inflammatory effects of baicalin, baicalein, and wogonin in vitro and in vivo. Inflammation 38, 110–125. doi: 10.1007/s10753-014-0013-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, S. L., Lai, S. F., Song, J. Z., Qiao, C. F., Liu, X., Zhou, Y., et al. (2010). Decocting-induced chemical transformations and global quality of du-shen-tang, the decoction of ginseng evaluated by UPLC-Q-TOF-MS/MS based chemical profiling approach. J. Pharm. Biomed. Anal. 53, 946–957. doi: 10.1016/j.jpba.2010.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, A. M., Ping, Y. H., Chang, G. F., Wang, J. Y., Chiu, J. H., Kuo, C. D., et al. (2011). Neuroprotective effect of oral S/B remedy (Scutellaria baicalensis Georgi and Bupleurum scorzonerifolfium Willd) on iron-induced neurodegeneration in the nigrostriatal dopaminergic system of rat brain. J. Ethnopharmacol. 134, 884–891. doi: 10.1016/j.jep.2011.01.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, C. C., and Shieh, D. E. (1996). The anti-inflammatory activity of Scutellaria rivularis extracts and its active components, baicalin, baicalein and wogonin. Am. J. Chin. Med. 24, 31–36. doi: 10.1142/S0192415X96000050

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, C. N., Yuan, Z. G., Zhang, X. L., Yan, R., Zhao, Y. Q., Liao, M., et al. (2012). Saikosaponin a and its epimer saikosaponin d exhibit anti-inflammatory activity by suppressing activation of NF-kappaB signaling pathway. Int. Immunopharmacol. 14, 121–126. doi: 10.1016/j.intimp.2012.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Massey, V. L., Stocke, K. S., Schmidt, R. H., Tan, M., Ajami, N., Neal, R. E., et al. (2015). Oligofructose protects against arsenic-induced liver injury in a model of environment/obesity interaction. Toxicol. Appl. Pharmacol. 284, 304–314. doi: 10.1016/j.taap.2015.02.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Mesko, M. F., Hartwig, C. A., Bizzi, C. A., Pereira, J. S. F., Mello, P. A., and Flores, E. M. M. (2011). Sample preparation strategies for bioinorganic analysis by inductively coupled plasma mass spectrometry. Int. J. Mass Spectr. 307, 123–136. doi: 10.1016/j.ijms.2011.03.002

CrossRef Full Text | Google Scholar

Niemela, O., Parkkila, S., Yla-Herttuala, S., Villanueva, J., Ruebner, B., and Halsted, C. H. (1995). Sequential acetaldehyde production, lipid peroxidation, and fibrogenesis in micropig model of alcohol-induced liver disease. Hepatology 22(4 Pt 1), 1208–1214. doi: 10.1016/0270-9139(95)90630-4

PubMed Abstract | CrossRef Full Text

Patwardhan, R. S., Sharma, D., Thoh, M., Checker, R., and Sandur, S. K. (2016). Baicalein exhibits anti-inflammatory effects via inhibition of NF-kappaB transactivation. Biochem. Pharmacol. 108, 75–89. doi: 10.1016/j.bcp.2016.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Rijpma, A., Meulenbroek, O., van Hees, A. M., Sijben, J. W., Vellas, B., Shah, R. C., et al. (2015). Effects of souvenaid on plasma micronutrient levels and fatty acid profiles in mild and mild-to-moderate Alzheimer’s disease. Alzheimers Res. Ther. 7:51. doi: 10.1186/s13195-015-0134-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Shimizu, T., Shibuya, N., Narukawa, Y., Oshima, N., Hada, N., and Kiuchi, F. (2018). Synergistic effect of baicalein, wogonin and oroxylin a mixture: multistep inhibition of the NF-kappaB signalling pathway contributes to an anti-inflammatory effect of Scutellaria root flavonoids. J. Nat. Med. 72, 181–191. doi: 10.1007/s11418-017-1129-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, S., Dai, J., Fang, J., Gou, X., Cao, H., Zheng, N., et al. (2013). Differences of excess and deficiency zheng in patients with chronic hepatitis B by urinary metabonomics. Evid. Based Complement. Alternat. Med. 2013:738245. doi: 10.1155/2013/738245

PubMed Abstract | CrossRef Full Text | Google Scholar

Sun, S., Dai, J., Wang, W., Cao, H., Fang, J., Hu, Y. Y., et al. (2012). Metabonomic evaluation of ZHENG differentiation and treatment by fuzhenghuayu tablet in hepatitis-B-caused cirrhosis. Evid. Based Complement. Alternat. Med. 2012:453503. doi: 10.1155/2012/453503

PubMed Abstract | CrossRef Full Text | Google Scholar

Tan, M., Schmidt, R. H., Beier, J. I., Watson, W. H., Zhong, H., States, J. C., et al. (2011). Chronic subhepatotoxic exposure to arsenic enhances hepatic injury caused by high fat diet in mice. Toxicol. Appl. Pharmacol. 257, 356–364. doi: 10.1016/j.taap.2011.09.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Tang, J. Y., He, A. H., Jia, G., Liu, G. M., Chen, X. L., Cai, J. Y., et al. (2018). Protective effect of selenoprotein X against oxidative stress-induced cell apoptosis in human hepatocyte (LO2) cells via the p38 pathway. Biol. Trace Elem. Res. 181, 44–53. doi: 10.1007/s12011-017-1025-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J. Y., Lee, C. Y., Pan, P. J., Chang, W. C., Chiu, J. H., Chen, W. S., et al. (2014). Herb-induced autoimmune-like hepatitis in C57BL/6J mice. Liver Int. 34, 583–593. doi: 10.1111/liv.12266

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, L., Zhou, G. B., Liu, P., Song, J. H., Liang, Y., Yan, X. J., et al. (2008). Dissection of mechanisms of chinese medicinal formula realgar-indigo naturalis as an effective treatment for promyelocytic leukemia. Proc. Natl. Acad. Sci. U.S.A. 105, 4826–4831. doi: 10.1073/pnas.0712365105

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Sun, H., Zhang, A., Jiao, G., Sun, W., and Yuan, Y. (2011). Pharmacokinetics screening for multi-components absorbed in the rat plasma after oral administration traditional chinese medicine formula Yin-Chen-Hao-Tang by ultra performance liquid chromatography-electrospray ionization/quadrupole-time-of-flight mass spectrometry combined with pattern recognition methods. Analyst 136, 5068–5076. doi: 10.1039/c1an15752c

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Zhang, A., Wang, P., Sun, H., Wu, G., Sun, W., et al. (2013). Metabolomics coupled with proteomics advancing drug discovery toward more agile development of targeted combination therapies. Mol. Cell Proteomics 12, 1226–1238. doi: 10.1074/mcp.M112.021683

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, Y., Cun, Y., Dong, J., Shao, J., Luo, S., Nie, S., et al. (2010). Polymorphisms in PPARD, PPARG and APM1 associated with four types of traditional Chinese medicine constitutions. J. Genet. Genomics 37, 371–379. doi: 10.1016/S1673-8527(09)60055-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Xi, Y., Wu, M., Li, H., Dong, S., Luo, E., Gu, M., et al. (2015). Baicalin attenuates high fat diet-induced obesity and liver dysfunction: dose-response and potential role of CaMKKbeta/AMPK/ACC pathway. Cell Physiol. Biochem. 35, 2349–2359. doi: 10.1159/000374037

PubMed Abstract | CrossRef Full Text | Google Scholar

Yan, S.-K., Liu, R.-H., Jin, H.-Z., Liu, X.-R., Ye, J., Shan, L., et al. (2015). “Omics” in pharmaceutical research: overview, applications, challenges, and future perspectives. Chin. J. Nat. Med. 13, 3–21. doi: 10.1016/s1875-5364(15)60002-4

CrossRef Full Text | Google Scholar

Ye, H. Z., Zheng, C. S., Xu, X. J., Wu, M. X., and Liu, X. X. (2011). Potential synergistic and multitarget effect of herbal pair chuanxiong rhizome-paeonia albifora pall on osteoarthritis disease: a computational pharmacology approach. Chin. J. Integr. Med. 17, 698–703. doi: 10.1007/s11655-011-0853-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Yohannes, E., Chang, J., Tar, M. T., Davies, K. P., and Chance, M. R. (2010). Molecular targets for diabetes mellitus-associated erectile dysfunction. Mol. Cell Proteomics 9, 565–578. doi: 10.1074/mcp.M900286-MCP200

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, H., Wang, X., Hu, P., Zhou, W., Zhang, M., Liu, J., et al. (2016). Serum metabolomic characterization of liver fibrosis in rats and anti-fibrotic effects of yin-chen-hao-tang. Molecules 21:E126. doi: 10.3390/molecules21010126

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Y., Lin, Y., Zhao, H., Guo, Q., Yan, C., and Lin, N. (2016). Revealing the effects of the herbal pair of euphorbia kansui and glycyrrhiza on hepatocellular carcinoma ascites with integrating network target analysis and experimental validation. Int. J. Biol. Sci. 12, 594–606. doi: 10.7150/ijbs.14151

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Z., Gao, L., Cheng, Y., Jiang, J., Chen, Y., Jiang, H., et al. (2014). Resveratrol, a natural antioxidant, has a protective effect on liver injury induced by inorganic arsenic exposure. Biomed. Res. Int. 2014:617202. doi: 10.1155/2014/617202

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, Z., Barcus, M., Kim, J., Lum, K. L., Mills, C., and Lei, X. G. (2016). High dietary selenium intake alters lipid metabolism and protein synthesis in liver and muscle of pigs. J. Nutr. 146, 1625–1633. doi: 10.3945/jn.116.229955

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, N., Dai, J., Cao, H., Sun, S., Fang, J., Li, Q., et al. (2013). Current understanding on antihepatocarcinoma effects of xiao chai hu tang and its constituents. Evid. Based Complement. Alternat. Med. 2013:529458. doi: 10.1155/2013/529458

PubMed Abstract | CrossRef Full Text | Google Scholar