- 1Department of Anatomy, Histology and Embryology, Jinzhou Medical University, Jinzhou, China
- 2Liaoning Key Laboratory of Diabetic Cognitive and Perceptive Dysfunction, Jinzhou Medical University, Jinzhou, China
- 3Institute of Neuroscience, Laboratory Animal Department, Kunming Medical University, Kunming, China
Purpose: Diabetic retinopathy (DR) is a serious complication of diabetes mellitus, which nearly happens to all the diabetic sufferers. This study aims to identify the preliminary molecular regulation involved in the therapeutic efficacy of astragaloside IV (AS- IV) for DR.
Methods: Diabetic rat models were established and treated with AS-IV. Optical coherence tomography (OCT) and Hematoxylin-eosin (HE) staining was employed to demonstrate the histopathological changes. The main targets of AS-IV were identified by searching from public databases of traditional Chinese medicine (GeneCards, PharmMapper and Swiss Target Prediction). Besides, disease targets of DR were also obtained by integrated data from GEO datasets and predicted from public databases. Protein-protein interaction (PPI) network was constructed by Cytoscape with overlapping genes and 10 core targets were selected, on which Gene Ontology (GO) along with Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis were conducted. The interaction between AS-IV and these crucial genes were analyzed using molecular docking. RT-qPCR and western blot were used to verify the expression variation of core targets.
Results: OCT imaging and HE staining demonstrated that AS-IV administration significantly increased retinal thickness in diabetic rats, obviously alleviating DR induced histopathological changes as well as elevated blood glucose levels. 107 common targets of AS-IV and DR were determined after intersection. PPI network analysis filtered 10 hub genes potentially targeted by AS-IV, including VEGFA, CASP3, HIF1α, STAT3, CTNNB1, SRC, AKT1, EGFR, IL1β and IL6. Enrichment analysis indicated that these genes were mainly enriched in biological processes like T cell activation, epithelial cell proliferation and protein kinase B signaling, and involved in oxidative stress, apoptosis and inflammation-related pathways. The molecular docking prediction suggested that AS-IV exhibited stable binding to these core targets. In addition, mRNA levels of core targets in diabetic rats were differentially expressed before and after AS-IV treatment. Western blot further revealed that AS-IV treatment elevated DR-depressed protein levels of PI3K and AKT.
Conclusion: Our study elucidated the effect of AS-IV in attenuating retinopathy induced by diabetes in rats and preliminarily unveiled the therapeutic efficacy of AS-IV in the treatment of DR might be attributed to activation of PI3K-AKT signaling pathway.
Introduction
Diabetic retinopathy (DR) represents as the most common and serious ocular complication of diabetes, with approximately 103 million DR patients worldwide as of 2020 and an estimated number of 160 million DR patients in 2045 (Teo, et al., 2021). DR causes irreversible visual impairment of patients, and presents with clinical manifestations such as visual abnormalities, blurred vision, and even blindness, which is the leading cause of blindness in adults in countries all over the world. In addition to its impact on vision, the presence of DR implies an increased risk of life-threatening systemic vascular complications (Cheung, et al., 2010). Also, DR leads to a poor quality of life and also increases the risk of other diabetic comorbidities and mortality, imposing a serious burden on society (Kramer, et al., 2011). Presently, the main treatments for DR include retinal laser photocoagulation, drug therapy, hormone therapy, and surgery. However, these treatments may result in undesirable effects such as increased angiogenesis, elevated intraocular pressure and retinal hemorrhage (Stitt, et al., 2016; Whitcup, et al., 2018). Studies have shown that age is a critical factor influencing DR, with a higher number of DR patients in older populations (Youngblood, et al., 2019). As global aging continues to develop, it is of great significance to actively explore the pathogenesis of DR and explore effective therapeutic drugs.
Traditional Chinese medicine (TCM) has been used in the prevention and treatment of chronic diseases for nearly 2,000 years and has made indelible contributions. TCM believes that the pathogenesis of DR can be derived from deficiency of both qi and yin. On the basis, people may develop the types of yin deficiency and fire, yin deficiency and stomach heat, yin deficiency in lung and kidney, yin deficiency in liver and kidney, yin and yang loss, etc. (Ai, et al., 2020). Although the pathogenesis and symptoms of DR varies, the key to the pathogenesis of DR lies in the stagnation of the eyelid channels and nutritional deficiencies of the eyes. Astragalus is the holy medicine for supplementing qi, which is listed as top grade by Shennong Materia Medica, and has the function of tonifying qi and dredging collaterals. Previous practice shows that “Huangqi Decoction” and “Huangqi injection” exert obvious effects in the prevention and treatment of diabetes induced neuropathy (Chang, et al., 2017; Zhang, et al., 2019). Astragaloside IV (AS-IV) is the main pharmacological active component of Astragalus membranaceus, which has anti-inflammatory, hypotensive, anti-diabetes and myocardial protective properties (Gui, et al., 2013; He, et al., 2013; Zhang, et al., 2020). In addition, AS-IV significantly suppresses iron death in retinal pigment epithelial cells (Tang, et al., 2022), inhibits the formation of late glycation end products (Motomura, et al., 2009), and also reduces apoptotic retinal ganglion cell by reducing phosphorylation of ERK1/2 and inhibiting activation of NF-κB and various cytokines (Hua, et al., 2019). Despite great potentials AS-IV exhibits in the treatment of DR, the protective mechanism of AS-IV involved in DR remains ambiguous.
With the rapid development of bioinformatics, system biology and comprehensive pharmacology, network pharmacology based on the concept of “drug-disease target-pathway” can explore the complex action mechanism of drugs on human body (Hopkins, 2008). Meanwhile, network pharmacology, based on system bioinformatics, can analyze the complex target and target signal network, which helps evaluate the feasibility and applicability of TCM in treating complex diseases. In this study, to investigate the molecular mechanism underlying the therapeutic effect of AS-IV in DR, multiple databases were used to screen potential drug targets of AS-IV and disease targets of DR, so as to obtain the AS-IV-DR composite targets. Protein-protein interaction (PPI) network is constructed to analyze connection among these targets and identify the hub genes in accordance with the topological interaction between these targets. Afterwards, enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) were employed to elucidate the biological functions involved by hub genes. Finally, molecular docking techniques along with in vivo experiments were applied to verify the reliability of these crucial targets (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/903485/fphar-13-903485-HTML/image_m/fphar-13-903485-g001.jpg)
FIGURE 1. Flow chart of study design. DR, diabetic retinopathy; AS-IV, astragaloside IV; RT-qPCR, real time quantitative polymerase chains reaction; WB, western blot.
Materials and Methods
Animal Grouping and Treatment
Male Sprague Dawley (SD) rats (180–220 g) were purchased from Experimental Animal Center of Kunming Medical University. The rats were kept in a controlled room at 22 ± 2.0°C and 50% ± 5% humidity, with light exposure for about 12 h a day and free access to food and water. One week after adaptive feeding, rats were injected intraperitoneally with streptozotocin (60 mg/kg, Merck reagent, lot V900890) dissolved in a freshly prepared sodiumcitrate buffer (50 mM, pH4.5) so as to induce diabetes mellitus (Du et al., 2017). Age-matched control rats received the equivalent volume of vehicles. Seventy-2 h after injection, blood samples were collected from the tail of the rats (fasted for 12 h) for the measurement of glucose levels. Rats with fasting blood glucose levels exceeding 16.7 mmol/L for three consecutive days (Sun, et al., 2020) were considered as diabetic rats and were used as animal models for subsequent experiments. AS-IV reduced the blood glucose levels of diabetes rats in a dose-dependent manner (Du et al., 2017), higher dose treatment exhibits better hypoglycemic effects. SD rats were divided into three groups (n = 8/group): 1) control group: normal rats were treated with normal saline (10 ml/kg); 2) diabetic group: diabetic rats received normal saline treatment (10 ml/kg); 3) AS-IV group: diabetic rats were intragastrically treated with AS-IV (100 mg/kg). Six weeks after diabetic insult, rats received AS-IV treatment consecutively for 28 days. The fasting blood glucose levels of rats were measured once a week.
Optical Coherence Tomography
Rats were anesthetized with intraperitoneal injection of 3% phenobarbital (50 mg/kg). Compound topicamide eye drops were applied externally to dilate the pupil, and flufloxacin protected the cornea during and after surgery. Rats in the control group, diabetic group, and AS-IV treatment group were recorded using the Micron IV Retinal imaging camera system (Phoenix Research Laboratory, Pleasanton, CA). Finally, the retinal thickness of OCT images was measured by ImageJ software.
Tissue Harvest
The blood glucose of rats was decreased on the 7th, 14th, 21st and 28th day after administration and the retinal thickness was increased on the 28th day after administration (Liu et al., 2010). Accordingly, rats in the AS-IV group were treated with 100 mg/kg AS-IV for 28 days, and then were killed with the eyeballs harvested. The collected eyeball was placed in a Petri dish with normal saline on ice for retinal separation. The eyeball wall was incised with a scalpel at the equivalent of the ciliary flattening, and ophthalmic scissors were used to expand the incision until the anterior segment was completely removed. Then, the posterior segment wall was reversed for vitreous tilting. After a fraction of retina was separated, ophthalmic forceps was applied to separate the retina to the optic papilla, and finally the root was cut out with ophthalmic scissors. The removed retina was placed into a new EP tube.
AS-IV Related Target Prediction
The two-dimensional molecular structure of AS-IV was obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/) (Liu et al., 2010), the largest free database of chemical information in the world. Firstly, the targets of AS-IV were predicted in Genecards (https://www.genecards.org/). Meanwhile, the two-dimensional molecules of AS-IV were uploaded into Swiss Target Prediction database (http://www.swisstargetprediction.ch/) as well as into PharmMapper database (http://lilab.ecust.edu.cn/pharmMapper/) respectively for target prediction (Liu et al., 2010), with species defined as “Homo sapiens”. Additionally, Uniprot protein database (https://www.uniprot.org/) (UniProt, 2021) was used to retrieve the predicted targets, and the protein targets and genes of AS-IV were standardized and integrated.
Integration of Affymetrix Microarray Data
We searched the raw data of DR from the Comprehensive Gene Expression Database (http://www.ncbi.nlm.nih.gov/geo). The GSE12610 dataset contained 5 RNA samples extracted from the retinas of musculus including 3 diabetic samples (streptozotocin induced diabetes at 5 weeks) and 2 normal samples. In GSE28831 dataset, the retinas of Evans rats were extracted and microarray experiments were performed on days 7, 28, and 84 after STZ-induced diabetes. GPL20710 [MoGene-2_0-ST] Affymetrix Mouse Gene 2.0 ST array (mogene20st_Mm_ENTREZG_19) and GPL7294 (Agilent-014879 whole rat genome microarray 4x44KG4131F) were used to analyze datasets GSE12610 and GSE28831. Both datasets identified genes and pathways involved in DR formation. The R programming language was used for statistical analysis and graphing, and Pearson correlation analysis was performed to verify the repeatability of the data within each group. Heat maps were generated to visualize correlations among all samples from the same dataset.
Identification of Differentially Expressed Genes and Construction of DR Targets
Microarray data linear module Bioconductor (Robinson, et al., 2010) identified DEGs by comparing their expression values in the retinas of normal and diabetic mice, with adjusted p < 0.05 and log2 fold change >1 or < −1 as selection criteria. The heat map of the DEGs was plotted using heat map package in R4.1.1. In addition, DR-related target genes were searched in five public databases: Genecards database (https://www.genecards.org/) (Fishilevich, et al., 2017), OMIM database (https://omim.org/) (Amberger, et al., 2015), PharmGkb database (https://www.pharmgkb.org/) (Barbarino, et al., 2018), TTD database (http://db.idrblab.net/ttd/) (Chen, et al., 2002) and DrugBank database, with “diabetic retinopathy” as the key word. Combined with results from Uniprot database (https://www.uniprot.org/) (UniProt, 2021) and GEO datasets, a disease target database of DR was established after elimination of redundant disease targets.
Construction of Protein-Protein Interaction Networks and Enrichment Analysis on Biological Function
In order to clarify the interaction between AS-IV related targets and DR targets, we screened the overlapping genes between AS-IV related targets and DR related targets, which were then imported into the STRING database (https://string-db.org/) (von Mering, et al., 2005) for interaction analysis. The analyzed results were saved as TSV files and imported into Cytoscape 3.8.2 for network generation. The plug-in “cytoHubba” was applied for network topology analysis in Cytoscape. In addition, these plug-ins of R programming “clusterProfiler”, “org.hs.eg.db”, “EnrichPlot” and “GGploT2” were used to carry out the GO and KEGG pathways analysis of those overlapping targets. GO analysis annotates genes and gene products based on molecular function (MF), biological process (BP), and cellular composition (CC). KEGG is a useful resource for systematic analysis of gene function and related high-level genomic function information.
Molecular Docking
Molecular docking is a computer structure-based method widely used in drug discovery, which can identify new compounds of therapeutic significance, predict ligand target interactions at molecular level, and depict the structure-activity relationship (SAR) in absence of chemical structure of other target regulators in advance (Pinzi and Rastelli, 2019; Wang et al., 2020). The high-resolution crystal structure of active components and their corresponding bioactive ligands can be obtained from the protein database (PDB), a single global repository of the 3D structure of biological macromolecules and their complexes determined experimentally (Burley, et al., 2017). We used PyMOL software to dehydrate the target protein structure and separate ligands and receptors, the target protein taken as the receptor and AS-Ⅳ as the ligand. The active site of molecular docking was determined according to the coordinates of the ligand in the target protein complex, and the Gridbox coordinates and size were set according to the active pocket of the target protein. Subsequently, molecular docking was performed in Autodock vina. Finally, the interaction between protein and ligand was studied by means of binding energy scores. The obtained docking protein structures were further optimized using PyMOL software.
Hematoxylin-Eosin Staining
The collected eyeballs were fixed with 4% paraformaldehyde for 1 h, and then dehydrated with sucrose solution. Subsequently, the eyeballs were embedded in the optimal cutting temperature compound medium (Sakura tissue TEK, West Chester, PA) and then stored at −80°C fridge overnight. Retinal frozen sections were prepared at thickness of 10 μm using Leica microtome (CM1860, Leica Microsystems, Germany). These sections were subsequently treated with hematoxylin solution (Servicebio, G1005) for 1 min, after washes under running water, and they were continuingly stained by eosin solution (Servicebio, G1005) for 2–3 min. After sections dehydrated and transparentized, the fields near the optic nerve (200 μm) were captured for analysis.
RT-qPCR
Total RNA was extracted from the harvested retinal tissues using the Trizol kit (Biosharp), and cDNA was reversely transcribed using reverse transcriptase kit (Promega). Later, amplification reaction was followed with 20 μL reaction system composed of 10 μL SYBR-Green qPCR Master Mix, 1 μL cDNA Template, 7.8 μL DEPC water and 0.6 μL each of forward and reverse primers. Primers sequences are shown in Table 1. The thermal reaction cycle of PCR reaction was: one cycle of initial denaturation at 95°C for 2 min; 40 cycles of denaturation at 95°C for 15 s, annealing at 59°C for 20 s, and elongation at 60°C for 40 s. GAPDH was used as an internal reference gene, and the relative expression were calculated using the 2−ΔΔCt method.
Western Blot
Retinal tissues were lysed with 100 μL of RIPA tissue lysates, and the supernatant was discarded after centrifugation (15,000 rpm) at 4°C for 15 min. Then, the protein concentration and purification were determined by BCA method. The protein extracts were separated with agarose gel electrophoresis, and then transferred to polyvinylidene fluoride (PVDF) membranes. The PVDF membranes were placed in solution containing 5% fat-free milk powder and blocked with gentle shaking at room temperature for 1 h. Afterwards, these membranes were incubated with primary polyclonal antibodies against p-PI3K(1:1,000, rabbit, Bioss, bs-20611R) and p-AKT (1:1,000, rabbit, Bioss, bs-2720R) and β-actin monoclonal antibodies (1:5,000, mouse, Proteintech, 66009-1-Ig) overnight at 4°C. After washing three times, secondary antibodies (IgG-HRP, 1:5,000, goat anti-rabbit, Abbkine, A21031; IgG-HRP, 1:5,000, goat anti-mouse, Abbkine, A21010) were added to incubate for 1 h. After washing, the protein bands were developed with ECL chemiluminescence substrate kit (BL520B, Biosharp). β-actin was used as a loading control. Finally, densitometry analysis was performed using ImageJ software, and the relative expression levels were calculated by using the gray value of the target strip than the gray value of β -actin.
Statistical Analysis
All data in this study were processed by SPSS21.0 statistical software (IBM) or R programming. Measurement data were expressed by mean ± standard deviation
Results
The Protective Effect of AS-IV on DR Retina
We firstly verified the effects of AS-IV treatment on retinal condition after DR by employing OCT on rats in the control group, the diabetic group and the treatment group. As shown in the representative images, detailed retinal layer and fundus photographs were observed in all groups, and no significant malformations were shown in any diabetic rats (Figure 2). However, the total retinal thickness of the diabetic model rats was significantly reduced compared with that of the control group, and the fundus vascular images of the diabetic model rats were blurrier as a result of vitreous opacity and corneal whiteness of the diabetic rats (Figure 2A). Meanwhile, a large number of new vessels appeared disorderly in the fundus of the diabetic rats, which were prominently decreased in the diabetic rats treated with AS-IV (Figure 2A). Additionally, the thinner retinas with different degrees of cavities and cell loss in diabetic rats became thicker after AS-IV treatment than those in the diabetic group (Figures 2B,C, p < 0.01).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/903485/fphar-13-903485-HTML/image_m/fphar-13-903485-g002.jpg)
FIGURE 2. The effect of AS-IV treatment on retinal vessels in diabetic rats. (A) Retinal vascular images of rats in control group, diabetic group and treatment groups. (B) The images of retina in control group, diabetic group and treatment groups. (C) Quantification of retinal thickness among three groups. DM, diabetes mellitus; AS-IV, astragaloside IV. **p < 0.01.
AS-IV Alleviated DR Induced Histopathological Changes in Rats
The morphological and histopathological changes of rat retina after DR were observed by HE staining. The retinal surface in the control group was smooth. Each layer presents a clear and complete structure, and the cells in each layer are arranged orderly, evenly and densely (Figure 3A). However in DM group, the cells in each layer of retina were arranged disorderly, and ganglion cells were denatured or necrotic. In addition, the thickened optic nerve fiber layer s, nerve fiber edema, inflammatory cell infiltration and vacuole formation were observed in the inner core and outer plexus layer (Figure 3A). These alterations corresponded to the basic pathological characteristics of DR. After treatment with AS-IV, the pathological changes induced by DR were largely alleviated (Figure 3A). The efficacy of AS-IV in DR rats was also demonstrated by long-term blood glucose levels monitoring. As revealed, the fasting blood glucose levels of rats in DM group and treatment group were much higher (>16.7 mmol/L) than that of control group at the 6th week prior to AS-IV treatment (Figure 3B, p < 0.01). There was no significant change in blood glucose of rats in control group as time prolonged. After AS-IV treatment, blood glucose levels decreased a lot within the later 4 weeks, which were significantly lower than that of DM group (Figure 3B, p < 0.01).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/903485/fphar-13-903485-HTML/image_m/fphar-13-903485-g003.jpg)
FIGURE 3. Effects of AS-IV on histopathological changes and blood glucose levels of rats after DR. (A) HE staining images of retinal tissues in rats among control, DM and AS-IV groups. Scale bar = 200 μm. (B) Blood glucose levels variation of rats in control, DM and AS-IV treatment groups. RPL, retinal photoreceptor layer; RPE, retinal pigment epithelium; ONL, outer nuclear layer; GCL, ganglion cell layer; INL, inner nuclear layer; OPL, outer plexiform layer; DM, diabetes mellitus; AS-IV, astragaloside IV; w, weeks. **p < 0.01.
Targets Prediction of AS-IV and DR
A total of 30 potential targets corresponding to AS-IV were identified from GeneCards, 265 corresponding AS-IV targets were identified from PharmMapper, and 108 potential AS-IV targets were predicted by using Swiss Target Prediction. After data deduplication, 398 potential targets were retained (Supplementary Table S1). Pearson correlation test was used to further verify the repeatability of intra-group data in two GEO datasets. Based on Pearson’s correlation test, there was a correlation among the samples of normal group and a strong correlation among the samples of DM group in the GSE12610 dataset (Figure 4A). Moreover, we found that for GSE28831, a correlation existed among the samples of normal group and a strong correlation among the samples of DM group (Figure 4B). Therefore, the Affymetrix microarray dataset GSE12610 and GSE28831 were normalized. 360 DEGs were screened out of GSE12610 dataset using limma packet (corrected p < 0.05, logFC>1/logFC <-1), among which 178 down-regulated genes and 182 up-regulated genes were identified (Figure 4C). In addition, 171 DEGs were identified from the GSE28831 dataset, including 111 up-regulated genes and 60 down-regulated genes (Figure 4D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/903485/fphar-13-903485-HTML/image_m/fphar-13-903485-g004.jpg)
FIGURE 4. DEGs identification in DR through data mining out of GEO datasets. Pearson correlation analysis of (A) GSE12610 dataset and (B) GSE28831 dataset. Color reflects the strength of the correlation. The correlation coefficient between 0 and 1 indicates a positive correlation, and a negative correlation between −1 and 0. The greater the absolute value of the number is, the stronger the correlation is. (C) The volcano map showed the differential expression between the normal group and the DM group in GSE12610 dataset, with red dots representing up-regulated genes and green dots representing down-regulated genes. (D) The volcano map exhibited differential expression between the control and DR groups in the GSE28831 dataset.
Identification and Analysis of Core Targets Between AS-IV and DR
A combined analysis of two datasets in GEO database (GSE12610, GSE28831) identified 527 differentially expressed genes associated with DR. In addition, we integrated disease targets from DrugBank, GeneCards, OMIM, PharmGKB, and TTD databases, with 3 genes identified in DrugBank, 4079 genes identified in GeneCards, 85 genes predicted in OMIM, 218 genes in PharmGKB and 27 genes in TTD. Together, 4703 disease targets of DR were determined after integration and de-duplication (Supplementary Table S2). We matched 398 potential targets of AS-IV with 4703 DR target databases, intersecting 107 shared targets (Figure 5A, Supplementary Table S3). The generated PPI network displayed that the network contains 100 nodes with 961 edges (Figure 5B). Hub genes were selected using plugin “cytoHubba” with target connectivity ≥30 (triple the median), and 10 core targets (VEGFA, CASP3, HIF1α, STAT3, CTNNB1, SRC, AKT1, EGFR, IL1β and IL6) out of the overlapping genes with the largest degree value stood out for later analysis (Figure 5C, Supplementary Table S4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/903485/fphar-13-903485-HTML/image_m/fphar-13-903485-g005.jpg)
FIGURE 5. Hub genes identification and biological function analysis. (A) Intersection of targets shared by AS-IV and DR. The blue and red ellipses represent the identified DR and AS-IV targets respectively, and the middle part represents the common targets of AS-IV and DR. (B) PPI networks were visually analyzed using Cityscape 3.8.2 software. The larger the node area, the darker the color, the more important it is in the network. (C) The top 10 key genes obtained by degree value method. Each node represents a protein, and each edge represents the relationship between two proteins. The larger the node in the PPI network the darker the color, the more the node protein will interact with the surrounding proteins. (D) GO enrichment analysis of hub genes, including the BP, CC and the MF items. (E) KEGG pathways terms involved by hub genes. DR, diabetic retinopathy; AS-IV, astragaloside IV; BP, biological processes; CC, cellular component; MF, molecular function.
Biological Functions Analysis of Identified Hub Genes
To explore the biological functions of the core targets, the GO and KEGG pathways were employed. GO analysis results showed a total of 2,269 entries involved by selected targets. In biological processes (BP) enrichment analysis, 2096 items were obtained, of which the top 10 include T cell activation, epithelial cell proliferation, protein kinase B signaling, regulation of protein kinase B signaling, leukocyte cell-cell adhesion, peptidyl-serine phosphorylation, negative regulation of catabolic process, positive regulation of protein kinase B signaling, tissue remodeling, regulation of reactive oxygen species metabolic process (Figure 5D). In the cell component (CC) enrichment analysis, the top 10 out of 68 items include cell leading edge, focal adhesion, cell-substrate junction, leading edge membrane, ruffle, extrinsic component of membrane, endocytic vesicle, ruffle membrane, integrin complex, protein complex involved in cell adhesion (Figure 5D). Molecular function (MF) enrichment analysis revealed 105 items, and the top 10 terms are phosphatase binding, protein Serine/Threonine kinase activity, integrin binding, protein phosphatase binding, growth factor receptor binding, nuclear receptor activity, ligand-activated transcription factor activity, extracellular matrix binding, fibronectin binding, chemokine binding (Figure 5D). KEGG pathway analysis showed that 150 enrichment pathways, of which the most significantly enriched 20 pathways present as proteoglycans in cancer, PI3K-AKT signaling pathway, focal adhesion signaling pathway, human cytomegalovirus infection, AGE-RAGE signaling pathway in diabetic complications, gastric cancer, Kaposi sarcoma-associated herpesvirus infection, lipid and atherosclerosis, colorectal cancer, endocrine resistance, Rap1 signaling pathway, chemical carcinogenesis-receptor activation, pancreatic cancer, EGFR tyrosine kinase inhibitor resistance, fluid shear stress and atherosclerosis, Yersinia infection, Hepatitis B, prostate cancer, Th17 cell differentiation, FoxO signaling pathway (Figure 5E, Supplementary Table S5).
Molecular Docking Analysis of Crucial Targets
In the present study, the possible interactions between 10 hub genes and AS-IV were investigated by molecular docking verification. In the process of molecular model construction, target genes were applied, VEGFA (PDB:1MKG), CASP3(PDB:1QX3), HIF1α (PDB:3HQU), STAT3 (PDB:6TLC), CTNNB1 (PDB:5IVN), SRC (PDB: 4U5J), AKT1 (PDB:5AAR), EGFR (PDB:3G5X), IL1β (PDB: 4X3A), IL6 (PDB:4O9H). Their interaction patterns with key amino acids and their binding at active sites are displayed in Figure 6. The results of binding free energy revealed that the core components of AS-IV have good binding activity with the core targets (Table 2). AS-IV was well connected to the corresponding protein, and the key amino acids around it mainly acted in the form of hydrogen bond. Hydrogen bond formation minimizes the energy of small molecule receptor complexes and is therefore the most stable. The consistent results of molecular docking further demonstrated the reliability of key targets predicated by network pharmacology.
Expression Changes of Key Target Genes in Rats with DR After AS-IV Treatment
RT-qPCR was used to detect mRNA expression changes of core targets identified above in DR rats with AS-IV treatment. Among these genes, levels of CASP3, HIF1α, VEGFA, IL6, IL1β, SRC and STAT3 were significantly elevated after DR insult, which were reversed after AS-IV treatment (Figure 7, p < 0.05, p < 0.01, p < 0.001). On the contrary, AKT1 and CTNNB1 were depressed in diabetic rats, but AS-IV treatment markedly upregulated their expression in DR rats (Figure 7, p < 0.001). These results suggest that the expression of AKT1, CASP3, HIF1α, VEGFA, CTNNB1, IL6, IL1β, SRC and STAT3 might play important roles in the treatment of DR with AS-IV.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/903485/fphar-13-903485-HTML/image_m/fphar-13-903485-g007.jpg)
FIGURE 7. Expression changes of 10 core targets. DM, diabetic group; AS-IV, treatment group. DM, diabetes mellitus; AS-IV, astragaloside IV. Data were presented as mean ± SD, *p < 0.05, **p < 0.01, ***p < 0.001.
AS-IV Treatment Might Enhance Retinal Functions of DR Rats via PI3K-AKT Signaling Pathway
AKT is a key downstream signal protein of PI3K, and PI3K/AKT pathway is an important signaling pathway involved in the progression and pathogenesis of DR (Abid, et al., 2004; Cui, et al., 2018). Thus we detected the protein levels of phosphorylated AKT (p-AKT) and p-PI3K to further investigate the therapeutic mechanism of AS-IV in DR. Results of western blot demonstrated that expression of p-AKT and p-PI3K were significantly reduced in DM group compared to that of control group (Figures 8A–C, p < 0.01). After AS-IV treatment, their levels exhibited prominent increase relative to that in DM group (Figures 8A–C, p < 0.01).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/903485/fphar-13-903485-HTML/image_m/fphar-13-903485-g008.jpg)
FIGURE 8. The role of PI3K-AKT signaling in the therapeutic efficacy of AS-IV in DR. (A) Protein bands of p-AKT, p-PI3K and β-actin in control, DM and AS-IV groups. Protein levels quantification of (B) p-AKT and (C) p-PI3K among control, DM and AS-IV groups. DM, diabetes mellitus; AS-IV, astragaloside IV. Data were presented as mean ± SD, **p < 0.01.
Discussion
At present, the pathogenesis of DR is complex and remains obscure. This paper mainly combined network pharmacology analysis with experimental verification to explore the molecular mechanism of AS-Ⅳ in the treatment of DR. Our results uncovered that AS-Ⅳ can protect retinal cells against diabetes induced retinopathy, which may act on target proteins AKT1, CASP3 and HIF1α, VEGFA, IL6, IL1β, SRC and CTNNB1, and might be involved in the regulation of PI3K-AKT signaling pathway.
There is increasing evidence that AS-IV (C41H68O14, molecular weight = 784.97) is a bioactive saponin extract of Astragalus membranaceus with a wide range of pharmacological effects, including anti-inflammatory and anti-tumor functions (Qu, et al., 2009; Wang, et al., 2019a). A recent study reported multiple protective effects of AS-IV in an animal model of diabetic retinopathy, including reduced apoptosis of retinal ganglion cells, reduced phosphorylation of ERK1/2, inhibition of NF-κB and various cytokines which induce anti-inflammatory effects, and down-regulation of the enzyme aldose reductase (an enzyme involved in the polyol pathway) (Ding, et al., 2014). In another study, AS-IV significantly increased the expression of miR-128 in retinal pigment epithelial cells (RPE), inhibited the apoptosis of RPE cells and alleviated the dysfunction of RPE (Wang, et al., 2020a). AS-IV has also reported to protect retinal capillary endothelial cells from inflammation and oxidative stress induced by high glucose concentrations by lowering GLUT1 and Nox4 expression levels (Qiao, et al., 2017). Consistent in this study, we elucidate the therapeutic efficacy of AS-IV in alleviating retinal damage induced by diabetes revealed by OCT imaging technology and HE staining. All of these highlight the great potentials of AS-IV in treating DR.
To further investigate the underlying mechanism of AS-IV treatment in DR, network-based pharmacology was applied to predict the potential targets of AS-IV and DR in this study. And 107 potential targets of AS-IV for DR were obtained after intersection, among which 10 core targets (VEGFA, CASP3, HIF1α, STAT3, CTNNB1, SRC, AKT1, EGFR, IL1β and IL6) stood out after more precise filtering. RT-qPCR validation of 10 core targets revealed that AKT1, CASP3, HIF1α, VEGFA, IL6, IL1β, SRC and CTNNB1 were modulated in the process of AS-IV treatment in DR, suggesting that AS-IV can be used a potential drug to treat DR through multi targets. It has been previously reported that in proliferative diabetic retinopathy, plasma and vitreous VEGFA levels are high (Han, et al., 2020), and VEGFA permeability is a possible factor in vascular leakage, which has a significant impact on the incidence of diabetic retinopathy (Xie, et al., 2020). VEGFA is not only the most important inducer of vascular leakage, but also a key angiogenic factor (Sharma, 2020). In a previous study, vegfa-flk-1 signaling has been shown to be the main mediator of endothelial cell mitosis, survival and microvascular permeability (Zhang, et al., 2021), suggesting that VEGF can become an important target for the treatment of ocular diseases. Inhibition of HIF-1α in DR can reduce VEGF expression and inhibit neovascularization in vivo. VEGF signaling pathway induces activation of PKC and AKT1, promotes endothelial cell proliferation, and increases neovascularization permeability, which is associated with retinal exudation and macular edema in diabetic retinopathy (Wang, et al., 2020b). In addition, HIF-1α is also associated with cell proliferation and inflammation in DR retinas (Min, et al., 2021), and specific knockout of HIF-1α in Muller cells attenuated retinal neovascularization, vascular leakage, and inflammation in streptozotocin induced DR mouse models (Lin, et al., 2011). For microvascular dysfunction of DR, its anti-angiogenesis effect is exerted by inhibiting the expression of VEGF, P-ERK, P-FAK and P-SRC crosstalk (Long, et al., 2019). Reduced mRNA and protein levels of FGF1 and HIF1α inhibit retinal neovascularization in diabetes mellitus (Guan, et al., 2020). Activated SRC is involved in signaling pathways that stimulate endothelial cell survival and angiogenesis (Li, et al., 2018). In hyperglycemic retinal endothelial cells, SRC-SHP2, β-catenin ubiquitination and endothelial permeability were significantly increased (Eshaq and Harris, 2020). In addition, the SRC kinase inhibitor PP2 was found to inhibit HIF1α expression (Jiang, et al., 1997), which suggests that SRC and β-catenin are also closely related to angiogenesis and inflammation.
Mounting evidences indicate that immune mechanism involves in the pathogenesis of DR (Mesquida, et al., 2019). AKT is a regulator of multiple signaling pathways involved in cytotoxic drug-induced inflammation, apoptosis, and autophagy (Yarmohammadi, et al., 2021). It has been reported that activation of PI3K/AKT signal plays an important role in blocking oxidative stress, inhibiting pro-inflammatory response and participating in alleviating the activation of DR microglia (Jacot and Sherris, 2011; Yang, et al., 2017). The PI3K/AKT pathway is also an important signaling pathway regulating cell survival and apoptosis (Franke, et al., 2003). In the PI3K/AKT pathway, phosphatidylinositol-dependent kinase 1 (PDK1) binds to AKT and phosphorylates AKT, thereby activating p70S6K1 and ultimately inhibiting apoptosis (Zhao, et al., 2010; Chai, et al., 2018). Recently, IL-6 has become an important factor in the pathogenesis of DR. The classical IL-6 signaling is associated with anti-inflammatory function, while trans-signaling correlates with pro-inflammatory response (Zegeye, et al., 2018). In a recent study, IL-6−/− mice showed significantly lower levels of superoxide production in the retina in response to angiotensin II stimulation than wild-type mice, suggesting a link between IL-6 and oxidative stress and inflammation in retinal pathology (Robinson, et al., 2020). IL-1 β is known for its role in inducing and amplifying neuroinflammation through its own stimulation (Liu, et al., 2012). Many of the molecular and cellular abnormalities detected in DR support the role of IL-1β in driving neuroinflammation. One separate study demonstrated that IL-1β may function in the development of diabetic retinal glial cell activation and endothelial dysfunction, thus interrupting the vicious cycle triggered by IL-1β auto-stimulation, which may be a mechanism for limiting the progression of DR (Kowluru and Odenbach, 2004).
Given current knowledge of these core targets, it can be inferred that AS-IV protects animals against DR through oxidation resistance, anti-inflammation, anti-apoptosis, inhibition of new blood vessel formation and vascular leakage. Meanwhile, the results of GO enrichment analysis on key targets, showed that the regulation of various kinases activity, insulin response and vascular endothelial cells play a major role in the AS-IV therapy of DR process. The enriched KEGG pathways of these crucial genes mainly include PI3K-AKT signaling pathway, AGE-RAGE signaling pathway in diabetic complications, and EGFR tyrosine kinase inhibitor resistance, suggesting that AS-IV may play an anti-inflammatory role by activating PI3K/AKT signaling pathway. Oxidative stress and inflammation may play an important role in delaying the progression of DR (Cecilia, et al., 2019). PI3K/AKT pathway is an important signaling pathway involved in the progression and pathogenesis of DR, endothelial cell dysfunction, and regulation of VEGF level and inflammation (Abid et al., 2004; Cui et al., 2018). AKT is a key downstream signaling protein of PI3K, and activated PI3K binds to AKT, which promotes the expression and secretion of pro-inflammatory cytokines through the activation of NF-κB pathway, resulting in imbalance of cytokine secretion and a series of inflammatory reactions (Dhandapani, et al., 2005; Koorella, et al., 2014; Shi, et al., 2016). A number of studies have shown that inhibition of PI3K signal can inhibit the secretion of pro-inflammatory factors and increase the secretion of anti-inflammatory factor IL-10 when toll-like receptor- (TLR-) inflammation occurs (Fallah, et al., 2011; Bai, et al., 2014). It has also been reported that blocking the PI3K/Akt/mTOR signaling pathway can reduce cell viability by destroying the functions of Akt and RhoB, which is conducive to the treatment of DR (Zhang, et al., 2019b). A study reported that PI3K is induced by VEGF/VEGFR2, which enhances PI3K activation through phosphorylation, and VEGF can activate AKT pathway in a PI3K-dependent manner (Wang, et al., 2019b). In fundamental agreement with the previous finding, we conducted western blot to detect the protein levels changes of AKT and PI3K, and the results revealed that AS-IV treatment reversed the DR-induced inhibition on levels of AKT and PI3K, implying AS-IV might play protective efficacy through activating PI3K-AKT signaling.
Conclusion
Collectively, this study identified core targets underlying the therapeutic effects of AS-IV in DR based on comprehensive network pharmacology, and preliminarily unveiled the therapeutic effects of AS-IV might be achieved via activation PI3K/AKT pathway, providing a strong theoretical and experimental basis for the follow-up clinical application of AS-IV in the treatment of DR.
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/Supplementary Material.
Ethics Statement
The animal study was reviewed and approved by Ethics Committee of Kunming Medical University.
Author Contributions
X-ZL, Z-FZ and T-HW contributed to the conceptual framework and revised the manuscript. Z-HZ and MX is responsible for the experimental design, operation, writing and revision. CF and YH collected experimental data and performed statistical analysis.
Funding
This study was supported by China Postdoctoral Science Foundation (No. 2017M612870) and Natural Science Foundation of Liaoning Province (No. 2019-ZD-0807).
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.903485/full#supplementary-material
References
Abid, M. R., Guo, S., Minami, T., Spokes, K. C., Ueki, K., Skurk, C., et al. (2004). Vascular Endothelial Growth Factor Activates PI3K/Akt/forkhead Signaling in Endothelial Cells. Arterioscler. Thromb. Vasc. Biol. 24 (2), 294–300. doi:10.1161/01.ATV.0000110502.10593.06
Ai, X., Yu, P., Hou, Y., Song, X., Luo, J., Li, N., et al. (2020). A Review of Traditional Chinese Medicine on Treatment of Diabetic Retinopathy and Involved Mechanisms. Biomed. Pharmacother. 132, 110852. doi:10.1016/j.biopha.2020.110852
Amberger, J. S., Bocchini, C. A., Schiettecatte, F., Scott, A. F., and Hamosh, A. (2015). OMIM.org: Online Mendelian Inheritance in Man (OMIM(R)), an Online Catalog of Human Genes and Genetic Disorders. Nucleic Acids Res. 43 (Database issue), D789–D798. doi:10.1093/nar/gku1205
Bai, W., Liu, H., Ji, Q., Zhou, Y., Liang, L., Zheng, R., et al. (2014). TLR3 Regulates Mycobacterial RNA-Induced IL-10 Production through the PI3K/AKT Signaling Pathway. Cell Signal 26 (5), 942–950. doi:10.1016/j.cellsig.2014.01.015
Barbarino, J. M., Whirl-Carrillo, M., Altman, R. B., and Klein, T. E. (2018). PharmGKB: A Worldwide Resource for Pharmacogenomic Information. Wiley Interdiscip. Rev. Syst. Biol. Med. 10 (4), e1417. doi:10.1002/wsbm.1417
Burley, S. K., Berman, H. M., Kleywegt, G. J., Markley, J. L., Nakamura, H., and Velankar, S. (2017). Protein Data Bank (PDB): The Single Global Macromolecular Structure Archive. Methods Mol. Biol. 1607, 627–641. doi:10.1007/978-1-4939-7000-1_26
Cecilia, O. M., José Alberto, C. G., José, N. P., Ernesto Germán, C. M., Ana Karen, L. C., Luis Miguel, R. P., et al. (2019). Oxidative Stress as the Main Target in Diabetic Retinopathy Pathophysiology. J. Diabetes Res. 2019, 8562408. doi:10.1155/2019/8562408
Chai, X., Sun, D., Han, Q., Yi, L., Wu, Y., and Liu, X. (2018). Hypoxia Induces Pulmonary Arterial Fibroblast Proliferation, Migration, Differentiation and Vascular Remodeling via the PI3K/Akt/p70S6K Signaling Pathway. Int. J. Mol. Med. 41 (5), 2461–2472. doi:10.3892/ijmm.2018.3462
Chang, Y., Xu, J., Yan, S., Liu, Z. P., Ren, W. C., Zhang, K. X., et al. (2017). Transcriptomics of Therapeutic Effect of Huangqi Liuyi Decoction in Treating Type 2 Diabetes. Zhongguo Zhong Yao Za Zhi 42 (14), 2760–2766. doi:10.19540/j.cnki.cjcmm.20170609.003
Chen, X., Ji, Z. L., and Chen, Y. Z. (2002). TTD: Therapeutic Target Database. Nucleic Acids Res. 30 (1), 412–415. doi:10.1093/nar/30.1.412
Cheung, N., Mitchell, P., and Wong, T. Y. (2010). Diabetic Retinopathy. Lancet 376 (9735), 124–136. doi:10.1016/S0140-6736(09)62124-3
Cui, J., Gong, R., Hu, S., Cai, L., and Chen, L. (2018). Gambogic Acid Ameliorates Diabetes-Induced Proliferative Retinopathy through Inhibition of the HIF-1α/VEGF Expression via Targeting PI3K/AKT Pathway. Life Sci. 192, 293–303. doi:10.1016/j.lfs.2017.11.007
Dhandapani, K. M., Wade, F. M., Wakade, C., Mahesh, V. B., and Brann, D. W. (2005). Neuroprotection by Stem Cell Factor in Rat Cortical Neurons Involves AKT and NFkappaB. J. Neurochem. 95 (1), 9–19. doi:10.1111/j.1471-4159.2005.03319.x
Ding, Y., Yuan, S., Liu, X., Mao, P., Zhao, C., Huang, Q., et al. (2014). Protective Effects of Astragaloside IV on Db/db Mice with Diabetic Retinopathy. PLoS One 9 (11), e112207. doi:10.1371/journal.pone.0112207
Du, H., Zhang, M., Yao, K., and Hu, Z. (2017). Protective Effect of Aster Tataricus Extract on Retinal Damage on the Virtue of its Antioxidant and Anti-inflammatory Effect in Diabetic Rat. Biomed. Pharmacother. 89, 617–622. doi:10.1016/j.biopha.2017.01.179
Eshaq, R. S., and Harris, N. R. (2020). Hyperglycemia-induced Ubiquitination and Degradation of β-catenin with the Loss of Platelet Endothelial Cell Adhesion Molecule-1 in Retinal Endothelial Cells. Microcirculation 27 (2), e12596. doi:10.1111/micc.12596
Fallah, M. P., Chelvarajan, R. L., Garvy, B. A., and Bondada, S. (2011). Role of Phosphoinositide 3-Kinase-Akt Signaling Pathway in the Age-Related Cytokine Dysregulation in Splenic Macrophages Stimulated via TLR-2 or TLR-4 Receptors. Mech. Ageing Dev. 132 (6-7), 274–286. doi:10.1016/j.mad.2011.05.003
Fishilevich, S., Nudel, R., Rappaport, N., Hadar, R., Plaschkes, I., Iny Stein, T., et al. (2017). GeneHancer: Genome-wide Integration of Enhancers and Target Genes in GeneCards. Database (Oxford), 2017, bax028. doi:10.1093/database/bax028
Franke, T. F., Hornik, C. P., Segev, L., Shostak, G. A., and Sugimoto, C. (2003). PI3K/Akt and Apoptosis: Size Matters. Oncogene 22 (56), 8983–8998. doi:10.1038/sj.onc.1207115
Guan, J. T., Li, X. X., Peng, D. W., Zhang, W. M., Qu, J., Lu, F., et al. (2020). MicroRNA-18a-5p Administration Suppresses Retinal Neovascularization by Targeting FGF1 and HIF1A. Front. Pharmacol. 11, 276. doi:10.3389/fphar.2020.00276
Gui, D., Huang, J., Guo, Y., Chen, J., Chen, Y., Xiao, W., et al. (2013). Astragaloside IV Ameliorates Renal Injury in Streptozotocin-Induced Diabetic Rats through Inhibiting NF-Κb-Mediated Inflammatory Genes Expression. Cytokine 61 (3), 970–977. doi:10.1016/j.cyto.2013.01.008
Han, N., Xu, H., Yu, N., Wu, Y., and Yu, L. (2020). MiR-203a-3p Inhibits Retinal Angiogenesis and Alleviates Proliferative Diabetic Retinopathy in Oxygen-Induced Retinopathy (OIR) Rat Model via Targeting VEGFA and HIF-1α. Clin. Exp. Pharmacol. Physiol. 47 (1), 85–94. doi:10.1111/1440-1681.13163
He, Y., Du, M., Gao, Y., Liu, H., Wang, H., Wu, X., et al. (2013). Astragaloside IV Attenuates Experimental Autoimmune Encephalomyelitis of Mice by Counteracting Oxidative Stress at Multiple Levels. PLoS One 8 (10), e76495. doi:10.1371/journal.pone.0076495
Hopkins, A. L. (2008). Network Pharmacology: the Next Paradigm in Drug Discovery. Nat. Chem. Biol. 4 (11), 682–690. doi:10.1038/nchembio.118
Hua, R., Qu, L., Ma, B., Yang, P., Sun, H., and Liu, L. (2019). Diabetic Optic Neuropathy and its Risk Factors in Chinese Patients with Diabetic Retinopathy. Invest. Ophthalmol. Vis. Sci. 60 (10), 3514–3519. doi:10.1167/iovs.19-26825
Jacot, J. L., and Sherris, D. (2011). Potential Therapeutic Roles for Inhibition of the PI3K/Akt/mTOR Pathway in the Pathophysiology of Diabetic Retinopathy. J. Ophthalmol. 2011, 589813. doi:10.1155/2011/589813
Jiang, B. H., Agani, F., Passaniti, A., and Semenza, G. L. (1997). V-SRC Induces Expression of Hypoxia-Inducible Factor 1 (HIF-1) and Transcription of Genes Encoding Vascular Endothelial Growth Factor and Enolase 1: Involvement of HIF-1 in Tumor Progression. Cancer Res. 57 (23), 5328–5335.
Koorella, C., Nair, J. R., Murray, M. E., Carlson, L. M., Watkins, S. K., and Lee, K. P. (2014). Novel Regulation of CD80/CD86-Induced Phosphatidylinositol 3-kinase Signaling by NOTCH1 Protein in Interleukin-6 and Indoleamine 2,3-dioxygenase Production by Dendritic Cells. J. Biol. Chem. 289 (11), 7747–7762. doi:10.1074/jbc.M113.519686
Kowluru, R. A., and Odenbach, S. (2004). Role of Interleukin-1beta in the Pathogenesis of Diabetic Retinopathy. Br. J. Ophthalmol. 88 (10), 1343–1347. doi:10.1136/bjo.2003.038133
Kramer, C. K., Rodrigues, T. C., Canani, L. H., Gross, J. L., and Azevedo, M. J. (2011). Diabetic Retinopathy Predicts All-Cause Mortality and Cardiovascular Events in Both Type 1 and 2 Diabetes: Meta-Analysis of Observational Studies. Diabetes Care 34 (5), 1238–1244. doi:10.2337/dc11-0079
Li, P., Chen, D., Cui, Y., Zhang, W., Weng, J., Yu, L., et al. (2018). Src Plays an Important Role in AGE-Induced Endothelial Cell Proliferation, Migration, and Tubulogenesis. Front. Physiol. 9, 765. doi:10.3389/fphys.2018.00765
Lin, M., Chen, Y., Jin, J., Hu, Y., Zhou, K. K., Zhu, M., et al. (2011). Ischaemia-induced Retinal Neovascularisation and Diabetic Retinopathy in Mice with Conditional Knockout of Hypoxia-Inducible Factor-1 in Retinal Müller Cells. Diabetologia 54 (6), 1554–1566. doi:10.1007/s00125-011-2081-0
Liu, X., Ouyang, S., Yu, B., Liu, Y., Huang, K., Gong, J., et al. (2010). PharmMapper Server: a Web Server for Potential Drug Target Identification Using Pharmacophore Mapping Approach. Nucleic Acids Res. 38, W609–W614. doi:10.1093/nar/gkq300
Liu, Y., Biarnés Costa, M., and Gerhardinger, C. (2012). IL-1β Is Upregulated in the Diabetic Retina and Retinal Vessels: Cell-specific Effect of High Glucose and IL-1β Autostimulation. PLoS One 7 (5), e36949. doi:10.1371/journal.pone.0036949
Long, L., Li, Y., Yu, S., Li, X., Hu, Y., Long, T., et al. (2019). Scutellarin Prevents Angiogenesis in Diabetic Retinopathy by Downregulating VEGF/ERK/FAK/Src Pathway Signaling. J. Diabetes Res. 2019, 4875421. doi:10.1155/2019/4875421
Mesquida, M., Drawnel, F., and Fauser, S. (2019). The Role of Inflammation in Diabetic Eye Disease. Semin. Immunopathol. 41 (4), 427–445. doi:10.1007/s00281-019-00750-7
Min, J., Zeng, T., Roux, M., Lazar, D., Chen, L., and Tudzarova, S. (2021). The Role of HIF1α-PFKFB3 Pathway in Diabetic Retinopathy. J. Clin. Endocrinol. Metab. 106 (9), 2505–2519. doi:10.1210/clinem/dgab362
Motomura, K., Fujiwara, Y., Kiyota, N., Tsurushima, K., Takeya, M., Nohara, T., et al. (2009). Astragalosides Isolated from the Root of astragalus Radix Inhibit the Formation of Advanced Glycation End Products. J. Agric. Food Chem. 57 (17), 7666–7672. doi:10.1021/jf9007168
Pinzi, L., and Rastelli, G. (2019). Molecular Docking: Shifting Paradigms in Drug Discovery. Int. J. Mol. Sci. 20 (18), 4331. doi:10.3390/ijms20184331
Qiao, Y., Fan, C. L., and Tang, M. K. (2017). Astragaloside IV Protects Rat Retinal Capillary Endothelial Cells against High Glucose-Induced Oxidative Injury. Drug Des. Devel Ther. 11, 3567–3577. doi:10.2147/DDDT.S152489
Qu, Y. Z., Li, M., Zhao, Y. L., Zhao, Z. W., Wei, X. Y., Liu, J. P., et al. (2009). Astragaloside IV Attenuates Cerebral Ischemia-Reperfusion-Induced Increase in Permeability of the Blood-Brain Barrier in Rats. Eur. J. Pharmacol. 606 (1-3), 137–141. doi:10.1016/j.ejphar.2009.01.022
Robinson, M. D., McCarthy, D. J., and Smyth, G. K. (2010). edgeR: a Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 26 (1), 139–140. doi:10.1093/bioinformatics/btp616
Robinson, R., Srinivasan, M., Shanmugam, A., Ward, A., Ganapathy, V., Bloom, J., et al. (2020). Interleukin-6 Trans-signaling Inhibition Prevents Oxidative Stress in a Mouse Model of Early Diabetic Retinopathy. Redox Biol. 34, 101574. doi:10.1016/j.redox.2020.101574
Sharma, T. (2020). Evolving Role of Anti-VEGF for Diabetic Macular Oedema: from Clinical Trials to Real Life. Eye (Lond) 34 (3), 415–417. doi:10.1038/s41433-019-0590-0
Shi, Z. M., Han, Y. W., Han, X. H., Zhang, K., Chang, Y. N., Hu, Z. M., et al. (2016). Upstream Regulators and Downstream Effectors of NF-Κb in Alzheimer's Disease. J. Neurol. Sci. 366, 127–134. doi:10.1016/j.jns.2016.05.022
Stitt, A. W., Curtis, T. M., Chen, M., Medina, R. J., McKay, G. J., Jenkins, A., et al. (2016). The Progress in Understanding and Treatment of Diabetic Retinopathy. Prog. Retin Eye Res. 51, 156–186. doi:10.1016/j.preteyeres.2015.08.001
Sun, Y., Zhu, Y., Liu, X., Chai, Y., and Xu, J. (2020). Morroniside Attenuates High Glucose-Induced BMSC Dysfunction by Regulating the Glo1/AGE/RAGE axis. Cell Prolif. 53 (8), e12866. doi:10.1111/cpr.12866
Tang, X., Li, X., Zhang, D., and Han, W. (2022). Astragaloside-IV Alleviates High Glucose-Induced Ferroptosis in Retinal Pigment Epithelial Cells by Disrupting the Expression of miR-138-5p/Sirt1/Nrf2. Bioengineered 13 (4), 8240–8254. doi:10.1080/21655979.2022.2049471
Teo, Z. L., Tham, Y. C., Yu, M., Chee, M. L., Rim, T. H., Cheung, N., et al. (2021). Global Prevalence of Diabetic Retinopathy and Projection of Burden through 2045: Systematic Review and Meta-Analysis. Ophthalmology 128 (11), 1580–1591. doi:10.1016/j.ophtha.2021.04.027
UniProt, C. (2021). UniProt: the Universal Protein Knowledgebase in 2021. Nucleic Acids Res. 49 (D1), D480–D489. doi:10.1093/nar/gkaa1100
von Mering, C., Jensen, L. J., Snel, B., Hooper, S. D., Krupp, M., Foglierini, M., et al. (2005). STRING: Known and Predicted Protein-Protein Associations, Integrated and Transferred across Organisms. Nucleic Acids Res. 33 (Database issue), D433–D437. doi:10.1093/nar/gki005
Wang, T., Zhang, Z., Song, C., Sun, L., Sui, X., Qu, Q., et al. (2020a). Astragaloside IV Protects Retinal Pigment Epithelial Cells from Apoptosis by Upregulating miR-128 E-xpression in D-iabetic R-ats. Int. J. Mol. Med. 46 (1), 340–350. doi:10.3892/ijmm.2020.4588
Wang, X., Gao, Y., Tian, N., Wang, T., Shi, Y., Xu, J., et al. (2019a). Astragaloside IV Inhibits Glucose-Induced Epithelial-Mesenchymal Transition of Podocytes through Autophagy Enhancement via the SIRT-NF-Κb P65 axis. Sci. Rep. 9 (1), 323. doi:10.1038/s41598-018-36911-1
Wang, Y., Lu, Y. H., Tang, C., Xue, M., Li, X. Y., Chang, Y. P., et al. (2019b). Calcium Dobesilate Restores Autophagy by Inhibiting the VEGF/PI3K/AKT/mTOR Signaling Pathway. Front. Pharmacol. 10, 886. doi:10.3389/fphar.2019.00886
Wang, Y., Zhai, W. L., and Yang, Y. W. (2020b). Association between NDRG2/IL-6/STAT3 Signaling Pathway and Diabetic Retinopathy in Rats. Eur. Rev. Med. Pharmacol. Sci. 24 (7), 3476–3484. doi:10.26355/eurrev_202004_20806
Wang, Y-B., Wang, Q-L., Yuan, H., Sai, M., Zuo, Z-F., and Wang, T-H. (2020). Construction of the Gene Network in the Spinal Cord Injury Treated by Coptidis Rhizoma Based on Network Pharmacological and Molecular Docking. Ibrain. 6 (3), 24–33.
Whitcup, S. M., Cidlowski, J. A., Csaky, K. G., and Ambati, J. (2018). Pharmacology of Corticosteroids for Diabetic Macular Edema. Invest. Ophthalmol. Vis. Sci. 59 (1), 1–12. doi:10.1167/iovs.17-22259
Xie, W., Zhou, P., Qu, M., Dai, Z., Zhang, X., Zhang, C., et al. (2020). Corrigendum: Ginsenoside Re Attenuates High Glucose-Induced RF/6A Injury via Regulating PI3K/AKT Inhibited HIF-1a/VEGF Signaling Pathway. Front. Pharmacol. 11, 1312. doi:10.3389/fphar.2020.01312
Yang, X., Huo, F., Liu, B., Liu, J., Chen, T., Li, J., et al. (2017). Crocin Inhibits Oxidative Stress and Pro-inflammatory Response of Microglial Cells Associated with Diabetic Retinopathy through the Activation of PI3K/Akt Signaling Pathway. J. Mol. Neurosci. 61 (4), 581–589. doi:10.1007/s12031-017-0899-8
Yarmohammadi, F., Hayes, A. W., and Karimi, G. (2021). Natural Compounds against Cytotoxic Drug-Induced Cardiotoxicity: A Review on the Involvement of PI3K/Akt Signaling Pathway. J. Biochem. Mol. Toxicol. 35 (3), e22683. doi:10.1002/jbt.22683
Youngblood, H., Robinson, R., Sharma, A., and Sharma, S. (2019). Proteomic Biomarkers of Retinal Inflammation in Diabetic Retinopathy. Int. J. Mol. Sci. 20 (19). doi:10.3390/ijms20194755
Zegeye, M. M., Lindkvist, M., Fälker, K., Kumawat, A. K., Paramel, G., Grenegård, M., et al. (2018). Activation of the JAK/STAT3 and PI3K/AKT Pathways Are Crucial for IL-6 Trans-signaling-mediated Pro-inflammatory Response in Human Vascular Endothelial Cells. Cell Commun. Signal 16 (1), 55. doi:10.1186/s12964-018-0268-4
Zhang, J., Wu, C., Gao, L., Du, G., and Qin, X. (2020). Astragaloside IV Derived from Astragalus Membranaceus: A Research Review on the Pharmacological Effects. Adv. Pharmacol. 87, 89–112. doi:10.1016/bs.apha.2019.08.002
Zhang, X., Qiu, B., Wang, Q., Sivaprasad, S., Wang, Y., Zhao, L., et al. (2021). Dysregulated Serum Lipid Metabolism Promotes the Occurrence and Development of Diabetic Retinopathy Associated with Upregulated Circulating Levels of VEGF-A, VEGF-D, and PlGF. Front. Med. 8, 779413. doi:10.3389/fmed.2021.779413
Zhang, Y., Gong, G., Zhang, X., Zhou, L., Xie, H., Tian, Y., et al. (2019a). Huangqi Guizhi Wuwu Decoction for Diabetic Peripheral Neuropathy: Protocol for a Systematic Review. Med. Baltim. 98 (31), e16696. doi:10.1097/MD.0000000000016696
Zhang, Z. Z., Qin, X. H., and Zhang, J. (2019b). MicroRNA-183 Inhibition Exerts Suppressive Effects on Diabetic Retinopathy by Inactivating BTG1-Mediated PI3K/Akt/VEGF Signaling Pathway. Am. J. Physiol. Endocrinol. Metab. 316 (6), E1050–E1060. doi:10.1152/ajpendo.00444.2018
Keywords: diabetic retinopathy, astragaloside IV, bioinformatics analysis, multiple targets, PI3K-AKT signaling pathway
Citation: Zhao Z-H, Xu M, Fu C, Huang Y, Wang T-H, Zuo Z-F and Liu X-Z (2022) A Mechanistic Exploratory Study on the Therapeutic Efficacy of Astragaloside IV Against Diabetic Retinopathy Revealed by Network Pharmacology. Front. Pharmacol. 13:903485. doi: 10.3389/fphar.2022.903485
Received: 28 March 2022; Accepted: 30 May 2022;
Published: 22 June 2022.
Edited by:
Fanfan Zhou, The University of Sydney, AustraliaReviewed by:
Mingzhu Wang, Shanghai University of Traditional Chinese Medicine, ChinaBo Chen, Shengjing Hospital of China Medical University, China
Wei Wei, Affiliated Hospital of Shandong University of Traditional Chinese Medicine, China
Copyright © 2022 Zhao, Xu, Fu, Huang, Wang, Zuo and Liu. 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: Ting-Hua Wang, Wangth_email@163.com; Zhong-Fu Zuo, zuozhongfu9807@163.com; Xue-Zheng Liu, 878599033@qq.com
†These authors have contributed equally to this work