Network pharmacology and molecular docking technology-based predictive study of the active ingredients and potential targets of rhubarb for the treatment of diabetic nephropathy
BMC Complementary Medicine and Therapies volume 22, Article number: 210 (2022)
Diabetic nephropathy (DN) is one of the most serious complications of diabetes and the main cause of end-stage renal failure. Rhubarb is a widely used traditional Chinese herb, and it has exhibited efficacy in reducing proteinuria, lowering blood sugar levels and improving kidney function in patients with DN. However, the exact pharmacological mechanism by rhubarb improves DN remain unclear due to the complexity of its ingredients. Hence, we systematically explored the underlying mechanisms of rhubarb in the treatment of DN. We adopted a network pharmacology approach, focusing on the identification of active ingredients, drug target prediction, gene collection, Gene Ontology enrichment and Kyoto Encyclopedia of Genes and Genomes enrichment. Molecular docking technology was used to verify the binding ability between the main active compounds and central therapeutic targets, and screen out the core active ingredients in rhubarb for the treatment of DN. Finally, molecular dynamics simulation was performed for the optimal core protein-ligand obtained by molecular docking using GROMACS software. The network analysis identified 16 active compounds in rhubarb that were linked to 37 possible therapeutic targets related to DN. Through protein–protein interaction analysis, TP53, CASP8, CASP3, MYC, JUN and PTGS2 were identified as the key therapeutic targets. By validation of molecular docking, finding that the central therapeutic targets have good affinities with the main active compounds of rhubarb, and rhein, beta-sitosterol and aloe-emodin were identified as the core active ingredients in rhubarb for the treatment of DN. Results from molecular dynamics simulations showed that TP53 and aloe-emodin bound very stably with a binding free energy of − 26.98 kcal/mol between the two. The results of the gene enrichment analysis revealed that the PI3K-Akt signalling pathway, p53 signalling pathway, AGE-RAGE signalling pathway and MAPK signalling pathway might be the key pathways for the treatment of DN, and these pathways were involved in podocyte apoptosis, glomerular mesangial cell proliferation, inflammation and renal fibrosis. Based on the network pharmacology approach and molecular docking technology, we successfully predicted the active compounds and their respective targets. In addition, we illustrated the molecular mechanisms that mediate the therapeutic effects of rhubarb against DN. These findings provided an important scientific basis for further research of the mechanism of rhubarb in the treatment of DN.
Because of the improvement of living conditions, increasingly sedentary lifestyles, increased rates of obesity and an ageing population, the global prevalence of diabetes is increasing at an alarming rate . Diabetic nephropathy (DN) is one of the most serious complications of diabetes and the main cause of end-stage renal failure . The incidence of DN in patients with diabetes is approximately 20–40% . Once DN enters the clinical phase, kidney damage progresses rapidly, which seriously affects patient quality of life . However, current treatments for DN are limited, and they mainly aim to control blood sugar, reduce blood pressure and regulate lipids, thereby improving hypercoagulability to delay the occurrence and development of proteinuria and protect kidney function [5, 6].
In recent years, some hospitals have assessed the efficacy of traditional rhubarb (Radix Rhei Et Rhizome), a widely used traditional Chinese herb, as a long-term treatment for DN [7, 8]. Clinically, rhubarb and its ingredients have been widely used in the treatment of DN with significant efficacy [9, 10]. Oral rhubarb reduced proteinuria, lowered blood sugar levels and improved kidney function in patients with DN . Rhein purified from rhubarb can decrease lipid levels and improve renal lesions in db/db mice with DN, an ideal animal model for type 2 diabetes research . However, the mechanism by which rhubarb alleviates DN has not yet been elucidated. Nevertheless, the integration of bioinformatics and network pharmacology provides a practical approach to explore the mechanisms of action . Network pharmacology can systematically reveal the active components in drugs and predict the relationships between drug components and gene targets [14, 15]. Molecular docking could verify the interaction between the active compounds and central therapeutic targets . Molecular dynamics relies on Newtonian mechanics to assess the binding stability and flexibility of ligand and receptor by simulating their movement . This study is the first to apply network pharmacology, molecular docking and molecular dynamics simulation technology in combination to fully reveal the underlying mechanisms by which rhubarb produces therapeutic effects for DN patients, and to predict the main active compounds and key targets of rhubarb for the treatment of DN, as well as the signalling pathways involved.
Materials and methods
Screening of the active compounds of rhubarb
The Traditional Chinese Medicine Systems Pharmacology (TCMSP) database is an authoritative Chinese herbal medicine database that captures the relationships among drugs, targets and diseases , whereas Traditional Chinese Medicines Integrated Database (TCMID) is a large repository of information pertaining to herbs and herbal ingredients . ‘Radix Rhei Et Rhizome’ was used as the keyword to retrieve all chemical compounds of rhubarb from the TCMSP (http://tcmspw.com/tcmsp.php, Ver.2.3) and TCMID databases (http://www.megabionet.org/tcmid/, Ver.1.0). Oral bioavailability (OB) refers to the percentage of an oral drug that reaches the systemic circulation. It reflects the degree of absorption and utilisation of drugs in the body . Drug similarity (DL) reflects the structural similarity between the compound and drug molecule. Thus, compounds with a high DL value are more likely to exhibit appropriate pharmacodynamic and pharmacokinetic properties . Therefore, OB and DL provide important references for screening active compounds and as the TCMSP database suggested OB ≥ 30% and DL ≥ 0.18 were used as the screening criteria to select the active compounds of rhubarb .
Prediction of the targets of rhubarb
DrugBank (https://www.drugbank.ca/) is a comprehensive bioinformatics and cheminformatics database that provides detailed information on drug types, chemical structures, drug targets and other information . The related targets of the active compounds of rhubarb were gathered from DrugBank under the condition of Homo sapiens. UniProt (http://www.uniprot.org/) is the most informative and extensive protein database, and it is used to standardise gene names . We obtained the target gene names corresponding to the target protein names of active chemical compounds using UniProt.
Collection of gene targets in DN
The human genes associated with DN were gathered from four databases, namely Online Mendelian Inheritance in Man (OMIM), GeneCards, Pharmacogenomics Knowledge Base (PharmGkb) and DrugBank. OMIM (https://omim.org/) is an authoritative and comprehensive database of human genes and genetic phenotypes . GeneCards (https://www.genecards.org/) is an integrative database that provides information on all predicted and annotated human genes . PharmGKB (https://www.pharmgkb.org/) is an integrated resource about how variation in human genetics leads to variation in responses to drugs . DrugBank (https://www.drugbank.ca/), as a comprehensive bioinformatics database, also provides relevant information on disease targets . The search term ‘diabetic nephropathy’ was used to retrieve the DN targets from the four databases. R package-VennDiagram was used to plot the venn diagram to visualize the relationship between the search results for each database .
Therapeutic targets of rhubarb for treating DN
We screened the active compounds of rhubarb and identified their target genes. We also gathered the DN-related genes. Taking the intersection of the aforementioned two groups of genes, the overlapped genes were considered potential therapeutic targets of rhubarb against DN.
Protein–protein interaction (PPI)
The PPIs of the therapeutic targets of rhubarb in the treatment of DN were gathered using STRING (https://string-db.org/, version 11.0), a database of known and predicted PPI that uses bioinformatic strategies to collect information . In this study, we limited the species to ‘Homo sapiens’, collected the PPIs with confidence scores > 0.4 and hid the disconnected nodes in the network.
The active compounds–therapeutic targets network was built by linking the active compounds with the therapeutic targets of rhubarb in the treatment of DN. The PPI network of therapeutic targets was established by connecting the therapeutic targets to their interacting targets. The two networks were visualised using Cytoscape (http://www.cytoscape.org/, version 3.7.2), which can graphically display and analyse networks .
In this study, to screen the central targets, we used the NetworkAnalyzer tool of Cytoscape to calculate three topological parameters for each node in the PPI network of therapeutic targets, including degree, closeness centrality and betweenness centrality [32, 33]. For each variable, higher values indicate greater importance of the node in the network . Nodes with values of all three topological parameters exceeding the median were selected to build sub-network, in which another selection was performed to finally obtain the central targets.
Molecular docking is the most widely used technique in drug design because it can predict the ability of ligands and proteins to bind, as well as the location of binding . In this study, to obtain the core active ingredients in rhubarb for the treatment of DN, molecular docking was performed on the key targets and the active compounds acting on them. The two-dimensional structures of the active compounds of rhubarb were gathered from PubChem. Then their three-dimensional structures were converted using ChemBio3D Ultra 17.0 and optimised by MMFF94 force field to obtain the three-dimensional conformation of the minimum free energy of each compound . The crystal structures of the key targets TP53, PTGS2, MYC, CASP3, CASP8, and JUN were all downloaded from PDB database (Protein Data Bank: www.rcsb.org), and their corresponding PDB IDs were as follows: 3DCY, 5IKR, 5I4Z, 6X8I, 5H31 and 5FV8 [37,38,39,40,41,42,43]. These crystal structures were imported into PyMOL 220.127.116.11 software (https://pymol.org/2/) for various optimisation operations such as dehydration, hydrogenation and isolation of original ligands, and the optimised targets were entered into AutoDockTools − 1.5.6 to construct a docking grid box. The active site of molecular docking was determined using the ligand coordinate in the target protein complex, and Autodock vina 1.1.2 was used for molecular docking research . Twenty conformations were generated for each molecular docking. The best affinity conformation was selected as the final docking conformation. The compound-target interactions in the docking conformation were analyzed by ligplot 2.2 and visualised in Pymol 2.3.
Molecular dynamics simulation
To further investigate the dynamic interaction process and the stability of binding between proteins and small ligands, molecular dynamics simulation was also performed. Molecular dynamics simulation is a popular technique to study protein motion by tracking their conformational changes over time . The GROMACS software package was used to perform molecular dynamics simulation in this study . The protein used the AMBER14SB force field parameter, while the ligand atomic charges were calculated by B3LYP/6-31G* basis set and the ligand topology was built by GAFF2 force field parameter. TIP3P water model was used to add the antagonist sodium ions and hydrogen atoms to neutralize the charge. Electrostatic interactions were treated separately using the Particle Mesh Ewald (PME) and Verlet algorithms. The heavy atoms of proteins were confined, and the steepest descent method was applied to minimize the energy at 50000 steps. The simulation system was equilibrated for 100 ps using canonical ensembles (NVT) and isothermal isobaric ensembles (NPT). Both Van-der-Waals and Coulomb interactions were calculated using a cut-off of 1.4 nm. Finally, the system performed a 100 ns molecular dynamics simulation at constant temperature (300 K) and constant pressure (1 bar) with a time step of 2 fs and the trajectory data was saved every 5 ps.
Gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway enrichment
We used the GO database (http://geneontology.org/) to clarify the possible biological mechanisms of therapeutic targets, including biological process (BP), cell component (CC) and molecular function (MF) terms deemed significant at P < 0.05 . We used the KEGG database (https://www.kegg.jp/) to further explore the vital biological relevance of therapeutic targets and verify the reliability of the integrated results . The enrichment analysis of GO and KEGG was performed using R package-Bioconductor clusterProfiler, which is widely used for gene clusters enrichment analysis .
The flowchart of the study based on network pharmacology is presented in Fig. 1.
Active compounds of rhubarb
In the search, 92 compounds were obtained from the TCMSP and TCMID databases, including flavonoids, steroids, alkaloids, glycosides and triterpenes. According to the OB and DL of the ingredients, 16 compounds were selected as potential active ingredients, and their characteristics are listed in Table 1.
Prediction of the targets of rhubarb
The targets of active compounds were predicted using the DrugBank database, and the target protein names were converted into target gene names using the UniProt database. Finally, 88 targets were predicted for these active compounds of rhubarb, and they were regarded as potential targets (Supplementary Table 1).
Collection of gene targets in DN
Regarding genes related to DN, 3319, 68, 79 and 43 genes were retrieved from GeneCards, OMIM, PharmGkb and DrugBank, respectively, (Supplementary Table 2). The genes commonly predicted by the four databases were considered the genes associated with DN, and the details are depicted in Supplementary Figure 1.
Therapeutic targets of rhubarb for DN
The target genes of rhubarb were predicted, and the genes associated with DN were obtained. Using the aforementioned commonly predicted genes, 37 possible therapeutic targets were gathered, and their features are listed in Table 2.
Active compounds–therapeutic targets network
The active compounds–therapeutic targets network is depicted in Fig. 2, and it revealed the relationships between the active compounds and therapeutic targets, including 46 total nodes (7 compound nodes, 37 therapeutic target nodes, 1 rhubarb node and 1 DN node) and 105 edges. This network fully reveals the characteristics of rhubarb for the treatment of DN through multi-components and multi-targets. The larger the Degree value of the compound node in the network, the more therapeutic targets this compound acts on. As presented in Fig. 2, the two nodes with the largest degree values were aloe-emodin and beta-sitosterol (degree = 17), and their degree values were much larger than those of the other nodes (mean, 6.8), indicating that aloe-emodin and beta-sitosterol may play significant roles in the effects of rhubarb against DN.
Therapeutic target PPI network
To identify the central therapeutic targets of the effects of rhubarb against DN, the PPI network of therapeutic targets is depicted in Fig. 3a, including 37 nodes and 223 edges. NetworkAnalyzer was used to calculate three topological features of the 36 targets, and the details are listed in Table 2. The median degree, betweenness and closeness values of these nodes were 12, 8.45993 and 0.57377, respectively. As shown in Fig. 3b and c, nodes with values of all three topological parameters exceeding the median were selected to build sub-network. In sub-network another selection was performed to finally obtain the central targets (Fig. 3c and d). Finally, 6 genes were identified as central therapeutic targets of the effects of rhubarb against DN, including TP53, CASP8, CASP3, MYC, JUN and PTGS2.
Results of molecular docking
In this study, molecular docking technology was used to obtain the core active ingredients in rhubarb for the treatment of DN. Based on the result of therapeutic target PPI network, the central therapeutic targets were identified, and they were selected for molecular docking with the active compounds of rhubarb acting on them, respectively. According to the receptor-ligand docking theory, docking energy is inversely correlated with binding affinity, and a more negative docking energy indicates stronger binding affinity between the protein and the ligand. The docking details are shown in Table 3. The six pairs of compounds and proteins with the best binding affinity were PTGS2 and rhein, TP53 and aloe-emodin, MYC and aloe-emodin, CASP3 and beta-sitosterol, CASP8 and beta-sitosterol, and JUN and rhein, respectively. The docking data of above six pairs of compounds and proteins were imported into Ligplot software to analyze the compound and protein interactions and visualised in Pymol software, which is depicted in Fig. 4. The binding energy values of the 6 central therapeutic targets with the active compounds acting on them were all less than − 5 kcal/mol, and at least one hydrogen bond formed between them, indicating these compounds have good binding affinity to the key therapeutic targets and can fully play the role of anti-MN. The screening criteria for core active ingredients were as follows: ① Compounds with the highest affinity for each key therapeutic target; ② Compounds that act on the key therapeutic targets must have a docking energy of less than − 5 kcal/mol. The core active ingredients in rhubarb for the treatment of DN were obtained by the above screening criteria were rhein, beta-sitosterol and aloe-emodin.
Results of molecular dynamics simulation
Since proteins and small ligands are in absolute motion, and the conformation involved in the molecular docking is only a relatively stationary state for them. Therefore, to get more confidence in the binding stability, we performed molecular dynamics simulation of TP53 and aloe-emodin, where TP53 was the most central therapeutic target with the highest Degree value in the PPI network, and aloe-emodin was the active ingredient of rhubarb that molecular docking results showed to bind most tightly to TP53. The RMSD curve reflects the fluctuation in protein conformation. As shown in Fig. 5a, the RMSD has some fluctuations in the early stage. However, the RMSD fluctuations of TP53 and aloe-emodin system stabilized after 55 nanoseconds (ns), indicating that the conformation of TP53 did not change significantly after the binding of aloe-emodin to it and the binding of the two was relatively stable. The RMSF curve reflects the fluctuation of the protein amino acid residues. As shown in Fig. 5b, the result of TP53 and aloe-emodin shows that the amino acids 114–134 and 287–354 of TP53 are highly volatile and has greater residue flexibility than other regions. The gyration radius curve reflects the compactness of the overall structure of protein. As shown in Fig. 5c, after a brief fluctuation at the beginning, the gyration radius curve gradually tends to balance, indicating that TP53 conformation is stable and is compactly folded. In addition, we also obtained the trajectory of the RMSD stationary phase (50–100 ns) of TP53 and aloe-emodin and calculated the binding free energy using gmx_MMPBSA. As shown in Fig. 5d, the binding free energy between the two is − 26.98 kcal/mol. Van der Waals energy (VDWAALS, − 30.81 kcal/mol), non-polar solvation energy (ESURF, − 4.07 kcal/mol) and total gas phase free energy (GGAS, − 30.81 kcal/mol) are conducive to the combination of both, while the polar solvation energy (EGB, 7.9 kcal/mol) and total solvation free energy (GSOLV, 3.83 kcal/mol) is not conducive to the interaction of both. As shown in Fig. 5e and f, the binding site of TP53 and aloe-emodin forms a relatively hydrophobic environment with strong hydrophobicity. The residues that interact with aloe-emodin are mainly hydrophobic residues, such as Thr312, Pro309, and Asn310.
GO and KEGG pathway enrichment
To illuminate the complex mechanisms of the effects of rhubarb against DN, we conducted GO analysis of BP, CC and MF for the 37 therapeutic targets. The top 10 most significant BP, CC and MF terms, as determined using p-values, are listed in Fig. 6 and the relationships between the therapeutic targets and the aforementioned entries are depicted in Fig. 7. The specific entries of GO enrichment analysis for BP, CC and MF are listed in Supplementary Table 3.
KEGG pathway enrichment analysis was performed to further clarify the underlying mechanisms of rhubarb in the treatment of DN. As presented in Supplementary Table 4, the 37 therapeutic targets were mapped into 106 KEGG pathways at P < 0.05. The top 30 matched KEGG pathways are depicted in Fig. 8. These 106 KEGG pathways involve human diseases, signalling pathways and pathophysiological mechanisms. The top 10 significant signalling pathways based on p-values were p53 signalling pathway, AGE-RAGE signalling pathway in diabetic complications, oestrogen signalling pathway, thyroid hormone signalling pathway, IL-17 signalling pathway, PI3K-Akt signalling pathway, TNF signalling pathway, VEGF signalling pathway, MAPK signalling pathway, ErbB signalling pathway and Wnt signalling pathway. The therapeutic targets of rhubarb and the associated DN genes in the AGE-RAGE signalling pathway are depicted in Supplementary Figure 2. These findings indicated that rhubarb may play a role in the treatment of DN by regulating the key targets in these signalling pathways, and most therapeutic targets participate in multiple signalling pathways.
DN is the main microvascular complication of diabetes. It significantly increases the morbidity and mortality of cardiovascular diseases, and it is one of the main complications prompting renal replacement therapy worldwide [50, 51]. Rhubarb is a well-known herbal medicine, and medicine preparations based on rhubarb have produced good results in the treatment of DN [52, 53]. However, the specific pharmacological mechanism of action of rhubarb is unclear. In this study, network pharmacology was used to explore the potential active compounds and underlying mechanisms of the effects of rhubarb against DN.
Through screening, we identified 16 active compounds in rhubarb, and after further molecular docking, rhein, beta-sitosterol and aloe-emodin were identified as the core active ingredients in rhubarb for the treatment of DN. Some of these compounds have been reported to have anti-DN effects. Rhein is a monomer of anthraquinone compounds mainly extracted from rhubarb . Many studies revealed that rhein has beneficial anti-oxidant and anti-inflammatory effects on DN [12, 55], and it can alleviate renal fibrosis . Emodin is a crystalline phenolic compound extracted from the root of Rheum palmatum . It can mitigate podocyte apoptosis , inhibit the proliferation of fibroblasts , and suppress inflammation to ameliorate renal dysfunction in DN . Catechin belongs to a group of polyphenolic compounds called flavanols, and it is widely distributed in plant foods [61, 62]. Some studies reported that catechin could ameliorate renal dysfunction in DN through its anti-oxidant and anti-inflammatory properties [63, 64]. These studies highlight the need for further research on the biological functions of the active compounds isolated from rhubarb.
From the therapeutic target PPI network, TP53, CASP8, CASP3, MYC, JUN and PTGS2 were selected as the core targets of rhubarb. These targets may play a major role in the therapeutic effects of rhubarb against DN. Caspases comprise a family of cysteine aspartate-specific proteases that mediate renal apoptosis . Active CASP3 degrades cell-stabilising proteins and other DNA repair enzymes, resulting in apoptotic cell death , and attenuation of CASP3 activity in db/db mice can inhibit the progression of DN . PPAR is a nuclear receptor that plays important roles in intermediary metabolism , and many studies suggested that the Pro12Ala polymorphism of the PPARG gene is a significant independent risk factor for DN [69, 70]. The transcription factor AP-1 (JUN) is widely involved in the transcriptional regulation of multiple genes involved in cell survival, proliferation and apoptosis . JUN activation is crucial for mesangial cell proliferation and extracellular matrix production, and mesangial expansion is a key pathologic feature of DN [72, 73]. These findings indicate that the therapeutic effect of rhubarb against DN is primarily mediated through the regulation of glucose and lipid metabolism, inhibition of glomerular mesangial cell infiltration and suppression of renal fibrosis.
Next, we performed GO and KEGG enrichment analyses of the therapeutic targets. According to the GO analysis, the therapeutic targets displayed strong correlations with BP (e.g., response to steroid hormone, transcription initiation from RNA polymerase II promoter, cellular response to drugs), CC (e.g., cyclin-dependent protein kinase holoenzyme complex, serine/threonine protein kinase complex, membrane raft) and MF terms (e.g., nuclear receptor activity, steroid hormone receptor activity, histone kinase activity). Hence, rhubarb may function through the aforementioned pathways. For example, aldosterone is an endogenous steroid hormone that is involved in the formation and development of DN . Aldosterone breakthrough is a clinical phenomenon that occurs upon treatment with the renin–angiotensin system (RAS) inhibitors angiotensin-converting enzyme inhibitors and angiotensin II receptor antagonists , and this phenomenon attenuates the organ-protective effects of RAS inhibitors against DN through the actions of mineralocorticoid receptors . Several preclinical studies revealed that mineralocorticoid receptor antagonists could ameliorate or cure kidney injury and dysfunction [77, 78]. Therefore, we speculate that suppressing the response to steroid hormone and the activity of steroid hormone receptor is beneficial for the treatment of DN. Tamadher et al. found that histone H3 serine 10 phosphorylation facilitated endothelial activation in DN , suggesting that histone kinase activity is related to gene activation in DN.
According to the KEGG terms, the therapeutic targets of rhubarb against DN were mainly associated with the PI3K-Akt signalling pathway, p53 signalling pathway, AGE-RAGE signalling pathway in diabetic complications and MAPK signalling pathway. Podocytes maintain the glomerular filtration barrier composed of the basement membrane and slit diaphragm , and podocyte apoptosis was associated with increased macroalbuminuria in DN . The PI3K-Akt signalling pathway is an apoptotic signalling transduction pathways, and decreased phosphorylation could give rise to podocyte apoptosis . Huang et al. reported that notoginsenoside R1 could activate the PI3K-Akt signalling pathway to exert anti-apoptotic and renal-protective effects in DN mice . p53 is a transcription factor, and the main outcomes of its activation are cell cycle arrest and apoptosis . Evidence suggests that p53 overexpression is associated with the progression of DN, as p53 mediates podocyte apoptosis related to DN and promotes the expression of pro-fibrotic genes such as plasminogen activator inhibitor-1 [85,86,87]. Inflammation is associated with the development and progression of DN, and many studies support that the MAPK signalling pathway plays a central role in high glucose-induced cell damage and the activation of inflammation [88, 89]. Inhibiting the MAPK signalling pathway can reduce the inflammatory response of DN and protect the kidneys [90, 91]. Studies illustrated that p38 MAPK might play an important role in high glucose-induced epithelial–mesenchymal transition in cultured human renal tubular epithelial cells . Zhong et al. reported that Cordyceps sinensis could mitigate epithelial–mesenchymal transition and the subsequent extracellular matrix deposition to exert a therapeutic effect on experimental diabetic renal fibrosis by inhibiting the p38 MAPK signalling pathway . These findings suggest that the MAPK signalling pathway plays a role in the fibrosis associated with DN. KEGG enrichment analysis demonstrated that the therapeutic effect mechanism of the effects of rhubarb against DN may be closely related to PI3K, P53 and MAPK, suggesting that the therapeutic effect of rhubarb may be mediated through the regulation of several important factors in these signalling pathways, and most therapeutic targets participate in multiple signalling pathways.
Finally, in order to explore the accurate therapeutic mechanism of rhubarb, molecular docking was performed. The results showed that the main active compounds of rhubarb have different affinities with the central therapeutic targets. These interactions are the main forces between proteins and compounds, which form stable complexes of proteins and compounds. Molecular docking studies provide an explanation of the protein-compound interactions, which laid the foundation for further research on the therapeutic mechanism of active compounds. Interestingly, Albersmeyer et al. reported a case of AKI induced by excessive ingestion of rhubarb in the treatment of DN . However, it has not been reported that regular consumption of Rhubarb caused secondary hyperoxaluria and renal failure. In addition, Liu et al. reported that the main active component of rhubarb, emodin, may provide clinical benefits in the treatment of AKI . Therefore, more clinical evidence needs to be collected to further investigate the relationship between Rhubarb and AKI.
Our study still has some limitations. Our study relied on data mining and analysis, and further validation using experimental evidence is required. More databases on the composition and targets information of TCM should be included, which can make the results of our results more reliable. In summary, a comprehensive understanding of rhubarb for the treatment of DN still depends on the common development of multi-disciplines.
We used network pharmacology analysis and molecular docking technology to explore the potential mechanism of the effects of rhubarb against DN, clarifying the core active ingredients in rhubarb were rhein, beta-sitosterol and aloe-emodin, and the key therapeutic targets were TP53, CASP8, CASP3, MYC, JUN and PTGS2. The therapeutic properties of rhubarb against DN arise from the regulation of biological pathways involved in podocyte apoptosis, inflammation and renal interstitial fibrosis. These findings demonstrate the importance of understanding traditional Chinese medicines and establish a base for further research on the pathogenesis of DN, as well as the development of new drugs.
Availability of data and materials
The datasets are included in this published article and its supplementary files. The algorithms used to process the data are available from the corresponding author on reasonable request. Please send the permissions document to email@example.com.
Dong Z, Sun Y, Wei G, Li S, Zhao Z. A nucleoside/nucleobase-rich extract from Cordyceps Sinensis inhibits the epithelial-mesenchymal transition and protects against renal fibrosis in diabetic nephropathy. Molecules (Basel, Switzerland). 2019;24(22):4119.
Kim SS, Kim JH, Kim IJ. Current challenges in diabetic nephropathy: early diagnosis and ways to improve outcomes. Endocrinol Metab (Seoul, Korea). 2016;31(2):245–53.
Bermejo S, Pascual J, Soler MJ. The current role of renal biopsy in diabetic patients. Minerva Med. 2018;109(2):116–25.
Alicic RZ, Rooney MT, Tuttle KR. Diabetic kidney disease: challenges, progress, and possibilities. Clin J Am Soc Nephrol. 2017;12(12):2032–45.
Tziomalos K, Athyros VG. Diabetic nephropathy: new risk factors and improvements in diagnosis. Rev Diabetic Stud. 2015;12(1–2):110–8.
Bhattacharjee N, Barma S, Konwar N, Dewanjee S, Manna P. Mechanistic insight of diabetic nephropathy and its pharmacotherapeutic targets: An update. Eur J Pharmacol. 2016;791:8–24.
Lu Z, Zhong Y, Liu W, Xiang L, Deng Y. The efficacy and mechanism of Chinese herbal medicine on diabetic kidney disease. J Diabetes Res. 2019;2019:2697672.
Wang L, Wang YH, Zhang XH, Yang XL, Wei HL, An ZC, et al. Effectiveness comparisons of traditional Chinese medicine on treating diabetic nephropathy proteinuria: a systematic review and meta-analysis. Medicine. 2019;98(43):e17495.
Guo J, Chen L, Chang B, Zheng M, Wen J, Sun H. Study on treatment of type 2 diabetic nephropathy with rhubarb-based Chinese medicine enema. Clin Focus. 2011;26(18):1595–8.
Li Q, Liu Y, Zhang X, Li T. Treating 50 cases of diabetic nephropathy with Bushen Quyu decoction Shaanxi. J Tradit Chin Med. 2011;32(08):967–9.
Qi Y. Observation on the intervention effect of rhubarb on patients with diabetic nephropathy. Health Protect Promot. 2015;22:65–7.
Gao Q, Qin WS, Jia ZH, Zheng JM, Zeng CH, Li LS, et al. Rhein improves renal lesion and ameliorates dyslipidemia in db/db mice with diabetic nephropathy. Planta Med. 2010;76(1):27–33.
Jain B, Raj U, Varadwaj PK. Drug target interplay: a network-based analysis of human diseases and the drug targets. Curr Top Med Chem. 2018;18(13):1053–61.
Li S, Zhang B. Traditional Chinese medicine network pharmacology: theory, methodology and application. Chin J Nat Med. 2013;11(2):110–20.
Zhou Z, Chen B, Chen S, Lin M, Chen Y, Jin S, et al. Applications of network pharmacology in traditional Chinese medicine research. Evid-based Complement Alternat Med. 2020;2020:1646905.
Luo Q, Shi X, Ding J, Ma Z, Chen X, Leng Y, et al. Network pharmacology integrated molecular docking reveals the Antiosteosarcoma mechanism of Biochanin A. Evid Based Complement Alternat Med. 2019;2019:1410495.
Nam KH. Molecular dynamics-from small molecules to macromolecules. Int J Mol Sci. 2021;22(7):3761.
Ru J, Li P, Wang J, Zhou W, Li B, Huang C, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. J Chem. 2014;6:13.
Huang L, Xie D, Yu Y, Liu H, Shi Y, Shi T, et al. TCMID 2.0: a comprehensive resource for TCM. Nucleic Acids Res. 2018;46(D1):D1117–20.
Xie Z, Tian Y, Lv X, Xiao X, Zhan M, Cheng K, et al. The selectivity and bioavailability improvement of novel oral anticoagulants: An overview. Eur J Med Chem. 2018;146:299–317.
Lin L, Wan L, He H, Liu W. Drug vector representation: a tool for drug similarity analysis. Mol Genet Genomics. 2020;295(4):1055–62.
He D, Huang JH, Zhang ZY, Du Q, Peng WJ, Yu R, et al. A network pharmacology-based strategy for predicting active ingredients and potential targets of LiuWei DiHuang pill in treating type 2 diabetes mellitus. Drug Des Devel Ther. 2019;13:3989–4005.
Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018;46(D1):D1074–d1082.
Breuza L, Poux S, Estreicher A, Famiglietti ML, Magrane M, Tognolli M, et al. The UniProtKB guide to the human proteome. Database. 2016;2016. https://doi.org/10.1093/database/bav120.
Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. OMIM.org: Online Mendelian inheritance in man (OMIM®), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2015;43(Database issue):D789–98.
Stelzer G, Rosen N, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, et al. The GeneCards suite: from gene data mining to disease genome sequence analyses. Curr Protoc Bioinformatics. 2016;54:1.30.31–31.30.33.
Barbarino JM, Whirl-Carrillo M, Altman RB, Klein TE. PharmGKB: a worldwide resource for pharmacogenomic information. Wiley Interdiscip Rev Syst Biol Med. 2018;10(4):e1417.
Probst D, Reymond JL. Exploring DrugBank in virtual reality chemical space. J Chem Inf Model. 2018;58(9):1731–5.
Chen H, Boutros PC. VennDiagram: a package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinform. 2011;12(1):1–7.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, 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–d613.
Otasek D, Morris JH, Bouças J, Pico AR, Demchak B. Cytoscape automation: empowering workflow-based network analysis. Genome Biol. 2019;20(1):185.
Batt G, Besson B, Ciron PE, de Jong H, Dumas E, Geiselmann J, et al. Genetic network analyzer: a tool for the qualitative modeling and simulation of bacterial regulatory networks. Methods Mol Biol (Clifton, NJ). 2012;804:439–62.
Liu X, Wu J, Zhang D, Wang K, Duan X, Zhang X. A network pharmacology approach to uncover the multiple mechanisms of Hedyotis diffusa Willd. On colorectal cancer. Evid Based Complement Alternat Med. 2018;2018:6517034.
Azuaje FJ, Zhang L, Devaux Y, Wagner DR. Drug-target network in myocardial infarction reveals multiple side effects of unrelated drugs. Sci Rep. 2011;1:52.
Plewczynski D, Łaźniewski M, Augustyniak R, Ginalski K. Can we trust docking results? Evaluation of seven commonly used programs on PDBbind database. J Comput Chem. 2011;32(4):742–55.
Tao Q, Du J, Li X, Zeng J, Tan B, Xu J, et al. Network pharmacology and molecular docking analysis on molecular targets and mechanisms of Huashi Baidu formula in the treatment of COVID-19. Drug Dev Ind Pharm. 2020;46(8):1345–53.
Burley SK, Berman HM, Kleywegt GJ, Markley JL, Nakamura H, Velankar S. Protein data Bank (PDB): the single global macromolecular structure archive. Methods Mol Biol (Clifton, NJ). 2017;1607:627–41.
McCoy JG, Bingman CA, Wesenberg GE, Phillips JGN. Crystal structure a TP53-induced glycolysis and apoptosis regulator protein from Homo sapiens. PDB. 2008. https://doi.org/10.2210/pdb3DCY/pdb.
Orlando BJ, Malkowski MG. The structure of Mefenamic acid bound to human Cyclooxygenase-2. PDB. 2016. https://doi.org/10.2210/pdb5IKR/pdb.
Koelmel W, Jung LA, Kuper J, Eilers M, Kisker C. Structure of apo OmoMYC. PDB. 2017. https://doi.org/10.2210/pdb5i4z/pdb.
Solania A, Xu JH, Wolan DW. Caspase-3 in complex with ketomethylene inhibitor reveals tetrahedral adduct. PDB. 2021. https://doi.org/10.2210/pdb6X8I/pdb.
Shen C, Pei J, Guo X, Quan J. Structural basis for dimerization of the death effector domains of Caspase-8. PDB. 2017. https://doi.org/10.2210/pdb5H31/pdb.
Zaccai NR, Mason JM, Brady RL. Structure of cJun-FosW coiled coil complex. PDB. 2017. https://doi.org/10.2210/pdb5FV8/pdb.
Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31(2):455–61.
Collier TA, Piggot TJ, Allison JR. Molecular dynamics simulation of proteins. In: Protein Nanotechnology. Berlin: Springer; 2020. p. 311–27.
Rakhshani H, Dehghanian E, Rahati A. Enhanced GROMACS: toward a better numerical simulation framework. J Mol Model. 2019;25(12):1–8.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The gene ontology consortium. Nat Genet. 2000;25(1):25–9.
Kanehisa M, Furumichi M, Tanabe M, Sato Y, Morishima K. KEGG: new perspectives on genomes, pathways, diseases and drugs. Nucleic Acids Res. 2017;45(D1):D353–d361.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012;16(5):284–7.
Wada T. Clinico-pathological progress in diabetic nephropathy. Rinsho byori Jpn J Clin Pathol. 2016;64(8):954–9.
Qi C, Mao X, Zhang Z, Wu H. Classification and differential diagnosis of diabetic nephropathy. J Diabetes Res. 2017;2017:8637138.
Xu L, Ma L. Research progress on the mechanism of rhubarb inhibiting renal fibrosis in diabetic nephropathy. Xinjiang J Tradit Chin Med. 2018;36(04):137–40.
Zhong J. Brief introduction of experimental research on treatment of diabetic nephropathy with single Chinese medicine and extract. J Pract Tradit Chin Internal Med. 2017;31(11):71–4.
Zhang Q, Liu L, Lin W, Yin S, Duan A, Liu Z, et al. Rhein reverses klotho repression via promoter demethylation and protects against kidney and bone injuries in mice with chronic kidney disease. Kidney Int. 2017;91(1):144–56.
Lin YJ, Zhen YZ, Wei JB, Wei J, Dai J, Gao JL, et al. Rhein lysinate protects renal function in diabetic nephropathy of KK/HlJ mice. Exp Ther Med. 2017;14(6):5801–8.
Chen Y, Mu L, Xing L, Li S, Fu S. Rhein alleviates renal interstitial fibrosis by inhibiting tubular cell apoptosis in rats. Biol Res. 2019;52(1):50.
Chen T, Zheng LY, Xiao W, Gui D, Wang X, Wang N. Emodin ameliorates high glucose induced-podocyte epithelial-mesenchymal transition in-vitro and in-vivo. Cell Physiol Biochem. 2015;35(4):1425–36.
Tian N, Gao Y, Wang X, Wu X, Zou D, Zhu Z, et al. Emodin mitigates podocytes apoptosis induced by endoplasmic reticulum stress through the inhibition of the PERK pathway in diabetic nephropathy. Drug Des Devel Ther. 2018;12:2195–211.
Ning Y, Wang J, Qu S. Effect of emodin on human kidney fibroblast proliferation. Zhongguo Zhong xi yi jie he za zhi Zhongguo Zhongxiyi jiehe zazhi = Chin J Integr Tradit Western Med. 2000;20(2):105–6.
Jing D, Bai H, Yin S. Renoprotective effects of emodin against diabetic nephropathy in rat models are mediated via PI3K/Akt/GSK-3β and Bax/caspase-3 signaling pathways. Exp Ther Med. 2017;14(5):5163–9.
Yokozawa T, Nakagawa T, Kitani K. Antioxidative activity of green tea polyphenol in cholesterol-fed rats. J Agric Food Chem. 2002;50(12):3549–52.
Bell JR, Donovan JL, Wong R, Waterhouse AL, German JB, Walzem RL, et al. (+)-Catechin in human plasma after ingestion of a single serving of reconstituted red wine. Am J Clin Nutr. 2000;71(1):103–8.
Chennasamudram SP, Kudugunti S, Boreddy PR, Moridani MY, Vasylyeva TL. Renoprotective effects of (+)-catechin in streptozotocin-induced diabetic rat model. Nutr Res (New York, NY). 2012;32(5):347–56.
Zhu D, Wang L, Zhou Q, Yan S, Li Z, Sheng J, et al. (+)-Catechin ameliorates diabetic nephropathy by trapping methylglyoxal in type 2 diabetic mice. Mol Nutr Food Res. 2014;58(12):2249–60.
Susztak K, Raff AC, Schiffer M, Böttinger EP. Glucose-induced reactive oxygen species cause apoptosis of podocytes and podocyte depletion at the onset of diabetic nephropathy. Diabetes. 2006;55(1):225–33.
Fuchs Y, Steller H. Live to die another way: modes of programmed cell death and the signals emanating from dying cells. Nat Rev Mol Cell Biol. 2015;16(6):329–44.
Ghosh S, Khazaei M, Moien-Afshari F, Ang LS, Granville DJ, Verchere CB, et al. Moderate exercise attenuates caspase-3 activity, oxidative stress, and inhibits progression of diabetic renal disease in db/db mice. Am J Physiol Renal Physiol. 2009;296(4):F700–8.
Elix C, Pal SK, Jones JO. The role of peroxisome proliferator-activated receptor gamma in prostate cancer. Asian J Androl. 2018;20(3):238–43.
Liu L, Zheng T, Wang F, Wang N, Song Y, Li M, et al. Pro12Ala polymorphism in the PPARG gene contributes to the development of diabetic nephropathy in Chinese type 2 diabetic patients. Diabetes Care. 2010;33(1):144–9.
Yu R, Bo H, Huang S. Association between the PPARG gene polymorphism and the risk of diabetic nephropathy: a meta-analysis. Genet Test Mol Biomarkers. 2012;16(5):429–34.
Karin M. The regulation of AP-1 activity by mitogen-activated protein kinases. Philos Trans R Soc Lond B Biol Sci. 1996;351(1336):127–34.
Micakovic T, Papagiannarou S, Clark E, Kuzay Y, Abramovic K, Peters J, et al. The angiotensin II type 2 receptors protect renal tubule mitochondria in early stages of diabetes mellitus. Kidney Int. 2018;94(5):937–50.
Ahn JD, Morishita R, Kaneda Y, Kim HJ, Kim YD, Lee HJ, et al. Transcription factor decoy for AP-1 reduces mesangial cell proliferation and extracellular matrix production in vitro and in vivo. Gene Ther. 2004;11(11):916–23.
Sato A. The necessity and effectiveness of mineralocorticoid receptor antagonist in the treatment of diabetic nephropathy. Hypertens Res. 2015;38(6):367–74.
Sato A, Saruta T. Aldosterone breakthrough during angiotensin-converting enzyme inhibitor therapy. Am J Hypertens. 2003;16(9 Pt 1):781–8.
Sato A, Fukuda S. Effect of aldosterone breakthrough on albuminuria during treatment with a direct renin inhibitor and combined effect with a mineralocorticoid receptor antagonist. Hypertens Res. 2013;36(10):879–84.
Barrera-Chimal J, Girerd S, Jaisser F. Mineralocorticoid receptor antagonists and kidney diseases: pathophysiological basis. Kidney Int. 2019;96(2):302–19.
Zuo C, Xu G. Efficacy and safety of mineralocorticoid receptor antagonists with ACEI/ARB treatment for diabetic nephropathy: a meta-analysis. Int J Clin Pract. 2019;73(12):e13413.
Alghamdi TA, Batchu SN, Hadden MJ, Yerra VG, Liu Y, Bowskill BB, et al. Histone H3 serine 10 phosphorylation facilitates endothelial activation in diabetic kidney disease. Diabetes. 2018;67(12):2668–81.
Reiser J, Kriz W, Kretzler M, Mundel P. The glomerular slit diaphragm is a modified adherens junction. J Am Soc Nephrol. 2000;11(1):1–8.
Beverly LJ, Starczynowski DT. IRAK1: oncotarget in MDS and AML. Oncotarget. 2014;5(7):1699–700.
Takano Y, Yamauchi K, Hayakawa K, Hiramatsu N, Kasai A, Okamura M, et al. Transcriptional suppression of nephrin in podocytes by macrophages: roles of inflammatory cytokines and involvement of the PI3K/Akt pathway. FEBS Lett. 2007;581(3):421–6.
Huang G, Lv J, Li T, Huai G, Li X, Xiang S, et al. Notoginsenoside R1 ameliorates podocyte injury in rats with diabetic nephropathy by activating the PI3K/Akt signaling pathway. Int J Mol Med. 2016;38(4):1179–89.
Chen Y, Yan R, Li B, Liu J, Liu X, Song W, et al. Silencing CCNG1 protects MPC-5 cells from high glucose-induced proliferation-inhibition and apoptosis-promotion via MDM2/p53 signaling pathway. Int Urol Nephrol. 2020;52(3):581–93.
Deshpande SD, Putta S, Wang M, Lai JY, Bitzer M, Nelson RG, et al. Transforming growth factor-β-induced cross talk between p53 and a microRNA in the pathogenesis of diabetic nephropathy. Diabetes. 2013;62(9):3151–62.
Samarakoon R, Overstreet JM, Higgins SP, Higgins PJ. TGF-β1 → SMAD/p53/USF2 → PAI-1 transcriptional axis in ureteral obstruction-induced renal fibrosis. Cell Tissue Res. 2012;347(1):117–28.
Jung DS, Lee SH, Kwak SJ, Li JJ, Kim DH, Nam BY, et al. Apoptosis occurs differentially according to glomerular size in diabetic kidney disease. Nephrol Dial Transplant. 2012;27(1):259–66.
Duyndam MC, Hulscher ST, van der Wall E, Pinedo HM, Boven E. Evidence for a role of p38 kinase in hypoxia-inducible factor 1-independent induction of vascular endothelial growth factor expression by sodium arsenite. J Biol Chem. 2003;278(9):6885–95.
Lu HJ, Tzeng TF, Liou SS, Da Lin S, Wu MC, Liu IM. Polysaccharides from Liriopes Radix ameliorate streptozotocin-induced type I diabetic nephropathy via regulating NF-κB and p38 MAPK signaling pathways. BMC Complement Altern Med. 2014;14:156.
Song Y, Wang X, Qin S, Zhou S, Li J, Gao Y. Esculin ameliorates cognitive impairment in experimental diabetic nephropathy and induces anti-oxidative stress and anti-inflammatory effects via the MAPK pathway. Mol Med Rep. 2018;17(5):7395–402.
Peng L, Li J, Xu Y, Wang Y, Du H, Shao J, et al. The protective effect of Beraprost sodium on diabetic nephropathy by inhibiting inflammation and p38 MAPK signaling pathway in high-fat diet/Streptozotocin-induced diabetic rats. Int J Endocrinol. 2016;2016:1690474.
Lv ZM, Wang Q, Wan Q, Lin JG, Hu MS, Liu YX, et al. The role of the p38 MAPK signaling pathway in high glucose-induced epithelial-mesenchymal transition of cultured human renal tubular epithelial cells. PLoS One. 2011;6(7):e22806.
Albersmeyer M, Hilge R, Schröttle A, Weiss M, Sitter T, Vielhauer V. Acute kidney injury after ingestion of rhubarb: secondary oxalate nephropathy in a patient with type 1 diabetes. BMC Nephrol. 2012;13(1):1–5.
Liu Y, Li M, Teh L, Zhao L, Ye N, Wu L, et al. Emodin-mediated treatment of acute kidney injury. Evid Based Complement Alternat Med. 2022;2022:5699615.
The authors express their gratitude to Chen Yiying from nephrology department of Jilin University, who helped us technologic support.
No funding was obtained for this study.
Ethics approval and consent to participate
This study is a bioinformatics analysis based on data mining, without the requirement of ethical review.
Consent for publication
This study is a bioinformatics analysis based on data mining without involving subjects.
The authors declare that they have no competing financial and non-financial interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Fu, S., Zhou, Y., Hu, C. et al. Network pharmacology and molecular docking technology-based predictive study of the active ingredients and potential targets of rhubarb for the treatment of diabetic nephropathy. BMC Complement Med Ther 22, 210 (2022). https://doi.org/10.1186/s12906-022-03662-6
- Diabetic nephropathy
- Network pharmacology
- Chinese traditional medicine
- Molecular docking