- Open Access
The potential effects and mechanisms of hispidulin in the treatment of diabetic retinopathy based on network pharmacology
BMC Complementary Medicine and Therapies volume 22, Article number: 141 (2022)
Diabetic retinopathy (DR), one of the most common and severe microvascular complication of diabetes mellitus (DM), is mainly caused by diabetic metabolic disorder. So far, there is no effective treatment for DR. Eriocauli Flos, a traditional Chinese herb, has been used in treating the ophthalmic diseases including DR. However, the active ingredients and molecular mechanisms of Eriocauli Flos to treat diabetic retinopathy remain elusive.
Here, the systems pharmacology model was developed via constructing network approach. 8 active components which were screened by oral bioavailability (OB ≥ 30%) and drug-likeness (DL ≥ 0.18) and 154 targets were selected from Eriocauli Flos through TCMSP database. Another 3593 targets related to DR were obtained from Genecards, OMIM, TTD, and Drugbank databases. The 103 intersecting targets of DR and Eriocauli Flos were obtained by Draw Venn Diagram. In addition, protein-protein interaction network was established from STRING database and the compound-target network was constructed by Cytoscape which screened top 12 core targets with cytoNCA module. Then the overlapping targets were analyzed by GO and KEGG enrichment. Moreover, two core targets were selected to perform molecular docking simulation. Subsequently, CCK8 assay, RT-PCR and Western blotting were applied to further reveal the mechanism of new candidate active component from Eriocauli Flos in high glucose-induced HRECs.
The results showed that the overlapping targets by GO analysis were enriched in cellular response to chemical stress, response to oxidative stress, response to reactive oxygen species, reactive oxygen species metabolic process and so on. Besides, the overlapping targets principally regulated pathways such as AGE-RAGE signaling pathway in diabetic complications, lipid atherosclerosis, fluid shear stress and atherosclerosis, and PI3K-Akt signaling pathway. Molecular docking exhibited that VEGFA and TNF-α, had good bindings to the great majority of compounds, especially the compound hispidulin. In vitro, hispidulin ameliorated high-glucose induced proliferation by down-regulating the expression of p-ERK, p-Akt, and VEGFA; meanwhile inhibited the mRNA levels of TNF-α.
In this study, through network pharmacology analysis and experimental validation, we found that hispidulin maybe has a potential targeted therapy effect for DR by decreasing the expression of p-Akt, p-ERK, and VEGFA, which resulted in ameliorating the proliferation in HRECs.
Diabetic retinopathy (DR) is a common complication of diabetes with microvascular lesions that is featured by vascular endothelial injury including pericyte loss, endothelial cell proliferation, basement membrane thickening, retinal ischemia and hypoxia, and angiogenesis [1, 2]. The major clinical pathogenesis now is mainly related to the activation of protein kinase C, the accumulation of terminal glycosylation products, and the hyperactivity of polyol pathway which may ultimately induce inflammatory reactions with upregulation of vascular endothelial growth factor (VEGF), tumor necrosis factor α (TNF-α), angiopoietin − 2 (Ang-2) and intercellular adhesion molecule-1 (ICAM-1), as well as result in the destruction of the blood-retinal barrier (BRB) and impairment of retinal vascular function [3, 4]. VEGFA is one of the potent pro-angiogenic factor that links angiogenesis and inflammation, which induces the proliferative changes associated with neovascularization in the progression of DR [5, 6]. More importantly, high glucose significantly increases the protein expression levels of VEGFA , and thereby, anti-VEGF therapy currently still serves as a prominent role in the management of DR . Eriocauli Flos, a Chinese herbal medicine, is traditionally used in ophthalmic diseases with the function of coursing wind, dissipating heat and improving eyesight in the theory of traditional Chinese medicine . The ingredients of Eriocauli Flos mainly belong to flavono, lavones, isoflavones, xanthones, naphthopyranones and phenolic acids, with pharmacology effects including anti-inflammatory and anti-bacterial [10, 11]. Hispidulin (Dinatin), one of the flavonoid compounds from Eriocauli Flos, has been reported had the anti-inflammatory, anti-cancer, anti-oxidative stress effects . For example, hispidulin inhibited HCC cell proliferation and metastasis by activating of peroxisome proliferator-activated receptor γ (PPARγ )  and another study demonstrated that hispidulin suppressed allergic inflammatory reaction by reducing histamine release and inflammatory cytokines such as tumor necrosis factor-α (TNF-α) and interleukin-4 (IL-4 ) . However, whether hispidulin has effects on high glucose-induced in vitro diabetic retina model by using HRECs cells is unclear.
Network pharmacology, combining system biology with computer technology, systematically observes the influence of drugs on diseases through network data analysis on the basis of “drug target of disease gene” interaction, reveals the molecular mechanism of drug action on diseases, as well as provides new ideas and directions for new drug research and development [15, 16]. Molecular docking, a method of simulating the interactions between proteins and small molecules, has become an important technology in the field of computer aided drug research [17, 18]. Therefore, the network pharmacology method was used to explore the targets, biological function and pathways of Eriocauli Flos against DR. Then, hispidulin was screened as the most effective compound from Eriocauli Flos by molecular docking. Furthermore, the effects of hispidulin against DR was further investigated by quantitative RT-PCR and Western blotting to test the mRNA and protein expression levels of p-Akt, p-ERK, and VEGFA. This study provided the theoretical foundation and laid the theoretical foundation for the subsequent experiment validation.
Material and methods
Hispidulin (1447-88-7) was purchased from Victory Biological Technology, China and D-glucose (G8270) was from Sigma, USA.
Screening active ingredients of Eriocauli Flos
All the compounds of Eriocauli Flos were gathered by using Traditional Chinese Medicine Systems Pharmacology Database (TCMSP) (http://sm.nwsuaf.edu.cn/lsp/tcmsp.phpphp),  and the active ingredients were filtered by oral bioavailability (OB ≥ 30%) and drug-likeness (DL ≥ 0.18) . The targets ‘s information of the active ingredients were also obtained from TCMSP and then affirmed through UniProt database (https://www.uniprot.org/ )  to be further analyzed.
Obtaining targets related to DR and intersection targets of Eriocauli Flos and DR
The targets associated with DR were predicted via GeneCards (https://www.genecards.org/),  Online Mendelian Inheritance in Man (OMIM) database (https://omim.org/),  Therapeutic Targets Database (TTD) (http://bidd.nus.edu.sg/group/ttd/ttd.asp)  and DrugBank database (https: www.drugbank.ca) . All targets were searched by key words of “diabetic retinopathy” and set specie as “homo sapiens”. At last, the total targets of DR and Intersection targets of Eriocauli Flos and DR were displayed by Draw Venn Diagram (http://bioinformatics.psb.ugent.be/webtools/Venn/) .
+Constructing “compound-target” network
To better acknowledge the interaction relationship between compounds and targets, the 8 active compounds and 103 intersection targets of Eriocauli Flos and DR were constructed into Cytoscape 3.8.0 to construct and analyze the network . Each node represents a compound or a target, and each line represents relationship between the active ingredients and targets.
Constructing protein-protein interaction (PPI) network of intersection targets of Eriocauli Flos and DR and screening key targets
Using STRING11.0 database (https://string-db.org/)  to predict the protein-protein interaction of 103 intersection targets of Eriocauli Flos and DR. To build the interaction network, minimum required interaction score > 0.9 was set and disconnected nodes in the network were hidden and 89 targets were obtained. Each node represents a protein and its structure, and each edge represents the associations between the different protein. And the PPI results were output as a TSV file which was then imported to Cytoscape 3.8.0 for further network analysis by using the cytoNCA plug-in,  and each filer got one score according to the degree, betweenness and closeness centrality. Degree is the most direct indicator to measure node centrality in network analysis and the higher the degree of a node reflects the more important it is in the network. Betweenness refers to the number of times that one node acts as the shortest bridge with two other nodes and the higher the number of times a node serves as an intermediary indicates the higher score of betweenness is. Closeness describes the average length of the shortest circuit of each node to other nodes and the closer a node is to other nodes means the higher score of closeness it is . Hence, 12 key targets were screened by the condition of the median value greater than that of three parameters in the score.
Gene ontology (GO) and the Kyoto encyclopedia of genes and gnomes (KEGG) enrichment analysis of the intersection targets
GO database (http://geneontology.org/),  including biological process, cell component, and molecular function terms, was used to describe gene functions and the results finally were input to excel for analysis. KEGG (www.kegg.jp/kegg/kegg1.html),  a knowledge database, was used to identify biological information and systematic function of target genes through obtaining significantly enriched biological pathways. Bioconductor clusterProfiler, an R package, was used to perform both enrichment analysis above by gene cluster .
The 2D structures of ligands (compounds) were downloaded from Pubchem database (https://pubchem.ncbi.nlm.nih.gov/)  and TCMSP database which then were imported to Chemoffice 2014 software for the sdf format switching to mol2 format (3D structure). The 3D structures of proteins (targets) were obtained from RCSB PDB database (http://www.rcsb.org/pdb/)  and were input to Pymol software to separate the original ligand out and remove the hydrone, phosphate radical and other inactive ligands from the proteins . The 3D structures of small molecules (ligands) and target proteins were imported into Auto Dock Vina software to acquire PDBQT format, and finally docking the molecules with the target proteins to observe the affinity and the hydrogen bond interaction was shown by using pymol .
Human retinal endothelial cells (HRECs) were obtained from BNCC and the cells previously were grown in Dulbecco’ Modified Eagle’s Medium (C11995500BT, Gibco, USA) supplemented with 10% fetal calf serum (10,270,106, Gibco, USA) and 1% (v/v) penicillin–streptomycin (SV30010, HyClone, China) at 37 °C in 5% CO2 incubator. For the hispidulin treatment study, HRECs were exposed to normal glucose (NG, 5 mmol/L) and high glucose (HG,30 mmol/L) with or without hispidulin for 48 h according to previous research about high glucose model establishment in HRECs. All the cells in the experiment were used from 3 to 6 passages.
HRECs viability assay
The viability of HRECs was detected by CCK-8 kit (C0038, Beyotime, China). Briefly, 8 × 103 cells per well were seeded into 96- well plate and cells were exposed to normal glucose (5 mM) with or without hispidulin (5-80 μM) for 24 h to assess cytotoxicity. Then 5 × 103 cells per well were seeded into 96- well plate and cells were treated with glucose (5 and 30 mM) with or without different concentrations of hispidulin (5-40 μM) for 48 h after the cells adhered for 24 h, Cell viability was measured following the manufacturer’s instructions.
RNA extraction and quantitative real-time polymerase chain reaction (RT-qPCR)
For the quantitative real-time PCR of VEGFA and TNF-a, 2.0× 105 HRECs were seeded into six-well plates. After treatment with Hispidulin for 48 h, total RNA was extracted using Trizol reagent (15596–026, Invitrogen, USA) and was reverse-transcribed into cDNA with ReverTra Ace qPCR RT Master Mix (FSQ-201, TOYOBO, Japan). Finally, 2 μl cDNA was used for fluorescence quantitative PCR analysis by using SYBR Green Realtime PCR Master Mix (QPK-201, TOYOBO, Japan). The relative CT method (△△CT) was used to calculate the relative transcriptome level, and the normalized target protein number was compared with the endogenous reference gene GAPDH. The specific primers used were as follows: human TNF-α (forward,5′-AGGCGCTCCCCAAGAAGACAG-3′; reverse,5′-AGCAGGCAGAAGAGAGCGTGGTG-3′), human VEGFA (forward,5′-GAACTTTCTGCTGTCTTGGG-3′;reverse,5′-CTTCGTGATGATTCTGCCC-3′), humanGAPDH (forward,5′-CAATGACCCCTTCATTGACC-3′;reverse,5′-GACAAGCTTCCCGTTCTCAG − 3′).
The cells were lysed with RIPA lysate (P0013C, Beyotime, China) containing a protease inhibitor, phosphatase inhibitor and PMSF. Then, lysates were centrifuged at 14,000 rpm, for 10 min at 4 C, and the supernatants were extracted, and resuspended in the 5× loading buffer (P0015, Beyotime, China) to prepare the protein samples. The electrophoresis was performed on SDS-PAGE and then the proteins were transferred to PDVF membranes (ISEQ00010, Millipore, USA). Following blocked of the membranes with 5% skim milk for 1 h and the membranes were then incubated overnight at 4 °C with primary antibodies as follows: VEGFA (ab214424, abcam, China), ERK (4695, Cell Signaling Technology, USA), p-ERK (4370, Cell Signaling Technology, USA), Akt (bsm-33,278 M, Bioss, China) and p-Akt (bsm-33,281 M, Bioss, China). Then membranes were incubated by the addition of a secondary antibody for 1 h at room temperature after another three times wash with TBST. Finally, the protein expression was detected by using ECL chemiluminescence substrate (P0018S, Beyotime, China).
All data were shown as mean ± SD and analyzed using one-way ANOVA and followed by Dunnett’ s multiple comparison test. P < 0.05 was considered as statistically significant.
Collection of targets of Eriocauli Flos and DR
A total of 8 compounds of Eriocauli Flos from TCMSP database, containing 154 targets, reached the criteria of OB ≥ 30% and DL ≥ 0.18. The goals (Table 1) including MOL001735 (Dinatin), MOL002281 (Toralactone), MOL002721 (Quercetagetin), MOL004112 (Patuletin), MOL006425 (2-benzo [1, 3] dioxol-5-yl-5,7-dimethoxy-chroman), MOL006426 (1,3,6-trihydroxy-2,5,7-trimethoxyxanthone), MOL006428 (7,30-dihydroxy-5,40,50-trimethoxyisoflavone) and MOL000098 (Quercetin). 3593 targets of diabetic retinopathy were obtained from Genecards, OMIM, TTD and Drugbank database, and 103 intersection targets were ultimately filtered from the targets of Eriocauli Flos and DR that the specific information were displayed by venn diagram in the Fig. 1.
Network between compounds and intersection targets of Eriocauli Flos and DR
The 8 compounds of Eriocauli Flos and 103 intersection targets were used to construct the compound-target network which consisted of 135 nodes (Fig. 2). Each edge represents the interaction of the compounds and targets. The total nodes and edges of the network were displayed in the Supplementary Table S1. And the network shows that one compound can associate with multiple compounds which accords with synthetic function of traditional Chinese medicine.
PPI network of intersection targets of Eriocauli Flos and DR
One hundred three intersection targets were input to STRING database to build the PPI network and then the results of PPI network were imported to cytoscape for further analysis. The network was composed of 89 nodes and 612 edges. The top 12 key targets (TP53, MAPK14, CCND1, MAPK1, AKT1, EGFR, EGF, FOS, MYC, VEGFA, TNF-α and JUN) were narrowed down by betweenness, closeness and degree parameter of cytoNCA module (Fig. 3) . The degree, betweenness and closeness of 12 core targets were listed in the Table 2.
GO and KEGG analysis of the intersection targets of Eriocauli Flos and DR
A total 2248 enrichment results were obtained from 103 intersection targets which contained 2095 biological process (BP), 49 molecular functions (MFs) and 153 cellular components (CCs). The top 10 enrichment results of BP, MFs and CCs were listed in the Fig. 4A, and the most notable biological processes were enriched in cellular response to chemical stress, response to oxidative stress, response to reactive oxygen species, reactive oxygen species metabolic process and so on. The KEGG enrichment result indicated that the genes were mostly enriched in AGE-RAGE signaling pathway in diabetic complications (hsa04933), lipid atherosclerosis (hsa05417), fluid shear stress and atherosclerosis (hsa05418) and PI3K-Akt signaling pathway (hsa04151), as shown in the Fig. 4B. And the detailed information about PI3K-Akt signaling pathway related to DR was displayed in the Fig. 5.
Molecular docking of the compounds and the core targets
We selected VEGFA and TNF-α, from the core targets, which was closely related to DR, to dock with the 8 compounds. Generally, the value of vina scores exhibits the binding activity between a protein and a compound and the lower value represents greater affinity. The results indicated that hispidulin, among the 8 compounds, had the highest affinity with VEGFA and the second highest affinity with TNF-α, and the affinity of VEGFA to each component was higher than that of the original ligand (HT7) while the affinity of TNF-α to the original ligand (VGY) was opposite (Table 3). In addition, we found that hispidulin formed hydrogen bonds with CYS-5, ASP-6 and GLU-12 in VEGFA while having two with TYR-227 and one with TYR-195 in TNF-α. The representative 3D diagram of the docking results was shown in the Fig. 6.
The effect of hispidulin on the proliferation in HRECs under high glucose condition
To evaluate the effect of hispidulin cytotoxicity on HRECs, we measured the cell viability of HRECs under normal glucose condition (5 mmol/L D-glucose) for 24 h. As shown in the Fig. 7A, the cell viability was decreased after treatment with hispidulin over 40 μM. Then, the HRECs were incubated with various concentrations of hispidulin for 48 h with normal glucose (NG, 5 mmol/L) or high glucose (HG,30 mmol/L) and the cell proliferation rate was significantly decreased from 10 μM to 40 μM (Fig. 7B). Thus, 10 μM and 20 μM of hispidulin treated in HRECs for 48 h was selected in the subsequent experiments.
The effect of hispidulin on the high glucose-induced VEGFA and TNF-α expression at the transcript level
To further examine the anti-inflammatory and anti-angiogenic effect of hispidulin on the regulation of VEGFA and TNF-α expression under high glucose conditions, we compared the transcription levels of VEGFA and TNF-α genes. The results demonstrated that the mRNA expression levels of VEGFA and TNF-α were increased after high glucose-stimulated HRECs culture for 48 h. Oppositely, the treatment of hispidulin significantly reduced the expression levels of VEGFA and TNF-α in a dose-dependent manner compared with high glucose group (Fig. 8).
The effect of hispidulin on the high glucose-induced VEGFA and phosphorylation of Akt and ERK1/2 on protein level
Through the network pharmacology analysis, the genes of VEGFA was predicted as the most crucial targets of Eriocauli Flos and DR. Furthermore, phosphorylation of Akt and ERK1/2 which play a vital role in the cell proliferation. Consequently, we found that treated with hispidulin showed down-regulating the expression of VEGFA only in the concentration of 20 μM, but decreasing phosphorylation of Akt and ERK1/2 in a dose-dependent manner (Fig. 9).
In this study, we employed a systematic network pharmacology method to explore the molecular mechanism of Eriocauli Flos in diabetic retinopathy. Through the PPI network and cytoNCA module analysis, two core genes VEGFA and TNF-α, were identified as major hub genes targets for DR, which has been paid to the key regulators in DR,  and also were validated as the biomarkers in our in vitro experiment. Furthermore, AKT1 and MAPK1 (ERK2), closely related to vascular endothelial proliferation, also existed in the core targets of PPI network which were tested in western blot experiment. Previous researches have proved that some proteins such as VEGFA, TNF-α, insulin-like growth factor 1 (IGF-1) and pigment epithelium-derived factor (PEDF) are responsible for the pathogenesis of DR, which contribute to retinal neovascularization in DR . In addition, the results of GO enrichment analysis indicated that the genes were enriched in the response to oxidative stress and inflammation, which was also consistent with the pathological mechanism of DR [40, 41]. According to KEGG enrichment, Eriocauli Flos may modulate DR by regulating multiple pathways such as the AGE-RAGE and the PI3K-Akt signaling pathway in diabetic complications. Previous reports have indicated that the interaction between AGE and its receptor RAGE leaded to loss of retinal pericytes, and triggered inflammatory response and blood vessel formation . In addition, VEGFA and VEGF receptor 2 (VEGFR2) have been implicated in direct activation through PI3K/Akt/eNOS pathway for endothelial cell proliferation, differentiation and angiogenesis . Moreover, the results of molecular docking more deeply displayed the close interaction relationship between effective components and overlap targets. Particularly, hispidulin had a higher affinity with VEGFA and TNF-α which implied that it may have a potential therapeutic effect on DR. Hence, we selected hispidulin as the underlying target and selected VEGFA and TNF-α as the biomarkers of DR to perform the experimental validation. In our work, we build model by the treatment of 30 μM D-glucose for 48 h in HRECs in the light of previous research,  but in the aspect of cell proliferation, there was no change in high glucose group compared with normal glucose group which was possibly caused by glucose toxicity and VEGFA maybe the main factor in the regulation of proliferation  which provided a clue that VEGF and its downstream may regulates vascular endothelial cell proliferation through the PI3K-Akt signaling pathway. Thus, we measured the expression of VEGFA at transcription level and protein level, and our data showed that the expression of VEGFA significantly increased compared with normal glucose group, and the treatment of hispidulin counteract the effect while the VEGFA expression on protein level was not significant change in the concentration of 10 μM hispidulin which indicated that 20 μM hispidulin may be the optimal concentration in high glucose- induced HRECs. And the phenomena above conformed with the network analysis (PPI) and results of molecular docking. For verification, we performed western blot analysis for several representative members (p-Akt and p-ERK) of the PI3K-Akt pathway which was selected from the top 10 pathways of preliminary KEGG analysis. And the results in our experiment showed that hispidulin affected HRECs proliferation and VEGFA expression by inhibiting expression of p-Akt and p-ERK, which suggested that the anti-proliferation and anti-angiogenic effect of hispidulin on HRECs under high glucose condition was mediated by the inactivation of Akt1 and Erk1/2 signaling and the results were also in accord with the results of KEGG enrichment analysis. All in all, based on the analysis of network pharmacology and molecular docking, our study revealed that hispidulin may have anti-proliferation and anti-angiogenic effect on DR.
Through network pharmacology analysis and by mining the data step by step, it was much clearer that Eriocauli Flos played a significant role in DR through multi-targets, multi-pathways approach. And hispidulin, a vital compound from Eriocauli Flos, exerts an anti-angiogenesis effect by down-regulating the expression of `p-AKT, p-ERK, and VEGF in vitro, therefore suggesting that hispidulin have potential as anti-VEGF agent in the therapy of DR. However, in our current study, it’s just simple verification of relevant targets and pathways in network pharmacology from the perspective of cell proliferation and the expression of VEGFA, and further in vivo studies are necessary to confirm the anti-angiogenic effect of hispidulin and to explore the underlying mechanisms of DR such as the effects on oxidative stress in terms of GO enrichment analysis.
Availability of data and materials
The datasets used or analyzed during the current study are available from the corresponding author.
Yang S, Zhang J, Chen L. The cells involved in the pathological process of diabetic retinopathy. Biomed Pharmacother. 2020;132:110818.
Aouiss A, Anka Idrissi D, Kabine M, Zaid Y. Update of inflammatory proliferative retinopathy: ischemia, hypoxia and angiogenesis. Curr Res Transl Med. 2019;67(2):62–71.
Cheung N, Mitchell P, Wong TY. Diabetic retinopathy. Lancet. 2010;376(9735):124–36.
Gui F, You Z, Fu S, Wu H, Zhang Y. Endothelial dysfunction in diabetic retinopathy. Front Endocrinol (Lausanne). 2020;11:591.
Zhao Y, Singh RP. The role of anti-vascular endothelial growth factor (anti-VEGF) in the management of proliferative diabetic retinopathy. Drugs Context. 2018;7:212532.
Behl T, Kotwani A. Exploring the various aspects of the pathological role of vascular endothelial growth factor (VEGF) in diabetic retinopathy. Pharmacol Res. 2015;99:137–48.
Wu MY, Yiang GT, Lai TT, Li CJ. The oxidative stress and mitochondrial dysfunction during the pathogenesis of diabetic retinopathy. Oxidative Med Cell Longev. 2018;2018:3420187.
Tan Y, Fukutomi A, Sun MT, Durkin S, Gilhotra J, Chan WO. Anti-VEGF crunch syndrome in proliferative diabetic retinopathy: a review. Surv Ophthalmol. 2021.
Wang M, Zhang Z, Cheang LC, Lin Z, Lee SM. Eriocaulon buergerianum extract protects PC12 cells and neurons in zebrafish against 6-hydroxydopamine-induced damage. Chin Med. 2011;6:16.
Qiao X, Ye G, Liu CF, Zhang ZX, Tu Q, Dong J, et al. Chemical analysis of Eriocaulon buergerianum and adulterating species by high-performance liquid chromatography with diode array detection and electrospray ionization tandem mass spectrometry. J Pharm Biomed Anal. 2012;57:133–42.
Fang JJ, Ye G, Chen WL, Zhao WM. Antibacterial phenolic components from Eriocaulon buergerianum. Phytochemistry. 2008;69(5):1279–86.
Patel K, Patel DK. Medicinal importance, pharmacological activities, and analytical aspects of hispidulin: a concise report. J Tradit Complement Med. 2017;7(3):360–6.
Han M, Gao H, Ju P, Gao MQ, Yuan YP, Chen XH, et al. Hispidulin inhibits hepatocellular carcinoma growth and metastasis through AMPK and ERK signaling mediated activation of PPARgamma. Biomed Pharmacother. 2018;103:272–83.
Kim DE, Min KJ, Kim MJ, Kim SH, Kwon TK. Hispidulin inhibits mast cell-mediated allergic inflammation through Down-regulation of histamine release and inflammatory cytokines. Molecules. 2019;24(11).
Boezio B, Audouze K, Ducrot P, Taboureau O. Network-based approaches in pharmacology. Mol Inform. 2017;36:10.
Zhang W, Huai Y, Miao Z, Qian A, Wang Y. Systems pharmacology for investigation of the mechanisms of action of traditional Chinese medicine in drug discovery. Front Pharmacol. 2019;10:743.
Dong D, Xu Z, Zhong W, Peng S. Parallelization of molecular docking: a review. Curr Top Med Chem. 2018;18(12):1015–28.
Saikia S, Bordoloi M. Molecular docking: challenges, advances and its use in drug discovery perspective. Curr Drug Targets. 2019;20(5):501–21.
Ru J, Li P, Wang J, Zhou W, Li B, Huang C, et al. TCMSP: a database of systems pharmacology for drug discovery from herbal medicines. Aust J Chem. 2014;6:13.
Li X, Xiang L, Lin Y, Tang Q, Meng F, Chen W. Computational analysis illustrates the mechanism of Qingfei Paidu decoction in blocking the transition of COVID-19 patients from mild to severe stage. Curr Gene Ther. 2021:1566–5232.
Breuza L, Poux S, Estreicher A, Famiglietti ML, Magrane M, Tognolli M, et al. The UniProtKB guide to the human proteome: Database (Oxford); 2016.
Safran M, Dalah I, Alexander J, Rosen N, Iny Stein T, Shmoish M, et al. GeneCards version 3: the human gene integrator. Database (Oxford). 2010; 2010:baq020.
Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. OMIM.org: online Mendelian inheritance in man (OMIM(R)), an online catalog of human genes and genetic disorders. Nucleic Acids Res. 2015;43(Database issue):D789–98.
Wang Y, Zhang S, Li F, Zhou Y, Zhang Y, Wang Z, et al. Therapeutic target database 2020: enriched resource for facilitating research and early development of targeted therapeutics. Nucleic Acids Res. 2020;48(D1):D1031–41.
Wishart DS, Feunang YD, Guo AC, Lo EJ, Marcu A, Grant JR, et al. DrugBank 5.0: a major update to the DrugBank database for 2018. Nucleic Acids Res. 2018;46(D1):D1074–82.
Huang XF, Zhang JL, Huang DP, Huang AS, Huang HT, Liu Q, et al. A network pharmacology strategy to investigate the anti-inflammatory mechanism of luteolin combined with in vitro transcriptomics and proteomics. Int Immunopharmacol. 2020;86:106727.
Doncheva NT, Morris JH, Gorodkin J, Jensen LJ. Cytoscape StringApp: network analysis and visualization of proteomics data. J Proteome Res. 2019;18(2):623–32.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019;47(D1):D607–13.
Tang Y, Li M, Wang J, Pan Y, Wu FX. CytoNCA: a cytoscape plugin for centrality analysis and evaluation of protein interaction networks. Biosystems. 2015;127:67–72.
Bringmann LF, Elmer T, Epskamp S, Krause RW, Schoch D, Wichers M, et al. What do centrality measures measure in psychological networks? J Abnorm Psychol. 2019;128(8):892–903.
The gene ontology C: expansion of the gene ontology knowledgebase and resources. Nucleic Acids Res. 2017;45(D1):D331–8.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Yu G, Wang LG, Han Y, He QY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS. 2012;16(5):284–7.
Wang Y, Xiao J, Suzek TO, Zhang J, Wang J, Zhou Z, et al. PubChem's BioAssay database. Nucleic Acids Res. 2012;40(Database issue):D400–12.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, et al. The Protein Data Bank. Nucleic Acids Res. 2000;28(1):235–42.
Seeliger D, de Groot BL. Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J Comput Aided Mol Des. 2010;24(5):417–22.
Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31(2):455–61.
Youngblood H, Robinson R, Sharma A, Sharma S. Proteomic biomarkers of retinal inflammation in diabetic retinopathy. Int J Mol Sci. 2019;20(19).
Capitao M, Soares R. Angiogenesis and inflammation crosstalk in diabetic retinopathy. J Cell Biochem. 2016;117(11):2443–53.
Calderon GD, Juarez OH, Hernandez GE, Punzo SM, De la Cruz ZD. Oxidative stress and diabetic retinopathy: development and treatment. Eye (Lond). 2017;31(8):1122–30.
Forrester JV, Kuffova L, Delibegovic M. The role of inflammation in diabetic retinopathy. Front Immunol. 2020;11:583687.
Zong H, Ward M, Stitt AW. AGEs, RAGE, and diabetic retinopathy. Curr Diab Rep. 2011;11(4):244–52.
Zhou L, Zhang S, Zhang L, Li F, Sun H, Feng J. MiR-199a-3p inhibits the proliferation, migration, and invasion of endothelial cells and retinal pericytes of diabetic retinopathy rats through regulating FGF7 via EGFR/PI3K/AKT pathway. J Recept Signal Transduct Res. 2021;41(1):19–31.
Biswas S, Chen S, Liang G, Feng B, Cai L, Khan ZA, et al. Curcumin analogs reduce stress and inflammation indices in experimental models of diabetes. Front Endocrinol (Lausanne). 2019;10:887.
Gao R, Zhu BH, Tang SB, Wang JF, Ren J. Scutellarein inhibits hypoxia- and moderately-high glucose-induced proliferation and VEGF expression in human retinal endothelial cells. Acta Pharmacol Sin. 2008;29(6):707–12.
We thank the Biological Laboratory of Basic School of Southwest Medical University for providing the platform to carry out the experiment and the assistance of professor Bo Chen.
This study was supported by Fund of Southwest Medical University in 2020 (Luzhou, China, Grant No:2020ZRZD013) and Fund of Key Laboratory of Medical Electrophysiology in 2021 (Luzhou, China, Grant No:KeyME-2021- 01).
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The original online version of this article was revised: the authors reported that certain authors' affiliations were assigned incorrectly.
About this article
Cite this article
Chen, Y., Sun, J., Zhang, Z. et al. The potential effects and mechanisms of hispidulin in the treatment of diabetic retinopathy based on network pharmacology. BMC Complement Med Ther 22, 141 (2022). https://doi.org/10.1186/s12906-022-03593-2
- Diabetic retinopathy
- Eriocauli Flos
- Network pharmacology
- Molecular docking