2.1 Network pharmacology analysis
2.1.1 Prediction of AD targets and DN disease genes
We searched the chemical structure and molecular properties of AD from the PubChem database (https://pubchem.ncbi.nlm.nih.gov/) and downloaded its SMILES format file. Using four online databases, namely PharmMapper (http://lilab.ecust.edu.cn/pharmmapper/), SwissTargetPrediction (http://www.swisstargetprediction.ch/), BATMAN-TCM (http://bionet.ncpsb.org/batman-tcm/) and SuperPred (http://prediction.charite.de/), we predicted the possible targets of AD based on its chemical structure, and integrated the results of the four databases, removed the duplicate targets, and obtained the potential target set of AD. To obtain the disease gene information of DN, we used four databases, DrugBank (https://www.drugbank.ca/), PharmGKB (https://www.pharmgkb.org/), GeneCards (https://www.genecards.org/) and OMIM (https://www.omim.org/) for analysis and summary. We entered the keywords “Diabetic Nephropathy” and “Diabetic Nephropathies” into these four databases respectively, obtained the related target information of DN, and used the Uniprot database to convert the protein names of the targets into the corresponding gene names, and performed deduplication and integration, and obtained the final target set of DN.
2.1.2 Intersection targets and network analysis of AD and DN
To obtain the intersection targets of AD and DN, we compared the drug target set of AD and the disease gene set of DN, and obtained the common targets of the two. To analyze the common targets of AD and DN, we used the VennDiagram package of R language, drew the Venn diagram of the targets of AD and DN, and obtained the set of common targets. To construct the protein-protein interaction (PPI) network of the intersection targets of AD and DN, we used the String database to input the intersection targets, queried the PPI information in the database, and set the minimum interaction score to 0.7, obtained the PPI network, and exported it as Cytoscape format, and used Cytoscape software for visualization and analysis. Using Cytoscape software, we used the intersection targets as nodes, and AD and DN as attributes, and drew the network graph. Using the cytoHubba plugin (https://apps.cytoscape.org/apps/cytohubba), according to the target interaction network of AD-DN, we used various algorithms, such as Degree, Betweenness, Closeness, EPC, MCC, etc., to calculate the centrality indicators of each target, and comprehensively considered various indicators, and screened out the core target set of AD-DN. And according to the degree value of the nodes, we used the cytoHubba function module to screen out the 7 core targets with the highest Degree value.
2.1.3 GO and KEGG enrichment analysis.
Using the DAVID database (https://david.ncifcrf.gov/), according to the common target set and core target set of AD-DN, we performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, and analyzed the biological processes (BP), cellular components (CC), molecular functions (MF) and signaling pathways (Pathway) involved by the common targets and core targets of AD-DN, and drew bar charts, bubble charts, mulberry charts and composite charts and other visualization forms to show the results of enrichment analysis.
2.2 In vitro experiment validation
2.2.1 Cell culture and treatment
Mouse glomerular podocyte MPC-5 cells were purchased from Shanghai Cell Bank and cultured in Dulbeccos Modified Eagle Medium (DMEM) containing 10% fetal bovine serum, and incubated in a 37 °C, 5% CO2 incubator. The cells were passaged to the 3rd to 5th generation for the experiment. Andrographolide AD was purchased from Shanghai Zhongyao Co., Ltd. of China Pharmaceutical Group, dissolved in dimethyl sulfoxide (DMSO) and diluted to different concentrations of reserve solution, and stored at 8 °C. WP1066 was purchased from Selleck Company in the United States, dissolved in DMSO and diluted to a 10 μM reserve solution, and stored in an 8 °C refrigerator. The cells were divided into the following 6 experimental groups: normal group (NC group), cultured with DMEM containing 5.5 mM glucose; diabetic nephropathy model group (M group), cultured with DMEM containing 30 mM glucose; M-AD 10 μM group (M-10 group), 10 μM of AD was added on the basis of the M group; M-AD 20 μM group (M-20 group), 20 μM of AD was added on the basis of the M group; M-AD 40 μM group (M-40 group), 40 μM of AD was added on the basis of the M group; M-AD 80 μM group (M-80 group), 80 μM of AD was added on the basis of the M group; M-AD 80 μM-WP1066 group (M-80-WP group), 10 μM of WP1066 was added on the basis of the M-80 group. After 48 h of cell culture, the cells and supernatant were collected and subjected to subsequent experiments. The WP1066 (10 μM, this concentration can effectively inhibit the phosphorylation and activation of STAT3, and will not have a significant effect on cell viability) used in this cell experiment is a novel, cell-permeable STAT3 inhibitor, which can inhibit the phosphorylation and activation of STAT3, and thus block the downstream signaling pathways, such as PI3K/AKT, MAPK, ERK, etc.
2.2.2 MTT assay to detect cell proliferation ability
The podocytes were passaged to a 96-well plate, and after 48 h of group intervention, 20 μL of MTT (5 mg/mL) solution was added, incubated at 37 °C for 4 h, the culture medium was removed, 150 μL of DMSO was added, and shaken for 10 min, and the absorbance (OD value) of each well was measured by an enzyme-linked immunosorbent assay (wavelength 490 nm), and the cell proliferation rate was calculated.
2.2.3 ELISA experiment
ELISA method to detect the levels of inflammatory factors (TNF-α, IL-1β), Advanced Glycation End Products (AGEs) and vascular active substances (PGE2, TXB2) in the cell supernatant. The podocytes were passaged to a 6-well plate, and after 48 h of group intervention, the cell supernatant was collected, and operated according to the instructions of the ELISA kit, and the OD value of each well was measured by an enzyme-linked immunosorbent assay (wavelength 450 nm), and the concentration of each indicator was calculated according to the standard curve.
2.2.4 Oxidative stress level detection experiment
DCFH-DA method to detect the level of ROS in cells, GSH/GSSG ratio to detect the level of GSH in cells. The podocytes were passaged to a 6-well plate, and after 48 h of group intervention, the original culture medium was replaced with serum-free culture medium containing 10 μM DCFH-DA, incubated at 37 °C for 30 min, and the original culture medium was replaced with serum-free culture medium containing 10 μM NEM, incubated at 37 °C for 15 min, and the cell lysate was collected, and operated according to the instructions of the GSH/GSSG detection kit, and the OD value of each well was measured by an enzyme-linked immunosorbent assay (wavelength 405 nm), and the data of each indicator was calculated according to the standard curve, and the GSH/GSSG ratio was calculated.
2.2.5 Detection of expression and activity of core targets in cells
Using laser confocal microscopy, we studied the effects of the expression changes of the core target STAT3 in cells on the expression of IL-6 and PIK3CA. The cells were seeded on glass coverslips, 2*104 per well, and cultured for 24 h. After treatment according to the previous method, on the third day after treatment, the cells were fixed with 4% polyformaldehyde for 15 min, permeabilized with 0.1% Triton X-100 for 10 min, blocked with 5% bovine serum for 1 h, and incubated with anti-STAT3, anti-IL-6, and anti-PIK3CA primary antibodies (1:200, Abcam, UK) overnight. Then, the cells were incubated with fluorescently labeled secondary antibodies (1:500, Abcam, UK) for 1 h, stained with DAPI for 10 min, washed three times with PBS, and observed and photographed with a laser confocal microscope (Leica, Germany). The expression and localization of STAT3, IL-6, and PIK3CA in the cells were analyzed. 2.2.6 Immunoprecipitation experiment
Combined with Western blot detection method, we studied whether the expression of STAT3 in cells was strongly correlated with IL-6 and PIK3CA. The cells were collected, lysed with RIPA buffer, centrifuged for 10 min, and the supernatant was taken. The protein concentration was determined by BCA method, and equal amounts of protein were taken. The proteins were incubated with anti-STAT3 primary antibody (1:200, Abcam, UK) overnight, and then protein A/G magnetic beads (Beyotime, China) were added and incubated at room temperature for 2 h. The magnetic beads were adsorbed with a magnetic rack, the supernatant was removed, the proteins were separated by SDS-PAGE, transferred to a membrane, and incubated with anti-IL-6 and anti-PIK3CA primary antibodies (1:1000, Abcam, UK) and corresponding secondary antibodies (1:5000, Abcam, UK). The membrane was developed with ECL luminescent liquid, and the gray value of the bands was analyzed with ImageJ software to reflect the interaction of STAT3 with IL-6 and PIK3CA in the cells.
2.2.7 Protein expression level detection
Western blot method was used to detect the expression changes of AGEs-RAGE signaling pathway-related proteins (RAGE, p38 MAPK, NF-κB), autophagy-related proteins (BINP3, NIX, FUNDC1, PINK1, Parkin) and core target proteins (STAT3, IL6, PIK3CA, JUN, CASP3, MAPK1 and ERBB2) in the cell lysate. The podocytes were passaged to a 6-well plate, and after 48 h of group intervention, the cell lysate was collected, and the protein concentration was measured by the BCA method, and the protein was separated by SDS-PAGE electrophoresis, transferred to the membrane, sealed with 5% skim milk, and added with the corresponding primary antibody (1:1000 dilution) and incubated overnight at 4 °C, and added with HRP-labeled secondary antibody (1:5000 dilution) and incubated at room temperature for 1 h, and colored with ECL luminescent liquid, and photographed with Gel Doc XR+ imaging system, and the gray value was analyzed by Image J software, and the relative protein expression level was calculated with β-actin as the internal reference.
2.2.8 RT-PCR experiment
RT-PCR method was used to detect the transcription level of AGEs-RAGE signaling pathway-related genes (RAGE, p38 MAPK, NF-κB), autophagy-related genes (BINP3, NIX, FUNDC1, PINK1, Parkin) and core targets (STAT3, IL6, PIK3CA, JUN, CASP3, MAPK1 and ERBB2) in the cells. The podocytes were passaged to a 6-well plate, and after 48 h of group intervention, the total RNA of the cells was collected, extracted by Trizol method, and the concentration and purity of RNA were measured by NanoDrop 2000, and reverse transcription was performed by PrimeScript RT reagent Kit to obtain cDNA. Real-time fluorescence quantitative PCR was performed with SYBR Premix Ex Taq II, and the fluorescence signal was detected by ABI 7500 system, and the relative gene expression level was calculated by 2-ΔΔCt method, with GAPDH as the internal reference. And the primer sequences in this study were shown in Table.1.
2.2.9 Statistical analysis
GraphPad Prism (version 22.0, Inc-La Jolla, USA) was used for statistical analysis. All experiments were performed at least three times, and the gray values of three independent experiments were read by ImageJ-pro software. All data were analyzed by Student’s t-test or one-way analysis of variance (one-way ANOVA), and were normally distributed, and passed the homogeneity of variance test. Data were expressed as mean ± SD of three repeated cell cultures, n=3. In addition, SPSS software (version 19, SPSS Inc., USA) was used for statistical analysis, and principal component analysis was used to define the most important parameters. *, indicates that the data of this group has statistical difference with another control group (P < 0.05); **, indicates that the data of this group has significant statistical difference with other groups (P < 0.01).