Bingyu Du1,2,3,4†
Yanyan Yin1,2,4†
Yuqing Wang5
Hui Fu6
Helin Sun5
Zhaodi Yue3,4
Shaohong Yu1,2,3*
Zhongwen Zhang1,2,5*- 1Department of Endocrinology and Metabology, Shandong University of Traditional Chinese Medicine, The First Affiliated Hospital of Shandong First Medical University and Shandong Provincial Qianfoshan Hospital, Jinan, China
- 2Department of Rehabilitation Medicine, The Second Clinical Medical College and Medical College, Shandong University of Traditional Chinese Medicine, Jinan, China
- 3Department of Rehabilitation Medicine, The Second Affiliated Hospital of Shandong University of Traditional Chinese Medicine, Jinan, China
- 4College of Rehabilitation Medicine, Shandong University of Traditional Chinese Medicine, Jinan, China
- 5Shandong Key Laboratory of Rheumatic Disease and Translational Medicine, Department of Endocrinology and Metabology, Shandong Institute of Nephrology, The First Affiliated Hospital of Shandong First Medical University and Shandong Provincial Qianfoshan Hospital, Jinan, China
- 6The Clinical Medical College, Cheeloo Medical College of Shandong University, Jinan, China
Aims: To evaluate the effectiveness and potential mechanism of calcium dobesilate (CaD) in diabetic kidney disease (DKD) patients.
Methods: We searched for available randomized controlled studies on DKD patients’ treatment with CaD through open databases. Continuous variables were expressed as standardized mean difference (SMD) with a 95% confidence interval (CI). The putative targets and possible pathways of CaD on DKD were analyzed by network pharmacology. Molecular docking was employed to verify the match between CaD and the target genes.
Results: In the meta-analysis, 42 trials were included, involving 3,671 DKD patients, of which 1,839 received CaD treatment in addition to conventional treatment, while 1,832 received conventional treatment. Compared with routine therapy, the levels of serum creatinine (Scr) and blood urea nitrogen (BUN) significantly decreased in the CaD treatment (early stage of DKD, Scr: p < 0.00001; BUN: p < 0.0001; clinical stage of DKD, Scr: p < 0.00001; BUN: p < 0.00001; kidney failure stage, Scr: p = 0.001; BUN: p = 0.004). The levels of serum cystatin C (Cys-C), urine levels of molecules reflecting kidney function (urinary albumin excretion rate (UAER) and micro glycoprotein), and inflammatory factors [hypersensitive c-reactive protein (hs-CRP)] were reduced compared with control groups, while glomerular filtration rate (GFR) was increased in patients treated with CaD for 12 weeks. CaD also showed a better effect on improving endothelial function. Network pharmacology results showed that the interaction pathway between CaD and DKD was mainly enriched in MAPK and chemokine signaling pathways. AKT1, CASP3, IGF1, MAPK8, and CCL5 might be the key targets for CaD in treating DKD.
Conclusion: Combination with CaD is effective and safe in patients with DKD. Inhibition of MAPK and chemokine signaling pathways might be vital in treating CaD in DKD patients.
1 Introduction
Diabetes mellitus (DM) has afflicted around 422 million people worldwide and has become a leading cause of morbidity and mortality (Zhang et al., 2017). Diabetic kidney disease (DKD) is one of the most severe microvascular consequences of diabetes, accounting for 30%–50% of all end-stage kidney disease patients (Qin et al., 2017). It also imposes a significant financial cost on patients, families, and society (Zhang, 2012). The severity of DKD was divided into five stages using an albuminuria-based methodology.
The interstitial space grows in the first two phases, the mesangial volume rises, and the glomerular basement membrane thickens. They are silent stages because there is no detectable microalbuminuria in clinical practice, and there are currently no effective indicators for their detection (Salem et al., 2020). The early and clinical stages are the third and fourth stages, respectively. Microalbuminuria to overt proteinuria is a sign of progressing from early stage to clinical overt diabetic kidney disease. End-stage renal failure (ESRD) is the ultimate stage (Cai et al., 2017). Because DKD is still reversible in its early phases, novel therapeutic medicines for aggressive therapy are urgently needed to avert progression to ESRD.
Calcium dobesilate (CaD) is a microcirculation-improving medication that improves hemodynamics, inhibits inflammatory responses, and suppresses interstitial fibrosis, among other things (Chen and Bai, 2017). CaD has long been utilized to treat diabetic retinopathy (DR) due to its potential to reduce oxidative stress by decreasing the activation of the p38MAPK and NF-B pathways (Liu et al., 2019; Ashraf et al., 2021). CaD’s efficacy in kidney illness (Zhang et al., 2021), chronic venous insufficiency (Ciapponi et al., 2004), thrombotic disorders (Michal and Giessinger, 1985), and various types of cardiac disease (Besirli et al., 2012) has received great attention recently. CaD considerably affects DKD in lowering the urine albumin excretion rate (Qin et al., 2017) and has been studied extensively in clinical trials (Zhou et al., 2018). There is little authoritative conclusion on the benefits and possible adverse effects due to the limited sample sizes and variable findings of currently available randomized controlled studies (RCTs) on CaD.
Some medications have a wide range of effects on humans, but they all point in the same therapeutic direction, implying that they can work on diverse targets in the same pathways. We have found that network-related approaches can be employed to emphasize our findings of drug action mechanisms in various data layers in the drug development process (Boezio et al., 2017). Based on the network pharmacology mechanism, we created a network between the putative CaD targets and the implicated gene targets of DKD. The PPI network is formed when common targets interact with each other. Furthermore, we may use GO and KEGG analyses on these proteins to identify the main pathway that plays a crucial role in CaD’s treatment of DKD.
2 Materials and methods
2.1 Meta-analysis
2.1.1 Randomized controlled study preparation
The PRISMA (preferred reporting items for systematic reviews and meta-analyses) guidelines were followed in this meta-analysis (Page et al., 2021). The PubMed Database, MED-LINE, Global Health, EMBASE, EBSCOhost, Cochrane Library, China National Knowledge Internet, VIP, Wanfang, and SinoMed databases were used to find the literature. Information from well-known registries, such as Current Controlled Trials, the World Health Organization International Clinical Trials Registry Platform, the Clinicaltrials.gov trials registry, and published review papers and editorials, was considered. The search terms were “Calcium Dobesilate” and “diabetic kidney disease” or “diabetic kidney disease” or “diabetic kidney disease” or “randomized” or “double-blind trial”, with no restrictions on subheadings or language. Other likely relevant citations were found in the reference lists of all included papers, and the literature not found in the abovementioned electronic databases was manually reviewed.
2.1.2 Inclusion and exclusion criteria of studies
The inclusive criteria were as follows: 1) study subjects were diagnosed as DKD according to the corresponding guidelines, 2) all patients were randomized to receive treatment with CaD and contemporary medications or contemporary medications alone, 3) sample size in each study group was ≥ 15cases, 4) follow-up in each study group was ≥ 8 weeks; and 5) the outcomes were quantitative to facilitate outcome analysis.
Exclusion criteria were as follows: 1) studies were nonrandomized or nonblinded, 2) patients enrolled had no definite diagnosis, 3) different medications were compared, 4) studies reported only symptomatic changes in patients without objective laboratory measurements; and 5) methodological quality was poor with a Jadad score <2.
2.1.3 Stage of diabetic kidney disease
DKD is categorized into five stages according to the albuminuria-based classification. There was a hemodynamic alteration at the start of the first stage, with increased glomerular capillary hydrostatic pressure but no abnormalities in the ultrastructure. Hyperglycemic effects begin in the second part of the first stage, with thickening of the glomerular basement membrane, increased mesangial volume, and interstitial expansion. Because microalbuminuria cannot be measured in clinical practice and no suitable test marker has yet been discovered, the second stage is quiet (Salem et al., 2020). The early stage of diabetic kidney disease is the third stage. Previous structural changes had become severe, and microalbuminuria had been diagnosed. The fourth stage is known as clinically severe diabetic kidney disease, and it is at this stage that these changes may progress to significant proteinuria, formerly known as “macroalbuminuria.” Microalbuminuria to overt proteinuria is a marker of progression from early stage to clinical overt diabetic kidney impairment. ESRD is the final stage (Cai et al., 2017).
2.1.4 Statistics analysis
RevMan version 5.3 was used to synthesize the data. Continuous variables were expressed as standardized mean difference (SMD) with a 95% confidence interval (CI). Chi-square and I2 tests were used to test heterogeneity. Nonheterogeneous results (P˃0.1, I2 < 50%) were adopted for the fixed effects model, and heterogenous results (p ≤ 0.1, I2 ˃ 50%) were adopted for the random effects model. Statistical significance was defined as a two-tailed p < 0.05. The fail-safe number was used to estimate the extent of publication bias (Cheng et al., 2016). The formula for the fail-safe number was Nfs0.05 = (ΣZ/1.64)2-S, where S is the number of the included trials.
2.2 Network pharmacology
2.2.1 Predicting potential targets of calcium dobesilate
The bioactive component of CaD was found in the published literature using the keyword “Calcium Dobesilate” in the PubMed database. The 2D structure of CaD was then retrieved from PubChem (Kim et al., 2021) and uploaded to Pharmmapper (Wang et al., 2017) to forecast the drug’s potential targets. All target names were entered into Uniprot sites (UniProt Consortium., 2019) and selected by Homo Saipan species to standardize the gene information.
2.2.2 Screening of targets for diabetic kidney disease
The target genes of DKD were found using the key phrases “diabetic kidney disease” or “diabetic kidney dis-ease” in the Online Mendelian Inheritance in Man database (OMIM) (Amberger et al., 2015), Gene Cards database (Stelzer et al., 2016), and DisGeNET database (Piñero et al., 2017).
2.2.3 Construction of protein–protein interaction network
To identify the interaction targets of CaD in the treatment of DKD, we selected the online drawing tool Interactive Venn (Heberle et al., 2015) to draw a Venn diagram, whose overlapping section represented the typical targets for CaD and DKD. These common targets were uploaded to the STRING 11.0 platform (Szklarczyk et al., 2019), and the PPI network was built using the STRING database and the Network Analyzer plugin of Cytoscape (Shannon et al., 2003). The stronger the interaction in a network, the higher the degree. The more significant degree nodes, which may play a critical role in the overall network, were screened based on the network’s topological qualities.
2.2.4 Enrichment analysis
Metascape1 (Zhou et al., 2019) was used to perform Gene Ontology (GO) functional analysis,Kyoto Encyclopedia of Genes, and Genomes (KEGG) pathway enrichment analysis, with p ≤ 0.01 as the cut-off threshold. Bioinformatics web tools2 and the EHBIO Gene Technology Platform3 were then used to display the top 10 GO items and 20 KEGG pathways that met the criteria.
2.2.5 Construction of the component–target–pathway network
The integrated network of component–target–pathway was constructed using Cytoscape 3.7.1. The topology parameters of the network were analyzed with the built-in tool Network Analyzer in Cytoscape to identify the relationships of protein targets with components and the involved pathways.
2.2.6 Molecular docking verification
Suitable 3D structure “pdb” files of the core targets were downloaded from RCSB Protein Data Bank (Berman et al., 2000). PyMol2.6.0 embellished the downloaded complexes to remove the original ligand and water molecules. The “sdf” file of the DKD bioactive component from the PubChem database was obtained, and its format was converted to a “pdb” file by Open Babel2.4.0 (O'Boyle et al., 2011). Then, we used AutoDockTools1.5.6 (Goodsell and Olson, 1990) to evaluate and verify the binding affinity of the compound–target relationship. The critical models were visualized by PyMol2.6.0 software and Discovery Studio4.5 software. A flow chart of the meta-analysis, network pharmacology, and molecular docking is shown in Figure 1.
FIGURE 1. Flow chart of the investigation of calcium dobesilate in the treatment of diabetic kidney disease.
3 Result
3.1 Meta-analysis
A total of 939 trials were retrieved through database searches—PubMed (48), MEDLINE (17), EMBASE (25), Cochrane library (5), China National Knowledge Internet (248), VIP (219), Wanfang (277), and SinoMed (100). Five hundred and thirty-seven duplicate records were eliminated using endnote software. Subsequently, among the remaining 402 pieces, 123 trials were excluded from the primary screening according to the inclusion and exclusion criteria. Of the remaining 279 trials, 273 were evaluated as eligible, and 42 shots were finally included for meta-analysis (Figure 2; Table 1).
FIGURE 2. Flow chart of the systematic search process.
TABLE 1. Meta-analysis of the randomized controlled study of DKD patients’ treatment with CaD.
3.1.1 Blood kidney function (Scr&BUN)
Among all of the 42 trials, 38 trials evaluated the Scr and BUN. Thirty-five of the 38 trials were focused on DKD in its early stages. Two trials (Zhang et al., 2013; Wang et al., 2015) focused on DKD patients in the clinical phase, while one (Wu et al., 2007) concentrated on DKD patients in the kidney failure stage. As for the studies performed on the early stage of DKD patients, the results showed that the function of the kidney was markedly better in the CaD treatment group than that of the control groups (Scr: SMD = −0.84; 95% CI, −1.05 to −0.62; p < 0.00001; I2 = 87%, BUN: SMD = −0.64; 95% CI, −0.96 to −0.32; p < 0.0001; I2 = 86%). The early phase was subdivided into 8-week and 12-week treatment cycles in most of the included literature. Furthermore, the results showed that CaD treatment for 8 weeks (Scr: SMD = −0.45; 95% CI, −0.73 to −0.17; p = 0.002; I2 = 78%, BUN: SMD = −0.48; 95% CI, −0.65 to −0.31; p < 0.00001; I2 = 0%) and 12 weeks (Scr: SMD = −0.97; 95% CI, −1.27 to −0.67; p < 0.00001; I2 = 89%, BUN: SMD = −0.55; 95% CI, −1.03 to −0.07; p = 0.03; I2 = 87%) could significantly reduce Scr and BUN levels compared to the control group. In the clinical stage of DKD, patients treated with CaD for 12 weeks also exhibited a statistically significant reduction in the expression level of the Scr and BUN compared with control groups (Scr: SMD = −2.08; 95% CI, −2.45 to −1.72; p < 0.00001; I2 = 99%; BUN: SMD = −0.70; 95% CI, −0.97 to −0.43; p < 0.00001; I2 = 19%). In the kidney failure stage of DKD, the data showed the same trend, that the groups treated with CaD for 12 weeks had a better effect on reducing Scr and BUN than comparators (Scr: SMD = −0.71; 95% CI, −1.15 to −0.27; p = 0.001; BUN: SMD = −0.63; 95% CI, −1.07 to −0.20; p = 0.004).
3.1.2 Serum cystatin C
Serum cystatin C (Cys-C) has been recognized as an ideal marker of kidney function (Dharnidharka et al., 2002). Nine trials (Hong, 2013; Wang et al., 2015; Cen et al., 2016; Liu et al., 2016; Chen and Bai, 2017; Fan and Ma, 2017; Qin et al., 2017; Li, 2018; and Zeng et al., 2018) reported serum Cys-C. The treatment durations were both 8 and 12 weeks. A total of 1,036 patients were included, 518 of whom were in the treatment group and 518 in the control group. The results proved that the CaD groups exhibited a statistically significant reduction in serum Cys-C compared with the control groups (SMD = −0.95; 95% CI, −1.29 to −0.61; p < 0.00001; I2 = 85%). Subgroup analysis was conducted to determine that the source of heterogeneity was from different treatment methods: one trial (Qin et al., 2017) reported the patients treated with CaD plus alprostadil and routine treatment, another (Wang et al., 2015) focused on the patients treated with CaD plus α-thioctic acid plus routine treatment, and the remaining studies (Zhang, 2012; Hong, 2013; Cen et al., 2016; Liu et al., 2016; Chen and Bai, 2017; Fan and Ma, 2017; Li, 2018) compared two treatment methods, including CaD, angiotensin-converting enzyme inhibitor (ACEI)/angiotensin II receptor blockers (ARBs) and routine treatment (CaD group) or ACEI/ARB and routine treatment (control group) (MD = −0.37; 95% CI, −0.43 to −0.30; p < 0.00001; I2 = 18%).
3.1.3 Glomerular filtration rate
Glomerular filtration rate (GFR) is widely recognized as a comprehensive measure of renal function, and its assessment is essential for clinical practice, research, and public health (Levey et al., 2020). GFR was reported in two trials (Wu et al., 2007; Zhang, 2018), and both treatment cycles were divided into 8 and 12 weeks. The results revealed that after 12 weeks of treatment, patients in the CaD group had considerably higher GFR levels than the control group (SMD = 1.66; 95% CI, 1.26 to 2.07; p < 0.00001; I2 = 22%), while patients who received CaD for 8 weeks had increased GFR levels, but there was no statistical difference (SMD = 0.17; 95% CI, −0.14 to 0.48; p = 0.27; I2 = 0%).
3.1.4 Urinary albumin excretion rate, 24 h urinary protein/24 h urinary albumin, α1-MG, and β2-MG
The kidney function can be quantified by detecting the level of some molecules in urine. These molecules are divided into primary outcomes and secondary outcomes. The primary outcomes include urine albumin excretion rate (UAER) and 24 h urine protein/24 h urine albumin, and the secondary outcomes consist of alpha-1-microglycoprotein (α1-MG) and β2-microglobulin (β2-MG).
3.1.4.1 Urinary albumin excretion rate
Evaluated UAER, as the primary symptom of DKD, leading to a decrease of GFR and proteinuria, is regarded as the essential criterion for diagnosing DKD in the early stage (Chen et al., 2017). UAER was reported in 15 trials (Hong, 2013; Li et al., 2013; Luo and Tao, 2014; Kang et al., 2014; Huang, 2015; Wang, 2015; Gao and Zhang, 2015; Liu et al., 2016; Yu et al., 2016; Deng and Yuan, 2016; Wang, 2016; Jiang, 2017; Yang, 2018; Zeng et al., 2018; Ma and Pang, 2018). A total of 1,288 patients were tested, 643 included in the treatment group and 645 in the control group. The results proved that UAER was reduced after CaD treatment compared with control groups (SMD = −1.29; 95% CI, −1.75 to −0.82; p < 0.00001; I2 = 93%).
3.1.4.2 24 h urinary protein/24 h urinary albumin
The data was reported in 12 trials (Wu et al., 2007; Sun, 2012; Zhang, 2012; Luo and Tao, 2014; Huang, 2015; Deng and Yuan, 2016; Liu et al., 2016; Zhu, 2016; Kuang, 2017; Deng, 2018; Yang, 2018; Zhang, 2008), which included a total of 891 patients, of which 448 were in the treatment group and 443 were in the control group. According to the results, we can see that the 24 h urinary protein/24 h urinary albumin of the CaD group is reduced compared with the control group (SMD = −1.95; 95% CI, −2.63 to −1.27; p < 0.00001; I2 = 94%)
3.1.4.3 α1-MG and β2-MG
α1-MG and β2-MG are tubular markers and could be predictive in diagnosing DKD patients due to their accuracy (Ferguson et al., 2015). Four trials (Hong, 2013; Zheng et al., 2013; Zhao et al., 2014; Liu et al., 2015) reported α1-MG. In trials targeting α1-MG, 284 patients were included, 142 in the treatment group and 142 in the control group. Control group patients were treated with ACEI, ARB, and routine treatment. The treatment duration was all less than or equal to 12 weeks. The result showed that the CaD group had a better curative effect (SMD = −2.32; 95% CI, −3.36 to −1.28; p < 0.0001; I2 = 91%). On the other hand, six trials (Wang et al., 2015; Jiang, 2016; Chen and Bai, 2017; Fan and Ma, 2017; Qin et al., 2017; Zhou and Chen, 2017) reported β2-MG. Seven hundred and thirty-six patients were included, half in the control and half in the CaD group. We conclude by analyzing the results: CaD can effectively improve the pathological changes of the kidney caused by DKD, making the MG in urine decrease significantly (SMD = −3.00; 95% CI, −4.41 to −1.58; p < 0.0001; I2 = 98%).
3.1.5 Endothelium function
Endothelin (ET) and nitric oxide (NO) were used to show endothelial function. Three trials (Wang et al., 2015; Wang, 2016; Kuang, 2017) reported NO, and four trials (Wang et al., 2015; Wang, 2016; Kuang, 2017; Zeng et al., 2018) reported ET. The treatment duration was all equal to or less than 12 weeks. The result showed that the CaD groups had a better effect on increasing NO (SMD = 0.68; 95% CI, 0.45 to 0.91; p < 0.00001; I2 = 0%) and decreasing ET (SMD = −0.81; 95% CI, −1.11 to −0.51; p < 0.00001; I2 = 51%) than the control groups.
3.1.6 Inflammation index
Inflammation response was one of the mechanisms related to DKD development. By modifying vascular permeability; releasing vasodilator and vasoconstrictor mediators; causing kidney fibrosis; and inducing cytotoxicity, apoptosis, and necrosis in the pathogenesis and progression of DKD, it may have an impact on the glomerular filtration function (Zhou et al., 2018). The enhanced inflammatory markers such as tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), and c-reactive protein (CRP) could be suggestive molecules in the development of diabetic vascular disease.
3.1.6.1 Hypersensitive c-reactive protein
Four trials (Chen and Bai, 2017; Fan and Ma, 2017; Qin et al., 2017; Jia et al., 2018) were collected to evaluate the Hs-CRP. Four hundred and ten patients were included in the total, half of which were in the CaD group and half in the control group. The results showed that CaD could effectively decrease the level of CRP compared with control groups (SMD = −1.43; 95% CI, −2.66 to −0.21; p = 0.02; I2 = 96%).
3.1.6.2 Interleukin-6
We collected two trials (Gao and Zhang, 2015; Qin et al., 2017) involving 160 DKD patients to assess the value of IL-6. In one trial, patients in the control group were treated with alprostadil, while in the other, the control group gave patients with telmisartan. According to this result, we can infer that no significant correlation was found between the CaD and the level of IL-6 (SMD = −6.07; 95% CI, −14.56 to 2.43; p = 0.16; I2 = 99%).
3.1.6.3 Tumor necrosis factor-α
Two studies (Qin et al., 2017; Jia et al., 2018) targeted the measure of TNF-α value. One hundred and ten patients were included in the total, 55 of which were in the control group. Although one trial (Qin et al., 2017) demonstrated a remarkable decrease in TNF-α in CaD patients compared to the control group, the other (Jia et al., 2018) did not exhibit a significant difference between the CaD group and control group (SMD = −8.07; 95% CI, −24.95 to 8.81; p = 0.35; I2 = 99%).
3.1.7 Hemodynamic index
Blood viscosity can reflect the state of hemodynamics, and alternation of hemodynamics is one of the critical pathways in DKD pathological development. We found two studies (Gao and Zhang, 2015; Ma and Pang, 2018) estimating the value of blood viscosity, and the result showed that CaD had no significant correlation with blood viscosity (SMD = −2.06; 95% CI, −5.35 to 1.23; p = 0.22; I2 = 98%).
3.2 Network pharmacology
3.2.1 Candidate targets of Calcium Dobesilate and diabetic kidney disease
For CaD, there are 293 related targets collected by searching the database. For DKD, we have retrieved 534 genes from the OMIM database, 1,678 genes from the Gene Cards database, and 1,189 genes from the DisGeNET database. After eliminating the redundancy, 2,771 known therapeutic targets were collected in this study. By drawing a Venn diagram to look for the intersection of CaD and DKD targets, 120 overlapping genes were obtained (Figure 3A).
FIGURE 3. Venn diagram and PPI work. (A) Venn diagram. The blue section indicates DKD-related targets, and the pink section indicates CaD-related targets. Twenty-five targets in the middle overlapping section are common targets of DKD and CaD. (B) PPI network. A total of 119 target proteins and 1,175 interacting edges are present in the network. Sizes and colors of the nodes are illustrated from big to small and blue to green in a descending order of degree values.
3.2.2 Protein–protein interaction network of diabetic kidney disease targets
To evaluate the role of potential targets in complicated diseases and discover connections, overlapping targets are submitted to the STRING11.0 platform to establish a PPI network. After analyzing the topology parameters of the PPI network, the 120 targets were sorted in descending order by degree and arranged in a concentric circle (Figure 3B). Degree reflected the importance of nodes by representing the number of connections between nodes and other nodes. According to the degree, the targets in the innermost circle were predicted as important targets.
3.2.3 Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment Analyses
We conducted GO functional analysis and KEGG pathway enrichment analysis to elucidate the biological impacts of CaD on gene functions and signaling pathways of relevant targets in treating DKD. As shown in Figure 4, the top 10 GO items and top 20 KEGG pathways were selected based on the p-value.
FIGURE 4. Enrichment analysis of the targets of CaD in treating DKD. (A) GO functional analysis. Top 10 items of each part are shown. (B) KEGG pathway enrichment analysis. The sizes of the bubbles are illustrated from big to small in a descending order of the number of potential targets involved in the pathways.
For biological processes, the targets were mainly enriched in the apoptotic signaling pathway, regulation of kinase activity, MAPK cascade, cellular response to chemical stress, and regulation of cytokine production. The top 10 items were mainly related to protein kinase activity, peptidase activity, and phosphotransferase activity (alcohol group as acceptor) in molecular functions. The cellular components were primarily concentrated in cytoplasmic vesicle lumen, secretory granule lumen, ficolin-1-rich granule lumen, membrane raft, vacuolar lumen, lytic vacuole, and lysosome; according to the results of the KEGG pathway enrichment analysis, most were involved in the MAPK signaling pathway, PI3K-Akt signaling pathway, the Ras signaling pathway, focal adhesion, fluid shear stress, atherosclerosis, Rap1 signaling pathway, and the chemokine signaling pathway. Interestingly, the typical targets of CaD and DKD are mainly enriched in the MAPK signaling pathway, and protein kinase activity was mentioned in the top 10 molecular function items. Additionally, biological processes involved apoptosis signaling pathways, regulation of kinase activity, and MAPK cascade. Therefore, we surmised that the MAPK signaling pathway might be critical. The number of targets concentrated in the chemokine pathway was relatively large in the figure, and the biological process highlighted cellular response to chemical stress and regulation of cytokine production. Meanwhile, the meta-analysis results showed that CaD significantly reduced the level of inflammatory factors. These imply that the chemokine signaling pathway cannot be overlooked.
3.2.4 Component–target–pathway network construction
A component–target–pathway network was constructed with Cytoscape3.7.1 based on the KEGG pathway enrichment analysis (Figure 5). Among the essential targets selected by the PPI network, AKT1, CASP3, IGF1, and MAPK8 were enriched in the MAPK signaling pathway, and CCL5, a typical inflammatory chemokine, was increased in the chemokine signaling pathway at the same time. We speculated that the five targets might be critical for CaD in treating DKD and used for further molecular docking.
FIGURE 5. Component–target–pathway network. A total of 77 nodes and 372 edges are present in the network. Orange diamond represents the bioactive component of CaD, 57 green squares represent targets, and 20 blue V-shapes represent pathways. Sizes of the green square node are illustrated from big to small in a descending order of degree values. A total of 372 edges represent the interaction relationship between components, targets, and pathways.
3.2.5 Molecular docking result analysis
In the present studies, the possible interaction activity between five key targets and their corresponding compounds of CaD was investigated with molecular docking verification. Among the docking results, most binding complexes possessed high binding affinity, averaging −5.12 kcal/mol. The modes of five binding complexes are displayed in Figure 6, including CaD-AKT1 docking (−5 kcal/mol), CaD-CASP3 docking (−5.55 kcal/mol), CaD-IGF1 docking (−4.62 kcal/mol), CaD-MAPK8 docking (−5.16 kcal/mol), and CaD-CCL5 docking (−5.26 kcal/mol). For concreteness, taking the CaD-AKT1 docking as an example, small-molecule ligand CaD may be embedded in the interfaced pocket formed by the interaction of amino acid residues in the protein (Figure 6Aa). Figure 6Ab shows three hydrogen bond formations between ligand and residues in SER 216, LEU202, and VAL201. The other essential residues (LYS214, TYR215, LEU213, GLN203, and ASN204) interacted with CaD through van der Waals forces. These forms of hydrogen bonds and interactions contribute to the stability of the binding of small molecules to the active sites of proteins.
FIGURE 6. Molecular docking diagram. (A) Five conformations of a molecular docking simulation. Diagrams (3D) represent that the molecular model of the compound is in the binding pocket of the protein. The compound is shown as a stick model in orange. The amino acid residues in the surrounding are represented by surface style. Diagrams (2D) show the interactions between the compound and surrounding residues. (B) 3D column diagram shows the affinity of six conformations. X-axis: bioactive component, y-axis: target names, and z-axis: docking affinity (absolute value of the binding energy).
4 Discussion
We gathered 42 trials of CaD intervention for DKD clinical treatment, including 3671 DKD patients in various phases of the disease, and divided them into the CaD experimental and control groups. The results showed that Scr, BUN, and Cys-C levels in the blood, as well as molecules including UAER, 24 h urine protein/24 h urine albumin, 1-MG, and 2-MG in the urine, were significantly lower in the CaD group compared to the control group, while GFR levels were substantially higher. These findings suggested that CaD could help DKD patients improve their glomerular filtration performance and kidney function. In addition, CaD could modulate endothelium contraction and relaxation by boosting NO and lowering ET, regulating microvasculature function. Inflammatory variables such as CRP levels also dropped in the CaD group, implying that CaD can help DKD patients lessen their inflammatory response. Compared to the control group, patients in the CaD group showed a trend of improvement in numerous indices, including glomerular filtering performance, endothelium function, and inflammatory function. As a result, we hypothesize that CaD may bind to multiple targets in vivo, activating or inhibiting multiple metabolic and disease-related pathways, reminding us that a network pharmacological approach can be used to investigate the mechanisms by which drugs improve multiple metrics in patients.
Then, using network pharmacology, we built a PPI network for the drug’s and disease’s shared targets and ran GO and KEGG analyses on these genes. According to BP enrichment, multiple genes are abundant in apoptosis, regulation of kinase activity, and cellular response to chemical stress, which could be CaD’s primary approach to suppressing the development of DKD. Based on the KEGG pathway enrichment analysis and literature research, we anticipated that CaD would play a therapeutic function in DKD primarily by modulating the MAPK signaling pathway and the chemokine signaling pathway (Figure 7). We identified AKT1, CASP3, IGF1, MAPK8, and CCL5 as primary targets based on the component–target–pathway network analysis, implying that the CaD-mediated DKD treatment is mainly related to the abovementioned targets.
FIGURE 7. MAPK and chemokine signaling pathways are influenced by dapagliflozin. The red nodes represent the key targets, the orange nodes represent common targets of CaD and DKD targets, and the green nodes represent the other targets of these two pathways. CaD affects the phosphorylation of MAPK14 in the p38/MAPK pathway. In the JNK pathway, CaD affects the phosphorylation of MAPK8 and MAPK10 and it also indirectly affects the activation of CASP3. As for the classical MAPK signaling pathway, CaD affects the activation of GRB2, HRAS, RAF1, and the phosphorylation of MEK and ERK. The chemokine signaling pathway is closely related to the classical MAPK signaling pathway. CaD inhibits the expression of CCL5. CaD inhibits apoptosis, oxidative stress, and inflammation by inhibiting the MAPK signaling pathway and the chemokine signaling pathway, thereby exerting a therapeutic effect on DKD.
Reactive oxygen species (ROS), the products of normal metabolism and xenobiotic exposure, are significant factors related to the pathogenesis of DKD (Bahmani et al., 2016; Miranda-Díaz et al., 2016), and their production can induce apoptosis (Xu et al., 2017). CaD has been shown to prevent cell apoptosis by lowering the number of ROS (Iriz et al., 2008). CASP3 is an essential component of apoptosis and is closely associated with the progression of DKD (Zhou et al., 2021). Hyperglycemia stimulates CASP3 cleavage and DNA breakage, and the resulting apoptosis leads to mesenchymal cell loss in DKD. (Mishra et al., 2005). One of the earliest events in DKD is glomerular thickening resulting from mesangial cell hypertrophy, and it has been determined that Akt signaling contributes to thylakoid hypertrophy in DKD. (Mahimainathan et al., 2006). DKD is associated with the loss of renal cells, particularly glomerular podocytes, which form the glomerular filtration barrier, and changes in Akt signaling as a critical event in podocyte loss during early DKD. In addition, inhibition of Akt and its downstream targets such as mTOR may provide future therapeutic benefits for treating DKD (Heljić and Brazil, 2011). IGF1 is an important growth factor that maintains the structure and function of nephritis and plays a vital role in the pathology of DKD (Bach and Hale, 2015). Previous studies have shown that IGF1 overexpression causes many histopathology changes such as kidney tissue hyperplasia, renal cell proliferation, nephromegaly, mesangial expansion, and increased inflammatory cytokines (Li et al., 2018). Donath also found that inhibition of IGF1R could reduce inflammation in DKD more effectively (Donath, 2013). Animal and human kidney biopsy studies have shown that activation of stress-activated protein kinases (p38 MAPK and JNK) is associated with the progression of inflammation and injury in multiple forms of kidney disease (Adhikary et al., 2004). In addition, pharmacological inhibitors of p38 MAPK or c-Jun N-terminal kinase (JNK) are effective in animal models of renal disease when used as prophylactic agents to prevent injury development or as interventional therapies to inhibit the progression of established injury (Tesch et al., 2016). These findings support the notion that p38 MAPK and JNK signaling are important therapeutic targets for preventing kidney damage. Besides, CCL5 is involved in the pathogenesis of diabetic kidney injury as an inflammatory chemokine upregulated in response to the metabolic and hemodynamic characteristics of the diabetic environment (Pérez-Morales et al., 2019). Subsequently, the desired molecular docking results strongly prove the criticality of the above five targets.
Podocytes are terminally differentiated cells, and podocyte damage is the critical event leading to proteinuria in DKD (Miranda-Díaz et al., 2016). ROS is essential for initiating podocyte apoptosis (Zhang et al., 2015). Activating p38 MAPK, a pro-apoptotic signaling factor downstream of ROS leads to cell apoptosis (Finkel and Holbrook, 2000). Evidence has shown that the increase of ROS could activate profibrotic factors, including TGF-β through P38-MAPK, which could promote the synthesis of type IV collagen (the main component of extracellular matrix (ECM)) and fibrin connection protein, thus causing an increase in ECM, and forming early DKD (Sakai et al., 2005). Previous studies have shown that leptin can signal through the leptin receptor isoform to stimulate glomerular endothelial cell proliferation, increase TGF-β1 synthesis, and type IV collagen production (Wolf and Ziyadeh, 2006). Leptin promotes fibrosis primarily on the glomerulus but can potentially prevent/reverse renal injury by normalizing metabolic disturbances, including hyperglycemia and hyperlipidemia (Suganami et al., 2005). Bulent O et al. suggested the potential importance of leptin in regulating glucose homeostasis and its possible direct application in treating disorders of glucose homeostasis (Yildiz and Haznedaroglu, 2006). These also provide help to broaden our treatment ideas. Furthermore, the MAPK pathway can be activated by chronic hyperglycemia, resulting in a local inflammatory response (Fang et al., 2012). Targeted therapy for inhibiting the p38 MAPK signaling pathway has shown preventive effects on streptozotocin-induced DKD (Lin et al., 2018). These findings indicate that CaD may significantly influence DKD by inhibiting the MAPK signaling pathway.
In addition tocell apoptosis, chemokine production has also been believed to play an essential role in the process of DKD (Dieter et al., 2019). The proinflammatory chemokine (C–C motif) ligand 12 caused glomerular sclerosis in T2D mice, and blocking its expression had a protective effect on DKD (Darisipudi et al., 2011). Inflammatory CCL5 is expressed in various cell types, including fibroblasts and renal tubular epithelial cells. Previously, upregulated CCL5 was found in the kidney, and its expression is directly related to the proteinuria concentration in kidney tubular cells (Navarro-González et al., 2011; Zhang et al., 2016). Consequently, it can be concluded from the existing research that the chemokine signaling pathway plays a vital role in CaD in treating DKD.
Although our results revealed the advantages of CaD for DKD, limitations did remain in the analysis. In the present study, we chose serum creatinine as a marker to assess renal function. At the same time, the literature has reported that CaD interferes to some extent with serum creatinine measured by sarcosine oxidase, both in vivo and in vitro (Zhou et al., 2018). However, we also evaluated serum Cys-C, which may be more sensitive for assessing renal function. Its levels are not interfered with by CaD, and many studies suggest using serum Cys-C to evaluate renal function in patients receiving CaD therapy (Guo et al., 2015). Other studies point out that microalbuminuria is the current gold standard for predicting and detecting diabetic kidney disease (Salem et al., 2020). The evaluation of urine albumin excretion rate (UAER) and 24 h urine protein/24 h urine albumin are covered in our article. Therefore, this issue is not sufficient to affect our findings.
5 Conclusion
In conclusion, our research systematically elucidated the underlying molecular mechanisms by which CaD interfered with DKD based on the meta-analysis, network pharmacology, and molecular docking. We predicted five key targets from complex networks and concluded that CaD exerts a therapeutic effect on DKD by inhibiting MAPK and the chemokine signaling pathway. We expect this research to provide additional reference directions for CaD as a medicine in treating diabetic microangiopathy. Further exploration can be done in the future based on this study.
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 authors.
Author contributions
ZZ and SY contributed to the study conception and design, literature review, and preparation of the manuscript. YW, HF, HS, and ZY contributed to the study conception and design. BD and YY drafted the manuscript, revised it critically for important intellectual content, and gave approval of the final version. Both authors read and approved the final manuscript.
Funding
The authors thank the Natural Science Foundation of China (No. 82000788), Chinese Postdoctoral Science Foundation (No. 2021M702040), Key Research and Development Plan of Shandong Province (No. 2018GSF118176), Guangdong Basic and Applied Basic Research Foundation (No. 2019A1515110023), Natural Science Foundation of Shandong Province (No. ZR2016HQ26), the Bethune-Merck’s Diabetes Research Foundation (Nos. B-0307-H-20200302 and G2016014), Medicine and Health Science and Technology Development Plan of Shandong Province (No. 2018WS201), and Science and technology plan of Shandong Provincial University (No. J16LK09) for funding.
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.
Footnotes
2http://www.bioinformatics.com.cn/.
3http://www.ehbio.com/ImageGP/.
References
Adhikary, L., Chow, F., Nikolic-Paterson, D. J., Stambe, C., Dowling, J., Atkins, R. C., et al. (2004). Abnormal p38 mitogen-activated protein kinase signalling in human and experimental diabetic nephropathy. Diabetologia 47 (7), 1210–1222. doi:10.1007/s00125-004-1437-0
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
Ashraf, H., Lotfi, M., Akbari, M., Jahanbani-Ardakani, H., Ghaffari, M., and Khalili, M. R. (2021). Effect of calcium dobesilate on retrobulbar blood flow and choroidal thickness in patients with non-proliferative diabetic retinopathy. Int. Ophthalmol. 41 (10), 3479–3486. doi:10.1007/s10792-021-01914-x
Bach, L. A., and Hale, L. J. (2015). Insulin-like growth factors and kidney disease. Am. J. Kidney Dis. 65 (2), 327–336. doi:10.1053/j.ajkd.2014.05.024
Bahmani, F., Kia, M., Soleimani, A., Mohammadi, A. A., and Asemi, Z. (2016). The effects of selenium supplementation on biomarkers of inflammation and oxidative stress in patients with diabetic nephropathy: A randomised, double-blind, placebo-controlled trial. Br. J. Nutr. 116 (7), 1222–1228. doi:10.1017/S0007114516003251
Berman, H. M., Westbrook, J., Feng, Z., Gilliland, G., Bhat, T. N., Weissig, H., et al. (2000). The protein Data Bank. Nucleic Acids Res. 28 (1), 235–242. doi:10.1093/nar/28.1.235
Besirli, K., Aydemir, B., Arslan, C., Kiziler, A. R., Canturk, E., and Kayhan, B. (2012). Calcium dobesilate may improve hemorheology in patients undergoing coronary artery bypass grafting. Rev. Bras. Cir. Cardiovasc. 27 (2), 260–266. doi:10.5935/1678-9741.20120057
Boezio, B., Audouze, K., Ducrot, P., and Taboureau, O. (2017). Network-based approaches in pharmacology. Mol. Inf. 36 (10), 1700048. doi:10.1002/minf.201700048
Cai, T., Wu, X. Y., Zhang, X. Q., Shang, H. X., Zhang, Z. W., Liao, L., et al. (2017). Calcium dobesilate prevents diabetic kidney disease by decreasing Bim and inhibiting apoptosis of renal proximal tubular epithelial cells. DNA Cell Biol. 36 (4), 249–255. doi:10.1089/dna.2016.3276
Cen, Q., Li, H. Q., Mo, J. T., Liu, P., and Zhou, Z. N. (2016). Clinical observation of calcium dobesilate combined with irbesartan in the treatment of early diabetic kidney disease. J. Guangxi Univ. Chin. Med. 1, 29–31.
Chen, C., Wang, C., Hu, C., Han, Y., Zhao, L., and Zhu, X. (2017). Normoalbuminuric diabetic kidney disease. Front. Med. 11 (3), 310–318. doi:10.1007/s11684-017-0542-7
Chen, P., and Bai, X. (2017). Clinical observation of benazepril combined with calcium dobesilate in the treatment of senile early diabetic kidney disease. China Dispens. 28 (3), 347–350.
Cheng, J., Liu, Y., Li, G., Zhang, Z., Ma, L., Yang, X., et al. (2016). Noninvasive ventilation improves cardiac function in patients with chronic heart failure. Oncotarget 7 (31), 48918–48924. doi:10.18632/oncotarget.10441
Ciapponi, A., Laffaire, E., and Roqué, M. (2004). Calcium dobesilate for chronic venous insufficiency: A systematic review. Angiology 55 (2), 147–154. doi:10.1177/000331970405500206
Darisipudi, M. N., Kulkarni, O. P., Sayyed, S. G., Ryu, M., Migliorini, A., Sagrinati, C., et al. (2011). Dual blockade of the homeostatic chemokine CXCL12 and the proinflammatory chemokine CCL2 has additive protective effects on diabetic kidney disease. Am. J. Pathol. 179 (1), 116–124. doi:10.1016/j.ajpath.2011.03.004
Deng, L. L. (2018). Safety and efficacy of calcium dobesilate in the treatment of diabetic kidney disease. Mod. Diagnosis Treat. 12, 1926–1927. CNKI: SUN: XDZD.0.2018-12-040.
Deng, X. Z., and Yuan, J. Y. (2016). Clinical study of shenyan kangfu tablet combined with calcium dobesilate in the treatment of diabetic kidney disease. Mod. Diagnosis Treat. 16, 3003–3004. CNKI: SUN: XDZD.0.2016-16-039.
Dharnidharka, V. R., Kwon, C., and Stevens, G. (2002). Serum cystatin C is superior to serum creatinine as a marker of kidney function: A meta-analysis. Am. J. Kidney Dis. 40 (2), 221–226. doi:10.1053/ajkd.2002.34487
Dieter, B. P., Meek, R. L., Anderberg, R. J., Cooney, S. K., Bergin, J. L., Zhang, H., et al. (2019). Serum amyloid A and Janus kinase 2 in a mouse model of diabetic kidney disease. PLoS One 14 (2), e0211555. doi:10.1371/journal.pone.0211555
Donath, M. Y. (2013). Targeting inflammation in the treatment of type 2 diabetes. Diabetes Obes. Metab. 15 (3), 193–196. doi:10.1111/dom.12172
Fan, H. H., and Ma, J. (2017). Clinical effect of calcium dobesilate combined with insulin in the treatment of early diabetic kidney disease in the elderly. Chin. J. Prim. Med. Pharm. 24, 3746–3749. doi:10.3760/cma.j.issn.1008-6706.2017.24.018
Fang, Y., Tian, X., Bai, S., Fan, J., Hou, W., Tong, H., et al. (2012). Autologous transplantation of adipose-derived mesenchymal stem cells ameliorates streptozotocin-induced diabetic nephropathy in rats by inhibiting oxidative stress, pro-inflammatory cytokines and the p38 MAPK signaling pathway. Int. J. Mol. Med. 30 (1), 85–92. doi:10.3892/ijmm.2012.977
Ferguson, T. W., Komenda, P., and Tangri, N. (2015). Cystatin C as a biomarker for estimating glomerular filtration rate. Curr. Opin. Nephrol. Hypertens. 24 (3), 295–300. doi:10.1097/MNH.0000000000000115
Finkel, T., and Holbrook, N. J. (2000). Oxidants, oxidative stress and the biology of ageing. Nature 408 (6809), 239–247. doi:10.1038/35041687
Gao, Z., and Zhang, Y. C. (2015). Clinical effect of telmisartan combined with calcium dobesilate in the treatment of senile diabetic kidney disease with hypertension. China Med. Her. 24, 136–138+145. CNKI: SUN: YYCY.0.2015-24-034.
Goodsell, D. S., and Olson, A. J. (1990). Automated docking of substrates to proteins by simulated annealing. Proteins 8 (3), 195–202. doi:10.1002/prot.340080302
Guo, X., Hou, L., Cheng, X., Zhang, T., Yu, S., Fang, H., et al. (2015). Strong negative interference by calcium dobesilate in sarcosine oxidase assays for serum creatinine involving the trinder reaction. Med. Baltim. 94 (23), e905. doi:10.1097/MD.0000000000000905
Heberle, H., Meirelles, G. V., da Silva, F. R., Telles, G. P., and Minghim, R. (2015). InteractiVenn: A web-based tool for the analysis of sets through Venn diagrams. BMC Bioinforma. 16 (1), 169. doi:10.1186/s12859-015-0611-3
Heljić, M., and Brazil, D. P. (2011). Protein kinase B/Akt regulation in diabetic kidney disease. Front. Biosci. 3 (1), 98–104. doi:10.2741/s135
Hong, Y. (2013). Effect of calcium dobesilate combined with olmesartan on proteinuria in diabetic kidney disease. Mil. Med. J. South China 12, 877–879. CNKI: SUN: HNGY.0.2013-12-005.
Huang, C. L. (2015). Clinical efficacy and safety of calcium dobesilate in the treatment of diabetic nephropathy. Chin. J. Clin. Pharmacol. 13, 1239–1241. doi:10.13699/j.cnki.1001-6821.2015.13.006
Iriz, E., Vural, C., Ereren, E., Poyraz, A., Erer, D., Oktar, L., et al. (2008). Effects of calcium dobesilate and diosmin-hesperidin on apoptosis of venous wall in primary varicose veins. Vasa. 37 (3), 233–240. doi:10.1024/0301-1526.37.3.233
Jia, M., Xie, Y. X., and Qiu, H. (2018). Effect of calcium dobesilate on hs-CRP, TNF- α in patients with stage 3diabetic kidney disease. Med. Innovation China 7, 5–8. CNKI: SUN: ZYCX.0.2018-07-004.
Jiang, H. (2016). Effect of calcium dobesilate combined with losartan on diabetic kidney disease. Health Res. 5, 569–570+573. CNKI: SUN: HZYG.0.2016-05-031.
Jiang, M. C. (2017). Therapeutic effect of calcium dobesilate on early diabetic kidney disease. Electron. J. Clin. Med. Literature 13, 2383–2392. doi:10.16281/j.cnki.jocml.2017.13.005
Kang, H., Ma, D., and Li, Y. (2014). Clinical observation of calcium dobesilate in the treatment of early diabetic kidney disease. Chin. J. Mod. Drug Appl. 1, 143–144. doi:10.14164/j.cnki.cn11-5581/r.2014.01.119
Kim, S., Chen, J., Cheng, T., Gindulyte, A., He, J., He, S., et al. (2021). PubChem in 2021: New data content and improved web interfaces. Nucleic Acids Res. 49 (D1), D1388–D1395. doi:10.1093/nar/gkaa971
Kuang, J. X. (2017). Study on the clinical therapeutic effect of calcium dobesilate on diabetic kidney disease. Strait Pharm. J. 04, 137–138. CNKI: SUN: HAIX.0.2017-04-069.
Levey, A. S., Coresh, J., Tighiouart, H., Greene, T., and Inker, L. A. (2020). Measured and estimated glomerular filtration rate: Current status and future directions. Nat. Rev. Nephrol. 16 (1), 51–64. doi:10.1038/s41581-019-0191-y
Li, J., Dong, R., Yu, J., Yi, S., Da, J., Yu, F., et al. (2018). Inhibitor of IGF1 receptor alleviates the inflammation process in the diabetic kidney mouse model without activating SOCS2. Drug Des. devel. Ther. 12, 2887–2896. doi:10.2147/DDDT.S171638
Li, L. B., Qin, A. P., Zhao, X. L., Li, W. J., Chen, K., and Tang, Z. (2013). Therapeutic effect of calcium dobesilate on nephropathy in the early aged patients with type 2 diabetes mellitus. Chin. J. Clin. Pharmacol. 12, 900–902. doi:10.13699/j.cnki.1001-6821.2013.12.011
Li, X. N. (2018). Clinical efficacy and safety of calcium dobesilate combined with candesartan in treatment of phase III diabetic kidney disease. J. Huaihai Med. 6, 649–650+654. doi:10.14126/j.cnki.1008-7044.2018.06.007
Lin, Y. C., Chang, Y. H., Yang, S. Y., Wu, K. D., and Chu, T. S. (2018). Update of pathophysiology and management of diabetic kidney disease. J. Formos. Med. Assoc. 117 (8), 662–675. doi:10.1016/j.jfma.2018.02.007
Liu, H. Z., Bai, J. M., and Zhao, P. (2015). Therapeutic effect of calcium dobesilate combined with enalapril on early type 2diabetic kidney disease in elderly patients. J. Pract. Diabetology 4, 22–24. CNKI: SUN: LNSY.0.2015-04-014.
Liu, J., Li, S., and Sun, D. (2019). Calcium dobesilate and micro-vascular diseases. Life Sci. 221, 348–353. doi:10.1016/j.lfs.2019.02.023
Liu, L. H., Peng, B., and Xu, Y. (2016). Clinical observation of calcium dobesilate combined with candesartan in the treatment of early diabetic kidney disease. Int. J. Urology Nephrol. 03, 408–410. CNKI: SUN: GYMY.0.2016-03-029.
Luo, S. H., and Tao, L. (2014). Therapeutic effect of calcium dobesilate combined with telmisartan on diabetic kidney disease. Zhejiang Med. J. 19, 1629–1631. CNKI: SUN: ZJYE.0.2014-19-017.
Ma, P., and Pang, X. F. (2018). Effect of irbesartan and calcium dobesilate on kidney hemodynamics in patients with early type 2 diabetes mellitus complicated with hypertension. Hainan Med. J. 8, 1055–1059. CNKI: SUN: HAIN.0.2018-08-006.
Mahimainathan, L., Das, F., Venkatesan, B., and Choudhury, G. G. (2006). Mesangial cell hypertrophy by high glucose is mediated by downregulation of the tumor suppressor PTEN. Diabetes 55 (7), 2115–2125. doi:10.2337/db05-1326
Michal, M., and Giessinger, N. (1985). Effect of calcium dobesilate and its interaction with aspirin on thrombus formation in vivo. Thromb. Res. 40 (2), 215–226. doi:10.1016/0049-3848(85)90332-9
Miranda-Díaz, A. G., Pazarín-Villaseñor, L., Yanowsky-Escatell, F. G., and Andrade-Sierra, J. (2016). Oxidative stress in diabetic nephropathy with early chronic kidney disease. J. Diabetes Res. 2016, 7047238. doi:10.1155/2016/7047238
Mishra, R., Emancipator, S. N., Kern, T., and Simonson, M. S. (2005). High glucose evokes an intrinsic proapoptotic signaling pathway in mesangial cells. Kidney Int. 67 (1), 82–93. doi:10.1111/j.1523-1755.2005.00058.x
Navarro-González, J. F., Mora-Fernández, C., Muros de Fuentes, M., and García-Pérez, J. (2011). Inflammatory molecules and pathways in the pathogenesis of diabetic nephropathy. Nat. Rev. Nephrol. 7 (6), 327–340. doi:10.1038/nrneph.2011.51
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
Page, M. J., McKenzie, J. E., Bossuyt, P. M., Boutron, I., Hoffmann, T. C., Mulrow, C. D., et al. (2021). The PRISMA 2020 statement: An updated guideline for reporting systematic reviews. BMJ 372, n71. doi:10.1136/bmj.n71
Pérez-Morales, R. E., Del Pino, M. D., Valdivielso, J. M., Ortiz, A., Mora-Fernández, C., and Navarro-González, J. F. (2019). Inflammation in diabetic kidney disease. Nephron 143 (1), 12–16. doi:10.1159/000493278
Piñero, J., Bravo, À., Queralt-Rosinach, N., Gutiérrez-Sacristán, A., Deu-Pons, J., Centeno, E., et al. (2017). DisGeNET: A comprehensive platform integrating information on human disease-associated genes and variants. Nucleic Acids Res. 45 (D1), D833–D839. doi:10.1093/nar/gkw943
Qin, L. L., Qin, W. J., Wang, J. F., and Lin, L. (2017). Combined treatment of diabetic nephropathy with alprostadil and calcium dobesilate. Exp. Ther. Med. 14 (5), 5012–5016. doi:10.3892/etm.2017.5115
Sakai, N., Wada, T., Furuichi, K., Iwata, Y., Yoshimoto, K., Kitagawa, K., et al. (2005). Involvement of extracellular signal-regulated kinase and p38 in human diabetic nephropathy. Am. J. Kidney Dis. 45 (1), 54–65. doi:10.1053/j.ajkd.2004.08.039
Salem, N. A., El Helaly, R. M., Ali, I. M., Ebrahim, H. A. A., Alayooti, M. M., El Domiaty, H. A., et al. (2020). Urinary Cyclophilin A and serum Cystatin C as biomarkers for diabetic nephropathy in children with type 1 diabetes. Pediatr. Diabetes 21 (5), 846–855. doi:10.1111/pedi.13019
Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303
Stelzer, G., Rosen, N., Plaschkes, I., Zimmerman, S., Twik, M., Fishilevich, S., et al. (2016). The GeneCards suite: From gene data mining to disease genome sequence analyses. Curr. Protoc. Bioinforma. 54, 1.30.1–1.30.33. doi:10.1002/cpbi.530.33
Suganami, T., Mukoyama, M., Mori, K., Yokoi, H., Koshikawa, M., Sawai, K., et al. (2005). Prevention and reversal of renal injury by leptin in a new mouse model of diabetic nephropathy. FASEB J. 19 (1), 127–129. doi:10.1096/fj.04-2183fje
Sun, L. J. (2012). Clinical observation on 81 cases of diabetic kidney disease treated with calcium dobesilate capsule and alprostadil. China Pract. Med. 35, 5–7. doi:10.14163/j.cnki.11-5547/r.2012.35.141
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 (D1), D607–D613. doi:10.1093/nar/gky1131
Tesch, G. H., Ma, F. Y., and Nikolic-Paterson, D. J. (2016). ASK1: A new therapeutic target for kidney disease. Am. J. Physiol. Ren. Physiol. 311 (2), F373–F381. doi:10.1152/ajprenal.00208.2016
UniProt Consortium (2019). UniProt: A worldwide hub of protein knowledge. Nucleic Acids Res. 47 (D1), D506–D515. doi:10.1093/nar/gky1049
Wang, J. (2016). Therapeutic effect of calcium dobesilate on early diabetic kidney disease. China Pract. Med. 22, 133–135. doi:10.14163/j.cnki.11-5547/r.2016.22.098
Wang, P., Wan, S., Du, X., Xue, Y., and Gao, H. (2015). Clinical study of calcium dobesilate dispersible tablets combined with α-lipoic acid in the treatment of diabetic kidney disease stage Ⅳ. Chin. J. Integr. Chin. West. Med. Nephrop. 16 (10), 898–900. CNKI: SUN: JXSB.0.2015-10-023.
Wang, X., Shen, Y., Wang, S., Li, S., Zhang, W., Liu, X., et al. (2017). PharmMapper 2017 update: A web server for potential drug target identification with a comprehensive target pharmacophore database. Nucleic Acids Res. 45 (W1), W356–W360. doi:10.1093/nar/gkx374
Wang, Y. (2015). Clinical effect of calcium and hydrochloric acid in the treatment of diabetic kidney disease patients with benner pury. China J. Pharm. Econ. 10, 35–36. CNKI: SUN: ZYWA.0.2015-10-011.
Wolf, G., and Ziyadeh, F. N. (2006). Leptin and renal fibrosis. Contrib. Nephrol. 1 51, 175–183. doi:10.1159/000095328
Wu, Q., Lei, Q. F., Zhou, Z. J., Jiang, H., and Cheng, J. Z. (2007). Clinical observation of calcium dobesilate in the treatment of chronic renal failure in diabetic kidney disease. J. Yangtze Univ. Nat. Sci. Ed. Med. V 1, 42–43+45. doi:10.3969/j.issn.1673-1409-B.2007.01.020
Xu, M., Sun, B., Li, D., Mao, R., Li, H., Li, Y., et al. (2017). Beneficial effects of small molecule oligopeptides isolated from panax ginseng meyer on pancreatic beta-cell dysfunction and death in diabetic rats. Nutrients 9 (10), 1061. doi:10.3390/nu9101061
Yang, L. (2018). Clinical efficacy and safety of calcium hydroxy benzenesulfonate in the treatment of diabetic kidney disease. China Rural. Health 10, 35. CNKI: SUN: NCWS.0.2018-10-035.
Yildiz, B. O., and Haznedaroglu, I. C. (2006). Rethinking leptin and insulin action: Therapeutic opportunities for diabetes. Int. J. Biochem. Cell Biol. 38 (5-6), 820–830. doi:10.1016/j.biocel.2005.09.013
Yu, R., Liang, G. Z., and Ji, C. Y. (2016). Effect of calcium dobesilate combined with losartan on elderly patients with type 2 diabetic kidney disease complicated with hypertension. Drugs & Clin. 2, 170–173. doi:10.7501/j.issn.1674-5515.2016.02.010
Zeng, L. J., Huang, Z. X., and Luan, L. L. (2018). Therapeutic effect of calcium dobesilate combined with valsartan on early diabetic kidney disease and its effect on serum ET and CysC levels. Hainan Med. J. 8, 1064–1066. doi:10.3969/j.issn.1003-6350.2018.08.008
Zhang, H., Guo, S. H., Xue, Z. K., Zhang, Y. R., Wang, J. R., Che, J. J., et al. (2021). Effect of calcium dobesilate in preventing contrast-induced nephropathy in patients with diabetes and chronic kidney disease. Clin. Sao Paulo. 76, e2942. doi:10.6061/clinics/2021/e2942
Zhang, L., Ji, L., Tang, X., Chen, X., Li, Z., Mi, X., et al. (2015). Inhibition to DRP1 translocation can mitigate p38 MAPK-signaling pathway activation in GMC induced by hyperglycemia. Ren. Fail. 37 (5), 903–910. doi:10.3109/0886022X.2015.1034607
Zhang, L., Li, Y. M., Wang, Z. F., Xing, X. Z., Zhu, L., Jiang, J. X., et al. (2017). Protective effect of calcium dobesilate combined with benazepril therapy on renal injury in patients with early diabetic kidney disease and the possible molecular mechanisms. J. hainan Med. Univ. 23 (6), 762–765. doi:10.13210/j.cnki.jhmu.20161223.007
Zhang, L. (2012). Observation on combined use of drugs in 64 cases of early diabetic kidney disease. Chin. J. Physicians' study 4, 57–58. doi:10.3760/cma.j.issn.1673-4904.2012.04.025
Zhang, M. X., Liu, S., Tan, D. M., and Tan, Z. H. (2013). Therapeutic effect of calcium dobesilate combined with Fosinopril on diabetic kidney disease in clinical stage. Chin. J. Gerontology 4, 927–928. doi:10.3969/j.issn.1005-9202.2013.04.091
Zhang, X. L. (2008). Clinical observation of calcium dobesilate in the treatment of 76 cases of early diabetic kidney disease. Med. J. Chin. People's Health 19, 2262–2264. CNKI: SUN: ZMYX.0.2008-19-022.
Zhang, Z., Zhang, X., Dong, J., Gao, W., Liu, F., Zhao, J., et al. (2016). Association of chemokine ligand 5/chemokine receptor 5 gene promoter polymorphisms with diabetic microvascular complications: A meta-analysis. J. Diabetes Investig. 7 (2), 212–218. doi:10.1111/jdi.12397
Zhao, P., Liu, H. Z., and Bai, J. M. (2014). Effect of calcium dobesilate combined with benazepril on early aged type 2 diabetic kidney disease. China J. Mod. Med. 34, 73–76. doi:10.3969/j.issn.1005-8982.2014.34.018
Zheng, X., Hu, Z. X., and Zhang, Y. H. (2013). Effect of calcium dobesilate combined with enalapril on proteinuria in diabetic kidney disease. Sichuan Med. J. 8, 1249–1250. doi:10.16252/j.cnki.issn1004-0501-2013.08.058
Zhou, R. Q., and Chen, Y. (2017). Efficacy of irbesartan combined with calcium dobesilate in the treatment of 30 cases of diabetic kidney disease. Shanghai Med. Pharm. J. 5, 15–17. doi:10.3969/j.issn.1006-1533.2017.05.006
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 (2), 98–113. doi:10.2478/jtim-2021-0020
Zhou, Y., Qi, C., Li, S., Shao, X., Mou, S., and Ni, Z. (2018). Diabetic nephropathy can Be treated with calcium dobesilate by alleviating the chronic inflammatory state and improving endothelial cell function. Cell. Physiol. biochem. 51 (3), 1119–1133. doi:10.1159/000495491
Zhou, Y., Zhou, B., Pache, L., Chang, M., Khodabakhshi, A. H., Tanaseichuk, O., et al. (2019). Metascape provides a biologist-oriented resource for the analysis of systems-level datasets. Nat. Commun. 10 (1), 1523. doi:10.1038/s41467-019-09234-6
Keywords: calcium dobesilate, diabetic kidney disease, network pharmacology, MAPK signaling pathway, chemokine signaling pathway
Citation: Du B, Yin Y, Wang Y, Fu H, Sun H, Yue Z, Yu S and Zhang Z (2022) Calcium dobesilate efficiency in the treatment of diabetic kidney disease through suppressing MAPK and chemokine signaling pathways based on clinical evaluation and network pharmacology. Front. Pharmacol. 13:850167. doi: 10.3389/fphar.2022.850167
Received: 18 January 2022; Accepted: 11 August 2022;
Published: 08 September 2022.
Edited by:
Jianping Chen, Shenzhen Traditional Chinese Medicine Hospital, ChinaReviewed by:
Norberto Perico, Mario Negri Pharmacological Research Institute (IRCCS), ItalyDingkun Gui, Shanghai Jiao Tong University, China
Ibrahim C. Haznedaroglu, Hacettepe University Hospital, Turkey
Copyright © 2022 Du, Yin, Wang, Fu, Sun, Yue, Yu and Zhang. 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: Shaohong Yu, sutcm2006@163.com; Zhongwen Zhang, zhangzhongwen@sdu.edu.cn
†These authors have contributed equally to this work and share first authorship