Skip to main content

Network pharmacology prediction and molecular docking-based strategy to discover the potential pharmacological mechanism of Huang–Qi–Gui–Zhi–Wu–Wu decoction against deep vein thrombosis

Abstract

Background

Huangqi Guizhi Wuwu decoction (HQGZWWD) has been used to treat and prevent deep vein thrombosis (DVT) in China. However, its potential mechanisms of action remain unclear. This study aimed to utilize network pharmacology and molecular docking technology to elucidate the molecular mechanisms of action of HQGZWWD in DVT.

Methods

We identified the main chemical components of HQGZWWD by reviewing the literature and using a Traditional Chinese Medicine Systems Pharmacology (TCMSP) database. We used GeneCards and Online Mendelian Inheritance in Man databases to identify the targets of DVT. Herb-disease-gene-target networks using Cytascape 3.8.2 software; a protein–protein interaction (PPI) network was constructed by combining drug and disease targets on the STRING platform. Additionally, we conducted Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses. Finally, molecular docking verification of active components and core protein targets was conducted.

Results

A total of 64 potential targets related to DVT were identified in HQGZWWD, with 41 active components; quercetin, kaempferol, and beta-sitosterol were the most effective compounds. The PPI network analysis revealed that AKT1, IL1B, and IL6 were the most abundant proteins with the highest degree. GO analysis indicated that DVT treatment with HQGZWWD could involve the response to inorganic substances, positive regulation of phosphorylation, plasma membrane protein complexes, and signaling receptor regulator activity. KEGG analysis revealed that the signaling pathways included pathways in cancer, lipid and atherosclerosis, fluid shear stress and atherosclerosis, and the phosphatidylinositol 3-kinases/protein kinase B(PI3K-Akt) and mitogen-activated protein kinase (MAPK) signaling pathways. The molecular docking results indicated that quercetin, kaempferol, and beta-sitosterol exhibited strong binding affinities for AKT1, IL1B, and IL6.

Conclusion

Our study suggests that AKT1, IL1B, and IL6 are promising targets for treating DVT with HQGZWWD. The active components of HQGZWWD likely responsible for its effectiveness against DVT are quercetin, kaempferol, and beta-sitosterol, they may inhibit platelet activation and endothelial cell apoptosis by regulating the PI3K/Akt and MAPK signaling pathways, slowing the progression of DVT.

Introduction

Deep vein thrombosis (DVT) is a disease in which blood in the deep veins of the limbs abnormally coagulates into clots due to damage to the venous wall and stagnation of blood flow, leading to narrowing or occlusion of the venous lumen [1, 2]. The incidence of DVT is high; for example, hospital-acquired DVT after major orthopedic surgery can reach 60% [3,4,5]. The danger lies in the fact that detachment of a DVT can cause fatal pulmonary embolism with a high mortality rate; it is a crucial cause of perioperative and unexpected hospital deaths [6, 7]. Low-molecular-weight heparin and vitamin K antagonists, such as warfarin and sulodexide, are routinely used in clinical practice to prevent DVT [8]; however, adverse reactions, mainly bleeding, hematoma formation, and decreased hemoglobin concentrations, have been discovered with the long-term, extensive clinical use. Furthermore, there is controversy on anticoagulants preventing DVT in patients with cerebral hemorrhage or combined blood system diseases [9]. Therefore, the effective prevention of DVT in long-term bedridden patients without causing new complications remains an urgent problem that needs to be solved.

Traditional Chinese medicine has a long history of use for the prevention and treatment of DVT [10]. Traditional Chinese medicine categorizes DVT as “pulse obstruction” and “blood stasis.” [11] Sun Simiao, a renowned Chinese physician over a thousand years ago, noted in Qian Jin Bei Ji Yao Fang that poor blood circulation is the root cause of thrombosis: “if qi and blood are stagnant, there will be pain; if the pulse is blocked, there will be swelling, and if stagnation persists for a long time, heat will arise.” HQGZWWD is composed of five traditional Chinese medicines: Huangqi (Hedysarum multijugum maxim, HM), Guizhi (Cinnamomi Ramulus, CR), Baishao (Paeoniae Radix Alba, PA), Shengjiang (Zingiber officinale Roscoe, ZR), and Dazao (Jujubae Fructus, JF) [12]. HQGZWWD is a medicinal formula that promotes blood circulation and nourishes the qi; this was reported by Zhang Zhongjing, a medical expert from the Eastern Han Dynasty (approximately 154–219 AD), in his book “Jin Gui Yao Lue” almost 2000 years ago [13]. Additionally, various medical books such as “San Yin Ji Yi Bing Zheng Fang Lun” from the Southern Song Dynasty (1174 AD), “Zheng Zhi Zhun Sheng: Lei Fang” from the Ming Dynasty (1602 AD) and “Yi Fang Ji Jie” from the Qing Dynasty recorded the effectiveness of HQGZWWD in treating DVT. Studies such as Yuebao have shown that using HQGZWWD can enhance hypercoagulability in patients who have undergone knee replacement surgery while decreasing the likelihood of developing postoperative DVT [14]. Furthermore, Zhao Zhili used HQGZWWD to treat lung cancer accompanied by DVT, significantly improving patients’ clinical symptoms and considerably shortening the acute or subacute period of DVT [15]. Despite its historical significance, our understanding of its mechanism primarily relies on the traditional Chinese medicine experience passed down through generations and clinical observations, without any scientific research regarding its molecular mechanisms.

The complexity of traditional Chinese medicine stems from the lack of adequate quantitative evidence to assess its therapeutic efficacy. Typically, traditional Chinese medicine involves intricate formulas comprising multiple herbs, with their composition and dosage often based on ancient texts and empirical knowledge rather than on contemporary scientific research [16]. Its peculiarity lies in its exclusive history within China, which has been subject to limited theoretical and fundamental investigations until recently. Consequently, many Western Scholars may not regard it as a science but rather an empirical practice lacking any theoretical foundation; these factors pose unique challenges to traditional Chinese medicine research [17].

Network pharmacology is a novel analytical approach that utilizes virtual computing and database retrieval to study the mechanisms of disease and drug action within a broader biological network. This method provides valuable insights into the pharmacological efficacy and mechanisms of the drugs [18]. The primary objective of network pharmacology research is to systematically address scientific challenges at multiple levels, which closely aligns with the fundamental concept of treating prediseases in traditional Chinese medicine, known as syndrome differentiation and treatment [19]. As an emerging research methodology, it offers a promising solution to overcome obstacles such as inadequate basic research on traditional Chinese medicine [20]. Based on this situation,We utilized network pharmacology methods to establish an herb-disease-gene-target network and constructed a protein–protein interaction (PPI) network combining drug and disease targets. We conducted a Gene Ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to identify the primary signaling pathways and biological functions. Finally, molecular docking was performed to verify functional components and their major targets. In summary, our study investigates the main components of HQGZWWD that act on DVT and their molecular biological mechanisms using network pharmacology and molecular docking methods. This study provides a reference for future research in this field.

Materials and methods

Screening of the active components in HQGZWWD

The TCMSP platform (http://tcmspw.com/tcmsp.php) was used to identify the active components of HQGZWWD (accessed on March 15, 2023). Our selection criteria were based on the absorption, distribution, metabolism, and excretion (ADME) processes. Oral bioavailability (OB) is the amount of a drug that enters the circulation after entering the human body. However, drug-likeness (DL) refers to the degree of similarity between a compound and known drugs. OB is a representative pharmacokinetic parameter in ADME, and DL is used as a qualitative concept in drug design to estimate the molecular characteristics of drugs [21]. Specifically, we screened for components with oral bioavailability (OB) ≥ 30% and drug-like properties (DL) ≥ 0.18 to obtain the appropriate active compounds [22,23,24,25].

Construction of drug active component target network and identification of DVT-predictive targets

We identified the active components of HQGZWWD by applying predefined thresholds (OB ≥ 30% and DL ≥ 0.18) to the TCMSP database and individually confirmed the respective targets. To standardize the target and gene symbols, we used the UniProt database (accessed on March 16, 2023) for conversion (https://www.uniprot.org/). To obtain disease targets from the databases, we searched for “deep venous thrombosis” as a keyword in the OMIM (https://omim.org/, accessed on March 16, 2023) and Gene cards databases (https://www.genecards.org/, accessed on March 16, 2023) [26, 27]. Next, we merged these targets, removed duplicates, and identified the remaining DVT targets. Using Cytascape (ver.3.8.2) network visualization software, we constructed a “herb-disease-gene-target” network of effective components and action targets. Finally, we analyzed network characteristics to elucidate the interactions between effective components and targets in herbs and diseases.

Construction and analysis of the PPI network

The Search Tool for Retrieving Interacting Genes/Proteins (STRING) (https://string-db.org/) was used to predict protein–protein interactions (accessed on March 18, 2023). The PPI network was created by introducing overlapping targets [28]. CytoScope (ver.3.8.2) was used to construct a core PPI network.

GO enrichment analysis and KEGG pathway analysis

To elucidate the biological functions and signaling pathways associated with DVT, we conducted GO and KEGG enrichment analyses using the Metascape database (https://metascape.org/, accessed on March 18, 2023) [29]. GO analysis identified relevant biological processes (BP), cellular components (CC), and molecular functions (MF). Additionally, KEGG enrichment analysis enabled us to identify significant signaling pathways involved in these biological processes. The P-value cutoff was established as P<0.01, where a lesser P-value signifies a higher likelihood of the current result being an authentic enrichment outcome rather than a random occurrence. The unit employed in this study was -log 10(P-value), and as this value increased, the reliability of the enrichment results increased.

Molecular docking

Molecular docking is a theoretical simulation method used in drug design to predict the binding modes and affinities of molecules by studying their interactions with receptors and ligands [30]. AutoDock Tools is a simulation software that facilitates the study of the interactions between biomolecules and small-molecule complexes; this enables researchers to accurately comprehend how protein targets interact with small-molecule compounds by simulating the ligand-receptor recognition [31]. In this study, we employed molecular docking to investigate whether the core components of HQGZWWD, identified through network pharmacology, could bind to core proteins. We selected the top three compounds based on degree value from the “herb-disease-gene-target” network's core components of HQGZWWD and chose the top three proteins based on degree value from the PPI network's core targets. The corresponding 3D structure files for proteins and small molecule compounds were downloaded from the RCSB database (https://www.rcsb.org/, accessed on March 19, 2023) and TCMSP database [32], followed by dehydration and hydrogenation before importing them into Autodock Tools (ver.1.5.6). The small molecule compound was then docked with the protein as a receptor using a molecular docking method to select the most suitable conformation. Finally, to generate a binding mode diagram, we imported the docked protein and small molecule files prepared using AutoDock Tools into PyMOL (ver.1.8.x). Using the 'find-polar contacts-to any atom' function, we can automatically identify hydrogen bonds or other interactions [33].

Results

Extraction of the active components in HQGZWWD

According to the screening criteria, OB ≥ 30% and DL ≥ 0.18 were set in the TCMSP database for screening; the results showed that 41 active components of HQGZWWD were obtained, as shown in Additional file 1.

Identification of core targets and construction of herb-disease-gene-target network

We searched the OMIM and GeneCards databases for potential targets related to DVT, removed duplicate targets, and identified 1731 potential targets. We then sorted the results from TCMSP, removed duplicate sites, matched the predicted drug targets with the corresponding disease targets, and obtained 64 core targets (Fig. 1). To further observe the compound-target interactions, we constructed a "herb-disease-gene-target" network (Fig. 2). The degree refers to the number of nodes that interact with a node, and its size represents its core degree. The higher the degree value, the more important this node is. Our findings indicate that three compounds in HQGZWWD are highly correlated with DVT: quercetin (C2; degree = 246), beta-sitosterol (C3; degree = 100), and kaempferol (C1; degree = 94).

Fig. 1
figure 1

Venn diagram showing the overlapping target genes for HQGZWWD against DVT. HQGZWWD, Huang-Qi-Gui-Zhi-Wu-Wu Decoction; DVT, Deep vein thrombosis

Fig. 2
figure 2

Herb-disease-gene-target network of HQGZWWD against DVT. The larger the font size, the more important its role in the compound. HQGZWWD Huang-Qi-Gui-Zhi-Wu-Wu Decoction, DVT Deep vein thrombosis, HM Hedysarum multijugum Maxim, CR Cinnamomi Ramulus, PA Paeoniae Radix Alba, ZR Zingiber officinale Roscoe, JF Jujubae Fructus. (C1)kaempferol, common components of HM and PA, (C2)quercetin, common components of HM and JF, (C3)beta-sitosterol, common components of CR, PA, ZR, and JF, (C4)sitosterol, common components of CR and PA, (C5)mairin, common components of HM, PA, and JF, (C6)( +)-catechin, common components of CR, PA, and JF, (C7)stigmasterol, common components of ZR and JF

Construction and analysis of the PPI network

The core overlapping diseases and drug component targets were imported into the STRING database to generate a PPI network diagram (Fig. 3). Molecular complex detection (MCODE) is a novel graph-theoretic clustering algorithm that calculates the information of each node in a PPI network graph using a k-means clustering algorithm to detect dense connection regions in large protein-protein interactions. Proteins within these regions may have similar structures or functions, providing a reference and guidance for further research on protein-disease interactions [34]. Target information from this network was further analyzed using MCODE, which revealed three central gene clusters; these clusters suggested that the proteins within them were closely connected and may share common functions or expression patterns (Fig. 4). Key target proteins were identified using the Cytoscape software, with higher degree values indicating stronger interactions with other targets. This analysis highlighted the potential central targets for HQGZWWD in DVT treatment, including AKT1 (degree = 51), IL6 (degree = 52), and IL1B (degree = 50) (Fig. 5).

Fig. 3
figure 3

PPI network. Empty nodes represent proteins of unknown 3D structures; filled nodes represent some 3D structures that are known or predicted. Edges represent protein–protein associations: the light blue edges represent from curated databases; the fuchsia edges represent experimentally determined; the green edges represent gene neighborhood; the red edges represent gene fusions; the dark blue edges represent gene co-occurrence; the light green edges represent text mining; the black edges represent co-expression; the light purple edges represent protein homology

Fig. 4
figure 4

Central gene cluster identified in the PPI network of HQGZWWD-DVT based on MCODE analysis. ac Was cluster 1, 2, and 3, respectively

Fig. 5
figure 5

PPI network diagram processed by Cytascape. The color of the target point changes gradually according to the degree value. The higher the degree value, the larger the circle. As the degree value changes, the color changes from light purple to deep purple

GO enrichment analysis and KEGG pathway analysis

GO and KEGG enrichment analyses were conducted to better understand the potential pathways and biological functions of HQGZWWD in DVT. Figure 6 displays the top five results of the GO analysis, which revealed that HQGZWWD primarily affected the response to inorganic substances (BP), positive regulation of phosphorylation (BP), plasma membrane protein complex (CC), and signaling receptor regulator activity (MF). Additionally, our KEGG analysis identified several signaling pathways involved in this treatment approach, including pathways in cancer, lipid and atherosclerosis, fluid shear stress and atherosclerosis, the PI3K-Akt pathway, and the MAPK signaling pathway (Fig. 7a, b).

Fig. 6
figure 6

GO analysis of key target genes. The top 5 items of biological function are listed on the vertical axis, including GOMF, GOCC, and GOBP, the horizontal axis in the figure represents the gene ratio. GO Gene ontology, BP biological process, CC cell composition

Fig. 7
figure 7

a KEGG analysis of key target genes, b Network of top 20 pathways. Red diamond represents gene, and green triangle represents pathway. The size of the nodes represents the value of the degree. The horizontal axis in the figure represents the gene ratio. KEGG Kyoto Encyclopedia of Gene and Genome

Molecular docking

We employed network pharmacology and PPI analyses to determine the top three compounds, quercetin, beta-sitosterol, and kaempferol, for docking with the top three proteins, AKT1, IL6, and IL1B. The best docking image of the receptor-ligand complex is presented in Fig. 8.

Fig. 8
figure 8

Molecular docking results of main chemical components of HQGZWWD and core proteins in PPI network a IL1B-beta-sitosterol, b IL1B-quercetin, c IL1B-kaempferol, d AKT1-beta-sitosterol, e AKT1-quercetin, f AKT1-kaempferol, g IL6-beta-sitosterol, h IL6-quercetin, i IL6-kaempferol

Our results indicate that beta-sitosterol can form hydrogen bonds with ALA-9 and GLY-101 through IL-1B protein (Fig. 8a), while quercetin can do so through GLU-25, TYR-24, ARG-91, and SER-30 (Fig. 8b), and kaempferol through SER-30, ARG-91, ASN-101, LEU-134, and PRO-78 (Fig. 8c). Similarly, for AKT1, beta-sitosterol formed hydrogen bonds via LYS-130 (Fig. 8d) and quercetin via ASP-348 and ALA-390 (Fig. 8e), while kaempferol formed hydrogen bonds via ARG-169/170 and ASP-103 (Fig. 8f). Finally, for the IL6 protein: beta-sitosterol forms hydrogen bonds via TRP-112 (Fig. 8g); quercetin does so through GLN-41/40, ASP-87, and LYS-170 (Fig. 8h), while kaempferol forms bonds between PHE-101/TR-P-49 and ALA-99 (Fig. 8i).

The tighter the binding between the ligands and receptors, the smaller the binding energy. Notably, binding energy < -5 kcal/mol indicates good affinity and binding activity between the receptor and ligand [35, 36]. Our molecular docking results showed that all three selected compounds (quercetin, beta-sitosterol, and kaempferol) had high affinities for the three core targets (AKT1, IL6, and IL1B), as their binding energies to proteins were < -5 kcal/mol (Table 1); this suggests that these compounds may play crucial roles in the treatment of DVT.

Table 1 Binding energy of molecular docking

Discussion

This study investigated the mechanism by which HQGZWWD treats DVT. By constructing and analyzing databases and network diagrams of traditional Chinese medicine prescriptions, we identified quercetin, beta-sitosterol, and kaempferol as highly correlated with DVT. Previous studies demonstrated that quercetin and its derivatives possess anticoagulant, antiplatelet, and antifibrinolytic activities [37, 38], effectively preventing pulmonary thromboembolism [39]. Furthermore, onion extracts rich in quercetin regulate MAPK under coagulation stimulation to prolong thrombosis time [40], consistent with our KEGG enrichment analysis results. Kaempferol and quercetin are flavonoids that inhibit prothrombin activity while regulating fibrinogen-thrombin interactions for antithrombotic formation in vivo and in vitro [41]. In 2018, researchers reported that plant-derived beta-sitosterol exhibits antithrombotic properties, such as in vivo anticoagulant effects and thromboprophylaxis [42]. These findings suggest that these components play an essential role in the treatment of DVT using HQGZWWD. However, research on the mechanisms of action of beta-sitosterol and kaempferol in thrombosis is still at an early stage and warrants further exploration.

The results of the PPI network analysis suggested that AKT1, IL6, and IL1B could be potential core targets for treating DVT with HQGZWWD; IL-1B and IL-6 are proinflammatory cytokines with a wide range of biological activities [43]. Initially, it was believed that IL6 did not play a role in venous thrombosis; however, as more research was conducted on its function, this statement was discarded [44]. Studies have revealed that patients with DVT have elevated levels of IL-6, and in vivo and in vitro studies have confirmed its involvement in DVT formation [45, 46]. Similar to psoriasis, studies have shown that IL-6 is involved in regulating inflammation-related thrombosis [47]. Regarding IL-1B, research has indicated increased body levels during acute pulmonary embolism, which can be reversed by atorvastatin treatment [48]. Currently, most studies on cytokines focus on their effects on atherosclerosis and cardiovascular disease [49,50,51]; however, research on their direct involvement in venous thrombosis started late. Further exploration is needed to determine the specific mechanism underlying these findings, which were significant in our study.

We conducted GO and KEGG analyses to investigate the mechanism of action of HQGZWWD in DVT. The results of the GO analysis showed that HQGZWWD treatment primarily involved responses to inorganic tolerance, positive regulation of phosphorylation, plasma membrane protein complex, and signaling receiver regulator activity. KEGG analysis revealed that multiple signaling pathways are involved in HQGZWWD treatment of DVT, including pathways in cancer, lipid and atherosclerosis, fluid shear stress, atherosclerosis, PI3K-Akt, and MAPK signaling pathways. The involvement of quercetin in cancer pathogenesis may be the main reason for its enrichment in cancer [52,53,54]. A previous study showed that the lack of protein tyrosine phosphatase non-receptor type 22 (PTPN22) significantly accelerates thrombus formation. Further research has shown that in platelets with PTPN22 defects, the phosphorylation of phosphodiesterase 5A (PDE5A) increases, and the cyclic guanosine monophosphate levels and vasodilator-stimulated phosphoprotein phosphorylation decrease [55]. Other studies have demonstrated that the platelet phenotypic response caused by Glycoprotein VI (GPVI) activation is linked to changes in protein kinase substrate phosphorylation. To investigate how adenosine diphosphate (ADP) secretion and thromboxane generation feedback affect GPVI activation and signal transduction, proteomic analysis was conducted; the results revealed significant alterations in the phosphorylation levels following ADP secretion, thromboxane feedback, and collagen-related peptide (CRP-XL) stimulation. The significantly altered proteins were primarily involved in biological processes such as vesicle-mediated transport and cell skeleton reconstruction, which cannot be separated from the role of quality model protein complexes and signal regulation receptor activity [56]. These identified biological processes, cellular components, and molecular functions aligned with our GO analysis results: response to inorganic tolerance, positive regulation of phosphorylation, plasma membrane protein complex, and signaling receiver regulator activity; therefore, it is worth investigating whether HQGZWWD exerts antagonistic effects on DVT formation via these biological processes. Our confidence was enhanced by the preliminary prediction results of the network pharmacology.

Although atherosclerosis and venous thrombosis were once considered distinct pathological conditions, recent studies have challenged this notion [57]. A study published in the NEJM suggested that atherosclerosis may contribute to the development of venous thrombosis [58]. Furthermore, Hong et al. discovered that patients with idiopathic DVT had significantly higher rates of coronary artery calcification than controls [59]. Another study reported that patients with acute DVT had a high incidence of acute myocardial infarction or ischemic stroke-related death during long-term follow-up [60]. These findings suggest a potential synergistic relationship between atherosclerosis and venous thrombosis. For instance, some studies have demonstrated that atherosclerosis can affect the adhesion of thrombotic substances to the peripheral venous system [61]. Specifically, atherosclerosis primarily affects the arterial tree system and is caused by an imbalance in lipid metabolism and maladaptive immune responses, resulting in chronic inflammation of the blood vessel walls [62,63,64]. Furthermore, chronic inflammatory diseases that cause changes in blood vessel walls are associated with a higher risk of DVT [65, 66]. Shared risk factors such as inflammation, hypercoagulability, and endothelial injury increase the likelihood of chronic inflammation induced by processes related to atherosclerosis, affecting the development of DVT and adhesion of DVT clots within veins [67, 68]. Traditional Chinese medicine treatment methods have shown that HQGZWWD is effective against coronary heart disease caused by arteriosclerosis, cerebral infarction, and peripheral vascular disease [69, 70]. Furthermore, literature search has revealed quercetin, a compound in HQGZWWD, can combat the progression of arteriosclerotic plaques through multiple mechanisms, such as regulating oxidized low-density lipoprotein-induced endothelial cell senescence [71] and modulating the autophagy of macrophages [72]. Since HQGZWWD has a positive therapeutic effect on arteriosclerosis, which is known to have a mutually promoting effect on DVT, it may be worthwhile to investigate whether HQGZWWD can simultaneously inhibit the progression of both conditions; this aligns with one of the objectives of network pharmacology, drug repositioning, and drug development.

The PI3K/Akt signaling pathway is a crucial regulatory pathway playing a significant role in pathophysiological processes such as cell growth, differentiation, and proliferation [73]. Recently, researchers have focused on its correlation with endothelial cell mobilization [74], differentiation [75, 76], and apoptosis [77]. Studies have suggested that this pathway regulates vascular endothelial growth factor (VEGF) secretion and controls apoptosis in blood vessels and vascular endothelial cells [78]. Additionally, microRNA-126 has been found to inhibit the PI3K/Akt signaling pathway to prevent vascular endothelial cell apoptosis for therapeutic effects in DVT treatment [79]. Furthermore, FXII activates PI3K/AKT signaling and promotes DVT progression by inducing inflammatory reactions [80]. However, kaempferol, the main component identified in this study, significantly reduced the phosphorylation of PI3K/AKT during the collagen/adrenaline-stimulated platelet activation tests. This delay ultimately resulted in a 34.6% reduction in clotting time. Additionally, animal experiments demonstrated the thromboprotective effect of resveratrol in collagen/adrenaline- and thrombin-induced acute thrombotic embolism models and FeCl3-induced carotid artery thrombosis models [41]. The target we selected was enriched in the PI3K/Akt signaling pathway, and kaempferol was one of the main components screened, which validated the presence of AKT1 in the PPI network analysis and suggested promising directions for future research. HQGZWWD can potentially inhibit endothelial cell apoptosis by activating the PI3K/Akt pathway, which could potentially slow the progression of DVT. Kaempferol, the primary constituent of HQGZWWD, may play a crucial role in this mechanism; however, additional cellular and animal experiments are necessary to further investigate this mechanism.

The MAPK signaling pathway is a crucial component of the eukaryotic signal transduction network; it plays a key role in regulating cell proliferation, differentiation, apoptosis, and stress responses under normal and pathological conditions [81, 82]. Studies have demonstrated that a sustained increase in reactive oxygen species can activate inflammatory reactions by activating MAPK signaling pathways. This activation leads to apoptosis of vascular endothelial cells and induces thrombosis [83]. Furthermore, studies have shown that ginsenoside-Rp3 can regulate MAPK signaling pathways to inhibit platelet activation and thrombosis [84]. These findings support our hypothesis that HQGZWWD may delay DVT progression by modulating the MAPK signaling pathway. Previous research has indicated that the activation of Raf1 in the MAPK pathway can regulate thromboxane production and platelet thrombosis [85]. However, quercetin, the primary compound in HQGZWWD, has been shown to inhibit IL-1 and chemokine production via the MAPK signaling pathway [86]. IL-1-induced inflammatory factors can activate proinflammatory thrombogenic processes in vascular endothelial cells, stimulate angiogenic mediator production, and promote thrombosis [87, 88]. Therefore, investigating whether HQGZWWD and its main components inhibit thrombus formation by regulating these processes is essential; the findings of this study reveal its potential.

Molecular docking experiments revealed that quercetin, beta sitosterol, and kaempferol, the main compounds obtained from HQGZWWD, have binding energies with core target proteins AKT1, TP53, and TNF of < -5 kcal/mol (Table 1). This indicated that these three small molecule compounds may bind to these proteins and regulate signal transduction; they play vital roles in the treatment of DVT. However, the therapeutic effects and molecular mechanisms of these agents active in DVT require further investigation; we believe they possess significant potential for the clinical prevention and treatment of DVT.

This study had some limitations. Although we used modern bioinformatics methods to explore the role of HQGZWWD in treating DVT through network pharmacology and molecular docking, there is still a lack of molecular dynamics(MD) simulations to supplement molecular docking in determining the dynamic stability of receptor ligands [89, 90]. Additionally, this study relied solely on conventional computational analysis without utilizing more comprehensive omics datasets for mining targets. Although quercetin, sitosterol, and kaempferol have been identified as the three most important bioactive components for treating DVT, they do not represent all the compounds in HQGZWWD. Therefore, further research involving MD simulations, deeper omics data mining, and molecular biology experiments are necessary to validate our findings.

Conclusions

In summary, we adopted a new approach that combined network pharmacology and molecular docking to predict AKT1, IL1B, and IL6 as the most likely targets of HQGZWWD for treating DVT. Quercetin, kaempferol, and beta-sitosterol are the main active components of HQGZWWD; they may inhibit platelet activation and endothelial cell apoptosis by regulating the PI3K/Akt and MAPK signaling pathways, slowing the progression of DVT. Additionally, our results offer new insights into the pathogenesis of DVT and suggest potential avenues for developing innovative treatment strategies. The initial findings from network pharmacology have increased our confidence. However, further clinical experiments and in-depth basic research are necessary to fully understand the therapeutic effects and molecular mechanisms of traditional Chinese medicine on DVT; this will be the next area of our research.

Availability of data and materials

The data used and/or analyzed during the current study are available from the corresponding author upon reasonable request. The datascxsets generated and/or analysed during the current study are available in the followed database: TCMSP database (http://tcmspw.com/tcmsp.php). UniProt database (https://www.uniprot.org/). OMIM database (https://omim.org/). Gene cards database (https://www.genecards.org/). STRING database (https://string-db.org/). Metascape database (https://metascape.org/). RCSB database (https://www.rcsb.org/).

Abbreviations

HQGZWWD:

Huang–Qi–Gui–Zhi–Wu–Wu decoction

DVT:

Deep vein thrombosis

GO:

Gene ontology

BP:

Biological process

CC:

Cell composition

KEGG:

Kyoto encyclopedia of gene and genome

VEGF:

Vascular endothelial growth factor

OB:

Oral bioavailability

DL:

Drug-likeness

References

  1. Line BR. Pathophysiology and diagnosis of deep venous thrombosis. Semin Nucl Med. 2001;31(2):90–101. https://doi.org/10.1053/snuc.2001.21406.

    Article  CAS  PubMed  Google Scholar 

  2. Stubbs MJ, Mouyis M, Thomas M, et al. Deep vein thrombosis. BMJ. 2018;360:k351. https://doi.org/10.1136/bmj.k351.

    Article  CAS  PubMed  Google Scholar 

  3. Park SJ, Kim CK, Park YS, et al. Incidence and factors predicting venous thromboembolism after surgical treatment of fractures below the hip. J Orthop Trauma. 2015;29(10):e349–54. https://doi.org/10.1097/BOT.0000000000000336.

    Article  PubMed  Google Scholar 

  4. Dixon J, Ahn E, Zhou L, et al. Venous thromboembolism rates in patients undergoing major hip and knee joint surgery at Waitemata District Health Board: a retrospective audit. Intern Med J. 2015;45(4):416–22. https://doi.org/10.1111/imj.12702.

    Article  CAS  PubMed  Google Scholar 

  5. Saleh J, El-Othmani MM, Saleh KJ, et al. Deep vein thrombosis and pulmonary embolism considerations in orthopedic surgery. Orthop Clin North Am. 2017;48(2):127–35. https://doi.org/10.1016/j.ocl.2016.12.003.

    Article  PubMed  Google Scholar 

  6. Di Nisio M, van Es N, Büller HR, et al. Deep vein thrombosis and pulmonary embolism. Lancet. 2016;388(10063):3060–73. https://doi.org/10.1016/S0140-6736(16)30514-1.

    Article  PubMed  Google Scholar 

  7. Bernardi E, Camporese G. Diagnosis of deep-vein thrombosis. Thromb Res. 2018;163:201–6. https://doi.org/10.1016/j.thromres.2017.10.006.

    Article  CAS  PubMed  Google Scholar 

  8. Chen B, Hu N. Low molecular weight heparin and aspirin for prevention of deep vein thrombosisafter orthopaedic surgery: a systematic review and meta-analysis. J Thromb Thrombolysis. 2021;52(2):553–9. https://doi.org/10.1007/s11239-020-02348-5.

    Article  CAS  PubMed  Google Scholar 

  9. Murthy SB, Gupta A, Merkler AE, et al. Restarting anticoagulant therapy after intracranial hemorrhage: a systematic review and meta-analysis. Stroke. 2017;48(6):1594–600. https://doi.org/10.1161/STROKEAHA.116.016327.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Chen C, Tang Q, Zhang W, et al. Combination of traditional chinese medicine and low-molecular-weight heparin prevents deep vein thrombosis after surgery: a meta-analysis. Clin Appl Thromb Hemost. 2019;25:1076029619890411. https://doi.org/10.1177/1076029619890411.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Zheng D, Qi G, Adu IK, et al. Efficacy of traditional Chinese medicine combined with rivaroxaban in the treatment of lower extremity deep vein thrombosis: a meta-analysis. Medicine (Baltimore). 2022;101(39):e29483. https://doi.org/10.1097/MD.0000000000029483.

    Article  CAS  PubMed  Google Scholar 

  12. Wang Y, Chen T, Yang C, et al. Huangqi Guizhi Wuwu decoction improves arthritis and pathological damage of heart and lung in TNF-Tg mice. Front Pharmacol. 2022;13:871481. https://doi.org/10.3389/fphar.2022.871481.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Li M, Li Z, Ma X, et al. Huangqi Guizhi Wuwu Decoction can prevent and treat oxaliplatin-induced neuropathic pain by TNFα/IL-1β/IL-6/MAPK/NF-kB pathway. Aging (Albany NY). 2022;14(12):5013–22. https://doi.org/10.18632/aging.203794.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Dai Y. Huang Qi Gui Zhi soup prevention of deep vein thrombosis after knee replacement clinical research [D].Anhui University of Chinese Medicine, 2014.

  15. Zhao Z, Zhao Y, Chen M, et al. Zhao Yuanhong’s differentiation and treatment on lower extremity venous thrombosis due to lung cancer with Huangqi Guizhi Wuwu Decoction: a case report. SHJTCM. 2018;52(12):29–30. https://doi.org/10.16305/j.1007-1334.2018.12.008.

    Article  Google Scholar 

  16. Zhou P. Traditional Chinese medicine. Comb Chem High Throughput Screen. 2010;13(10):836. https://doi.org/10.2174/138620710793360329.

    Article  CAS  PubMed  Google Scholar 

  17. Zhu S, Song Y, Chen X, et al. Traditional Chinese and western medicine for the prevention of deep venous thrombosis after lower extremity orthopedic surgery: a meta-analysis of randomized controlled trials. J Orthop Surg Res. 2018;13(1):79. https://doi.org/10.1186/s13018-018-0785-2.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Niu M, Zhang S, Zhang Bo, et al. Interpretation of network pharmacology evaluation method guidance. Chin Trad Herbal Drugs. 2021;52(14):4119–29. https://doi.org/10.7501/j.issn.0253-2670.2021.14.001.

    Article  Google Scholar 

  19. Hopkins AL. Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. 2008;4(11):682–90. https://doi.org/10.1038/nchembio.118.

    Article  CAS  PubMed  Google Scholar 

  20. Li X, Wei S, Niu S, et al. Network pharmacology prediction and molecular docking-based strategy to explore the potential mechanism of Huanglian Jiedu Decoction against sepsis. Comput Biol Med. 2022;144:105389. https://doi.org/10.1016/j.compbiomed.2022.105389.

    Article  CAS  PubMed  Google Scholar 

  21. Zhang J, Wang R, Liang X, et al. Computation and molecular pharmacology to trace the anti-rheumatoid activity of Angelicae Pubescentis Radix. BMC Complem Med Ther. 2022;22(1):312. https://doi.org/10.1186/s12906-022-03769-w.

    Article  CAS  Google Scholar 

  22. Xu X, Zhang W, Huang C, et al. A novel chemometric method for the prediction of human oral bioavailability. Int J Mol Sci. 2012;13(6):6964–82. https://doi.org/10.3390/ijms13066964.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Tian S, Wang J, Li Y, et al. The application of in silico drug-likeness predictions in pharmaceutical research. Adv Drug Deliv Rev. 2015;86:2–10. https://doi.org/10.1016/j.addr.2015.01.009.

    Article  CAS  PubMed  Google Scholar 

  24. Gao Y, Ji W, Lu M, et al. Systemic pharmacological verification of Guizhi Fuling decoction in treating endometriosis-associated pain. J Ethnopharmacol. 2022;297:115540. https://doi.org/10.1016/j.jep.2022.115540.

    Article  CAS  PubMed  Google Scholar 

  25. Wu Y, You X, Lin Q, et al. Exploring the pharmacological mechanisms of xihuang pills against prostate cancer via integrating network pharmacology and experimental validation in vitro and in vivo. Front Pharmacol. 2022;12:791269. https://doi.org/10.3389/fphar.2021.791269.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Amberger JS, Bocchini CA, Schiettecatte F, et al. OMIM.org: Online Mendelian Inheritance in Man (OMIM®), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2015;43:D789–98. https://doi.org/10.1093/nar/gku1205.

    Article  CAS  PubMed  Google Scholar 

  27. Rebhan M, Chalifa-Caspi V, Prilusky J, et al. GeneCards: a novel functional genomics compendium with automated data mining and query reformulation support. Bioinformatics. 1998;14(8):656–64. https://doi.org/10.1093/bioinformatics/14.8.656.

    Article  CAS  PubMed  Google Scholar 

  28. Szklarczyk D, Gable AL, Lyon D, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13. https://doi.org/10.1093/nar/gky1131.

    Article  CAS  Google Scholar 

  29. Zhou Y, Zhou B, Pache L, et al. Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat Commun. 2019;10(1):1523. https://doi.org/10.1038/s41467-019-09234-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Chen G, Seukep AJ, Guo M, et al. Recent advances in molecular docking for the research and discovery of potential marine drugs. Mar Drugs. 2020;18(11):545. https://doi.org/10.3390/md18110545.

    Article  CAS  PubMed  Google Scholar 

  31. Cosconati S, Forli S, Perryman AL, et al. Virtual screening with autodock: theory and practice. Expert Opin Drug Discov. 2010;5(6):597–607. https://doi.org/10.1517/17460441.2010.484460.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Burley SK, Bhikadiya C, Bi C, et al. RCSB Protein Data Bank: powerful new tools for exploring 3D structures of biological macromolecules for basic and applied research and education in fundamental biology, biomedicine, biotechnology, bioengineering and energy sciences. Nucleic Acids Res. 2021;49(D1):D437–51. https://doi.org/10.1093/nar/gkaa1038.

    Article  CAS  PubMed  Google Scholar 

  33. Seeliger D, de Groot BL. Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J Comput Aided Mol Des. 2010;24(5):417–22. https://doi.org/10.1007/s10822-010-9352-6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinform. 2003;4:2. https://doi.org/10.1186/1471-2105-4-2.

    Article  Google Scholar 

  35. Zhu W, Li Y, Zhao J, et al. The mechanism of triptolide in the treatment of connective tissue disease-related interstitial lung disease based on network pharmacology and molecular docking. Ann Med. 2022;54(1):541–52. https://doi.org/10.1080/07853890.2022.2034931.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. An W, Huang Y, Chen S, et al. Mechanisms of Rhizoma Coptidis against type 2 diabetes mellitus explored by network pharmacology combined with molecular docking and experimental validation. Sci Rep. 2021;11(1):20849. https://doi.org/10.1038/s41598-021-00293-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Quintal Martínez JP, Segura Campos MR. Flavonoids as a therapeutical option for the treatment of thrombotic complications associated with COVID-19. Phytother Res. 2023;37(3):1092–114. https://doi.org/10.1002/ptr.7700.

    Article  CAS  PubMed  Google Scholar 

  38. Guglielmone HA, Nuñez-Montoya SC, Agnese AM, et al. Quercetin 3,7,3’,4’-tetrasulphated isolated from Flaveria bidentis inhibits tissue factor expression in human monocyte. Phytomedicine. 2012;19(12):1068–71. https://doi.org/10.1016/j.phymed.2012.06.013.

    Article  CAS  PubMed  Google Scholar 

  39. Guglielmone HA, Agnese AM, Nuñez-Montoya SC, et al. Antithrombotic “in vivo” effects of quercetin 3,7,3’,4’-tetrasulfate isolated from Flaveria bidentis in an experimental thrombosis model in mice. Thromb Res. 2020;195:190–2. https://doi.org/10.1016/j.thromres.2020.07.040.

    Article  CAS  PubMed  Google Scholar 

  40. Lee SM, Moon J, Chung JH, et al. Effect of quercetin-rich onion peel extracts on arterial thrombosis in rats. Food Chem Toxicol. 2013;57:99–105. https://doi.org/10.1016/j.fct.2013.03.008.

    Article  CAS  PubMed  Google Scholar 

  41. Choi JH, Park SE, Kim SJ, et al. Kaempferol inhibits thrombosis and platelet activation. Biochimie. 2015;115:177–86. https://doi.org/10.1016/j.biochi.2015.06.001.

    Article  CAS  PubMed  Google Scholar 

  42. Gogoi D, Pal A, Chattopadhyay P, et al. First report of plant-derived β-sitosterol with antithrombotic, in vivo anticoagulant, and thrombus-preventing activities in a mouse model. J Nat Prod. 2018;81(11):2521–30. https://doi.org/10.1021/acs.jnatprod.8b00574.

    Article  CAS  PubMed  Google Scholar 

  43. Schultheiß C, Willscher E, Paschold L, et al. The IL-1β, IL-6, and TNF cytokine triad is associated with post-acute sequelae of COVID-19. Cell Rep Med. 2022;3(6):100663. https://doi.org/10.1016/j.xcrm.2022.100663.

    Article  CAS  PubMed Central  Google Scholar 

  44. Ridker PM, Cushman M, Stampfer MJ, et al. Inflammation, aspirin, and the risk of cardiovascular disease in apparently healthy men. N Engl J Med. 1997;336(14):973–9. https://doi.org/10.1056/NEJM199704033361401.

    Article  CAS  PubMed  Google Scholar 

  45. Salemi R, Gattuso G, Tomasello B, et al. Co-occurrence of interleukin-6 receptor Asp358Ala variant and high plasma levels of IL-6: an evidence of IL-6 trans-signaling activation in deep vein thrombosis (DVT) patients. Biomolecules. 2022;12(5):681. https://doi.org/10.3390/biom12050681.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Zhang Y, Zhang Z, Wei R, et al. IL (Interleukin)-6 contributes to deep vein thrombosis and is negatively regulated by miR-338-5p. Arterioscler Thromb Vasc Biol. 2020;40(2):323–34. https://doi.org/10.1161/ATVBAHA.119.313137.

    Article  CAS  PubMed  Google Scholar 

  47. Wang Y, Golden JB, Fritz Y, et al. Interleukin 6 regulates psoriasiform inflammation-associated thrombosis. JCI Insight. 2016;1(20):e89384. https://doi.org/10.1172/jci.insight.89384.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Sun TW, Zhang JY, Li L, et al. Effect atorvastatin on serum tumor necrosis factor alpha and interleukin-1β following acute pulmonary embolism. Exp Lung Res. 2011;37(2):78–81. https://doi.org/10.3109/01902148.2010.514647.

    Article  CAS  PubMed  Google Scholar 

  49. Libby P. Interleukin-1 beta as a target for atherosclerosis therapy: biological basis of CANTOS and beyond. J Am Coll Cardiol. 2017;70(18):2278–89. https://doi.org/10.1016/j.jacc.2017.09.028.

    Article  CAS  PubMed  Google Scholar 

  50. Qiao L, Ma J, Zhang Z, et al. Deficient chaperone-mediated autophagy promotes inflammation and atherosclerosis. Circ Res. 2021;129(12):1141–57. https://doi.org/10.1161/CIRCRESAHA.121.318908.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Ridker PM, Rane M. Interleukin-6 signaling and anti-interleukin-6 therapeutics in cardiovascular disease. Circ Res. 2021;128(11):1728–46. https://doi.org/10.1161/CIRCRESAHA.121.319077.

    Article  CAS  PubMed  Google Scholar 

  52. Reyes-Farias M, Carrasco-Pozo C. The anti-cancer effect of quercetin: molecular implications in cancer metabolism. Int J Mol Sci. 2019;20(13):3177. https://doi.org/10.3390/ijms20133177.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Guo H, Ding H, Tang X, et al. uercetin induces pro-apoptotic autophagy via SIRT1/AMPK signaling pathway in human lung cancer cell lines A549 and H1299 in vitro. Thorac Cancer. 2021;12(9):1415–22. https://doi.org/10.1111/1759-7714.13925.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Ward AB, Mir H, Kapur N, et al. Quercetin inhibits prostate cancer by attenuating cell survival and inhibiting anti-apoptotic pathways. World J Surg Oncol. 2018;16(1):108. https://doi.org/10.1186/s12957-018-1400-z.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Wang X, Wei G, Ding Y, et al. Protein tyrosine phosphatase PTPN22 negatively modulates platelet function and thrombus formation. Blood. 2022;140(9):1038–51. https://doi.org/10.1182/blood.2022015554.

    Article  CAS  PubMed  Google Scholar 

  56. Babur Ö, Melrose AR, Cunliffe JM, et al. Phosphoproteomic quantitation and causal analysis reveal pathways in GPVI/ITAM-mediated platelet activation programs. Blood. 2020;136(20):2346–58. https://doi.org/10.1182/blood.2020005496.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Mozzini C, Garbin U, Fratta Pasini AM, et al. An exploratory look at NETosis in atherosclerosis. Intern Emerg Med. 2017;12(1):13–22. https://doi.org/10.1007/s11739-016-1543-2.

    Article  Google Scholar 

  58. Prandoni P, Bilora F, Marchiori A, et al. An association between atherosclerosis and venous thrombosis. N Engl J Med. 2003;348(15):1435–41. https://doi.org/10.1056/NEJMoa022157.

    Article  PubMed  Google Scholar 

  59. Hong C, Zhu F, Du D, et al. Coronary artery calcification and risk factors for atherosclerosis in patients with venous thromboembolism. Atherosclerosis. 2005;183(1):169–74. https://doi.org/10.1016/j.atherosclerosis.2005.03.047.

    Article  CAS  PubMed  Google Scholar 

  60. Prandoni P, Villalta S, Bagatella P, et al. The clinical course of deep-vein thrombosis: prospective long-term follow-up of 528 symptomatic patients. Haematologica. 1997;82(4):423–8.

    CAS  Google Scholar 

  61. Keller K, Hobohm L, Münzel T, et al. Impact of systemic atherosclerosis on clinical characteristics and short-term outcomes in patients with deep venous thrombosis or thrombophlebitis. Am J Med Sci. 2022;363(3):232–41. https://doi.org/10.1016/j.amjms.2021.09.002.

    Article  PubMed  Google Scholar 

  62. Malek AM, Alper SL, Izumo S. Hemodynamic shear stress and its role in atherosclerosis. JAMA. 1999;282(21):2035–42. https://doi.org/10.1001/jama.282.21.2035.

    Article  CAS  PubMed  Google Scholar 

  63. Moore KJ, Tabas I. Macrophages in the pathogenesis of atherosclerosis. Cell. 2011;145(3):341–55. https://doi.org/10.1016/j.cell.2011.04.005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Libby P. Inflammation in atherosclerosis. Arterioscler Thromb Vasc Biol. 2012;32(9):2045–51. https://doi.org/10.1161/ATVBAHA.108.179705.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  65. Grainge MJ, West J, Card TR. Venous thromboembolism during active disease and remission in inflammatory bowel disease: a cohort study. Lancet. 2010;375(9715):657–63. https://doi.org/10.1016/S0140-6736(09)61963-2.

    Article  PubMed  Google Scholar 

  66. Fox EA, Kahn SR. The relationship between inflammation and venous thrombosis: a systematic review of clinical studies. Thromb Haemost. 2005;94(2):362–5. https://doi.org/10.1160/TH05-04-0266.

    Article  CAS  PubMed  Google Scholar 

  67. Goldhaber SZ, Bounameaux H. Pulmonary embolism and deep vein thrombosis. Lancet. 2012;379(9828):1835–46. https://doi.org/10.1016/S0140-6736(11)61904-1.

    Article  PubMed  Google Scholar 

  68. Keller K, Hobohm L, Munzel T, et al. Impact of symptomatic atherosclerosis in patients with pulmonary embolism. Int J Cardiol. 2019;278:225–31.

    Article  PubMed  Google Scholar 

  69. Huang L. Clinical Study of Huangqi Guizhi Wuwu Decoction in the Treatment of Coronary Heart Disease Stable Angina Pectoris (Yang Deficiency Pulse Obstruction Syndrome). Chin For Med Res. 2020;18(17):15–7. https://doi.org/10.14033/j.cnki.cfmr.2020.17.006.

    Article  Google Scholar 

  70. Bai-rong Hu. Effect of Huangqi Guizhi Wuwu Decoction on patients with cerebral infarction in convalescence. Chin J Mod Drug Appl. 2021;15(22):188–91. https://doi.org/10.14164/j.cnki.cn11-5581/r.2021.22.074.

    Article  Google Scholar 

  71. Jiang YH, Jiang LY, Wang YC, et al. Quercetin Attenuates Atherosclerosis via Modulating Oxidized LDL-Induced Endothelial Cellular Senescence [published correction appears in Front Pharmacol. 2020 May 29;11:772]. Front Pharmacol. 2020;11:512. https://doi.org/10.3389/fphar.2020.00512.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  72. Cao H, Jia Q, Yan L, et al. Quercetin suppresses the progression of atherosclerosis by regulating MST1-mediated autophagy in ox-LDL-Induced RAW264.7 macrophage foam cells. Int J Mol Sci. 2019;20(23):6093. https://doi.org/10.3390/ijms20236093.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  73. Luo Q, Du R, Liu W, et al. PI3K/Akt/mTOR signaling pathway: role in esophageal squamous cell carcinoma, regulatory mechanisms and opportunities for targeted therapy. Front Oncol. 2022;12:852383. https://doi.org/10.3389/fonc.2022.852383.

    Article  PubMed  PubMed Central  Google Scholar 

  74. Bouvard C, Galy-Fauroux I, Grelac F, et al. Low-molecular-weight fucoidan induces endothelial cell migration via the PI3K/AKT pathway and modulates the transcription of genes involved in angiogenesis. Mar Drugs. 2015;13(12):7446–62. https://doi.org/10.3390/md13127075.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Li T, Zhang L, Jin C, et al. Pomegranate flower extract bidirectionally regulates the proliferation, differentiation and apoptosis of 3T3-L1 cells through regulation of PPARγ expression mediated by PI3K-AKT signaling pathway. Biomed Pharmacother. 2020;131:110769. https://doi.org/10.1016/j.biopha.2020.110769.

    Article  CAS  PubMed  Google Scholar 

  76. He J, Zhang N, Zhu Y, et al. MSC spheroids-loaded collagen hydrogels simultaneously promote neuronal differentiation and suppress inflammatory reaction through PI3K-Akt signaling pathway. Biomaterials. 2021;265:120448. https://doi.org/10.1016/j.biomaterials.2020.120448.

    Article  CAS  PubMed  Google Scholar 

  77. Li Y, Xia J, Jiang N, et al. Corin protects H2O2-induced apoptosis through PI3K/AKT and NF-κB pathway in cardiomyocytes. Biomed Pharmacother. 2018;97:594–9. https://doi.org/10.1016/j.biopha.2017.10.090.

    Article  CAS  PubMed  Google Scholar 

  78. Fang J, Ding M, Yang L, et al. PI3K/PTEN/AKT signaling regulates prostate tumor angiogenesis. Cell Signal. 2007;19(12):2487–97. https://doi.org/10.1016/j.cellsig.2007.07.025.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  79. Chen L, Wang J, Wang B, et al. MiR-126 inhibits vascular endothelial cell apoptosis through targeting PI3K/Akt signaling. Ann Hematol. 2016;95(3):365–74. https://doi.org/10.1007/s00277-015-2567-9.

    Article  CAS  PubMed  Google Scholar 

  80. Meng Y, Yin Q, Ma Q, et al. FXII regulates the formation of deep vein thrombosis via the PI3K/AKT signaling pathway in mice. Int J Mol Med. 2021;47(5):87. https://doi.org/10.3892/ijmm.2021.4920.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Sun Y, Liu WZ, Liu T, et al. Signaling pathway of MAPK/ERK in cell proliferation, differentiation, migration, senescence and apoptosis. J Recept Signal Transduct Res. 2015;35(6):600–4. https://doi.org/10.3109/10799893.2015.1030412.

    Article  CAS  PubMed  Google Scholar 

  82. Wagner EF, Nebreda AR. Signal integration by JNK and p38 MAPK pathways in cancer development. Nat Rev Cancer. 2009;9(8):537–49. https://doi.org/10.1038/nrc2694.

    Article  CAS  PubMed  Google Scholar 

  83. Wen C, Ying Y, Zhao H, et al. Resistance exercise affects catheter-related thrombosis in rats through miR-92a-3p, oxidative stress and the MAPK/NF-κB pathway. BMC Cardiovasc Disord. 2021;21(1):440. https://doi.org/10.1186/s12872-021-02233-w.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Irfan M, Jeong D, Kwon HW, et al. Ginsenoside-Rp3 inhibits platelet activation and thrombus formation by regulating MAPK and cyclic nucleotide signaling. Vascul Pharmacol. 2018;109:45–55. https://doi.org/10.1016/j.vph.2018.06.002.

    Article  CAS  PubMed  Google Scholar 

  85. Manne BK, Münzer P, Badolia R, et al. PDK1 governs thromboxane generation and thrombosis in platelets by regulating activation of Raf1 in the MAPK pathway. J Thromb Haemost. 2018;16(6):1211–25. https://doi.org/10.1111/jth.14005.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  86. Cheng SC, Huang WC, Pang JH, et al. Quercetin Inhibits the Production of IL-1β-Induced Inflammatory Cytokines and Chemokines in ARPE-19 Cells via the MAPK and NF-κB Signaling Pathways. Int J Mol Sci. 2019;20(12):2957. https://doi.org/10.3390/ijms20122957.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  87. Cavalli G, Colafrancesco S, Emmi G, et al. Interleukin 1α: a comprehensive review on the role of IL-1α in the pathogenesis and treatment of autoimmune and inflammatory diseases. Autoimmun Rev. 2021;20(3):102763. https://doi.org/10.1016/j.autrev.2021.102763.

    Article  CAS  PubMed  Google Scholar 

  88. Yadav V, Chi L, Zhao R, et al. Ectonucleotidase tri(di)phosphohydrolase-1 (ENTPD-1) disrupts inflammasome/interleukin 1β-driven venous thrombosis. J Clin Invest. 2019;129(7):2872–7. https://doi.org/10.1172/JCI124804.

    Article  PubMed  PubMed Central  Google Scholar 

  89. Rehman HM, Sajjad M, Ali MA, et al. Identification of NS2B-NS3 protease inhibitors for therapeutic application in ZIKV infection: a pharmacophore-based high-throughput virtual screening and MD simulations approaches. Vaccines (Basel). 2023;11(1):131. https://doi.org/10.3390/vaccines11010131.

    Article  CAS  PubMed  Google Scholar 

  90. Hosen MA, El Bakri Y, Rehman HM, et al. A computational investigation of galactopyranoside esters as antimicrobial agents through antiviral, molecular docking, molecular dynamics, pharmacokinetics, and bioactivity prediction [published online ahead of print, 2023 Apr 7]. J Biomol Struct Dyn. 2023. https://doi.org/10.1080/07391102.2023.2198606.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to thank all authors of references.

Funding

This research is supported by This work is supported by the research project of Science and Technology Bureau of Luzhou (2021-JYJ-60).

Author information

Authors and Affiliations

Authors

Contributions

WF and SL contributed equally to this work, they led the conception and design of the research and wrote the article. YY climbed to make critical revisions to important knowledge content, while JL ultimately approved the version to be submitted. All authors approved the final manuscript as submitted and agree to be accountable for all aspects of the work.

Corresponding authors

Correspondence to Yunkang Yang or Jie Liang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file: 1.

Effective compounds of Huang-Qi-Gui-Zhi-Wu-Wu Decoction.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Fan, W., Lan, S., Yang, Y. et al. Network pharmacology prediction and molecular docking-based strategy to discover the potential pharmacological mechanism of Huang–Qi–Gui–Zhi–Wu–Wu decoction against deep vein thrombosis. J Orthop Surg Res 18, 475 (2023). https://doi.org/10.1186/s13018-023-03948-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13018-023-03948-6

Keywords