Network Pharmacology Combined with Bioinformatics to Investigate the Mechanisms and Molecular Targets of Astragalus Radix-Panax notoginseng Herb Pair on Treating Diabetic Nephropathy

Background Astragalus Radix (AR)-Panax notoginseng (PN), a classical herb pair, has shown significant effects in treating diabetic nephropathy (DN). However, the intrinsic mechanism of AR-PN treating DN is still unclear. This study aims to illustrate the mechanism and molecular targets of AR-PN treating DN based on network pharmacology combined with bioinformatics. Materials and Methods The Traditional Chinese Medicine Systems Pharmacology database was used to screen bioactive ingredients of AR-PN. Subsequently, putative targets of bioactive ingredients were predicted utilizing the DrugBank database and converted into genes on UniProtKB database. DN-related targets were retrieved via analyzing published microarray data (GSE30528) from the Gene Expression Omnibus database. Protein-protein interaction networks of AR-PN putative targets and DN-related targets were established to identify candidate targets using Cytoscape 3.8.0. GO and KEGG enrichment analyses of candidate targets were reflected using a plugin ClueGO of Cytoscape. Molecular docking was performed using AutoDock Vina software, and the results were visualized by Pymol software. The diagnostic capacity of hub genes was verified by receiver operating characteristic (ROC) curves. Results Twenty-two bioactive ingredients and 189 putative targets of AR-PN were obtained. Eight hundred and fifty differently expressed genes related to DN were screened. The PPI network showed that 115 candidate targets of AR-PN against DN were identified. GO and KEGG analyses revealed that candidate targets of AR-PN against DN were mainly involved in the apoptosis, oxidative stress, cell cycle, and inflammation response, regulating the PI3K-Akt signaling pathway, cell cycle, and MAPK signaling pathway. Moreover, MAPK1, AKT1, GSK3B, CDKN1A, TP53, RELA, MYC, GRB2, JUN, and EGFR were considered as the core potential therapeutic targets. Molecular docking demonstrated that these core targets had a great binding affinity with quercetin, kaempferol, isorhamnetin, and formononetin components. ROC curve analysis showed that AKT1, TP53, RELA, JUN, CDKN1A, and EGFR are effective in discriminating DN from controls. Conclusions AR-PN against DN may exert its renoprotective effects via various bioactive chemicals and the related pharmacological pathways, involving multiple molecular targets, which may be a promising herb pair treating DN. Nevertheless, these results should be further validated by experimental evidence.


Introduction
Diabetes mellitus (DM) is a common metabolic disorder characterized by chronic hyperglycemia due to defective insulin secretion and/or utilization, which has been considered an enormous global health concern [1]. According to the International Diabetes Federation data, 352 million people between the ages of 20 and 65 years are suffering from DM worldwide, which will occur in up to 486 million by 2045 [2]. Diabetic nephropathy (DN), a major microvascular complication of DM, is a leading cause of end-stage renal disease [3]. Research studies indicate that 20%-40% of patients suffering from DM develop DN and patients with DN have a higher risk of death with 3 to 12 times than those with DM alone [4,5]. Advances have been made in the treatment of DN in recent years with the introduction of lifestyle changes, including exercise and diet, and therapies such as angiotensin-converting enzyme inhibitors, angiotensin receptor blockers, sodium-glucose cotransporter 2 inhibitors, and glucagon-like peptide-1 receptor [3,6]. Unfortunately, these drugs were reported to have adverse effects, such as diabetic ketoacidosis, gastrointestinal, and hepatotoxicity [7,8], and there are still no pharmacotherapies available to cure or reverse disease progression. Highlighting the requirement for novel therapeutic with better efficacy and safe treating DN is still necessary.
Chinese herbal medicine with various bioactive chemicals and the related pharmacological pathways, involving multiple cellular and molecular targets, and fewer side effects has been frequently and extensively used for thousands of years in China [9], which may be to provide beneficial and effective therapy for DN. According to the theory of traditional Chinese medicine (TCM), the basic pathology of DN is qi deficiency and blood stasis [10]. Hence, tonifying qi and activating blood circulation are the basic therapeutic principles for DN. Astragalus Radix (AR, Huang qi in Chinese) was initially recorded in Shennong Bencao Jing, which has been extensively known to exert the tonifying Qi effect. Panax notoginseng (PN, San qi in Chinese) was found in Bencao Gangmu with multiple beneficial effects, especially activating blood circulation and decreasing swelling to relieving pain. e herb pair is considered the most basic form of TCM therapy containing two special Chinese herbal medicines, which plays a vital role in connecting mutual enhancement, assistance, and restraint [11]. AR-PN, a classic highly valued and widely used herb pair, is well known to possess functions with tonifying qi and activating blood circulation. Modern pharmacological research shows that AR treatment ameliorated the severity of DN by inhibiting inflammation-related gene IL-1β and IL-18 expression and fibrosis indexes in DN rats [12]. Likewise, it has been confirmed that AR possesses antioxidative stress effects [13]. Previous research demonstrated that PN can significantly reduce albuminuria, proteinuria, and serum creatinine and improve the metabolism of serum lipids in patients with DN [14]. In addition, both in vivo and in vitro research studies confirmed that the combination of AR and PN could strengthen the renoprotection effects by inhibiting inflammation, alleviating oxidative stress, and upregulating autophagy [15]. Taken together, the AR-PN herb pair playing a renoprotection role in DN is mainly aimed at anti-inflammation and antioxidative stress. However, the underlying pharmacological mechanism of the AR-PN herb pair in the treatment of DN remains to be illustrated. A novel approach to clarify the pharmacological mechanisms of the efficacy of the AR-PN herb pair on DN is necessary.
Network pharmacology was first proposed by Hopkins in 2007, which provides a convenient and systematic method to reveal the complex mechanisms of drug treating disease through identifying core targets based on the theory of systems biology and network pharmacology [16]. Moreover, network pharmacology is also contributed to integrate and extract possible signaling pathways that drugs prevent and treat diseases to improve the therapeutic effect of drugs, and the pattern is "drugtarget-disease." It has been applied to explore the mechanisms of TCM treating some diseases, including diabetes mellitus [17] and ischemic stroke [18], showing that applying network pharmacology to elucidate the mechanisms of bioactive components and the potential of TCM is an emerging trend. Bioinformatics is the application of tools of computation and analysis to the capture and interpretation of biological data, which is an interdisciplinary field of computer science, mathematics, physics, and biology [19]. Currently, bioinformatics has been increasingly recommended in the study of biological fields, particularly in genomics and proteomics. miRNAs involve the fine-tuning of mRNA abundance via binding to the 3′ untranslated region of a target mRNA and results in translational repression. Increasing evidence suggests that the downregulation or upregulation of miRNAs is closely associated with podocyte injury, inflammation, accumulation of the extracellular matrix, fibrosis, and so on in DN, concerning their potential as biomarkers and miRNA modulation as a therapeutic option for DN [20,21]. Besides, some microarray data analysis studies related to DN have been performed and numerous DEGs have been identified, which may be potential biomarkers and target candidates in DN [22]. Except for DN, bioinformatics strategies are also involved in identifying the core genes of other diseases, such as diabesity [23], familial hypercholesterolemia [24], lung squamous cell carcinoma [25], and prostate cancer [26], indicating that bioinformatics plays vital roles in predicting potential biomarkers and providing therapeutic options. Collectively, network pharmacology integrated with bioinformatics is a novel strategy to study ingredient identification, target prediction, and therapeutic options.
Consequently, we integrated network pharmacology with Gene Expression Omnibus (GEO) database to further identify the comprehensive pharmacological mechanisms, screen for potential molecular targets, and explore underlying pathways of AR-PN herb pair on treating DN, which may provide a guideline in the further research and development. e flowchart of this study is shown in Figure 1.

Screening of Active Ingredients of AR-PN Herb Pair.
All active ingredients of the AR-PN herb pair were collected from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database (http://tcmspw.com/tcmsp.php), which is a unique systematic pharmacology platform to display the relationships between components, targets, and diseases [27]. e absorption, distribution, metabolism, and excretion (ADME) model [28], which was an important contributor to predict the pharmacokinetic properties of chemical ingredients, was used to screen the compounds of the AR-PN herb pair. Two vital parameters, i.e., oral bioavailability (OB) and druglikeness (DL), among all ADMErelated properties were used to identify chemical ingredients of the AR-PN herb pair. e OB refers to the relative amount of a drug absorbed into the systemic circulation after extravascular administration [29]. DL is used to estimate the similitude physiochemical properties or/and structural features between an ingredient and certified drug candidates [30], which is calculated by the following formula: where "X" represents the active ingredients of the AR-PN herb pair and "Y" stands for the average DL index of all molecules that are collected from the DrugBank database (https://www.drugbank.ca) [31]. In our present study, en, probe IDs were applied to identify the corresponding genes. Subsequently, differentially expressed genes (DEGs) between DN glomeruli tissues and controlled normal glomeruli tissues were screened by the limma package of R software. DEGs with P < 0.05 and |log2 fold change (FC)| > 1 were selected as genes of DN and were visualized using a volcano plot.

Construction of Protein-Protein Interaction Network.
Protein-protein interaction (PPI) is the core of almost all biological processes, which present the interaction of proteins and may provide potential targets for therapeutic intervention [34]. e PPI networks of the AR-PN herb pair putative targets as well as DN-related DEGs were constructed and visualized using the BisoGenet [35], a plugin of Cytoscape 3.8.0, which were established based on the currently available PPI databases, including Biological General Subsequently, a merged network was performed according to the intersection data of the two networks.

Network Topological Characteristics Analysis.
e three crucial topological properties, including betweenness centrality (BC), closeness centrality (CC), and degree centrality (DC), were used to calculate and analyze topological properties of every node in the merged interaction network using the CytoNCA [36], a plugin of Cytoscape 3.8.0. e node is more important, and the received quantitative value is higher. First, the degree of centrality was calculated and the network was extracted according to more than twice the median degree of DC. en, more than the median values of DC, BC, and CC as the threshold values were used to further extract hub nodes in the network. Ultimately, a core subnetwork was obtained for the ensuing analysis.

GO and KEGG Pathway Enrichment
Analysis. GO, a freely available public resource, provides a standardized language to systematically describe the functional properties of gene products from all species, which consists of three ontologies, namely, biological process (BP), cellular component (CC), and molecular functional (MF) [37]. BP describes the functioning of biological objectives to which they contribute, CC describes the place of a gene product within the cell and their extracellular environments, and MF represents the biochemical activities, such as binding and catalysis, of gene products [38]. KEGG is a highly complex network structure, which has been an exceptionally comprehensive site for describing metabolic pathways, gene signaling networks, and practical application of genomic information [39]. e biological interpretations of the hub genes in the core network were reflected by GO and KEGG pathway enrichment analysis using ClueGO, a plugin of Cytoscape 3.8.0, which is utilized to visualize nonrepetitive biological terms as functionally grouped networks. It is worth noting that the P value less than or equal to 0.01 was set as the threshold in both the GO or KEGG functional categories. Furthermore, the component-target-pathway network was established and visualized using Cytoscape 3.8.0 software.

Molecular
Docking. Molecular docking, a method of analyzing the binding site and affinity between the conformation of molecules and macromolecular targets, has been considered a crucial technique for structure-based drug discovery [40]. e core DN-related genes targeted by the AR-PN herb pair bioactive components was identified according to the component-target-pathway network, which were further validated by molecular docking with the experimentally verified AR-PN herb pair bioactive components possessing renoprotective effects on DN, including quercetin [41], kaempferol [42], formononetin [43], and isorhamnetin [44]. First, the mol2 formats of these bioactive components were obtained from the TCMSP database and PubChem database, which were converted to the format of pdb and pdbqt via Pymol software and AutoDock Tools 1.5.6, respectively. Consequently, the 3-dimensional structures of potential receptors saved in pdb format were downloaded from the PDB database (https://www.wwpdb.org/). e removal of the ligand and water in the protein was performed via AutoDock Tools 1.5.6., which was saved in pdbqt format. Subsequently, molecular docking was performed to evaluate the combination mode and affinity between the ligand-receptor and molecules by AutoDock Vina, which was run with default settings. e smaller the lowest energy value is, the more stable and the stronger affinity between ligand and receptor binding. Ultimately, we selected the top two receptor proteins with the lowest energy value and the ligand for visualization and construction by Pymol software. 4 Evidence-Based Complementary and Alternative Medicine

Receiver Operating Characteristic (ROC) Curves for Hub
Genes. To evaluate the diagnostic accuracy of each hub gene for DN, ROC curves were plotted and area under the curve (AUC) values were calculated, which were performed with "pROC" package of R software 4.0.3. P < 0.05 was considered to be significant.  Figure 2.

Identification of DN-Related DEGs.
Ultimately, 850 DNrelated DEGs were identified between DN glomeruli tissues and controlled normal glomeruli tissues with |log2 FC| > 1 and P < 0.05 by analyzing the sequencing data downloaded from the GEO database. A volcano plot of the distribution of DEGs is exhibited in Figure 3, which includes 253 upregulated genes dotted in red, 597 downregulated genes dotted in blue, and other genes without significant differences dotted in black.

PPI Network Construction, Merging, and Topology
Analysis. To illustrate the interactions and underlying mechanisms between AR-PN herb pair putative targets and DN-related DEGs, we visually constructed PPI networks of AR-PN herb pair putative targets and DN-related DEGs, respectively. e PPI network of AR-PN herb pair putative targets contained 6063 nodes and 149,352 edges ( Figure 4(a)), while the PPI network of DN-related DEGs contained 8073 nodes and 181,168 edges ( Figure 4(b)). In the PPI network graph, nodes represent interacting proteins, and edges represent interactions. A new PPI network was generated after merging these two PPI networks to identify the candidate targets for the AR-PN herb pair in the treatment of DN, which comprised 4390 nodes and 123,047 edges ( Figure 4(c)). Subsequently, we analyzed the topological properties of the merged PPI network based on BC, CC, and DC key parameters and extracted targets above twofold median values of DC as well as more than the median values of BC and CC. erefore, nodes with DC that were above twofold median values (DC > 72) of all nodes were the first extracted, which included 997 nodes and 44,709 edges ( Figure 4(d)). Whereafter, nodes with DC > 226, BC > 0.000528, and CC > 0.45 of these 997 nodes were the second identified core targets containing 115 nodes and 2,334 edges ( Figure 4(e)). e information of 115 candidate targets sorted in descending order based on the value of the degree is presented in Table 2.

GO and KEGG Pathway Enrichment Analysis.
To further evaluate better the molecular mechanism of the AR-PN herb pair on DN, GO and KEGG pathway enrichment analyses associated with 115 candidate targets were performed using the Cytoscape plugin ClueGO, which included GO biological process and the signaling pathway. e GO biological processes were mainly involved in regulating DNA-binding transcription factor activity, regulating protein modification by small protein conjugation or removal, rhythmic process, proteasome-mediated ubiquitin-dependent protein catabolic process, and cellular response to heat, and so on ( Figure 5(a)). According to the KEGG pathway analysis, the signaling pathways were mainly focused on the PI3K-Akt signaling pathway, viral carcinogenesis, cell cycle, micro-RNAs in cancer, and hepatitis B, and so on ( Figure 5(b)).

Molecular Docking Analysis.
According to the network of component-target-pathway, the core targets (MAPK1, AKT1, GSK3B, CDKN1A, TP53, RELA, MYC, GRB2, JUN, and EGFR) were selected for molecular docking with quercetin, kaempferol, isorhamnetin, and formononetin bioactive components verified experimentally via AutoDock Vina 1.1.2 software. As shown in Table 3, all the core targets mentioned above had a great binding affinity with quercetin, kaempferol, isorhamnetin, and formononetin. Furthermore, the MAPK1 and JUN target proteins were molecularly docked with quercetin, kaempferol, isorhamnetin, and formononetin active components with the lowest energy value that was less than -8.0 kcal·mol −1 . erefore, we selected the MAPK1 and JUN receptor proteins with the lowest energy value and the ligands for visualization and construction by Pymol software (Figure 7).

Discussion
Microarrays are a powerful and effective high-throughput technology and provide information on gene profiling data to researchers to explore biomarker, molecular pathways, target selectivity, and development of prognostic, and so on involved in complex disorders [45]. ey have been widely applied to investigate the key molecule biomarkers and pathways of various diseases, including ovarian cancer [46], multiple myeloma [47], systemic lupus erythematosus (SLE) [48], and idiopathic pulmonary fibrosis [49]. DN is the most common microvascular complication of DM, which has been considered the major cause of ESRD worldwide [3].
Nevertheless, effective treatments remain scarce. Previous studies have been demonstrated that AR and PN play a critical role in treating DN [12][13][14]50], suggesting AR-PN may be a potential herb pair for the treatment of DN. For all we know, this is the first study to perform a comprehensive analysis of network pharmacology combined with gene expression profiling to further explore the underlying pharmacological mechanisms and therapeutic targets of the AR-PN herb pair on DN. Our findings demonstrated that a total of 22 bioactive compounds and 115 core targets for the AR-PN herb pair in the treatment of DN were identified. According to the KEGG pathway enrichment, the AR-PN herb pair played a therapeutic role in DN through regulating the PI3K-Akt signaling pathway, cell cycle, and MAPK signaling pathway, and so on. Furthermore, MAPK1, AKT1, GSK3B, CDKN1A, TP53, RELA, MYC, GRB2, JUN, and EGFR were regarded as core targets of the AR-PN herb pair against DN according to the component-target-pathway network. What is more, molecular docking analysis validated that these core target genes showed strong binding interactions and affinity with quercetin, kaempferol, isorhamnetin, and formononetin bioactive compounds of the AR-PN herb pair.
As shown in the component-target network (Figure 2), 22 compounds in the AR-PN herb pair with possible efficacy against DN were identified, many of which shared targets and had a synergistic effect, indicating that the AR-PN herb   Evidence-Based Complementary and Alternative Medicine pair exerts effects in various bioactive components and multiple cellulars and molecular targets way. Quercetin, a natural flavonoid, is well known to possess multiple pharmacological functions, including but not limited to antioxidant, anti-inflammatory, anticancer, antidiabetic, neuroprotective, and antiobesity [51,52]. e use of quercetin contributes to a meaningful improvement in disease symptoms, renal hypertrophy index, renal histopathology, oxidative stress, and blood glucose in a rat model of DN [41]. Furthermore, quercetin treatment is beneficial to antioxidant capacity and vasoprotective effects against STZ-induced ER stress of diabetic rats via reducing the lipid peroxidation and increasing the expression of vascular endothelial growth factor (VEGF) and its receptor, VEGFR2, and nitric oxide (NO) [53]. Accumulating evidence also demonstrated that the nephroprotective effects of quercetin are related to inhibiting renal tubular epithelial-mesenchymal transition and renal fibrosis and reducing the production of cytokines in DN, including IL-6, TNF-α, and IL-1β [54][55][56]. What is more, quercetin inhibited glomerular mesangial cell proliferation and epithelial-mesenchymal transition in DN via reactivation of the Hippo pathway and inhibition of the TGF-β/PI3K/Akt pathway, respectively [57,58]. Kaempferol, a natural polyflavonol, also exhibits nephroprotective effects in addition to its antioxidant and anti-inflammatory properties [59]. Studies demonstrated that kaempferol plays vital nephroprotective roles associating with inhibiting oxidative stress, proinflammatory cytokines (TNF-α and IL-1β), and fibrosis in DN via inhibiting hyperglycemia-induced activation RhoA/Rho-kinase [42]. Moreover, kaempferol exerts antioxidant potential in DN rats by the upregulation of the Nrf-2/HO-1 axis [60]. Formononetin possesses multiple biological activities, such as antioxidant, anti-inflammatory, and regulating immune responses [61,62]. Formononetin improved oxidative stress and prevented the progression of renal fibrosis in DN by activating the Nrf2/ARE signaling pathway [43]. Experimental studies have shown that isorhamnetin has renoprotective effects on DN through improving fasting blood glucose, increasing autophagosomes in renal tissues, and suppressing miRNA regulation of autophagy genes [44]. Besides, isorhamnetin inhibits oxidative stress, inflammation, and apoptosis in STZ-induced diabetic rats [63].       Figure 6: Component-target-pathway network of AR-PN herb pair against DN. The larger size and darker the color of nodes were shown, the higher the degree value was represented. Rhombus represents targets, triangles represent ingredients of AR-PN herb pair, and arrowhead represents signal pathways.  which were considered as core targets that the AR-PN herb pair plays therapeutic roles in DN. To further illustrate the mechanism of AR-PN herb pair on DN, we conducted GO and KEGG pathway enrichment analyses of 115 core targets. As shown in Figure 5(a), the BP terms of AR-PN herb pair on DN were mainly related to DNA synthesis, transcription, protein catabolism, cell cycle, apoptosis, oxidative stress, and inflammation regulation. Interestingly, apoptosis, oxidative stress, and inflammation are deemed as the important pathogenesis of DN [3]. KEGG enrichment analysis showed that the AR-PN herb pair may alleviate the development of DN through the action of core targets in multiple signaling pathways ( Figure 5(b)), including PI3K-Akt signaling pathway, cell cycle, and MAPK signaling pathway. Among them, the enriched targets of the PI3K-Akt signaling pathway are markedly higher than others. Phosphoinositide-3-kinase (PI3K) can promote the phosphorylation of protein kinase B (PKB, Akt), which involves diversified cell activities, e.g., cell proliferation, protein synthesis, glycolysis, autophagy, and apoptosis [64]. Previous studies demonstrated that restoring podocyte autophagy, attenuating podocyte apoptosis, and inhibiting MC proliferation of DN were closely associated with inhibiting the PI3K/AKT signaling pathway [65][66][67]. Furthermore, notoginsenoside R1, which is extracted from PN, increased the phosphorylation of both PI3K and Akt and ameliorated podocyte injury in DN rats by activating the PI3K/Akt signaling pathway [68]. Cell cycle refers to order progression from one cell cycle phase to the next, involving a family of cyclins and their associated kinases, and mediates cell proliferation, which will be ceased if its checkpoint secures, such as Rad3-related protein kinase and checkpoint kinase 1, causing severe problems [69]. Moreover, the cell cycle is regulated by several cellular pathways, e.g., PI3K/Akt signaling pathway that is famous for regulating cell proliferation [64]. MAPK pathway acts on inhibiting cell apoptosis and promoting proliferation and anti-inflammatory effects in DN rats, which is associated with MAPK1 [70]. Viral carcinogenesis accelerates DNA damage and virus integration, which are closely related to the progression of human cancers [71], as well as micro-RNAs in cancer [72], suggesting AR-PN herb pair also exerts anticancer effects. us, our identified core targets are highly related to the signaling pathways associated with cell proliferation, apoptosis, and cell cycle. Furthermore, the network of component-target-pathway showed that MAPK1, AKT1, GSK3B, CDKN1A, TP53, RELA, MYC, GRB2, JUN, and EGFR were core targets in the network due to possessing the higher degree ( Figure 6), suggesting that these core targets were identified as core genes that were strongly associated with the pathogenesis of AR-PN herb against DN. MAPK1, also known as an extracellular signal-regulated protein kinase, is a member of the MAPKs family that is contributed to enhancing cellular proliferation and participating in cell metabolism, proliferation, apoptosis, differentiation, and survival in physiological and pathological processes [73]. It has been demonstrated that downregulated expression of MAPK1 plays a protective role in lipopolysaccharide-induced podocyte damage in DN [74]. Many studies have observed that the expression of MAPK1 was increased in various cancers, such as breast cancer [75], hepatocellular carcinoma [76], and gastric cancer [77], which are involved in inhibiting cell proliferation, cell cycle, and cell apoptosis. Furthermore, they found that the upregulation of MAPK1 was negatively regulated by miR-20a [75], miR-193b [76], and miR-378 [77], respectively. is suggests that MAPK1 vital gene plays essential roles in the development and progression of DN, which could participate in cell metabolism, cell proliferation, and apoptosis of podocyte, mesangial cell, etc. AKT1 is an important factor to regulate inflammation and apoptosis in kidney disease [78]. e phosphorylation of AKT1 can directly promote mesangial cell (MC) proliferation, basement membrane thickening, renal tubular epithelial cell transdifferentiation, and podocyte injury that are well known to play a critical role in the pathogenic mechanisms of DN [67,79]. Evidence revealed that GSK3B is phosphorylated by AKT at serine 9, resulting in facilitating the expression of inflammatory cytokines [80]. Moreover, the inhibition of GSK3B induces human cell proliferation [81]. Kumar et al. found that AKT1 was one of the essential genes in the development and progression of SLE, with regard to increasing autoimmune responses via different signaling pathways [48]. TP53, a tumor suppressor, mediates glucose and lipotoxicity and causes a stop of the cell cycle, apoptosis, DNA repair, and metabolic changes, which is inhibited by microRNA-770-5p through regulating podocyte apoptosis [82]. is indicates that AKT1, GSK3B, TP53, and RELA might be associated with cell proliferation, cell cycle, and apoptosis of DN. e activation of JUN N-terminal kinase is associated with inflammatory factor secretion during DN [83]. Furthermore, JUN is a highly unstable protein due to polyubiquitination by GSK3 at Ser-243, and its level is inversely associated with GSK3 activity in mammalian cells that have entered the cell cycle [84], showing that JUN interacts with GSK3 via polyubiquitination. Nevertheless, the roles of JUN in DN require further research. Epidermal growth factor receptor (EGFR) belongs to the family of receptor tyrosine kinase ErbB receptors, and the activation of which can result in phosphorylation on specific tyrosine residues, such as ErbB2, ErbB3, or ErbB4, and further regulate cell proliferation, differentiation, and apoptosis [85]. e activation of EGFR leads to podocyte injury and loss in DN, while EGFR deletion in podocytes ameliorates glomerular injury and attenuates the progression of DN [86]. Previous studies also demonstrated that EGFR inhibition plays a nephroprotective role in DN via decreasing reactive oxygen species (ROS), inhibiting endoplasmic reticulum(ER) stress, alleviating podocyte injury, and increasing autophagy [87,88]. RELA is one of the transcription factor genes in the NF-κB signaling pathway that contributes to human disease processes, notably inflammatory diseases and cancer [89]. Site-specific phosphorylation has been proven to be a prerequisite for additional modifications regulating RELA activity, especially acetylation and ubiquitination [90]. A previous study showed that the phosphorylation of RELA at S536 may be related to regulating NF-kB activity, which was further corroborated in a mouse model that knocked down RELA S534A, a nonphosphorylatable mutant version of the human RELA S536 homolog in murine [91]. Besides, RELA also participates in the process of cell proliferation [92]. However, the role of RELA in DN remains unclear and requires further research. To explore the direct or indirect interactions between core target genes, an interrelation analysis for them was performed. As shown in Figure 9, the core target genes interacted with each other via multiple signaling pathways mainly associated with the PI3K-Akt signaling pathway, MAPK signaling pathway, cell cycle, apoptosis, microRNAs in cancer, and viral carcinogenesis. e core genes MAPK1, AKT1, TP53, RELA, MYC, GRB2, and EGFR are involved directly or indirectly in the PI3K-Akt signaling pathway regulating cell proliferation and MAPK signaling pathway resulting in apoptosis [93][94][95][96]. e cell cycle interacts with cell proliferation, differentiation, DNA damage, and cell apoptosis through TP53, GSK3B, CDKN1A, and MYC [97,98]. Collectively, these findings indicate that the core genes could be critical factors for AR-PN herb pair against the development of DN, which require to be further experimentally validated to prove this hypothesis.
Molecular docking was performed to validate the binding interactions and affinity between active components (quercetin, kaempferol, isorhamnetin, and formononetin) and the core targets (MAPK1, AKT1, GSK3B, CDKN1A, TP53, RELA, MYC, GRB2, JUN, and EGFR). It has been demonstrated that a docking score between bioactive components and targets structures less than −7, −5, and −4.25 kcal·mol −1 represents a strong binding affinity, a good binding affinity, and binding affinity, respectively [99], which means the smaller the docking score, the more stable and stronger the affinity between components and targets (ligand and receptor) binding. As shown in Table 3, the binding energies of all docking are less than −5 kcal·mol −1 , indicating that these target genes have great binding interactions and affinity with quercetin, kaempferol, formononetin, and isorhamnetin bioactive components. Moreover, the energy value of MAPK1 and JUN docked with quercetin, kaempferol, isorhamnetin, and formononetin was less than −8.0 kcal·mol −1 among the core target proteins. MAPK1 had the least binding affinity with quercetin, kaempferol, and formononetin at −8.4 kcal/mol. JUN docked with kaempferol and isorhamnetin had the lowest energy value at 9.7 kcal/mol. e docking results provide an efficient method to estimate the binding modes of AR-PN herb pair with DN-related core target proteins and provide us with sound directions for further validated.
Subsequently, ROC curves were applied to evaluate the diagnostic value of the core genes for DN. TP53, AKT1, RELA, and JUN with an AUC of 0.97436, 0.93162, 0.92308, and 0.86325, respectively, indicate that these hub genes exhibit high diagnostic values to distinguish DN patients and controls. Of course, their expression and related function also need to be further validated by mutation analysis, immune infiltration analysis, and RT-PCR.
Our study still has some limitations. Firstly, we only explored the potential functional mechanism of AR-PN herb pair on DN and have not identified the pharmacologically mutual interferences in bioactive components of AR-PN herb pair. Secondly, the DEGs of DN are acquired by analyzing published microarray data from the GEO database. Nevertheless, the sample size of DN-related DEGs is small, which needs to be enhanced in the future. Furthermore, the application of network pharmacology and bioinformatics in this study is just initial screening of the mechanisms and molecular targets of AR-PN herb pair against DN, and the results of this study have not been validated through experimental data. Hence, further studies, especially pharmacological experiments and clinical validations, are still required to be carried out to explain the complex bioactive chemical interactions. Moreover, we should provide larger groups and experimental evidence to validate these results with the interdisciplinary research works integrating computer science, network science, mathematics, and pharmacology in the future.

Conclusion
In summary, we utilized network pharmacology combined with bioinformatics to explore the mechanisms and molecular targets of AR-PN herb pair against DN in the present study. Our findings revealed that AR-PN herb pair exerts pharmacological effects on DN with various bioactive chemicals and the related pharmacological pathways, involving molecular targets and multiple cellulars, including antioxidant, anti-inflammatory, and antiapoptosis, and regulating the cell cycle.
e AR-PN herb pair and its components may be promising drugs in the treatment of DN in the future.

Data Availability
e data used to support the findings of this study are included within the article.

Disclosure
Jie Zhao and Chao Mo are joint first authors of this article.