- 1Department of Endocrinology, Shandong Provincial Hospital Affiliated to Shandong First Medical University, Jinan, China
- 2Shandong Key Laboratory of Endocrinology and Lipid Metabolism, Shandong Provincial Hospital, Jinan, China
- 3Department of Endocrinology, Shandong Provincial Hospital, Shandong University, Jinan, China
- 4Shandong Provincial Institute of Dermatology and Venereology, Shandong University, Jinan, China
Background: Diabetic kidney disease (DKD) becomes the leading cause of death for end-stage renal disease, whereas the potential mechanism is unclear and effective therapy is still rare. Our study was designed to investigate the cellular mechanism of Forsythiaside against DKD.
Materials and Methods: The targets of Forsythiaside and the DKD-related targets were obtained from databases. The overlapping targets in these two sets were regarded as potential targets for alleviation of DKD by Forsythiaside. The targets of diabetic podocytopathy and tubulopathy were also detected to clarify the mechanism of Forsythiaside ameliorating DKD from the cellular level.
Results: Our results explored that PRKCA and RHOA were regarded as key therapeutic targets of Forsythiaside with excellent binding affinity for treating DKD podocytopathy. Enrichment analysis suggested the underlying mechanism was mainly focused on the oxidative stress and mTOR signaling pathway. The alleviated effects of Forsythiaside on the reactive oxidative species accumulation and PRKCA and RHOA proteins upregulation in podocytes were also confirmed.
Conclusion: The present study elucidates that Forsythiaside exerts potential treatment against DKD which may act directly RHOA and PRKCA target by suppressing the oxidative stress pathway in podocytes. And Forsythiaside could be regarded as one of the candidate drugs dealing with DKD in future experimental or clinical researches.
1 Introduction
As the population who suffers from diabetes surges rapidly, the prevalence rate of diabetic kidney disease (DKD) significantly increases and occurs in approximately 40% of patients with type 2 diabetes mellitus (de Boer et al., 2011). DKD as one of the most deleterious complications of diabetes, has become the leading cause of end-stage renal disease (ESRD), which considerably does harm for the life quality of patients and lays serious economic burden in the society and families (Doshi and Friedman, 2017). A sufficient understanding the underlying mechanisms of renal injury under the condition of high glucose facilitates the advancement of prevention and treatment of DKD.
A variety of cells were involved in the pathogenesis of DKD and the clinical effects and mechanism of DKD occurrence and development have been preliminarily elucidated. Studies showed that glomerular changes were undoubtedly major characteristics of DKD (Gilbert, 2017). Traditional study revealed that glomerular podocyte injury played an important role in the process of DKD development (Chen et al., 2020). Furthermore, the role of tubular abnormality in promoting DKD progression has been gradually illustrated. High glucose can do harm to the kidney through various pathways, including renal glomerular (Ren et al., 2017), tubular (Zhao et al., 2021), vascular (Zhan et al., 2021) and renal interstitial injuries (Cheng et al., 2020). Intensive glucose control in DKD could slow down but not significantly prevent the progression of DKD (MacIsaac et al., 2017). Other conventional therapy which was restricted to renin-angiotensin aldosterone system inhibitor also showed limited benefit for delaying the disease progression (Coleman et al., 2020). Thus, conventional western therapies fail to provide favorable utility for DKD.
Traditional Chinese Medicine (TCM) shows remarkable therapeutic effect on chronic diseases, and has been widely supplemented to intervene DKD in China and other Asian countries (Lovisa et al., 2016; Chen D. Q. et al., 2018; Chen H. et al., 2018; Wang et al., 2018). Chinese botanical drug medicine as a treasure trove of TCM is rich in lots of natural compounds, which are effective in treating diseases for their achievement of the multi-target, and multi-mechanism system regulation. However, there is less research of natural compounds from the Chinese botanical drug medicine against renal injury in DKD (Zhang et al., 2019).
According to the theory of TCM, Forsythiaside, as a main constituent of Forsythia suspensa (Thunb.) Vahl, has long been thought to have benefits involving clearing heat and detoxifying, dispelling carbuncle and is commonly used in treating infections and pyrexia (Guo et al., 2007; Lu et al., 2010). Previous pharmacologic studies indicated that the Forsythiaside has multiple biological functions such as anti-inflammatory (Wang et al., 2016), antioxidant (Lu et al., 2010), and hepatoprotection (Zhang et al., 2018). Nowadays, Forsythiaside has a wide range of applications in various diseases, for instance, inflammation (Zeng et al., 2017), virus infection (Law et al., 2017), neurodegeneration (Chen L. et al., 2019), oxidative stress (Huang et al., 2015), and liver injury (Gong et al., 2021). Recently, study revealed that Forsythiaside A showed protective effects against adriamycin-induced nephropathy in rat by reducing proteinuria, serum creatinine and urea nitrogen levels and apoptotic cells (Lu et al., 2020). Moreover, Forsythiaside alleviated adriamycin-triggered renal injury and enhanced superoxide dismutase activities, inhibited malondialdehyde and lactate dehydrogenase productions (Lu et al., 2020). However, the influence and underlying mechanism of Forsythiaside on the alleviation of DKD is yet to be completely elucidated. Therefore, the nephro-protective activities of Forsythiaside in DKD were further investigated in this study.
With the continuous progress and rapid development of various high-throughput technologies, a large number of scientific researches have surged, which further contribute to an enormous and novel bioinformation network. Network analysis as an emerging bioinformation approach for drug discovery and development shows the marked advantages of convenience of acquiring information, low cost, and a flood of valuable resource. Network analysis emerges as the new mode of large-scale data processing develops and investigates the relationship between compounds and medicines, diseases and targets by calculating the synergistic probability of multi-compound, multi-target and multi-pathway (Hopkins, 2007). And network analysis is useful to explain the mechanism of diseases from the systemic view and guides the discovery of new medicines for elucidating the potential interaction of the body and medicines. It has gradually become a new trend to apply the methods of network analysis to clarify the possibility and mechanism of Chinese botanical drug medicine and compound in treating disease.
Nowadays, there have been major advances in the understanding and treatment of DKD, whereas few researches about natural compounds of Chinese botanical drug medicine relieving the development of DKD from evidence-based single and specific cellular facet were available. Therefore, the objective of this study was to reveal the effects of Forsythiaside on the DKD alleviation and its mechanism in the treatment of DKD from the specific cellular level. The target proteins of Forsythiaside and DKD, diabetic podocytopathy and diabetic tubulopathy were retrieved and established respectively by public database. Then the enrichment analyses were further carried out. The molecular docking simulation based on “lock-key principle” and energy matching was performed to calculate the binding affinity between macromolecule and ligand for virtually prediction of drug targets. Finally, the potential mechanism of Forsythiaside alleviating renal damage was preliminarily confirmed and Forsythiaside-targets-pathways network was constructed and visualized by Cytoscape 3.8.2 software (Pang et al., 2018). Our research could provide preliminarily theoretical basis for the curative effects of natural compounds of TCM on the treatments of DKD and broaden the horizon for clinical integration of traditional Chinese and Western medicine on DKD treatment in future.
2 Materials and methods
2.1 Analysis of predicted targets of forsythiaside
The drug-like properties of Forsythiaside are summarized in Table 1. The potential intervention targets of Forsythiaside were extracted from the ChEMBL database (Gaulton et al., 2017) (https://www.ebi.ac.uk/chembl/g/), the PubChem database (Kim et al., 2019) (http://pubchem.ncbi.nlm.nih.gov), the Swiss Target Prediction database (Gfeller et al., 2014) (http://www.swisstargetprediction.ch/), the Binding database (Chen et al., 2001) (https://www.bindingdb.org/bind/index.jsp), and the PharmMapper database (Wang et al., 2017) (http://lilab-ecust.cn/pharmmapper/index.html). All potential targets were summarized and their names standardized using the universal Protein Resource site (RRID: SCR_002380) (http://www.uniprot.org/).
2.2 Prediction of diabetic kidney disease-related target proteins
The target proteins involved in the pathogenesis of DKD were predicted and acquired from a variety of databases. The terms “Diabetic kidney disease” and “Diabetic nephropathy” were used as keywords across multiple databases. We used Therapeutic Target (Wang et al., 2020) (http://bidd.nus.edu.sg/group/cjttd/), GeneCards (Rebhan et al., 1998) (https://www.genecards.org/), DisGeNET (RRID: SCR_006178) (https://www.disgenet.org/) (Piñero et al., 2020), DrugBank (Wishart et al., 2008) (https://go.drugbank.com/) and Online Mendelian Inheritance in Man (OMIM, RRID: SCR_006437) (https://omim.org/) (Boyadjiev and Jabs, 2000) databases to screen and identify potential targets of DKD.
2.3 Identification of overlapping targets of forsythiaside and diabetic kidney disease
2.3.1 Acquiring candidate target molecules by target mapping
To clarify the candidate target molecules via which Forsythiaside could intervene in DKD, the predicted targets of Forsythiaside were mapped onto the DKD-related proteins and overlapping targets were extracted as key genes. Overlapping targets were visualized using a Venn diagram.
2.3.2 Construction of a protein-protein interaction network
Once mapped, overlapping target molecules between Forsythiaside and DKD were inserted into the STRING website (Szklarczyk et al., 2017) (https://string-db.org/), a useful tool for the discovery of functional protein interactions and further generation of a PPI network.
2.3.3 Network construction based on overlapping targets of forsythiaside and diabetic kidney disease
Cytoscape software (CluePedia Cytoscape plugin, RRID: SCR_015784) (https://cytoscape.org/), an open source software platform for visualizing intricate networks and integrating networks with attribute data (Shannon et al., 2003), was applied to establish the networks of Forsythiaside and DKD’s respective targets, the PPI network of overlapping targets between Forsythiaside and DKD, and overlapping target-pathway networks. The properties of each PPI network were calculated and analyzed using the Cytoscape plug-in Network Analyzer. Higher parameter values indicate that a node is more significant.
2.3.4 KEGG enrichment analysis based on overlapping targets between forsythiaside and diabetic kidney disease
To further clarify the potential role of Forsythiaside in alleviating DKD, overlapping targets were inserted into the Database for Annotation, Visualization and Integrated Discovery site (Huang, et al., 2009) (version 6.8; https://david.ncifcrf.gov/) to carry out KEGG enrichment analysis. Bubble diagrams of KEGG enrichment analysis were drawn using R Project for Statistical Computing (RRID: SCR_001905) (https://www.r-project.org/).
2.4 Analysis of diabetic podocytopathy and tubulopathy-related target molecules
In order to specify the cellular mechanism of Forsythiaside in the treatment of DKD, target molecules which were related to diabetic podocytopathy and diabetic tubulopathy were predicted and acquired from multiple databases. The terms “Diabetic podocytopathy” or “Diabetic tubulopathy” were inserted into the GeneCards and OMIM databases as keywords.
2.5 Establishment of a protein-protein interaction network, and network analysis based on overlapping targets of forsythiaside and diabetic podocytopathy or tubulopathy
Overlapping targets between Forsythiaside and diabetic podocytopathy were acquired by drawing a Venn diagram, and candidate diabetic tubulopathy target genes acted on by Forsythiaside were obtained by using similar methods.
Upon target mapping, overlapping target molecules between Forsythiaside and diabetic podocytopathy or tubulopathy were used to generate a PPI network using the STRING database (RRID: SCR_005223) (https://cn.string-db.org/) following the methods described in Section 2.3.2.
Furthermore, PPI networks were developed for the following: diabetic podocytopathy and targets; diabetic tubulopathy and targets; overlapping targets between Forsythiaside and diabetic podocytopathy; overlapping targets between Forsythiaside and diabetic tubulopathy; and overlapping target-pathway network were constructed using Cytoscape software (RRID: SCR_003032) (https://cytoscape.org/). The values of degree were calculated using the Cytoscape plug-in Network Analyzer, and the nodes were ranked by these values.
2.6 GO pathway and KEGG enrichment based on overlapping proteins between forsythiaside and diabetic podocytopathy or tubulopathy
In order to elucidate the pathways through which Forsythiaside may intervene in diabetic podocytopathy or tubulopathy, the overlapping targets of Forsythiaside and diabetic podocytopathy or tubulopathy were inserted into the Database for Annotation, Visualization and Integrated Discovery (version 6.8; https://david.ncifcrf.gov/) and the Metascape database (RRID: SCR_016620) (http://metascape.org/) (Zhou et al., 2019). The top 10 pathways were displayed.
Upon three enrichment analyses, the key pathways through which Forsythiaside alleviates DKD, diabetic podocytopathy, and tubulopathy, respectively were identified. The same pathways were extracted as key ones that Forsythiaside might alter to treat DKD, by ameliorating targets in diabetic podocytopathy or tubulopathy. The overlapping Forsythiaside and diabetic podocytopathy or tubulopathy target proteins that were identified as critical proteins through which Forsythiaside may intervene in order to treat DKD.
2.7 Molecular docking simulation
2.7.1 Acquirement and preparation of protein receptors
The crystal structures of critical proteins were acquired from RCSB Protein Data Bank database (Goodsell et al., 2020) (http://www.pdb.org/), and PDB codes of the protein receptors were summarized in Table 2. PyMOL software (version 2.5.0) was used to edit the protein complexes by removing the original ligand and water molecules from them, and the modified protein was saved in pdb format (Lu, 2020). Autodock software (version 4.2.6) (Morris et al., 2009), a suite of automated docking tools, and Autodock Tools (version 1.5.6) (http://mgltools.scripps.edu/documentation/links/autodock) were used to prepare the appropriate protein receptor by adding hydrogen and computing charges. The protein macromolecule was saved in pdbqt format.
2.7.2 Ligand preparation
Before docking, the structure of Forsythiaside was acquired using PubChem and the Traditional Chinese Medicine System Pharmacology Database (http://lsp.nwu.edu.cn/tcmsp.php) (Ru et al., 2014), which was saved in Mol2 format and further analyzed. The conformation of Forsythiaside (Pubchem CID: 5281773) was displayed in Supplementary Figure S1. Then, Autodock software (version 4.2.6) was applied to transfer the structure in Mol2 format to pdbqt format to allow energy minimization for the docking.
2.7.3 Molecular docking
Upon the preparation of ligand and receptor, Autodock software (version 4.2.6) was used to evaluate and confirm the binding affinity between ligand and receptor and further construct the binding model. Binding models were visualized using PyMOL software (RRID: SCR_000305) (https://pymol.org/2/), and Discovery Studio Visualizer software (version 21.1.0.20298) was applied to transform the 3D structure of the complex to 2D.
2.8 Reagent, cell culture and treatment
Forsythiaside A powder (FA, the main component of Forsythiaside, purity: 99.43%) was purchased from MedChemExpress (cat. HY-N0028; Monmouth Junction, NJ, United States), and dissolved in DMSO at various concentrations (2.5, 5, and 10 μg/mL). Conditionally immortalized human podocytes were kindly supplied by Dr. Yi Fan (Department of Pharmacology, Shandong University School of Medicine, Jinan, China) and cultured in RPMI-1640 medium (cat no. 11875093; Gibco, Grand Island, NY, United States) supplemented with 10% fetal bovine serum (cat no. 10099141; Gibco), 100 units/mL penicillin and 100 mg/mL streptomycin (1% P/S concentration, cat no. 10378016; Gibco) until cells became confluent.
Podocytes were incubated in various media containing 5.5 mM D-glucose (NG), 5.5 mM D-glucose +34.5 mM D-mannitol (MG) (Chen Y. et al., 2019), 40 mM D-glucose (HG), or 40 mM D-glucose + Forsythiaside at various concentrations (HG + FA 2.5, 5, or 10 μg/mL). D-glucose and D-mannitol were purchased from Sigma-Aldrich Canada Co. (cat nos. 50-99-7 and #69-65-8; Oakville, ON, Canada). After treatment, podocytes were harvested for protein extraction.
2.9 Cell viability assessment by Cell Counting Kit-8 assay
Human podocytes were seeded at a density 3×104 cells/cm2 in RPMI 1640 medium on 96-well plates. After cell adherence, RPMI-1640 medium was replaced with other media, including NG, MG, HG, or HG + FA at concentrations of 2.5, 5, or 10 μg/mL. Cells were subsequently incubated for 48 h. Cell viability was detected using the CCK-8 assay kit (cat no. CK04-11; Dojindo Molecular Technologies, Inc., Kumamoto, Japan) according to the manufacturer’s instructions. The absorbance of each sample was measured at a wavelength of 450 nm.
2.10 Western blotting
Proteins were isolated from cultured podocytes after various treatments. Protein levels were detected using a BCA Protein Assay Kit (cat no. P0012; Beyotime, Shanghai, China). Proteins were subsequently separated by SDS-PAGE electrophoresis. Antibodies against protein kinase C alpha (PRKCA) (cat no. 21991-1-AP; Proteintech Group Inc., Wuhan, China, 1:1,000) and RHOA (cat no. 10749-1-AP; Proteintech, 1:1,000) were used for western blotting. An antibody against β-actin (cat no. 60008-1-Ig; Proteintech, 1:5,000) was used as a reference.
2.11 Reactive oxidative species measurement
ROS were measured in podocytes using a ROS detection kit (cat no. S0033S; Beyotime, Shanghai, China) according to the manufacturer’s instructions. Briefly, podocytes were seeded on sterile glass coverslips in a 6-well plate and incubated with different media for 48 h. Thereafter, slides were incubated with a 10 μM DCFH-DA working solution for 30 min at 37°C. After washing three times with warm PBS, ROS intensity was detected under a fluorescent microscope (Olympus FSX100).
2.12 Statistical analysis
Technical and biological replicates were performed three times. Data were processed using SPSS software (RRID: SCR_002865) (https://www.ibm.com/products/spss-statistics) and presented as means ± standard error of the mean. Student’s t-test was used to evaluate the significance of difference between two groups. One-way ANOVA and two-way ANOVA tests were used for multiple statistic comparisons. p < .05 was considered as statistically significant.
3 Results
3.1 Forsythiaside-target network and diabetic kidney disease-target network construction
A flowchart summarizing the results of the network analysis study is shown in Figure 1. We used the ChEMBL, PubChem, Swiss Target Prediction, Binding, and PharmMapper databases to screen and acquire the putative targets of Forsythiaside. 534 predicted targets were identified as the potential target proteins of Forsythiaside (Figure 2A). The target proteins were standardized using universal Protein Resource site, and detailed information is given in Supplementary Table S1.
FIGURE 1. Flowchart of exploring the mechanisms of Forsythiaside in DKD from the specific cellular level.
FIGURE 2. Compound-target network and disease-target network. (A) Forsythiaside-target network. Green diamond node represents Forsythiaside and red circular nodes represent corresponding targets. (B) DKD-target network. Yellow circular node represents DKD and green circular nodes represent corresponding targets. FT: Forsythiaside; DKD: diabetic kidney disease.
By using the Therapeutic Target, GeneCards, DisGeNET, DrugBank, and OMIM databases, a total of 1,554 target proteins were identified as potential targets of DKD (Figure 2B). Information on target proteins is given in Supplementary Table S2.
3.2 Protein-protein interaction network of Forsythiaside-diabetic kidney disease targets
To acquire the targets through which Forsythiaside acted on DKD, the 534 putative targets of Forsythiaside were mapped to the 1554 DKD-related targets to obtain overlapping targets of Forsythiaside and DKD. The overlapping targets were visualized using a Venn diagram shown in Figure 3A. Consequently, 155 targets were identified as candidate targets through which Forsythiaside might alleviate DKD. Cytoscape software (version 3.8.2) was applied to construct the PPI network, which could evaluate the action of the target proteins in disease, and reveal their potential interactions. The PPI network of overlapping proteins between Forsythiaside and DKD was established, and 155 targets were sorted in descending order by degree. As shown in Figure 3B, targets were presented as circles and ordered by the degree. Hub targets were presented in the innermost center circle, and the core targets were AKT1, ALB, TNF, MAPK3, MAPK1, MMP9, and EGFR. The interactive relationship was analyzed using the STRING tool; results are displayed in Figure 3C.
FIGURE 3. Venn diagram, PPI network and KEGG enrichment of Forsythiaside-DKD targets. (A) Venn diagram of overlapping targets of Forsythiaside and DKD. (B) PPI network of Forsythiaside-DKD targets. A total of 155 overlapping targets were presented as circle and ordered by the degree level. The hub targets were presented in the innermost center circle, and the core targets were AKT1, ALB, TNF, MAPK3, MAPK1, MMP9, and EGFR. (C) The interactive relationship was analyzed and displayed as PPI network by STRING website. (D) KEGG enrichment analysis of potential targets that Forsythiaside acted for alleviating DKD. DKD: diabetic kidney disease.
3.3 KEGG pathway enrichment analysis
To systematically elucidate the mechanism of Forsythiaside treatment, KEGG enrichment analysis of 155 overlapping genes between Forsythiaside and DKD was performed. The top 10 pathways were selected and drawn according to counts of hub genes and p values (Figure 3D). Several pathways were selected as the key action mechanism underlying DKD treatment with Forsythiaside, including Type II diabetes mellitus and mTOR signaling pathways.
3.4 Diabetic podocytopathy and tubulopathy target network
The GeneCards and OMIM databases were used to acquire the targets of diabetic podocytopathy. A total of 120 proteins were documented as potential targets of diabetic podocytopathy (Figure 4A). Detailed information on these proteins was recorded in Supplementary Table S3. The same analysis was conducted, and 633 predicted target proteins of diabetic tubulopathy were identified from the databases (Figure 4B). Relative information of proteins is presented in Supplementary Table S4.
FIGURE 4. Disease-target network. (A) Diabetic podocytopathy-target network. Yellow triangle node represents diabetic podocytopathy and blue circular nodes represent corresponding targets. A total of 120 proteins were documented as the potential targets of diabetic podocytopathy. (B) Diabetic tubulopathy-target network. Yellow triangle node represents diabetic tubulopathy and blue circular nodes represent corresponding targets. Total 633 predicted target proteins of diabetic tubulopathy were identified.
3.5 Protein-protein interaction networks of forsythiaside-diabetic podocytopathy and tubulopathy targets
A Venn diagram was drawn to identify the overlapping targets between Forsythiaside and diabetic podocytopathy (Figure 5A). Nine overlapping targets were identified and arranged based on the degree level, including PPARG, INSR, PRKCA, SOD2, GCK, ACE, RHOA, HDAC1, and CCR5 (Figure 5B). Additionally, the interaction relationship between nine overlapping targets were further investigated. And the PPI network which reflected the interaction between these targets was displayed in Figure 5C.
FIGURE 5. Venn diagram and PPI network of Forsythiaside-disease targets. (A) Venn diagram of overlapping targets of Forsythiaside and diabetic podocytopathy. (B) PPI network of Forsythiaside-diabetic podocytopathy targets. Nine overlapping targets were presented as circle and ordered by the degree level. (C) The PPI network which reflected the interaction between nine overlapping targets was displayed by STRING website. (D) Venn diagram of overlapping targets of Forsythiaside and diabetic tubulopathy. (E) PPI network of Forsythiaside-diabetic tubulopathy targets. Thirty-six overlapping targets were presented as circle and ordered by the degree level. Five genes were of significance in the network, including ALB, PPARG, IGF1, F2 and ACE. (F) The PPI network analysis of 36 overlapping targets between Forsythiaside and diabetic tubulopathy was performed by STRING website.
Similarly, the overlapping genes between Forsythiaside and diabetic tubulopathy were selected. A total of 36 genes were identified, as shown in a Venn diagram (Figure 5D). A circle network was drawn and genes were ordered by their degree levels. Five genes were of significance in this network, including ALB, PPARG, IGF1, F2, and ACE (Figure 5E). PPI network analysis was also performed; the result is shown in Figure 5F.
3.6 GO and KEGG pathway enrichment analysis of overlapping genes of forsythiaside-diabetic podocytopathy and forsythiaside-diabetic tubulopathy
For clarifying the potential mechanism of Forsythiaside on diabetic podocytopathy and tubulopathy systematically, enrichment analyses were conducted with the 9 and 36 discovered overlapping genes, respectively. For diabetic podocytopathy, the top GO items were the response to oxidative stress, cellular components of side of membrane, and molecular function of kinase activity (Figure 6A). The top KEGG pathways were the mTOR signaling pathway, Type II diabetes mellitus, and Rap1 and Ras signaling pathway (Figure 6B). For diabetic tubulopathy, the overlapping targets were primarily enriched in the positive regulation of ion transport, cellular components of side of membrane, and molecular function of lipid binding (Figure 6C). Additionally, the top KEGG enrichments were associated with the longevity-regulating and AMPK signaling pathways (Figure 6D).
FIGURE 6. GO and KEGG enrichment of Forsythiaside-disease targets. (A) GO and (B) KEGG enrichment analysis of potential targets that Forsythiaside acting on diabetic podocytopathy. (C) GO and (D) KEGG enrichment analysis of potential targets that Forsythiaside acting on diabetic tubulopathy.
3.7 Molecular docking analysis
After conducting enrichment analysis, we found that the enrichment pathways of the overlapping targets of Forsythiaside and DKD were consistent with those of Forsythiaside and diabetic podocytopathy, suggesting that Forsythiaside alleviated DKD mainly through its action on podocytes. Thus, molecular docking simulation was performed to verify the possibility of Forsythiaside binding to the targets in diabetic podocytopathy. The potential interaction between 9 overlapping genes and Forsythiaside was explored with molecular docking verification. And the complexes of retrieved targets and Forsythiaside were filtered by the docking affinity values which were reported by AutoDock. The greater absolute value of minus value of the docking affinity showed the stronger binding affinity of Forsythiaside and the targets. Total nine docking complexes were analyzed and the models of binding complexes were displayed. The number of conformations analyzed for docking was 50. And the corresponding free energy of binding and the putative affinity were summarized in Table 2. Docking data suggested strong binding affinity of RHOA- Forsythiaside docking (−1.14 kcal/mol, Figures 7A–C) and PRKCA- Forsythiaside docking (−.15 kcal/mol, Figures 7D–F), while other 7 targets showed weak binding affinities with Forsythiaside. To be specific, for RHOA- Forsythiaside docking complex, small molecule ligand Forsythiaside could interact with the residue of RHOA protein via a wide range of interaction forces. As shown in Figure 7A, Forsythiaside could insert into the interface pocket of RHOA protein. Docking complex displayed that five hydrogen bonds were formed between Forsythiaside and RHOA’s residues including PHE-39, GLU-40, and GLY-62 residue of RHOA. And the distances of hydrogen bonds and Forsythiaside were 2.8, 2.4, and 2.7, 1.9 and 1.9 Å, respectively (Figure 7B). Other interaction forces were also found in the Forsythiaside -RHOA docking complex. ALA-61 and VAL-38 residue contributed to form pi-alkyl interactions with Forsythiaside. GLN-63, GLU-64, GLY-14, THR-37, and THR-19 residue interacted with Forsythiaside via van der Waals forces. The carbon hydrogen bond and unfavorable donor-donor forces were also found in the interaction of Forsythiaside with TYR-66 and PHE-39 residue, respectively (Figure 7C).
FIGURE 7. Molecular docking data of RHOA target and PRKCA target and compound-target-pathway network of Forsythiaside and diabetic podocytopathy. (A) Schematics (3D) showed that molecular model of Forsythiaside fitting in the binding pocket of the RHOA protein. (B) Schematics (3D) represented the interactions between Forsythiaside and surrounding residues of RHOA. The yellow dashed lines showed hydrogen bonds and the interaction distances were indicated. (C) Schematics (2D) revealed the interaction forces between Forsythiaside and RHOA. (D) Schematics (3D) represented that molecular model of Forsythiaside fitting in the binding pocket of the PRKCA protein. (E) Schematics (3D) showed the interactions between Forsythiaside and surrounding residues of PRKCA. The yellow dashed lines showed hydrogen bonds and the interaction distances were indicated. (F) Schematics (2D) revealed the interaction forces between Forsythiaside and PRKCA. (G) The compound-target-pathway network was established to explore the potential mechanisms of Forsythiaside treating diabetic podocytopathy.
In the PRKCA-Forsythiaside docking complex, Forsythiaside could potentially fit into the interface pocket of the PRKCA protein (Figure 7D). The complex showed that three hydrogen bonds were formed between Forsythiaside and GLN-513 and ALA-511 of PRKCA. Interactions between the three hydrogen bonds are shown in Figure 7E with a distance of 2.5, 2.2, and 1.8 Å, respectively. Moreover, the GLU-543 residue contributed to form pi-anion, and the ALA-511 residue contributed to form pi-alkyl, interactions with Forsythiaside. The other residues, such as TYR-512, ILE-510 and MET-551, PHE-547, and ASP-544 interacted with Forsythiaside via van der Waals forces (Figure 7F). Consequently, Forsythiaside stably bound to PRKCA and RHOA proteins through multiple interactive forces.
3.8 Compound-target-pathway network
As shown in Figure 7G compound-target-pathway network was established to explore the potential mechanism of Forsythiaside treatment of diabetic podocytopathy. Two targets, PRKCA and RHOA proteins, and four signaling pathways, including the mTOR signaling pathway, Type II diabetes mellitus, and Ras and Rap1 signaling pathway, were predicted and considered as key nodes and pathways through which Forsythiaside ameliorated DKD podocytopathy.
3.9 Forsythiaside relieved impaired viability of podocytes induced by high glucose treatment
To investigate the direct effect and drug toxicity of Forsythiaside on podocytes, we performed a CCK-8 assay. Firstly, cell viability of podocytes was examined after high glucose stimulation. As shown in Figure 8A, cell viability declined in the HG podocyte group compared to the NG podocyte group (p < .05). Cell viability was further assessed in the HG + FA group. Results showed that the viability of podocytes in this group were not significantly different when compared with podocytes in the NG group (p > .05, Figure 8A), indicating that the concentrations of Forsythiaside used in the experiment were non-toxic. Additionally, we found that there was a significant increase in viability in media with 2.5 and 5 μg/mL Forsythiaside compared with the HG group (p < .05). However, treatment with 10 μg/mL Forsythiaside did not show significant increase of cell viability when compared with HG group (p > .05).
FIGURE 8. Assessment of podocyte viability by Cell Counting Kit-8 assay and reactive oxidative species measurement by fluorescence. (A) Podocyte viability was examined by Cell Counting Kit-8 assay after various treatments for 48 h. (B) Reactive oxidative species which were labeled with green fluorescence were detected after podocytes were exposed to normal glucose (NG), isosmotic mannitol (MG), or high glucose (HG) for 48 h with or without Forsythiaside (FA 2.5 μg/mL, FA 5 μg/mL, and FA 10 μg/mL) treatment. Magnification: ×200. Scale bar = 50 µm. Data are mean ± S.E.M. *p < .05. FA: Forsythiaside.
3.10 Forsythiaside significantly inhibited ROS production in podocytes under high glucose conditions
To clarify a potential mechanism for the renal protective effect of Forsythiaside, potential changes in ROS were measured with a reagent solution, DCFH-DA, which fluoresced upon oxidation by ROS. As shown in Figure 8B, increased ROS production was observed in podocytes in the HG group compared to those in the NG and MG groups, manifesting as green fluorescence. However, Forsythiaside significantly reduced the HG-induced increase in ROS. These results suggest a potential nephroprotective effect of Forsythiaside related to the inhibition of oxidative stress inhibition in podocytes, consistent with our network analysis; findings are described in Section 3.6.
3.11 Forsythiaside suppressed the expression of PRKCA and RHOA in podocytes treated with high glucose
As shown in Figures 9A,B marked increase in expression of PRKCA protein in podocytes was detected in the HG group when compared to the NG and MG groups, as detected by western blotting (p < .05), while treatment with Forsythiaside in HG medium visibly suppressed the expression of PRKCA (p < .05). Additionally, western blotting results showed that the expression of RHOA was significantly upregulated in podocytes in the HG group compared with those in the NG and MG groups (p < .05), which was in accordance with the results of PRKCA protein. It is worth noting that 2.5 μg/mL Forsythiaside treatment markedly inhibited the expression of RHOA protein compared with the HG group (p < .05, Figures 9A,C). However, there was no significant decrease of RHOA protein expression in 5 μg/mL and 10 μg/mL Forsythiaside groups (Figure 9C). The results of western blotting confirmed the network analysis and molecular docking findings described in Sections 3.5, 3.7. Taken together, our results demonstrate the potential protective mechanism of Forsythiaside under high glucose conditions may be attributed to the suppression of a PRKCA and RHOA-related pathway.
FIGURE 9. Examination of PRKCA and RHOA protein expression levels after Forsythiaside treatment. (A) The protein levels of PRKCA and RHOA in podocytes were examined by Western blotting. (B) Quantification of PRKCA protein expression in podocytes after various treatments. (C) Quantification of RHOA protein expression in podocytes after different treatments. Values are expressed as mean ± S.E.M. *p < .05. NG: normal glucose; MG: isosmotic mannitol; HG: high glucose; FA: Forsythiaside.
4 Discussion
Of the common chronic complications of diabetes, DKD in one of the heaviest burdens for patients (Thomas et al., 2015). Conventional therapies such as intensive control of blood glucose, pressure, and lipids have not effectively inhibited DKD progression (Ahlqvist et al., 2015). Additional factors including oxidative stress, inflammation, and metabolic alteration have also been verified to influence the occurrence and development of DKD (Lehtonen, 2020). Current advances in treatment for DKD show few benefits in patient prognosis and survival. Thus, key pathways that may be involved in the progression of DKD must be investigated in order to identify effective treatments for DKD.
Bioactive constituents extracted from a variety of medicinal plants are considered more effective than synthetic medications for the management of multiple chronic diseases (Gorji et al., 2018; Mao et al., 2019). Therefore, natural products are regarded to have potential applications in the treatment of DKD. Forsythiaside has been shown to be an effective compound for the prevention and treatment of chronic disease (Cheng et al., 2015). In recent years, studies on the pharmacological effects of Forsythiaside have gradually increased. These show that Forsythiaside has anti-oxidation (Huang et al., 2015), anti-inflammation (Chen L. et al., 2018), and immunoregulatory (Cheng et al., 2014) properties. Studies have shown that one of the core mechanisms of DKD development is the accumulation of ROS, inflammatory factors, and immune complexes in the kidney, which further activate cytotoxic signaling in the kidney, leading to renal dysfunction (Huang et al., 2020). Forsythiaside may protect against DKD by intervening these mechanisms. However, to our knowledge there are few studies on the material basis and mechanism of the prevention and treatment of DKD by Forsythiaside.
In this study, we found that Forsythiaside was effective in preventing and treating DKD by performing network analysis. The pathways of Forsythiaside in the treatment of DKD were mTOR signal pathway and Type II diabetes mellitus. Our findings revealed that Forsythiaside could not only intervene various targets to regulate the same pathway, but also show the ability to act on the same target to control the multiple pathways for management of DKD, preliminarily indicating that the natural compound exerts multi-target and multi-pathway regulating effect on the improvement of the chronic and complicated disease.
A great number of cell types are involved in the pathogenesis of DKD. Deciphering the molecular mechanisms that regulate the occurrence and development of DKD, and exploring effective therapeutic measures for delaying renal impairment are essential. As the two inherent components of the kidney, the renal tubule and glomerular podocytes are the most frequently studied biological units in renal dysfunction. The proximal renal tubule is sensitive to various pathogenic factors, including metabolic and hemodynamic stimuli which induce the apoptosis of tubular cells and further lead to the progression of tubulo-interstitial lesions in DKD (Habib, 2013). Glomerular podocytes are the primary mediators of the integrity and normal function of the glomerular filtration barrier. Damage to podocytes causes detachment and dysfunction, which further promotes DKD progression (Pagtalunan et al., 1997; Steffes et al., 2001).
For clarifying the exact mechanism from the cellular aspect, the genes and pathways that Forsythiaside intervened were mapped to the targets and pathways of diabetic podocytopathy and tubulopathy by construction of network. For exploring how Forsythiaside treated DKD from the cellular level, we analyzed the overlapping pathways and related genes that Forsythiaside played an important role in the pathogenesis of DKD and diabetic podocytopathy or tubulopathy simultaneously. And the overlapping pathways, including mTOR signaling pathway and Type II diabetes mellitus were identified as the key pathways that Forsythiaside ameliorated DKD by improving podocytopathy. Few overlapping pathways were observed between DKD and diabetic tubulopathy in which Forsythiaside worked, indicating the mechanism the Forsythiaside alleviated DKD was mainly involved in the improvement of podocytopathy other than tubulopathy.
Further molecular docking simulation was performed to elucidate the specific target protein in podocytopathy that Forsythiaside worked to intervene for the amelioration of DKD, and we found that Forsythiaside could play a key role in the treatment of DKD podocytopathy by targeting RHOA and PRKCA protein.
Ras homolog family member A (RHOA) belongs to the family of Small GTPase. The primary function of RHOA is mainly related with cytoskeleton organization when binding to a variety of effector proteins, and then is responsible for control of cellular responses (Kloc et al., 2019). Researches revealed that RHOA exerted the detrimental effects on the progression of DKD. The suppression of the RHOA/ROCK pathway also possessed a protective capability for glomerular inflammation and fibrosis in diabetic rats (Hirose et al., 2010). Another research indicated that catalpol rescued disrupted cytoskeleton in podocytes under high glucose, which was mainly attributable to the inhibition of RHOA activities (Rao et al., 2017). Consistent with previous studies, our research revealed that Forsythiaside inhibited the expression of RHOA protein in podocytes. Therefore, we inferred that RHOA may be involved in the dysfunction of podocytes in DKD pathogenesis and Forsythiaside could target and intervene RHOA for management of diabetic podocytopathy.
Protein kinase C alpha (PRKCA), another target of Forsythiaside we acquired from network analysis and molecular docking in this study, is a serine/threonine-protein kinase which is involved in the regulation of a variety of cellular processes by activating signaling cascade involving mTOR pathway (Bhalla et al., 2017). Previous studies revealed that hyperglycemia induced the activation of PKC system which played a vital role in DKD glucotoxicity (Craven and DeRubertis, 1989; Lee et al., 1989) and the abnormal modulation of renal function (Chen et al., 2004). Activation of PRKCA has been implicated as a central mediator in the pathogenesis of DKD and the development of albuminuria (Kang et al., 1999). The inhibition of PRKCA exerted the beneficial effect in DKD (Parameswaran et al., 2000; Koya, 2014). Combined with our findings, Forsythiaside ameliorated the detrimental effects of podocyte damage under HG condition possibly by preventing PRKCA upregulation.
The targets that Forsythiaside bound and intervened were primarily identified, whereas the underlying mechanism is still unclear. Further Go enrichment analysis was performed and the results showed that Forsythiaside could relieve oxidative stress and improve kinase activity in DKD podocytes. Dysfunction of podocytes has been verified as a critical driving factor for DKD occurrence and development (Bid et al., 2013; Reiser and Sever, 2013). Podocyte loss which induced glomerular damage and glomerulosclerosis may result from excess generation of ROS which further caused podocyte apoptosis (Susztak et al., 2006). Thus, prevention of redox reactions is probably beneficial for podocyte survival in DKD. Forsythiaside as the active component of could inhibit the increased ROS induced by hydrogen peroxide (Huang et al., 2015). RHOA was reported to play a vital role of induing excess oxidative stress (Bhalla et al., 2017) and Sanggenon C alleviated cerebral injury by suppressing oxidative stress via controlling RHOA signaling pathway (Shimokawa et al., 2016). Additionally, the role for PRKCA which controlled oxidative phosphorylation was identified in nephrotoxic nephritis mice (Kvirkvelia et al., 2018), suggesting that PRKCA was essential in the regulation of oxidative stress. Combined with the findings which Forsythiaside inhibited the expressions of RHOA and PRKCA proteins, and markedly alleviated the ROS accumulation, we supposed that Forsythiaside could act on the renal podocytes and suppress oxidative stress in diabetic podocytopathy via intervening RHOA and PRKCA.
Furthermore, our study clarified that Forsythiaside alleviated the detrimental effect of hyperglycemia on podocytes in diabetic podocytopathy possibly by inhibiting the mTOR signaling pathway. Research revealed that Everolimus, an inhibitor of mTOR remodeled aberrant podocyte behavior by substantial increase of RHOA activity and enhanced actin stress fibers (Jeruschke et al., 2013), suggesting that RHOA was involved in the regulation of mTOR signaling pathway, and the potential roles in DKD pathogenesis was preliminarily elucidated in our network analysis study. And PRKCA could also modulate a wide range of cellular processes by activating mTOR signaling pathways. These studies suggested that the pretreatment of Forsythiaside could efficiently act on the kidney to keep podocytes functioning normally through targeting RHOA and PRKCA and could be developed as a potential natural antioxidant and therapy for DKD podocyte injury by controlling mTOR signaling pathways.
In conclusion, by using network analysis, we found that Forsythiaside intervened DKD by mainly targeting diabetic podocytopathy. Furthermore, our study elucidated that Forsythiaside exerted potential treatment against DKD which may act directly RHOA and PRKCA target by suppressing the oxidative stress pathway in podocytes, which preliminarily reflected the characteristics of multi-target and multi-action pathway of Forsythiaside in treating DKD from the point of view of cellular precision therapy. Our findings could offer a reference for investigation of therapeutic benefits of Forsythiaside in DKD from the cellular facet. However, systematic studies of Forsythiaside intervention on DKD are still insufficient. Furthermore, limitations existed in the study. The lack of comprehensiveness in current network information, and the delayed update of database may become the inevitable limitations. The explanation of data should be more careful. Therefore, it is of great clinical significance to study the exact mechanism of Forsythiaside against DKD podocytopathy and to further develop related pharmaceutical preparations in future, and future research study should focus on the therapeutic interventions of Forsythiaside for DKD patients in clinic.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Materials and methods section.
Author contributions
All authors listed have made substantial contributions to the study. CX conceived and designed the study, HM searched the databases and conducted the in vitro experiments, XC analyzed the data and wrote the paper, and HZ reviewed and corrected the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China Grants (No. 82200916), Shandong Provincial Natural Science Foundation of China Grants (No. ZR2022QH068), Shandong Provincial Medicine and Health Science and Technology Development Program (No. 202103060766), and Shandong Provincial Hospital Research incubation Fund (No. 2021FY059).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2022.1096536/full#supplementary-material
SUPPLEMENTARY FIGURE S1 | The conformation of Forsythiaside which was acquired from PubChem database (Pubchem CID: 5281773).
References
Ahlqvist, E., van Zuydam, N. R., Groop, L. C., and McCarthy, M. I. (2015). The genetics of diabetic complications. Nat. Rev. Nephrol. 11, 277–287. doi:10.1038/nrneph.2015.37
Bhalla, M., Law, D., Dowd, G. C., and Ireton, K. (2017). Host serine/threonine kinases mTOR and protein kinase C-α promote InlB-mediated entry of Listeria monocytogenes. Infect. Immun. 85, 000877–e117. doi:10.1128/IAI.00087-17
Bid, H. K., Roberts, R. D., Manchanda, P. K., and Houghton, P. J. (2013). RAC1: An emerging therapeutic option for targeting cancer angiogenesis and metastasis. Mol. Cancer Ther. 12, 1925–1934. doi:10.1158/1535-7163.MCT-13-0164
Boyadjiev, S. A., and Jabs, E. W. (2000). Online Mendelian Inheritance in Man (OMIM) as a knowledgebase for human developmental disorders. Clin. Genet. 57, 253–266. doi:10.1034/j.1399-0004.2000.570403.x
Chen, D. Q., Feng, Y. L., Cao, G., and Zhao, Y. Y. (2018). Natural products as a source for antifibrosis therapy. Trends Pharmacol. Sci. 39, 937–952. doi:10.1016/j.tips.2018.09.002
Chen, H., Yang, T., Wang, M. C., Chen, D. Q., Yang, Y., and Zhao, Y. Y. (2018). Novel RAS inhibitor 25-O-methylalisol F attenuates epithelial-to-mesenchymal transition and tubulo-interstitial fibrosis by selectively inhibiting TGF-β-mediated Smad3 phosphorylation. Phytomedicine 42, 207–218. doi:10.1016/j.phymed.2018.03.034
Chen, J. S., Lee, H. S., Jin, J. S., Chen, A., Lin, S. H., Ka, S. M., et al. (2004). Attenuation of mouse mesangial cell contractility by high glucose and mannitol: Involvement of protein kinase C and focal adhesion kinase. J. Biomed. Sci. 11, 142–151. doi:10.1007/BF02256557
Chen, L., Lin, L., Dong, Z., Zhang, L., and Du, H. (2018). Comparison of neuroprotective effect of Forsythia suspensa leaf extract and forsythiaside, one of its metabolites. Nat. Prod. Res. 32, 2705–2708. doi:10.1080/14786419.2017.1374266
Chen, L., Yan, Y., Chen, T., Zhang, L., Gao, X., Du, C., et al. (2019). Forsythiaside prevents β-amyloid-induced hippocampal slice injury by upregulating 2-arachidonoylglycerol via cannabinoid receptor 1-dependent NF-κB pathway. Neurochem. Int. 125, 57–66. doi:10.1016/j.neuint.2019.02.008
Chen, S. J., Lv, L. L., Liu, B. C., and Tang, R. N. (2020). Crosstalk between tubular epithelial cells and glomerular endothelial cells in diabetic kidney disease. Cell Prolif. 53, e12763. doi:10.1111/cpr.12763
Chen, X., Liu, M., and Gilson, M. K. (2001). BindingDB: A web-accessible molecular recognition database. Comb. Chem. High. Throughput Screen 4, 719–725. doi:10.2174/1386207013330670
Chen, Y., Liu, Q., Shan, Z., Zhao, Y., Li, M., Wang, B., et al. (2019). The protective effect and mechanism of catalpol on high glucose-induced podocyte injury. BMC Complement. Altern. Med. 19, 244. doi:10.1186/s12906-019-2656-8
Cheng, G., Zhao, Y., Li, H., Wu, Y., Li, X., Han, Q., et al. (2014). Forsythiaside attenuates lipopolysaccharide-induced inflammatory responses in the bursa of Fabricius of chickens by downregulating the NF-κB signaling pathway. Exp. Ther. Med. 7, 179–184. doi:10.3892/etm.2013.1378
Cheng, L., Li, F., Ma, R., and Hu, X. (2015). Forsythiaside inhibits cigarette smoke-induced lung inflammation by activation of Nrf2 and inhibition of NF-κB. Int. Immunopharmacol. 28, 494–499. doi:10.1016/j.intimp.2015.07.011
Cheng, Y., Wang, D., Wang, F., Liu, J., Huang, B., Baker, M. A., et al. (2020). Endogenous miR-204 protects the kidney against chronic injury in hypertension and diabetes. J. Am. Soc. Nephrol. 31, 1539–1554. doi:10.1681/ASN.2019101100
Coleman, C. I., Weeda, E. R., Kharat, A., Bookhart, B., and Baker, W. L. (2020). Impact of angiotensin-converting enzyme inhibitors or angiotensin receptor blockers on renal and mortality outcomes in people with type 2 diabetes and proteinuria. Diabet. Med. 37, 44–52. doi:10.1111/dme.14107
Craven, P. A., and DeRubertis, F. R. (1989). Protein kinase C is activated in glomeruli from streptozotocin diabetic rats. Possible mediation by glucose. J. Clin. Invest. 83, 1667–1675. doi:10.1172/JCI114066
de Boer, I. H., Rue, T. C., Hall, Y. N., Heagerty, P. J., Weiss, N. S., and Himmelfarb, J. (2011). Temporal trends in the prevalence of diabetic kidney disease in the United States. JAMA 305, 2532–2539. doi:10.1001/jama.2011.861
Doshi, S. M., and Friedman, A. N. (2017). Diagnosis and management of type 2 diabetic kidney disease. Clin. J. Am. Soc. Nephrol. 12, 1366–1373. doi:10.2215/CJN.11111016
Gaulton, A., Hersey, A., Nowotka, M., Bento, A. P., Chambers, J., Mendez, D., et al. (2017). The ChEMBL database in 2017. Nucleic Acids Res. 45, D945–D954. doi:10.1093/nar/gkw1074
Gfeller, D., Grosdidier, A., Wirth, M., Daina, A., Michielin, O., and Zoete, V. (2014). SwissTargetPrediction: A web server for target prediction of bioactive small molecules. Nucleic Acids Res. 42, W32–W38. doi:10.1093/nar/gku293
Gilbert, R. E. (2017). Proximal tubulopathy: Prime mover and key therapeutic target in diabetic kidney disease. Diabetes 66, 791–800. doi:10.2337/db16-0796
Gong, L., Zhou, H., Wang, C., He, L., Guo, C., Peng, C., et al. (2021). Hepatoprotective effect of forsythiaside a against acetaminophen-induced liver injury in zebrafish: Coupling network pharmacology with biochemical pharmacology. J. Ethnopharmacol. 271, 113890. doi:10.1016/j.jep.2021.113890
Goodsell, D. S., Zardecki, C., CostanzoDi, L., Duarte, J. M., Hudson, B. P., Persikova, I., et al. (2020). RCSB Protein Data Bank: Enabling biomedical research and drug discovery. Protein Sci. 29, 52–65. doi:10.1002/pro.3730
Gorji, N., Moeini, R., and Memariani, Z. (2018). Almond, hazelnut and walnut, three nuts for neuroprotection in alzheimer's disease: A neuropharmacological review of their bioactive constituents. Pharmacol. Res. 129, 115–127. doi:10.1016/j.phrs.2017.12.003
Guo, H., Liu, A. H., Ye, M., Yang, M., and Guo, D. A. (2007). Characterization of phenolic compounds in the fruits of Forsythia suspensa by high-performance liquid chromatography coupled with electrospray ionization tandem mass spectrometry. Rapid Commun. Mass Spectrom. 21, 715–729. doi:10.1002/rcm.2875
Habib, S. L. (2013). Diabetes and renal tubular cell apoptosis. World J. Diabetes 4, 27–30. doi:10.4239/wjd.v4.i2.27
Hirose, A., Tanikawa, T., Mori, H., Okada, Y., and Tanaka, Y. (2010). Advanced glycation end products increase endothelial permeability through the RAGE/Rho signaling pathway. FEBS. Lett. 584, 61–66. doi:10.1016/j.febslet.2009.11.082
Hopkins, A. L. (2007). Network pharmacology. Nat. Biotechnol. 25, 1110–1111. doi:10.1038/nbt1007-1110
Huang, C., Lin, Y., Su, H., and Ye, D. (2015). Forsythiaside protects against hydrogen peroxide-induced oxidative stress and apoptosis in PC12 cell. Neurochem. Res. 40, 27–35. doi:10.1007/s11064-014-1461-5
Huang, da.W., Sherman, B. T., and Lempicki, R. A. (2009). Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 4, 44–57. doi:10.1038/nprot.2008.211
Huang, Y., Chi, J., Wei, F., Zhou, Y., Cao, Y., and Wang, Y. (2020). Mitochondrial dna: A new predictor of diabetic kidney disease. Int. J. Endocrinol. 2020, 3650937. doi:10.1155/2020/3650937
Jeruschke, S., Büscher, A. K., Oh, J., Saleem, M. A., Hoyer, P. F., Weber, S., et al. (2013). Protective effects of the mTOR inhibitor everolimus on cytoskeletal injury in human podocytes are mediated by RhoA signaling. PLoS. One 8, e55980. doi:10.1371/journal.pone.0055980
Kang, N., Alexander, G., Park, J. K., Maasch, C., Buchwalow, I., and Luft, F. C. (1999). Differential expression of protein kinase C isoforms in streptozotocininduced diabetic rats. Kidney Int. 56, 1737–1750. doi:10.1046/j.1523-1755.1999.00725.x
Kim, S., Chen, J., Cheng, T., Gindulyte, A., He, J., He, S., et al. (2019). PubChem 2019 update: Improved access to chemical data. Nucleic Acids Res. 47, D1102–D1109. doi:10.1093/nar/gky1033
Kloc, M., Uosef, A., Wosik, J., Kubiak, J. Z., and Ghobrial, R. M. (2019). RhoA pathway and actin regulation of the golgi/centriole complex. Results Probl. Cell. Differ. 67, 81–93. doi:10.1007/978-3-030-23173-6_5
Koya, D. (2014). Dual protein kinase C alpha and beta inhibitors and diabetic kidney disease: A revisited therapeutic target for future clinical trials. J. Diabetes Investig. 5, 147–148. doi:10.1111/jdi.12154
Kvirkvelia, N., McMenamin, M., Warren, M., Jadeja, R. N., Kodeboyina, S. K., Sharma, A., et al. (2018). Kidney-targeted inhibition of protein kinase C-α ameliorates nephrotoxic nephritis with restoration of mitochondrial dysfunction. Kidney Int. 94, 280–291. doi:10.1016/j.kint.2018.01.032
Law, A. H., Yang, C. L., Lau, A. S., and Chan, G. C. (2017). Antiviral effect of forsythoside A from Forsythia suspensa (Thunb.) Vahl fruit against influenza A virus through reduction of viral M1 protein. J. Ethnopharmacol. 209, 236–247. doi:10.1016/j.jep.2017.07.015
Lee, T. S., Saltsman, K. A., Ohashi, H., and King, G. L. (1989). Activation of protein kinase C by elevation of glucose concentration: Proposal for a mechanism in the development of diabetic vascular complications. Proc. Natl. Acad. Sci. U. S. A. 86, 5141–5145. doi:10.1073/pnas.86.13.5141
Lehtonen, S. (2020). Metformin protects against podocyte injury in diabetic kidney disease. Pharm. (Basel) 13, 452. doi:10.3390/ph13120452
Lovisa, S., Zeisberg, M., and Kalluri, R. (2016). Partial epithelial-to-mesenchymal transition and other new mechanisms of kidney fibrosis. Trends Endocrinol. Metab. 27, 681–695. doi:10.1016/j.tem.2016.06.004
Lu, C., Zheng, S. F., and Liu, J. (2020). Forsythiaside A alleviates renal damage in adriamycin-induced nephropathy. Front. Biosci. (Landmark Ed. 25, 526–535. doi:10.2741/4818
Lu, T., Piao, X. L., Zhang, Q., Wang, D., Piao, X. S., and Kim, S. W. (2010). Protective effects of Forsythia suspensa extract against oxidative stress induced by diquat in rats. Food Chem. Toxicol. 48, 764–770. doi:10.1016/j.fct.2009.12.018
Lu, X. J. (2020). DSSR-enabled innovative schematics of 3D nucleic acid structures with PyMOL. Nucleic Acids Res. 48, e74. doi:10.1093/nar/gkaa426
MacIsaac, R. J., Jerums, G., and Ekinci, E. I. (2017). Effects of glycaemic management on diabetic kidney disease. World J. Diabetes 8, 172–186. doi:10.4239/wjd.v8.i5.172
Mao, Q. Q., Xu, X. Y., Cao, S. Y., Gan, R. Y., Corke, H., Beta, T., et al. (2019). Bioactive compounds and bioactivities of ginger (zingiber officinale roscoe). Foods 8, 185. doi:10.3390/foods8060185
Morris, G. M., Huey, R., Lindstrom, W., Sanner, M. F., Belew, R. K., Goodsell, D. S., et al. (2009). AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 30, 2785–2791. doi:10.1002/jcc.21256
Pagtalunan, M. E., Miller, P. L., Jumping-Eagle, S., Nelson, R. G., Myers, B. D., Rennke, H. G., et al. (1997). Podocyte loss and progressive glomerular injury in type II diabetes. J. Clin. Invest. 99, 342–348. doi:10.1172/JCI119163
Pang, X. C., Kang, D., Fang, J. S., Zhao, Y., Xu, L. J., Lian, W. W., et al. (2018). Network pharmacology-based analysis of Chinese herbal Naodesheng formula for application to Alzheimer's disease. Chin. J. Nat. Med. 16, 53–62. doi:10.1016/S1875-5364(18)30029-3
Parameswaran, N., Disa, J., Spielman, W. S., Brooks, D. P., Nambi, P., and Aiyar, N. (2000). Activation of multiple mitogen-activated protein kinases by recombinant calcitonin gene-related peptide receptor. Eur. J. Pharmacol. 389, 125–130. doi:10.1016/s0014-2999(99)00874-2
Piñero, J., Ramírez-Anguita, J. M., Saüch-Pitarch, J., Ronzano, F., Centeno, E., Sanz, F., et al. (2020). The DisGeNET knowledge platform for disease genomics: 2019 update. Nucleic Acids Res. 48, D845–D855. doi:10.1093/nar/gkz1021
Rao, J., Ye, Z., Tang, H., Wang, C., Peng, H., Lai, W., et al. (2017). The RhoA/ROCK pathway ameliorates adhesion and inflammatory infiltration induced by AGEs in glomerular endothelial cells. Sci. Rep. 7, 39727. doi:10.1038/srep39727
Rebhan, M., Chalifa-Caspi, V., Prilusky, J., and Lancet, D. (1998). GeneCards: A novel functional genomics compendium with automated data mining and query reformulation support. Bioinformatics 14, 656–664. doi:10.1093/bioinformatics/14.8.656
Reiser, J., and Sever, S. (2013). Podocyte biology and pathogenesis of kidney disease. Annu. Rev. Med. 64, 357–366. doi:10.1146/annurev-med-050311-163340
Ren, X. D., Zhang, Y. W., Wang, X. P., and Li, Y. R. (2017). Effects of Dangguibuxue decoction on rat glomerular mesangial cells cultured under high glucose conditions. BMC Complement. Altern. Med. 17, 283. doi:10.1186/s12906-017-1774-4
Ru, J., Li, P., Wang, J., Li, B., Huang, C., Li, P., et al. (2014). TCMSP: A database of systems pharmacology for drug discovery from herbal medicines. J. Cheminform. 6 (13), 13. doi:10.1186/1758-2946-6-13
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, 2498–2504. doi:10.1101/gr.1239303
Shimokawa, H., Sunamura, S., and Satoh, K. (2016). RhoA/rho-kinase in the cardiovascular system. Circ. Res. 118, 352–366. doi:10.1161/CIRCRESAHA.115.306532
Steffes, M. W., Schmidt, D., McCrery, R., and Basgen, J. M.International Diabetic Nephropathy Study Group (2001). Glomerular cell number in normal subjects and in type 1 diabetic patients. Kidney Int. 59, 2104–2113. doi:10.1046/j.1523-1755.2001.00725.x
Susztak, K., Raff, A. C., Schiffer, M., and Böttinger, E. P. (2006). Glucose-induced reactive oxygen species cause apoptosis of podocytes and podocyte depletion at the onset of diabetic nephropathy. Diabetes 55, 225–233. doi:10.2337/diabetes.55.01.06.db05-0894
Szklarczyk, D., Morris, J. H., Cook, H., Kuhn, M., Wyder, S., Simonovic, M., et al. (2017). The STRING database in 2017: Quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 45, D362–D368. doi:10.1093/nar/gkw937
Thomas, M. C., Brownlee, M., Susztak, K., Sharma, K., Jandeleit-Dahm, K. A., Zoungas, S., et al. (2015). Diabetic kidney disease. Nat. Rev. Dis. Prim. 1, 15018. doi:10.1038/nrdp.2015.18
Wang, M., Chen, D. Q., Chen, L., Cao, G., Zhao, H., Liu, D., et al. (2018). Novel inhibitors of the cellular renin-angiotensin system components, poricoic acids, target Smad3 phosphorylation and Wnt/β-catenin pathway against renal fibrosis. Br. J. Pharmacol. 175, 2689–2708. doi:10.1111/bph.14333
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, W356–W360. doi:10.1093/nar/gkx374
Wang, Y., Zhang, S., Li, F., Zhou, Y., Zhang, Y., Wang, Z., et al. (2020). Therapeutic target database 2020: Enriched resource for facilitating research and early development of targeted therapeutics. Nucleic Acids Res. 48, D1031–D1041. doi:10.1093/nar/gkz981
Wang, Y., Zhao, H., Lin, C., Ren, J., and Zhang, S. (2016). Forsythiaside A exhibits anti-inflammatory effects in LPS-stimulated BV2 microglia cells through activation of Nrf2/HO-1 signaling pathway. Neurochem. Res. 41, 659–665. doi:10.1007/s11064-015-1731-x
Wishart, D. S., Knox, C., Guo, A. C., Cheng, D., Shrivastava, S., Tzur, D., et al. (2008). DrugBank: A knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 36, D901–D906. doi:10.1093/nar/gkm958
Zeng, X. Y., Yuan, W., Zhou, L., Wang, S. X., Xie, Y., and Forsythoside, Y. J. Fu. (2017). Forsythoside A exerts an anti-endotoxin effect by blocking the LPS/TLR4 signaling pathway and inhibiting Tregs in vitro. Int. J. Mol. Med. 40, 243–250. doi:10.3892/ijmm.2017.2990
Zhan, X., Chen, W., Chen, J., Lei, C., and Wei, L. (2021). Telmisartan mitigates high-glucose-induced injury in renal glomerular endothelial cells (rGECs) and albuminuria in diabetes mice. Chem. Res. Toxicol. 34, 2079–2086. doi:10.1021/acs.chemrestox.1c00159
Zhang, L., Yang, L., Shergis, J., Zhang, L., Zhang, A. L., Guo, X., et al. (2019). Chinese herbal medicine for diabetic kidney disease: A systematic review and meta-analysis of randomised placebo-controlled trials. BMJ Open 9, e025653. doi:10.1136/bmjopen-2018-025653
Zhang, Y., Miao, H., Yan, H., Sheng, Y., and Ji, L. (2018). Hepatoprotective effect of Forsythiae Fructus water extract against carbon tetrachloride-induced liver fibrosis in mice. J. Ethnopharmacol. 218, 27–34. doi:10.1016/j.jep.2018.02.033
Zhao, N., Luo, Q., Lin, R., Li, Q., and Ma, P. (2021). MiR-142-3p ameliorates high glucose-induced renal tubular epithelial cell injury by targeting BOD1. Clin. Exp. Nephrol. 25, 1182–1192. doi:10.1007/s10157-021-02102-y
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, 1523. doi:10.1038/s41467-019-09234-6
Nomenclature
DKD diabetic kidney disease
FA Forsythiaside A powder
HG podocytes incubated in media with 40 mM D-glucose for 48 h
HG + FA podocytes incubated in media with 40 mM D-glucose and varying concentrations of Forsythiaside (2.5, 5, and 10 μg/mL) for 48 h
MG podocytes incubated in media with 5.5 mM D-glucose and 34.5 mM D-mannitol for 48 h
NG podocytes incubated in media with 5.5 mM D-glucose for 48 h
OMIM Online Mendelian Inheritance in Man database
PKC protein-kinase C
PPI protein-protein interaction
PRKCA protein kinase C alpha
ROS reactive oxygen species
TCM traditional Chinese medicine
Keywords: forsythiaside, diabetic kidney disease, network analysis, oxidative stress, podocytopathy, traditional Chinese medicine
Citation: Xu C, Miao H, Chen X and Zhang H (2023) Cellular mechanism of action of forsythiaside for the treatment of diabetic kidney disease. Front. Pharmacol. 13:1096536. doi: 10.3389/fphar.2022.1096536
Received: 12 November 2022; Accepted: 30 December 2022;
Published: 13 January 2023.
Edited by:
Claudio Ferrante, University of Studies G. d’Annunzio Chieti and Pescara, ItalyReviewed by:
Güneş Ak, Selçuk University, TürkiyeXianju Huang, South-Central Minzu University, China
Copyright © 2023 Xu, Miao, Chen 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: Chunmei Xu, m18253161803@163.com; Haiqing Zhang, haiqingzhang7576@163.com