Core Bioactive Components Promoting Blood Circulation in the Traditional Chinese Medicine Compound Xueshuantong Capsule (CXC) Based on the Relevance Analysis between Chemical HPLC Fingerprint and In Vivo Biological Effects

Compound xueshuantong capsule (CXC) is an oral traditional Chinese herbal formula (CHF) comprised of Panax notoginseng (PN), Radix astragali (RA), Salvia miltiorrhizae (SM), and Radix scrophulariaceae (RS). The present investigation was designed to explore the core bioactive components promoting blood circulation in CXC using high-performance liquid chromatography (HPLC) and animal studies. CXC samples were prepared with different proportions of the 4 herbs according to a four-factor, nine-level uniform design. CXC samples were assessed with HPLC, which identified 21 components. For the animal experiments, rats were soaked in ice water during the time interval between two adrenaline hydrochloride injections to reduce blood circulation. We assessed whole-blood viscosity (WBV), erythrocyte aggregation and red corpuscle electrophoresis indices (EAI and RCEI, respectively), plasma viscosity (PV), maximum platelet aggregation rate (MPAR), activated partial thromboplastin time (APTT), and prothrombin time (PT). Based on the hypothesis that CXC sample effects varied with differences in components, we performed grey relational analysis (GRA), principal component analysis (PCA), ridge regression (RR), and radial basis function (RBF) to evaluate the contribution of each identified component. Our results indicate that panaxytriol, ginsenoside Rb1, angoroside C, protocatechualdehyde, ginsenoside Rd, and calycosin-7-O-β-D-glucoside are the core bioactive components, and that they might play different roles in the alleviation of circulation dysfunction. Panaxytriol and ginsenoside Rb1 had close relevance to red blood cell (RBC) aggregation, angoroside C was related to platelet aggregation, protocatechualdehyde was involved in intrinsic clotting activity, ginsenoside Rd affected RBC deformability and plasma proteins, and calycosin-7-O-β-D-glucoside influenced extrinsic clotting activity. This study indicates that angoroside C, calycosin-7-O-β-D-glucoside, panaxytriol, and protocatechualdehyde may have novel therapeutic uses.


Introduction
Traditional Chinese medicine (TCM) employs compounds called Chinese herbal formulas (CHFs). In the past three decades, these compounds have been attracting increasing attention for their complementary therapeutic effects to western medicines [1,2]. CHFs are comprised of multiple components and affect numerous targets, yet the bioactive components of most CHFs have not been elucidated. In the Pharmacopoeia of the People's Republic of China (2010), only a few components, or even a single component, based on their content in herbs instead of their clinical activities, were controlled in current CHF quality standards. In order to better guarantee the clinical safety and effectiveness of CHF, the bioactive components evidenced by suitable methods should be controlled.
Bioactivity control, which reflects information directly related to safety and effectiveness, is indispensable in CHF quality control [3]. Therefore, a change from traditional evaluation methods that emphasize chemical characteristics to a method that quantifies biological effects is urgently required [4]. Also, the Guidance for Industry Botanical Drug Products from the U.S. Food & Drug Administration (FDA) (2004) clearly indicated that the functional study of known components in herbal medicine should be improved. The recognition and control of bioactive components has thus become a key technology of CHF development.
According to the theory of modern chemical biology, drug effects are achieved by the regulation of biomacromolecules by small molecules in it. Hence, there must be a group of bioactive components responsible for a drug's clinical effects. The concept of relevance between drug components and effects was first introduced in 2002 [5], and the earliest methodology to match drug effects with prominent peaks in high-performance liquid chromatography (HPLC) fingerprints was reported in 2000 [6]. Kong et al. studied the relationship between the HPLC fingerprint and the antibacterial effects of artificial Calculus bovis with the help of chemometric methods [7]. Lu et al. examined the corresponding relationship between the HPLC fingerprint and the effects of Houttuynia cordata injection using multivariate statistical analysis [8]. In the past decade, studies on the relevance between drug components and effects have provided valuable information to improve TCM development, including quality control, classification, CHF prescription analysis and optimization, the identification of new medicinal plant sources, and drug design. The concept of blood stasis is described in TCM theory as a slowing or pooling of the blood due to disruption of the heart's Qi [9]. It is often understood in biomedical terms as hematological disorders, such as hemorrhage, congestion, thrombosis, and local ischemia [9]. In 1998, Mchedlishvili et al reported that decreased blood flow velocity disturbs normal blood flow, which indicates hemorheological abnormalities [10]. It was explained by pathology as a state resulting from a sluggish or impeded blood flow in the body, or abnormal blood outside the vessels that remains in the body and fails to disperse [11]. Once blood stasis develops, the circulation will be further affected, resulting in new pathological changes [11,12].
The core bioactive components of CXC that improve blood circulation impairment remain unknown. The present study was designed to assess the relationship between CXC components and circulatory effects and the specific research process was shown in Fig. 1.

Ethics statement
We obtained 130 specific-pathogen free (SPF) male Sprague-Dawley rats, weighing 180-220 g, from the Guangdong medical experimental animal center (certification No: SCXK-(Yue) 2008-0002) and raised in an SPF laboratory (license: SCXK-(Yue) 2009-0020) at Sun Yat-sen University. This experiment was approved by Animal Ethics Committee of the Life Sciences School, Sun Yatsen University (permission No: 2012082502). The temperature of the SPF laboratory was 20-23uC, and the relative humidity was 50-65%. Rats were fed with standard pellet feed. Experiments  began after the rats had adapted to the new environment for one week. The harm to rats was minimized during the experimental process.

CXC sample preparation
Based on the original formula, CXC samples were prepared by changing the content of 4 herbs. The specific weights of 4 herbs were designed according to a four-factor, nine-level uniform design (UD) ( Table 1). For each CXC sample, PN was soaked in 50% ethanol twice. The soaking liquids were filtered separately, and the two filtrates were mixed together. The ethanol was recovered, and the water was evaporated to produce the PN extraction. Then RA, SM, and RS were mixed together and refluxed twice with 50% ethanol. The reflux extraction liquid was filtered separately, and both filtrates were mixed together. The ethanol was recovered, and the water was evaporated to produce the three-herb extraction, which was combined with the PN extraction to form the CXC sample [13]. Each CXC sample was made into the same dose of 152 mg/ml by diluting with normal saline. The 4 herbs (batch No: 120523) were provided by Guangdong Zhongsheng Pharmaceutical Co., Ltd. (Dongguan, China).

HPLC and cluster analysis of CXC samples
Liang et al have previously established the standard HPLC fingerprint of CXC and identified 21 components [14]. The HPLC fingerprints of nine CXC samples were obtained using the established method, and a Dionex Acclaim 120 C 18 column (3 mm, 150 mm64.6 mm) was used with acetonitrile (A) and 0.05% orthophosphoric acid (B) in gradient elution mode. The elution profile was 0-50 min (15%R34%A), 50-95 min (34%R75%A); the detection wavelengths were 203 nm and 270 nm; the column temperature was set at 25uC, and the flow rate was 1.0 ml/min. Peak areas were calculated with a Chromeleon 6.8 Chromatography Data System (Dionex). The peak areas of 21 identified components were adopted as independent component variables ( Table 2). Cluster analysis of the nine CXC samples was conducted in SPSS 18.0. The clustering method was between-groups linkage, and the distance calculating method was Pearson's correlation.

Experimental animal model
Rats were divided randomly into 13 groups: control, model, Asp (100 mg/kg/d), CDDP (800 mg/kg/d), and 1-9 CXC samples (1,520 mg/kg/d), with 10 rats in each group. Rats received treatment (10 ml/kg/d) via gavage once daily for 10 consecutive days. The control and model groups received the same volume of normal saline. Thirty minutes after the tenth administration, all rats (except those in the control group) were subcutaneously injected with Adr (0.8 mg/kg), while the control group was subcutaneously injected with the same volume of normal saline. Two hours later, all the rats (except those in the control group) were kept in ice-cold water (0-4uC) for 5 minutes. Two hours later, they received a second subcutaneous injection of Adr (0.8 mg/kg) [15]. After the re-injection, all the rats were fasted for 12 hours and then given the last drug treatment.

Blood sample collection
Chloral hydrate (10%, 3.5 ml/kg) was administered via intraperitoneal injection for anesthesia. Blood was drawn from the abdominal aorta minutes after the last administration and collected into plastic tubes with 3.8% sodium citrate (citrate/ blood: 1/9, v/v). All the blood samples were strictly processed and examined using standard operating procedures.

Hemorheology and blood clotting activity assessments
Blood plasma obtained from 1.5 ml blood (3,820 r/min, 20uC, 15 min) was put into a Sysmex CA-510 fully automatic blood coagulation analyzer to determine activated partial thromboplastin time (APTT), prothrombin time (PT), and plasma viscosity (PV). Blood (0.9 ml) was put directly into a fully automatic selfcleaning hemorheology analyzer (LBY-N6B) to measure the erythrocyte aggregation index (EAI) and the red corpuscle electrophoresis index (RCEI), as well as whole-blood viscosity (WBV) at 5/s, 30/s, 50/s, 150/s, and 200/s shear rates. Blood (0.9 ml) was centrifuged (3,820 r/min, 20uC, 15 min) for hematocrit detection by an automatic dynamic blood sedimentation tester (LBY-XC40). Platelet-rich plasma (PRP) was separated from 3.0 ml blood sample (500 r/min, 20uC, 10 min), then the remaining without PRP was used to obtain platelet-poor plasma (PPP) (3,000 r/min, 20uC, another 10 min). Platelet counts were kept approximate 200610 9 /L in each PRP sample. Platelet aggregation was induced by 5 ml adenosine diphosphate (ADP). PRP (300 ml) and PPP (300 ml) were used to detect maximum platelet aggregation rate (MPAR) on a platelet aggregation analyzer (LBY-NJ4). All the tests were completed in the drug non-clinical evaluation research center of Guangzhou General Pharmaceutical Research Institute. 8 The relevance analysis between HPLC fingerprint and drug effects: mathematical approaches The peak areas of 21 identified components were regarded as 21 independent component variables (p1-p21). For each component, the corresponding peak area was divided by its area average in 9 CXC samples for nondimensionalization (Table 3). For each measurement index, the effect of each CXC sample was represented by the average value in the corresponding group. The negative indices (the larger the value, the smaller the effect) were switched to positive by taking the reciprocal. Then, the average value in each group was divided by the corresponding average in 9 CXC sample groups for nondimensionalization. The dimensionless data are shown in Table 4.
Factor analysis (FA): five factors with no mutual correlation were extracted from 11 intercorrelated indices. The relevance between 11 indices and 5 factors was calculated with a rotation matrix method to explain the clinical meaning of the factors ( Table 5). The scores of five factors were then calculated to be the new dependent variables F1-F5 ( Table 6). The processes above were implemented in SPSS 18.0 (Factor Analysis, Data Reduction, Analyze).
In grey relational analysis (GRA), the grey relational degree (GRD) is used to rank the influence of compared series, which can be represented by the relative distance between the compared series and reference series in an imaging grey space without making any prior assumptions about the distribution type [16]. In the experiment, five factors (F1-F5) were chosen as 5 reference series, and 21 component variables (p1-p21) were chosen as 21 compared series. Then, the GRD between the compared and reference series was calculated with a resolution ratio of 0.5 [16]. The higher the GRD, the greater the effect of the corresponding component.
For principal component analysis (PCA), 21 intercorrelated original independent component variables were recombined into two mutual independent principal components (Table 7), which were regarded as new independent variables. Then, the regression equations between the two principal components and the factors were constructed by a stepwise regression analysis approach. Once a strict regression equation was established (P,0.05), two principal components would be replaced by the 21 original component variables (Table 7)   In ridge regression-radial basis function (RR-RBF) analysis, every ridge parameter K corresponded to a standard ridge regression coefficient (SRRC) for each independent component variable. The ridge trace of each independent component variable, a curve of SSRC on different K (0.02-1, increasing progressively by 0.02), was drawn on a coordinate system. Next, the variables with stable but small absolute values of SRRC and the variables with unstable SRRCs verging infinitely to zero were removed [17]. After that, the remaining variables were used to draw a new ridge trace (Fig. 2), in which the minimum K stabilizing all the ridge traces of independent component variables could be identified. Finally, with this minimum K (0.18), the proper SRRCs of the remaining variables were obtained to evaluate the contribution of each component. Meanwhile, the RBF network between the same remaining variables and each factor (F1-F5) was constructed. Then, the bond value (BV) and standard importance of variable (SIV) of the remaining variables were calculated. BV was used to evaluate the contributions of the components, and SIV represented the reliability of independent component variables in the established RBF network. All the processes above were implemented in SPSS 18.0 (Ridge, Catreg, Regression, Analyze; RBF, Neural Networks, Analyze).

Statistical analysis
The data are presented as mean 6 SD and were evaluated by one-way analysis of variance (ANOVA), Student's t-tests, and Dunnett's multiple comparisons in SPSS 18.0 (SPSS Inc., Chicago, IL, USA). P-values less than 0.05 or 0.01 were considered statistically significant.

CXC sample generation
UD was first proposed by Fang in 1978 and is based on either a quasi-Monte Carlo method or a number-theoretic method. UD is capable of producing samples with high representativeness in the studied experimental domain [18]. Most notably, four-factor, nine-level UD was adopted to guarantee the component difference within 9 CXC samples while reducing workload.

HPLC and cluster analyses of CXC samples
The HPLC fingerprint of CXC was established in 2012 by our team; we detected 42 common peaks covering 4 herbs, and 21 unique components were identified by RRLC/MS/MS and diode array detector (DAD) spectra (Table 2) [14]. In this experiment, the HPLC fingerprint of 9 CXC samples was constructed as shown in Fig. 3. According to the peak areas of 21 identified components, 9 CXC samples could be divided into 7 categories: S2 and S3 belonged to a class, S7 and S8 belonged to a class, and the rest of the samples each represented a class.

CXC sample effects in vivo
As an acute stressor, Adr injections can cause various hemorheological disorders, including blood hypercoagulability and increased blood clotting activity [19]. In addition, skin blood flow reduces rapidly when humans are exposed to ice-cold water [20]. It was previously reported that the combination of Adr injection and exposure to ice-cold water could induce poor blood circulation [21]. Therefore, in this experiment, the rats were placed in ice-cold water during the time interval between Adr injections being administered to generate an appropriate model. Table 7. The scores of two principal components (extracted from 21 component variables). Rats with blood stasis exhibit sluggish circulation, leading to changes in hemorrheologic parameters [12]. In recent years, hemorheology has commonly been used to diagnose and prevent diseases and to evaluate the therapeutic effects of drugs on blood circulation [22]. WBV reflects the intrinsic resistance of blood flow in the vasculature, and abnormal WBV is associated with increased cardiovascular risk. The reduction of blood fluidity can lead to tissue ischemia more rapidly in patients with atherosclerotic diseases [23]. WBV at a low shear rate reflects red blood cell (RBC) aggregation and WBV at a high shear rate reflects RBC deformability [24]. WBV is determined not only by RBCs but also by plasma viscosity (PV), which depends on the type and concentration of plasma proteins. Both EAI and RCEI reflect the extent of aggregation among RBCs. Blood platelet aggregation induced by ADP plays a major role in the development and extension of arterial thrombosis [25]. APTT is a test of the intrinsic clotting activity and PT is used to evaluate the overall efficiency of the extrinsic clotting pathway [26]. A prolonged APPT or PT indicates a deficiency in coagulation factors or the presence of inhibitors of coagulation [27]. Fig. 4 and Fig. 5 show the different effects observed in all 13 treatment groups, with particularly noticeable differences in the 9 CXC groups. There were significant differences between the control rats and model rats for all 11 indices. After treatment, the results showed that indices differed greatly among the 9 CXC samples, and we plan to assess possible reasons for this in subsequent analyses. In addition, the original data that were used to generate Fig. 4 and Fig. 5 were uploaded as Table S1. 4 Five factors (F1-F5) extracted from the original data FA simplifies multivariate data consisting of a large number of intercorrelated variables by grouping them into a smaller set of independent factors or clusters according to the basic underlying relationships between them. Factors are fewer in number than the original variables, account for a significant proportion of data variance, and are useful in predictive regression models [28]. The results in Table 5 show that five mutually independent factors reflected 98% of the original data. Factors 1-5 represented five aspects of the vascular system with clear clinical significance: RBC aggregation (F1), platelet aggregation (F2), intrinsic clotting activity (F3), RBC deformability and plasma proteins (F4), and extrinsic clotting activity (F5).

Selection and combination of suitable mathematical approaches
Every method has a different theoretical basis and suitable scope, hence a combination of different methods is necessary for a comprehensive and accurate interpretation of the relationship between HPLC data and experimental data. In fact, regarding 21 components from 4 herbs as independent variables probably led to collinearity. Therefore, PCA, which is often used for data reduction by structuring many variables into a much smaller number of components to eliminate variable collinearity, was chosen for this study. Two principal components contributing to 96.67% of the total variance were extracted, and the scores between them and the 21 component variables were calculated for the transformation (Table 6). Ngo et al. demonstrated that RR is clearly superior in that it provides estimated parameters much closer to the true values than those obtained by the ordinary least squares regression method [29]. The RBF network, a powerful type of artificial neural network for forecasting and modeling a multivariable system, is one of the artificial learning methods that describe complex linear or non-linear relationships between inputs and outputs [30]. GRA is suitable for solving problems with  Table 2. C shows the cluster analysis result of 9 CXC samples. doi:10.1371/journal.pone.0112675.g003 complicated interrelationships between multiple factors and variables and has been demonstrated as a useful approach to solve multiple-attribute decision-making problems [31]. Generally speaking, one of the main purposes of relation analysis is to identify the importance order of factors, so the GRD rankings are more significant [32]. Overall, GRA, PCA, RR, and RBF network were combined to control the mutual interference among the 21 component variables to better assess the data.

The core bioactive components promoting blood circulation
Under the hypothesis that the effects of CXC samples varied with component ratio, the relationship between component and effect differences was analyzed to evaluate the contribution of each component using the methods of GRA, PCA, and RR-RBF. The computed results are shown in Table 8. In addition, a dynamic component-effect analysis bubble chart based on the analysis results is provided as Component-effect relevance bubble chart S1 to provide a clear picture of the comprehensive effects of all 21 components.
F1 represents RBC aggregation in the vascular system. The results show that panaxytriol and ginsenoside Rb 1 were the major effective components to suppress RBC aggregation. Ginsenoside Rb 1 , with high PN saponin content, can ameliorate lipopolysaccharide (LPS)-induced microcirculatory disturbance in the rat mesentery [33]. The ability of ginsenoside Rb 1 to suppress RBC aggregation, revealed in this study, might be the foundation of its effects on blood circulatory disturbances, as evidenced by increasing studies. Panaxytriol has a relatively low content in PN and has only been reported to exhibit potent anti-tumor properties [34]. Until now, there have been no descriptions of the pharmacological effects of panaxytriol against vascular dysfunction. The results of this study suggest that it could be useful to further investigate panaxytriol's novel therapeutic effects.
F2 represents platelet aggregation in the vascular system. Angoroside C was more closely related to platelet aggregation than the other components. Angoroside C is one of the main compounds involved in the anti-inflammatory effect of RS and has a significant effect on thromboxane B 2 (TXB 2 ) release and nitric oxide (NO) activity in LPS-induced macrophages [35]. Thromboxane A 2 (TXA 2 ) has been demonstrated to induce platelet aggregation, and TXB 2 is the stable hydrolysate of TXA 2 . Therefore, angoroside C's significant suppression of platelet aggregation might be partly attributed to the effect on TXB 2 release, which ameliorates blood hypercoagulability.
F3 represents intrinsic clotting activity in the vascular system. All 21 components exhibited relatively low relevance to F3, suggesting that the contribution of every single component to F3 was small. Nevertheless, the results showed that protocatechualdehyde had the closest relevance to F3 and might be the bioactive component. Protocatechualdehyde is known for its antiproliferative and antioxidant activities, as well as its ability to inhibit aldose WBV at all shear rates, EAI, and PV increased significantly and RCEI decreased significantly in model rats. After treatment, Asp was significantly effective in decreasing WBV, EAI, and MPAR. CDDP was significantly effective in reducing WBV and EAI. For CXC samples, WBV decreased significantly at 5/s, 30/s, 50/s and 150/s in S1-S5 group. EAI decreased significantly in all groups except for S7 and S9 group. RCEI increased significantly in S2, S4, and S5 group. PV decreased in S1, S3, and S5 group. * P,0.05 and ** P,0.01 vs control group, # P,0.05 and ## P,0.01 vs model group, n = 10. doi:10.1371/journal.pone.0112675.g004 reductase in rat lens epithelial cells [36]. This study indicated that protocatechualdehyde might influence intrinsic clotting activity and could be of therapeutic interest with respect to promoting blood circulation.
F4 represents RBC deformability and plasma proteins in the vascular system, which is closely related to blood microcirculatory disturbance. In this study, ginsenoside Rd had a significantly higher influence on F4 than the other components. More than 20 basic research papers on ginsenoside Rd have been recently published in many journals. In these studies, ginsenoside Rd was reported to have many different effects as a bioactive component: it could attenuate early oxidative damage and sequential inflammatory responses after transient focal ischemia in rats [37], relieve redox imbalance and improve stroke outcome after focal cerebral ischemia in aged mice [38], and prevent the development of atherosclerosis by blocking Ca 2+ influx through receptor-and store-operated Ca 2+ channels in vascular smooth muscle cells [39,40]. The present study provides evidence for ginsenoside Rd's possible regulation of RBC deformability and plasma proteins, yet elucidation of the underlying mechanism requires further research.
F5 represents extrinsic clotting activity in the vascular system. We found that calycosin-7-O-b-D-glucoside was most closely related to F5. It has been reported that calycosin-7-O-b-Dglucoside can increase endogenous antioxidant levels [41] and improve cell membrane fluidity in the brains of rats subjected to experimental stroke [42]. However, there is still no direct evidence showing calycosin-7-O-b-D-glucoside's activity on the extrinsic Figure 5. The animal experiment results of CXC samples 1-9. PT and APTT decreased significantly and MPAR increased significantly in model rats. After treatment, APTT was significantly influenced in S5 group. MPAR decreased significantly in S1 and S5 group. PT was prolonged in S3, S5, and S6 group. * P,0.05 and ** P,0.01 vs control group, # P,0.05 and ## P,0.01 vs model group, n = 10. doi:10.1371/journal.pone.0112675.g005 For GRD, RC, PCC, SRRC and BV, positive and negative of value means positive and negative effect contribution respectively. The larger the absolute value, the higher the effect contribution. SIV represents the relative importance and reflects the reliability of corresponding variables. clotting pathway. The present study indicates a new potential pharmacological activity of calycosin-7-O-b-D-glucoside. We determined that panaxytriol, ginsenoside Rb 1 , angoroside C, protocatechualdehyde, ginsenoside Rd, and calycosin-7-O-b-D-glucoside were the core bioactive components of CXC and that they acted on different aspects of the vascular system. Panaxytriol and ginsenoside Rb 1 were relevant to RBC aggregation, angoroside C to platelet aggregation, protocatechualdehyde to intrinsic clotting activity, ginsenoside Rd to RBC deformability and plasma proteins, and calycosin-7-O-b-D-glucoside to extrinsic clotting activity. Moreover, this study provides an important reference for the development of novel therapeutics on promoting blood circulation that employ angoroside C, calycosin-7-O-b-D-glucoside, panaxytriol, and protocatechualdehyde.
Besides, according to the attribution of 6 bioactive components, PN might affect two aspects of the vascular system, including F1 (RBC aggregation) and F4 (RBC deformability and plasma proteins), while RA, RS, and SM might influence the F5 (extrinsic coagulation system), F2 (platelet aggregation), and F3 (intrinsic coagulation activity), respectively. Consequently, this work also partly revealed the mutual promotion among 4 herbs and demonstrated the multicomponent and multitarget properties of CXC.
However, the results also showed that rosmarinic acid, ononin, calycosin-7-O-b-D-glucoside, and panaxytriol might have negative effects on RBC aggregation, RBC deformability, intrinsic clotting activity, and platelet aggregation (Table 8). Both CHFs and organisms are very complex systems and the specific interactions between them are still unclear; therefore, further exploration and verification are needed to understand the pharmacological mechanisms of the 4 components listed above. Nevertheless, this study demonstrated that one specific aspect of the vascular system might be affected by two or more CXC components, and the same component might have a positive effect on one aspect yet a negative effect on another, which reflects both the synergistic and the antagonistic actions of CHF components.
Comprehensive qualitative or quantitative detection based on the low or high content of core bioactive components is necessary to improve CXC quality control. Therefore, we established the bioactive HPLC fingerprint of CXC based on its core bioactive components (Fig. 6). In summary, we explored the core bioactive components of CXC with a series of HPLC and animal experiments. The results suggested that the core bioactive components of CXC promoting blood circulation were panaxytriol, ginsenoside Rb 1 , angoroside C, protocatechualdehyde, ginsenoside Rd, and calycosin-7-O-b-D-glucoside. Overall, these findings provide a novel perspective to pinpoint the core bioactive components in a CHF, which will be helpful for improving quality control and inspiring further clinical studies of CHFs. It might also provide guidance for drug discovery. Table S1 The results of the animal experiments in 13 treatment groups. Control group and model group received the same volume of normal saline (NS) for the treatment. * P,0.05 and ** P,0.01 vs control group, # P,0.05 and ## P,0.