Deciphering the Active Compounds and Mechanisms of Qixuehe Capsule on Qi Stagnation and Blood Stasis Syndrome: A Network Pharmacology Study

Background Qixuehe capsule (QXH), a Chinese patent medicine, has been demonstrated to be effective in the treatment of menstrual disorders. In traditional Chinese medicine (TCM) theory, qi stagnation and blood stasis syndrome (QS-BSS) is the main syndrome type of menstrual disorders. However, the pharmacodynamic effect of QXH in treating QS-BSS is not clear, and the main active compounds and underlying mechanisms remain unknown. Methods A rat model of QS-BSS was established to evaluate the pharmacodynamic effect of QXH. Thereafter, a network pharmacology approach was performed to decipher the active compounds and underlying mechanisms of QXH. Results QXH could significantly reduce the rising whole blood viscosity (WBV) and plasma viscosity (PV) but also normalize prothrombin time (PT), activated partial thromboplastin time (APTT), thrombin time (TT), and fibrinogen (FIB) content in QS-BSS rats. Based on partial least-squares-discriminant analysis (PLS-DA), the low-dose QXH-intervened (QXH-L) and the high-dose QXH-intervened (QXH-H) groups seemed the most effective by calculating the relative distance to normality. Through network pharmacology, QXH may improve hemorheological abnormality mainly via 185 compounds-51 targets-28 pathways, whereas 184 compounds-68 targets-28 pathways were associated with QXH in improving coagulopathy. Subsequently, 25 active compounds of QXH were verified by UPLC-Q/TOF-MS. Furthermore, 174 active compounds of QXH were shared in improving hemorheological abnormality and coagulopathy in QS-BSS, each of which can act on multiple targets to be mainly involved in complement and coagulation cascades, leukocyte transendothelial migration, PPAR signaling pathway, VEGF signaling pathway, and arachidonic acid metabolism. The attribution of active compounds indicated that Angelicae Sinensis Radix (DG), Paeoniae Radix Rubra (CS), Carthami Flos (HH), Persicae Semen (TR), and Corydalis Rhizoma (YHS) were the vital herbs of QXH in treating QS-BSS. Conclusion QXH can improve the hemorheology abnormality and coagulopathy of QS-BSS, which may result from the synergy of multiple compounds, targets, and pathways.


Introduction
Qixuehe capsule (QXH) is a Chinese patent medicine and used clinically for the treatment of menstrual disorders for a dozen years. Menstrual disorders are known to be a common gynecological disease, which is characterized by abnormal menstrual cycle, dysmenorrhea, amenorrhea, uterine bleeding, postpartum abdominal pain, lochia, breast swelling, infertility, and chloasma, seriously affecting female's health [1,2]. Although the symptoms of menstrual disorders are complex and diverse, qi stagnation and blood stasis syndrome (QS-BSS) is the main syndrome type of menstrual disorders in traditional Chinese medicine (TCM) theory. QS-BSS is defined as due to the poor stagnation of the operation of the qi, the blood is running impediment and the pathological state of blood stasis occurs, which is mainly related to menstrual disorders, hyperlipidemia, liver injury, ischemic brain injury, and hypertension [3][4][5].
e formation of QS-BSS is accompanied by the following pathological changes: (1) abnormal hemorheology; (2) increased blood coagulability or reduced fibrinolytic activity; (3) microcirculation disturbance; (4) increased platelet aggregation or decreased release function; (5) hemodynamic disturbance; (6) manifestations of static blood in pathological section; (7) vessel obstruction [6]. Numerous studies have shown that the pathological changes of hemorheological abnormality and coagulopathy are more important during the formation of QS-BSS [7,8]. TCM used to treat syndrome by compound formulas, which are composed of several or even tens of herbal medicines in different quantities, following the composition theory "Monarch, Minister, Assistant and Guide" [9]. QXH is based on the addition and subtraction of Xuefu Zhuyu decoction (a classic blood-activating formula) and Chaihu Shugan powder (a classic qi-regulating formula), which consists of Angelicae Sinensis Radix (DG), Paeoniae Radix Rubra (CS), Chuanxiong Rhizoma (CX), Persicae Semen (TR), Bupleuri Radix (CH), Carthami Flos (HH), Salviae Miltiorrhizae Radix et Rhizoma (DS), etc., so as to regulate liver-qi stagnation, promote blood circulation, and remove blood stasis in the treatment of QS-BSS [10,11]. However, the potential active compounds and underlying mechanisms of QXH in treating QS-BSS remain unclear.
In recent years, network pharmacology has been proposed as a promising approach for understanding compound formulas [12,13]. Studies have shown that Xuefu Zhuyu decoction and compound Danshen formula have been analyzed by network pharmacology, and the main active compounds in the formulas and the treatment mechanisms have been elucidated [14,15]. erefore, network pharmacology is promising to decipher the active compounds and possible mechanisms of QXH.
In this study, we evaluated the pharmacodynamic effect of QXH in the treatment of QS-BSS based on animal experiments from the two aspects: hemorheology and coagulation function. Hemorheological abnormality is reflected in changes in blood fluidity, deformability, and viscosity, which can be measured by whole blood viscosity (WBV) and plasma viscosity (PV). Meanwhile, coagulopathy is reflected in the abnormality of coagulation factor function; thrombin time (TT), prothrombin time (PT), activated partial thromboplastin time (APTT), and fibrinogen content (FIB) are the important monitoring index of coagulopathy in clinical aspects [16]. Subsequently, a network pharmacology approach was performed to decipher the active compounds of QXH and its regulation mechanisms of QS-BSS. By constructing two herb-compound-pathological process relationship diagrams, the possible mechanisms of QXH in the treatment of QS-BSS were systematically elucidated. A workflow of the entire study procedure is depicted in Figure 1.

Animal
Handling. Female Sprague-Dawley rats (200 ± 20 g) were provided by Dossy Experimental Animals Co., Ltd., Chengdu, China. All procedures were conducted in accordance with the Regulations of Experimental Animal Administration issued by the State Committee of Science and Technology of the People's Republic of China and approved by the Animal Ethics Committee of Shaanxi University of Chinese Medicine.

Sample Collection and Indicator Measurements.
On the 14 th day, after subcutaneous injection of Adr, the rats were managed to be fasting but had free access to clean water. e next morning at 8 : 00, rats were anesthetized with 10% chloral hydrate (350 mg/kg i.p.). ereafter, 4 mL blood samples of rats were collected by the abdominal aorta intubation approach into a 5 mL disposable human venous blood sample collection container (heparin sodium anticoagulation).

QXH Database Construction.
e chemical compounds of 15 herbs of QXH were data-mined from relevant databases including Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform (TCMSP, http://lsp.nwu.edu.cn/index.php) [19] and TCM-MESH (http://mesh.tcm.microbioinformatics.org), and all candidate compounds were also related to the literature to determine the attribution relationship with herbs. Latin names, the number of active compounds of each herb in QXH involving in mediating hemorheological abnormality and coagulopathy, and the corresponding abbreviations of each herb are displayed in Table 1 Figure 1: e workflow of QXH in treating QS-BSS combined with pharmacodynamic evaluation and network pharmacology prediction.
Evidence-Based Complementary and Alternative Medicine 3 contribute to its therapeutic effects. In order to obtain the potential active compounds from these herbs of QXH, we applied a method incorporating oral bioavailability (OB), Caco-cell permeability (Caco-2), and drug likeness (DL) evaluation in this work [20,21]. OB is an important pharmacokinetic parameter which describes the ratio of the amount absorbed into systemic circulation when dosed orally compared to the amount in blood upon intravenous dosing. In the present study, the OB values were evaluated using the robust in silico tool OBioa-vail1.1 [22]. e compounds with OB ≥ 30% were assessed as active compounds for further analysis. Caco-2 is another vital parameter frequently used as a model for studying the passive diffusion of drugs across the intestinal epithelial barrier. us, the in silico Caco-2 permeability model which was constructed by 100 drug molecules with satisfactory statistical results (R 2 > 0.8) was employed to select compounds that are more likely to possess good permeability [23]. Given those compounds with Caco-2 value less than − 0.4 not permeable, the threshold of Caco-2 permeability is set at − 0.4.
DL is a qualitative property of chemicals that describes how the pharmacokinetic and pharmaceutical properties of compounds, such as solubility and potency, correspond in the majority of known drugs. In order to identify drug-like compounds, we applied a database-dependent model using the Tanimoto coefficient [24] to calculate the DL (see equation (1)) of each compound in QXH: In equation (1), A represents the molecular descriptors of herbal compounds and B displays the average molecular properties of all compounds in Drugbank.
ose compounds with DL ≥ 0.18 were preserved. e compounds were adopted as the candidate compounds for further analysis when they met all these three criteria.
Besides, due to the profound pharmacological effects, biological activity, and high contents, those compounds with low OB, Caco-2, or DL values were also selected for further analysis (Table S1). Last but not least, all active compounds were imported into PubChem (https://pubchem.ncbi.nlm. nih.gov) [25] and ChemSpider (http://www.chemspider. com) [26] to standardize their English names.

Sample Preparation and UPLC-Q/TOF-MS Verification.
Reference standards were accurately weighed and directly prepared in methanol as individual stock solutions; subsequently, the stock solutions of twelve standards were continuously diluted with methanol to obtain a suitable concentration for further analysis, respectively. QXH samples were divided into two. 1.0000 g QXH powder was accurately weighed for each and extracted with methanol (25 mL) and water (25 mL) in ultrasonic water bath (40 kHz, 500 W) for 30 min at room temperature, respectively. Subsequently, the QXH methanol and water extracts were centrifuged at 14,000 rpm for 10 min, and the supernatant was stored at 4°C until analysis.   e source temperature and desolvation gas temperature were set to 100 and 350°C, respectively, and the cone gas flow and desolvation gas flow were set to 50 L/h and 600 L/h, respectively. e collision energy was set to 4 eV for positive ion mode and 2 eV for negative ion mode. Sodium formate solution was used to calibrate the mass spectrometer prior to the experiment. Leucine-enkephalin was used as an external reference (Lock-Spray ™ ) to correct the mass during data

Target Prediction.
In order to identify the corresponding targets of active compounds in QXH, several methods combined with chemometric strategy, database dependent, information integration, and literature-mining were implemented [27]. Firstly, all active compounds were imported into TCMSP and TCM-MESH again to obtain relevant targets. Subsequently, all active compounds were imported into STITCH (http://stitch.embl.de/) [28], erapeutic Targets Database (TTD, http://bidd.nus.edu.sg/ group/ttd/) [29], Comparative Toxicogenomics Database (CTD, http://ctdbase.org/) [30], and Google Scholar to obtain H-C-T. In view of the preliminary mining results, the active compounds with less targets were introduced into the Similarity Ensemble Approach (SEA, http://sea.bkslab.org/) [31] and Swiss Target Prediction (http://www. swisstargetprediction.ch) [32] to predict their possible targets. Given the prediction results with Max TC < 0.57 and Possibility < 0 not considered as probably having affinity to the appropriate target, the targets with Max TC > 0.57 in SEA and Possibility > 0 in Swiss Target Prediction were kept in the result, respectively. ereafter, for the purpose of dissecting the mechanisms of the candidate targets in regulating the pathological processes of QS-BSS, the target information was submitted to TTD, CTD, Online Mendelian Inheritance in Man Database (OMIM, http://www.omim. org/) [33], and Uniprot (https://www.uniprot.org) [34] to determine the targets which were related to hemorheological abnormality and coagulopathy, respectively. Finally, all the targets from the previous mining were imported into Uniprot to obtain their corresponding gene symbol, and only the targets of Homo sapiens were kept for further analysis.

Gene Ontology Enrichment and KEGG Pathway Analysis.
e Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov) web server was employed to perform gene ontology enrichment and KEGG pathway analysis for the targets related to hemorheological abnormality and coagulopathy, respectively [35]. Meanwhile, gene ontology enrichment and KEGG pathway analysis were also performed for the overlapped targets involving in the pathological processes of hemorheological abnormality and coagulopathy.

Network Construction.
e herb-active compound-target (H H -C H -T H ) and herb-active compound-target (H C -C C -T C ) networks of QXH were constructed in the hemorheological abnormality module and coagulopathy module, respectively. All networks were generated and analyzed by Cytoscape 3.6.0 (http://www.cytoscape.org/), an open-source software package project for visualizing, integrating, modeling, and analyzing the interaction networks [36].

Data Processing and Statistical Analysis. All LC-MS and
MS/MS data were processed with MassLynx ™ (V4.1, Waters, Manchester, UK) software. Molecular formula speculations of the compounds were determined with Elemental Composition software. en, the acquired raw data matrices on hemorheology and coagulation function were imported to SIMCA-P (version 14.1, Umetrics, Sweden) for principal component analysis (PCA) and partial least-squares-discriminant analysis (PLS-DA). PCA is an unsupervised multivariate statistical approach used for variable reduction and separation into classes. To maximize class discrimination, the data were further analyzed utilizing PLS-DA [16,37]. e parameters of PCA (R 2 X and Q 2 ) and PLS-DA (R 2 X, R 2 Y and Q 2 ) were used to evaluate the model, and Q 2 is an appraisal criterion of the predictive ability of a model. R 2 X represents the capability wherein the primary components of variables are used to build the model and sample, whereas R 2 Y represents the conformity level of the sample with the model. en, the PLS-DA model was validated by 200-repeated permutation tests and CV-ANOVA (analysis of variance testing of cross-validated predictive residuals) tests. A model with the p value (obtained by CV-ANOVA) below 0.5 or with R 2 and Q 2 values (obtained by permutation tests) lower than the values of the model itself was considered to be valid. e value of variable importance in the projection (VIP) from the PLS-DA model represents the contribution of variables to the classification of the model which was analyzed. e differential parameters between the control and the model groups were selected by VIP >1. All data were expressed as mean ± standard deviation (SD) and analyzed by one-way ANOVA in SPSS 16.0 (SPSS Inc., Chicago, IL, USA). e criterion for statistical significance was p < 0.05.

Effects of QXH on Hemorheology in Rats with QS-BSS.
QS-BSS mainly occurs with increased blood viscosity and poor blood fluidity. WBV and PV are the vital hemorheological parameters and commonly used to evaluate the clinical therapeutic efficacy of drugs. In Table 2, the levels of WBV in low and high shear rates (1 s − 1 , 5 s − 1 , 50 s − 1 , 100 s − 1 , Evidence-Based Complementary and Alternative Medicine and 200 s − 1 ) and PV were significantly elevated in the model group in comparison with the control group (p < 0.01), indicating a successful modeling. Compared with the model group, the WBV at the shear rate of 1 s − 1 , 5 s − 1 , 50 s − 1 , 100 s − 1 , and 200 s − 1 in QXH and CDDP groups significantly decreased near to normality (p < 0.05 or p < 0.01). Meanwhile, QXH and CDDP groups could significantly decrease the rising PV in QS-BSS rats (p < 0.01), suggesting that both QXH and CDDP treatment could significantly improve the hemorheological abnormality in rats with QS-BSS, and QXH is no less effective than CDDP in improving hemorheology in QS-BSS rats.

Effects of QXH on Coagulation Function in Rats with QS-BSS.
APTT, PT, TT, and FIB, as coagulation function test items in clinics, are crucial quantitative indexes of pathological changes in the coagulation system of the body [36]. As is illustrated in Table 3, compared with the control group, TT, PT, and APTT levels were significantly shortened (p < 0.05 or p < 0.01) and FIB content was significantly elevated (p < 0.01) in the model group. Compared with the model group, the level of TT, PT, APTT and FIB in QXH groups had different degrees of callback. Specifically, QXH-L treatment could significantly rectify the perturbed TT and PT after QS-BSS modeling (p < 0.05); the QXH-M group significantly improved the abnormality of TT and APTT (p < 0.05 or p < 0.01); the QXH-H group showed a significantly effect in improving TT (p < 0.01). CDDP significantly brought the disrupted TT, APTT, and FIB after QS-BSS modeling to near normality (p < 0.05 or p < 0.01). All above results fully prove that QXH is comparable with CDDP on improving coagulation function in QS-BSS rats.

Evaluation of the Effects of QXH on Improving Hemorheological Abnormality and Coagulopathy Based on Multivariate Statistical Analysis.
In this study, PCA was used firstly to determine the general interrelation between control, model, QXH-L, QXH-M, QXH-H, and CDDP groups. As is shown in Figure S1, the model group and other groups were separated clearly, indicating that the rats of the model group have significantly altered on hemorheology and coagulation function.
e parameters of the PCA model: Figure S1) suggesting good ability of prediction and reliability of the model, and the variances of Comp [1] and Comp [2]  Subsequently, the second round of the PLS-DA method was used to analyze the VIP value of each index between control and model groups to obtain the most contribution indexes to the two groups' differences. e parameters of PLS-DA between control and model groups:  Figure 2(b), the value of VIP from the PLS-DA model represents the contribution of variables of the indexes to the classification between control and model groups. e larger the VIP value, the greater the contribution of those indexes to the difference between control and model groups. e VIP value of WBV5 s − 1 , WBV1 s − 1 , WBV200 s − 1 , FIB, and WBV50 s − 1 are all >1, indicating WBV and FIB are vital to distinguish the difference between normal and QS-BSS levels and can be used as the essential indexes to evaluate the pharmacodynamic effect of QXH in the treatment of QS-BSS. e VIP value of TT and APTT are about 0.9, suggesting the poor contribution to the difference between control and model groups.
Single-nucleotide polymorphisms in the ESR1 gene have been associated with vascular diseases and conditions, such as arterial hypertension, CVD, cerebrovascular disease, and serum lipid level alterations [43]. In the coagulopathy module, besides PTGS1, PTGS2, and PIK3CG, F2, F10, PRKACA, NOS3, PTPN2, NQO1, and HMOX1 were also included. F2 and F10 cleave bonds after Arg and Lys, convert fibrinogen to fibrin and activate factors V, VII, VIII, and XIII, in combination with thrombomodulin and protein C, and has function in blood homeostasis, inflammation, and wound healing [44]. NOS3 promotes the production of NO; NO can mediate vascular endothelial growth factor-(VEGF-) induced angiogenesis in coronary vessels and promote blood clotting through the activation of platelets [45].

H H -C H -T H Network of QXH (Hemorheological Abnormality Module).
Hemorheological abnormality is an important part of the formation of QS-BSS, which is characterized by changes in blood fluidity, viscosity, deformability, and coagulability. us, we constructed a H H -C H -T H network through network analysis to illuminate the relationship among herbs, active compounds, and candidate targets (Figure 3 (for original picture, see Figure S7)). is network consists of 251 nodes (15 herbs, 185 active compounds, and 51 candidate targets; Tables S2 and S3). As is depicted in the H H -C H -T H network, red, green, and blue nodes represent herbs, active compounds, and candidate targets, respectively. e size of the nodes is positively related to the value of degree. Degree, which is the number of edges connected to the node, and the betweenness, another centrality index defined by the number of times a node acts a bridge along the shortest path between two other nodes, was calculated. Actually, betweenness reflects the fraction of the shortest paths in the network that pass through any particular node and a measure of the importance of a node as a hub in a network [46]. To gain insights into the pharmacological mechanisms of QXH on hemorheological abnormality, we performed gene ontology enrichment and KEGG pathway analysis on the 51 relevant targets, and the top 10 of biological process, molecular function, and cellular component and 28 KEGG pathways were predicted ( Figure 4). e results of biological process were mainly enriched in regulation of cholesterol storage, lipoprotein metabolic process, blood circulation, and regulation of lipid storage, which were consistent with the pathologic mechanisms of hemorheological abnormality. e degree and betweenness of the top 10 candidate compounds and targets involvied in hemorheological abnormality module are shown in Tables 4 and 5. Among the active compounds, quercetin, luteolin, isorhamnetin, butylated hydroxytoluene, and isoliquiritigenin were the top ones, interacting with 30, 13, 13, 12, and 12 targets, respectively. Quercetin acted on PTGS1, PTGS2, ALOX15, and PPARA and participated in the PPAR signaling pathway, prolactin signaling pathway, and VEGF signaling pathway. Numerous studies have shown that quercetin plays an important role in anti-inflammatory, antiproliferative, and antiatherosclerotic, especially in reducing blood fat and increasing coronary blood flow, which was contained by HH, TR, JG, WY, CH, and ZQ [40,47]. Luteolin acting on PTGS1, PTGS2, AR, and APP could depress hypertensive aorta remodeling in spontaneous hypertensive rats, relax porcine coronary and aortic arteries, and improve cardiac function after myocardium ischemia/reperfusion injury [41,42], which is vital to the treatment of cardiovascular diseases contained by HH, JG, XF, and ZQ. Butylated hydroxytoluene acting on PTGS1, PTGS2, AR, ESR1, and CYP3A4 is an antioxidant that has been shown to have an antiatherosclerotic effect in normal cholesterol-fed rabbits, which was contained by CS [48]. Kaempferol acting on AR, PIK3CG, ESR1, ESR2, CYP3A4, and APOE was involved in regulation of lipolysis in adipocytes, prolactin signaling pathway, estrogen signaling pathway, arachidonic acid metabolism, and VEGF signaling pathway, which is an important active member of dietary flavonoids discovered in many kinds of fruits and vegetables. Che et al. [49] demonstrated that kaempferol could effectively protect human umbilical vein endothelial cells against oxidative stress induced by ox-LDL, which indicated that kaempferol is capable of preventing atherosclerotic vascular disease. Rosmarinic acid acting on PTGS2, AR, ESR1, IL4, and NR1H4 was the active compound in DS and HH, which could inhibit the formation of malondialdehyde in human platelets in vitro, and its IC50 value is 3.37 nmol/L, indicating that rosmarinic acid has antiplatelet aggregation activity [50]. Berberine is a new cholesterol-lowering drug, and studies have shown that oral administration of berberine in 32 hypercholesterolemia patients for 3 months reduces serum cholesterol and triglycerides by 29% and 35%, respectively [51,52].
Obviously, though the analysis of hemorheological abnormality module revealed that quercetin, luteolin, isorhamnetin, butylated hydroxytoluene, isoliquiritigenin, nicotinic acid, kaempferol, rosmarinic acid, campesterol, and berberine may play an important role in improving hemorheological abnormality, the hub targets of QXH in the treatment of QS-BSS were PTGS2, PTGS1, AR, ESR1, CHRM5, ESR2, MMP2, PIK3CG, CYP3A4, and HMGCR. 8 Evidence-Based Complementary and Alternative Medicine Quercetin was the active compound of HH, WY, CH, and ZQ; luteolin was the active compound of HH, JG, XF, and ZQ; isorhamnetin was only contained by CH; butylated hydroxytoluene was contained by CS; nicotinic acid was the active compound in DG; rosmarinic acid was contained by DS and HH; and campesterol and berberine were contained by TR and YHS, respectively. Kaempferol was the active compound of HH, CH, and ZQ. Based on this, we can infer that HH, CS, DG, TR, YHS, CH, and ZQ were exerting effects for improving hemorheological abnormality.
KEGG analysis showed that the targets (related to hemorheological abnormality) mainly involved in PPAR signaling pathway, fat digestion and absorption, arachidonic acid metabolism, prolactin signaling pathway, inflammatory mediator regulation of TRP channels, regulation of lipolysis in adipocytes, and TNF signaling pathway.

H C -C C -T C Network of QXH (Coagulopathy Module).
Interestingly, quercetin, isorhamnetin, luteolin, kaempferol, berberine, wogonin, and isoliquiritigenin were not only playing a significant role in regulating hemorheological abnormality but also exerting effects for improving coagulopathy (Table 4). In Figure 5 (original picture see Figure S8), this network consisted of 267 nodes (15 herbs, 184 candidate compounds, and 68 candidate targets, Tables S2 and S3). e degree and betweenness of the top 10 candidate compounds and targets involved in improving coagulopathy are shown in Tables 4 and 5. Tangeretin was one of the active compounds of ZQ, which interacted with the hub targets (such as PTGS1, PTGS2, IL6, F2, F10, NQO1, and HMOX1) in coagulation function module. Columbamine can significantly reduce the levels of TC, TG, and LDL-c in the serum of golden hamsters (HFHC) induced by high fat and high cholesterol, while increasing the level of HDL-c [53]. Ellagic acid was one of the main active compounds of CS, which acted on HMOX1, NQO1, F12, SYK, and NOS3, and was involved in PI3K-Akt signaling pathway and platelet activation. Tanshinone IIB acted on F2, PTGS2, MPI ,and other targets with antiplatelet aggregation activity [54]. Betaine was proved to protect against coagulation events in vivo and vitro, decreasing lipid peroxidation in plasma [55].       are involved in the generation of thromboxane A2, which promote platelet activation and aggregation, vasoconstriction, and proliferation of vascular smooth muscle cells [56].
Tangeretin was the active compound in ZQ. Ferulic acid was shared in DG, CX, and SM. Ellagic acid was only contained by CS. Columbamine was the active compound in YHS, and wogonin was contained by NX. Hence, we can infer that ZQ, CX, and SM also contributed to improving coagulopathy in addition to HH, DG, CS, TR, and YHS.
KEGG analysis ( Figure 6) showed that the targets (related to coagulopathy) were mainly involved in complement and coagulation cascades, VEGF signaling pathway, vascular smooth muscle contraction, calcium signaling pathway, arachidonic acid metabolism, and leukocyte transendothelial migration. Over recent years, several studies have linked complement and thrombosis to the complement and coagulation cascades [57]. e thrombosis consists of insoluble fibrin, deposited platelets, accumulated white blood cells, and trapped red blood cells, which lead to the formation of blood stasis syndrome. VEGF signaling pathway leads to endothelial cell apoptosis, and a study has shown that the genetic repression of mouse VEGF expression regulates coagulation cascade [58]. Moreover, calcium signaling pathway plays a vital role in the development of platelet activation [59].

Analysis of the Overlapped Compounds, Targets, and Pathways in the Pathological Processes of Hemorheological
Abnormality and Coagulopathy. By analysis of the overlapped compounds, targets, and pathways involved in hemorheological abnormality and coagulopathy, 174 of 195 active compounds, 17 of 102 targets (Table S4), and 6 of 50 pathways associated with QXH were shared in the improvement of hemorheological abnormality and coagulopathy, indicating QXH in treating QS-BSS by the synergy of multiple compounds, targets, and pathways ( Figure 7). en, the biological process of overlapped targets was analyzed, which indicated that the pathological processes of hemorheological abnormality and coagulopathy were involved in inflammatory response, response to oxidative stress, cyclooxygenase pathway, angiogenesis, immune system process, regulation of blood pressure, and response to hypoxia of total 6 biological processes. Hence, in addition to hemorheological abnormality and coagulopathy, the participation of inflammation and immunity in the formation of QS-BSS is also worthy of our attention. Complement and coagulation cascades, arachidonic acid metabolism, VEGF signaling pathway, leukocyte transendothelial migration, Fc epsilon RI signaling pathway, and pathways in cancer might be the hub pathways of QXH in treating QS-BSS, which were shared in improving hemorheological abnormality and coagulopathy.

Discussion
TCM is a precious heritage which has been practiced for thousands of years in China and surrounding countries. Nowadays, herbal medicine has been widely used for complicated disease treatment and is gradually becoming acceptable as an alternative medicine worldwide [60]. TCM network pharmacology was a new research strategy for translating TCM from an experience-based medicine to an evidence-based medicine system, which predicted the target profiles and pharmacological actions of herbal compounds and revealed drug-gene-disease comodule associations to interpret the combinatorial rules and network regulation effects of herbal formula [61]. erefore, we herein utilized a combination method that included pharmacodynamics evaluation and network pharmacology to elucidate the pharmacodynamic effect, substances, and potential mechanisms of QXH through improving hemorheology and coagulation function in treating QS-BSS.
Hemorheology and coagulation function provide an important approach for studies on the mechanisms of QS-BSS and the therapeutic effects of Chinese herbal compounds that promote blood circulation and remove blood stasis [16]. In our work, we established the QS-BSS model in rats to evaluate the pharmacodynamic effect of QXH in improving hemorheological abnormality and coagulopathy. CDDP, one Chinese patent drug, was used as the positive control drug in this study. Preclinical and clinical studies have suggested that CDDP can expand blood vessel, promote blood circulation, relieve blood stasis, improve microcirculation, and improve hemorheological property [62]. e results of pharmacodynamic evaluation showed that QXH and CDDP groups significantly regulated the hemorheological abnormality and coagulopathy in rats with QS-BSS. Especially, QXH-L and QXH-H groups could significantly improve the abnormal increase of WBV and PV in rats with QS-BSS and had different degrees of callback on APTT, TT, PT, and FIB, which are comparable with CDDP in improving hemorheology and coagulation function. In terms of PCA and PLS-DA, after the preadministration of QXH, the same parameters in the QS-BSS model were close to those in the control and CDDP groups, also indicating that QXH has significant pharmacodynamic effect on QS-BSS. All of the above results were consistent with those of previous studies [63][64][65]. Hence, QXH, through promoting blood circulation and removing blood stasis in treating QS-BSS, gives us a glance at TCM compound formula in the treatment of the disease with a variety of pathogenesis by the synergy of various herbal medicines. ereby, we performed a network pharmacology approach which contains OB, Caco-2, and DL screening, active compound verification, target prediction, and network analysis to dissect the active compounds and possible mechanisms of QXH on regulating hemorheological abnormality and coagulopathy. e main results are as follows: (1) 195 active compounds in 15 herbs and 102 targets were identified as candidate targets related to QS-BSS. ese active compounds provided clues to our study of the substance on the basis of the pharmacodynamic effect of QXH.
(2) Meanwhile, 25 active compounds were verified by UPLC-Q/TOF-MS on QXH methanol and water extracts, which are corresponding to the results of network pharmacology. (3) rough the network analysis of QXH on the two vital pathological processes of QS-BSS, the mechanism of QXH in the treatment of QS-BSS is gradually clear.
In the hemorheological abnormality network module, 51 targets participated in signaling pathways such as PPAR signaling pathway, prolactin signaling pathway, fat digestion and absorption, estrogen signaling pathway, VEGF signaling pathway, and leukocyte transendothelial migration. A recent study on antihyperlipidemic activity of quercetin (included in HH, NX, WY, YHS, CH, XF, and ZQ) revealed that it may be involved in regulating the PPAR signaling pathway [66], which is comparable to our findings. e PPAR signal pathway is concerned with the formation of QS-BSS. PPARc is highly expressed in macrophage-derived foam cells of both early and advanced atherosclerotic lesions [67]. PPARc ligands may attenuate inflammation and hence,    Figure 5: e H C -C C -T C network of QXH (coagulopathy module, red node: herbs, green node: compounds, and blue node: targets). atherosclerosis in the vessel wall. But the PPARc expression in monocytes increased the expression of CD36, leading to an enhanced macrophage uptake of oxidized LDL, thus enhancing foam cell formation [68].
In the coagulopathy network module, 68 targets participated in signaling pathways such as complement and coagulation cascades, VEGF signaling pathway, vascular smooth muscle contraction, calcium signaling pathway, arachidonic acid metabolism, and long-term depression. F2, F3, F7, F9, F10, SERPINE1, and SERPIND1 were included by complement and coagulation cascades, which acted on blood coagulation, platelet activation, blood homeostasis, inflammation, and wound healing. Studies have shown that APTT and PT reflect changes in the coagulation factor of the endogenous and exogenous coagulation pathways, respectively. TT and FIB are related to the common coagulation pathway in plasma [36]. However, the analysis results of coagulopathy network module fully support the rationality of using APTT, TT, PT, and FIB to evaluate QXH on improving coagulation function in QS-BSS. Another study on identifying the compounds of DS on intervening thrombotic diseases revealed that tanshinone IIB could improvemicrocirculation and prevent thrombosis involved in complement and coagulation cascades based on UPLC-LTQ-Orbitrap MS [69], which is also comparable to our findings.
Analysis of the targets overlapped such as PTGS1, PTGS2, PIK3CG, and PECAM1 in the pathological processes of hemorheological abnormality and coagulopathy implied that inflammation and immunity were also important for QS-BSS which need further study. Complement and coagulation cascades, arachidonic acid metabolism, and VEGF signaling pathway were shared in hemorheological abnormality and coagulopathy, which need more attention, and further experimental validation will be carried out.
However, although we have predicted the main active compounds of QXH through network pharmacology and some of which has been verified by UPLC-Q/TOF-MS, to clarify the substance basis of the pharmacodynamic effect of QXH, further study on the migrating components absorbed into blood by serum pharmacochemistry based on animal experiments is needed [70]. Meanwhile, the network analysis  Evidence-Based Complementary and Alternative Medicine 13 also provided a vital reference for experimental investigations of hub targets and signaling pathways of QXH in treating QS-BSS.

Conclusions
rough the pharmacodynamic evaluation of QXH based on animal experiments and the underlying pharmacological mechanism prediction of QXH, our study has proved that the effect of QXH in the treatment of QS-BSS is significant, and the underlying pharmacological mechanisms might be related to complement and coagulation cascades, leukocyte transendothelial migration, PPAR signaling pathway, VEGF signaling pathway, arachidonic acid metabolism, inflammatory mediator regulation of TRP channels, and TNF signaling pathway, for which the further experimental validation is essential.

QXH:
Qixuehe capsule QXH-L: Low-dose QXH-intervened QXH-M: Middle-dose QXH-intervened QXH-H: High-dose QXH-intervened TCM: Traditional Chinese medicine Adr: Adrenaline hydrochloride CDDP: Compound danshen dripping pills PT: Prothrombin time APTT: Activated partial thromboplastin time TT: rombin time FIB: Fibrinogen content WBV: Whole blood viscosity PV: Plasma viscosity PCA: Principal component analysis PLS: Partial least squares CV-ANOVA: Analysis of variance testing of cross-validated predictive residuals UPLC: Ultra-performance liquid chromatography MS: Mass spectrometer TCMSP: Traditional Chinese medicine systems pharmacology database and analysis platform TTD: erapeutic targets database CTD: Comparative toxicogenomics database SEA: Similarity ensemble approach OMIM: Online Mendelian inheritance in man database OB: Oral bioavailability Caco-2: Caco-cell permeability DL: Drug-likeness KEGG: Kyoto Encyclopedia of Genes and Genomes   Table S1: the OB, DL, and Caco-2 of the active compounds in QXH. Table S2: the active compounds involved in hemorheological abnormality module and coagulopathy module. Table S3: the candidate targets involved in hemorheological abnormality module and coagulopathy module. Table S4: the overlapped targets in hemorheological abnormality module and coagulopathy module. Figure S1: PCA score plots of hemorheology and coagulation function indexes among all groups. Figure S2: the permutation tests (n = 200) (A) and CV-AVNOA tests (B) for the PLS-DA model among all groups. Figure S3: PLS-DA score plots of hemorheology and coagulation function indexes between control and model groups. Figure S4: the permutation tests (n = 200) (A) and CV-AVNOA tests (B) for the PLS-DA model between control and model groups. Figure S5: base peak intensity chromatograms of 25 standards in positive ion mode or negative ion mode. Figure S6: base peak intensity chromatograms of QXH methanol extracts (A) in positive ion mode and (C) in negative ion mode; base peak intensity chromatograms of QXH water extracts (B) in positive ion mode and (D) in negative ion mode by UPLC-Q/TOF-MS. Figure S7: the HH-CH-TH network of QXH (hemorheological abnormality module; red node: herbs, green node: compounds, and blue node: targets) (the original picture for Figure 3). Figure S8: the HC-CC-TC network of QXH (coagulopathy module, red node: herbs, green node: compounds, and blue node: targets) (the original picture for Figure 5). (Supplementary Materials)