Abstract Polydatin (PD), a small natural compound originally extracted from Polygonum cuspidatum, exerts distinct biological functions in a variety of diseases. However, the action mechanism of PD has yet to be systematically explored. In this study, we firstly corroborated the druggability of PD by evaluating the medicinal properties of PD using a TCMSP server. We then conducted in silico target-prediction for PD using PharmMapper and ChemMapper, which led to the identification of 15 potential targets overlapping in both approaches. These 15 targets were subsequently evaluated by GeneMANIA, GO biological process and KEGG pathway analysis, which finally contribute to the construction of a drug-target-pathway network for PD. The network analysis revealed that these targets were mainly associated with cancer, cell growth and apoptosis, hormones and other physiological processes, outlining the pharmacological influences of PD on multiple integrated pathways involved in a particular network. __________________________________________________________________ Polydatin (PD), a small natural compound originally extracted from Polygonum cuspidatum, exerts distinct biological functions in a variety of diseases. graphic file with name c8ra03124j-ga.jpg Introduction Nowadays, Chinese herb medicine (CHM) has been widely recognized as a valuable resource for innovative drug discovery, since more and more compounds identified from CHM have shown distinct therapeutic activities in various diseases.^[28]1,2 Polydatin (PD; 3,4′,5-trihydroxystibene-3-β-mono-d-glucoside; [29]Fig. 1A) is one of the major active compounds originally extracted from the root and rhizome of Hu zhang (Polygonum cuspidatum), a typical CHM. PD could be found in red wine, peanuts, hop cones, cocoa-containing products and many daily diets,^[30]3 exhibiting a series of pharmacological activities, such as anti-oxidative, anti-inflammatory, and protective effects in multiple organs (nephron-, hepatic- and cardiovascular system).^[31]4–10 Furthermore, recent studies have shown that PD may act as an anti-tumor agent by inhibiting cell proliferation and inducing apoptosis in a variety of tumors.^[32]11,12 However, current knowledge about the medicinal characteristics as well as the action mechanism of PD is still fragmented, and its system-level mechanism warrants further investigation. Fig. 1. Chemical structure of polydatin (PD) and roadmap of PD target-identification. (A) Chemical structure of PD downloaded from the ZINC database (ZINC:04098633); (B) technology roadmap for the identification of potential PD targets using ADME evaluation, reverse docking, chemical similarity, GO processes & KEGG pathways analysis, which finally led to a PD-target-pathway network. [33]Fig. 1 [34]Open in a new tab With the advance in computational methods, bioinformatics has been applied to the target prediction for herb medicines,^[35]13 which has significantly promoted drug discovery and corresponding drug design processes, while saving money, time and manpower. Thus, as a natural active compound from CHM, PD is equally suitable for deciphering its system-level mechanism by using in silico target-prediction assisted with network pharmacological analysis. Here, we started the study with evaluating the medicinal property of PD using TCMSP server. Then, having validated the druggability of PD, we performed compound-target prediction by reverse docking and chemical similarity analysis, respectively. The overlapped targets identified by both approaches were chosen for gene ontology (GO) biological process and KEGG pathway analysis. Finally, we constructed a drug-target-pathway network to provide an overview on mechanisms of PD action, thus achieving a comprehensive understanding of PD at system-level. The technical route for identifying the potential targets of PD is shown in [36]Fig. 1B. Results ADME-related parameters of polydatin (PD) ADME-related parameters of PD were explored by TCMSP, which provides the details on important ADME-related parameters such as human oral bioavailability (OB), drug-likeness (DL), Caco-2 permeability (Caco-2), blood–brain barrier permeability (BBB) and Lipinski's rule of five (MW, Alog P, TPSA, Hdon, Hacc).^[37]14 Generally, the concept of DL has been widely used to filter out drug compounds with unsatisfactory characteristics, especially those with poor ADMET-related profiles.^[38]15 The average value of DL for Drugbank compounds ≥0.18 has been reported as a standard for screening bioactive drug compounds in systematic analysis of TCM.^[39]14,16,17 As shown in [40]Table 1, the DL value of PD was calculated to be 0.5 by TCMSP, indicating PD possessed a good drug-likeness. Pharmacological and molecular parameters of PD. MOL_ID Name MW Alog P Hdon Hacc OB (%) Caco-2 BBB DL TPSA RBN MOL013289 Polydatin 390.42 1.11 6 8 21.44 −0.9 −1.81 0.5 139.84 5 [41]Open in a new tab Lipinski's rule of five identifies the important characteristics requested to oral delivery for drug compounds, which include the molecular weight (MW) < 500 daltons (Da), the calculated log P (Clog P) < 5 (or Mlog P > 4.15), number of hydrogen-bond donors < 5 and number of hydrogen-bond acceptors < 10, respectively.^[42]18,19 Currently, the “rule of five” is commonly referred to as a “druglike” guideline in drug lead optimisation.^[43]20 As listed in [44]Table 1, our result showed that the PD's characteristics met most of the requirements, suggesting PD is a good candidate for drug development. Identification of potentially drug targets We used two different approaches as described in the method to predict the potentially drug targets. [45]Fig. 1B showed the 293 potentially protein targets of PD from over 7000 receptor-based pharmacophore models obtained using PharmMapper. Likewise, the potential 314 targets were predicted using ChemMapper. In a ChemMapper job, upon selection of the Homo sapiens, 135 protein targets met this requirement were identified. Finally, a total of 15 potentially protein targets overlapped in both sets of results were selected for further analysis ([46]Table 2). Potentially drug targets of PD identified by PharmMapper and ChemMapper. Rank PDB ID Name Target gene 1 1IVO Epidermal growth factor receptor EGFR 2 1P49 Steryl-sulfatase STS 3 1FDT Estradiol 17-beta-dehydrogenase 1 HSD17B1 4 3OS8 Estrogen receptor ESR1 5 1NDE Estrogen receptor beta ESR2 6 1H10 RAC-alpha serine/threonine-protein kinase AKT1 7 2XYG Caspase-3 CASP3 8 4GUB Casein kinase II subunit alpha CSNK2A1 9 3IXK Beta-secretase 1 BACE1 10 1H92 Tyrosine-protein kinase Lck LCK 11 1CW3 Aldehyde dehydrogenase, mitochondrial ALDH2 12 2QM9 Fatty acid-binding protein, adipocyte FABP4 13 5EA2 NAD(P)H dehydrogenase [quinone] 1 NQO1 14 1WSZ Histo-blood group ABO system transferase ABO 15 2QMY Ribosyldihydronicotinamide dehydrogenase [quinone] NQO2 [47]Open in a new tab Analysis by GeneMANIA To find additional genes that may have missed, the 15 candidate targets identified in above screening were further analyzed using GeneMANIA, which offers a simple, intuitive web interface to show the composite gene–gene functional interactions from a input gene list, and extends the gene list by finding more genes and/or gene products that are most related to them, and shows how they are networked.^[48]21 Our result revealed that among the 15 drug targets and 20 interacting proteins, 43.32% shared similar protein domains and 20.09% shared the same reaction within a pathway. Other results characteristics were co-expression, predicted, physical interactions and co-localization, respectively ([49]Fig. 2). All above results suggested that GeneMANIA represents a high accuracy of prediction algorithm. Fig. 2. Network of potential PD targets. Black protein nodes refer to protein targets, and different connecting colors represent different correlations. Functional association of protein targets were analyzed by GeneMANIA. [50]Fig. 2 [51]Open in a new tab GO and KEGG pathway analysis, and network construction To further explore the action mechanism of PD, analysis on interaction network for its 15 potential protein targets was performed by DAVID Bioinformatics Resources 6.8. According to the significance of P value, the top five processes involving in the PD action include positive regulation of nitric oxide biosynthetic process (GO:0045429), positive regulation of neuron apoptotic process (GO:0043525), regulation of phosphatidylinositol 3-kinase (PI3-K) signaling (GO:0014066), positive regulation of cell growth (GO:0030307) and response to UV-A (GO:0070141) ([52]Table 3 and [53]Fig. 3A). Some of the mechanisms mentioned above have been reported in other studies, such as the PI3-K signaling dependent mechanisms, through which a prominent nephroprotective effects against renal I/R injury has been documented previously.^[54]22 Meanwhile, these 15 targets participate in five KEGG pathways including estrogen signaling pathway (hsa04915), proteoglycans in cancer (hsa05205), prolactin signaling pathway (hsa04917), MAPK signaling pathway (hsa04010) and endometrial cancer (hsa05213) ([55]Table 4 and [56]Fig. 3B). Notably, all of these pathways were associated with cancer-related diseases, in line with the notion that PD showed the potential as anti-tumor agent.^[57]11,12 GO analysis of PD potentially targets. Term P value GO:0045429 – positive regulation of nitric oxide biosynthetic process 5.72 × 10^−4 GO:0043525 – positive regulation of neuron apoptotic process 5.72 × 10^−4 GO:0014066 – regulation of phosphatidylinositol 3-kinase signaling 0.0018695 GO:0030307 – positive regulation of cell growth 0.002164 GO:0070141 – response to UV-A 0.00249926 GO:0032355 – response to estradiol 0.0025336 GO:0051091 – positive regulation of sequence-specific DNA binding transcription factor activity 0.0033557 GO:0048015 – phosphatidylinositol-mediated signaling 0.0034186 GO:0043066 – negative regulation of apoptotic process 0.00575773 GO:0031659 – positive regulation of cyclin-dependent protein serine/threonine kinase activity involved in G1/S transition of mitotic cell cycle 0.00665179 GO:0006809 – nitric oxide biosynthetic process 0.01078828 GO:0007165 – signal transduction 0.01299481 GO:0097194 – execution phase of apoptosis 0.01326248 GO:0030520 – intracellular estrogen receptor signaling pathway 0.0173734 GO:0051000 – positive regulation of nitric-oxide synthase activity 0.01819367 GO:0050999 – regulation of nitric-oxide synthase activity 0.0214684 GO:0071392 – cellular response to estradiol stimulus 0.02554755 GO:0045907 – positive regulation of vasoconstriction 0.02636148 GO:0071364 – cellular response to epidermal growth factor stimulus 0.02717479 GO:0007611 – learning or memory 0.02961091 GO:0038128 – ERBB2 signaling pathway 0.03123184 GO:0030163 – protein catabolic process 0.03365853 GO:0048146 – positive regulation of fibroblast proliferation 0.04410908 GO:0071407 – cellular response to organic cyclic compound 0.04810051 GO:0048661 – positive regulation of smooth muscle cell proliferation 0.04889694 GO:0050727 – regulation of inflammatory response 0.0512825 GO:0006468 – protein phosphorylation 0.05398018 GO:0043154 – negative regulation of cysteine-type endopeptidase activity involved in apoptotic process 0.05603698 GO:0071260 – cellular response to mechanical stimulus 0.05761688 GO:0050729 – positive regulation of inflammatory response 0.05919432 GO:0031295 – T cell costimulation 0.06312722 GO:0042060 – wound healing 0.0646961 GO:0006919 – activation of cysteine-type endopeptidase activity involved in apoptotic process 0.06704484 GO:0046854 – phosphatidylinositol phosphorylation 0.07561013 GO:0007169 – transmembrane receptor protein tyrosine kinase signaling pathway 0.07715958 GO:0006469 – negative regulation of protein kinase activity 0.07947925 GO:0055114 – oxidation-reduction process 0.08536772 GO:0006979 – response to oxidative stress 0.08793844 GO:0030168 – platelet activation 0.09175962 GO:1901796 – regulation of signal transduction by p53 class mediator 0.0986003 [58]Open in a new tab Fig. 3. Enrichment analysis. (A) GO enrichment analysis map of potential PD targets by DAVID Bioinformatics Resources 6.8; (B) KEGG enrichment analysis map of potential PD targets by DAVID Bioinformatics Resources 6.8. [59]Fig. 3 [60]Open in a new tab KEGG analysis of PD potentially targets. Term P value hsa04915: estrogen signaling pathway 7.35 × 10^−4 hsa05205: proteoglycans in cancer 0.00551241 hsa04917: prolactin signaling pathway 0.007546708 hsa04010: MAPK signaling pathway 0.08096943 hsa05213: endometrial cancer 0.093609348 [61]Open in a new tab Based on target predicting and pathway enrichment analysis, a systematic network was constructed using Cytoscape 3.2.1. As shown in [62]Fig. 4, the interaction network contains 61 nodes and 125 edges. Interestingly, we found that EGFR, CASP3, AKT1, and ESR1 were identified as potential hub proteins in this network. Notably, the functions and pathways associated with these hub nodes play important roles in different diseases, such as cardiovascular disease and cancer. Fig. 4. Polydatin-target-pathway network. Red rectangle = polydatin; green diamond = protein targets; blue oval = pathways. [63]Fig. 4 [64]Open in a new tab Discussion Nowadays, by using the methods of network pharmacology and bioinformatics, the strengths to speculate and model pharmacokinetic, toxicity and metabolic endpoints have been greatly improved, leading to an accelerated process of drug discovery. As drug toxicity and poor pharmacokinetics are the major causes of failures during drug development, it has been widely accepted that these technical aspects should be prioritized in the drug discovery process.^[65]18 Here, we started the study with evaluating the medicinal property of PD using TCMSP server, suggesting PD is a good candidate for drug development. So far, a myriad of in silico target predicting methods have been developed and widely used in the field of drug research.^[66]23 In the present study, we employed PharmMapper and ChemMapper bioinformatic tools for target prediction, through which we identified 15 potential targets of PD overlapped by using both algorithms. Among these 15 putative targets, some have been well documented previously. For instance, the epidermal growth factor receptor (EGFR) is a transmembrane receptor that belongs to the family of receptor tyrosine kinases (RTK).^[67]24 It could be activated by a series of ligands, leading to the activation of many signaling pathways related to control proliferation, differentiation, and survival.^[68]25 In particularly, EGFR plays an important role in carcinogenesis.^[69]26,27In vitro experimental assays revealed that PD could reduce the phosphorylation of EGFR in ovarian cancer cell lines.^[70]28 Moreover, Caspase-3 (CASP3), as a member of Caspases, is an executioner that implements apoptosis.^[71]29,30 It mediates apoptosis by destroying a series of cellular substrates, and activates the degradation of DNA.^[72]30,31 E2 binds to its classical receptors-estrogen receptor 1 (ESR1) and estrogen receptor 2 (ESR2) and triggers genomic and non-genomic actions.^[73]32,33 The detailed mechanisms of E2/ERs interactions contain ligand binding, receptor dimerization and activation of estrogen response elements (ERES) in target genes, leading to activation or repression of related genes.^[74]34,35 Resveratrol, as an analogue of PD, can be used to treat infertility at low concentration, while exhibiting potential benefit to reduce proliferation of human endometrium at high doses.^[75]36 Furthermore, the GeneMANIA results ([76]Fig. 2) provided additional information on shared protein domains, pathway, co-expression, indicating these targets and their interacting proteins may have the same or similar functions. In recent years, studies on PD have revealed that it plays important roles in a wide range of cellular activities, such as anti-thrombotic formation and anti-neoplastic, strengthening the systolic and diastolic function of cardiac myocytes, improving microcirculation of important organs postshock.^[77]37–45 Indeed, our drug-target-pathway network based on the target-prediction and GO/KEGG pathway analysis showed that PD possesses multiple protein targets that involve in an array of interactive pathways. Thus, we have reasons to believe that PD is associated with a variety of pharmacological activities against many diseases. Materials and methods Evaluation of drug-likeness (DL) by TCMSP server The TCMSP server ([78]http://ibts.hkbu.edu.hk/LSPtcmsp.php) is a unique systems-level pharmacology platform for TCM that can calculate ADME-related characteristics for naturally occurring compounds.^[79]14 It contains chemicals, targets, drug-target networks and associated drug-target-disease networks. Meanwhile, it also provides a systematic evaluation about pharmacokinetic properties for natural compounds, which included drug-likeness (DL), oral bioavailability (OB), Caco-2 permeability (Caco-2), blood–brain barrier permeability (BBB) and Lipinski's rule of five (MW, Alog P, TPSA, Hdon, Hacc).^[80]14,16,17 DL is a qualitative notion used in drug design, and it also helps to optimize pharmacokinetic and pharmaceutical properties, such as chemical stability and solubility. A database-dependent model is constructed on account of the molecular descriptors and Tanimoto coefficient.^[81]16,17 Oral bioavailability (OB) is also one of the most important pharmacokinetic properties of medicine by orally because it plays an key role for the efficiency of the drug delivery to the systemic circulation in human body. Its value was evaluated by OBioavail 1.1 model in the TCMSP database.^[82]14,16 Furthermore, the Caco-2 permeability prediction model preCaco-2 could solve a problem for oral medicine, that was movement across the intestinal epithelial barrier, which determines the extent and rate of human drug absorption and eventually affects its bioavailability. This model was used in the TCMSP database.^[83]14,16 In this research, the drug name “Polydatin (PD)” was entered as the search term, and its druggability was analyzed at the molecular level. Identification of potentially drug targets by PharmMapper and ChemMapper PharmMapper, as a reverse docking server, can identify the potentially protein target candidates for the given probe small molecules using pharmacophore mapping approach.^[84]46 For a given compound, it can provide the top 300 targets, sorted by fit score in descending order. Meanwhile, the ChemMapper server can speculate potentially protein targets for small molecules by chemical similarity analysis.^[85]47 Both are powerful tools for potential drug targets prediction. Analysis by GeneMANIA GeneMANIA uses a database of organism-specific weighted networks to construct the resulting composite network, and it can analyze gene lists, generating hypotheses about gene function and prioritising genes for functional assays.^[86]21,48 After selecting Homo sapiens from the nine optional organisms, the candidate gene was entered into the search bar and the results were emerged. GO and KEGG pathway analysis, and network construction The functional annotation tool DAVID Bioinformatics Resources 6.8, as a web-based software, can be utilized to help comprehend relationships within gene expression data and provide visual systematization information on the candidate gene. Potentially protein targets were uploaded to the DAVID Bioinformatics Resources 6.8 server ([87]https://david.ncifcrf.gov/summary.jsp) following the online instructions, and GO and KEGG pathway detailed information for PD was generated and analyzed. For a deeper understanding of the complex relationships among drug, protein targets and pathways, we constructed a PD-target-pathway network, and analyzed them using Cytoscape 3.2.1. Conclusions In the present study, the DL of PD was calculated by TCMSP server, and our result showed that PD was a good drug candidate with a DL value of 0.5. Next, we performed targets prediction by using PharmMapper and ChemMapper, and the resultant 15 overlapped targets identified by using both algorithms were chosen for further analysis. Further GO and KEGG pathway analysis were performed and a drug-target-pathway network of PD action was finally constructed. Our results showed that PD exhibited many functions in anticancer, metabolic, endocrine and other physiological activities. Together, this work provides a visual systematization overview of PD involved signaling pathways and molecular mechanisms of action. Conflicts of interest The authors declare no conflicts of interest in association with this manuscript. List of abbreviations PD Polydatin CHM Chinese herb medicine OB Oral bioavailability DL Drug-likeness Caco-2 Caco-2 permeability BBB Blood–brain barrier permeability MW Molecular weight Clog P Calculated log P PI3-K Phosphatidylinositol 3-kinase EGFR Epidermal growth factor receptor RTK Receptor tyrosine kinases Supplementary Material RA-008-C8RA03124J-s001 [88]RA-008-C8RA03124J-s001.xlsx^ (15.4KB, xlsx) Acknowledgments