Original Research ARTICLE
Quantitative Evaluation of the Compatibility Effects of Huangqin Decoction on the Treatment of Irinotecan-Induced Gastrointestinal Toxicity Using Untargeted Metabolomics
- 1MOE Key Laboratory of Drug Quality Control and Pharmacovigilance, China Pharmaceutical University, Nanjing, China
- 2State Key Laboratory of Natural Medicine, China Pharmaceutical University, Nanjing, China
- 3State Key Laboratory for Quality Research in Chinese Medicines, Macau University of Science and Technology, Macau, China
Huangqin decoction (HQD), a traditional Chinese medicine (TCM), has been widely used to treat gastrointestinal syndrome in China for thousands of years. Chemotherapy drug irinotecan (CPT-11) is used clinically to treat various kinds of cancers but limited by its side effects, especially delayed diarrhea. Nowadays, HQD has been proved to be effective in attenuating the intestinal toxicity induced by CPT-11. HQD consists of four medicinal herbs including Scutellaria baicalensis Georgi, Glycyrrhiza uralensis Fisch, Paeonia lactiflora Pall, and Ziziphus jujuba Mill. Due to its complexity, the role of each herb and the multi-herb synergistic effects of the formula are poorly understood. In order to quantitatively assess the compatibility effects of HQD, mass spectrometry-based untargeted metabolomics studies were performed. The serum metabolic profiles of rats administered with HQD, single S. baicalensis decoction, S. baicalensis-free decoction and baicalin/baicalein combination were compared. A time-dependent trajectory upon principal component analysis was firstly used to visualize the overall differences. Then metabolites deregulation score and relative area under the curve were calculated and used as parameters to quantitatively evaluate the compatibility effects of HQD from the aspect of global metabolic profile and the specifically altered metabolites, respectively. The collective results indicated that S. baicalensis played a crucial role in the therapeutic effect of HQD on irinotecan-induced diarrhea. Both HQD and SS decoction regulated glycine, serine and threonine pathway. This study demonstrated that metabolomics was a promising tool to elucidate the compatibility effects of TCM or combinatorial drugs.
Huangqin decoction (HQD) has been widely used to treat gastrointestinal syndrome in China for thousands of years. Irinotecan (CPT-11) is used clinically to treat various kinds of cancers (Ikuno et al., 1995; Souglakos et al., 2006; Kummar et al., 2011) but limited by its side effects like diarrhea, including short duration early secretory diarrhea and life-threatening delayed diarrhea (Sparreboom et al., 1998; Santos et al., 2000; Slatter et al., 2000; Xie et al., 2002). Nowadays, HQD (Wang et al., 2015) and its modified pharmaceutical preparation PHY906 (Farrell and Kummar, 2003; Kummar et al., 2011) have been proved to be effective in attenuating the intestinal toxicity induced by CPT-11 in preclinical and clinical trials.
Huangqin decoction consists of four medicinal herbs including Scutellaria baicalensis Georgi, Glycyrrhiza uralensis Fisch, Paeonia lactiflora Pall, and Ziziphus jujuba Mill. It is believed that the interaction between herbs could lead to synergistic effects (Cheng, 2015). For HQD, researches have revealed that Angelica dahurica extracts had synergistic action on S. baicalensis by enhancing the intestinal absorption of its major bioactive component baicalin in vivo (Liang et al., 2012). This herbal pair also showed synergistic effect on each other in relieving oxidative stress (Yang et al., 2011). Furthermore, S. baicalensis could work synergistically with Coptis chinensis or Pinellia ternata to get more rapid dissolution of the effective ingredients (Takasuna et al., 1995; Feng et al., 2008).
Although it has been confirmed that multi-components would interacted with multi-targets and finally led to synergistic effects in most formulas, the detailed potential mechanisms are still ambiguous. Due to the complexity of chemical composition, only using a single index or a simple superposition method could not be feasible to unveil the multi-herb synergistic characteristics of TCM (Jin et al., 2016). Thus, recently, researchers tried to use knockout/knockin chromatography, bioactive equivalence oriented feedback screening method (Song et al., 2016) and system pharmacology strategy (Huang et al., 2014; Zhang W.J. et al., 2016) to recognize the bioactive equivalent combinatorial components (BECCs) or the most effective chemicals in one formula. As a novel holistic approach to study metabolic alternations, metabolomics was also introduced to study the complex herbal prescriptions (Pelkonen et al., 2012). However, how to quantitatively assess the efficacy of herbal drugs based on metabolomics is still challenging. Most metabolomics studies only compared the separation among control, disease and therapy groups based on PCA score plots visually, and without quantitative basis (Zhang et al., 2010; Li et al., 2016). In order to translate the direct phenotypic graph to quantifiable data, relative mean distance values (RDV) of different groups in PLS-DA score plots were calculated and used to quantitatively assess drug toxicity (Aa et al., 2010). However, the RDV are largely dependent on the number of principal components. A recently proposed method, Pathifier (Drier et al., 2013), which has been successfully used in genomics to quantify the pathway deregulation level in a sample for cancer early diagnosis and stratify (Drier et al., 2013; Livshits et al., 2015; Huang et al., 2016) showed its practicability.
In this study, in order to assess the compatibility effects of HQD on the treatment of irinotecan-induced diarrhea quantitatively, liquid chromatograph-mass spectrometer and gas chromatograph-mass spectrometer based untargeted metabolomics studies were performed. The efficacy effects of HQD, single S. baicalensis decoction, S. baicalensis-free decoction or baicalin/baicalein combination on the treatment of irinotecan-induced diarrhea were evaluated. The overall perturbations observed in serum metabolic profiles were compared. Two parameters, MDS and RAUC, were used to quantitatively evaluate the compatibility effects of HQD from the aspect of global metabolic profile and the specifically altered metabolites, respectively.
Materials and Methods
Chemicals and Reagents
Irinotecan hydrochloride (CPT-11) was obtained from Hengrui (Jiangsu, China). Four component herbs of HQD, Scutellaria baicalensis Georgi (Hebei province), G. uralensis Fisch (Inner Mongolia of China), P. lactiflora Pall (Anhui province), and Z. jujuba Mill (Henan province) were purchased from Tongrentang drugstore (Nanjing, China) and authenticated by Dr. Wei Zhang (State Key Laboratory for Quality Research in Chinese Medicines, Macau University of Science and Technology, China). Baicalin, baicalein, methoxyamine hydrochloride, N-methyl-N-trifluoroacetamide (MSTFA) and pyridine were purchased from Sigma–Aldrich (St. Louis, MO, USA). Methanol, acetonitrile and ethyl acetate for high-performance LC-grade were obtained from Honeywell (Burdick and Jackson, NJ, USA). Deionized water was purified by a Milli-Q system (Millipore, Bedford, MA, USA).
The dry herb pieces of S. baicalensis Georgi (dried roots) (9 g), P. lactiflora Pall (dried roots) (6 g), G. uralensis Fisch (dried and honey-fried roots and rhizomes) (6 g), and Z. jujuba Mill (dried fruits) (6 g) were extracted twice with boiling water (1:15 and then 1:10) for 1.5 and 1 h and filtered through gauze. Then the merged mixtures were concentrated to dryness by a 2.5 L freeze dry system (Labconco, USA). Finally, suspension was prepared by adding 27 ml water (1 g/mL) for intragastric administration (i.g.). Single S. baicalensis and S. baicalensis-free decoctions were prepared through the same procedure with the equal weights of medicine. Besides, according to the ratio in HQD, BB were mixed (13.6:2.5) mg/mL and suspended in water to prepare for the bioactive constituents.
Animal Study and Sampling
Sixty-three male Sprague-Dawley rats (weighing 200 ± 20g) were purchased from the Sino-British SIPPR/BK Lab Animal Ltd. (Shanghai, China). All animals were kept in a standard animal laboratory with regulating temperature and humidity (20–22°C, 45 ± 5%), and 12/12 h light-dark cycle, with free access to chow and water. All experiments were authorized by the Animal Ethics Committee of China Pharmaceutical University and carried out under the Guidelines for Care and Use of Laboratory Animals.
After 1-week accommodation, the rats were randomly divided into six groups: CPT-11 alone (T), HQD treatment (T/HQD), Single S. baicalensis decoction treatment (T/SS), BB combination treatment (T/BB), S. baicalensis-free decoction treatment (T/SF) and control (C). Normal saline or CPT-11 diluted with it was administered intravenously (i.v.) to the animals via the tail vein with 150 mg/kg/day of body weight from days 1 to 2. Treatment groups were intragastrically administrated with HQD, SS decoction, BB combination, or SF decoction from day 0 (about 0.5 h before administration of CPT-11), twice each day for continuous 7 days. Oral and intravenous administration schedules and dosages were shown in Table 1.
Body weight and diarrhea score (DS) were monitored twice per day (before oral gavage) from days -1 to 13. Rats feces were scored using the following scale: 0, normal (normal stool or absent); 1, slight diarrhea (wet and soft stool); 2, moderate diarrhea (wet and unformed stool with moderate perianal staining of the coat); 3, severe diarrhea (watery stool with severe perianal staining of the coat) (Kurita et al., 2000). As shown in Supplementary Figure S1, blood samples from retinal venous plexus were collected at days -1, 1, 4, 7, and 10, then clotted, centrifuged at 8000 rpm for 10 min, and stored at -80°C until GC/MS and LC/MS analysis. Two rats from each group were euthanized by day 4. The jejunum, ileum, cecum and colon were collected to prepare paraffin section and perform hematoxylin-eosin staining for histological examination.
Sample Preparation and Metabolomic Analysis
Sample preparation and analytical method of instruments were based on our previous studies (as shown in Supplementary Information) (Wang et al., 2015; Dai et al., 2016). All tested samples were analyzed randomly. LC/MS analysis was performed on a Shimadzu Prominence series ultrafast liquid chromatography (UFLC) system coupled with ion trap/time-of-flight hybrid mass spectrometry (IT-TOF/MS) (Shimadzu Inc., Kyoto, Japan). The column for separation was Phenomenex Kinetex C18 (100 mm × 2.1 mm, 2.6 μm) (Phenomenex, Torrance, CA, USA). GC/MS analysis was performed on Shimadzu GCMS-QP2010 Ultra (Ultra GC-Q/MS; Shimadzu Inc., Kyoto, Japan) equipped with an Rtx-5MS capillary column (30.0 m × 0.25 mm ID, 0.25 μm).
Data Processing and Statistical Analysis
Peak detection and alignment were performed by profiling solution version 1.1 (Shimadzu, Kyoto, Japan) to get the matrix of peak intensities from the original LC-MS and GC-MS chromatographs, respectively. Following background-peak-filtering, 80 rule, limitation of QCs, normalization by the total ion intensity, ion fusion, and Pareto scaling (Dai et al., 2016; Zhang P. et al., 2016), PCA and OPLS-DA (orthogonal projection to latent structures-discriminant analysis) were conducted by SIMCA-P (version 13.0, Umetrics, Sweden). Parameters of PCA (R2X and Q2) and OPLS-DA (R2X, R2Y and Q2) were used to evaluate the models. Besides, variable importance in the projection (VIP) values in the OPLS-DA model were used to select potential biomarkers. Additionally, Mann–Whitney U test was performed in MeV (version 4.6.11) and p-values were adjusted by false discovery rate correction based on the Benjamini-Hochberg method (Broadhurst and Kell, 2006).
Biomarkers Screening and Identification
For each time point, OPLS-DA models were constructed between group T and other groups, respectively. Features with VIP > 1.0 and p < 0.05 between the two groups were regarded as potential biomarkers and selected for identification.
For LC/MS datasets, peaks were preliminarily identified by referring to online databases (HMDB, METLIN, and Lipid MAPS) or related literatures (ppm < 30). Furthermore, the preliminarily identified metabolites were determined by comparing with the reference standards available in our lab or the fragment models in databases. For GC/MS datasets, metabolites were located in the original chromatographs and preliminarily identified by NIST 11 (National Institute of Standards and Technology). Some important biomarkers were further confirmed by reference standards available in our lab. The potential metabolites were connected by MetaboAnalyst for pathway analysis2.
Data Quality Evaluation in Metabolomics Analysis
For validating the stability of the analysis process, quality control samples (QCs) were first prepared by mixing serum from each test sample (10 μL aliquot). QCs were injected ten times before the batch process and injected one time every eight samples during the analysis process, to monitor the stability of sample preparation and instrument.
The procedure of removing substances originating from the decoctions is as follows. First, combining groups C and T into one group, called CT group. The OPLS-DA model between T/HQD and CT was constructed. Then, trend plots of features in S-plot were shown, features absent in CT were the components from decoction that need to be removed (Liu and Yan, 2013).
Compatibility Effects Evaluation Methods
After serum metabolomics analysis, GC/MS data, LC/MS (+) and LC/MS (-) data were integrated as the whole dataset, all the efficacy evaluations were based on whole dataset. Scheme of compatibility effects evaluation was showed in Figure 1.
First, PCA model based on whole dataset was established. Average value of the first or the second principle components (PC1 and PC2) of every group at each time point was calculated. Each dot represents the average metabolic status at a certain time. They were plotted on one graph with PC1 as X-axis, PC2 as Y-axis and then different group was curved on separate graph so that the results could be displayed more clearly.
Second, R package pathifier was applied to perform metabolites set analysis (Drier et al., 2013). All of the samples dataset were looked as a cloud. We built up the “principal curve” through healthy samples in the cloud. The “principal curve” was firstly put forward as a non-parametric non-linear extension of the linear PCA. It passed through the middle of the data to provide a non-linear summary of the data (Hastie and Stuetzle, 1989). Due to the limited metabolic pathway coverage of untargeted metabolomics approach, we chose to calculated the overall MDS instead of pathway deregulation score (PDS) using the R package pathifier (Huang et al., 2016). Every sample was projected onto the principle curve and it’s MDS was the projection distance value of each sample. The level of sample metabolites deregulation was quantified by MDS. Area under the MDS curves were also calculated.
Third, common re-regulation metabolites in T/HQD, T/SS, and T/BB three groups, T/HQD and T/SF two groups were under investigation. Fold change values relative to group control at each time point compose the relative curves. RAUC from days 4 to 10 was calculated. The lower the RAUC level, the sample is more similar to control.
Body Weight and Histopathology Examination
Comparing to the group C, the relative body weight (RBW), DS and its AUC, reduction of food intake and histopathology examination of group T indicated that CPT-11 induced serious delayed diarrhea which mainly lasted from days 4 to 5. The results were consistent with our previous research (Wang et al., 2015) (Figures 2A–D, 3).
FIGURE 2. Effect of TCM (HQD, SS, BB, or SF decoctions) on CPT-11-induced toxicity in rats. (A) Effects of TCM in protecting against weight loss induced by CPT-11. The change in body weight was calculated on the basis of that on day –1. Each point represents the mean of 8–10 animals. (B) Effects of TCM in protecting against food-intake reduce induced by CPT-11. (C) Effects of TCM in protecting against diarrhea induced by CPT-11. (D) Area under the curve (AUC) of delayed diarrhea score (DS) (days 3.5–8.5). ∗∗∗p < 0.001 vs. T group. Details of experimental procedures are given in section “Materials and Methods.” SS, single S. baicalensis; BB, baicalin and baicalein; SF, S. baicalensis free; T, CPT-11.
FIGURE 3. Effect of TCM (HQD, SS, BB, or SF decoctions) on CPT-11 induced intestinal damage in rats. Hematoxylin and eosin (H&E) staining was used to visualize the formalin-fixed sections of jejunum, ileum, cecum, and colon 3 days after treatment with CPT-11. SS, single S. baicalensis; BB, baicalin and baicalein; SF, S. baicalensis free; T, CPT-11.
Meanwhile, the serious delayed diarrhea was partly alleviated in decoction treated groups (T/HQD, T/SS, T/BB, and T/SF) with different degrees. Taken the DS, RBW, incidence of severe diarrhea and histological examination in day 4 and/or 5 into account, the treatment effects of these groups were sorted as HQD > SS decoction > BB combination > SF decoction (Table 2). These results were indicative, but not conclusive since the difference of treatment effect among these decoctions were not significant. Additionally, DS and histological examination was not so objective and RBW was not so closely related with serious delayed diarrhea.
Data Quality Evaluation in Metabolomics Analysis
Aiming to detect as many endogenous metabolites as possible, full scan mode was performed in sample analysis which meant all substances in serum were detected. So, the exogenous components derived from HQD were also included in the dataset. In order to eliminate the influence of those exogenous components, chemical constituents originating from the HQD were recognized and removed first. Taken the constituents in HQD as an example, OPLS-DA was conducted between non-decoction treated groups (C and T) and HQD treated group (R2X = 0.767, R2Y = 0.971, Q2 = 0.766). S-plot and trend plot were used to screen out decoction-originated compounds. For example, peak 102 was filtered out by S-plot, and trend plot demonstrated it was originated from HQD (Supplementary Figure S2). Finally, 53 and 159 peaks from GC/MS and LC/MS were retained, respectively.
As shown in Supplementary Figure S3, QC samples were close in PCA score plots. Furthermore, the retention time shift was less than 0.1 min and the RSD values of all peaks in QCs were lower than 15%, which demonstrated the whole analytical process was stable and reliable (Naz et al., 2014).
Time-dependent Trajectory Based Qualitative Evaluation
To construct a model containing overall metabolites information, the datasets from LC/MS (ESI+), LC/MS (ESI-) and GC/MS were combined and analyzed by a PCA model. In the PCA model of all test samples (R2X = 0.919, Q2 = 0.591), the trajectory change of group C in different days was relatively small, while the deviation in group T was obvious. Comparing with day -1 in group T, the distances of other days were correlated to DS and RBW: the distance of day 4 (i.e., the day of delayed diarrhea) was the longest, then it was greatly reduced at day 7 and finally it was close to zero at day 10 (Figure 4). Also, the results of PCA indicated the differences between group C and group T at days 4, 7, and 10 were related to the delayed diarrhea (Supplementary Figure S3). Meanwhile, the trajectory analysis indicated the differences between group C and decoction-treated groups during delayed diarrhea were less obvious than that between group C and group T.
FIGURE 4. The time-dependent trajectory of the metabolic profiles of control, T and TCM (HQD, SS, BB, or SF decoctions) treatment groups. Each dot represents the average metabolic status at a certain time. (A) Summary graph. (A-1) to (A-6) The trajectory change for each group. SS, single S. baicalensis; BB, baicalin and baicalein; SF, S. baicalensis free; T, CPT-11.
Particularly at day 4, the deviations of group T/HQD and T/SS were much less than that in group T. OPLS-DA score plots of day 4 showed the metabolic profiles of group T/HQD, T/SS and T/BB were separated from group C and group T, while group T/SF were mixed with group T (Supplementary Figure S4). These results visualized the time-dependent trajectory of CPT-11 exposure and decoction treatment effect from the aspect of global metabolic alternation.
Quantitative Evaluation of the Compatibility Effects
Using the biomarkers selection method described above, totally 47 discriminating metabolites (VIP > 1 and p < 0.05) were screened out and identified. Based on Spearman correlation analysis between altered metabolites and DS/RBW, four metabolites were removed because of the weak correlation (-0.5 < Spearman correlation coefficients < 0.5) (Supplementary Figure S5). The remaining 43 metabolites were used to calculate MDS.
Principal curve and data-cloud of samples were shown in Supplementary Figure S6. In general, the MDS demonstrated that the CPT-11 induced metabolites deregulation on day 1, further (the most) on day 4, and then the abnormalities were almost recovered on days 7 and 10. While decoction-treated groups (T/HQD, T/SS, T/BB, and T/SF) were closer to control. Furthermore, MDS curves and the areas under the curves (Figures 5A,B) evaluated the treatment effects of decoctions quantitatively. At day 4, group T, T/SF, or T/BB were significantly different (p < 0.05) from group T/HQD and T/SS. However, groups T/HQD and T/SS had no significant difference between each other (Supplementary Figure S7).
FIGURE 5. Relative quantitation of TCM efficacy based on MDS. (A) Time dependent MDS curves of CPT-11 and TCM (HQD, SS, BB, or SF decoctions) treatment groups. (B) Area under the curve (AUC) of MDS curves. ∗∗p < 0.01, ∗∗∗p < 0.001 vs. T group; #p < 0.05 vs. SF group. SS, single S. baicalensis; BB, baicalin and baicalein; SF, S. baicalensis free; T, CPT-11.
On the overall MDS basis, we further compared difference between decoction treatment groups and group T based on re-regulation metabolites. The relative concentration of creatine, creatinine, glycine, LysoPE (20:4), proline, valine, nonanedioic acid, acetylcarnitine, citric acid and methionine were commonly re-regulated in T/HQD, T/SS, and T/BB treatment groups. Creatine, creatinine, glycine, LysoPE (20:4), proline, valine, nonanedioic acid and alanine were re-regulated in both T/HQD and T/SF treatment groups. Creatine, creatinine, glycine, LysoPE (20:4), proline and valine are the common re-regulation metabolites.
Relative area under the curve results demonstrated that groups T/HQD, T/SS and T/BB had obvious variance. For group T/HQD, RAUC of creatine, creatinine, LysoPE (20:4), proline, valine, citric acid and methionine were significantly lower than that of group T. As to group T/SS, RAUC of creatinine, proline, valine, and methionine were significantly lower than group T. Except for valine, other re-regulated metabolites in the T/BB group were the same with T/SS group (Figure 6A). In addition, five common re-regulated metabolites were significantly restored in group T/HQD (Figure 6B). However, only LysoPE (20:4) was significantly restored in group T/SF.
FIGURE 6. Re-regulated metabolites’ RAUC from days 4 to 10. (A) T/HQD, T/SS, T/BB vs. T. (B) T/HQD, T/SF vs. T. (Control is the zero level of Y-axis). ∗p < 0.05, ∗∗p < 0.01 vs. T group [not false discovery rate (FDR)-adjusted p-values]. SS, single S. baicalensis; BB, baicalin and baicalein; SF, S. baicalensis free; T, CPT-11.
All above results indicated that CPT-11 induced distinct change of metabolites in days 4–10, and groups T/HQD, T/SS, and T/BB showed retracement trend. Moreover, group T/HQD showed the most obvious trend of retracement in common metabolites. Using RAUC as an indicative parameter, the difference in treatment effect of each TCM on CPT-11 induced intestinal toxicity could be quantified from the aspect of the specific change of individual differential metabolite.
After biomarkers identification and spearman correlation analysis, 43 metabolites were reserved. All of the reserved 43 metabolites were significantly perturbed at least in 1 day during the process of CPT-11 exposure and listed in Supplementary Tables S1, S2. Heat-maps showed the time-dependent fluctuations (red = increased concentration, green = decreased concentration) of metabolites as a result of CPT-11 perturbation and decoction deregulation (Supplementary Figure S8). Except for leucine, lysine, threonine, tryptophan, glutamic acid, uric acid, glycocholic acid, linoleic acid, LysoPC (18:0), LysoPC (20:3) and LysoPE (16:0), other metabolites were re-regulated at least 1 day in decoction treatment groups (T/HQD, T/SS, T/BB, or T/SF).
It has been reported that five metabolic pathways in serum were identified associated with CPT-11 exposure (phenylalanine, tyrosine and tryptophan biosynthesis; phenylalanine metabolism; glycine, serine and threonine metabolism; alanine, aspartate and glutamate metabolism; and primary bile acid biosynthesis) and HQD attenuated the side-effects with reversed glutamine, tryptophan metabolism, and lipid metabolism (Wang et al., 2015). In the present study, D-glutamine and D-glutamate metabolism, arginine and proline metabolism, and valine, leucine and isoleucine biosynthesis pathways are also considered to be relevant pathways influenced by CPT-11. These metabolism changes are related to energy, amino acid and lipids metabolism. Furthermore, we also found that decoction treatment groups achieved a comprehensive efficacy via multiple pathways and targets (Figure 7). The number and extend of re-regulated metabolites in groups T/SS, T/BB, and T/SF are all less than those in group T/HQD.
FIGURE 7. The perturbed metabolites and metabolic pathways related to CPT-11 treatment as well as TCM (HQD, SS, BB, or SF decoctions) modification. Elevation (up arrows) and reduction (down arrows) of the levels of metabolites observed in rats are indicated. SS, single S. baicalensis; BB, baicalin and baicalein; SF, S. baicalensis free; T, CPT-11.
It is noted that HQD alleviated the delayed diarrhea induced by CPT-11 mainly via modulating glycine, serine and threonine metabolism, and phenylalanine, tyrosine and tryptophan biosynthesis pathways. Besides, glycine, serine and threonine metabolism, and arginine and proline metabolism pathways play an important role in SS decoction mitigates the delayed diarrhea.
Qualitative Analysis of TCM Efficacy
Traditional Chinese medicine formulas are an important source for the development of new drugs to prevent and treat gastrointestinal toxicity induced by chemotherapeutics (Liu and Cheng, 2012). In studies of HQD, its component S. baicalensis was described in Bencao Gangmu (Compendium of Materia Medica) about 420 years ago for the treatment of diarrhea, dysentery, and inflammation (Li, 1593/2012). Moreover, as the major bioactive constituents in S. baicalensis, baicalin and its aglycones baicalein can protect injured intestinal epithelial cells (Yao et al., 2016) and active autophagy in cancer cells (Li-Weber, 2009). The mechanism research further demonstrated that baicalin might be able to promote the polarization of macrophages to an M2 phenotype to relieve ulcerative colonic inflammation (Zhu et al., 2016). Afterward, other studies also shown that BB could protect intestinal epithelial cells and inhibit mTORC1 in human cancer cells (Aryal et al., 2014). And so, we think it’s reasonable to speculate that S. baicalensis and its active component BB (BB) have ameliorative effect on CPT-11-induced gastrointestinal toxicity. Because baicalin is absorbed as the form of baicalein and they have “enterohepatic circulation” phenomenon in vivo (Li et al., 2012), rats in group T/BB were administrated the mixture of BB. Mean metabolic position had been used in PCA trajectory plots to map the metabolic changes in chemicals exposure by Nicholson et al. (2002). So, we studied the serum metabolome to reveal trajectories of PCA score plot derived from C, T and TCM treatment groups.
The results showed that HQD achieved the best efficacy with alleviated intestinal tissue pathological damage, SF decoction was the lowest one, and SS decoction and BB combination were in between. The results not only supported that CPT-11 induced metabolic profiles changes in a time-dependent manner but also revealed the overall time-dependent dynamic changes of difference groups. These observations were in agreement with the theory of the compatibility of HQD.
In the development of herbal therapy, single herb formulae were gradually evolved into the multi-herb formulae. Modernization of TCM focuses on multi-components formulae (Wang et al., 2012). Consequently, the results of this study may provide some heuristic guidance for TCM formulation rationalization and optimization.
Relatively Quantitative Analysis of TCM Efficacy Based on MDS
We developed a new method to quantitatively evaluate the efficacy of TCM referring to PDS. PDS was first proposed in genomics study, and was used to transforms gene-level information into pathway-level information for cancer stratification and diagnosis (Drier et al., 2013; Huang et al., 2016). In the light of this idea, we calculated the overall MDS instead of PDS to evaluate the level of global metabolic alteration quantitatively. Using the method presented above, we selected the best relevant features, screened significantly changed metabolites, and then calculated the overall MDS at each time points to assess the compatibility effects dynamic changes of serum profile. It is especially interesting that MDS of groups T/HQD and T/SS have no significant difference. However, the MDS of groups T/SF and T are significantly higher at day 4.
These results demonstrated that, S. baicalensis played a major role in anti-CPT-11 toxicity effects among the four herb constituents of HQD formula. The reason leading to this results might be that HQD’s targeted facilitated intestinal recovery function could not be reflected completely in blood (Lam et al., 2010).
As to the RAUCs of these common differential metabolites, it is especially interesting that proline, methionine and creatinine are related to T/HQD, T/SS, and T/BB three treatment groups. A recent study demonstrated that CPT-11 triggered inflammation, induced apoptosis and antioxidant, changed the profile of intestinal bacterial (Lam et al., 2010, 2014). Methionine is an endogenous antioxidant which plays an essential role in protection against oxidative damage (Luo and Levine, 2009). Other studies also revealed that methionine made a contribution to the stability of the Escherichia coli (Oshima and Miwa, 2016). Proline and transmembrane protein were both involved in the intestinal mucosa epithelium mucin secretion, to ensure the integrity of the mucous barrier (Wu, 1998). The gastrointestinal mucosal barrier played a crucial role in the regulation of intestinal secretion, absorption (Wu, 1998) and intestinal immunity (Wu, 2009). Serum creatinine is a traditional indicator of kidney injury (Wyss and Kaddurah-Daouk, 2000). The up-regulated of creatinine in the group T was in consistent with our previous research that CPT-11 brought out renal injury (Yao et al., 2017). Thus, we speculated that HQD, SS decoction and BB combination could contribute to maintain the stability and integrity of the intestinal tract, and mitigate the damage of tissues.
Altered Pathways Related to HQD and SS Decoction Efficacy
Intriguingly, both HQD and SS decoction elevated the level of glycine, which belongs to glycine, serine, and threonine pathway and it is highly associated with intestinal mucosal and barrier function. It has been demonstrated that CPT-11 caused destruction of mucosal architecture, increased the number of apoptotic cell, and triggered inflammation (Lam et al., 2010).
It has been established that intestinal inflammation increased threonine utilization and caused threonine deficiency (Remond et al., 2009). Furthermore, threonine is an necessary amino acid for the mucous layer (Faure et al., 2005; Mao et al., 2011). Thus, threonine is particularly important for the maintenance of gut integrity. Glycine-serine conversion is strictly related to N5N10methylentetrahydrofolate synthesis (N5N10-methylen-FH4 is a coenzyme in the purine and pyrimidine synthesis process). Purine and pyrimidine are essential substances for intestinal epithelial cell proliferation and protein synthesis (Tedeschi et al., 2013; Wang et al., 2014). The content of threonine and glycine were decreased in group T and re-regulated in groups T/HQD and T/SS. These observations supported the result showed here that HQD and SS decoction alleviated the delayed diarrhea via modulating glycine, serine and threonine metabolism.
The previous literature has showed that both phenylalanine, tyrosine, and tryptophan biosynthesis pathway and arginine and proline metabolism pathway are related to inflammation reaction. Significant increase in serum phenylalanine and the phenylalanine-tyrosine ratio and inducible nitric oxide synthase (iNOS could catalyze the synthesis of arginine) mRNA are indicators of the inflammatory disease magnitude (Wannemacher et al., 1976; Lam et al., 2010). A significant down-regulation of phenylalanine was found only in HQD treatment group, implying the efficacy specificity of HQD probably resulted from regulated phenylalanine, tyrosine, and tryptophan biosynthesis pathway. It is worth pointing out here that, HQD also elevated content of proline although the arginine and proline metabolism pathway is one of the most important pathways regulated by SS. We speculate that SS decoction might re-regulated other metabolites in this pathway.
From the above, we could draw the conclusion that both HQD and SS decoction could protect the gastrointestinal tract by altering the glycine, serine and threonine metabolism pathway, and HQD altering the phenylalanine, tyrosine and tryptophan biosynthesis pathway specifically. Due to their multi-component and multi-targets, the potential mechanism of HQD, SS decoction in treating severe diarrhea induced by CPT-11 deserves further investigation and discussion.
In this report, we proposed a method to quantitatively assess the efficacy of HQD, which would be further used to elucidate the compatibility effects of complex formulae. Time-dependent trajectory and MDS combine RAUC based quantitative evaluation all showed distinct changes in various compatibility. Collectively, these results indicated that S. baicalensis played a crucial role in HQD on the treatment of irinotecan-induced diarrhea and exerted nearly the same efficacy effect with HQD in the aspect of global metabolic alteration. Both HQD and SS decoction could protect the gastrointestinal tract mainly by altering the glycine, serine and threonine metabolism pathway. Our study also demonstrated that metabolomics could be served as a promising tool to elucidate the compatibility effects of TCM or combinatorial drug.
D-NC carried out most of the studies, performed statistical analysis and wrote the manuscript. XW and J-QC participated in the animal experiments. BL participated in the data processing work. WZ provided professional advices. Z-JZ and F-GX designed the study and revised the manuscript. All authors have read and approved the final version.
This study was financially supported by the National Science Foundation of China (No.81302733, 81573385), the research project of Chinese Ministry of education (No.113036A), the Program for Jiangsu province Innovative Research Team, the Program for New Century Excellent Talents in University (No.NCET-13-1036), a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD).
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.
We would like to thank Mr. Yi-Qiao Gao for his warm encouragement and scientific support.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fphar.2017.00211/full#supplementary-material
AUC, area under the curve; BB, baicalin and baicalein; C, control; GC/MS, gas chromatography-mass spectrometry; HQD, Huangqin decoction; i.g., intragastric; i.v., intravenous; LC/MS, liquid chromatography-mass spectrometry; MDS, metabolites deregulation score; OPLS-DA, orthogonal partial least squares-discriminant analysis; PCA, principle component analysis; QC, quality control; RAUC, relative area under the curve; SF, S. baicalensis free; SS, single S. baicalensis; T, CPT-11; TCM, traditional Chinese medicine.
Aa, J., Shao, F., Wang, G., Huang, Q., Zha, W., Yan, B., et al. (2010). Gas chromatography time-of-flight mass spectrometry based metabolomic approach to evaluating toxicity of triptolide. Metabolomics 7, 217–225. doi: 10.1007/s11306-010-0241-8
Aryal, P., Kim, K., Park, P. H., Ham, S., Cho, J., and Song, K. (2014). Baicalein induces autophagic cell death through AMPK/ULK1 activation and downregulation of mTORC1 complex components in human cancer cells. FEBS J. 281, 4644–4658. doi: 10.1111/febs.12969
Dai, D., Gao, Y., Chen, J., Huang, Y., Zhang, Z., and Xu, F. (2016). Time-resolved metabolomics analysis of individual differences during the early stage of lipopolysaccharide-treated rats. Sci. Rep. 6:34136. doi: 10.1038/srep34136
Farrell, M. P., and Kummar, S. (2003). Phase I/IIA randomized study of PHY906, a novel herbal agent, as a modulator of chemotherapy in patients with advanced colorectal cancer. Clin. Colorectal Cancer 2, 253–256. doi: 10.3816/CCC.2003.n.007
Faure, M., Moennoz, D., Montigon, F., Mettraux, C., Breuille, D., and Ballevre, O. (2005). Dietary threonine restriction specifically reduces intestinal mucin synthesis in rats. J. Nutr. 135, 486–491.
Feng, Y.-B., Luo, W.-Q., and Zhu, S.-Q. (2008). [Explore new clinical application of Huanglian and corresponding compound prescriptions from their traditional use]. Zhongguo Zhong Yao Za Zhi 33, 1221–1225.
Huang, C., Zheng, C. L., Li, Y., Wang, Y. H., Lu, A. P., and Yang, L. (2014). Systems pharmacology in drug discovery and therapeutic insight for herbal medicines. Brief. Bioinform. 15, 710–733. doi: 10.1093/bib/bbt035
Huang, S., Chong, N., Lewis, N. E., Jia, W., Xie, G., and Garmire, L. X. (2016). Novel personalized pathway-based metabolomics models reveal key metabolic pathways for breast cancer diagnosis. Genome Med. 8, 34. doi: 10.1186/s13073-016-0289-9
Ikuno, N., Soda, H., Watanabe, M., and Oka, M. (1995). Irinotecan (Cpt-11) and characteristic mucosal changes in the mouse ileum and cecum. J. Natl. Cancer Inst. 87, 1876–1883. doi: 10.1093/jnci/87.24.1876
Jin, Y., Qu, C., Tang, Y. P., Pang, H. Q., Liu, L. L., Zhu, Z. H., et al. (2016). Herb pairs containing Angelicae Sinensis Radix (Danggui): a review of bio-active constituents and compatibility effects. J. Ethnopharmacol. 181, 158–171. doi: 10.1016/j.jep.2016.01.033
Kummar, S., Copur, M. S., Rose, M., Wadler, S., Stephenson, J., O’rourke, M., et al. (2011). A phase I study of the Chinese herbal medicine PHY906 as a modulator of irinotecan-based chemotherapy in patients with advanced colorectal cancer. Clin. Colorectal Cancer 10, 85–96. doi: 10.1016/j.clcc.2011.03.003
Kurita, A., Kado, S., Kaneda, N., Onoue, M., Hashimoto, S., and Yokokura, T. (2000). Modified irinotecan hydrochloride (CPT-11) administration schedule improves induction of delayed-onset diarrhea in rats. Cancer Chemother. Pharmacol. 46, 211–220. doi: 10.1007/s002800000151
Lam, W., Bussom, S., Guan, F., Jiang, Z., Zhang, W., Gullen, E. A., et al. (2010). The four-herb Chinese medicine PHY906 reduces chemotherapy-induced gastrointestinal toxicity. Sci. Transl. Med. 2:45ra59. doi: 10.1126/scitranslmed.3001270
Lam, W., Jiang, Z. L., Guan, F. L., Hu, R., Liu, S. H., Chu, E., et al. (2014). The number of intestinal bacteria is not critical for the enhancement of antitumor activity and reduction of intestinal toxicity of irinotecan by the Chinese herbal medicine PHY906 (KD018). BMC Complement. Altern. Med. 14:490. doi: 10.1186/1472-6882-14-490
Li, C. R., Zhang, L., Wo, S. K., Zhou, L. M., Lin, G., and Zuo, Z. (2012). Pharmacokinetic interactions among major bioactive components in Radix Scutellariae via metabolic competition. Biopharm. Drug Dispos. 33, 487–500. doi: 10.1002/bdd.1815
Li, J., Yang, H., Meng, S., Zhou, J., and Ding, Y. (2016). Metabonomic profiles reveal dose-dependent effects of Bu-Shen-Gu-Chi-Wan on the serum in experimental periodontitis of rat model. J. Ethnopharmacol. 193, 248–254. doi: 10.1016/j.jep.2016.07.051
Liang, X. L., Liao, Z. G., Zhu, J. Y., Zhao, G. W., Yang, M., Yin, R. L., et al. (2012). The absorption characterization effects and mechanism of Radix Angelicae dahurica extracts on baicalin in Radix Scutellariae using in vivo and in vitro absorption models. J. Ethnopharmacol. 139, 52–57. doi: 10.1016/j.jep.2011.10.001
Li-Weber, M. (2009). New therapeutic aspects of flavones: the anticancer properties of Scutellaria and its main active constituents Wogonin, Baicalein and Baicalin. Cancer Treat. Rev. 35, 57–68. doi: 10.1016/j.ctrv.2008.09.005
Pelkonen, O., Pasanen, M., Lindon, J. C., Chan, K., Zhao, L., Deal, G., et al. (2012). Omics and its potential impact on R&D and regulation of complex herbal products. J. Ethnopharmacol. 140, 587–593. doi: 10.1016/j.jep.2012.01.035
Remond, D., Buffiere, C., Godin, J. P., Mirand, P. P., Obled, C., Papet, I., et al. (2009). Intestinal inflammation increases gastrointestinal threonine uptake and mucin synthesis in enterally fed minipigs. J. Nutr. 139, 720–726. doi: 10.3945/jn.108.101675
Slatter, J. G., Schaaf, L. J., Sams, J. P., Feenstra, K. L., Johnson, M. G., Bombardt, P. A., et al. (2000). Pharmacokinetics, metabolism, and excretion of irinotecan (CPT-11) following i.v. infusion of C-14 CPT-11 in cancer patients. Drug Metab. Dispos. 28, 423–433.
Song, H. P., Wu, S. Q., Hao, H., Chen, J., Lu, J., Xu, X., et al. (2016). A chemical family-based strategy for uncovering hidden bioactive molecules and multicomponent interactions in herbal medicines. Sci. Rep. 6:23840. doi: 10.1038/srep23840
Souglakos, J., Androulakis, N., Syrigos, K., Polyzos, A., Ziras, N., Athanasiadis, A., et al. (2006). FOLFOXIRI (folinic acid, 5-fluorouracil, oxaliplatin and irinotecan) vs FOLFIRI (folinic acid, 5-fluorouracil and irinotecan) as first-line treatment in metastatic colorectal cancer (MCC): a multicentre randomised phase III trial from the Hellenic Oncology Research Group (HORG). Br. J. Cancer 94, 798–805. doi: 10.1038/sj.bjc.6603011
Sparreboom, A., De Jonge, M. J. A., De Bruijn, P., Brouwer, E., Nooter, K., Loos, W. J., et al. (1998). Irinotecan (CPT-11) metabolism and disposition in cancer patients. Clin. Cancer Res. 4, 2747–2754.
Takasuna, K., Kasai, Y., Kitano, Y., Mori, K., Kobayashi, R., Hagiwara, T., et al. (1995). Protective effects of kampo medicines and baicalin against intestinal toxicity of a new anticancer camptothecin derivative, irinotecan hydrochloride (cpt-11), in rats. Jpn. J. Cancer Res. 86, 978–984. doi: 10.1111/j.1349-7006.1995.tb03010.x
Tedeschi, P. M., Markert, E. K., Gounder, M., Lin, H. X., Dolfi, S. C., Chan, L. L. Y., et al. (2013). Contribution of serine, folate, and glycine metabolism to the ATP, NADPH, and purine requirements of cancer cells. Cell Death Dis. 4:e877. doi: 10.1038/cddis.2013.393
Wang, J., Fan, H., Wang, Y., Wang, X., Zhang, P., Chen, J., et al. (2015). Metabolomic study of Chinese medicine Huang Qin decoction as an effective treatment for irinotecan-induced gastrointestinal toxicity. RSC Adv. 5, 26420–26429. doi: 10.1039/C5RA02581H
Wang, S. P., Hu, Y. Y., Tan, W., Wu, X., Chen, R. E., Cao, J. L., et al. (2012). Compatibility art of traditional Chinese medicine: from the perspective of herb pairs. J. Ethnopharmacol. 143, 412–423. doi: 10.1016/j.jep.2012.07.033
Wang, W. W., Wu, Z. L., Lin, G., Hu, S. D., Wang, B., Dai, Z. L., et al. (2014). Glycine stimulates protein synthesis and inhibits oxidative stress in pig small intestinal epithelial cells. J. Nutr. 144, 1540–1548. doi: 10.3945/jn.114.194001
Wannemacher, R. W. Jr., Klainer, A. S., Dinterman, R. E., and Beisel, W. R. (1976). The significance and mechanism of an increased serum phenylalanine-tyrosine ratio during infection. Am. J. Clin. Nutr. 29, 997–1006.
Xie, R., Mathijssen, R. H., Sparreboom, A., Verweij, J., and Karlsson, M. O. (2002). Clinical pharmacokinetics of irinotecan and its metabolites: a population analysis. J. Clin. Oncol. 20, 3293–3301. doi: 10.1200/JCO.2002.11.073
Yang, P., Kong, X., Cheng, C., Li, C., Yang, X., and Zhao, X. E. (2011). Synthesis and biological evaluation of 8-substituted and deglucuronidated scutellarin and baicalin analogues as antioxidant responsive element activators. Sci. China Chem. 54, 1565–1575. doi: 10.1007/s11426-011-4361-4
Yao, J., Cao, X., Zhang, R., Li, Y. X., Xu, Z. L., Zhang, D. G., et al. (2016). Protective effect of baicalin against experimental colitis via suppression of oxidant stress and apoptosis. Pharmacogn. Mag. 12, 225–234. doi: 10.4103/0973-1296.186342
Yao, Y., Zhang, P., Wang, J., Chen, J., Wang, Y., Huang, Y., et al. (2017). Dissecting target toxic tissue and tissue specific responses of irinotecan in rats using metabolomics approach. Front. Pharmacol. 8:122. doi: 10.3389/fphar.2017.00122
Zhang, P., Chen, J. Q., Wang, Y., Huang, Y., Tian, Y., Zhang, Z. J., et al. (2016). Discovery of potential biomarkers with dose- and time-dependence in cisplatin-induced nephrotoxicity using metabolomics integrated with a principal component-based area calculation strategy. Chem. Res. Toxicol. 29, 776–783. doi: 10.1021/acs.chemrestox.5b00519
Zhang, W. J., Bai, Y. F., Wang, Y. H., and Xiao, W. (2016). Polypharmacology in drug discovery: a review from systems pharmacology perspective. Curr. Pharm. Des. 22, 3171–3181. doi: 10.2174/1381612822666160224142812
Zhu, W., Jin, Z. S., Yu, J. B., Liang, J., Yang, Q. D., Li, F. J., et al. (2016). Baicalin ameliorates experimental inflammatory bowel disease through polarization of macrophages to an M2 phenotype. Int. Immunopharmacol. 35, 119–126. doi: 10.1016/j.intimp.2016.03.030
Keywords: Huangqin decoction, delayed diarrhea, metabolomics, metabolites deregulation score, relative area under the curve
Citation: Cui D-N, Wang X, Chen J-Q, Lv B, Zhang P, Zhang W, Zhang Z-J and Xu F-G (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
Received: 25 January 2017; Accepted: 05 April 2017;
Published: 21 April 2017.
Edited by:Thomas Efferth, Johannes Gutenberg-Universität Mainz, Germany
Reviewed by:Pierluigi Caboni, University of Cagliari, Italy
Feng Li, Baylor College of Medicine, USA
Jeremy Everett, University of Greenwich, UK
Copyright © 2017 Cui, Wang, Chen, Lv, Zhang, Zhang, Zhang and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.