- Hunan Key Laboratory of Kidney Disease and Blood Purification, Department of Nephrology, The Second Xiangya Hospital of Central South University, Changsha, Hunan, China
Danggui-Shaoyao-San (DSS) is one of traditional Chinese medicine, which recently was found to play a protective role in diabetic kidney disease (DKD). However, the pharmacological mechanisms of DSS remain obscure. This study would explore the molecular mechanisms and bioactive ingredients of DSS in the treatment of DKD through network pharmacology. The potential target genes of DKD were obtained through OMIM database, the DigSee database and the DisGeNET database. DSS-related targets were acquired from the BATMAN-TCM database and the STITCH database. The common targets of DSS and DKD were selected for analysis in the STRING database, and the results were imported into Cytoscape to construct a protein-protein interaction network. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis and Gene Ontology (GO) enrichment analysis were carried out to further explore the mechanisms of DSS in treating DKD. Molecular docking was conducted to identify the potential interactions between the compounds and the hub genes. Finally, 162 therapeutic targets of DKD and 550 target genes of DSS were obtained from our screening process. Among this, 28 common targets were considered potential therapeutic targets of DSS for treating DKD. Hub signaling pathways including HIF-1 signaling pathway, TNF signaling pathway, AMPK signaling pathway, mTOR signaling pathway, and PI3K-Akt signaling pathway may be involved in the treatment of DKD using DSS. Furthermore, TNF and PPARG, and poricoic acid C and stigmasterol were identified as hub genes and main active components in this network, respectively. In this study, DSS appears to treat DKD by multi-targets and multi-pathways such as inflammatory, oxidative stress, autophagy and fibrosis, which provided a novel perspective for further research of DSS for the treatment of DKD.
Introduction
Diabetic kidney disease (DKD) is one of the most common microvascular complications of diabetes and is the leading cause of end-stage renal disease (Tuttle et al., 2014). Pathologically, DKD is characterized by glomerular basement membrane thickening, mesangial expansion, glomerulosclerosis and tubulointerstitial changes (Tervaert et al., 2010). A great number of studies have shown that the pathogenesis of DKD involves multiple mechanisms, including oxidative stress, inflammation, fibrosis and lipid accumulation (Liang et al., 2018; Tanaka et al., 2018; Chen X et al., 2019). Although existing treatments can slow down the progression of DKD in some patients, additional potent treatments need to be explored to fight against DKD.
Traditional Chinese medicine has been widely applied in the treatment of disease for thousands of years, and the application scope of many traditional Chinese medicine has been expanding (Tong et al., 2018; Peng et al., 2020). Danggui-Shaoyao-San (DSS) consists of Paeonia lactiflora Pall., Alisma orientalis (Sam.) Juzep., Angelica sinensis (Oliv.) Diels., Poria cocos (Schw.) Wolf., Atractylodes macrocephala Koidz. and Ligusticum chuanxiong Hort. In 2009, Chen et al. (2009) reported that they identified 41 components in DSS extract and determined 14 bioactive components. In addition, Yin et al. (2021) found that different solvents yield different extracts with different compounds and the extracts from DSS improved metabolic syndrome symptoms and the biochemical factors by regulating intestinal floras and improving hepatic gene expressions and metabolites. Actually, DSS used to be primarily used to treat gynecological diseases (Kotani et al., 1997; Lee et al., 2016). Subsequently, researchers have demonstrated that DSS can be applicable to ameliorate neurological disorders including Alzheimer’s disease (Huang et al., 2020; Kim and Cho, 2020; Li et al., 2021). Fu et al. (2022) recently reviewed that polysaccharide extracted from DSS play significant roles in the treatment of Alzheimer’s disease through various mechanisms, such as anti-inflammatory, anti-neuronal injury, and immune regulation. Furthermore, recent studies have showed that DSS may play an important role in the treatment of diabetes and its complications (Xiaobing et al., 2020; Shi et al., 2021). Especially, Liu et al. (2012) found that DSS has the renal protective effects on STZ-induced diabetic rats. And Yang et al. (2022) screened bioactive compounds from DSS for treating sodium retention in nephrotic syndrome and found that 1,2,3,4,6-O-pentagalloylglucose is a potential bioactive ingredient responsible for the effect of DSS on natriuresis. However, the bioactive components and molecular mechanisms of DSS as related to DKD remained unclear due to the limited research. Therefore, the aim of this study was to explore the potential key pharmacologic mechanisms of DSS on DKD.
Network pharmacology is an emerging discipline that integrate virtual computing technology, public databases, and high-throughput data for revealing the complicated network relationship between drugs and diseases. Nowadays, many researchers have used network pharmacology to study the therapeutic targets of traditional Chinese medicine, laying a foundation for the clinical application of traditional Chinese medicine (Qin et al., 2020; Zhou et al., 2021). DSS, as a traditional Chinese medicine, have multiple active compounds, however, its bioactive components and molecular mechanisms for the treatment of DKD needs to be further explored. In this study, the network pharmacology based on a compound-target-disease interaction network and biologic analysis were used to explore the potential mechanisms of DSS in the treatment of DKD and provide a basis for subsequent pharmacologic experimental research.
Materials and methods
Identifying diabetic kidney disease-Related Targets
Human genes associated with DKD were gathered from the OMIM database (http://omim.org) (Amberger et al., 2015), the DigSee database (http://210.107.182.61/geneSearch/) (Kim et al., 2013; Kim et al., 2017) and the DisGeNET database (https://www.disgenet.org/) (Piñero et al., 2015; Piñero et al., 2020). Terms like (diabetic and nephropathy) or (“diabetic nephropathy”) or (“diabetic kidney disease”) were used as baits to search DKD-related targets from the OMIM database. By removing duplicate items, 37 genes were obtained. The keywords “Diabetic nephropathy” or “Diabetic nephropathies” were used in the DisGeNET database or the DigSee database to screen targets related to DKD. Genes with ‘Score_gda’ > 0.25 in the DisGeNET database or number of evidence articles >5 in the DigSee database were selected for further analysis. Seventy-four DKD-related genes and 91 DKD-related genes were collected from these two databases, respectively. All genes were then inputted into UniProt (https://www.uniprot.org/) and converted into a unified gene name (UniProt, 2019). After deleting duplicates, 162 genes were identified as DKD targets finally.
Active ingredients
DSS prescription consists of the following herbs: Paeonia lactiflora Pall., Alisma orientalis (Sam.) Juzep., Angelica sinensis (Oliv.) Diels., Poria cocos (Schw.) Wolf., Atractylodes macrocephala Koidz. and Ligusticum chuanxiong Hort. Active ingredients of every herb were mainly accessed through the TCMSP database (https://tcmspw.com/tcmsp.php), a platform provides information about Chinese herbal medicine, related targets and ingredients-targets-diseases networks (Ru et al., 2014). Names of herbs were inputted into the search box to obtain the ingredients of DSS. And the ingredients meeting following criterions were identified as pivotal active ingredients: 1) oral bioavailability ≥30%; 2) drug-likeness ≥ 0.18; 3) number of rotatable bonds ≤10; 4) molecular weight: 180-500.
Predicting Danggui-Shaoyao-San-Related Targets
Targets of DSS active ingredients were acquired from the BATMAN-TCM database (http://bionet.ncpsb.org/batman-tcm/) (Liu et al., 2016) and the STITCH database (http://stitch.embl.de/) (Szklarczyk et al., 2016). Results were limited to “Homo sapiens.” Targets predicted from the BATMAN-TCM database were accepted for subsequent analysis under the retrieval of score >10. All obtained targets were standardized through the UniProt (https://www.uniprot.org/) (UniProt, 2019). Finally, 550 targets were predicted as DSS corresponding targets.
Identifying ingredients-disease common targets
The screened DKD targets and DSS active ingredients related targets were imported into TBtools (version 1.0686) (Chen C et al., 2020) and the overlaid targets were identified as potential DKD targets that may be targeted by DSS.
Enrichment analysis and network construction
Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enrichment analysis and Gene Ontology (GO) enrichment analysis were conducted by online function annotation tool DAVID (https://david.ncifcrf.gov/, version 6.8) (Huang da et al., 2009a; Huang da et al., 2009b) and visualized by a data visualization website (http://www.bioinformatics.com.cn/). Protein-protein interaction (PPI) networks for DKD targets, DSS targets and common targets were established by STRING database (https://string-db.org/, version 11.0) (Szklarczyk et al., 2019). PPI network of DKD targets, DSS-active ingredients-targets network, DSS-active ingredients-common targets network and DSS-active ingredients-common targets-enriched KEGG pathways network were further constructed by Cytoscape (version 3.7.2). Cytohubba and MCODE tools were used to screen hub genes.
Molecular docking simulation
Crystal structures of screened hub macromolecule receptors were obtained from PDB database (http://www.pdb.org) (Berman et al., 2002) and were then imported into Pymol software (version 2.4.0) for removal of water molecules, co-crystallized ligands and ions. Subsequently, hydrogens were added and gasteiger charges were computed by AutoDockTools (version 1.5.6). The structures were saved as PDBQT format files after AD4 type was assigned. 3D-structures of corresponding ingredients were downloaded from pubchem database (https://pubchem.ncbi.nlm.nih.gov/) and were converted into Mol2 format files by Open Babel GUI software (version 2.3.1 http://openbabel.org) (O'Boyle et al., 2011). Then the small molecule ligand files were imported into AutoDock Tools (version 1.5.6) and torsions were set automatically. The structures were then saved as PDBQT format files, too. Subsequently, PDBQT files of macromolecule receptors and corresponding small molecule ligands were imported into AutoDockTools (version 1.5.6) for construction of mating pockets. Further molecular docking simulation was conducted in AutoDockTools (version 1.5.6) by using Genetic algorithm. Conformations with the lowest docking binding free energy were selected as the most possible docking modes and visualized using PyMOL (version 2.4.0).
3 Results
Diabetic kidney disease therapeutic targets, enriched pathways and PPI network
Flow chart of the network pharmacology is given in Figure 1. Briefly, DKD therapeutic targets and DSS targets were obtained from online databases, and the shared targets of these two parts were then identified as latent therapeutic targets of DSS against DKD. Function analysis and PPI network construction were then implemented to elucidate the potential mechanism of DSS in the treatment of DKD. Subsequent molecular docking was performed to further clarify the drug-target interactions.
FIGURE 1. Schema of network pharmacology analysis of Danggui-Shaoyao-San (DSS) against diabetic kidney disease (DKD).
An aggregate of 162 targets were screened as DKD therapeutic targets from OMIM, DigSee and DigSee database (Figure 2A). KEGG enrichment analysis of these genes showed that pathways, including HIF-1 signaling pathway, PI3K-AKT signaling pathway, VEGF signaling pathway, TNF signaling pathway and NOD-like receptor signaling pathway, participated in the development of DKD (Figure 2B). A PPI network of these targets were then constructed. One hundred and sixty DKD potential therapeutic targets nodes and 2,465 protein-protein edges were present in Figure 3. According to the number of interactions between points, we sorted these disease targets. The points closing to the center and being redder indicate more numbers of interactions between points and other points and its importance. As presented in Figure 3, targets like INS, ALB, AKT1, IL6, VEGFA and TNF might be key targets of DKD.
FIGURE 2. DKD targets and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. (A) Venn diagram composed of the OMIM database retrieved DKD targets, DigSee database retrieved DKD targets and DisGeNET database retrieved DKD targets. (B) Dot bubble plot of the top 20 KEGG pathways enriched by the targets of DKD. The size of these dots represents the number of enriched genes, and the color represents the −log10 (FDR).
FIGURE 3. Protein-protein interactions network of DKD targets. Color and location of one protein represent the number of connections with other proteins. The redder the color and the more central the position, the more interaction with other proteins. FDR: false discovery rate. DKD: diabetic kidney disease.
Danggui-Shaoyao-San targets, function, and protein-protein interaction network
Twenty-four active components of 6 Chinese herbal medicines were recognized according to the screening conditions. And 550 targets of these active components were further predicted by the BATMAN-TCM database and the STITCH database. KEGG enrichment analysis suggested that these targets of DSS are involved in a variety of metabolic pathways, such as glycine, serine and threonine metabolism, histidine metabolism, and fatty acid degradation (Figure 4A). Besides, PPAR signaling pathway, cGMP-PKG signaling pathway, and cAMP signaling pathway were also been identified (Figure 4A). Other important terms like inflammatory mediator regulation of TRP channels, insulin resistance, apoptosis, HIF-1 signaling pathway, Rap1 signaling pathway, AMPK signaling pathway, and mTOR signaling pathway are also enriched by these targets (data not show). GO enrichment analysis for molecular function indicated that targets of DSS possess protein homodimerization activity, heterodimerization activity and oxidoreductase activity (Figure 4B). Furthermore, these targets were also enriched in enzyme binding, transcription factor binding and other molecular functions (Figure 4B). GO enrichment analysis for cellular component implied that most of these targets were located in cytosol and mitochondrion (Figure 4B). GO enrichment analysis for biological process showed that these active components of DSS were possibly associated with oxidation-reduction process, response to drug, apoptotic process and so on (Figure 4B). Moreover, we constructed a network composed of traditional Chinese herbal medicines of DSS, active ingredients, targets and interactions between targets (Figure 5).
FIGURE 4. KEGG enrichment, Gene Ontology (GO) enrichment of the DSS targets. (A) Dot bubble plot of the top 20 KEGG pathways enriched by the predicted targets of DSS. The size of these dots represents the number of enriched genes, and the color represents -log10 (FDR). (B) Histogram plot of the top 10 enriched molecular function, cellular component and biological process of DSS targets. BP: biological processes; CC: cellular components; MF: molecular functions.
FIGURE 5. Protein-protein interaction network of proteins targeted by 24 active chemical ingredients of DSS. The green rectangles represent herbs of DSS, the purple triangles represent screened active ingredients of these herbs, and the blue circles represent their targets. The red lines represent the interactions between herbs and the chemical ingredients, the grass green lines represent the interactions between ingredients and the predicted targets, while the protein-protein interactions were represented as prey lines.
Danggui-Shaoyao-San-diabetic kidney disease common targets
As presented in Figure 6A, 28common targets were identified as DKD related proteins targeted by DSS. Subsequently, we constructed a network consisting of DSS herbs, active ingredients, common targets, interactions between targets and disease (Figure 6B). Based on the PPI network of common targets, 16 core target genes including HIF1A, TNF, IL1B, IL10, ICAM1, NOTCH1, INS, AKT1, AGTR1, SERPINE1, PTGS2, NOS3, IGF1, ADIPOQ, PPARG, and APOE were recognized by Cytoscape MCODE tool (Figure 6C; Table 1). In order to further clarify the function of these common genes, KEGG pathway enrichment analysis was carried out. As shown in Figure 7A, DSS may ameliorate renal damage of DKD through HIF-1 signaling pathway, TNF signaling pathway, mTOR signaling pathway and so on. GO enrichment analysis indicated that these common targets may be associated with nitric oxide biosynthetic, inflammatory response, NF-kappa B transcription factor activity and other important biological processes (Figure 7B). Besides, these common targets located in extracellular space, plasma membrane and cytosol most (Figure 7B). And GO enrichment analysis for molecular function showed that these targets may be related to hormone activity, protease binding, insulin receptor binding and so on (Figure 7B). To better understand the relationship between these drug-disease common targets and KEGG pathways, we constructed a network composed of herbs, chemical ingredients, common targets and KEGG pathways Figure 8. And the relationships among these chemical ingredients, targets and KEGG pathways were also summarized in Supplemenatry Table S1.
FIGURE 6. Matching of target genes between DKD and DSS, and network analysis of the common targets. (A) Venn diagram composed of the ingredients targets and DKD targets. (B) Herbs-active ingredients-common targets-Disease network. The green rectangles represent herbs of DSS, the purple triangles represent screened active ingredients of these herbs, and the blue circles represent the common targets of DSS and DKD. The green lines represent the interactions between herbs and the chemical ingredients, the orange lines represent the interactions between ingredients and the common targets, the yellow lines represent interactions between DKD and the common targets, while the protein-protein interactions were represented as prey lines. (C) Protein-protein interaction network constructed of 16 core targets screened by MCODE tools.
FIGURE 7. Functional enrichment analysis of common targets. (A) Dot bubble plot of the top 20 KEGG pathways enriched by the common targets of DKD and DSS. Size of the dots represents the number of enriched genes, and the color represents −log10 (p value). (B) Histogram plot of the top 10 biological processes, cellular components and molecular functions of common targets by GO enrichment analysis. BP: biological processes; CC: cellular components; MF: molecular functions.
FIGURE 8. Herbs-active ingredients-common targets-enriched KEGG pathways network. Green hexagons represent herbs of DSS, the purple triangles represent screened active ingredients of these herbs, the red circles represent the 16 core common targets of DSS and DKD, pink circles represent other common targets, and the orange rectangles represent enriched KEGG pathways. The green lines represent the interactions between herbs and the chemical ingredients, the blue lines represent the interactions between ingredients and the common targets, the red lines represent interactions between the common targets and KEGG pathways, while the protein-protein interactions were represented as prey lines.
Pivotal intersection targets and molecular docking
Subsequently, the pivotal intersection targets were screened by Cytohubba tool in Cytoscape software. As showed in Figure 9A and Table 1, TNF, PPARG, ATK1 and IL1B were identified as hub genes in this network, and may play a critical role in the therapeutic effect of DSS against DKD, especially TNF and PPARG. Molecular docking of these 2 proteins and their predicted interacting DSS ingredients were performed then. The binding energy and binding activity of these interacting pairs were summarized in Table 2. The interaction between poricoic acid C and TNF requires the lowest docking binding free energy (−6.29 Kcal/Mol), followed by poricoic acid C and PPARG (−5.53 Kcal/Mol), and stigmasterol and TNF (−5.46 Kcal/Mol), indicating a good binding activity in the molecular docking, while the interaction between folic acid and PPARG was not certain (-2.6 Kcal/Mol). The interaction mode of stigmasterol and TNF was displayed in Figures 9B, C. As shown in Figure 9B, stigmasterol was positioned in the small groove on the surface of TNF. And Figure 9C suggested that stigmasterol may interact with TNF through ALA-35 residue and LEU-37 residue. Stigmastrol may be pulled to TNF through 2 hydrogen bonds and thus exerts the therapeutic potential (Figure 9C). Similarly, Figures 9D–I showed the molecular docking results of poricoic acid C with TNF, eburicoic acid with PPARG, and poricoic acid C with PPARG.
FIGURE 9. Screening of pivotal common targets and molecular docking of these targets binding to DSS ingredients. (A) Protein-protein interaction network of 28 common targets, ranked by Cutohubba tool. The redder the color, the more likely it is to be pivotal targets. (B,C) Molecular docking of stigmasterol with TNF. (D,E) Molecular docking of poricoic acid C with TNF. (F,G) Molecular docking of eburicoic acid with PPARG. (H,I) Molecular docking of poricoic acid C with PPARG. The blue small molecular structures represent the ingredients. And the white macromolecules represent screened important macromolecular receptors in panel (B,D,F,H). The green structures in panel (C,E,G,I) represent screened macromolecular receptor structures, while the red sticks represent the interacted residuals, the yellow dotted line represent intermolecular hydrogen bonding.
Discussion
Although there are many ways to slow the progression of DKD, there is still a lack of an effective way to prevent or cure it due to its complicate pathogenic mechanism. Therefore, it is urgent to find new and multi-effective drugs for the treatment of DKD from a broader perspective. DSS is one of promising drugs for DKD. Recent studies demonstrated that DSS plays a crucially reno-protective role in STZ-induced diabetic rats (Liu et al., 2012; Xiaobing et al., 2020). Unfortunately, the underlying mechanisms of DSS against DKD is still unknown. In this study, we investigate the potential molecular mechanism and active ingredients of DSS in the treatment of DKD based on network pharmacology.
In the current study, 28 DSS-DKD-common-targets were obtained and 16 core targets were acquired based on PPI network of these common targets by MCODE. The core targets include HIF1A, TNF, IL1B, IL10, ICAM1, NOTCH1, INS, AKT1, AGTR1, SERPINE1, PTGS2, NOS3, IGF1, ADIPOQ, PPARG, and APOE. Besides, KEGG pathways enrichment analysis of these 28 common targets found that DSS might ameliorate diabetic kidney injury through HIF-1 signaling pathway, TNF signaling pathway, mTOR signaling pathway and so on. Subsequent analysis by Cytohubba further ranked the DSS-DKD-common-targets. TNF and PPARG were identified as the most core targets.
TNF is a pro-inflammatory cytokine that promotes the expression of a variety of transcription factors, cytokines, cell adhesion molecules, mediators of inflammatory processes, and acute-phase proteins by binding to its receptors (Ortiz et al., 1995; Navarro and Mora-Fernández, 2006). TNF was also reported to be involved in the development and progression of DKD (Navarro-González and Mora-Fernández, 2008; Park et al., 2019). And TNF is able to induce diabetic renal injury by inducing apoptosis, the production of reactive oxygen species and so on (Laster et al., 1988; Radeke et al., 1990; Boyle et al., 2003). Other inflammatory factors like IL1B and IL10, as well as acute-phase protein marker of inflammation, ICAM1, were also found to be critical core targets of DSS against DKD. These targets have been proved to play an important role in the occurrence and development of DKD (Sangoi et al., 2016; Tang and Yiu, 2020). These proteins also participant in inflammation-related signaling pathways such as TNF signaling pathway and T cell receptor signaling pathway. Yin et al. (2021) also reported that DSS could reduce the expression of TNF in fructose-fed rats and Bi et al. (2021) demonstrated that DSS alleviates atherosclerosis through decreasing the expression of molecule including TNF, IL-1B and ICAM1. Hence, in this study, we found that DSS may improve DKD through inflammation-related signaling pathways by targeting TNF and other related targets, which was consistent with previous research.
PPARG is another important DKD-related gene. Previous study identified PPARG as a transcription factor modulating DKD target genes (Chen L et al., 2020). And Sharaf et al. (2018) found that PPARG showed lower expression in DKD patients and the abundance of PPARG were negatively correlated with microalbumin. Nevertheless, few studies have explored the specific role of PPARG in DKD. PPARG, a type II nuclear receptor, participates in fatty acid storage and glucose metabolism (Kökény et al., 2021). Recently, more and more researchers pay attention to the important role of abnormal lipid metabolism and renal ectopic fat deposition in the pathogenesis of DKD (Chen X et al., 2019; Yoshioka et al., 2021). In this study, we reported that DSS regulates a variety of metabolic pathways including PPAR signaling pathway. Consistently, Liu et al. (2022) also demonstrated that DSS regulates lipid metabolism though PPARG in scopolamine-induced amnesia. These findings indicated that DSS may play a reno-protective role in DKD through metabolism pathways by targeting PPARG, which is worthy to be validated by experiments in the future.
Further molecular docking of TNF and PPARG with their predicted active ingredients of DSS was conducted. Among them, the docking of poricoic acid C and TNF had the lowest binding energy (−6.29 Kcal/Mol), followed by poricoic acid C and PPARG (−5.53 Kcal/Mol), and stigmasterol and TNF (−5.46 Kcal/Mol), mirroring a good binding activity in the molecular docking. Poricoic acid C is a kind of triterpenoid compounds. Previous study found that poricoic acid A could protect against acute kidney injury-to-chronic kidney disease transition in rats and attenuate fibroblast activation and abnormal extracellular matrix remodeling in renal fibrosis (Chen DQ et al., 2019; Chen DQ et al., 2020). Poricoic acid B could decrease the production of cytokines (TNF, IL-1B, and IL-6) (Zhang et al., 2021). And our data showed that poricoic acid C might interact with TNF and thus ameliorate renal injury of DKD. On the other hand, Chen et al. (2009) demonstrated that poricoic acid could significantly decrease free fatty acid-induced intracellular triglyceride accumulation (Kim et al., 2019). And consistently, we predicted that poricoic acid C could interact with PPARG. Meanwhile, Stigmasterol is a natural with anti-inflammatory activity. Khan et al. (2020) found that stigmasterol significantly suppresses the inflammation by decreasing the expression of proinflammatory mediators (TNF, IL-6, IL-1B) and increasing the expression of anti-inflammatory cytokine (IL-10) in collagen-induced arthritis. And Kishore demonstrated that stigmasterol might inhibit the progression of DKD by ameliorating oxidative stress (Kishore et al., 2016). In general, our analysis shows that poricoic acid C and stigmasterol may play a protective role in DKD through anti-inflammatory and anti-lipid deposition via TNF and PPARG. However, there is no research on the role of poricoic acid C in DKD at present. Subsequent experiments can explore the reno-protective function of these active ingredients.
We also explore the potential mechanisms of DSS against DKD through functional enrichment analysis. KEGG pathways enrichment analysis suggests that DSS participates in the regulation of several important pathways, such as HIF-1 signaling pathway, TNF signaling pathway, mTOR signaling pathway, AMPK signaling pathway, NF-kappa B signaling pathway and so on, some of which were also key pathways of DKD. These pathways are closely related to hypoxia injury, inflammation and autophagy. Previous studies also reported that DSS alleviated kidney injury through its anti-oxidant effect (Liu et al., 2012). In cerebral ischemic-reperfusion injury, Luo et al. (2021) also found that DSS could attenuate oxidative stress and inhibit neuronal apoptosis. Besides, in lipopolysaccharide-induced human umbilical vein endothelial cells injury model, Bi et al. (2021) suggested that DSS could inhibit the phosphorylation of NF-kappa B. In addition, Zhang and others demonstrated that DSS improves renal fibrosis through attenuating hypoxia and regulating autophagy in unilateral ureteral obstruction rat model (Zhang et al., 2019). The above results are consistent with the pathways of DSS against DKD, strongly indicating that DSS is a promising drug for DKD treatment through regulating multiple cellular biological processes, such as inflammatory, oxidative stress, autophagy and fibrosis.
Conclusion
The present study provides new insight into the pharmacology of DSS and the therapeutic mechanism for DKD, which also indicates that network pharmacology could be an effective way to explore the potential mechanism of intricate prescription for diseases. Our results using network pharmacology revealed that DSS play critically protective role in the treatment of DKD through multi-targets and multi-pathways (Figure 10).
FIGURE 10. Schematic diagram of the network pharmacology for the pharmacological mechanisms of DSS for treatment of DKD.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
All authors contributed to the article and approved the submitted version. JY and CL participated in data collection and analysis. YL, YH, HZ, SL, CZ, NJ, and MY participated in manuscript writing. LS is responsible for the integrity and accuracy of the data.
Funding
This work was supported by grants from the National Natural Science Foundation of China (No. 81730018) and the National Key R&D Program of China (No. 2018YFC1314002).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.832299/full#supplementary-material
References
Ahmad Khan, M., Sarwar, A., Rahat, R., Ahmed, R. S., and Umar, S. (2020). Stigmasterol protects rats from collagen induced arthritis by inhibiting proinflammatory cytokines. Int. Immunopharmacol. 85, 106642. doi:10.1016/j.intimp.2020.106642
Amberger, J. S., Bocchini, C. A., Schiettecatte, F., Scott, A. F., and Hamosh, A. (2015). OMIM.org: online Mendelian Inheritance in Man (OMIM®), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 43, D789–D798. doi:10.1093/nar/gku1205
Berman, H. M., Battistuz, T., Bhat, T. N., Bluhm, W. F., Bourne, P. E., Burkhardt, K., et al. (2002). The protein data bank. Acta Crystallogr. D. Biol. Crystallogr. 58, 899–907. doi:10.1107/s0907444902003451
Bi, Y., Han, X., Lai, Y., Fu, Y., Li, K., Zhang, W., et al. (2021). Systems pharmacological study based on UHPLC-Q-Orbitrap-HRMS, network pharmacology and experimental validation to explore the potential mechanisms of Danggui-Shaoyao-San against atherosclerosis. J. Ethnopharmacol. 278, 114278. doi:10.1016/j.jep.2021.114278
Boyle, J. J., Weissberg, P. L., and Bennett, M. R. (2003). Tumor necrosis factor-alpha promotes macrophage-induced vascular smooth muscle cell apoptosis by direct and autocrine mechanisms. Arterioscler. Thromb. Vasc. Biol. 23, 1553–1558. doi:10.1161/01.Atv.0000086961.44581.B7
Chen, C., Chen, H., Zhang, Y., Thomas, H. R., Frank, M. H., He, Y., et al. (2020). TBtools: an integrative toolkit developed for interactive analyses of big biological data. Mol. Plant 13, 1194–1202. doi:10.1016/j.molp.2020.06.009
Chen, D. Q., Feng, Y. L., Chen, L., Liu, J. R., Wang, M., Vaziri, N. D., et al. (2019). Poricoic acid A enhances melatonin inhibition of AKI-to-CKD transition by regulating Gas6/AxlNFκB/Nrf2 axis. Free Radic. Biol. Med. 134, 484–497. doi:10.1016/j.freeradbiomed.2019.01.046
Chen, D. Q., Wang, Y. N., Vaziri, N. D., Chen, L., Hu, H. H., and Zhao, Y. Y. (2020). Poricoic acid A activates AMPK to attenuate fibroblast activation and abnormal extracellular matrix remodelling in renal fibrosis. Phytomedicine 72, 153232. doi:10.1016/j.phymed.2020.153232
Chen, L., Qi, J., Chang, Y. X., Zhu, D., and Yu, B. (2009). Identification and determination of the major constituents in Traditional Chinese Medicinal formula Danggui-Shaoyao-San by HPLC-DAD-ESI-MS/MS. J. Pharm. Biomed. Anal. 50, 127–137. doi:10.1016/j.jpba.2009.03.039
Chen, L., Wu, B., Wang, S., Xiong, Y., Zhou, B., Cheng, X., et al. (2020). Identification of cooperative gene regulation among transcription factors, LncRNAs, and MicroRNAs in diabetic nephropathy progression. Front. Genet. 11, 1008. doi:10.3389/fgene.2020.01008
Chen, X., Han, Y., Gao, P., Yang, M., Xiao, L., Xiong, X., et al. (2019). Disulfide-bond A oxidoreductase-like protein protects against ectopic fat deposition and lipid-related kidney damage in diabetic nephropathy. Kidney Int. 95, 880–895. doi:10.1016/j.kint.2018.10.038
Fu, X., Liu, Q., Sun, X., Chang, H., Liu, Y., and Han, J. (2022). Research advances in the treatment of alzheimer's disease with polysaccharides of danggui-shaoyao-san. J. Alzheimers Dis. 85, 7–19. doi:10.3233/jad-210656
Huang da, W., Sherman, B. T., and Lempicki, R. A. (2009a). Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 37, 1–13. doi:10.1093/nar/gkn923
Huang da, W., Sherman, B. T., and Lempicki, R. A. (2009b). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi:10.1038/nprot.2008.211
Huang, J., Wang, X., Xie, L., Wu, M., Zhao, W., Zhang, Y., et al. (2020). Extract of Danggui-Shaoyao-San ameliorates cognition deficits by regulating DHA metabolism in APP/PS1 mice. J. Ethnopharmacol. 253, 112673. doi:10.1016/j.jep.2020.112673
Kim, Y., and Cho, S. H. (2020). Danggui-Shaoyao-San for dementia: a PRISMA-compliant systematic review and meta-analysis. Medicine 99, e18507. doi:10.1097/md.0000000000018507
Kim, J., So, S., Lee, H. J., Park, J. C., Kim, J. J., and Lee, H. (2013). DigSee: Disease gene search engine with evidence sentences (version cancer). Nucleic Acids Res. 41, W510–W517. doi:10.1093/nar/gkt531
Kim, J., Kim, J. J., and Lee, H. (2017). An analysis of disease-gene relationship from Medline abstracts by DigSee. Sci. Rep. 7, 40154. doi:10.1038/srep40154
Kim, J. H., Sim, H. A., Jung, D. Y., Lim, E. Y., Kim, Y. T., Kim, B. J., et al. (2019). Poria cocus Wolf extract ameliorates hepatic steatosis through regulation of lipid metabolism, inhibition of ER stress, and activation of autophagy via AMPK activation. Int. J. Mol. Sci. 20, E4801. doi:10.3390/ijms20194801
Kishore, L., Kaur, N., and Singh, R. (2016). Renoprotective effect of Bacopa monnieri via inhibition of advanced glycation end products and oxidative stress in STZ-nicotinamide-induced diabetic nephropathy. Ren. Fail. 38, 1528–1544. doi:10.1080/0886022x.2016.1227920
Kökény, G., Calvier, L., and Hansmann, G. (2021). PPARγ and TGFβ-major Regulators of metabolism, inflammation, and fibrosis in the lungs and kidneys. Int. J. Mol. Sci. 22, 10431. doi:10.3390/ijms221910431
Kotani, N., Oyama, T., Sakai, I., Hashimoto, H., Muraoka, M., Ogawa, Y., et al. (1997). Analgesic effect of a herbal medicine for treatment of primary dysmenorrhea--a double-blind study. Am. J. Chin. Med. 25, 205–212. doi:10.1142/s0192415x9700024x
Laster, S. M., Wood, J. G., and Gooding, L. R. (1988). Tumor necrosis factor can induce both apoptic and necrotic forms of cell lysis. J. Immunol. 141, 2629–2634.
Lee, H. W., Jun, J. H., Kil, K. J., Ko, B. S., Lee, C. H., and Lee, M. S. (2016). Herbal medicine (danggui shaoyao san) for treating primary dysmenorrhea: a systematic review and meta-analysis of randomized controlled trials. Maturitas 85, 19–26. doi:10.1016/j.maturitas.2015.11.013
Li, H., Li, D., Zhao, G., Gao, Y., and Ke, J. (2021). Effects of Danggui-Shaoyao-San on depression- and anxiety-like behaviors of rats induced by experimental tooth movement. J. Orofac. Orthop. 83, 23–33. doi:10.1007/s00056-021-00323-0
Liang, G., Song, L., Chen, Z., Qian, Y., Xie, J., Zhao, L., et al. (2018). Fibroblast growth factor 1 ameliorates diabetic nephropathy by an anti-inflammatory mechanism. Kidney Int. 93, 95–109. doi:10.1016/j.kint.2017.05.013
Liu, I. M., Tzeng, T. F., Liou, S. S., and Chang, C. J. (2012). Beneficial effect of traditional chinese medicinal formula danggui-shaoyao-san on advanced glycation end-product-mediated renal injury in streptozotocin-diabetic rats. Evidence Based complementary Altern. Med. 2012, 140103. doi:10.1155/2012/140103
Liu, Z., Guo, F., Wang, Y., Li, C., Zhang, X., Li, H., et al. (2016). BATMAN-TCM: a bioinformatics analysis tool for molecular mechANism of traditional Chinese medicine. Sci. Rep. 6, 21146. doi:10.1038/srep21146
Liu, P., Zhou, X., Zhang, H., Wang, R., Wu, X., Jian, W., et al. (2022). Danggui-Shaoyao-San attenuates cognitive impairment via the microbiota-gut-brain Axis with regulation of lipid metabolism in scopolamine-induced amnesia. Front. Immunol. 13, 796542. doi:10.3389/fimmu.2022.796542
Luo, Y., Chen, H., Tsoi, B., Wang, Q., and Shen, J. (2021). Danggui-Shaoyao-San (DSS) ameliorates cerebral ischemia-reperfusion injury via activating SIRT1 signaling and inhibiting NADPH oxidases. Front. Pharmacol. 12, 653795. doi:10.3389/fphar.2021.653795
Navarro, J. F., and Mora-Fernández, C. (2006). The role of TNF-alpha in diabetic nephropathy: pathogenic and therapeutic implications. Cytokine Growth Factor Rev. 17, 441–450. doi:10.1016/j.cytogfr.2006.09.011
Navarro-González, J. F., and Mora-Fernández, C. (2008). The role of inflammatory cytokines in diabetic nephropathy. J. Am. Soc. Nephrol. 19, 433–442. doi:10.1681/asn.2007091048
O'Boyle, N. M., Banck, M., James, C. A., Morley, C., Vandermeersch, T., and Hutchison, G. R. (2011). Open Babel: an open chemical toolbox. J. Cheminform. 3, 33. doi:10.1186/1758-2946-3-33
Ortiz, A., Bustos, C., Alonso, J., Alcázar, R., López-Armada, M. J., Plaza, J. J., et al. (1995). Involvement of tumor necrosis factor-alpha in the pathogenesis of experimental and human glomerulonephritis. Adv. Nephrol. Necker Hosp. 24, 53–77.
Park, J., Guan, Y., Sheng, X., Gluck, C., Seasock, M. J., Hakimi, A. A., et al. (2019). Functional methylome analysis of human diabetic kidney disease. JCI insight 4, 128886. doi:10.1172/jci.insight.128886
Peng, Y., Ren, D., Song, Y., Hu, Y., Wu, L., Wang, Q., et al. (2020). Effects of a combined fucoidan and traditional Chinese medicine formula on hyperglycaemia and diabetic nephropathy in a type II diabetes mellitus rat model. Int. J. Biol. Macromol. 147, 408–419. doi:10.1016/j.ijbiomac.2019.12.201
Piñero, J., Queralt-Rosinach, N., Bravo, À., Deu-Pons, J., Bauer-Mehren, A., Baron, M., et al. (2015). DisGeNET: a discovery platform for the dynamical exploration of human diseases and their genes. Database 2015, bav028. doi:10.1093/database/bav028
Piñero, J., Ramírez-Anguita, J. M., Saüch-Pitarch, J., Ronzano, F., Centeno, E., Sanz, F., et al. (2020). The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 48, D845–d855. doi:10.1093/nar/gkz1021
Qin, T., Wu, L., Hua, Q., Song, Z., Pan, Y., and Liu, T. (2020). Prediction of the mechanisms of action of shenkang in chronic kidney disease: a network pharmacology study and experimental validation. J. Ethnopharmacol. 246, 112128. doi:10.1016/j.jep.2019.112128
Radeke, H. H., Meier, B., Topley, N., Flöge, J., Habermehl, G. G., and Resch, K. (1990). Interleukin 1-alpha and tumor necrosis factor-alpha induce oxygen radical production in mesangial cells. Kidney Int. 37, 767–775. doi:10.1038/ki.1990.44
Ru, J., Li, P., Wang, J., Zhou, W., Li, B., Huang, C., et al. (2014). TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform. 6, 13. doi:10.1186/1758-2946-6-13
Sangoi, M. B., de Carvalho, J. A., Tatsch, E., Hausen, B. S., Bollick, Y. S., Londero, S. W., et al. (2016). Urinary inflammatory cytokines as indicators of kidney damage in type 2 diabetic patients. Clin. Chim. Acta. 460, 178–183. doi:10.1016/j.cca.2016.06.028
Sharaf, S. A., Kantoush, N. A., Ayoub, D. F., Ibrahim, A. A., Abdelaal, A. A., Aziz, R. A., et al. (2018). Altered expression of WFS1 and NOTCH2 genes associated with diabetic nephropathy in T2DM patients. Diabetes Res. Clin. Pract. 140, 304–313. doi:10.1016/j.diabres.2018.03.053
Shi, J. J., Liu, H. F., Hu, T., Gao, X., Zhang, Y. B., Li, W. R., et al. (2021). Danggui-Shaoyao-San improves cognitive impairment through inhibiting O-GlcNAc-modification of estrogen α receptor in female db/db mice. J. Ethnopharmacol. 281, 114562. doi:10.1016/j.jep.2021.114562
Szklarczyk, D., Santos, A., von Mering, C., Jensen, L. J., Bork, P., and Kuhn, M. (2016). STITCH 5: augmenting protein-chemical interaction networks with tissue and affinity data. Nucleic Acids Res. 44, D380–D384. doi:10.1093/nar/gkv1277
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–d613. doi:10.1093/nar/gky1131
Tanaka, S., Sugiura, Y., Saito, H., Sugahara, M., Higashijima, Y., Yamaguchi, J., et al. (2018). Sodium-glucose cotransporter 2 inhibition normalizes glucose metabolism and suppresses oxidative stress in the kidneys of diabetic mice. Kidney Int. 94, 912–925. doi:10.1016/j.kint.2018.04.025
Tang, S. C. W., and Yiu, W. H. (2020). Innate immunity in diabetic kidney disease. Nat. Rev. Nephrol. 16, 206–222. doi:10.1038/s41581-019-0234-4
Tervaert, T. W., Mooyaart, A. L., Amann, K., Cohen, A. H., Cook, H. T., Drachenberg, C. B., et al. (2010). Pathologic classification of diabetic nephropathy. J. Am. Soc. Nephrol. 21, 556–563. doi:10.1681/asn.2010010010
Tong, X., Xu, J., Lian, F., Yu, X., Zhao, Y., Xu, L., et al. (2018). Structural alteration of gut microbiota during the amelioration of human type 2 diabetes with hyperlipidemia by metformin and a traditional Chinese herbal formula: a multicenter, randomized, open label clinical trial. mBio 9, e02392. doi:10.1128/mBio.02392-17
Tuttle, K. R., Bakris, G. L., Bilous, R. W., Chiang, J. L., de Boer, I. H., Goldstein-Fuchs, J., et al. (2014). Diabetic kidney disease: a report from an ADA consensus conference. Am. J. Kidney Dis. 64, 510–533. doi:10.1053/j.ajkd.2014.08.001
UniProt (2019). UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 47, D506–d515. doi:10.1093/nar/gky1049
Xiaobing, L., Chunling, N., Wenyu, C., Yan, C., and Zhenzhen, L. (2020). Effect of danggui-shaoyao-san-containing serum on the renal tubular epithelial-mesenchymal transition of diabetic nephropathy. Curr. Pharm. Biotechnol. 21, 1204–1212. doi:10.2174/1389201021666200416094318
Yang, M., Ni, L., Wang, Y., Xuan, Z., Wu, H., Zhan, W., et al. (2022). Screening bioactive compounds from Danggui-shaoyao-san for treating sodium retention in nephrotic syndrome using bio-affinity ultrafiltration. J. Ethnopharmacol. 292, 115171. doi:10.1016/j.jep.2022.115171
Yin, J., Lu, J., Lei, P., He, M., Huang, S., Lv, J., et al. (2021). Danggui-Shaoyao-San improves gut microbia dysbiosis and hepatic lipid homeostasis in fructose-fed rats. Front. Pharmacol. 12, 671708. doi:10.3389/fphar.2021.671708
Yoshioka, K., Hirakawa, Y., Kurano, M., Ube, Y., Ono, Y., Kojima, K., et al. (2021). Lysophosphatidylcholine mediates fast decline in kidney function in diabetic kidney disease. Kidney Int. 101, 510–526. doi:10.1016/j.kint.2021.10.039
Zhang, M. Y., Chen, H. H., Tian, J., Chen, H. J., Zhu, L. L., Zhao, P., et al. (2019). Danggui shaoyao san ameliorates renal fibrosis via regulation of hypoxia and autophagy. Evidence-based complementary Altern. Med. 2019, 2985270. doi:10.1155/2019/2985270
Zhang, L., Yin, M., Feng, X., Ibrahim, S. A., Liu, Y., and Huang, W. (2021). Anti-inflammatory activity of four triterpenoids isolated from poriae cutis. Foods (Basel, Switz. 10, 3155. doi:10.3390/foods10123155
Zhou, X. F., Zhou, W. E., Liu, W. J., Luo, M. J., Wu, X. Q., Wang, Y., et al. (2021). A network pharmacology approach to explore the mechanism of HuangZhi YiShen capsule for treatment of diabetic kidney disease. J. Transl. Int. Med. 9, 98–113. doi:10.2478/jtim-2021-0020
Glossary
ADIPOQ adiponectin
AGTR1 angiotensin II receptor type 1
AKT protein kinase B
AKT1 AKT serine/threonine kinase 1
ALB albumin
AMPK adenosine monophosphate-activated protein kinase
APOE apolipoprotein E
cAMP cyclic adenosine monophosphate
cGMP cyclic guanosine monophosphate
DKD diabetic kidney disease
DSS Danggui-Shaoyao-San
GO Gene Ontology
HIF-1 hypoxia-inducible factor 1
HIF1A hypoxia inducible factor 1 subunit alpha
ICAM1 intercellular adhesion molecule 1
IGF1 insulin like growth factor 1
IL10 interleukin-10
IL1B interleukin-1 beta
IL6 interleukin-6
INS insulin
KEGG Kyoto Encyclopedia of Genes and Genomes
mTOR mechanistic target of rapamycin kinase
NF-kappa B nuclear factor-kappa B
NOD-like receptor nucleotide-binding oligomerization domain-like receptors
NOS3 nitric oxide synthase 3
NOTCH1 notch receptor 1
PI3K phosphatidylinositol 3-kinase
PKG protein kinase G
PPAR peroxisome proliferator-activated receptor
PPARG peroxisome proliferator-activated receptor gamma
PPI Protein-protein interaction
PTGS2 prostaglandin-endoperoxide synthase 2
Rap1 Ras-related protein 1
SERPINE1 serpin family E member 1
TNF tumor necrosis factor
VEGF vascular endothelial growth factor
VEGFA vascular endothelial growth factor A
Keywords: Danggui-Shaoyao-San, diabetic kidney disease, network pharmacology, inflammation, oxidative stress
Citation: Yang J, Li C, Liu Y, Han Y, Zhao H, Luo S, Zhao C, Jiang N, Yang M and Sun L (2022) Using network pharmacology to explore the mechanism of Danggui-Shaoyao-San in the treatment of diabetic kidney disease. Front. Pharmacol. 13:832299. doi: 10.3389/fphar.2022.832299
Received: 12 January 2022; Accepted: 29 July 2022;
Published: 19 August 2022.
Edited by:
Quan Hong, Chinese PLA General Hospital, ChinaCopyright © 2022 Yang, Li, Liu, Han, Zhao, Luo, Zhao, Jiang, Yang and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lin Sun, sunlin@csu.edu.cn
†These authors have contributed equally to this work