Abstract

Objective. To explore the main bioactive compounds and investigate the underlying mechanism of Pollen Typhae (PT) against diabetic retinopathy (DR) by network pharmacology and molecular docking analysis. Methods. Bioactive ingredients and the target proteins of PT were obtained from TCMSP, and the related target genes were acquired from the SwissTargetPrediction database. The target genes of DR were obtained from GeneCards, TTD database, DisGeNET database, and DrugBank. The compound-target interaction network was established based on Cytoscape 3.7.2. The protein-protein interaction (PPI) network was constructed via STRING database and Cytoscape 3.7.2. Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were visualized through DAVID database and Bioinformatics. Ingredient-gene-pathway network analysis was conducted to further screen the ingredients, target proteins, and pathways closely related to the biological mechanism on PT for DR, and molecular docking analysis was performed by SYBYL-X 2.1.1 software. Finally, the mechanism and underlying targets of PT in the treatment of DR were predicted. Results. A total of 8 compounds and 171 intersection targets were obtained based on the online network database. 7 main compounds were screened from compound-target network, and 53 targets including the top six key targets (PTGS2, AKT1, VEGFA, MAPK3, TNF, and EGFR) were further acquired from PPI analysis. The 53 key targets covered 80 signaling pathways, among which PI3K-Akt signaling pathway, focal adhesion, Rap1 signaling pathway, VEGF signaling pathway, and HIF-1 signaling pathway were closely connected with the biological mechanism involved in the alleviation of DR by PT. Ingredient-gene-pathway network shows that AKTI, EGFR, and VEGFA were core genes, kaempferol and isorhamnetin were pivotal ingredients, and VEGF signaling pathway and Rap1 signaling pathway were closely involved in anti-DR. The docking results indicated that five main compounds (arachidonic acid, isorhamnetin, quercetin, kaempferol, and (2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one) had good binding activity with EGFR and AKT1 targets. Conclusion. The active ingredients in PT may regulate the levels of inflammatory factors, suppress the oxidative stress, and inhibit the proliferation, migration, and invasion of retinal pericytes by acting on PTGS2, AKT1, VEGFA, MAPK3, TNF, and EGFR targets through VEGF signaling pathway, PI3K-Akt signaling pathway, Rap1 signaling pathway, and HIF-1 signaling pathway to play a therapeutic role in diabetic retinopathy.

1. Introduction

Diabetic retinopathy (DR), which results from chronic high blood glucose levels and almost occurs in medium and late stage of type 1 and type 2 diabetes, is one of the most serious microvascular complications of diabetes mellitus [1]. Diabetic retinopathy can result in severe visual impairment, vitreous hemorrhage, and even blindness [2]. The World Health Organization estimated that more than 93 million people worldwide are affected by DR. The rapidly rising incidence of DR is a serious threat to public health which has been put into much attention. Moreover, DR is the main cause of adult-acquired blindness [3, 4], which dramatically affects a patient’s life and spirit and leads to a massive social economic burden. Currently, a variety of pharmaceutical and surgical treatments have been developed for DR [5]; however, their adverse side effects and limited therapeutic effect have contributed to making DR a challenge for clinical practice. Therefore, it is crucial to develop potential and effective therapeutic strategies to treat DR based on its pathogenesis.

Traditional Chinese medicine (TCM), as a type of alternative drug, is considered a multicomponent and multitarget treatment for diseases corresponding to complex organic targets. Thus, it is necessary to study pharmacological effects and molecular mechanisms of TCMs. Pollen Typhae (PT) is a classical and vital Chinese medicine with well-defined phytochemicals, which has been used for the treatment of hyperlipidemia [6], nonproliferative diabetic retinopathy [7], and vitreous hemorrhage [8]. Previous reports have shown that active ingredients of PT can inhibit the apoptosis [9], improve insulin resistance [10], and increase insulin sensitivity [11]. PT is an effective TCM for the treatment of DR. However, due to the complexity of the active ingredients of PT, its molecular mechanism has not been elucidated clearly. Therefore, the objective of this study is to elucidate the active ingredients of PT correlated to the molecular mechanism of PT on the targets of DR.

Network pharmacology, an emerged discipline of Bioinformatics, is increasingly applied in TCM research in recent years. This approach is in accordance with the holistic theory that TCM emphasizes the treatment of diseases with multiple components and multiple targets [12]. Consequently, systematic pharmacology has been recognized as a suitable method to achieve the interaction between the active ingredients and various targets and explore the underlying mechanisms against complicated disease [13, 14]. Molecular docking technology, a key-and-lock principle, is a theoretical simulation approach to investigate the binding capacity of protein receptors to small medicine molecules using computers to study the receptors’ characteristics as well as the interaction between receptors and medicinal molecules. We have referenced the previously published article [15] to investigate the mechanisms of PT on DR.

In this study, the bioactive components of PT and the specific target proteins of PT and DR were combined to construct the pharmacological networks of “components-targets-diseases” and PPI. Moreover, several biologic pathways have been proposed to be involved in DR, for example, PI3K-Akt signaling pathways, VEGF signaling pathways, HIF-1 signaling pathways, and Rap1 and TNF signaling pathways. This suggests that it is of great significance to explore the characteristics of ingredients and targets of PT on DR in biological networks and molecular docking analysis. The workflow is summarized in Figure 1.

2. Materials and Methods

2.1. Collection of Chemical Ingredients and Related Targets in PT

The chemical ingredients of PT were gathered from TCMSP database (http://old.tcmsp-e.com/tcmsp.php). In order to select the ingredients, which have better pharmacokinetic properties and oral bioavailability in vivo, the ingredients were filtrated by suggested criterion in TCMSP database, and the ingredients whose drug-likeness (DL) ≥0.18 and oral bioavailability (OB) ≥30% were regarded as putative major ingredients and retained [16]. Additionally, the SMILES number of ingredients obtained from PubChem database was imported into the SwissTargetPrediction database through Probability >0 to select the predicted targets.

2.2. Collection of the Targets of Diabetic Retinopathy

The DR-related targets were searched in the DrugBank (https://http://www.drugbank.ca/), GeneCards platform (https://www.genecards.org), Therapeutic Target Database (TTD), and DisGeNET database with “diabetic retinopathy” as the keywords [17, 18]. These collected targets were merged and duplicates were removed.

2.3. Establishment of Compound-Target Network and Analysis

To better demonstrate the anti-DR ingredients and targets of PT, Cytoscape 3.7.2 software was used to visualize the results and analyze the topological properties of the network. In the network, nodes represent ingredients or targets, and edges indicate compound-target gene interactions. In the topological analysis, the topological parameters were significant indexes to judge the importance of nodes, mainly including the degree of nodes and the betweenness centrality. Afterwards, the main active components were screened.

2.4. Protein-Protein Interaction Analysis

In order to explore the underlying mechanism of PT on DR, protein-protein interaction (PPI) networks were constructed using the STRING database and Cytoscape 3.7.2 software [19, 20]. The 171 intersection target genes were imported into the STRING database and the species were selected as “Homo sapiens” to obtain the interaction relationship between the targets, and it was saved in TSV format. Then, the file with TSV format was set into Cytoscape 3.7.2 software and the topology parameters were calculated to build an advanced protein-protein interaction network diagram and filter the core targets.

2.5. Gene Ontology and KEGG Pathway Enrichment Analysis

Gene ontology (GO) analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) are important methods to screen the key targets and enriched pathways in DAVID database (https://david.ncifcrf.gov/tools.jsp) setting  < 0.05 as the threshold [21, 22]. In this study, we chose the top 10 functional categories in biological process (BP), cellular component (CC), and molecular function (MF) and the top 20 DR-related KEGG terms, by using Bioinformatics (http://www.bioinformatics.com.cn/), to explore the key targets and the related pathways [23].

2.6. Construction of Ingredient-Gene-Pathway Network and Analysis

The ingredient-gene-pathway network was established by importing top 5 intersecting ingredients, 7 intersecting genes, and the leading 20 KEGG pathways into Cytoscape software. In this network, nodes represent ingredients, targets, and pathways, and edges indicate compound-target and target-pathway interactions. The topological parameters, mainly including the degree of nodes and the betweenness centrality, were utilized to evaluate the centrality features of nodes in ingredient-gene-pathway network.

2.7. Molecular Docking Analysis

Top five molecular compounds based on compound-target network were screened, and their molecular structures were downloaded from TCMSP. The top two key targets were selected by ingredient-gene-pathway network, and their PDB was obtained from Protein Data Bank (PDB) database. Then the PDB formats of proteins and the compounds were docked and visualized in the SYBYL-X 2.1.1 software, and the docking score was calculated to evaluate the binding degree between the bioactive components and the targets of PT. It is generally believed that the CSCORE >4.0 indicates that the docking molecules have good binding activity with the target.

3. Results

3.1. Screening of Active Ingredients and Related Targets

With the values of OB ≥30% and DL ≥0.18 properties, 8 active compounds, (2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one, arachidonic acid, isorhamnetin, quercetin, kaempferol, beta-sitosterol, testosterone palmitate, and kaempferol-3-O-α-L-rhamnosyl(1 ⟶ 2)-β-D-glucoside qt were selected from TCMSP database and 283 related targets of the compounds were screened from SwissTargetPrediction database. In addition, a total of 4077 DR-related targets were identified through the screening of the disease target database after removing the duplicate data. Then, the intersection of 171 targets was obtained.

3.2. The Compound-Target Network

In order to identify the interaction relationship between the ingredients and related targets, a compound-target network and topological analysis were made by using Cytoscape 3.7.2 and the degree values of compounds-targets were calculated. The higher the degree is, the closer the relationship between compounds and targets is, and the more important the compounds are in this network. Figure 2 shows that five compounds probably were considered as significant therapeutic compounds in DR, namely, arachidonic acid, isorhamnetin, quercetin, kaempferol, and (2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one; and the seven ingredients are listed in Table 1.

3.3. Construction and Topological Analysis of the PPI Network

The intersection targets PT and DR were inputted into the STRING database to obtain the PPI network, with 171 nodes and 1649 edges. Cytoscape 3.7.2 was used for merging and calculating topological analysis. The higher the degree is, the more important the targets are in this network. As is shown in Figure 3, the 53 key targets including the top six key targets (PTGS2, AKT1, VEGFA, MAPK3, TNF, and EGFR) were acquired, and the details of the 53 key targets are listed in Table 2. The PPI analysis demonstrated that these targets may be the key targets for PT at treating DR. In addition, the related abbreviations are provided in Supplementary Table S1.

3.4. GO Enrichment Analysis of Key Targets

GO enrichment analysis process consists of biological process (BP), cellular component (CC), and molecular function (MF). A total of GO terms, including 227 biological processes, 27 cellular components, and 55 molecular functions with  < 0.05, were obtained through inputting the 53 key genes into the DAVID website. Among them, the top ten most significantly enriched GO terms, respectively, were input into Bioinformatics (http://www.bioinformatics.com.cn/) to obtain the GO analysis diagram. As is shown in Figure 4, these targets were mainly enriched in signal transduction, cell proliferation, and apoptosis in BP, including positive regulation of transcription from RNA polymerase II promoter (GO:0045944), signal transduction (GO:0007165), negative regulation of apoptotic process (GO:0043066), positive regulation of cell proliferation (GO:0008284), and positive regulation of apoptotic process (GO:0043065). In CC analysis, the cellular components of these key targets mainly were nucleus (GO:0005634), plasma membrane (GO:0005886), extracellular space (GO:0005615), cytosol (GO:0005829), and extracellular region (GO:0005576). The MF analysis indicates that these key targets are closely related to protein binding (GO:0005515). These analyses suggested that PT may modulate angiogenesis by regulating the proliferation and apoptosis of retinal vascular cells, as well as the binding of core targets to its receptor, to treat DR. The details of GO enrichment are provided in Supplementary Table S2.

3.5. KEGG Enrichment Analysis

There were 80 KEGG signaling pathways obtained by the 53 key genes which were inputted into the DAVID database. Analysis of the KEGG pathway was carried out via Bioinformatics (http://www.bioinformatics.com.cn/), and a bubble diagram was obtained by listing the top 20 signaling pathways according to the enrichment score. As shown in Figure 5, the results of KEGG pathway enrichment included VEGF signaling pathway, PI3K-Akt signaling pathway, Rap1 signaling pathway, and HIF-1 signaling pathway, which were mainly related to oxidative stress, inflammation response, and retinal angiogenesis and proliferation. It is indicated that PT may act through these signaling pathways in treating PT. Furthermore, the top 20 signaling pathways and their specific related genes information are shown in Table 3.

3.6. The Ingredient-Gene-Pathway Network

The ingredient-gene-pathway network consists of 29 nodes, including 5 active ingredients, 7 intersecting genes, and 20 pathways. As shown in Figure 6, AKTI, EGFR, and VEGFA are the core genes. According to the topological analysis, kaempferol and isorhamnetin are pivotal ingredients and could be the key effective components of PT for the treatment of DR. In this network, several pathways are involved in anti-DR, mainly including VEGF signaling pathway and Rap1 signaling pathway. It is indicated that PT may act through AKTI, EGFR, and VEGFA via VEGF signaling pathway and Rap1 signaling pathway in treating PT.

3.7. Molecular Docking Results

Five compounds had a docking score higher than 4.0 with EGFR and AKT1, suggesting good binding activity. All docking scores are enlisted in Table 4. As shown in Figure 7(A), (2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one stably bound to AKT1 through ARG23, ASN53, LYS14, ARG86, TYR18, and ILE19; arachidonic acid stably bound to the active site of AKT1 through hydrogen bonding interactions with ASN53, ARG86, and LYS14 on the AKT1 target protein (Figure 7(B)); isorhamnetin bound to the active site of AKT1 stably through hydrogen bond interactions with LYS14, ARG25, ASN53, GLU17, TYR18, and ILE19 on the AKT1 target protein (Figures 7(C)); quercetin and kaempferol stably bound to AKT1 through ASN54, ASN53, ARG86, GLU17, and LYS14 (Figures 7((D) and (E))); (2R)-5,7-dihydroxy-2-(4-hydroxyphenyl)chroman-4-one stably bound to the active site of EGFR through THR360 and ARG40 on EGFR target protein (Figure 7(a)); arachidonic acid stably bound to the active site of EGFR through hydrogen bonding interactions with THR358 and THR330 on the EGFR target protein (Figure 7(b)); isorhamnetin bound to EGFR stably through hydrogen bond interactions with ASP323, SER324, GLU320, THR358, ASP355, THR360, and ASN328 (Figure 7(c)); quercetin stably bound to EGFR through ASN331, SER324, THR358, HIS359, and PHE351 on the EGFR target protein (Figure 7(d)); kaempferol stably bound to the active site of AKT1 through PHE357, HIS359, ASN331, and ASN328 (Figure 7(e)).

4. Discussion

Diabetic retinopathy progresses through a variety of pathophysiological pathways, including the stimulation of VEGF in the eyes, oxidative stress, inflammation response, and the activation of the hexosamine pathway. The current treatments are monotherapy with limited efficacy, such as intravitreal injection of anti-VEGFs, laser and surgery, and other methods. Pollen Typhae (PT), a multicomponent and multitarget traditional Chinese medicine (TCM), has been clinically used in the treatment of DR and has a significant effect. The previous pharmacological research demonstrated that the components of PT have good antioxidant and anti-inflammatory effects [24]. Therefore, it is of great significance to further explore the active ingredients and molecular mechanism anti-DR. In our study, combining compound-target and ingredient-gene-pathway networks, we predicted that active constituents (arachidonic acid, isorhamnetin, quercetin, and kaempferol) are closely related to PT for DR. Based on the PPI analysis and ingredient-gene-pathway network, core target proteins (PTGS2, AKT1, VEGFA, MAPK3, TNF, and EGFR) were obtained. Considering the results of the KEGG pathway enrichment analysis and ingredient-gene-pathway network, the key signaling pathways (VEGF signaling pathway, PI3K-Akt signaling pathway, Rap1 signaling pathway, and HIF-1 signaling pathway) were selected. Arachidonic acid is one of the essential polyunsaturated fatty acids in human body and plays a key role in regulating inflammatory responses. Previous report indicated that 5,6-dihydroxytrienoic lactone, a stable metabolite of arachidonic acid, could mediate microvascular dilation to contribute to the anti-DR activity [25]. Furthermore, lots of experiments have shown that arachidonic acid and its metabolites can suppress inflammation [26] and reduce free radical generation to participate in the prevention of diabetes and its microvascular complications. Pollen Typhae is particularly rich in flavonoids, especially flavanol triglycosides including derivatives of isorhamnetin, quercetin, and kaempferol. Modern pharmacological studies have confirmed that flavonoids have good effect on the antioxidant and anti‐inflammation activities, regulation of blood glucose and lipids, and regulation of cell cycle and vascular integrity. Several studies have demonstrated that isorhamnetin has good effects on attenuation of diabetes-induced oxidative stress and protection of nerves [27, 28]. Studies have certified that quercetin can protect the nerves in diabetic rat retina and prevent diabetic retinopathy in rats by inducing heme oxygenase-1 expression [29, 30]. Researchers have found that quercetin has significant inhibitory activity against the expressions of MCP-1, MMP9, and VEGF and may have a promising role on the treatment of DR [31], and it also can relieve inflammatory responses by suppression of nuclear factor kappa B in mouse. In this study, ingredient-gene-pathway network shows that kaempferol is one of the core ingredients. Many experiments and findings have confirmed that kaempferol could protect retinal ganglion cells through regulating vasohibin-1 from high-glucose-induced injury [32] and prevent retinal pigment epithelial cell damage from oxidative stress [32]. Another study has indicated that kaempferol can inhibit VEGF and PGF expression and in vitro angiogenesis under high-glucose environment [33]. Therefore, these might demonstrate that the active ingredients of PT play a crucial role in developing microvascular events and it might be promising in the prevention and treatment of DR.

The PPI network analysis illustrates that gene targets (PTGS2, AKT1, VEGFA, MAPK3, and TNF) may be the main targets of PT to DR. The ingredient-gene-pathway network shows that AKT1 and EGFR are the core target genes. AKT1 is a suppressor gene and plays a key role in DR, and pharmacological studies have demonstrated that high glucose promotes retinal endothelial cell migration through activation of AKT1 [34, 35]. The previous studies have indicated that EGFR inhibitors played the important role in attenuating inflammatory infiltration and angiogenesis in mice with diabetic retinopathy [36]. Another in vivo study indicated that it can prevent pathological retinal neovascularization via inhibiting EP (4)R-EGFR-Gab1-AKT signaling pathway [37]. Studies and experiments have certified that regulating the EGFR pathway also could inhibit the proliferation, migration, and invasion of endothelial cells and retinal pericytes of diabetic retinopathy rats [38]. Lots of studies have suggested that PTGS2 could reduce the inflammatory response by downregulating the PTGS2 level [39, 40]. The study suggests that inhibition of p-c-Jun expression could prevent retinal neuronal cell death in NMDA-induced excitotoxicity [41]. Furthermore, inhibiting the activation of the stress response kinase JNK can improve insulin action in retina [42]. Jun was determined as a target of c-Jun N-terminal kinase and plays a critical role in regulating angiogenesis in vitro [43]. Jun was proved to be involved in improving hypoxia in retinal vascular endothelial cells by regulating cooperation of c-Jun/AP-1 and HIF-1alpha [44], which is potentially important for treatment of DR. Less relevant studies have shown that ESR1 has a direct effect on DR; however, some articles previously reported that ESR1 plays a critical role in glucose transporter GLUT4 regulation and maintains glycemic homeostasis [45, 46]. A previous research indicated that MAPK3 is a key target for treatment of rats with diabetic retinopathy [47].

Functional and enrichment analysis shows that these pathways (VEGF signaling pathway, PI3K-Akt signaling pathway, Rap1 signaling pathway, and HIF-1 signaling pathway) are significantly involved in antioxidative stress and anti-inflammation and suppress the angiogenesis of retina. Study suggested that modulating the PI3K/AKT signaling pathway could ameliorate diabetes-induced proliferative retinopathy [48, 49]. In vitro evidence showed that the activation of PI3K/Akt signaling pathway was closely related to oxidative stress inhibition and proinflammatory response of microglial cells associated with diabetic retinopathy [50]. Focal adhesion signaling pathway was important for maintaining the integrity of retinal vascular structure and function [51, 52]. According to the study results, those genes enriched in focal adhesion signaling pathway were considerably upregulated [5355]. Moreover, proteomic bioinformatic analysis in DR indicated that DR development was closely related to focal adhesion and PI3K-Akt signaling pathway [56]. For the VEGF signaling pathway, one study’s results suggested that the inhibition of the VEGF signaling pathway could protect diabetic retina, and this protective effect might be related to anti-inflammatory and antioxidative activities [57]. In vivo studies indicated that downregulation of the VEGF and HIF1 signaling pathway could inhibit the proliferation and angiogenesis of high-glucose-induced human retinal endothelial cells [5860]. Additionally, lots of studies and experiments results have confirmed that downregulation of the VEGF and HIF1 signaling pathway could also prevent hypoxia-induced retinal angiogenesis [6163]. One experiment demonstrated that HIF-1 signaling pathway had a protective effect on DR [64], and it is suggested that targeting the HIF-1 signaling pathway could be a therapeutic approach to prevent development of DR.

5. Conclusion

In summary, this study is aimed at illuminating the therapeutic effect on DR through network pharmacology and molecular docking analysis. We predicted that these active ingredients in PT, including arachidonic acid, isorhamnetin, quercetin, and kaempferol, and core targets, such as PTGS2, AKT1, VEGFA, MAPK3, TNF, and EGFR, mainly intervene in core targets through VEGF signaling pathway, PI3K-Akt signaling pathway, Rap1 signaling pathway, HIF-1 signaling pathway, and other signaling pathways, thereby exerting their roles in the treatment DR.

Data Availability

The datasets in this study can be obtained from the corresponding author upon reasonable request. The authors are responsible for providing the final supplementary material files that will be published along with the article.

Disclosure

Rongrong zhou and De Jin are the co-first authors.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

Acknowledgments

This work was funded by the 2015 Traditional Chinese Medicine Scientific Research (no. 201507001-11).

Supplementary Materials

Table S1: the abbreviations. Table S2: the GO enrichment with related genes. (Supplementary Materials)